key: cord-0017958-0e8gpu9u
authors: Mugunthan, Susithra Priyadarshni; Harish, Mani Chandra
title: Multi-epitope-Based Vaccine Designed by Targeting Cytoadherence Proteins of Mycoplasma gallisepticum
date: 2021-05-17
journal: ACS Omega
DOI: 10.1021/acsomega.1c01032
sha: c4f9faab6d8efb3feb7b74c423c69735b0346e71
doc_id: 17958
cord_uid: 0e8gpu9u

[Image: see text] Mycoplasma gallisepticum causes chronic respiratory disease in chickens leading to large economic losses in the poultry industry, and the impacts remain to be a great challenge for a longer period. Among the other approaches, a vaccine targeting the adhesion proteins of M. gallisepticum would be a promising candidate in controlling the infection. Thus, the present study is aimed to design a multi-epitope vaccine candidate using cytoadhesion proteins of M. gallisepticum through an advanced immunoinformatics approach. As a result, the multi-epitope vaccine was constructed, which comprised potential T-cell and B-cell binding epitopes with appropriate adjuvants. The designed multi-epitope vaccine represented high antigenicity with viable physiochemical properties. The prospective three-dimensional structure of the epitope was predicted, refined, and validated. The molecular docking analysis of multi-epitope vaccine candidates with the chicken Toll-like receptor-5 predicted effective binding. Furthermore, codon optimization and in silico cloning ensured high expression. Thus, the present finding indicates that the engineered multi-epitope vaccine is structurally stable and can induce a strong immune response. Furthermore, the multi-epitope vaccine is suggested to be a suitable vaccine candidate for the M. gallisepticum infection due to its effective binding capacity and precise specificity.

Mycoplasma gallisepticum is an avian pathogen causing chronic respiratory disease in chickens 1−4 and infectious sinusitis in turkeys. 5 M. gallisepticum are cell wall-less bacteria belonging to the family Mycoplasmataceae and class Mollicutes 6 and are usually characterized by small cells and a genome with an absence of a peptidoglycan layer. 7 The remarkable characteristic of M. gallisepticum is colonization of host cells by cytoadhesion, followed by cilial shedding and inflammation. 8 Initial adhesion of M. gallisepticum to the host cells is a crucial step in infection and pathogenesis and along with a family of variable lipoproteins (vlhA) it plays an important role in pathogenesis by immune evasion. 9 The M. gallisepticum cytoadherence is a multistep progression involving the primary cytoadherence proteins GapA and CrmA 10 and cytoadherence accessory proteins Hlp3 and PlpA. 11 The cytoadherence proteins and cytoadherence accessory proteins are virulence factors in M. gallisepticum infection, and prophylactic strategies targeting these proteins may be an important preventive measure against M. gallisepticum infection.

Avian mycoplasmosis has caused a substantial economic loss globally by significantly reducing egg production, meat quality, and hatchability rate. 12 Control of M. gallisepticum by antibiotics is practiced in most poultry farms, and the widely used antibiotics are tetracyclines, macrolides (tylosin and tilmicosin), fluoroquinolones (enrofloxacin and difloxacin), and pleuromutilins (tiamulin). 13, 14 The widespread use of antibiotics has led to the rise and spread of resistant bacteria over time. 15 The quinolone resistance in M. gallisepticum has become common, which hinders treatment and control efforts. 16 The application of vaccinology might be a potential approach in the prevention and control of M. gallisepticum infections in layer and breeder flocks. The currently available vaccines are live-attenuated vaccines, live vaccines, and bacterins. 17 Live-attenuated vaccines such as the F-strain, ts-11 strain, and 6/85 strain and bacterins are widely used, but a sudden onset of M. gallisepticum infections can pose a challenge for the suppression of infection. 18, 19 Liveattenuated vaccines may get reverted back and cause side effects, whereas bacterins are highly expensive and involve repetitive doses. Thus, the development of highly efficient vaccines at low cost is of immense importance, which can be effectively achieved through new recombinant strategies. Reports on in silico approaches for the vaccine development for mycoplasma species are not observed, except the identification of some virulence genes. 20 Since a large population of birds is to be administered concurrently, effortless administration with a least amount of labor is a vital requisite to consider the vaccine desirable for poultry. 21 To fix the current state of mycoplasma infections in chickens, an effective vaccine is therefore necessary.

