Introduction

Transmissible gastroenteritis (TGE) is a highly contagious acute digestive system disease caused by transmissible gastroenteritis virus (TGEV) infection1. This infection manifests in viral enteritis, severe watery diarrhea, loss of appetite, vomiting, dehydration and often, high mortality. Pigs of different ages can be infected with TGEV, especially piglets under two weeks of age are the most susceptible to infection, with mortality rates as high as 100%2. This disease has become a worldwide epidemic and is one of the important diseases that cause early death of piglets in various pig-raising countries, causing heavy economic losses to the global pig breeding industry3. Currently, there are no effective drugs for treating TGEV infection, and vaccination is considered to be the most effective strategy for control and prevention, whereby the vaccine stimulates the body to develop an immune response and produce protective antibodies, which in turn produce antiviral functions4. Conventional vaccines (attenuated vaccines, inactivated vaccines or triple vaccines) are mostly used in China, but they have defects such as unreliable immune effects and atavism5. Due to its frequent genetic variation, the circulating strain does not match the vaccine strain, and vaccine development lags behind2, Therefore, an efficient and rapid development platform combined with a broad spectrum vaccine development strategy will hopefully solve this problem.

TGEV is an enveloped, single-stranded positive RNA virus with a full genome length of approximately 28.5 kb organized into a 5′UTR-ORF1a/lb-S-3a-3b-E-M-N-ORF7-3′UTR structure, containing nine open reading frames (ORFs), with ORF2, ORF4, ORF5 and ORF6 encoding four structural proteins, namely, Spike (S), Nucleocapsid (N), Membrane (M), Envelope glycoprotein (E) respectively6. Among these, the S protein is the main structural protein of TGEV, located on the virus surface, functioning as the main antigen inducing the body to produce neutralizing antibodies, and plays an important biological role in virus virulence, histophilicity, receptor-binding sites, and hemagglutination, which is the main target for vaccine development7,8. However, full-length S proteins pose challenges due to potential antigenic interference and toxicity concerns when used in vaccine formulations9. In addressing this challenge and achieving immune homeostasis at the same time, studies have used immunoinformatics tools to examine the T- and B-cell dominant epitopes of the S-protein. Combining computer science and immunology, this approach leverages advanced neural networks to screen S gene sequences using various algorithms thoroughly, elucidating key epitope information for vaccine design10. Yet epitope-based vaccines encounter difficulty in eliciting a robust immune response, necessitating delivery platforms to enhance antigen immunogenicity11. Considering ferritin’s widespread use as a nanoparticle delivery platform, its 24 subunits are arranged symmetrically around a hollow core in an octahedral manner, providing a highly ordered, repetitive structure with abundant surface area for multivalent display of target sequences to enhance immune responses and vaccine effectiveness12, this study explored the tandem fusion of ferritin’s N-terminal end with epitopes to construct candidate nanoparticle vaccine. These candidate nanoparticle vaccine were constructed with a view to providing a practical new approach for the creation of new anti-TGEV vaccines and to provide theoretical and data support for the improvement of antigen design in nanoparticle vaccine research.

Materials and methods

The research employs theoretical guidance to design candidate nanoparticle vaccines, primarily involving the screening of T-cell and B-cell dominant epitopes. Ferritin self-assembled nanoparticles serve as the delivery system for constructing these candidate nanoparticle vaccines. Subsequently, the vaccine properties are analyzed, followed by immune simulation and computerized cloning. As shown in Fig. 1 (Draw with the “Biorender” website), the immune response mechanism caused by the candidate nanoparticle vaccine is briefly described.

Fig. 1
figure 1

Schematic diagram of immune response mechanism of candidate nanoparticle vaccine against TGEV S protein.

T cell epitope prediction

