A Computational Reverse Vaccinology Approach for the Design and Development of Multi-Epitopic Vaccine Against Avian Pathogen Mycoplasma gallisepticum

Introduction

Mycoplasma gallisepticum is cell-wall-less bacteria, the key avian respiratory pathogen which is the causative agent for chronic respiratory disease (CRD) in chickens and infectious sinusitis in turkeys and house finches (1). M. gallisepticum infection tends to intensify the pathology of accompanying viral infections like Newcastle disease virus (NDV), infectious bronchitis (IB), and CRD in chickens. M. gallisepticum is the crucial pathogen in poultry industries leading to significant economic loss by infecting layers and breeder and broiler poultry flocks (2). The reason for CRD outbreak and spread includes heat/cold stress, excessive ammonia, and unhygienic conditions (3, 4). The affected chickens manifest clinical symptoms like respiratory distress, nasal and ocular discharge loss of appetite, and decreased egg production.

M. gallisepticum are classified under the class Mollicutes and family Mycoplasmataceae. M. gallisepticum has a small genome with 996 kb responsible for its limited biosynthetic capabilities and low replication rates and tends to infect the respiratory tract and colonize in the mucosal surface tract of chicken. Transmission of M. gallisepticum in respiratory tract mucosa induces local inflammation in the trachea, lungs, and air sacs and sporadically in the conjunctiva and oviduct. The most crucial step for initial infection is the adhesion of M. gallisepticum to the host cell surface which is achieved by cytoadherence proteins GapA and cytoadherence-related molecule A (CrmA) (5). Previous reports suggest that the absence of GapA and CrmA in mutant strains of M. gallisepticum does not cause infection, proving the importance of these two cytoadherence proteins in infection. Further, the fibronectin-binding proteins M. pneumoniae-like protein A (PlpA) and HMW3-like protein (Hlp3) are known as cytoadherence accessory proteins because they are also involved in cytoadherence (6–8); these two protein-encoding genes are absent in avirulent strain M. gallisepticum R_high (9). After initial interaction with host cells, the variable lipoprotein hemagglutinin A (vlhA) causes extensive phase variation, which makes it possible for M. gallisepticum to escape from host immune response (10, 11). The vlhA gene family consists of 43 genes in R low strain, which plays a prominent role in M. gallisepticum pathogenesis. Among these, vlhA 1.07, vlhA 4.01, vlhA 2.02, vlhA 3.03, and vlhA 5.13 were predominantly expressed through the course of infection (12, 13). For the development of recombinant vaccine against M. gallisepticum infection in chickens, these proteins could be the effective target.

M. gallisepticum infection is controlled by antibiotics, which is a regular practice in poultry farms; antibiotics like fluoroquinolones (enrofloxacin, difloxacin), pleuromutilins (tiamulin), tetracyclines, and macrolides (tylosin, tilmicosin) were widely used (14, 15). The extensive usage of antibiotics has paved a way in the emergence and increase of antibiotic resistance bacteria (16). One such example is quinolone resistance in M. gallisepticum which happens to be frequent, dampening the treatment and control methods (17). The antibiotics can reduce the population of M. gallisepticum in the respiratory tract, thereby reducing the risk of spread to neighboring flocks (18). Currently available vaccines for M. gallisepticum infections are killed bacterins, live attenuated strains, and recombinant protein; of these, the commercially approved live attenuated vaccine includes 6/85 strain, ts-11, and F strains which are widely used in layers and commercial flocks. Live attenuated vaccines often display pathogenicity as well as side effects, while bacterins are extremely costly and repeated doses are needed to improve the avian immune system. Additionally, vaccination stress and post-vaccination reaction can also affect chickens (19).

Thus, novel recombinant strategies are required for developing more efficient and less expensive vaccines. So far, no in silico approaches had been reported for the development of vaccine; however, using in silico approaches few virulence genes have been identified for mycoplasma species (8, 20). Moreover, in silico-based multi-epitopic vaccines for cancer were tested in animals and showed promising results (21–23). Thus, there is scope for novel in silico-based multi-epitopic vaccines to be successful for infectious diseases. Therefore, a successful vaccine is essential to addressing the current state of Mycoplasma infections in chickens.

Immunoinformatics is a rapidly developing field in bioinformatics, which uses up-to-date technologies to predict potential peptides for vaccine development against various diseases. Recent advances in immunoinformatics allow us to exploit the available data to predict the most efficient epitopic regions in antigenic proteins, which further led to epitope-based vaccine development. The B-cells and T-cells play a central role in the activation of immune response against various viral and bacterial infections. Immunization based on epitopic vaccine has high potential in stimulation of humoral and cellular immune response. These vaccines are composed of highly immunogenic B and T cell epitopes which intensify the immune response through B-cell and T-cell activation. The setback of epitopic vaccines is as follows: due to its smaller size, it may not be recognized as an immunogenic substance by the host immune system, but it can be overcome by the usage of multiple epitopes and adjuvant. The prospective advantages of epitope-based peptide vaccine are as follows: it decreases the risk associated with other types of vaccinations; the epitopes can be wisely engineered and optimized to increase its potency in eliciting stronger immune response; and it has chemical stability due to its smaller size. In addition to that, epitope-based vaccines also offer cost effectiveness, safety, and selection-based immunity.

Hence, in this study we employed immunoinformatics approaches for the development of epitope-based peptide vaccine from cytoadherence proteins, fibronectin-binding proteins, and variable lipoprotein hemagglutinin A. The selected proteins have a quintessential role in infection and pathogenesis of M. gallisepticum into the host cells based upon earlier studies (5, 6, 9, 10, 12, 13). Further, we used reliable computational approaches to predict the highest T-cell as well as B-cell epitopes for these proteins. A multi-epitope vaccine was constructed, validated, and docked against chicken TLR-2 and TLR-5. The pipeline of the computational methods used for vaccine development is depicted in Figure 1. Eventually, this study is aimed at assisting the future laboratory efforts in developing potential vaccine for M. gallisepticum. The designed multi-epitopic vaccine could act as a legitimate prophylactic therapy for M. gallisepticum and its associated illness in chickens.

FIGURE 1

Figure 1. Illustration of steps involved in the development of multi-epitopic vaccine candidate from antigenic protein of M. gallisepticum.