In both the preclinical and clinical phases, the development of a vaccine candidate is expensive, complex, and time-consuming, requiring in vitro and in vivo protocols for assessing efficacy. The computational methods have facilitated the in silico evaluations for an intended vaccine candidate with the advancement of bioinformatics, significantly reducing the needed time for preclinical and clinical procedures. 22, 23 The recent advancement in the field of bioinformatics and immunoinformatics had led to the efficient screening and validation of proteins, which can be a potential vaccine target. 24 Earlier studies reported the use of in silico assays to investigate the efficiency of the designed vaccine for human diseases. 25 The advantage of this approach is the identification of a larger range of vaccine candidates without the need to propagate the pathogenic microorganism in a laboratory while delivering a multivalent, highly purified, and sustainable alternative to whole-cell vaccines. 26 Furthermore, the molecular docking approach predicts the binding domains and molecular interactions between the receptor and vaccine candidate. 27, 28 To date, a very few in silico studies were carried out to design a vaccine in veterinary diseases and no studies have been reported for M. gallisepticum. 29, 30 The major limitations of the in silico approach is that it cannot predict polysaccharides or lipids which may be active components in proteins, 31 the requirement to identify suitable antigen-processing sites for the predicted epitopes which were not addressed by many epitope prediction tools. 32 In this study, we used robust immunoinformatics tools to predict MHC-I and MHC-II binding T-cell-specific epitopes and B-cell epitopes from four M. gallisepticum antigenic proteins and developed a base to design a candidate for a multi-epitopebased vaccine against M. gallisepticum infection in chickens ( Figure 1 ).

Vaccines are critically responsible for providing host organisms with protection from a specific disease, helping to save lives worldwide. The process for the production of vaccines is typically highly painstaking and costly and takes a long period of time for their achievement. The design of vaccines can be accomplished by using different immunoinformatics methods in less time. They include the detection of a pathogen's latent antigenic proteins, accompanied by the determination of different immune-dominant epitopes that are responsible for the development of immune responses against the pathogen, both humoral and cell-mediated. Therefore, multi-epitopebased vaccines can be engineered for the probable antigenic proteins of a pathogen by identifying the B-and T-cell epitope regions. Currently, numerous peptide vaccines are under development, and most of them are for human diseases. 33 Limited research studies have been carried out in the field of in silico vaccines for poultry and other animals. Among those, few studies focused on the prediction of epitopes for animal diseases such as foot and mouth disease 33 and animal trypanosomiasis 34 in designing the multi-epitope vaccine, which offered effective immunity when compared with the available vaccines. The immunoinformatics approach is a valid method for designing multi-epitope vaccines against infectious diseases in poultry. 35−38 The immunoinformatics analysis and evaluation in the present study show that the vaccine construct has all the factors within the most favorable range for it to be considered as a potential vaccine candidate, which further needs to be experimentally validated.

Antigenicity and Allergenicity Prediction. Initially, the GapA, PlpA, Hlp3, and CrmA proteins were subjected to antigenicity and allergenicity prediction. Verifying the antigenic and allergic nature of the selected proteins is crucial in the selection of epitopes for the further study. 39 The VaxiJen v 2.0 server was used for the evaluation of the antigenicity of the proteins. The threshold score for antigenicity is 0.4. All protein sequences selected are nonallergens that are predicted using AllerTop v.2.0.0. Table 1 summarizes both the proteins' antigenic score and allergenicity.

T-Cell Epitope Prediction. Due to the lack of chicken MHC alleles, human MHC alleles were used to predict the epitopes. Earlier studies used human alleles because of the lack of chicken alleles in immunoinformatics online tools. 30, 40, 41 Studies have shown that chicken B−F alleles can stimulate an immune response similar to human class I homologous alleles, particularly in antigen presentation. 42, 43 Due to the lack of chicken MHC-I and MHC-II alleles in the T-cell epitope prediction database, the closest similar human alleles (HLA*B 40:06, HLA*B 41:04, and HLA*B 41:03 for MHC I and DRB1:1482, DRB1:1366, DRB1:1310, and DRB1:1445 for MHC II) were used. The NetMHCcons 1.1 server predicted MHC-I binding epitopes. The epitopes with % rank <0.50 were preferred as strong binders for MHC I (Tables 2 and S1). The NetMHCpan server predicted MHC-II binding epitopes. The epitopes with % rank <2 were preferred as strong binders for MHC II (Tables 3 and S2) .