The TGEV-S protein cellular epitope prediction template sequence (GenBankID: ACN71196.1) was acquired in FASTA format from the NCBI database. T cell epitopes are divided into cytotoxic T cell epitopes (CTL) and helper T cell epitopes (Th). CTL epitopes were analyzed using the NetMHCpan4.1 server (https://services.healthtech.dtu.dk/service.php? NetMHCpan-4.1) and the NetMHCpanEL4.1 server (http://tools.iedb.org/mhci/)13,14. Both servers use the same screening method, employing the consensus recommendation method composed of ANN, SMM, and CombLib. These methods use artificial neural networks to predict the binding of known epitopes to MHC molecules, identifying strong binding sequences (5%) within the same allele as potential candidates for sequences of length 14. Th epitopes were analyzed using the NetMHCIIpan-4.0 server (https://services.healthtech.dtu.dk/service.php? NetMHCIIpan-4.0)15, which also uses an artificial neural network to predict potential binding peptide segments of any known MHC-II molecule sequence and outputs the likelihood of the peptide segment naturally existing with MHC-II alleles, identifying strong binding epitopes (10%) of length 15. MHC-II molecules activate macrophages to recognize and clear intracellular pathogens by releasing IFN-γ. The secretion of cytokine IL-4 promotes the proliferation and differentiation of antigen-presenting cells. Therefore, Epitope evaluation included IFN epitope analysis (https://webs.iiitd.edu.in/raghava/ifnepitope/index.php) and IL4pred (http://crdd.osdd.net/raghava/il4pred/) for IFN⁃γ and IL inducibility16,17.

B cell epitope prediction

Linear B-cell epitopes were analyzed by the ABCpred server (https://webs.iiitd.edu.in/raghava/abcpred/ABC_submission.html)18, a trained recurrent neural network server that grades predicted B-cell epitopes based on their quality. The analysis was conducted with a sequence length set to 16, and a prediction threshold of 0.8 was applied to ensure robustness. Additionally, to ensure high confidence in the results, the entire amino acid sequence was analyzed using the PEPTIDES server (http://imed.med.ucm.es/Tools/antigenic.html)19. The PEPTIDES server employs the Kolaskar and Tongaonkar methods to predict potential antibody responses by simulating 3D epitopes of immunogenicity recognition, reporting fragments only when the minimum number of residues is 8.

Screening and conservativeness analysis of candidate epitopes

After applying the filters above, the outcomes from various servers were compared, focusing on the overlapping regions of epitopes. The immunogenicity of these overlapping epitopes was assessed using the IEDBMHC-I server (http://tools.iedb.org/immunogenicity/)20, while toxicity was evaluated via the ToxinPred server (https://webs.iiitd.edu.in/raghava/toxinpred/ design.php) 21. Epitopes with an immunogenicity score > 0.25 and deemed non-toxic were designated as candidate epitopes. Candidate epitopes were evaluated for conservation with 45 TGEV-S protein amino acid sequences downloaded from the NCBI database via the IEDB server (http://tools.iedb.org/conservancy/)22. At the same time, the conservation of amino acid residues in dominant epitopes is plotted through Weblogo.

Vaccine design and evaluation

CTL, Th, and B cell dominant epitopes were tandem linked using AAY, GPGPG, and KK flexlinker, respectively. AAY provides flexibility to protect proteins from degradation, GPGPG facilitates immunoprocessing and presentation, while KK, a dual lysyl linker, retains epitope-independent immunoreactivity. Ferritin is a self-assembled nanoparticle whose highly ordered and repetitive symmetric structure serves as a nanoplatform for antigen delivery and immunostimulation, contributing to the induction of efficient immune responses and enhancing vaccine efficacy. The N-terminal of ferritin (GenBank: WP_000949190) was linked to the dominant epitope construct via EAAAK to generate a candidate nanoparticle vaccine. The antigenicity of candidate nanoparticle vaccines was assessed using the ANTIGENpro server (http://scratch.proteomics.ics.uci.edu/)23, while solubility was evaluated using the SOLpro overexpression server (http://scratch.proteomics.ics.uci.edu/index.html)24. AllerTop Sensitization Server (http://www.ddg-pharmfac.net/AllerTOP/)25, ProtScale Hydrophobicity Server (https://web.expasy.org/cgi-bin/protscale/protscale.pl#opennewwindow)26, TMHMM-2.0 Transmembrane Structural Domain Server (https://services.healthtech.dtu.dk/services/TMHMM-2.0/) for comprehensive analysis27, and the ExPASy Physicochemical Properties Server (https://web.expasy.org/protparam/) for evaluating various computational parameters including relative molecular mass, theoretical pI, atomic composition, estimated half-life, extinction coefficient, instability index, fat index and GRAVY28.

Vaccine secondary structure prediction and tertiary structure modeling and refinement, disulfide engineering and validation

Vaccine secondary structures were analyzed using the SOPMA server (https://npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl?page=npsa_sopma.html) and PSIPRED 4.0 server (http://bioinf.cs.ucl.ac.uk /psipred) to analyze transmembrane helices, folding and domain identification29,30. Tertiary structure modeling was conducted using the homology modeling trRosetta server (https://yanglab.nankai.edu.cn/trRosetta/) to generate primary 3D models31. Subsequently, the primary 3D model pdb files were subjected to molecular dynamics simulations performing repetitive structural perturbations32 through the GalaxyRefine server (http://galaxy.seoklab.org/cgi-bin/submit.cgi?), aimed at refining errors in the primary initial models. The refined model was further stabilized by disulfide engineering through the DbD2 server (http://cptweb.cpt.wayne.edu/DbD2/index.php) to enhance structural stability33. The models were evaluated before and after refinement using the CHECK Serve (https://saves.mbi.ucla.edu/)34, with the final models subjected to validation through Z-score generation via the ProSA-web server (https://prosa.services.came.sbg.ac.at/prosa.php)35.

Analysis of conformational B-cell epitopes and molecular docking

Conformational B-cell epitopes, comprising both discontinuous and continuous epitopes formed through protein folding, were evaluated in the final refined candidate nanoparticle vaccine tertiary structure using the ElliPro server (http://tools.iedb.org/ellipro/)36. A threshold score > 0.5 and a maximum distance of 6 Angstrom were employed to assess the presence of these epitopes. Also, in order to properly stimulate an immune response, interactions between antigenic and immune receptor molecules are essential. TLR3, known for regulating immune responses, promoting immune cell maturation and differentiation, and participating in pro-inflammatory factor responses, plays a vital role in intrinsic antiviral immune responses. Conversely, TLR4 acts as a key pattern recognition receptor, inducing cytokines and bridging innate and acquired immunity, contributing to the host cell’s resistance against viral infections. Protein–protein docking studies were conducted using the pyDockWEB server (https://life.bsc.es/pid/pydockweb) to analyze the interaction patterns between the candidate nanoparticle vaccine and TLR-3 (PDB ID: 2A0Z) and TLR-4 (PDB ID: 4G8A)37. Docking complexes were visualized using the PyMOL tool, while complex 2D interworking interfaces were visualized by the LigPLot + tool.

Molecular dynamics simulation

Molecular dynamics simulations were performed using the Gromacs 2022.3 software to examine the stability of the candidate vaccine docking complex38,39. The GAFF force field and hydrogenation operation were added to each complex and RESP potential was calculated. The Tip3p water model was used for step energy minimization and balance protocol. The simulation environment of 300 K and 1 Bar was selected and the trajectory analysis was carried out in 100 ns time to calculate the physical quantities of RMSD, RMSF, Rg, SASA and hydrogen bond of the complex. In addition, the Gibbs free energy topography was evaluated by using the built-in “gmx_sham” function and "xpm2txt.py" foot.

Immunosimulation and computerized cloning

The immune simulation of the candidate nanoparticle vaccine was performed using the C-ImmSim server (https://150.146.2.1/C-IMMSIM/index.php)40, an agent-based computational immune response simulator that utilizes position-specific scoring matrices (PSSM) and machine learning methods, respectively, to predict epitope immune interactions. To mimic commercial vaccine dosing intervals, the immunization simulation retained default parameters, with the number of simulation steps set to 1050 and time steps adjusted to 1, 84, and 170 (8 h per step) for three immunizations. In order to express foreign genes in host organisms, codon optimization is necessary depending on the particular host organism; therefore, candidate nanoparticle vaccine sequences were submitted to the JCat server (http:/jcat.de/)41 for codon optimization using E. coli K12 as the host. ExpOptimizer server (https://www.novopro.cn/tools/codon-optimization.html) (NOVOPRO, China) validated Codon Adaptation Index (CAI) values against GC content. Finally, the adapted nucleotide sequences were cloned into the pET28-a( +) expression vector using SnapGene software.

Result

Results of the TGEV S protein dominant epitope screening

From the 1448 amino acids of S protein obtained from NCBI (The complete amino acid sequence is shown in supplemental table S1), CTL dominant epitopes were identified using the NetMHCpan-4.1 and NetMHCpanEL4.1 dual servers, revealing three overlapping dominant epitope regions. The dominant epitopes were analyzed via the NetMHCIIpan-4.0 server and further evaluated by IFNepitope and IL4pred server, identifying three overlapping dominant epitope regions. Screening with the ABCpred and PEPTIDES servers revealed six overlapping linear B-cell epitope regions. All identified dominant epitopes were confirmed as non-toxic peptides by the ToxinPred server and highly immunogenic (immunogenicity score > 0.25) by the IEDBMHC-I server. Conservation analysis conducted by the IEDB server indicated high conservation of linear B-cell dominant epitopes with Th-2 and Th-3 across all 45 sequences (The S protein GenBankID of 45 strains was shown in supplemental table S2), while CTL dominant epitopes showed less conservation with Th-1 (Detailed scores for conservative analysis are shown in supplementary table S3). Nonetheless, considering their significance in vaccine molecular design, they were included in the study. The results are shown in Table 1. The analysis of amino acid residues using Weblogo showed that most of the dominant epitopes were conserved, with only a few amino acid residues experiencing mutations, as shown in Fig. 2. In the analysis of epitope hiding, it was found that all the selected epitopes except B-6 epitopes were located externally (The complete extracellular region of S protein is shown in supplementary figure S8).

Table 1 TGEV S protein B-cell, CTL, and Th dominant epitope screening results.
Fig. 2
figure 2

Weblogo analysis of dominant epitope amino acid residue conservation results (Polar: green; Neutral: purple; Basic: blue; Acidic: red; Hydrophobic: black).

TGE S protein vaccine construction

Based on the above analytical screening, the candidate nanoparticle vaccine was designed by selecting the above T-cell and B-cell potential epitopes for vaccine construction. Figure 3 illustrates the selection of 12 dominant epitopes. Th epitopes were connected using GPGPG flexlinker, CTL epitopes using AAY flexlinker, and B-cell epitopes using KK flexlinker. Additionally, the ferritin N-terminus was tandemly linked with the C-terminus of the dominant epitope B-6 using EAAAK flexlinker to form the complete candidate nanoparticle vaccine (The amino acid sequence of the candidate vaccine is shown in supplemental Table S4).

Fig. 3
figure 3

TGEV S protein multi-epitope vaccine kit.

Assessment of candidate nanoparticle vaccine

The antigenicity evaluation of the candidate nanoparticle vaccine yielded a score of 0.752474 using the ANTIGENpro server. Solubility analysis during overexpression indicated a score of 0.647714 via the SOLpro server. Additionally, the AllerTop server confirmed the non-sensitizing nature of the candidate nanoparticle vaccine, indicating its antigenicity, solubility, and non-sensitizing properties. ProtScale server analysis of hydrophobicity revealed a large number of negative regions for candidate nanoparticle vaccines, and the results are shown in Fig. 4A. The TMHMM-2.0 server analyzed the transmembrane structural domains in the candidate nanoparticle vaccine were all extracellular regions, and the results are shown in Fig. 4B, indicating that the candidate nanoparticle vaccine has hydrophilic and extracellular secretion functions. Further characterization of physical and chemical properties via the ExPASy server revealed essential parameters: the candidate nanoparticle vaccine comprised 368 amino acids, with a molecular weight of 41814.09 and a theoretical pI of 6.69. The instability index (II) computed as 28.75 classifies the protein as stable—aliphatic index as 91.47 and the Grand average of hydropathicity (GRAVY) as -0.153. The secondary structure analysis, performed using the SOPMA and PSIPRED servers, depicted a composition of Helix, sheet, Turn, and Coil, with PSIPRED indicating percentages of Alpha helix (43.75%), Extended strand (22.01%), Beta turn (8.42%), and Random coil (25.82%), as shown in Fig. 4C,D, respectively.

Fig. 4
figure 4

Results of candidate nanoparticle vaccine evaluation. (A) candidate nanoparticle vaccine hydrophobicity analysis; (B) candidate nanoparticle vaccine transmembrane domain analysis; (C) candidate vaccine secondary structure analyzed by SOPMA server; (D) candidate vaccine secondary structure analyzed by PSIPRED server.

Modeling, refinement and evaluation of tertiary structures of candidate nanoparticle vaccine

The tertiary structure of the candidate nanoparticle vaccine was generated using the trRosetta server via homology modeling, resulting in a model with a TM-score of 0.632. Visualization of the structure was performed using Pymol software, as depicted in Fig. 5F. Additionally, the trRosetta server generated a 2D view of the tertiary structure, and the results are shown in Fig. 5A–E. To ensure that the protein resembled the natural structure, the original tertiary structure model of the generated candidate vaccine was later refined using the Galaxy-Refine server (The detailed scores of the tertiary structure refinement model are shown in supplementary table S5). Structural comparisons were made using Pymol software, with results illustrated in Fig. 6B. The level of refinement pre- and post model refinement was assessed by generating Ramachandran plots of the protein’s tertiary structure using the PROCHECK server. In the crude model, 88.4% of residues fell within the most favoured regions, 10.0% in additional allowed regions, 0.6% in generously allowed regions, and 0.9% in disallowed regions, as depicted in Fig. 6A. Following refinement, the percentages shifted to 93.6%, 4.9%, 0.9%, and 0.6% for residues in the respective regions, as illustrated in Fig. 6C (The detailed comparison results before and after the tertiary structure refinement of candidate vaccines are shown in supplementary Figure S1-S5). The results demonstrate that the model quality was improved after refinement. However, further stabilization of the protein tertiary structure through disulfide engineering was deemed necessary for the refined model. This engineering was executed via the DbD2 server. Notably, the refined model displayed a B-factor of 0 in the analysis results, indicating that no sequence mutations were required to introduce disulfide bonds, as shown in Fig. 6D. The ProSA web server performed the final refined model scoring and validation, and the results are shown in Fig. 6E,F. The model exhibited a Z-Score of -4.83, followed by predominantly negative local energy scores. These findings affirm that the refined model conforms to natural protein conformation, rendering it an optimal candidate for subsequent analytical steps.

Fig. 5
figure 5

Modeling of the tertiary structure of the candidate nanoparticle vaccine. (A)–(E): Are 2D structural parameters (including Omega, Phi, Distance, Theta, and Contact); (F) The tertiary structure of the candidate nanoparticle vaccine as visualized by the Pymol software (Th-dominant epitopes are in blue, CTL-dominant epitopes are in purple, B-cell-dominant epitopes are in yellow, and ferritin is in red).

Fig. 6
figure 6

Refinement of the tertiary structure of the candidate nanoparticle vaccine. (A) Ramachandran graphic validation of the tertiary structure; (B) cartoon comparison of the before and after model structure (green for the coarse model, red for the refined model); (C) Ramachandran graphic validation of the refined tertiary structure; (D) disulfide engineering of the refined model (the yellow part is the disulfide bond); (E) ProSA SEB server Evaluating the refined model; (F) ProSA SEB server evaluating the local model quality.

Analysis of structural B cell epitopes

Ellipro server identified a total of 6 Linear epitopes and 6 discontinuous epitopes within the 368 amino acids of the candidate nanoparticle vaccine. As presented in Table 2, all predicted scores surpassed the server’s threshold of 0.5.

Table 2 Prediction results of candidate vaccine epitopes by Ellipro server.

Docking of candidate nanoparticle vaccine molecules

The ability of candidate nanoparticle vaccines to interact with TLRs is an important indication of defense against pathogens. To evaluate this, protein–protein docking was conducted using the pyDockWEB server. The server comprehensively ranked the complexes, with the top-ranked docked complexes selected for subsequent analysis. From the docking analysis, it is obvious that the candidate nanoparticle vaccine can be stably docked with TLR-3 and TLR-4. The composite scores (kcal/mol) for docking with TLR-3 were as follows: Electrostatics: -10.549, Desolvation: -47.666, Vdw: 86.453, with a Total score of -49.569. Pymol software was utilized to visualize the local interaction sites, as depicted in Fig. 7A–C, while LigPlot + software generated the 2D interaction interface, illustrated in Fig. 7D, There were 6 amino acid residues with hydrogen bond and 23 amino acid residues with hydrophobic interaction. The composite scores (kcal/mol) for docking with TLR-4 were as follows: Electrostatics: -9.895, Desolvation: -47.614, Vdw: 63.102, resulting in a Total score of -51.199. Pymol software was used to visualize the local interaction sites, depicted in Fig. 7E–G, while LigPlot + software generated the 2D interaction interface, illustrated in Fig. 7H, There were 6 amino acid residues with hydrogen bond and 23 amino acid residues with hydrophobic interaction. Overall, the docking analysis revealed total energy consumption scores around − 50 kcal/mol, with the two-dimensional interaction interface displaying numerous interacting amino acid residues, indicating a robust interaction between the candidate vaccine and TLR-3 and TLR-4.

Fig. 7
figure 7

Results of molecular docking of candidate nanoparticle vaccine with TLR-3 and TLR-4. (A) and (E): Docking of candidate nanoparticle vaccine with TLR-3 and TLR-4 (candidate vaccine structure is in red, TLR-3 and TLR-4 are in green); (B) and F: Magnification of local interaction position; (C) and (G) Specific interaction amino acid residues; (D) and (H): Two-dimensional interaction interface.

Molecular dynamics simulation analysis

In MD simulation, Fig. 8A–G shows the results of docking complex analysis between TLR-3 and candidate vaccine, and Fig. 8H–N shows the results of docking complex analysis between TLR-4 and candidate vaccine. It can be seen from the figure that RMSD curves remain between 4–5 nm and 12–14 nm respectively after 45 ns, indicating that the complex tends to be stable during the simulation process and the structure is relatively stable. The RMSF curve Chain A showed different behaviors up to 1.5 and 4.0 nm, respectively. Most of the regions in TLR-3 showed relative stability, while some regions in TLR-4 showed higher fluctuations, which proved that the dynamic behavior of amino acid residues was more stable in the TLR-3 complex. Rg curve analysis further demonstrated the spatial conformation of the complex, and both complexes fluctuated within 100 ns. In the analysis of hydrogen bond curves, about 10 and 20 hydrogen bonds were found in 20 to 60 ns, respectively, showing strong interactions. SASA analysis showed a decline from 560 to 500 nm2 and 810 nm2 to 750 nm2 within the first 40 ns, respectively, and then remained stable, reflecting the tight binding and stability of the complex. In Gibbs free energy analysis, the 3D free energy landscape map is drawn, and the RMSD is in the range of 0.0–2.5. The Gyrate fell to 0.0–4.8 and 0.0–5.8 respectively, which further confirmed the relatively good stability of the two docking complexes.

Fig. 8
figure 8

Molecular dynamics simulation results, TLR-3 docking complex (A)–(G), TLR-4 docking complex (H)–(N).

Immune simulation of candidate nanoparticle vaccine

The immune simulation reaction was performed on the C-ImmSimm server, and the results are shown in Fig. 9. The candidate nanoparticle vaccine was injected three times and significantly induced primary immune responses after three injections. And three peak curves were formed, which proved that immune memory was successfully created and memory T and B lymphocytes produced a rapid response. Notably, secondary immune responses were stimulated, with gradual increments observed in primary immune responses after each subsequent injection. A marked elevation in secondary and tertiary antibody levels (IgG1 + IgG2, IgM and IgG + IgM) was observed. Furthermore, there was a continual increase in the concentration of active B cells, plasma B cells, helper T cells, and cytotoxic T cells, peaking during antigen injection, indicating effective immune responses, establishment of immune memory, and efficient antigen clearance. Increases in dendritic cells and macrophages pointed to excellent antigen presentation by antigen-presenting cells (APCs), while the candidate nanoparticle vaccine resulted in significant increases in cytokine (IFN-γ, IL-4, and IL-10) levels. Therefore, the simulated immune response demonstrated that the candidate nanoparticle vaccine elicited substantial production of immunoglobulins, APCs, cytokines, and activated B and T cells, showcasing its ability to induce potent immunogenic responses post-vaccination within the host.

Fig. 9
figure 9

Results of immune response induced by three doses of the vaccine. (A) Overall change trend of B lymphocytes; (B) Plasma B lymphocytes count; (C) B lymphocytes population per entity-state; (D) and (E) Helper t lymphocyte population; (F) Regulatory T cell population; (G) Cytotoxic T lymphocyte population; (H) Cytotoxic T-cell population per state; (I) Natural killer cell population; (J) Dendritic cell population per state; (K) Macrophages population per state; (L) EP cell population; (M) The change trend of antigen and immunoglobulin; (N) Trends of cytokines and interleukin.

Codon optimization and in-silico cloning of vaccine

The efficient expression of candidate nanoparticle vaccines is a key step in silic cloning. Therefore, the amino acid sequences of the candidate nanoparticle vaccine were codon-optimized for prokaryotic organisms to improve their translational efficiency based on E. coli K12 in the JCat server. The full length of the nucleotide sequence of the candidate nanoparticle vaccine was 1104 bp after reverse translation (The optimization results of Escherichia coli favoritism are shown in supplementary Figure S6, The complete nucleotide sequence is shown in supplemental table S6). The CAI-Value of the improved sequence reached 1.0, while the GC-Content of E coli was at 47.92%. Upon validation using the ExpOptimizer server, the CAI was determined to be 0.74, with a corresponding GC content of 47.92%. These results, presented in Fig. 10A,B, indicated that all CAI values exceeded the threshold of 0.5, ensuring abundance in the DNA sequences. Additionally, the GC content fell within the optimal scoring range of 30 to 70%. To facilitate expression and purification, XhoI and BamHI restriction endonuclease sites were selected using SnapGene software, and the adapted nucleotide sequences were inserted into the pET28-a(+) expression vector with a 6× his tag sequence. The prokaryotic expression plasmids were constructed, and the results are shown in Fig. 10C. To verify successful insertion, the constructed plasmids underwent double enzyme digestion and PCR simulation, followed by analysis using 1% agarose gel electrophoresis. The results, shown in Fig. 10D,E, confirmed the successful insertion of the candidate vaccine, with bands observed at the correct positions (See supplementary Figure S7 for the complete agarose gel electrophoresis simulation).

Fig. 10
figure 10

Construction of the prokaryotic expression system of the candidate nanoparticle vaccine. (A) ExpOptimizer server to verify the CAI value; (B) ExpOptimizer server to verify the GC content; (C) Candidate nanoparticle vaccine (red) inserted into pET-28a(+); (D) Double enzymatic digestion of pET-28a-Vaccine plasmid (Cropped from supplemental Figure S7); (E) PCR simulation of the candidate nanoparticle vaccine Analysis (Cropped from supplemental Figure S7).

Discussion

In recent years, coronaviruses have caused serious public health problems with an alarming frequency of recombination and mutation, posing a great threat to human and animal health42. Porcine enteric coronaviruses have led to more frequent infections in swine herds and have caused significant losses to the global farming industry due to their potential for cross-species transmission1. Transmissible gastroenteritis virus causes acute, high-contact enteric infections in neonatal piglets, which are often clinically mixed with other enteroviruses or secondary infections, posing a serious challenge to disease prevention and control in swine farms3. Vaccination is known as a highly effective public health measure for disease prevention, aiding in immune response stimulation and pathogen containment, yet the scarcity of approved and efficacious viral vaccines remains a problem4. At present, many achievements in basic virology, vaccine technology, immunoinformatics and synthetic biology have opened up new opportunities for vaccine development against TGEV43. Strategies for the development of new vaccines as appropriate solutions for the prevention and treatment of viral infections are a primary need and will help to find solutions to the current biosecurity problems on swine farms.

However, the preparation and development of new vaccines is a complex and systematic process constrained by uncertainties related to safety, efficacy, and production feasibility, which limits the advancement of traditional vaccinology and the widespread use of existing vaccines44. One promising method to overcome these challenges is the development of subunit vaccines using immunoinformatics, which relies on advanced computational prediction of protein information45. However, discrepancies between predicted and actual results, along with the intricate and dynamic nature of the animal immune system, can lead to failures in later-stage experiments46. An effective strategy to mitigate these issues is to perform multiple server comparison screenings for the same sequence and use different alleles of animals for a comprehensive evaluation47. Since the epitope of TGEV S protein has not been fully resolved, predictive screening is a necessary method to rapidly enrich its epitope landscape and discard redundant sequence interference. In this study, immunoinformatics-related servers and software were employed for repeated comparison, and screening of TGEV S proteins, and the dominant epitopes identified were combined with a ferritin self-assembled nanoparticle delivery platform to construct candidate vaccines. This tandem polyvalent presentation of dominant epitopes enhances the activation of antigen-presenting cells (APCs) and the maturation of dendritic cells (DCs), thereby amplifying B cell signaling and effectively triggering both humoral and cellular immune responses, It will be beneficial for candidate nanoparticle vaccines to produce protection against TGEV5. Ultimately boosting the host’s protective immunity, the relatively conserved dominant epitope reduces the antigen sequence, alleviating excessive antigen load and allergic reactions of the host, without reversing viral pathogenesis and effectively combating viral genetic variation and antigen drift. The dominant Th-1 and CTL-1 epitopes in the study overlap with the reported neutralizing antibody sites at 538–591 aa and 378–395 aa at site D, providing robust support for the vaccine’s practical efficacy in virus neutralization48. Additionally, to ensure enhanced immune response and higher affinity, the candidate nanoparticle vaccine underwent systematic evaluation (antigenicity, allergenicity, immunogenicity, hydrophobicity, solubility, and toxicity). Most of the secondary structure converted to antigen mosaicism, forming a solid structural basis, with B-cell epitope predictions exceeding server thresholds, indicating a potential for inducing antibody formation49. The docking of TLR-3 and TLR-4 molecules revealed numerous interacting amino acid residues and low binding free energy, MD simulation shows that the docking complex and its conformation are relatively stable at the microscopic level, elucidating the molecular mechanism of antiviral activity and suggesting that the vaccine may induce transcriptional up-regulation of type I interferon, pro-inflammatory cytokines, and chemokines, thereby promoting the host’s antiviral response50. In the immune simulations, Although C-ImmSimm server cannot simulate the ferritin self-assembly process, the secondary and tertiary immune responses following vaccination with candidate nanovaccines exhibited higher levels compared to primary immune responses. This resulted in elevated immunoglobulin production, IFN-γ secretion, and sustained cellular and CTL responses. These findings align with previous observations, such as the results seen in the simulation of COVID-19 S protein47, the research results of tuberculosis vaccine MP3RT in animal experiments51, and the self-assembly results of ferritin and IAV main proteins can trigger a powerful immune response52. Moreover, to address the challenges associated with protein expression purification for large-scale industrial vaccine production, the sequence was incorporated into the pET28-a(+) expression vector with His tags added at both ends. However, despite the systematic design and evaluation of the constructed nanoparticle vaccine candidate, the study’s limitations necessitate future focus on in vitro and in vivo testing to ascertain its practical application.

Conclusions

Multi-epitope vaccines and nanoparticle vaccines have received much attention and progressed to clinical trials. In this study, we successfully combined the two using immunoinformatics leveraging a comprehensive analysis of the S protein’s full-length gene sequence to identify immunologically superior B-cell and T-cell epitopes tandemly linked with ferritin. The resulting candidate nanoparticle vaccine exhibits high antigenicity, along with non-allergenic and non-toxic properties. Through stable binding to immune receptors, TLR-3 and TLR-4, it triggers strong cellular and humoral immunity in the body as shown in the immune simulation and efficient expression in computerized clones. We believe that the candidate nanoparticle vaccine could be a potent vaccine against the TGEV pathogen and, at the same time, provide a new approach and an efficient research platform for accelerating the creation of new vaccines against porcine enteric coronaviruses.