Methodology

Dataset Construction and Initial Analysis

A dataset was constructed with four cytoadherence proteins and five vlhA proteins of M. gallisepticum strain R low from UniProt Knowledgebase. The dataset was further improved by incorporating the protein name and UniProt accession number. Further initial antigenicity and allergenicity analyses for all the proteins were carried out using ANTIGENpro (http://scratch.proteomics.ics.uci.edu/) (24), Vaxign-ML (25), and AllerTOP v. 2.0 (26), respectively.

Prediction of T Cell Epitope

The lack of chicken MHC alleles in the immunoinformatics database was a major drawback; thus, human MHC alleles were substituted for chicken alleles to predict the epitopes. Previous studies used human alleles to predict epitopes against poultry pathogens (27, 28). Reports have revealed that chicken B–F alleles have the capacity to induce immune response comparable with human class I alleles predominantly during antigen presentation (29, 30). For this reason, the most appropriate human alleles for chicken were chosen from a previous study (31). Alleles HLA^*B 40:06, HLA^*B 41:04, and HLA^*B 41:03 for MHC class I and alleles DRB1:1482, DRB1:1366, DRB1:1310, and DRB1:1445 for MHC class II were chosen.

MHC I Binding Epitope Prediction

To predict MHC I binding affinity, the NetMHCcons 1.1 server was used (http://www.cbs.dtu.dk/services/NetMHCcons) (32); the epitopes were screened based on threshold IC50 values: IC50 <50 nM indicates strong binders, whereas IC50 > 500 nM indicates weak binders. The percentage rank for strong binding epitopes was set as 0.5, and that for weak binding epitopes was 2. The length of the epitopes was set as 9 mers. The class I immunogenicity scores for the selected peptides were calculated by the IEDB Analysis Resource (http://tools.iedb.org/immunogenicity) (33).

MHC II Binding Epitope Prediction

MHC II binding epitopes were predicted using the NetMHCIIpan 3.1 Server (http://www.cbs.dtu.dk/services/NetMHCIIpan-3.1/) (34). Among the available options, the peptide length was set to 15 amino acids. A number of peptides binding to MHC alleles were then predicted by the server, from which peptides were selected based on percentile rank: the percentage rank for strong binding epitopes was set as 0.5, and that for weak binding epitopes was 2.

B-Cell Epitope Prediction

The potential linear B cell epitopes were predicted by BepiPred (http://tools.iedb.org/bcell/) from the IEDB research resource (35). BepiPred predicts continuous epitopes by combining two residue properties with the hidden Markov model.

Construction of Multi-Epitopic Vaccine

A multi-epitopic vaccine against M. gallisepticum was designed from CTL, HTL, and B cell epitopes from the antigenic proteins. The adjuvant improves the immunogenicity of the vaccine candidate; here, we used the avian β defensin 8, chicken TLR-5 agonist in the N-terminal of the multi-epitopic vaccine construct followed by EAAAK linker and CTL epitopes. The intra-CTL, intra-HTL, and intra-B-cell epitopes were linked using AAY, GPGPG, and KK linkers, respectively.

Determination of Antigenicity, Allergenicity, and Solubility Profile

ANTIGENpro (http://scratch.proteomics.ics.uci.edu/) (24) predicted the antigenicity of the multi-epitope vaccine, whereas AllerTOP v2.0 (http://www.ddg-pharmfac.net/AllerTOP/) (26) is employed to analyze the allergic nature of the designed vaccine construct. The SOLpro online tool was used for the calculation of the solubility Escherichia coli (E. coli) host (http://scratch.proteomics.ics.uci.edu/) (36).

Determination of Physiochemical Properties

The ExPASy-ProtParam server (http://web.expasy.org/protparam/) was employed to calculate the physicochemical properties of the designed multi-epitopic vaccine. This tool computes several parameters like molecular weight, isoelectric point (pI), half-life, atomic composition, aliphatic index, instability index, and grand average hydropathicity (GRAVY) (37).

Secondary Structure Analysis

A secondary structure analysis of the multi-epitope vaccine construct was performed using the Garnier–Osguthorpe–Robson (GOR IV) online server with mean accuracy of 64.4% (https://npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl?page=npsa_gor4.html) (38) and position specific iterated prediction (PSIPRED) analysis on outputs from PSI-BLAST (http://bioinf.cs.ucl.ac.uk/psipred/) (39).

Tertiary Structure Prediction, Refinement, and Validation

The I-TASSER server was used in tertiary-structure modeling of the vaccine construct (https://zhanglab.ccmb.med.umich.edu/I-TASSER/) (40). The tertiary structure modeled by I-TASSER was refined by the GalaxyRefine server (http://galaxy.seoklab.org/cgi-bin/submit.cgi?type=REFINE) (41). The galaxy server applies molecular dynamic simulation for overall structural relaxation and swaps amino acids with high-probability rotamers to refine the structures. The output generally contains five refined models, with GDT-HA, RMSD, MolProbity, Clash score, Poor rotamers, and Rama favored scores (42). The output was authenticated by PROCHECK (https://servicesn.mbi.ucla.edu/PROCHECK/), which verifies the stereochemical quality of refined protein structures by examining residue by residue and the overall structural geometry (43).

Immunological Analysis

The presence of linear (continuous) and conformational (discontinuous) B cell epitopes in the multi-epitopic vaccine construct was analyzed. The BcePred online server (https://webs.iiitd.edu.in/raghava/bcepred/index.html) was employed in the prediction of continuous B-cell epitopes (44). ElliPro from IEDB (http://tools.iedb.org/ellipro/) was used in the prediction of discontinuous B-cell epitopes, which predicts epitopes based on a protein's 3D structure (45).

Multi-Epitope Vaccine Protein Disulfide Engineering

Disulfide bonds offer considerable stability and strengthen the geometric conformation of proteins. The online DbD2 server, available at http://cptweb.cpt.wayne.edu/DbD2/index.php, was used in this purpose. The pair of residues able to form a disulfide bond, if every amino acid residue is mutated to cysteine, is identified by this web server (46).

Molecular Docking

In this study, interactions between the vaccine construct with chicken TLR-2 and TLR-5 receptors was studied as they are bacterial sensing TLR (47). The active and passive residues in both the receptor and vaccine construct were predicted by CPORT (48). The HADDOCK 2.4 (https://www.bonvinlab.org/software/haddock2.4/) online server performed docking between multi-epitopic vaccine construct and chicken TLR-2 and TLR-5 receptors (49). HADDOCK Refinement Interface was used to refine the best cluster which has the lowest HADDOCK score. The best structure after refinement from the docked complex was chosen, and their binding affinity was calculated using the PRODIGY web server (50, 51). Lastly, PDBsum (https://www.ebi.ac.uk/thorntonsrv/databases/pdbsum/Generate.html) mapped the interacting residues between the docked complexes (52).

Codon Adaptation and in silico Cloning of a Multi-Epitopic Vaccine Construct

The JAVA Codon Adaptation Tool (JCAT) server (http://www.jcat.de) accomplished codon optimization as per E. coli codon usage (53). Certain factors that are important to the efficacy of gene expression including codon adaptation index (CAI) and GC content adjustment were enhanced. Cloning and expression of the multi-epitopic vaccine construct in a suitable expression vector were achieved using SnapGene in silico cloning.

In silico Immune Simulation

In silico immune simulation was performed by the C-ImmSim server (https://kraken.iac.rm.cnr.it/C-IMMSIM/) (54) to analyze the immunogenicity and immune response of the designed multi-epitopic vaccine. The entire simulation ran for 1,100 time steps, where one time step is 8 h. Three peptide injections were given at day 1, day 84, and day 168, and other parameters were kept default.

Results

Dataset Construction and Initial Analysis

The dataset for this study was constructed by including the protein name, UniProt accession number, and antigenicity and allergenicity of the selected proteins (Table 1). The threshold score for antigenicity is 0.4; all the proteins obtained scores higher than the threshold score and thus has a potential to induce immune response. All the selected protein sequences are non-allergens.

TABLE 1

Table 1. Dataset and antigenicity and allergenicity prediction.

T-Cell Epitope Prediction

Highly ranked MHC-I-binding peptides were selected based on IC50; IC50 <50 nM were identified as strong binders by using the NetMHCcons 1.1 server (Supplementary Table S1). MHC-II binding peptides were predicted using the NetMHCIIpan 3.1 server, and the epitopes with %rank <2 were regarded as strong binders (Supplementary Table S2) and preferred for further studies.

B-Cell Epitope Prediction

The BepiPred server predicted the B-cell binding epitopes with the threshold score above 0.350 and considered them as probable epitopes. The epitopes from the antigenic proteins of M. gallisepticum with length ranging from 10 to 15 peptides are given in Supplementary Table S3.

Multi-Epitopic Vaccine Design

The selected MHC-I and MHC-II-binding epitopes and B-cell epitopes derived from M. gallisepticum GapA, plpA, Hlp3, CrmA, vlhA 1.07, vlhA 2.02, vlhA 3.03, vlhA 4.01, and vlhA 5.13 were used to design the multi-epitope vaccine construct. Avian β-defensin 8 sequences were added as an adjuvant followed by the epitopes and linked by specific linkers (Figure 2).

FIGURE 2

Figure 2. (A) Three-dimensional structure of the engineered vaccine construct. (B) Graphical representation of the final multi-epitopic vaccine construct along with CTL, HTL, and B-cell epitopes.

Antigenicity and Allergenicity Prediction of Multi-Epitopic Vaccine

The foremost decisive factor needs to be guaranteed while designing a vaccine, which is the antigenic potential of the vaccine candidate to induce a humoral and/or cell-mediated immune response against the targeted virus and allergenicity of the constructed vaccine. The antigenicity score for multi-epitope vaccine was 0.683753 as predicted by the ANTIGENpro server. The final vaccine construct was predicted to be non-allergenic by AllerTOP v2.0. The solubility prediction using the SOLpro server illustrated that the vaccine construct is soluble with a probability of 0.749662.

Physiochemical Characterization of the Designed Vaccine

The constructed multi-epitope vaccine contains 498 amino acids and a molecular weight of ~54.15 kDa. The estimated pI was 9.75, which points out the alkaline nature of the multi-epitopic vaccine. The estimated instability index, aliphatic index, and grand average of hydropathicity (GRAVY) were 26.09, 79.64, and −0.386, respectively, suggesting that the vaccine is stable, thermostable, and hydrophilic in nature. The total number of positively charged residues (Asp+ Glu) was predicted to be 65 while the total number of negatively charged residues (Arg+ Lys) was predicted to be 33. The extinction coefficient at 280 nm in H₂O was 47,595. The estimated half-lives in mammalian reticulocytes (in vitro), yeast cells (in vivo), and E. coli (in vivo) were 30 h, <20 h, and <10 h, respectively.

Secondary and Tertiary Structure Analysis

The secondary structure of the multi-epitopic vaccine consists of 48.80% (243 aa) of random coils, 34.74% (173 aa) alpha-helices, and 16.47% (82 aa) betasheets. Figure 3 represents the graphical image of the GOR IV server and PSIPRED server. The tertiary structure of the vaccine was modeled using the I-TASSER server. The I-TASSER server modeled five tertiary structures of the designed multi-epitopic vaccine protein based on 10 threading templates. All the 10 preferred templates showed high-quality alignment as per their Z-score values (ranging from 0.38 to 2.40). The predicted model has a C-score value of −2.31. The preferred model had an estimated TM score of 0.44 ± 0.14 and RMSD of 12.9 ± 4.2Å. The TM score indicates the structural similarity between two structures. A TM score >0.5 specifies a model of correct topology, and a TM score <0.17 means random similarity.

FIGURE 3

Figure 3. Predicted secondary structure of the multi-epitope vaccine construct. (A) Graphical illustration of the secondary structure by the PSIPRED server. (B) GOR IV secondary structure prediction-alpha helix (blue), random coils (yellow), and extended strand (red).

Multi-Epitopic Vaccine Structure Refinement and Validation

The 3D structure modeled by the I-TASSER server was validated by the PROCHECK server and Ramachandran plot analysis was performed, and the output revealed that 52.6% of residues fall in the favored region. For further refinement, the GalaxyRefine server (http://galaxy.seoklab.org/cgi-bin/submit.cgi?type=REFINE) was used; the output provided five refined models for the multi-epitopic vaccine construct (Supplementary Table S4). Among them, model four was considered to be the best-refined structure with parameters such as GDT-HA of 0.8886, RMSD of 0.576, MolProbity of 2.565, Clash score of 19.2, Poor rotamers of 1.2, and favored Rama of 81.7%. Thus, this model (Figure 4) was finally selected for further investigations.

FIGURE 4

Figure 4. (A) Tertiary structure model of the multi-epitope vaccine construct and (B) Ramachandran plot validation.

Immunological Property Assessment

The continuous B-cell epitopes in multi-epitope vaccine was predicted by BcePred with default parameters. Vital properties of epitopes like hydrophilicity, antigenicity, flexibility, accessibility, polarity, and exposed surfaces were predicted (Table 2). Conformational B cell epitopes in multi-epitope vaccine were predicted by ElliPro. In total, seven conformational epitopes were predicted; the details of the prediction are given in Supplementary Table S5 and Figure 5.

TABLE 2

Table 2. Results of linear B-cell epitopes in multi-epitope vaccine by the BcePred server.

FIGURE 5

Figure 5. The immunogenic conformational B-cell epitopes are represented as yellow globules. The details of (A–G) are given in Supplementary Table S5.

Vaccine Protein Disulfide Engineering

In total, 55 pairs of amino acid residues were predicted to form a disulfide bond using the DbD2 server. Just two residues, including GLU 145-SER 148 and ILE 152-ALA 156 were substituted by cysteine residues and found to fulfill the disulfide bond formation with residue assessment by chi3 and B-factor energy parameters (Figure 6). The residue screening was performed based on chi3 value (−87 to +97) and <2.5 energy value.

FIGURE 6

Figure 6. (A) Initial model without disulphide bonds. (B) Mutant model; the yellow stick represents the disulphide bond formation.

Molecular Docking

HADDOCK clustered 133 structures in 15 clusters, showing 66% of the models formed by water-refined HADDOCK for multi-epitopic vaccine-Chicken TLR-2 complex and 78 structures in 11 clusters, showing 39% of the models formed by water-refined HADDOCK for multi-epitopic vaccine-Chicken TLR-5 complex. The most precise predicted cluster is the cluster with the smallest HADDOCK ranking. The HADDOCK refinement server refined the top cluster where 20 structures were grouped into one cluster, resulting in 100% of the HADDOCK water-refined version for both the complexes. In Table 3, the results of the refined models were presented. A strong binding affinity between the multi-epitopic vaccine and chicken TLR-2 receptor and chicken TLR-5 receptor is indicated by the Haddock scores of −225.0 ± 4.2 and −230.9 ± 0.0, respectively; a negative score suggests better docking. The docked complex along with interacting residues between docked complexes is shown in Figure 7. A total of 19 hydrogen bonds, 6 salt bridges, and 160 non-bonded contacts were analyzed between the multi-epitopic vaccine and chicken TLR-2 receptor, and 13 hydrogen bonds, 3 salt bridges, and 150 non-bonded contacts were analyzed between the multi-epitopic vaccine and the chicken TLR-5 receptor. The Gibbs-free energy (ΔG) is a decisive parameter for deciding that the probability of an interaction can occur under specific conditions in the cell. The binding affinity of the docked complex was analyzed by the PRODIGY web server. The G-values were −14.4 and −14.6 kcal mol-1, and the dissociation constant was 2.9 E-11 and 2.0 E-11 for the vaccine construct and chicken TLR-2 and TLR-5, respectively. The results indicated by negative Gibbs-free energy (ΔG) shows that the docked complex is energetically viable.

TABLE 3

Table 3. Docking statistics of best-refined docked chicken TLR-2, TLR-5, and vaccine construct.

FIGURE 7

Figure 7. Molecular docking of multi-epitopic vaccine with chicken immune receptors. (A) Chicken TLR-2 (green) and multi-epitopic vaccine (red). (B) Chicken TLR-5 (cyan) and multi-epitopic vaccine (red).

Codon Adaptation and in silico Cloning

The CAI in the adapted sequence was 0.9713, and the GC content of the enhanced sequence was expected to be 51.74%, whereas the GC content of E. coli was 50.7. To perform in silico cloning, the vaccine construct sequence was examined for restriction enzyme sites; the Xho I and Bam HI restriction enzyme sites were not present in the vaccine sequence, so these enzymes were used for in silico cloning in the pET30a (+) vector. A clone of 6,882 bp was achieved subsequent to introduction of the vaccine construct into the pET30a (+) vector (Figure 8).

FIGURE 8

Figure 8. In silico restriction cloning of the gene sequence of the vaccine construct into the pET30a(+) expression vector; the red-colored region represents the gene coding the multi-epitopic vaccine, and the black-colored region points out the vector backbone.

In silico Immune Simulation

In silico immune simulation was performed to characterize the immune profile of the designed multi-epitopic vaccine construct. The immune response was extensively activated, and there is an increase in antibody titer after injection Figure 9A. The cytokine response was also analyzed; the results show reliable and vigorous response following injection (Figure 9B). The B and T cell populations were also increased considerably (Figures 9C,D). These results confirm that the designed multi-epitopic vaccine effectively elicits an immune response.

FIGURE 9

Figure 9. In silico immune simulation of multi-epitopic vaccine by the C-ImmSim server. (A) Antibody production in response to antigen injections (vertical black lines). (B) Expression of cytokine level. (C) Level of B-cell population. (D) Level of T-cell population.

Discussion

Vaccines are significantly accountable for providing host organisms with defense from a specific disease. The conventional vaccine development process is extremely meticulous, expensive, and time consuming. Immunoinformatics tools are a boon to design and develop vaccine in a short period of time with specificity and ease. It allows detecting pathogens underlying immunogenic proteins, together with prediction of diverse immune-dominant epitopes which are involved in the development of humoral and cell-mediated immune responses against the pathogen. For that reason, a multi-epitope-based peptide vaccine can be designed with immunogenic proteins of a pathogen. In the present scenario, a large number of peptide vaccines are in progress, majority of which are for human infectious diseases and tumor (55–58). Very few studies are reported in the field of in silico vaccine for poultry and animals. Previous reports focused on the prediction of epitopes for animal diseases like foot and mouth disease (59) and animal trypanosomiasis (60) in designing the multi-epitopic vaccine and offered effective immunity as evaluated with currently available vaccines. The following studies validate the immunoinformatics method to design multi-epitopic vaccines against infectious diseases in poultry (61–64). The advantages of this method include evaluation of complete antigenic epitopes and molecular modeling to study the probable binding with host proteins (65, 66). The in silico validations like molecular docking and in silico cloning and immune simulation (67–70) were included in this study to make the constructed multi-epitope as functionally suitable as a vaccine candidate against M. gallisepticum infection.

The purpose of multi-epitopic vaccine is to identify the specific CTL, HTL, and B-cell epitopes from the antigenic proteins that induce specific and effective immune response. Initial antigenicity and allergenicity analyses revealed that selected proteins are antigenic and non-allergic. Vaccine design requires an adjuvant that is necessary for increasing the effectiveness of the vaccine (71); hence, avian β-defensin was used as an adjuvant at the N-terminal followed by a sequence of various CTL, HTL, and B-cell epitopes present in antigenic proteins (72). The selected epitopes together with the adjuvant have been merged together with the help of suitable linkers. A good vaccine candidate is supposed to be competent in initiating an immune response without causing an allergic reaction; hence, the antigenicity and allergenicity of the multi-epitopic vaccine construct were analyzed. The vaccine construct is antigenic and non-allergic in nature. The molecular weight was 54.15 kDa, and the theoretical pI was 9.75, which signifies the basic nature of the vaccine construct. The aliphatic index and GRAVY score indicate that the vaccine is thermostable and hydrophilic. The secondary structure of the vaccine construct revealed that the random coil dominated the structure followed by alpha-helices and beta sheets. The tertiary structure was selected based on the highest c-score. The structure was further refined, and a better structure was obtained after refinement. The prediction of discontinuous and continuous B-cell epitopes in the vaccine construct shows that it can possibly interact with antibodies and is flexible. Further disulfide engineering was performed to stabilize the vaccine construct. The interaction prototype of the vaccine construct with chicken TLR-2 and TLR-5 receptors was analyzed by molecular docking studies. The results indicated a better interaction. The negative Gibbs-free energy (ΔG) indicates that the docked complex is energetically viable. Codon optimization was done for the E. coli K12 strain; the estimated GC content and CAI indicate enhanced the transcriptional and translational efficiency. Restriction enzymes Xho I and Bam HI were selected for cloning in the pET30a(+)vector. The expression of codon was adequate. In addition to this, in silico immune simulation results revealed that upon injection of the designed multi-epitopic vaccine, a strong immune response was triggered.

The currently available live attenuated vaccine and bacterins are commonly used in commercial birds; these vaccines cannot aid in control during the sudden onset of M. gallisepticum infection; strict biosecurity has to be followed to control and eradicate the infection. The live attenuated vaccine shows adverse side effects and pathogenicity, while the bacterins are associated with high cost and repeated dosage. The limitations of the currently available vaccine can be addressed by the present study which employs robust immunoinformatics tools to design multi-epitopic vaccine, and the result demonstrates that the designed multi-epitopic vaccine construct have every aspect be developed as an effective vaccine candidate, such as immune specificity, small size, and absence of adverse effects.

Conclusion

The current study listed out the applications of immunoinformatics to predict the epitopes of M. gallisepticum as a potential strategy to speed up vaccine development. M. gallisepticum infection has emerged as a severe problem in the poultry industry, resulting in huge economic loss. The present treatment strategies could not effectively eradicate the M. gallisepticum infection. Recently, there have been numerous attempts in developing vaccines against M. gallisepticum, which are mostly limited to live attenuated vaccines. This study predicted the potential CTL, HTL, and B-cell epitopes from the antigenic proteins against M. gallisepticum and designed a multi-epitope vaccine candidate. This approach of vaccine development will have advantages over the classical vaccine development methods by using small epitopes that it will bind specifically to the paratopic region of the host immune cell and elicit higher immune response. The results reported in this study could be supported by in vivo and in vitro laboratory evaluations in the future. Nevertheless, these findings provide a novel way to develop vaccine against M. gallisepticum.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

Author Contributions

SPM and HMC designed and performed the experimental studies. SPM carried out the in silico experiments. The manuscript was written by SPM and HMC.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

The authors thank the management of Thiruvalluvar University for providing the infrastructural facilities to conduct this research work.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2021.721061/full#supplementary-material

References

1. Levisohn S, Kleven SH. Avian mycoplasmosis (Mycoplasma gallisepticum). Revue Scientifique Technique. (2000) 19:425–42. doi: 10.20506/rst.19.2.1232

CrossRef Full Text | Google Scholar

2. Yoder HWJ. Mycoplasma gallisepticum Infection. Ames, IA: Iowa State University Press (1991).

3. Hochachka WM, Dhondt AA. Density-dependent decline of host abundance resulting from a new infectious disease. Proc Natl Acad Sci USA. (2000) 97:5303–6. doi: 10.1073/pnas.080551197

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Ley DH. Mycoplasma gallisepticum infection. In: Saif YM, editor. Diseases of Poultry. 11th ed. Ames, IA: Iowa State Press (2003). p. 722–43.

5. Browning GF, Marenda MS, Markham PF, Noormohammadi AH, Whitear KG. Mycoplasma. In: Gyles CL, Prescott JF, Songer JG, Thoen CO, editors. Pathogenesis of Bacterial Infections in Animals. 4th ed Ames, IA: Blackwell Publishing Professional (2010). p. 549–73. doi: 10.1002/9780470958209.ch29

CrossRef Full Text | Google Scholar

6. Majumder S, Silbart LK. Interaction of Mycoplasma gallisepticum with chicken tracheal epithelial cells contributes to macrophage chemotaxis and activation. Infect Immun. (2015) 84:266–74. doi: 10.1128/IAI.01113-15

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Mizutani M, Miyata M. Behaviors and energy source of Mycoplasma gallisepticum gliding. J Bacteriol. (2019) 201:e00397-9. doi: 10.1128/JB.00397-19

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Shahbaaz M, Bisetty K, Ahmad F, Hassan MI. In silico approaches for the identification of virulence candidates amongst hypothetical proteins of Mycoplasma pneumoniae. Comput Biol Chem. (2015) 59:67–80. doi: 10.1016/j.compbiolchem.2015.09.007

PubMed Abstract | CrossRef Full Text | Google Scholar

9. May M, Papazisi L, Gorton TS, Geary SJ. Identification of fibronectin-binding proteins in Mycoplasma gallisepticum strain R. Infect Immun. (2006) 74:1777–85. doi: 10.1128/IAI.74.3.1777-1785.2006

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Pflaum K, Tulman ER, Beaudet J, Canter J, Geary SJ. Variable lipoprotein hemagglutinin A gene (vlhA) expression in variant Mycoplasma gallisepticum strains in vivo. Infect Immun. (2018) 86:e00524-18. doi: 10.1128/IAI.00524-18

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Beaudet J, Tulman ER, Pflaum K, Canter JA, Silbart LK, Geary SJ. Immunologic pathways in protective versus maladaptive host responses to attenuated and pathogenic strains of Mycoplasma gallisepticum. Infect Immun. (2019) 87:e00613-8. doi: 10.1128/IAI.00613-18

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Pflaum K, Tulman ER, Beaudet J, Liao X, Geary SJ. Global changes in Mycoplasma gallisepticum phase-variable lipoprotein gene vlhA expression during in vivo infection of the natural chicken host. Infect Immun. (2015) 84:351–5. doi: 10.1128/IAI.01092-15

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Ron M, Gorelick-Ashkenazi A, Levisohn S, Nir-Paz R, Geary SJ, Tulman E, et al. Mycoplasma gallisepticum in vivo induced antigens expressed during infection in chickens. Vet Microbiol. (2015) 175:265–74. doi: 10.1016/j.vetmic.2014.12.007

PubMed Abstract | CrossRef Full Text | Google Scholar

14. El Gazzar M, Laibinis VA, Ferguson-Noel N. Characterization of a ts-1-like Mycoplasma galisepticum isolate from commercial broiler chickens. Avian Dis. (2011) 55:569–74. doi: 10.1637/9689-021711-Reg.1

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Kleven SH. Control of avian mycoplasma infections in commercial poultry. Avian Dis. (2008) 52:367–74. doi: 10.1637/8323-041808-Review.1

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Chernova OA, Medvedeva ES, Mouzykantov AA, Baranova NB, Chernov VM. Mycoplasmas and their antibiotic resistance: the problems and prospects in controlling infections. Acta Naturae. (2016) 8:24–34. doi: 10.32607/20758251-2016-8-2-24-34

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Jumbe N, Louie A, Leary R, Liu W, Deziel MR, Tam VH, et al. Application of a mathematical model to prevent in vivo amplification of antibiotic-resistant bacterial populations during therapy. J Clin Invest. (2003) 112:275–85. doi: 10.1172/JCI200316814

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Mavromati J, Mavromati E, Gjeta Z. The Effect of a macrolid antibiotic on the control of mycoplasmas and production efficiency in broilers. Biotechnol Anim Husb. (2011) 27:721–31. doi: 10.2298/BAH1103721M

CrossRef Full Text | Google Scholar

19. Ishfaq M, Zhang W, Ali Shah SW, Wu Z, Wang J, Ding L, et al. The effect of Mycoplasma gallisepticum infection on energy metabolism in chicken lungs: through oxidative stress and inflammation. Microb Pathog. (2020) 138:103848. doi: 10.1016/j.micpath.2019.103848

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Gaurivaud P, Baranowski E, Pau-Roblot C, Sagné E, Citti C, Tardy F. Mycoplasma agalactiae secretion of β-(1 → 6)-Glucan, a rare polysaccharide in prokaryotes, is governed by high-frequency phase variation. Appl Environ Microbiol. (2016) 82:3370–83. doi: 10.1128/AEM.00274-16

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Safavi A, Kefayat A, Sotoodehnejadnematalahi F, Salehi M, Modarressi MH. Production, purification, and in vivo evaluation of a novel multiepitope peptide vaccine consisted of immunodominant epitopes of SYCP1 and ACRBP antigens as a prophylactic melanoma vaccine. Int Immunopharmacol. (2019) 76:105872. doi: 10.1016/j.intimp.2019.105872

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Kardani K, Hashemi A, Bolhassani A. Comparative analysis of two HIV-1 multiepitope polypeptides for stimulation of immune responses in BALB/c mice. Mol Immunol. (2020) 119:106–22. doi: 10.1016/j.molimm.2020.01.013

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Safavi A, Kefayat A, Mahdevar E, Ghahremani F, Nezafat N, Modarressi MH. Efficacy of co-immunization with the DNA and peptide vaccines containing SYCP1 and ACRBP epitopes in a murine triple-negative breast cancer model. Hum Vaccin Immunother. (2021) 17:22–34. doi: 10.1080/21645515.2020.1763693

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Magnan CN, Zeller M, Kayala MA, Vigil A, Randall A, Felgner PL, et al. High-throughput prediction of protein antigenicity using protein microarray data. Bioinformatics. (2010) 26:2936–43. doi: 10.1093/bioinformatics/btq551

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Ong E, Wang H, Wong MU, Seetharaman M, Valdez N, He Y. Vaxign-ML: supervised machine learning reverse vaccinology model for improved prediction of bacterial protective antigens. Bioinformatics. (2020) 36:3185–91. doi: 10.1093/bioinformatics/btaa119

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Dimitrov I, Flower DR, Doytchinova I. AllerTOP–a server for in silico prediction of allergens. BMC Bioinform. (2013) 14:4. doi: 10.1186/1471-2105-14-S6-S4

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Ali SA, Almofti YA, Abd-Elrahman KA. Immunoinformatics approach for multiepitopes vaccine prediction against glycoprotein b of avian infectious laryngotracheitis virus. Adv Bioinformatics. (2019) 2019:1270485. doi: 10.1155/2019/1270485

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Aziz F, Tufail S, Shah MA, Salahuddin Shah M, Habib M, Mirza O, et al. In silico epitope prediction and immunogenic analysis for penton base epitope-focused vaccine against hydropericardium syndrome in chicken. Virus Res. (2019) 273:197750. doi: 10.1016/j.virusres.2019.197750

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Vainio O, Koch C, Toivanen A. B-L antigens (class II) of the chicken major histocompatibility complex control T-B cell interaction. Immunogenetics. (1984) 19:131–40. doi: 10.1007/BF00387856

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Hala K, Boyd R, Wick G. Chicken major histocompatibility complex and disease. Scand J Immunol. (1981) 14:607–16. doi: 10.1111/j.1365-3083.1981.tb00602.x

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Valdivia-Olarte H, Requena D, Ramirez M, Saravia LE, Izquierdo R, et al. Design of a predicted MHC restricted short peptide immunodiagnostic and vaccine candidate for Fowl adenovirus C in chicken infection. Bioinformation. (2015) 11:460–5. doi: 10.6026/97320630011460

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Karosiene E, Lundegaard C, Lund O, Nielsen M. NetMHCcons: a consensus method for the major histocompatibility complex class I predictions. Immunogenetics. (2012) 64:177–86. doi: 10.1007/s00251-011-0579-8

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Calis JJ, Maybeno M, Greenbaum JA, Weiskopf D, De Silva AD, Sette A, et al. Properties of MHC class I presented peptides that enhance immunogenicity. PLoS Comput Biol. (2013) 9:e1003266. doi: 10.1371/journal.pcbi.1003266

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Andreatta M, Karosiene E, Rasmussen M, Stryhn A, Buus S, Nielsen M. Accurate pan-specific prediction of peptide-MHC class II binding affinity with improved binding core identification. Immunogenetics. (2015) 67:641–50. doi: 10.1007/s00251-015-0873-y

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Larsen JE, Lund O, Nielsen M. Improved method for predicting linear B-cell epitopes. Immunome Res. (2006) 2:2. doi: 10.1186/1745-7580-2-2

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Cheng J, Randall AZ, Sweredoski MJ, Baldi P. SCRATCH: a protein structure and structural feature prediction server. Nucleic Acids Res. (2005) 33:W72–6. doi: 10.1093/nar/gki396

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Gasteiger E, Hoogland C, Gattiker A, Wilkins MR, Appel RD, Bairoch A. Protein identification and analysis tools on the ExPASy server. In: The Proteomics Protocols Handbook Springer. Humana Press (2005). p. 571–607. doi: 10.1385/1-59259-890-0:571

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Garnier J, Gibrat J-F, Robson B. GOR method for predicting protein secondary structure from amino acid sequence. Methods Enzymol. (1996) 266:540–53. doi: 10.1016/S0076-6879(96)66034-0

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Buchan DW, Jones DT. The PSIPRED protein analysis workbench: 20 years on. Nucleic Acids Res. (2019) 47:W402–7. doi: 10.1093/nar/gkz297

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Yang J, Yan R, Roy A, Xu D, Poisson J, Zhang Y. The I-TASSER Suite: protein structure and function prediction. Nat methods. (2015) 12:7–8. doi: 10.1038/nmeth.3213

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Heo L, Park H, Seok C. GALAXYREFINE: protein structure refinement driven by side-chain repacking. Nucleic Acids Res. (2013) 41:W384–8. doi: 10.1093/nar/gkt458

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Ko J, Park H, Heo L, Seok C. GalaxyWEB server for protein structure prediction and refinement. Nucleic Acids Res. (2012) 40:W294–7. doi: 10.1093/nar/gks493

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Laskowski RA, MacArthur MW, Moss DS, Thornton JM. PROCHECK: a program to check the stereochemical quality of protein structures. J Appl Cryst. (1993) 26:283–91. doi: 10.1107/S0021889892009944

CrossRef Full Text | Google Scholar

44. Saha S, Raghava GPS. BcePred:prediction of continuous B-cell epitopes in antigenic sequences using physico-chemical properties. In: Nicosia G, Cutello V, Bentley PJ, Timis J, editors. ICARIS 2004, LNCS 3239. Springer (2014). p. 197–204. doi: 10.1007/978-3-540-30220-9_16

CrossRef Full Text | Google Scholar

45. Ponomarenko JV, Bui H, Li W, Fusseder N, Bourne PE, Sette A, et al. ElliPro: a new structure-based tool for the prediction of antibody epitopes. BMC Bioinform. (2008) 9:514. doi: 10.1186/1471-2105-9-514

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Craig DB, Dombkowski AA. Disulfide by design 2.0: a web-based tool for disulfide engineering in proteins. BMC Bioinform. (2013) 14:346. doi: 10.1186/1471-2105-14-346

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Vinkler M, Bainová H, Bryja J. Protein evolution of Toll-like receptors 4, 5 and 7 within Galloanserae birds. Genet Sel Evol. (2014) 46:72. doi: 10.1186/s12711-014-0072-6

PubMed Abstract | CrossRef Full Text | Google Scholar

48. de Vries SJ, Bonvin AM. CPORT: a consensus interface predictor and its performance in prediction-driven docking with HADDOCK. PLoS ONE. (2011) 6:e17695. doi: 10.1371/journal.pone.0017695

PubMed Abstract | CrossRef Full Text | Google Scholar

49. van Zundert G, Rodrigues J, Trellet M, Schmitz C, Kastritis PL, Karaca E, et al. The HADDOCK2.2 Web Server: user-friendly integrative modeling of biomolecular complexes. J Mol Biol. (2016) 428:720–5. doi: 10.1016/j.jmb.2015.09.014

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Vangone A, Bonvin AM. Contacts-based prediction of binding affinity in protein–protein complexes. Elife. (2015) 4:e07454. doi: 10.7554/eLife.07454.018

CrossRef Full Text | Google Scholar

51. Xue LC, Rodrigues JP, Kastritis PL, Bonvin AM, Vangone A. PRODIGY: a web server for predicting the binding affinity of protein–protein complexes. Bioinformatics. (2016) 32:3676–8. doi: 10.1093/bioinformatics/btw514

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Laskowski RA, Jabłońska J, Pravda L, Vareková RS, Thornton JM. PDBsum: structural summaries of PDB entries. Protein Sci. (2018) 27:129–34. doi: 10.1002/pro.3289

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Grote A, Hiller K, Scheer M, Munch R, Nortemann B, Hempel DC, et al. JCat: a novel tool to adapt codon usage of a target gene to its potential expression host. Nucleic Acids Res. (2005) 33:526–31. doi: 10.1093/nar/gki376

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Castiglione F, Mantile F, De Berardinis P, Prisco A. How the interval between prime and boost injection affects the immune response in a computational model of the immune system. Comput Math Methods Med. (2012) 2012:842329. doi: 10.1155/2012/842329

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Rapin N, Lund O, Bernaschi M, Castiglione F. Computational immunology meets bioinformatics: the use of prediction tools for molecular binding in the simulation of the immune system. PLoS ONE. (2010) 5:e9862. doi: 10.1371/journal.pone.0009862

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Kolesanova EF, Sanzhakov MA, Kharybin ON. Development of the schedule for multiple parallel “difficult” Peptide synthesis on pins. Internat J peptides. (2013) 2013:197317. doi: 10.1155/2013/197317

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Kar T, Narsaria U, Basak S, Deb D, Castiglion F, Mueller DM, et al. A candidate multi-epitope vaccine against SARS-CoV-2. Sci Rep. (2020) 10:10895. doi: 10.1038/s41598-020-67749-1

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Kolesanova EF, Melnikova MV, Bolshakova TN, Rybalkina EY, Sivov IG. Bacteriophage MS₂ as a tool for targeted delivery in solid tumor chemotherapy. Acta naturae. (2019) 11:98–101. doi: 10.32607/20758251-2019-11-2-98-101

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Liu YF, Lin CY, Hong HM. In silico design, synthesis and potency of an epitope-based vaccine against foot-and-mouth disease virus. Int. J. Pharmacol. (2017) 13:122–33. doi: 10.3923/ijp.2017.122.133

CrossRef Full Text | Google Scholar

60. Michel-Todó L, Bigey P, Reche PA, Pinazo MJ, Gascón J, Alonso-Padilla J. Design of an epitope-based vaccine ensemble for animal trypanosomiasis by computational methods. Vaccines. (2020) 8:130. doi: 10.3390/vaccines8010130

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Osman MM, ElAmin EE, Al-Nour MY, Alam SS, Adam RS, Ahmed AA, et al. In silico design of epitope based peptide vaccine against virulent strains of (HN)-Newcastle Disease Virus (NDV) in poultry species. Int J Multidiscip Cur Res. (2016) 868–78.

62. Hasan M, Ghosh P, Azim KF, Mukta S, Abir RA, Nahar J, et al. Reverse vaccinology approach to design a novel multi-epitope subunit vaccine against avian influenza A (H7N9) virus. Microb Pathog. (2019) 130:19–37. doi: 10.1016/j.micpath.2019.02.023

PubMed Abstract | CrossRef Full Text | Google Scholar

63. Ingale AG, Goto S. Prediction of CTL epitope, in silico modeling and functional analysis of cytolethal distending toxin (CDT) protein of Campylobacter jejuni. BMC Res Notes. (2014) 7:92. doi: 10.1186/1756-0500-7-92

PubMed Abstract | CrossRef Full Text | Google Scholar

64. Unni PA, Ali AMMT, Rout M, Thabitha A, Vino S, Lulu SS. Designing of an epitope-based peptide vaccine against walking pneumonia: an immunoinformatics approach. Mol Biol Rep. (2019) 46:511–27. doi: 10.1007/s11033-018-4505-0

PubMed Abstract | CrossRef Full Text | Google Scholar

65. Jabbar B, Rafique S, Salo-Ahen OMH, Ali A, Munir M, Idrees M, et al. Antigenic peptide prediction from E6 and E7 oncoproteins of HPV types 16 and 18 for therapeutic vaccine design using immunoinformatics and MD simulation analysis. Front Immunol. (2018) 9:3000. doi: 10.3389/fimmu.2018.03000

PubMed Abstract | CrossRef Full Text | Google Scholar

66. Tosta SFO, Passos MS, Kato R, Salgado A, Xavier J, Jaiswal AK, et al. Multi-epitope based vaccine against yellow fever virus applying immunoinformatics approaches. J Biomol Struct Dyn. (2021) 39:219–35. doi: 10.1080/07391102.2019.1707120

PubMed Abstract | CrossRef Full Text | Google Scholar

67. Behbahani M, Moradi M, Mohabatkar H. In silico design of a multi-epitope peptide construct as a potential vaccine candidate for Influenza A based on neuraminidase protein. In Silico Pharmacol. (2021) 9:36. doi: 10.1007/s40203-021-00095-w

PubMed Abstract | CrossRef Full Text | Google Scholar

68. Sanches Rodrigo CO, Tiwari S, Ferreira Laís CG, Oliveira Flávio M, Lopes Marcelo D, Passos Maria JF, et al. Immunoinformatics design of multi-epitope peptide-based vaccine against schistosoma mansoni using transmembrane proteins as a target. Front Immunol. (2021) 12:621706. doi: 10.3389/fimmu.2021.621706

PubMed Abstract | CrossRef Full Text | Google Scholar

69. Ashfaq UA, Saleem S, Masoud MS, Ahmad M, Nahid N, et al. Rational design of multi epitope-based subunit vaccine by exploring MERS-COV proteome: reverse vaccinology and molecular docking approach. PLoS ONE. (2021) 16:e0245072. doi: 10.1371/journal.pone.0245072

PubMed Abstract | CrossRef Full Text | Google Scholar

70. Hossain MS, Hossan MI, Mizan S, Moin AT, Yasmin F, Akash A, et al. Immunoinformatics approach to designing a multi-epitope vaccine against Saint Louis Encephalitis Virus. Informatics Med Unlocked. (2021) 22:100500. doi: 10.1016/j.imu.2020.100500

CrossRef Full Text | Google Scholar

71. Shey RA, Ghogomu SM, Esoh KK, Nebangwa ND, Shintouo CM, Nongley NF, et al. In-silico design of a multi-epitope vaccine candidate against onchocerciasis and related filarial diseases. Sci Rep. (2019) 9:4409. doi: 10.1038/s41598-019-40833-x

PubMed Abstract | CrossRef Full Text | Google Scholar

72. Henderson GP, Jensen GJ. Three-dimensional structure of Mycoplasma pneumoniae's attachment organelle and a model for its role in gliding motility. Mol Microbiol. (2006) 60:376–85. doi: 10.1111/j.1365-2958.2006.05113.x

PubMed Abstract | CrossRef Full Text | Google Scholar