B-Cell Epitope Prediction. An important step in the vaccine design includes the prediction and identification of Bcell epitopes. 39 The B-cell epitopes in the cytoadherence proteins were predicted using the BepiPred server. The threshold was 0.350, and the peptides above this threshold are considered as potential B-cell epitopes ( Figure 2 ). The epitopes from all the selected proteins with length varying from 10 to 15 peptides are listed in Table 4 .

Multi-epitope Vaccine Design. The initiation of adaptive immunity largely depends on the regulation of T cells in poultry. 44 Avian β-defensin, an adjuvant, was added at the Nterminal to enhance the efficacy of the vaccine, followed by a sequence of T-cell and B-cell epitopes present on the cytoadhesion and cytoadhesion accessory proteins, viz., GapA, PlpA, Hlp3, and CrmA ( Figure 3 ). 45 These epitopes along with the adjuvant have been merged together with EAAAK, which enhances the bifunctional catalytic activity and firmness of the vaccine construct. 46 The AAY and GPGPG linkers were inserted between the T-cell and B-cell epitopes, respectively, which allows dissociation and identification of individual epitopes. 47, 48 mRNA Structure Prediction. Prediction of the mRNA secondary structure is essential, which is a key factor in translation initiation, elongation, and mRNA biogenesis. 49, 50 Free energy related to the whole mRNA structure was obtained using Mfold online service. The minimum free energy of the secondary RNA structure was ΔG = −156.50 kcal/mol ( Figure  4 ). This score indicates the stability and translation efficiency of protein in hosts. The increase in mRNA stability is directly related to the increases in the expression rate.

Antigenicity, Allergenicity, and Solubility Profile Prediction of the Multi-epitope Vaccine. A creditable vaccine candidate should be capable of initiating an immune response without causing allergic reactions, thus the antigenicity and allergenicity of the vaccine construct were predicted. 39, 41 The main criteria that have to be ensured while designing a vaccine is the antigenicity to induce a humoral and/or cellmediated immune response against the targeted microorganism and allergenicity of the constructed vaccine. The antigenic scores for the multi-epitope vaccine are 0.5669 and 0.5532 as predicted by VaxiJen v2.0 and ANTIGENpro servers, respectively. The vaccine was found to be nonallergenic using AllerTOP v2.0. The predicted solubility upon over expression by the SOLpro server showed the vaccine construct as soluble with probability 0.636137. The designed vaccine construct is antigenic and non-allergenic and can be endorsed as a suitable candidate.

Physiochemical Characterization of the Designed Vaccine. The assessment of physiochemical properties of the designed multi-epitope vaccine is essential for validating the safety and efficacy of the vaccine. 51 Thus, various physical and chemical properties were analyzed. The designed multi-epitope vaccine is composed of 196 amino acids with a molecular weight of∼ 21.40 kDa. The theoretical pI was calculated as 9.41, which indicate that the multi-epitope vaccine is significantly basic in nature. The instability index was estimated to be 35.78, implying a stable protein. The computed aliphatic index and grand average of hydropathicity (GRAVY) were found to be 92.14 and −0.094, respectively, signifying that the vaccine is thermostable and hydrophilic (Table 5) .

