Abstract
Lymphoid leukosis is a poultry neoplastic disease caused by avian leukosis virus (ALV) and is characterized by high morbidity and variable mortality rates in chicks. Currently, no effective treatment and vaccination is the only means to control it. This study exploited the immunoinformatics approaches to construct multi-epitope vaccine against ALV. ABCpred and IEDB servers were used to predict B and T lymphocytes epitopes from the viral proteins, respectively. Antigenicity, allergenicity and toxicity of the epitopes were assessed and used to construct the vaccine with suitable adjuvant and linkers. Secondary and tertiary structures of the vaccine were predicted, refined and validated. Structural errors, solubility, stability, immune simulation, dynamic simulation, docking and in silico cloning were also evaluated.The constructed vaccine was hydrophilic, antigenic and non-allergenic. Ramchandran plot showed most of the residues in the favored and additional allowed regions. ProsA server showed no errors in the vaccine structure. Immune simulation showed significant immunoglobulins and cytokines levels. Stability was enhanced by disulfide engineering and molecular dynamic simulation. Docking of the vaccine with chicken’s TLR7 revealed competent binding energies.The vaccine was cloned in pET-30a(+) vector and efficiently expressed in Escherichia coli. This study provided a potent peptide vaccine that could assist in tailoring a rapid and cost-effective vaccine that helps to combat ALV. However, experimental validation is required to assess the vaccine efficiency.
Similar content being viewed by others
Introduction
Avian leukosis virus (ALV) is one of the most prevalent retroviruses causing neoplasms in chickens1,2. ALV related to the genus Alpharetrovirus and subfamily Orthoretrovirinae3. ALVs were shown causing avian leukemia/sarcoma in birds with major economic losses in the global poultry industry1,4. ALV causes high morbidity rates than mortality rates In the infected flocks which negatively impacts hatchability and egg productivity5,6,7. In addition to the neoplasia, other factors contributed to mortality, including delayed maturation and diminished fertility rate, which have serious economic implications4.The co-infections of the virulent ALV isolates with reticuloendotheliosis virus (REV) and marek’s disease virus (MDV) are the main factor for prevalence of avian neoplastic diseases in China that hampering the future control of disease3. Therefore the affected flocks showed lower egg production and hemorrhages in the skin surrounding the feather follicles and phalanges8,9. As a result of loss of appetite, chickens characterized by weight loss, abnormal feathers, pale combs and wattles, depression, paralysis and finally death8. Postmortem studies revealed enlarged livers, spleens, hearts, kidneys, and congested lungs in the affected chickens7,9,10.
Horizontal, vertical, and genetic transmissions are the main three ways of ALV transmission8. The horizontal transmission occurs through direct contact with infected chickens and materials. In vertical infection, the disease is transmitted from hens to their off-springs via the eggs. Genetic transmission is through endogenous retroviruses occurred through a mendelian manner8. Another possible transmission route is through contaminated live and attenuated vaccines as some ALV strains were isolated from contaminated live vaccines11.
ALV is a member of the leukosis/sarcoma (L/S) group of avian retroviruses. It is an RNA virus family that uses the enzyme reverse transcriptase which originates the generation of the DNA provirus that is combined with the host genome during viral replication12. In terms of structure, ALV is a simple virus and has three chief genes located from its 5′-to-3′ end of its genome. These genes are gag/pro, pol and env genes, which code for the antigens and proteases that distinguish a different group of viron, as well as, the reverse transcriptase enzyme and the glycoprotein or the envelope protein12. The genomic sequence associated with an ordinance of the viral replications that flank these structural genes in the DNA provirus, forming long terminal repeats (LTRs) for the virus with promoters and enhancers13. Based on their viral envelope characteristics, chicken leukosis virus envelope subgroups are categorized into seven subgroups identified as A, B, C, D, E, J and K. The exogenous subgroups A, B, C, D, and J are highly pathogenic, while those of subgroup E viruses are endogenous and less pathogenic10,14. Virus subgroups are mainly differentiated by the receptors they use, which determine their host range and tropism. The highly pathogenic ALV-K subgroup has recently been detected in China based solely on its genomic sequence, posing a major threat to Asia's poultry industry14,15. In ALV genomes, subgroup differences can be primarily detected in the env gene, which binds to cell-specific receptors and dictates viral pathogenicity and tropism16.
Despite more than a century of ALV-related research, no vaccine and no specific treatment exist to prevent chickens from the disease. Limited anti-ALV vaccines are being developed. Among them four targeting ALV strain J, one targeting strain B, and another one targeting strain A17. Also, a limited protection and less immunogenicity were observed in some anti-ALV vaccine trials17.However defeating and preventing the ALV disease from spreading has been accomplished by eliminating affected birds, managing hatcheries and farms and implementing biosecurity measures1,4.
With the proliferation of big genomic and proteomic data and computing power, a successful vaccine could be developed by accessing vast information about the virus structure, protein structures, and predicted epitopes. A combination of computational methods and immunoinformatics had formed the bases for the development of multi-epitopes vaccines18,19. In reverse vaccinology, which analyzes a pathogen's genome and identifies its prospective antigenic proteins, multi-epitope vaccines are produced in a shorter period of time and at a lower cost, with safety and highest therapeutic efficiency and minimal adverse reactions18,19. In addition, to the boosting of immunogenicity, vaccines with multi-epitopes may improve immune responses because they contain epitopes representative of different target genes20. Moreover, this novel computational science can restrict the use of animal models and replaces some clinical trials21. Therefore, designing and developing multi-epitope-based peptide vaccines via computational science holds great promise for controlling and preventing the occurrence of infectious diseases in the future. Despite the advances in the field of in silico vaccines formulation, some drawbacks were emerged. For instance, the in silico approaches cannot predict lipids or polysaccharides as they are significant active constituents in proteins22. Also predictions of optimum antigen-processing sites are not well addressed and clarified by many epitope prediction tools23. On the other hand, common traditional vaccines might use inactivated or attenuated forms of microorganism, which might provide adverse effects and autoimmune diseases in infected hosts. Also the development of traditional vaccines candidate such as chimeric virus vaccines, live attenuated vaccines, virus-like particles (VLPs), subunit vaccines and nucleic acid vaccines is complex process, expensive, time consuming, need complicated in vitro and in vivo studies for evaluating their potential efficacy24,25.
Hence, this research targeted the use of immunoinformatics tools to predict multi epitopes vaccine that could efficiently combat the ALV infection. The proteome of ALV was investigated, evaluated and analyzed for epitopes prediction. The results demonstrated that the constructed vaccine has successfully comprised all the most favorable factors to be a potential vaccine candidate, which further requires experimental validation. This study provided good theoretical basis essential for developing a new vaccine strategy and extending the current knowledge in controlling ALV.
Materials and methods
The immunoinformatics steps for the in silico vaccine design were visualized in the flow chart presented in Fig. 1.
ALV protein’s sequences retrieval
The ALV demonstrated three proteins: polymerase protein, envelope protein, and transacting factor protein with the following accession numbers NP_040550.1, NP_040548.1, and NP_040549.1, respectively. The sequences of these three proteins were retrieved from the National Center for Biotechnology Information (NCBI) at (https://www.ncbi.nlm.nih.gov/protein)26.
Physiochemical properties, antigenicity, allergenicity, and toxicity of the retrieved proteins
A physiochemical profile of each retrieved protein was assessed using the Protparam server (http://web.expasy.org/protparam/)27. The potent antigenicity was determined using the VaxiJen v2.0 server (http://www.ddg-pharmfac.net/vaxijen/) based on adefault viral threshold of 0.4 for each antigenic protein28,29. Allergenicity and toxicity of the proteins were assessed using the AllerTOP server (http://www.ddg-pharmfac.net/AllerTop/)30 and the ToxinPred server (http://crdd.osdd.net/raghava/toxinpred/)31, respectively. The same servers were later used to assess the physiochemical properties, antigenicity, allergenicity and toxicity properties of the constructed vaccine.
Subcellular localization and transmembrane topologies of virus proteins
Subcellular localization of viral proteins is considered as an important clue to the function of the immune cells and judging the potential efficacy of vaccine targets32. In addition, surface-localized proteins are among the best candidates for the recombinant vaccine, since they are the first molecular patterns of pathogens that contacted by the host immune system32. For the detection of the viral protein subcellular localization, the Phobius server (https://phobius.sbc.su.se/index.html) was used33. The server provided a combination of a transmembrane topology (TMHs) and a signal peptide predictor.
Epitopes prediction and conservancy
A total of 50, 13, and 3 strains sequences were retrieved for the polymerase, envelope, and transacting factor protein, respectively. These strains were used for epitopes conservancy and were presented in Table 1. BioEdit program version 7.2.5, is a multiple sequence alignment (MSA) tool, was used to align each protein strains sequences34. The analysis of the aligned sequences was conducted in order to identify the conserved epitopes that effectively act against B and T lymphocytes. Epitopes that had 100% conservancy (no mutations) among the strains were selected for further analysis, while non-conserved epitopes were excluded.
B cell epitopes prediction
To identify the B cell epitopes from the three ALV proteins, the ABCpred server (https://webs.iiitd.edu.in/raghava/abcpred/ABC_submission.html) was utilized35,36. The reference sequence of each protein was used as an input into ABCpred server. The default thresholds (0.51), as well as 12mers epitope length were used for the prediction strategy. VaxiJen v2.0, AllerTOP, and ToxinPred servers were further used to assess the antigenic, allergenic and toxic potentiality of the conserved B cell epitopes, respectively.
T cell epitopes prediction
Based on Immune Epitope Database (IEDB) analysis resources at (https://www.iedb.org/), different T cell epitope prediction tools were analyzed37,38. The reference sequence was used as an input for each protein analysis. The data for epitopes that interacted with the major histo-compatibility complex class I and II (MHC-I and MHC-II) is not yet organized in the IEDB resources for chicken alleles. Accordingly, the human alleles were exploited to predict epitopes from the ALV-retrieved proteins interacting with T cell epitope as previously described39,40.
Cytotoxic T cell epitopes prediction
The IEDB prediction method at (http://tools.iedb.org/mhci/) provided a number of MHC-I binding prediction methods. In this study, the prediction of the MHC-I interacted alleles were obtained by Artificial Neural Network, NetMHC (ANN)41. The human reference alleles sets (HLA-A, HLA-B, and HLA-C) were used for the prediction process. Conserved epitopes with a score equal to or less than 1 (≤ 1) percentile rank with nine amino acids length bound to alleles were only analyzed. The conserved cytotoxic T-cell epitopes were further assessed for antigenicity, allergenicity and toxicity predictions.
Helper T cell epitopes prediction
The IEDB MHC-II binding prediction tool (http://tools.iedb.org/mhcii/) was used to investigate the reference sequence of ALV proteins for epitope prediction against MHC-II37,38. The human alleles reference sets (HLA-DP, HLA-DQ, and HLA-DR) were employed to search for promising epitopes. The analysis comprises the Neural Networks Align method, NetMHCII, version 2.2 (NN-align)37,38,41. The approach was used to find potential epitopes having a percentile rank score equal to or less than 10 (≤ 10). The core sequence and peptide lengths were set to 9 and 18 amino acids, respectively. The antigenic, allergenic, and toxic evaluation of the conserved helper T cell epitopes was carried out using VaxiJen v2.0, AllerTOP, and ToxinPred servers, respectively.
Assembly of the multi-epitope vaccine
The primary assembly of the vaccine sequence was accomplished fusing the B and T cells predicted epitopes that demonstrate conservancy, antigenicity score more than 1 and were shown to be non-allergenic and non-toxic. The elected B and T helper epitopes were fused by the GPGPG linkers while the T cytotoxic epitopes were fused by the YAA linkers42. The 5′-amino terminal of the vaccine was supported by the β-defensin 3 (Q5U7J2) as an adjuvant after separation with the EAAAK linker. Moreover, the sequence was provided with the 6His-tag for purification and identification of the vaccine upon expression42,43,44,45.
Secondary and tertiary structures prediction of the vaccine
Predicting the secondary structure of the vaccine is available for free on the raptor X server (http://raptorx.uchicago.edu/)46. The secondary structure (SS), disorder regions (DISO), and solvent accessibility (ACC) played crucial roles in predicting α-helices, β-pleated sheets, and coiled structures, respectively. For the tertiary structure prediction, the vaccine sequence was submitted to the same raptor X server46. The obtained results were received in the form of a PDB file that was further used for vaccine sequence refinement and adaptation.
Vaccine tertiary structure refinement and adaptation
Refinement of the vaccine's tertiary structure was carried out using the GalaxyWEB server (http://galaxy.seoklab.org/)47,48. The structured model is made more dynamic and effective by utilizing the side chains, establishing them and repacking these side chains to achieve relaxation of the structured model47,49. As a result, the refinement and the structure's physical quality were improved. In the SAVES server, the Ramachandran plot (https://saves.mbi.ucla.edu/results?job=869538&p=procheck) was used to analyze the model validation50, and the ProsA server (https://prosa.services.came.sbg.ac.at/prosa.php) was used to analyze the errors in the protein structure51,52.
Determination of the solubility of the vaccine
Protein-sol server (https://protein-sol.manchester.ac.uk/) consists of a set of predictive algorithms and a theoretical computation was used to measure the solubility of the vaccine protein53. The vaccine solubility was assessed in comparison with Escherichia coli proteins, which have an estimated solubility of 0.45. Thus, proteins with solubility score greater than 0.45 are therefore considered soluble53. To confirm the solubility of the vaccine, the SOLpro server (http://scratch.proteomics.ics.uci.edu/) was further employed54. For a protein to be soluble in the SOLpro server the probability should be ≥ 0.5 and insoluble protein scores < 0.5.
Determination of the stability of the vaccine
The disulfide bonding in a given protein between the cysteine residues plays an important role in strengthening of the protein’s geometric conformation and enhances its extensive stability55. Disulfide-by-Design 2.0 (DbD2) (http://cptweb.cpt.wayne.edu/DbD2/) is a web-based tool that facilitates designing disulfide bonds in vaccine construct by substituting particular amino acid with cysteine in high-mobility and unstable regions of proteins55. This was followed by formation of disulfide bonds between cysteine residues. The parameters such as the intra-chain, inter-chain and build Cβ for Gly were chosen. The angle (− 87° or + 97° ±) was set to 30 and Cα–Cβ–Sγ angles (114.6° ±) was set to 10 for proper prediction of the bonds55.
Immune simulation
For mimicry of the immune response and immunogenicity of the ALV vaccine in the host, C-ImmSim server (https://kraken.iac.rm.cnr.it/C-IMMSIM/) was used56. Two injections were given with time step set at 1 and 90 (the server provided each time step as 8 h while the time step 1 represents the injection at time = zero). The other simulation parameters were set to default. The measure of diversity (Simpson index, D) was interpreted from the plot56.
Molecular dynamic simulation (MD)
iMODS server (https://chaconlab.org/multiscale-simulations/imod) was used to analyze the collective motions of protein vaccine57,58. A normal modes analysis (NMA) in internal coordinates is conducted by the server to determine the stability of the vaccine protein. This server structured the dynamics of the protein complex and provided various results data, such as deformability, eigenvalues, B-factors, variance maps, co-variances, elastic networks in the atoms, and residue indexes in terms of magnitude and direction57,58.
Prediction of discontinuous B-cell epitopes
The ElliPro in the IEDB (http://tools.iedb.org/ellipro/) was used to predict the discontinuous B cell epitopes59. ElliPro tool predicts discontinuous and linear antibody epitopes based on the protein 3D structure. The prediction method was based on the default parameters of the sever59. For instance, the minimum score and the maximum distance (Angstrom) of the selected epitopes prediction were set to 0.5 and 6, respectively.
Active sites detection in the vaccine structure
Searching for a ligand-binding region on a protein is an essential step prior to molecular docking process. The process primarily based on multiple factors such as detection of hydrophobic or hydrophilic interactions, salt bridges and electrostatic and hydrogen bonding interactions. The computed atlas of surface topography of proteins (CASTp 3.0) website (http://sts.bioe.uic.edu/castp/index.html?3igg) was exploited to determine the vaccine’s active regions60,61. The default probe radius of 1.4 ˚A was used.
Molecular docking of the vaccine protein with chicken TLR7
Protein–protein interaction is essential for functioning of many biological molecules62. Analyzing the complex structures formed between these molecules is of great importance to assess the molecular interactions or the affinities between these molecules. Toll-like receptors (TLRs) are considered as recognition receptors that play a paramount role in recognition of pathogen. In birds, there are ten genes encoding for TLRs, among them, TLR7 was chosen for the docking with the vaccine construct since it is a viral-sensing TLR63. Thus the designed vaccine was docked against the chicken TLR7 using the HADDOCK 2.4 server (https://www.bonvinlab.org/software/haddock2.4/)62. Refinement interface in HADDOCK server was used to provide the accurate cluster. PRODIGY web server (https://wenmr.science.uu.nl/prodigy/)64,65 was used to calculate the binding affinities of the best chosen clusters at 25 °C. Finally, the interaction between the vaccine and the chicken TLR7 was visualized by PDBsum server (https://www.ebi.ac.uk/thorntonsrv/databases/pdbsum/Generate.html)66.
In silico molecular cloning and codon adaptation
The in silico cloning ensures that a particular host would express the vaccine protein upon cloning in suitable vector67. To facilitate successful cloning, optimization process and cloning of the vaccine construct in the expression vector were performed. The optimization comprises the elimination of different restriction enzymes cleavage sites, prokaryotic ribosomal binding sites, and rho-independent terminators of transcription in the sequence of the vaccine67. A reverse translation of the vaccine protein sequence into a DNA sequence was performed with the Java Codon Adaptation Tool (JCAT) (http://www.prodoric.de/JCat) because cloning uses DNA rather than proteins67. The codon adaptation index and the GC content were in ranges of 0.8–1.0 and 30–70%, respectively. The sequences of the restriction enzymes Xho1 (5-CTCGAG-3) and BamHI (5-GGATCC-3) were added at the 5’and the 3’ ends of the DNA, respectively. A restriction cloning module from SnapGene (https://www.snapgene.com/)67 was used to clone the DNA sequence located between the restriction sites of BamHI and Xho1 in the pET-30a(+) vector.
Results
Characteristics of the virus proteome
Polymerase, envelope, and transacting factor proteins from the ALV were retrieved from the NCBI database. These three proteins were found to be stable and hydrophilic using the ProtParam server. The VaxiJen server was used to determine and prove their antigenicity. The three proteins were used as inputs to predict B and T cell epitopes for designing the vaccine against ALV. All the physical and chemical features of the three proteins were provided in Table 2.
Multiple sequence alignment and epitopes conservancy
The ClustalW program provided in the Bioedit tool was used for multiple sequence alignment (MSA) of all retrieved strains. MSA was exploited to search for conserved epitopes among the retrieved stains from polymerase, envelope, and transacting factor proteins. Epitopes length that was not broken by mutated amino acids from other strains is considered conserved epitope. During the MSA, the retrieved strain sequences demonstrated high epitopes conservancy.
Linear B-cell epitopes prediction
The ABCpred server received the reference sequences from each protein. In the server, a trained recurrent neural network provided the predicted B-cell epitopes based on their scores. Generally, an epitope passing the threshold of 0.51 is more likely to have a higher peptide score. Based on the ABCpred server, 39, 29, and 10 epitopes were predicted from the polymerase, envelope, and transacting factor proteins, respectively. After assessing the antigenicity, allergenicity and the toxicity of the predicted epitopes from each protein, 11, 10 and 6 epitopes from polymerase, envelope and transacting proteins were chosen as B cell epitopes, respectively. These epitopes were provided in Table 3.
Cytotoxic T lymphocyte epitopes prediction
Based on the reference sequences of polymerase, envelope, and transacting factor, multiple epitopes were predicted against human alleles (HLA-A, HLA -B, HLA-C) using IEDB MHC-1 binding prediction tools. Antigenic, allergenic, and toxic effects were then assessed for the predicted epitopes. A total of 6, 11, and 15 epitopes were obtained from the polymerase, envelope, and transacting factor proteins, respectively, and were elected as T cytotoxic cell epitopes due to their high antigenicity scores, non-allergenicity, non-toxicity and the allelic interactions. These epitopes were provided in Table 4.
Helper T lymphocyte epitopes prediction
The reference sequence of each of the three proteins (polymerase, envelope, and transacting factor) was analyzed against the human alleles (HLA-DR, DQ, DP) using IEDB MHC-1I binding prediction tools with a percentile rank of (≤ 10). A vast amount of epitopes were predicted from the three proteins. The predicted epitopes were analyzed for antigenic, allergenic, and toxic outcomes. A total of 21, 6, and 7 epitopes were obtained from the polymerase, envelope, and transacting factor proteins, respectively. They were elected as T helper cell epitopes due to their high antigenicity scores, non-allergenicity, non-toxicity, and allelic interactions. These epitopes were provided in Table 5.
Structure of the assembled vaccine
The entire number of predicted B cell, T cytotoxic, and T helper epitopes from the three proteins of ALV were used in the construction of the vaccine. Adjuvant, linkers, and 6-His-tags were also embedded in the final structure of the vaccine. Thus the final vaccine structure comprised 738 amino acids. The antigenicity score of the assembled vaccine was 0.8535 when examined in the VaxiJen server. Also, the vaccine protein was non-allergic in the AllerTOP server.
Physiochemical properties of the assembled vaccine
ProtParam server was used to examine the physiochemical properties of the assembled vaccine. The predicted vaccine weighed 77.121 kilo Dalton (kd) and possessed a theoretical isoelectric point of 9.81, indicating the proposed vaccine had an alkaline pH. Negatively and positively charged residues were 33 and 79 respectively. The Extinction coefficient at 280 nm measured in water was shown to be 132,125 assuming all pairs of Cys residues forming cystines. The instability index score (II) was 38.24, indicating a stable vaccine protein, while the aliphatic index score was 78.73, indicating a hydrophilic vaccine. The grand average water affinity was -0.130, suggesting a hydrophilic vaccine.
Secondary and tertiary structures prediction of the assembled vaccine
The SS3, ACC, and DISO for the secondary structure were predicted using the Raptor X server. The SS3 showed 23%, 15%, and 61% of the residues as α-helix, β-sheets and coiled, respectively. The ACC provided 49% as exposed residues, 21% as medium residues and 29% as buried residues. The DISO (disordered predicted regions) was 43 (5%). Figure 2 showed the primary sequence, the tertiary and the refined structures of the vaccine construct.
Vaccine tertiary structure refinement and validation
The vaccine's stability was assessed via the Ramachandran plot after refinementt. In the plot, 90.9% of residues were located in the most favored region. While regions of additional allowed, generously allowed, and disallowed comprised residues of 6.1%, 1.9%, and 1.0%, respectively (Fig. 3a). The ProsA server provided a Z score of -5.68 demonstrating a favorable model structure (Fig. 3b).
Solubility of the assembled vaccine
Based on the Protein-Sol server, a scaled solubility score of 0.499 was obtained for the vaccine construct, competed with 0.45 for the population solubility of E. coli (Fig. 4a). As a confirmation, SOLpro was further used to predict the solubility. The probability of the proposed vaccine upon expression on SOLpro was 0.9843, greater than 0.5, provided by the server.
Stability of the assembled vaccine
By engineering disulfide bonds into the structure of the proposed vaccine, the structural stability of the vaccine was improved. The improvement in stability was made possible by substituting the amino acids in the highly mobile regions in the sequence of the vaccine by cysteine residues. As per the Disulfide by Design 2.0 server, 94 amino acid pairs were identified to form disulfide bonds. However based on the Chi3 angle between + 97 and − 87 and a tolerance of 30 and a maximum Ca–Cb–S angle of 114.60 in the server, five pairs of residues (amino acids) were unstable regions and were replaced by cysteine-cysteine residues. The position and the replaced residues in the vaccine structure were A107-R127; I150-G179; P210-P280; P278-P312 and G500-L538 and were shown in Fig. 4b,c.
Immune simulation
The obtained immune simulation results were coincided with actual immune responses. This was proved by marked increase in the primary, secondary and tertiary immune responses accompanied by drop in the antigen concentration (Fig. 5a). The cytokines and interleukins (IL) levels during the injections showed that the IL-2 level was compatible with the measure of diversity (Simpson index, D) (Fig. 5b). The elevation of the measure of diversity over time is considered as danger signal together with leukocyte growth factor. Therefore, the lower the measure of diversity value, the lower the diversity. In addition the primary response, for instance, was featured by augmented IgM level, while, secondary and tertiary responses provided marked elevation in the population of B-cells and the antibodies level (Fig. 5c). This showed the development of immune memory accompanied by rapid clearance of the antigen upon subsequent exposures. Moreover the population of T-cytotoxic (TC) (Fig. 5d) and T-helper (TH) (Fig. 5e) lymphocytes showed high response level coincided with memory development. The natural killer cells maintained high levels throughout the duration of exposure (Fig. 5f).
Molecular dynamic simulation (MD)
A Normal mode analysis (NMA) was performed on the MD of the vaccine protein using the iMODS server. As shown in Fig. 6a, the arrows indicated the direction in which each vaccine protein residue moves. Deformability was also demonstrated with hinges in the chief chain, as a result of an individual distortion of the residues (Fig. 6b). Experimental B-factors were calculated on the basis of the PDB field and the NMA data (Fig. 6c). A normal mode of deformability of the vaccine structure was shown by the eigenvalue, which directly correlated to the energy required with the deformability. The obtained eigenvalue (7.182836e-07) demonstrated the stiffness of the motion (Fig. 6d). The lower eigenvalue is always associated with the easier deformation of the protein structure. The normal mode variance is inversely related to the eigenvalue. Figure 6e illustrated the cumulative variance and individual variance as green and purple bars, respectively. It was possible to determine the correlations between proteins by examining the covariance matrix (Fig. 6f). Thus, red identified correlated motions, white indicated uncorrelated motions, and blue indicated anti-correlated motions. Spring-connected or joined atom pairs were demonstrated in the elastic network model. A single-atom pair spring was represented as a dot, and colored according to its stiffness, with darker dots denoting stiffer strings, and vice versa (Fig. 6g).
Discontinuous B-cell epitopes prediction
Table 6 and Fig. 7 demonstrated six discontinuous B cell epitopes. The scores of these epitopes were ranged from 0.996 to 0.615 with a total of 405 predicted residues. The size of the conformational epitopes ranged from 4 to 108 residues.
Active site detection in the vaccine structure
Prior to docking the active binding pocket of the vaccine construct was predicted by CASTp. The highest pocket area and volume were 7086.109 Å2 and 18,104.055 Å3, respectively. This result indicated pocket area and volume were suitable for binding during the molecular docking. The surface binding pocket of the vaccine and the sequence were shown in Fig. 8.
Molecular docking of the vaccine protein with chicken TLR7
The interaction between the vaccine construct and chicken TLR7 was assessed by HADDOCK software. HADDOCK clustered 13 structures in 3 cluster(s), which represents 6% of the water-refined models. Upon refinement, 20 structures were grouped into one cluster, resulting in 100% of the HADDOCK water-refined version. The binding affinity between the vaccine and the chicken TLR7 was − 263.0 ± 3.1 demonstrating the strong binding between the molecules. As shown in Fig. 9, this binding was evident by 20 hydrogen bonds, 2 salt bridge, and 184 non-bonded contacts. These bonding events between the amino acids of the molecules were provided in Table 7. Additionally PRODIGY web server showed binding affinity in terms of Gibbs free energy (ΔG) and thermodynamics (dissociation constant) between the docked molecules. Such kind of binding affinity decided the real interaction between the docked molecules under certain circumstances within the cell. The server showed ΔG values − 21.1 kcal/mol for the vaccine construct and chicken TLR7 and the dissociation constant was 3.1e−16 indicating the docked molecules were energetically viable.
In silico molecular cloning
The potential host expression of the target protein was performed by JCAT. The protein sequence of the vaccine was reversibly translated into DNA sequence. The index codon adaptation value of the DNA sequence was equal to 1.0, indicating a higher proportional abundance of the common codons. The improved GC content value of the DNA sequence was 59.53%, providing significant GC content. As shown in Fig. 10 the DNA sequence was cloned into the pET-30a(+) vector utilizing Snapgene software between the restriction enzyme cutting locations of BamH1 and Xho1.
Discussion
The most common avian retrovirus that causes a variety of neoplastic diseases in chicken is the avian leukosis virus (ALV)2. Globally, the ALV morbidity and mortality rates contributed to the poultry industry’s economic decline1,3. This is accompanied by subsequent adverse effects on the food supply worldwide. Preventing and controlling viral infection in avian industry is always via mass vaccination means. Therefore vaccines designed to combat avian viral diseases will significantly alleviate selection pressure on the virus and on the field strains68. Concerning ALV infection in poultry, many anti-ALV vaccines were developed, but they targeted only specific strains. Also some of the vaccine trials had less immunogenicity and limited protection69. Currently, neither known treatment nor vaccination against ALV is available. Multiple studies used the multi-epitope vaccine prediction against ALV and evaluated their possibility as effective vaccine candidate via challenging in chickens70,71,72. For instance, one study provided a novel oral vaccine of recombinant gp85 protein in L. plantarum with a significant increase in antibodies post-inoculation72. The study demonstrated a protection against ALV-J and showed protective immune response against early ALV-J infection based on viremia analysis73. Another study showed the impact of polysaccharides from Ulvapertusa as anti-ALV-J. The polysaccharides demonstrated strongest suppression of the ALV-J activity as they bound with the viral particles and obstacle ALV-J adsorption by host cells accompanied by significant reduction of gp85 protein expression74. However these studies reported partial immune protections against ALV-J infections in chickens.
In this study a vaccine with multi-epitopes was designed and showed increased immunogenicity and enhanced immune responses as a result of the existence of epitopes from various target genes. Also the designed vaccine activated the humoral and cell-mediated immunity as previously described75. These could solve limitations occurred during controlling the ALV infection75,76. Most importantly, the safety and effectiveness, allergenicity and the immunogenicity of the predicted vaccine were also taken into consideration to ensure the safety of the designed epitopes77. In addition, the toxic effect, the solvent accessibility of the amino acids, the identification of B cells, and MHCmolecules were also contemplated to ensure the effectiveness of the predicted epitope vaccine. All these measures give the predicted vaccine an advantage over the traditional ones for controlling the ALV infection.
Thus the conserved predicted epitopes from ALV proteins were submitted to the ABCpred server. Based on ANN, Hidden Markov model (HMM) and support vector machine (SVM) in the ABCpred server the B cell epitopes were predicted78. Furthermore, the predicted epitopes were subjected to antigenic, allergenic and toxic analysis to confirm their suitability as B cell epitopes. Also T cell epitopes were predicted from their reference sequences using the IEDB server. In addition to their high binding affinity to MHC alleles, the predicted epitopes demonstrated high antigenicity score in VaxiJen server, and they revealed no allergic or toxic characteristics. Therefore they were picked to enter the vaccine protein structure. With the aid of expedient linker sequences (protein spacers), the generated B- and T-cell epitopes were fused together49,52. Linkers are crucial to the assembly of stable, bioactive fusion proteins. Essentially, linkers reduced the likelihood of junctional antigen formation as well as enhancing antigen processing and presentation52. They are also important to construct and facilitate structural flexibility and reduced rigidity52. A sequence with the least junctional immunogenicity was generated in this study using the linkers GPGPG and YAA. The GPGPG linkers were applied to facilitate immune processing and merge the B-cells and T-helper cell epitopes. The YAA linkers ameliorated the immunogenicity of a vaccine by impacting protein stability and epitope presentation capacity and were used to fuse the cytotoxic T-cells24,79. As an adjuvant, the β-defensin was added via an EAAAK linker at the N terminus of the vaccine construct to improve the immunogenicity of the vaccine. EAAAK are helical linkers used to control the distance and decrease the interference between the domains24,79. As a 45amino acids peptide with a relatively small size, the β-defensin was used for its immune modulation and antimicrobial features44. To facilitate purification and downstream testing, a small 6His tag was added to the proposed vaccine at the C-terminal to prevent protein structure from being altered80.
The stability of the vaccine was confirmed by the ProtParam server based on its physiochemical properties. VaxiJen and AllerTOP servers were used to assess the antigenic and allergenic features of the vaccine. The results indicated that the vaccine was antigenic without causing any allergic reactions. In order to select the best score of the model generated by the 3D structure of the vaccine protein, the secondary and tertiary structures of the vaccine construct were analyzed. The Ramachandran plot showed favorable results in the distribution of the vaccine residues and provided a stable structure. The ProSA server indicating that the overall model is suitable for acceptance as a potential ALV vaccine51,52.
The solubility of the designed vaccine in this study was calculated with the protein-sol and SOLpro servers. As a comparison with the solubility of E coli, Protein-sol presented the vaccine as a soluble protein and predicted a scaled solubility of 0.499, an increase over 0.45 from the average solubility of the E. coli population. According to the SOLpro server, the predicted solubility was 0.9843, which confirmed this result. To obtain disulfide bonds between the vaccine residues, the proximity and geometry composition of the residue pairs were evaluated for the formation of disulfide bonds. Five unstable regions in the vaccine structure were replaced by the formation of disulfide bonds. Disulfide bonding increases the stability of the vaccine protein as previously stated51,52.
Immune simulation demonstrated results that consistent with the real immune responses. Generally there were elevated levels of the immune responses after repeated exposure to the vaccine (antigen). In addition, there was marked development in the memory cells of B and T lymphocytes. Most importantly, IL-2 and IFN-γ were elevated following the initial injection and provided peak levels after antigen repeated exposures, showing the high levels of T-helper lymphocytes and efficient immunoglobulin production. The Simpson index demonstrated a possible different immune response, indicating the vaccine structure contains multiple B and T cells epitopes44. A study by Landman et al., demonstrated the interaction of the NK cells during ALV infection81. They showed that during ALV infection in immunosuppressed chicken, the NK cells provided reduced killing activity than the NK cells of the uninfected controls. Natural killer cells play a paramount defense mechanism in host and surveillance of tumor, resulting in cell death and secretion of cytokines and chemokine. Moreover, NK cells have a significant role in immune regulation of T cells and DC functions during viral infection in mouse models81. In addition to that, there is scarcity in ALV vaccine researches concerning the immune system of chickens. Thus vaccine research studies may provide a realistic importance and increase the knowledge of immune system of ALV with great insights into human retroviral diseases69.
Molecular dynamics simulation (SD) was used to assess the complex stability of the vaccine protein. In previous studies, macromolecule stability was associated and correlated with the fluctuations of atoms82,83. Therefore MD was performed to evaluate the essential dynamics and complex stability of the vaccine based on the protein normal modes in the iMODS server. The analysis showed that no atoms had a significant distortion in the vaccine protein structure indicating less chance of deformability with proper stiffness motion.
It is noteworthy that bioinformatics and immunologic analysis tools provided that the chimeric vaccine should comprises linear and discontinuous B-cell epitopes in addition to MHC-I and MHC-II epitopes84. Our predicted vaccine was shown comprising all these epitopes which strongly facilitate the interaction against the humoral and adaptive immunity of the host84.
The geometry and topology features of protein structures, such as interior cavities, pockets in the structure surface and the cross channels prior to the docking process are essential to study the function of proteins. The vaccine construct showed surface binding pocket suitable for docking with chicken TLR7. Based on the molecular docking, the constructed vaccine and the TLR7 demonstrated a good binding affinity. The vaccine strongly bound to the chicken TLR7 revealed by the negative values of the docking process62. Among the ten chicken TLRs, TLR7 has a propensity to recognize the viral constituents located on the extracellular surfaces63, thus has the advantages to be elected for docking against ALV predicted vaccine.
Molecular cloning is an important step to produce recombinant vaccines. Prior to cloning into the pET-30a [+] vector, reverse transcription and adaptation of the vaccine protein to DNA by the JCAT were performed on E. coli strain K12. The DNA sequence showed a Cal-Value of 1.00 and a GC ratio of 59.57%, demonstrating a high expression in bacteria. Cloning of the vaccine construct gene in the vector was typically carried out in multiple cloning sites. This result provided prolific cloning of the vaccine protein.
Conclusion
This study demonstrated the urgency need for effective vaccine strategy against ALV due to the lack of treatment or approved antiviral drugs. This study inclusively exploited the computational and immunoinformatics approaches to design and evaluate a multi-epitope vaccine candidate against ALV. Constructing abjunctive vaccine with antigenic characteristics, devoid of allergenicity and toxicity is a crucial footstep to combat ALV. This study provided a potential vaccine epitopes with immunogenic adjuvant and suitable linkers. The vaccine was stable and provokes strong immune response interactions. Moreover the vaccine showed favorable interaction with the chicken immune receptor as confirmed by molecular docking analysis. However, validation of this vaccine via experimental studies is essential to guarantee the immunogenicity and protective efficacy of the vaccine.
Data availability
The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.
References
Fadly, A. M. Isolation and identification of avian leukosis viruses: A review. Avian Pathol.: J. W.V.P.A 29(6), 529–535. https://doi.org/10.1080/03079450020016760 (2000).
Fandiño, S., Gomez-Lucia, E., Benítez, L. & Doménech, A. Avian leukosis: Will we be able to get rid of it?. Anim.: Open Access J. MDPI 13(14), 2358. https://doi.org/10.3390/ani13142358 (2023).
Zheng, L. P. et al. Current epidemiology and co-infections of avian immunosuppressive and neoplastic diseases in chicken flocks in central China. Viruses 14(12), 2599. https://doi.org/10.3390/v14122599 (2022).
Borodin, A. M. et al. Eradication of avian leukosis virus subgroups J and K in broiler cross chickens by selection against infected birds using multilocus PCR. PloS One 17(6), e0269525. https://doi.org/10.1371/journal.pone.0269525 (2022).
Zavala, G. & Cheng, S. Detection and characterization of avian leukosis virus in Marek’s disease vaccines. Avian Dis. 50(2), 209–215. https://doi.org/10.1637/7444-092405R.1 (2006).
Gao, Y. L. et al. Avian leukosis virus subgroup J in layer chickens, China. Emerg. Infect. Dis. 16(10), 1637–1638. https://doi.org/10.3201/eid1610.100780 (2010).
Khan, A. et al. Lymphoid leukosis in fayoumi birds reared in countryside. Intl. J. Agric. Biol. 22, 620–626. https://doi.org/10.17957/IJAB/15.1107 (2019).
Li, H. et al. Aberrant expression of liver microRNA in chickens infected with subgroup J avian leukosis virus. Virus Res. 169(1), 268–271. https://doi.org/10.1016/j.virusres.2012.07.003 (2012).
Payne, L. N. & Nair, V. The long view: 40 years of avian leukosis research. Avian Pathol: J. W.V.P.A 41(1), 11–19. https://doi.org/10.1080/03079457.2011.646237 (2012).
Abdel-Latif, M. M. & Khalafalla, A. I. Detection by PCR of multiple subgroups of avian leukosis virus (ALV) in broilers in the Sudan. J. Anim. Vet. Adv. 4(3), 407–413 (2005).
Li, Q. et al. Recombinant subgroup B avian leukosis virus combined with the subgroup J env gene significantly increases its pathogenicity. Vet. Microbiol. 250, 108862. https://doi.org/10.1016/j.vetmic.2020.108862 (2020).
Cui, N., Su, S., Chen, Z., Zhao, X. & Cui, Z. Genomic sequence analysis and biological characteristics of a rescued clone of avian leukosis virus strain JS11C1, isolated from indigenous chickens. J. Gen. Virol. 95(Pt 11), 2512–2522. https://doi.org/10.1099/vir.0.067264-0 (2014).
Parvizpour, S., Pourseif, M. M., Razmara, J., Rafi, M. A. & Omidi, Y. Epitope-based vaccine design: A comprehensive overview of bioinformatics approaches. Drug Discov. Today 25(6), 1034–1042. https://doi.org/10.1016/j.drudis.2020.03.006 (2020).
Li, H. et al. Diversity of avian leukosis virus subgroup J in local chickens, Jiangxi, China. Sci. Rep. 11(1), 4797. https://doi.org/10.1038/s41598-021-84189-7 (2021).
Kheimar, A. et al. A genetically engineered commercial chicken line is resistant to highly pathogenic avian leukosis virus subgroup. J. Microorgan. 9(5), 1066. https://doi.org/10.3390/microorganisms9051066 (2021).
Payne, L. N. et al. A novel subgroup of exogenous avian leukosis virus in chickens. J. Gen. Virol. 72(Pt 4), 801–807. https://doi.org/10.1099/0022-1317-72-4-801 (1991).
Wang, X., Zhao, P. & Cui, Z. Z. Identification of a new subgroup of avian leukosis virus isolated from Chinese indigenous chicken breeds]. Bing Du XueBao 28(6), 609–614 (2012) (Chinese).
Patronov, A. & Doytchinova, I. T-cell epitope vaccine design by immunoinformatics. Open Biol. 3(1), 120139. https://doi.org/10.1098/rsob.120139 (2013).
Vij, S., Thakur, R. & Rishi, P. Reverse engineering approach: A step towards a new era of vaccinology with special reference to Salmonella. Exp. Rev. Vacc. 21(12), 1763–1785. https://doi.org/10.1080/14760584.2022.2148661 (2022).
Zhang, J. et al. The immunogenicity and protective immunity of multi-epitopes DNA prime-protein boost vaccines encoding Amastin-Kmp-11, Kmp11-Gp63 and Amastin-Gp63 against visceral leishmaniasis. PloS One 15(3), e0230381. https://doi.org/10.1371/journal.pone.0230381 (2020).
Saldanha, L., Langel, Ü. & Vale, N. In Silico studies to support vaccine development. Pharmaceutics 15(2), 654. https://doi.org/10.3390/pharmaceutics15020654 (2023).
Jorge, S. & Dellagostin, O. A. The development of veterinary vaccines: A review of traditional methods and modern biotechnology approaches. Biotechnol. Res. Innov. 1, 6–13. https://doi.org/10.1016/j.biori.2017.10.001 (2017).
Silva-Arrieta, S., Goulder, P. J. R. & Brander, C. In silico veritas? Potential limitations for SARS-CoV-2 vaccine development based on T-cell epitope prediction. PLoS Pathog. 16(6), e1008607. https://doi.org/10.1371/journal.ppat.1008607 (2020).
Shawan, M. M. A. K. et al. Advances in computational and bioinformatics tools and databases for designing and developing a multi-epitope-based peptide vaccine. Int. J. Pept. Res. Ther. 29(4), 60. https://doi.org/10.1007/s10989-023-10535-0 (2023).
Dey, J. et al. Designing a novel multi-epitope vaccine to evoke a robust immune response against pathogenic multidrug-resistant Enterococcus faecium bacterium. Gut Pathog. 14(1), 21. https://doi.org/10.1186/s13099-022-00495-z (2022).
NCBI Resource Coordinators. Database resources of the National Center for Biotechnology Information. Nucl. Acids Res. 44(D1), D7–D19. https://doi.org/10.1093/nar/gkv1290 (2016).
Gasteiger, E., Hoogland, C., Gattiker, A., Wilkins, M. R., Appel, R.D. & Bairoch, A. Protein Identification and Analysis Tools on the ExPASy Server. The Proteomics Protocols Handbook, pp. 571–607 (2005). https://doi.org/10.1385/1-59259-890-0:571
Doytchinova, I. A. & Flower, D. R. VaxiJen: A server for prediction of protective antigens, tumour antigens and subunit vaccines. BMC Bioinform. 8, 4. https://doi.org/10.1186/1471-2105-8-4 (2007).
Zaharieva, N., Dimitrov, I., Flower, D. & Doytchinova, I. Immunogenicity prediction by VaxiJen: A 10 year overview. J. Proteom. Bioinform. https://doi.org/10.4172/jpb.1000454 (2017).
Dimitrov, I., Bangov, I., Flower, D. R. & Doytchinova, I. AllerTOP vol 2–a server for prediction of allergens. J. Mol. Model. 20(6), 2278. https://doi.org/10.1007/s00894-014-2278-5 (2014).
Gupta, S., Open Source Drug Discovery Consortium et al. Approach for predicting toxicity of peptides and proteins. PloS One 8(9), e73957. https://doi.org/10.1371/journal.pone.0073957 (2013).
Mao, Y., Su, Q., Li, J., Jiang, T. & Wang, Y. Avian leukosis virus contamination in live vaccines: A retrospective investigation in China. Vet. Microbiol. 246, 108712. https://doi.org/10.1016/j.vetmic.2020.108712 (2020).
Käll, L., Krogh, A. & Sonnhammer, E. L. Advantages of combined transmembrane topology and signal peptide prediction–the Phobius web server. Nucl. Acids Res. 35(Web Server issue), W429–W432. https://doi.org/10.1093/nar/gkm256 (2007).
Hall, T. A. BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. in Nucleic Acids Symposium Series, Jan 1 (Vol. 41, No. 41, pp. 95–98). Information Retrieval Ltd., c1979-c2000 (1999)
Ansari, H. R. & Raghava, G. P. Identification of conformational B-cell Epitopes in an antigen from its primary sequence. Immun. Res. 6, 6. https://doi.org/10.1186/1745-7580-6-6 (2010).
Saha, S. & Raghava, G. P. Prediction of continuous B-cell epitopes in an antigen using recurrent neural network. Proteins 65(1), 40–48. https://doi.org/10.1002/prot.21078 (2006).
Vita, R. et al. The immune epitope database (IEDB): 2018 update. Nucl. Acids Res. 47(D1), D339–D343. https://doi.org/10.1093/nar/gky1006 (2019).
Dhanda, S. K. et al. IEDB-AR: Immune epitope database-analysis resource in 2019. Nucl. Acids Res. 47(W1), W502–W506. https://doi.org/10.1093/nar/gkz452 (2019).
Mugunthan, S. P. & Harish, M. C. Multi-epitope-based vaccine designed by targeting cytoadherence proteins of mycoplasma gallisepticum. ACS Omega 6(21), 13742–13755. https://doi.org/10.1021/acsomega.1c01032 (2021).
Andongma, B. T. et al. design of a promiscuous chimeric multi-epitope vaccine against Mycobacterium tuberculosis. Comput. Struct. Biotechnol. J. 21, 991–1004. https://doi.org/10.1016/j.csbj.2023.01.019 (2023).
Nielsen, M. & Lund, O. NN-align An artificial neural network-based alignment algorithm for MHC class II peptide binding prediction. BMC Bioinform. 10, 296. https://doi.org/10.1186/1471-2105-10-296 (2009).
Omoniyi, A. A. et al. design and analyses of a multi-epitope vaccine against Crimean-Congo hemorrhagic fever virus through reverse vaccinology and immunoinformatics approaches. Sci. Rep. 12(1), 8736. https://doi.org/10.1038/s41598-022-12651-1 (2022).
Pandey, R. K., Ojha, R., Aathmanathan, V. S., Krishnan, M. & Prajapati, V. K. Immunoinformatics approaches to design a novel multi-epitope subunit vaccine against HIV infection. Vaccine 36(17), 2262–2272. https://doi.org/10.1016/j.vaccine.2018.03.042 (2018).
Almofti, Y. A., Abd-Elrahman, K. A. & Eltilib, E. E. M. Vaccinomic approach for novel multi epitopes vaccine against severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2). BMC Immunol. 22(1), 22. https://doi.org/10.1186/s12865-021-00412-0 (2021).
Mahmoud, N. A., Elshafei, A. M. & Almofti, Y. A. A novel strategy for developing vaccine candidate against Jaagsiekte sheep retrovirus from the envelope and gag proteins: An in-silico approach. BMC Vet. Res. 18(1), 343. https://doi.org/10.1186/s12917-022-03431-0 (2022).
Wang, S., Li, W., Liu, S. & Xu, J. RaptorX-Property: A web server for protein structure property prediction. Nucl. Acids Res. 44(W1), W430–W435. https://doi.org/10.1093/nar/gkw306 (2016).
Ko, J., Park, H., Heo, L. & Seok, C. GalaxyWEB server for protein structure prediction and refinement. Nucl. Acids Res. 40(Web Server issue), W294–W297. https://doi.org/10.1093/nar/gks493 (2012).
Seok, C. et al. Accurate protein structure prediction: What comes next?. Bio Des. 9, 47–50. https://doi.org/10.34184/kssb.2021.9.3.47 (2021).
Srivastava, S. et al. Structural Basis for designing multiepitope vaccines against COVID-19 infection: Vaccine design and validation. JMIR Bioinform. Biotechnol. 1(1), e19371. https://doi.org/10.2196/19371 (2020).
Lovell, S. C. et al. Structure validation by Calpha geometry: phi, psi and Cbeta deviation. Proteins 50(3), 437–450. https://doi.org/10.1002/prot.10286 (2003).
Wiederstein, M. & Sippl, M. ProSA-web: Interactive web service for the recognition of errors in three-dimensional structures of proteins. Nucl. Acids Res. 35, W407–W410. https://doi.org/10.1093/nar/gkm290 (2007).
Shawan, M. M. A. K. et al. Insilico modeling and immunoinformatics probing disclose the epitope based peptide vaccine against Zika virus envelope glycoprotein. Indian J. Pharm. Biol. Res. 2, 44–57. https://doi.org/10.30750/ijpbr.2.4.10 (2014).
Hebditch, M., Carballo-Amador, M. A., Charonis, S., Curtis, R. & Warwicker, J. Protein-Sol: A web tool for predicting protein solubility from sequence. Bioinformatics (Oxford, England) 33(19), 3098–3100. https://doi.org/10.1093/bioinformatics/btx345 (2017).
Magnan, C. N., Randall, A. & Baldi, P. SOLpro: Accurate sequence-based prediction of protein solubility. Bioinformatics (Oxford, England) 25(17), 2200–2207. https://doi.org/10.1093/bioinformatics/btp386 (2009).
Craig, D. B. & Dombkowski, A. A. Disulfide by design 2.0: A web-based tool for disulfide engineering in proteins. BMC Bioinform. 14, 346. https://doi.org/10.1186/1471-2105-14-346 (2013).
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 5(4), e9862. https://doi.org/10.1371/journal.pone.0009862 (2010).
López-Blanco, J. R., Aliaga, J. I., Quintana-Ortí, E. S. & Chacón, P. iMODS: Internal coordinates normal mode analysis server. Nucl. Acids Res. 42(Web Server issue), W271–W276. https://doi.org/10.1093/nar/gku339 (2014).
Lopéz-Blanco, J. R., Garzón, J. I. & Chacón, P. iMod: multipurpose normal mode analysis in internal coordinates. Bioinformatics (Oxford, England) 27(20), 2843–2850. https://doi.org/10.1093/bioinformatics/btr497 (2011).
Ponomarenko, J. et al. ElliPro: A new structure-based tool for the prediction of antibody epitopes. BMC Bioinform. 9, 514. https://doi.org/10.1186/1471-2105-9-514 (2008).
Tian, W., Chen, C., Lei, X., Zhao, J. & Liang, J. CASTp 3.0: Computed atlas of surface topography of proteins. Nucl. Acids Res. 46(W1), W363–W367. https://doi.org/10.1093/nar/gky473 (2018).
Halder, S. K. et al. In silico identification and analysis of potentially bioactive antiviral phytochemicals against SARS-CoV-2: A molecular docking and dynamics simulation approach. BioMed Res. Int. 2023, 5469258. https://doi.org/10.1155/2023/5469258 (2023).
Honorato, R. V. et al. Structural biology in the clouds: The WeNMR-EOSC ecosystem. Front. Mol. Biosci. 8, 729513. https://doi.org/10.3389/fmolb.2021.729513 (2021).
Vinkler, M., Bainová, H. & Bryja, J. Protein evolution of Toll-like receptors 4, 5 and 7 within Galloanserae birds. Genet. Select. Evolut.: GSE 46(1), 72. https://doi.org/10.1186/s12711-014-0072-6 (2014).
Vangone, A. & Bonvin, A. M. Contacts-based prediction of binding affinity in protein-protein complexes. eLife 4, e07454. https://doi.org/10.7554/eLife.07454 (2015).
Xue, L. C., Rodrigues, J. P., Kastritis, P. L., Bonvin, A. M. & Vangone, A. PRODIGY: A web server for predicting the binding affinity of protein-protein complexes. Bioinformatics (Oxford, England) 32(23), 3676–3678. https://doi.org/10.1093/bioinformatics/btw514 (2016).
Laskowski, R. A., Jabłońska, J., Pravda, L., Vařeková, R. S. & Thornton, J. M. PDBsum: Structural summaries of PDB entries. Prot. Sci.: Publ. Prot. Soc. 27(1), 129–134. https://doi.org/10.1002/pro.3289 (2018).
Grote, A. et al. JCat: A novel tool to adapt codon usage of a target gene to its potential expression host. Nucl. Acids Res. 33(Web Server issue), W526–W531. https://doi.org/10.1093/nar/gki376 (2005).
Ravikumar, R., Chan, J. & Prabakaran, M. Vaccines against major poultry viral diseases: Strategies to Improve the breadth and protective efficacy. Viruses 14(6), 1195. https://doi.org/10.3390/v14061195 (2022).
Feng, M. & Zhang, X. Immunity to Avian Leukosis Virus: Where are we now and what should we do?. Front. Immunol. 7, 624. https://doi.org/10.3389/fimmu.2016.00624 (2016).
Xu, Q., Ma, X., Wang, F., Li, H. & Zhao, X. Evaluation of a multi-epitope subunit vaccine against avian leukosis virus subgroup J in chickens. Virus Res. 210, 62–68. https://doi.org/10.1016/j.virusres.2015.06.024 (2015).
Xu, Q. et al. Design and construction of a chimeric multi-epitope gene as an epitope-vaccine strategy against ALV-J. Prot. Expr. Purif. 106, 18–24. https://doi.org/10.1016/j.pep.2014.10.007 (2015).
Wang, X. et al. A novel multi-variant epitope ensemble vaccine against avian leukosis virus subgroup. J. Vacc. 35(48 Pt B), 6685–6690. https://doi.org/10.1016/j.vaccine.2017.10.019 (2017).
Wang, S. et al. Oral immunization of chickens with recombinant Lactobacillus plantarum vaccine against early ALV-J infection. Front. Immunol. 10, 2299. https://doi.org/10.3389/fimmu.2019.02299 (2019).
Sun, Y. et al. Antiviral activity against avian leucosis virus subgroup J of degraded polysaccharides from Ulva pertusa. BioMed Res. Int. 9415965, 2018. https://doi.org/10.1155/2018/9415965 (2018).
Barman, A., Deb, B. & Chakraborty, S. Prediction of potential epitopes for peptide vaccine formulation against teschovirus a using immunoinformatics. Int. J. Pept. Res. Ther. 26(2), 1137–1146. https://doi.org/10.1007/s10989-019-09916-1 (2020).
Kar, T. et al. A candidate multi-epitope vaccine against SARS-CoV-2. Sci. Rep. 10(1), 10895. https://doi.org/10.1038/s41598-020-67749-1 (2020).
Zepp, F. Principles of vaccine design-Lessons from nature. Vaccine 28(Suppl 3), C14–C24. https://doi.org/10.1016/j.vaccine.2010.07.020 (2010).
Ojha, R., Pareek, A., Pandey, R. K., Prusty, D. & Prajapati, V. K. Strategic development of a next-generation multi-epitope vaccine to prevent Nipah virus zoonotic infection. ACS Omega 4(8), 13069–13079. https://doi.org/10.1021/acsomega.9b00944 (2019).
Nain, Z. et al. Proteome-wide screening for designing a multi-epitope vaccine against emerging pathogen Elizabethkingia anophelis using immunoinformatic approaches. J. Biomol. Struct. Dyn. 38(16), 4850–4867. https://doi.org/10.1080/07391102.2019.1692072 (2020).
Fan, J. et al. Engineering his-tagged senecavirus a for one-step purification of viral antigens. Vaccines 10(2), 170. https://doi.org/10.3390/vaccines10020170 (2022).
Landman, W. J. et al. Effect of an in ovo infection with a Dutch avian leukosis virus subgroup J isolate on the growth and immunological performance of SPF broiler chickens. Avian Pathol: J. W.V.P.A 31(1), 59–72. https://doi.org/10.1080/03079450120106633 (2002).
Cheatham, T. I., Miller, J. L., Fox, T., Darden, T. A. & Kollman, P. A. Molecular dynamics simulations on solvated biomolecular systems: the particle mesh Ewald method leads to stable trajectories of DNA, RNA, and proteins. J. Am. Chem. Soc. 117(14), 4193–4194. https://doi.org/10.1021/ja00119a045 (1995).
Tehei, M. & Zaccai, G. Adaptation to high temperatures through macromolecular dynamics by neutron scattering. FEBS J. 274(16), 4034–4043. https://doi.org/10.1111/j.1742-4658.2007.05953.x (2007).
Ali, M. et al. Exploring dengue genome to construct a multi-epitope based subunit vaccine by utilizing immunoinformatics approach to battle against dengue infection. Sci. Rep. 7(1), 9232. https://doi.org/10.1038/s41598-017-09199-w (2017).
Acknowledgements
Special acknowledge to Lamees Tarig Mohamed for writing assistance
Funding
This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors Ethics declarations.
Author information
Authors and Affiliations
Contributions
S.O.E.: conceived the study, designed the study and wrote and drafted the manuscript. N.A.M.: conceived the study, designed the study, did data analysis, and wrote and drafted the manuscript. Y.A.A.: Conceived the study and data analysis, coordinated the study, and drafted the manuscript. All authors gave final approval for publication.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Elshafei, S.O., Mahmoud, N.A. & Almofti, Y.A. Immunoinformatics, molecular docking and dynamics simulation approaches unveil a multi epitope-based potent peptide vaccine candidate against avian leukosis virus. Sci Rep 14, 2870 (2024). https://doi.org/10.1038/s41598-024-53048-6
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-024-53048-6
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.