Secondary and Tertiary Structure Prediction. Evaluating the secondary structure of the multi-epitope vaccine revealed Figure 5 shows the graphical descriptions of the GOR IV server and the performance of the secondary structure from the PSIPRED server. I-TASSER servers have been used to model the tertiary structure of the vaccine. Based on 10 threading templates, the I-TASSER server predicted five tertiary structure models of the designed chimeric protein. All the 10 templates chosen showed good alignment according to their Z-score values (ranging from 1.71 to 2.25). The C-score value of the selected model is −4.83. The C-score value is directly proportional to the quality of the structure. With an estimated root-mean-square deviation (rmsd) of 17.3 ± 2.7 Å, this model had an estimated TM-score of 0.22 ± 0.06. As a scale for evaluating the structural similarity between two structures, the TM-score has been proposed. 52 A TM-score >0.5 reveals a correct topology model, and a TM-score <0.17 implies random similarity. These cut-off values are protein length-independent. Structure Refinement and Validation. I-TASSERmodeled three-dimensional (3D) structure was subjected to the PROCHECK server, where the Ramachandran plot was generated and the output revealed 37.1% residues that were present in the favored region. The Ramachandran plot allows one to visualize energetically allowed and disallowed dihedral angles psi (ψ) and phi (φ). Hence, the 3D model was subjected to further refinement by the GalaxyRefine server (http://galaxy. (Table  S3) . Thus, this model ( Figure 6 ) was finally selected for further investigations.

Immunological Property Assessment. The presence of B-cell epitopes in the multi-epitope vaccine construct plays an important role in triggering humoral immune responses. 53 The BcePred server was used for the prediction of continuous B-cell epitopes in the chimeric multi-epitope vaccine with default parameters. Important properties of epitopes such as hydrophilicity, antigenicity, flexibility, accessibility, polarity, and exposed surfaces were predicted (Table 6 ). Ellipro was used to predict the conformational B-cell epitopes in the multiepitope vaccine. Only one conformational epitope region was predicted (Table 7 and Figure 7 ). The analysis of continuous and discontinuous B-cell epitopes in the designed vaccine construct revealed that they can interact with antibodies and are flexible.

Vaccine Protein Disulfide Engineering. Disulfide engineering was performed to stabilize the vaccine construct by particular geometric confirmations. 36, 54 It was predicted that a total number of 60 pairs of amino acid residues may form a disulfide bond through the DbD2 server. Among them, only three residues, including LYS 2-TYR 5, ARG 49-ASP 66, and TYR 88-TYR 95, which were replaced by cysteine residues, were able to create disulfide bond formation following the residue assessment by chi3 and B-factor energy parameters (Figure 8 ). On the basis of −87 to +97 chi3 values and <2.5 energy values, residue screening was performed.

Molecular Docking. Molecular docking was carried out between the multi-epitope vaccine construct and chicken Tolllike receptor-5 (TLR5). Among the chicken TLRs, TLR5 has a propensity toward the recognition of bacterial components present on the extracellular surfaces, 71 thus chicken TLR5 was selected for docking. HADDOCK uses biochemical and biophysical interaction data to perform docking. 55 HADDOCK clustered 104 structures in 13 clusters, representing 52% of the models produced by water-refined HADDOCK. The most accurate cluster of all is the top cluster with the smallest HADDOCK ranking. Using the HADDOCK refinement server, a representative model of the top cluster was subjected to further refinement, where 20 structures were grouped into one cluster, Table 8 . A strong binding affinity between the vaccine and the receptor is indicated by a Haddock score of 297.1 ± 3.1; a negative score suggests better docking (Figure 9 ). A total of 19 hydrogen bonds, 1 salt bridge, and 182 nonbonded contacts were analyzed ( Table 9 ). The binding affinity of a complex, or in terms of thermodynamics, the Gibbs free energy (ΔG), is a critical quantity for deciding whether or not an interaction can actually occur under particular conditions in the cell. Using PRODIGY web servers, the binding affinity of the docked complex was analyzed. The G-values were found to be −14.4 kcal/mol for the vaccine construct and chicken TLR5, and the dissociation constant was 3.0 × 10 −11 . The results show that, as indicated by negative Gibbs free energy (ΔG), the docked complex is energetically viable. Codon Optimization and In Silico Cloning. The initial step in validating a vaccine candidate is to carry out serological analysis, 56 which requires the expression of the vaccine in a suitable expression system. A similar approach has been used in earlier studies for in silico designed vaccines prior to in vitro expression. 57 We chose E. coli as the expression system. The Java codon adaptation tool (JCAT) aids in optimized cloning and expression in the E. coli K12 strain as a host. The codon adaptation index (CAI) in the adapted sequence was 0.92, and the GC content of the improved sequence was estimated to be 54.08% and the GC content of E. coli was 50.7. The graphical representation of codon usage in the optimized multi-epitope vaccine gene is shown in Figure 10 . In order to perform in silico cloning, the codon-optimized vaccine construct sequence was inspected for restriction enzyme sites; HindIII and BamHI enzymes were not found in the vaccine sequence, so these enzymes were used for in silico cloning purpose. Finally, a successful clone of 5938 bp was obtained following the insertion of the vaccine construct into the pET28a(+) vector ( Figure 11 ).

A multi-epitope-based vaccine was developed against M. gallisepticum through an immunoinformatics approach, which is advantageous over traditional methods in many aspects. The application of computational methods can be used to design and develop an effective vaccine in a short time and low cost. The amino acid sequences of cytoadherence and cytoadherence accessory proteins were retrieved, their potential B-cell and Tcell binding epitopes were predicted, and a vaccine was constructed, including epitopic sequences with an immunogenic adjuvant merged together with the suitable linkers that can trigger strong immunogenic responses. The constructed multiepitope vaccine was found to be nonallergenic and -antigenic. Disulfide engineering further increased the stability of the multiepitope vaccine construct. A stable and strong interaction of the vaccine with the chicken immune receptor was confirmed by molecular docking analysis. The multi-epitope cDNA in silico Evaluation of Protein Antigenicity and Allergenicity. The selected protein sequences were submitted to the VaXiJen v 2.0 server to identify antigenicity, assisting in the prediction of potential antigens with a threshold score of 0.4. The models were tested using leave-one-out cross-sectional models; the validations were with 70−89% accuracy. 58 Similarly, the prediction of allergenicity was made using AllerTop v. 2.0. 59 

http://pubs.acs.org/journal/acsodf Article Therefore, allergens and nonallergens with possible sensitivity and protein specificity have been predicted. Prediction of T-Cell Epitopes. Human MHC-I and MHC-II alleles were used as alternates for chicken MHC alleles due to the nonavailability of chicken alleles for peptide-binding predictions in immunoinformatics tools. 40 Hence, the best suitable human alleles for chicken were chosen from a previous study. 41 MHC-I and -II Binding Prediction. The NetMHCcons 1.1 server was used to predict MHC-I binding affinity (http://www. cbs.dtu.dk/services/NetMHCcons); 60 the thresholds were selected for strong binding peptides with IC 50 values of 50 nM for strong binders and 500 nM for weak binders. Thus, IC 50 < 50 nM peptides were categorized as strong binders and 50 < IC 50 < 500 nM peptides were considered weak binders. Similarly, in terms of the percentage (percent) rank, the thresholds for strong binding peptides were set at 0.5 while those for weak binders were maintained at 2. The protein sequence was set to a length of nine amino acids. Using the IEDB analysis resource (http:// tools.iedb.org/immunogenicity), 61 the class I immunogenicity score was predicted. The NetMHCIIpan 3.1 server was used to predict the MHC-II binding epitope (http://www.cbs.dtu.dk/ services/NetMHCIIpan-3.1/); 62 the peptide length was set to 15 amino acids. The server then predicted a number of peptides binding to MHC alleles, from which only peptides which were recognized as strong binders (% rank <2) for all the four MHC alleles were selected.

B-Cell Epitope Prediction. BepiPred was used to predict linear B-cell epitopes (http://tools.iedb.org/bcell/) from the IEDB research resource. 63 By combining two residue properties with the hidden Markov model, BepiPred predicts continuous epitopes. BepiPred was assessed on the data set of epitopes extracted from the literature, AntiJen, and HIV databases.

Construction of the Multi-epitope Vaccine. The multiepitope vaccine against M. gallisepticum was designed using the MHC-I and MHC-II binding epitopes from the GapA, PlpA, Hlp3, and CrmA proteins. Furthermore, to improve the immunogenicity of the multi-epitope vaccine, the TLR5 agonist, avian β defensin 8, was selected as an adjuvant. The EAAAK linker was used to link the adjuvant with CTL epitopes. The intra-CTL and intra-HTL epitopes were linked using AAY and GPGPG linkers, respectively.

mRNA Secondary Structure Prediction. The Mfold online server was used to predict the RNA secondary structure of the multi-epitope vaccine (http://unafold.rna.albany.edu/ ?q=mfold). The Mfold server provides the prediction of true positive base pairs and structures within minimum ΔG thermodynamically. 64 Antigenicity, Allergenicity, and Solubility Profile Prediction. The antigenicity of the multi-epitope vaccine was predicted using the VaxiJen server (http://www.ddgpharmfac. net/vaxijen/VaxiJen/VaxiJen.html) and ANTIGENpro (http://scratch.proteomics.ics.uci.edu/); 65 the allergenicity of the vaccine was verified using AllerTOP v2.0 (http://www.ddgpharmfac.net/AllerTOP/). The AllerTOP v2.0 server utilizes machine learning techniques to sort the allergens, such as nearest neighbors, auto-and cross-variance transformation, and amino acid descriptors. The accuracy of this server is 85.3% at fivefold cross validation. 59 In addition, the SOLpro online tool was employed for the prediction of solubility upon overexpression of the protein construct in the E. coli host (http:// , S17, G20, F21, M22, R23, V24, P25, N26, N27, E28, A29, Q30, C31, E32, Q33, A34, I37, C38, S39, K40, D41, H42, C43, H47, T48, R49, A50,  F51, G52, H53, C54, Q55, R56, G57, V58, P59, C60, C61, V64, Y65, D66, E67, A68, A69, A70, K71, Q72, E73, F74, T75, G76, F77, D78, A79, A82,  Y83, E85, N86, Q87, Y88, Y89, P90, P91, A92, A93, A94, Y95, Q96, E97, Y98, D99, Y100, Y101, P102, P103, A104, A105, Y107, N108, E109, L116,  P118, G119, P120, G121, I122, D130, R131, S132, R133, A134, T135, N136, G137, P138, G139, P140, G141, K142, I143, H144, N145, K146, L148,  K149, N150, L151, A152, K153, L154, K155, K156, G157, P158, G159, P160, G161, F162, L165, Y166, L167, N168 , S169, N171, S172, L173, S174, F175, I176, G177, P178, G179, P180, G181, K182, T183, F184, S185, L186, N187, K188, G189, L190, N191, K192, V193, I194 ACS Omega http://pubs.acs.org/journal/acsodf Article scratch.proteomics.ics.uci.edu/), with 74% accurate prediction and a 10-fold cross-validation method. 66 Determination of Physiochemical Properties. The ExPASy-ProtParam server, accessible at (http://web.expasy. org/protparam/), has calculated the physiochemical properties of the vaccine construct. Multiple parameters, aliphatic index, instability index, half-life, isoelectric point (pI), molecular weight, and atomic composition, including grand average hydropathicity (GRAVY), were analyzed. The half-life of the protein represents the time that the molecule disappears after its synthesis in the cell. The instability index indicates the in vitro stability of a protein molecule. 67 Secondary Structure Prediction. The secondary structure of the multi-epitope vaccine construct was predicted using the Garnier-Osguthorpe-Robson (GOR IV) online server (https:// npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl?page=npsa_gor4. html) with a mean accuracy of 64.4% 68 and position-specific iterated prediction (PSIPRED) analysis on outputs from PSI-BLAST (http://bioinf.cs.ucl.ac.uk/psipred/). 69 Tertiary Structure Prediction, Refinement, and Validation. The tertiary structure was modeled using the I-TASSER server (https://zhanglab.ccmb.med.umich.edu/I-TASSER/). 70 The 3D structure modeled by I-TASSER was subjected to refinement by the GalaxyRefine server (http://galaxy.seoklab. org/cgi-bin/submit.cgi?type=REFINE). 71 This server replaces amino acids with high-probability rotamers and applies molecular dynamic simulation for overall structural relaxation. The output usually includes five refined models, with different parameter scores, comprising GDT-HA, rmsd, MolProbity, Clash score, Poor rotamers, and Rama favored. 72 The refined structure was validated by PROCHECK (https://servicesn.mbi. ucla.edu/PROCHECK/) which analyzes the stereochemical quality of a protein structure by analyzing the residue-by-residue geometry and overall structure geometry. 73 Immunological Analysis for the Chimeric Multiepitope Vaccine. The multi-epitope vaccine was analyzed for linear (continuous) and conformational (discontinuous) Bcell epitopes. Continuous B-cell epitopes were predicted by BcePred (https://webs.iiitd.edu.in/raghava/bcepred/index. html). 74 Discontinuous B-cell epitopes were predicted by ElliPro, which gives insights about linear and discontinuous antibody epitopes based on a protein antigen's 3D structure (http://tools.iedb.org/ellipro/). 75 Multi-epitope Vaccine Protein Disulfide Engineering. Disulfide engineering is an important biotechnological tool to design new disulfide bonds in the target protein through cysteine mutation of residues in the highly mobile region of protein.

Disulfide bonds provide substantial stability and reinforce the protein geometric conformation. For this purpose, the DbD2 online server was used (http://cptweb.cpt.wayne.edu/DbD2/ index.php). The pair of residues capable of forming a disulfide Molecular Docking. To generate a stable immune response, it is necessary for the vaccine to interact with the target immune cell receptor. In this study, interactions between the vaccine 

http://pubs.acs.org/journal/acsodf Article construct and chicken TLR5 were studied as they are bacterialsensing TLR. 77 CPORT 78 was used in predicting the active and passive residues for the interactions. The docking of the vaccine with chicken TLR5 was performed by HADDOCK 2.4 (https:// www.bonvinlab.org/software/haddock2.4/). 79 To refine the chosen cluster, HADDOCK Refinement Interface was used. The best structure after refinement from the docked complex was chosen and its binding affinity was calculated using the PRODIGY web server. 80, 81 Last, the interacting residues between the vaccine construct and chicken TLR5 were mapped using PDBsum (https://www.ebi.ac.uk/thorntonsrv/ databases/pdbsum/Generate.html). 82 Codon Optimization and In Silico Cloning of the Vaccine Construct. The JCAT server (http://www.jcat.de) was used to execute codon optimization according to E. coli codon usage. 83 Parameters that are significant to the effectiveness of gene expression including CAI and GC content adjustment were optimized. The assessment of cloning and expression of the vaccine construct in an appropriate expression vector was carried out using the SnapGene restriction cloning tool (GSL Biotech LLC 2015).

The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acsomega.1c01032.

List of predicted MHC-I binding epitopes, list of MHC-II binding epitopes, and (five refined models from the Galaxy refine server (PDF) 

Avian mycoplasmas, The mycoplasmas

Avian mycoplasmosis ( Mycoplasma gallisepticum)

Mycoplasma gallisepticum Infection

Infectious sinusitis in wild turkeys

Molecular biology and pathogenicity of mycoplasmas. Microbiol

Highlights of mycoplasma research: an historical perspective

Transcription and Translation

Role of phenotypic diversity in pathogenesis of avian mycoplasmosis

Role of the GapA and CrmA cytadhesins ofMycoplasma gallisepticum in promoting virulence and host colonization

Identification offibronectin-binding proteins in Mycoplasma gallisepticum strain R

The Effect of a macrolid antibiotic on the control of mycoplasmas And production efficiency in broilers

Characterization of a ts-1-like Mycoplasma galisepticum isolate from commercial broiler chickens

Control of avian mycoplasma infections in commercial poultry

Mycoplasmas and Their Antibiotic Resistance: The Problems and Prospects in Controlling Infections

Application of a mathematical model to prevent in vivo amplification of antibiotic-resistant bacterial populations during therapy

Effects of time-specific F-strain Mycoplasma gallisepticum inoculation overlays on prelay ts11-strain Mycoplasma gallisepticum inoculation on performance characteristics of commercial laying hens

Evaluation of Antibiotics to Control Mycoplasma gallisepticum in Broiler Breeder Chickens

Vaccine and vaccination for Mycoplasma gallisepticum (MG)

Mycoplasma agalactiae Secretion of β-(1→6)-Glucan, a Rare Polysaccharide in Prokaryotes, Is Governed by High-Frequency Phase Variation

Introduction to poultry vaccines and immunity

The process development challenge for a new vaccine

Vaccine safety evaluation: practical aspects in assessing benefits and risks

Molecular Binding Mechanism and Pharmacology Comparative Analysis of Noscapine for Repurposing against SARS-CoV-2 Protease

Multiepitope Subunit Vaccine to Evoke Immune Response against Acute Encephalitis

Vaccine Formulation and Optimization for Human Herpes Virus-5 through an Immunoinformatics Framework

Antimicrobial Peptide Designing and Optimization Employing Large-Scale Flexibility Analysis of Protein-Peptide Fragments

High biorecognizing aptamer designing and optimization against human herpes virus-5

In silico Prediction of Peptide based Vaccine against Fowlpox Virus (FPV)

Immunoinformatics Approach for Multiepitopes Vaccine Prediction against Glycoprotein B of Avian Infectious Laryngotracheitis Virus

The development of veterinary vaccines: a review of traditional methods and modern biotechnology approaches

In silico veritas? Potential limitations for SARS-CoV-2 vaccine development based on T-cell epitope prediction

Synthesis and Potency of an Epitope-based Vaccine Against Foot-and-mouth Disease Virus

Design of an Epitope-Based Vaccine Ensemble for Animal Trypanosomiasis by Computational Methods. Vaccines

In Silico Design of Epitope Based Peptide Vaccine against Virulent Strains of (HN)-Newcastle Disease Virus (NDV) in Poultry Species

Reverse vaccinology approach to design a novel multi-epitope subunit vaccine against avian influenza A (H7N9) virus

Prediction of CTL epitope, in silico modeling and functional analysis of cytolethal distending toxin (CDT) protein of Campylobacter jejuni

T cell epitope engineering: an avian H7N9 influenza vaccine strategy for pandemic preparedness and response

Designing of an epitope-based peptide vaccine against walking pneumonia: an immunoinformatics approach

In silico epitope prediction and immunogenic analysis for penton base epitope-focused vaccine against hydropericardium syndrome in chicken

B-L antigens (class II) of the chicken major histocompatibility complex control T-B cell interaction

Chicken major histocompatibility complex and disease

Cell-mediated immunity in poultry

Chicken avian β-defensin 8 modulates immune response via the mitogen-activated protein kinase signaling pathways in a chicken macrophage cell line

Structural basis and designing of peptide vaccine using PE-PGRS family protein of Mycobacterium ulceransAn integrated vaccinomics approach

Exploring dengue genome to construct a multi-epitope based subunit vaccine by utilizing immunoinformatics approach to battle against dengue infection

Exploring Leishmania secretory proteins to design B and T cell multi-epitope subunit vaccine using immunoinformatics approach

Immunoinformatics based approaches to design a novel multi epitope-based vaccine for immune response reinforcement against Leptospirosis

A multi-method and structure-based in silico vaccine designing against Echinococcus granulosus through investigating enolase protein

A candidate multi-epitope vaccine against SARS-CoV-2

Scoring function for automated assessment of protein structure template quality

Novel Immunoinformatics Approaches to Design Multi-epitope Subunit Vaccine for Malaria by Investigating

HADDOCK: a protein-protein docking approach based on biochemical or biophysical information

Peptides for immunological purposes: design, strategies and applications

2020, 69, 101413. (58) Doytchinova, I. A.; Flower, D. R. VaxiJen: a server for prediction of protective antigens, tumour antigens and subunit vaccines

AllerTOP-a server for in silico prediction of allergens

NetMHCcons: a consensus method for the major histocompatibility complex class I predictions

Properties of MHC class I presented peptides that enhance immunogenicity

Accurate pan-specific prediction of peptide-MHC class II binding affinity with improved binding core identification

Improved method for predicting linear B-cell epitopes

Mfold web server for nucleic acid folding and hybridization prediction

SCRATCH: a protein structure and structural feature prediction server

Highthroughput prediction of protein antigenicity using protein microarray data

Protein identification and analysis tools on the ExPASy server

GOR method for predicting protein secondary structure from amino acid sequence

The PSIPRED Protein Analysis Workbench: 20 years on

The I-TASSER Suite: Protein structure and function prediction

Protein structure refinement driven by side-chain repacking

GalaxyWEB server for protein structure prediction and refinement

PROCHECK: a program to check the stereochemical quality of protein structures

BcePred: Prediction of Continuous B-Cell Epitopes in Antigenic Sequences Using Physico-chemical Properties

ICARIS 2004

ElliPro: a new structure-based tool for the prediction of antibody epitopes

Disulfide by Design 2.0: a webbased tool for disulfide engineering in proteins

Protein evolution of Toll-like receptors 4, 5 and 7 within Galloanserae birds

CPORT: A consensus interface predictor and its performance in prediction-driven docking with HADDOCK

The HADDOCK2.2 web server: User-friendly integrative modeling of biomolecular complexes

PRODIGY: A web server for predicting the binding affinity of protein−protein complexes

PDBsum: Structural summaries of PDB entries

JCat: a novel tool to adapt codon usage of a target gene to its potential expression host