Skip to main content

Insilco vaccine design of spike and hemagglutinin esterase proteins of bovine coronavirus

Abstract

Background

Bovine coronavirus (BCoV) is a widespread, fatal illness in cattle that has a large economic impact in particularly disease-prone hosts. BCoV does not have an effective vaccine. As a result, the objective of this study was to use immunoinformatics and computational tools to design a multi-epitope vaccine against Spike and haemagglutinin proteins of Bovine coronavirus. BCoV proteins were also subjected to protein analysis.

Methods

A variety of tools of bioinformatics were used for data analysis. Conserved B and T cell epitopes against BCoV target proteins were predicted using the Immune Epitope Database (IEDB). Proteins were characterized utilizing a variety of servers, including Protparam, PSIPRED and the GOR IV servers. The antigenicity, allergenicity, and toxicity of the anticipated epitopes were assessed as well.

Results

Several MHCI epitopes were predicted from S and HE proteins. As top epitopes, the peptides 77NMALKGTLL85 and 56SYMDLNPAL65 were proposed from Spike and hemagglutinin proteins, respectively. These epitopes exhibited high scores of antigenicity, no allergenicity, no toxicity, and a strong connection to Bola alleles. Moreover, three epitopes (1204YYYPE1208, 379TCQPQ384, and 720QLQPINY726) from Spike glycoprotein were selected as surface, linear, and antigenic epitopes using B cell scales. The methods dropped to anticipate effective and safe epitope(s) to cover all B cell scales from HE protein.

Conclusion

Three B cell epitopes (1204YYYPE1208, 379TCQPQ384, and 720QLQPINY726) were predicted from Spike protein (S) of BCoV only. MHC I epitopes of S and HE proteins of BCoV predicted two epitopes (77NMALKGTLL85 and 56SYMDLNPAL65 respectively) to have a strong link to Bola alleles, as well as high antigenicity and safety. The predicted epitopes’ activity should be tested experimentally as a multi-epitope vaccine against BCoV using in vitro and in vivo trials.

Introduction

Bovine respiratory disease (BRD) is caused by pathogenic organisms and is typically exacerbated by stress factors such as the environment, nutrition, transportation, and commingling with cattle from various herds [1]. Bovine coronavirus (BCoV) is now generally recognized as a major cause of neonatal calf diarrhea [2]. The disease affects the respiratory tract of cattle and has been linked with different syndromes in cattle: calf diarrhea, winter dysentery in adult cattle, and respiratory infection in cattle of different age groups [3,4,5,6,7,8]. BCoV can significantly economically impact the veterinary industry [9]. The bovine coronavirus (BCoV) is a member of a family of microorganisms that causes intestinal and respiratory disease in a variety of mammalian and avian species. BCoV’s function in calf hood diarrhea is well-established, and it remains a problem in calf-rearing enterprises [10]. Along with wild ruminant CoVs, porcine hemagglutinating encephalomyelitis virus, equine coronavirus, HCoV-OC43, HECoV-44, and canine respiratory coronavirus, bovine coronavirus (BCoV) belongs to the species Betacoronavirus 1 (subgenus Embecovirus) of the Betacoronavirus genus. In the Bovine Respiratory Disease Complex, BCoV is a major pathogen. Viral particles are big (100–150 nm), pleomorphic, and contain four key structural proteins: a membrane (M) glycoprotein, an envelope (E) glycoprotein, a spike (S) glycoprotein, and a hemagglutinin-esterase (HE) glycoprotein [5, 11, 12].

Bovine coronavirus BCoV is pleomorphic and enclosed, with a diameter varying from 65 to 210 nm with a double layer of short (hemagglutinin) and long (spike) surface projections. The big genome is made up of 27 to 32 kb of single-stranded positive-sense RNA that codes for five main structural proteins. The 50-kDa nucleocapsid (N) is the most conserved among them, making it a popular target for viral RNA detection techniques. 29 The presence of a surface hemagglutinin-esterase (HE) glycoprotein (120–140 kDa) is unique to some group 2 CoVs, including BCoV and wild ruminant CoVs. To reverse hemagglutination, the HE operates as a receptor-destroying enzyme (esterase) [11, 12].

Bovine coronavirus (BCoV), like other coronavirus (CoVs), has an outer-surface S glycoprotein (190 kDa). It is made up of two subunits, one of which carries the dominant neutralizing epitopes and the other of which mediates viral membrane fusion. Because HE and S proteins elicit neutralizing antibodies that can prevent viral attachment and infectivity, they are significant for immunity and vaccine development [13, 14].

Bovine coronavirus (BCoV) respiratory illness has little control methods. The vaccinations for BCoV are approved for the prevention of newborn intestinal illness [9, 15, 16]. There are three inactivated vaccines approved for the treatment of neonatal enteric disease, and they are given to pregnant cows/heifers throughout pregnancy to increase humoral immunity in the newborn calf. A modified live virus vaccination against BCoV, the vaccine was found to be safe as four colostrum-deprived newborn calves remained healthy after oronasal administration of ten doses of the vaccine that can induce an active antibody and cell-mediated immunity at the same time [9, 15].

Three antigenic groups of coronaviruses have been identified, and all BCoV strains classified around the world belonged to the 2a subgroup [9]. The International Commission for Virus Taxonomy (ICTV) has recommended a revision of the Coronaviridae family to include the Alpha, Beta, and Gammacorona-virus genera in a new subfamily Coronavirinae. BCoV is classified as a member of the Betacoronavirus genus, which is part of the Corona-virinae subfamily, Coronaviridae family, and Nidovirales order [17,18,19,20].

Recently, prevention and control of the infectious diseases have been among the top public health priorities. However, controlling disease due to pathogens that move between animals and humans has been challenging. Such zoonotic pathogens have been responsible for the majority of new human disease threats and a number of recent international epidemics [21].

Developing a medication or vaccine using traditional methods was a labor-intensive and time-consuming process that also carried the risk of toxicity and safety. By combining computational and experimental approaches, significant progress in building insilico tools not only accelerates drug discovery and vaccine design, but also decreases cost and time [22]. There are various possible advantages of peptide vaccinations over conventional organism vaccines. Most crucially, it permits the immune response to concentrate solely on relevant epitopes while avoiding those that trigger non-protective responses, immune evasion, or undesirable side effects such autoimmunity [23]. Furthermore, peptide vaccines are safe and cost-effective approach. Because of its ease of manufacture and simple composition, peptide vaccines are generally inexpensive to produce [24].

The main objective of this study was to develop a multi-epitope vaccine against Bovine coronavirus (BCoV) glycoprotein Spike (S) and hemagglutinin-esterase (HE) glycoproteins, using several immunoinformatics tools.

Materials and methods

Retrieval of protein sequences

The sequences of spike (S) and hemagglutinin esterase (HE) proteins of bovine coronavirus (BCoV) were obtained from the National Center for Biotechnology Information’s (NCBI) GeneBank database in October 2021 (http://www.ncbi.nlm.nih.gov/protein/).

Multiple sequence alignment and construction of phylogenetic tree

Multiple sequence alignment (MSA) was achieved on the recovered spike (S) and hemagglutinin proteins (HE) bovine coronavirus (BCoV) using ClustalW program in BioEdit software version 7.2.5 [25]. The maximum likelihood parameter was used to create phylogenetic trees of S and HE protein sequences from bovine coronavirus (BCoV) using MEGA7.0.26 (7170509) software [26].

Structural analysis, physicochemical and antigenic assessments of the target proteins

The antigenicity qualities of the reference sequences S (NP_150077.1) and HE (NP_150076.1) proteins of BCoV were assessed using the VaxiJen v2.0 server [27]. The physico-chemical properties of the target protein sequences of BCoV were evaluated using the online tool ProtParam, which included GRAVY (Grand average of hydropathicity), half-life, molecular weight, stability index, and amino acids atomic composition [28]. Secondary structures of BCoV, S and HE proteins were predicted using PSIPRED (http://bioinf.cs.ucl.ac.uk/psipred/) and the GOR IV servers (https://npsa-prabi.ibcp.fr/cgi-bin/secpred sopma.pl) [29] The existence of disulphide bonds was expected using DIANNA v1.1. (http://clavius.bc.edu/ ~clotelab/DiANNA/). TMHMM, an online tool (http://www.cbs.dtu.dk/services/TMHMM/). Existed used to predict the trans-membrane topologies of S and HE proteins [30]. DD-BLAST (http://www.ncbi.nlm.nih.gov/BLAST/) and PFAM (https://pfam.xfam.org/) [31,32,33,34] were used to search the defined conserved domains in the targeted protein sequences. Blastp in NCBI database (https://blast.ncbi.nlm.nih.gov/Blast.cgi) using reference sequence (refseq - protein) database was used to compare spike reference sequences of different coronaviruses in human and animals against IBV spike protein sequence. COBALT multiple alignments were also used to construct the phylogenetic tree (https://www.ncbi.nlm.nih.gov/blast/treeview/treeView.cgi) [31] 2.4 Prediction of B-cell epitopes The free online server, Immune Epitope Database (IEDB) (http://tools.iedb.Org/mhci/) [35] was used to predict B epitopes in S and HE proteins. Linear B-cell epitopes were forecasted using BepiPred from IEDB [36]. Emini surface accessibility prediction tool was used to predict surface epitopes [37]. However, the antigenic epitopes were inspected using kolaskar and Tongaonkar antigenicity method [38]. 2.5 Prediction of MHC I The MHC Class-I prediction tool was utilized using Immune Epitope Database (IEDB) (http://tools.iedb.org/mhci/). Cow alleles (BoLAHD6, BoLA-JSP.1, BoLA-T2a, BoLA-T2b, and BoLA-T2c) were used to predict CTL epitopes of MHC Class-I target proteins [39]. The reference sequence of the S and HE proteins of bovine coronaviruses (BCoV) were submitted in FASTA format; the Percentile ranks of binding affinity to specific cow MHC-I molecules were taken within the range of 0–2.

Assessment of B and T-cell epitopes

Conservation analysis

The conservation of anticipated epitopes was determined using Epitope Conservancy Analysis in the Immune Epitope Database (IEDB) (http://tools.iedb.org/conservancy/) [40]. The conservation score was calculated using a sequence identity criterion of 100%. 2.6.2Antigenicity, allergenicity and Toxicity The antigenicity qualities of the conserved predicted epitopes of S and HE proteins of BCoV were assessed using the VaxiJen v2.0 server (http://www.jenner.ac.uk/VaxiJen) [27]. Based on a training set containing known allergens and non-allergens from various species, AllerTop v.2.0 categorizes the allergens of a protein utilizing amino acid E-descriptors, k nearest neighbor’s machine learning algorithms, and the auto- and cross-covariance (ACC) transformation (http://www.ddg-pharmfac.net/AllerTOP) [41]. The ToxinPred server (http://crdd.osdd.net/raghava/toxinpred/) was also utilized to assess epitopes toxicity as well as mutation and physico-chemical characteristics. ToxinPred predicts the peptides’ toxicity based on the patterns found in hazardous peptides. As a result, motifs were utilized in the prediction algorithms [42]. 2.6.3 Structural modelling Three-dimensional (3D) structures of the reference Spike S protein were modeled using phyer 2 server [43]. Whereas, Raptor X online server [44, 45] was used to model the reference HE protein of BCoV. Visualization of the result was performed using Chimera software 1.8 [46]. 3. Results: 3.1 Multiple sequence alignment Multiple sequence alignment of all retrieved sequences of S and HE proteins using ClustalW through BioEdit software showed high conservancy between the aligned sequences.

Antigenicity, physicochemical properties and structural analysis of the target proteins of BCoV

VaxiJen predicted scores of 0.4786 and 0.5367 for S and HE proteins respectively, indicating that the sequences were effectively antigenic. Using Expasy tools in protparam, physical and chemical properties of S and HE protein sequences from bovine coronavirus (BCoV) revealed that both targeted proteins were stables (Table 1).

The spike S protein, showed that it contained.

1363 amino acids (aa) with a molecular weight of 150235.48 kDa, while, HE protein contained 424 amino acids (aa) with a molecular weight of 47617.94 kDa (Table 1).

Table 1 Physical and chemical properties of S and HE protein sequences of BCoV using protparam server

The secondary structure of BCoV spike S and HE proteins were analyzed through PSIPRED and GOR IV server Fig. 1.

Fig. 1
figure 1

(a): Secondary structure of Spike (S); (b) Secondary structure of Heamagglutinase Estrase (HE) proteins using GOR IV server

DiANNA1.1 tool calculated 28 and 7 disulphides bond (S–S) positions in S and HE proteins respectively. The trans-membrane protein topology was investigated via online Tool TMHMM. Residues from 1 to 1307 were found to be exposed to the surface, residue from 1308 to 1330 were found inside trans-membrane-region and residues from 1331 to 1363 were buried within the core-region of the S protein. In HE protein, residues from 1 to 392 were found to be exposed to the surface, residue from 383 to 415 were found inside trans-membrane-region and residues from 416 to 424 were buried within the core-region. The closest homologs for S and HE proteins, when comparing different coronaviruses in humans and animals with S and HE protein sequences of BCoV, are Rabbit coronavirus and Rabbit coronavirus HKU14 protein, respectively, followed by spike surface glycoprotein of Human coronavirus and hemagglutinin-esterase of Human coronavirus OC43. Based on multiple alignments of COBALT, a phylogenetic tree of BCoV viruses and other coronaviruses on humans and animals was constructed (Tables 2 and 3).

Table 2 Blastp similarity search of Spike S protein of BCoV against other spike proteins of coronaviruses in human and animals
Table 3 Blastp similarity search of HE protein of BCoV against other HE proteins of other coronaviruses in human and animals

Prediction of B-cell epitopes

Utilizing the Bepipred Linear Epitope Prediction program, several epitopes were predicted using B cell prediction methods. The conservation of predicted epitopes was determined using Epitope Conservancy Analysis in the Immune Epitope Database (IEDB) (http://tools.iedb.org/conservancy/). The conservation scores were calculated using a sequence identity criterion of 100%.

Twenty-five linear conserved epitopes from spike S and fifteen epitopes from HE proteins were identified when the predicted epitopes were shortened. Between locations 379 and 1208 of spike S protein, three epitopes (379TCQPQ384, 720QLQPINY726 and 1204YYYPE1208) with varying lengths were identified as linear, surface, and antigenic epitopes (Table 4; Fig. 2).

In HE protein, however, no epitopes have been discovered as linear, surface, or antigenic.

Table 4 List of B cell epitopes of BCoV proteins predicted by different B cell methods
Fig. 2
figure 2

The position of proposed B cell epitopes in the 3D structure of reference spike (S) protein of BCoV using Chimera software 1.8

Prediction of MHC1 epitopes

Only seven epitopes showed a high linkage with MHCI bovine alleles for spike glycoprotein (77NMALKGTLL85, 11TAFAVIGDL19, 986AAAGVPFYL994, 865TQLQVANSL873, 1103SQQLSDSTL1111, 1300YEYYVKWPW1308, and 972TLAATSASL980) and three epitopes in HE protein showed a high linkage with MHCI bovine allele (56SYMDLNPAL65, 186YKVEADFYL194 and 360SSVWPLYPY368).

Assessment of B-cell and MHC I epitopes

Antigenicity, allergenicity and toxicity of MHCI epitopes

VaxiJen v2.0 server, AllerJen v2.0, and ToxiPred were used to estimate the possible antigenicity, allergenicity, and toxicity of MHCI predicted epitopes. The S protein, B cell epitope 720QLQPINY726, and MHC1 epitope 77NMALKGTLL85 showed significant antigenicity, no allergencity, and no toxicity. While using the VaxiJen v2.0 server, the HE protein MHC1 epitope 56SYMDLNPAL65 showed a high antigenicity score as well as no allergencity or toxicity (Tables 5 and 6; Figs. 3 and 4.)

Table 5 List of MHC1 epitopes of S protein of BCoV
Table 6 List of MHC1 epitopes of HE protein of BCoV
Fig. 3
figure 3

The position of proposed MHC1 epitopes in the 3D structure of reference spike (S) protein of BCoV using Chimera software 1.8

Fig. 4
figure 4

The position of proposed MHC1 epitopes in the 3D structure of reference HE protein of BCoV using Chimera software 1.8

Discussion

To improve protection against viral infections, a single antigen or a group of antigens must be presented to the immune system [23]. In silico vaccine design can be thought of as picking good virus protein fragments and then putting them together to make a final vaccine. As a subunit of the final vaccine, a fragment with numerous merits can be picked. An ideal subunit, for example, would have several B-cell and T-cell epitopes, as well as a high antigenicity to elicit human and animals protective responses [47]. An antigenic epitope is a basic unit that triggers a cellular or humoral immune response. As a result, a multi-epitope vaccination composed of a series of or overlapping peptides is a highly effective strategy to prevent and treat malignancies and viral infections [23, 48]. Epitopes capable of inducing both B-cell and T-cell immunity are known to be effective vaccine candidates [49]. However, in this study no overlapping between B-cell and T-cell Epitopes was detected. Peptide-based vaccinations are developed using an emerging computational method that combines immunoinformatic prediction with rigorous experimental validation, making epitope identification within protein antigens much easier [50]. These vaccines are designed to accurately induce immune responses against antigens using essential epitope peptides rather than the full antigen, and hence have various advantages over traditional vaccines, including improved specificity, increased safety, reduced costs, and fewer adverse events [51].

In the present study, spike (S) and hemagglutinin-esterase (HE) glycoproteins of BCoV were analyzed using a variety of online and offline immunoinformatics tools. Protparam protein characterization of the BCoV spike S and HE glycoproteins indicated that both proteins are positive in nature and stables.

A protein whose instability index is smaller than 40 is predicted as stable, a value above 40 predicts that the protein may be unstable.

The proteins also showed good antigenic characteristics using the Vaxijen 2.0v server.

ClustalW was used to align sequences. The aligned sequences of the two analyzed proteins (S and HE) showed strong conservancy between sequences. On the other hand, a few regions had different amino acids in various sequences. We analyzed the evolutionary divergence of each protein.

B-cell epitope prediction is critical for the development of vaccine components and immunodiagnostic reagents. Continuous epitopes are used in the majority of epitope prediction algorithms Linear B cell epitopes have been shown to play an important function in viral neutralization [52, 53]. Based on amino acid qualities such as hydrophilicity, surface accessibility, flexibility, and antigenicity, the IEDB prediction algorithm was utilized to predict linear, surface, and antigenic epitopes [23]. It is known that both cellular and humoral immune responses are essential against coronaviruses infection [54]. The development of neutralizing antibodies, in addition to T cell-mediated immunity, is required to protect against virus infection, such as Coronaviruses. While cytotoxic lymphocytes can kill infected cells, antibodies have the ability to destroy infected cells while also preventing the infectious virus from infecting the cell (neutralization) [54, 55]. In this study, three shortened conserved epitopes 379TCQPQ384, 720QLQPINY726 and 1204YYYPE1208 were predicted as surface, linear, and antigenic epitopes using B cell prediction algorithms. The top epitope was 720QLQPINY726 due to high safety and antigenicity.

BoLA is called bovine leukocyte antigen (the MHC in cattle) and is located on chromosome 23; it shares the overall structure of other mammals’ MHC. MHC is an adaptive immune system component that addresses both immunological and evolutionary biological issues at the same time [56].

In the bovine genome project, the bovine MHC-II locus was not sequenced completely [57, 58]. As a result, the analysis was completed with the BoLA MHC-I alleles Only. Cell-mediated immunity induced by cytotoxic T lymphocytes (CTLs) is vital for the defense against viral diseases. CTLs are responsible for the immune elimination of intracellular, pathogens such as viruses because these cells recognize the presented endogenous antigenic peptides by the MHC class I molecules.

The induction of CTL is highly dependent on how effectively the antigenic peptide is delivered into the major histocompatibility complex (MHC) class I presentation pathway in antigen presenting cells (APCs) [24]. MHCI prediction methods on the other hand, revealed seven, high linkage epitopes with MHCI allele from S protein (77NMALKGTLL85, 11TAFAVIGDL19, 986AAAGVPFYL994, 865TQLQVANSL873, 1103SQQLSDSTL1111, 1300YEYYVKWPW1308, and 972TLAATSASL980). In HE protein, only three epitopes showed a high linkage with MHCI allele (56SYMDLNPAL65, 186YKVEADFYL194 and 360SSVWPLYPY368). The predicted epitopes exhibited high levels of antigenicity, no allergenicity, and no toxicity for both spike S and HE proteins.

Conclusion

Three B cell epitopes were predicted from spike glycoprotein of BCoV in this investigation, but the best epitope was 720QLQPINY726 of BCoV because of high antigenicity. However, the HE did not predict B cell epitopes. Several T cell epitopes were predicted from S and HE glycoproteins of BCoV. The predicted epitopes exhibited high safety and antigenicity as well as strong binding with bola alleles. The proposed epitopes should be included into a multi-epitopes vaccine against BCoV and evaluated using in vivo and in vitro trails. By generating humoral and cellular responses, this vaccine could be used as a peptide vaccine to reduce BCoV infection in bovines. The peptide vaccination against BCoV’s spike protein (S) and (HE) which designed to cover all strains in different serotypes, might be effectively replace existing vaccines, reducing recurrent outbreaks and their associated significant economic losses.

Data availability

All the data supporting the findings are contained within the manuscript.

References

  1. Peel DS. The effect of market forces on bovine respiratory disease. Veterinary Clinics: Food Anim Pract. 2020;36(2):497–508.

    Google Scholar 

  2. Mebus CA, Stair EL, Rhodes MB, Twiehaus MJ. Neonatal calf diarrhea: propagation, attenuation, and characteristics of a coronavirus-like agent. Amer J Vet Res. 1973;34:145–50.

    PubMed  CAS  Google Scholar 

  3. Saif LJ. A review of evidence implicating bovine coronavirus in the etiology of winter dysentery in cows: an enigma resolved? Cornell Vet. 1990;80:303–11.

    PubMed  CAS  Google Scholar 

  4. Saif LJ. Coronaviruses of domestic livestock and poultry: interspecies transmission, pathogenesis, and immunity. In: Perlman S, Gallagher T, Snijder EJ, editors. The nidoviruses. Washington, DC: ASM; 2007. pp. 279–96.

    Google Scholar 

  5. Saif LJ. Bovine respiratory coronavirus. Vet Clin North Am Food Anim Pract. 2010;26:349–64.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Fulton RW. Bovine respiratory disease research (1983–2009). Anim Health Res Rev. 2009;10:130–9.

    Article  Google Scholar 

  7. Fulton RW. Viral diseases of bovine respiratory tract: bovine herpesvirus-1, parainfluenza – 3 virus, bovine respiratory syncytial virus, bovine adenoviruses, bovine coronavirus, and bovine viral diarrhea viruses. In: Anderson DE, Rings DM, editors. Current veterinary therapy-food animal practice. St. Louis, MO: Saunders Elsevier; 2008. pp. 171–91.

    Google Scholar 

  8. Fulton RW, Herd HR, Sorensen NJ, Confer AW, Ritchey JW, Ridpath JF, Burge LJ. Enteric disease in postweaned beef calves associated with bovine coronavirus clade 2. J Veter Diagn Investig. 2015;27:97–101.

    Article  Google Scholar 

  9. Cho KO, Hasoksuz M, Nielsen PR, et al. Cross-protection studies between respiratory and calf diarrhea and winter dysentery coronavirus strains in calves and RT-PCR and nested PCR for their detection. Arch Virol. 2001;146:2401–19.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  10. Brownlie J. Coronaviridae. In: MacLachlan NJ, Dubovi EJ, editors. Fenner’s Veterinary Virology. 5th ed. San Diego, California: Elsevier; 2017. pp. 435–59.

    Google Scholar 

  11. Saif LJ, Alhamo M. Bovine coronavirus infection. In: Coetzer JAW, Thomson GR, Maclachlan NJ, Penrith ML, editors. Infectious diseases of livestock. 3rd ed. Oxford: Oxford University Press; 2018. https://anipedia.org/resources/bovine-coronavirus-infection/1033., United Kingdom.

    Google Scholar 

  12. Andersen KG, Rambaut A, Lipkin WI, Holmes EC, Garry RF. The proximal origin of SARS-CoV-2. Nat Med. 2020;26:450–2.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  13. Lai MHKV. Coronaviridae: the viruses and their replication. In: Fields BN, Howley PM, DMK, editors. Fields Virology. Philadelphia: Lippincott - Raven; 2001. pp. 1163–86.

    Google Scholar 

  14. Deregt D, Gifford GA, Ijaz MK, et al. Monoclonal antibodies to bovine coronavirus glycoproteins E2 and E3: demonstration of in vivo virus-neutralizing activity. J Gen Virol. 1989;70(Pt 4):993–8.

    Article  PubMed  CAS  Google Scholar 

  15. Compendium of veterinary products. 12th ed. Huron P, MI: North American Compendiums; 2010. p. 1–1848.

  16. Hasoksuz M, Sreevatsan S, Cho KO, Hoet AE, Saif LJ. Molecular analysis of the S1 subunit of the spike glycoprotein of respiratory and enteric bovine coronavirus isolates. Virus Res. 2002;84:101–9.

    Article  PubMed  CAS  Google Scholar 

  17. Adams MJ, Carstens EB. Ratification vote on taxonomic proposals to the International Committee on Taxonomy of Viruses (2012). Arch Virol. 2012;157:1411–22.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  18. Woo PC, Lau SK, Lam CS, Lau CC, Tsang AK, Lau JH, Bai R, Teng JL, Tsang CC, Wang M, et al. Discovery of seven Novel mammalian and avian coronaviruses in the Genus Deltacoronavirus supports Bat coronaviruses as the Gene source of Alphacoronavirus and Betacoronavirus and Avian coronaviruses as the Gene source of Gammacoronavirus and Deltacoronavirus. J Virol. 2012;86:3995–4008.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  19. ICTV Virus Taxonomy. 2019 Release. Available online: https://talk.ictvonline.org/ictv-reports/ictv_9th_report/positive-sense-rna-viruses-2011/w/posrna_viruses/223/coronaviridae-figures (accessed on 23 March 2021).

  20. McDaniel CJ, Cardwell DM, Moeller RB, Gray GC. (2014). Humans and cattle: a review of bovine zoonoses. Vector Borne Zoonotic Dis. (2014);14(1):1-19.https://doi.org/10.1089/vbz.2012.1164

  21. Usmani SS et al. In silico tools and databases for designing peptide-based vaccines and drugs. 2018;112:221-63.

  22. Ali SA, Almofti YA. and K.A.J.A.i.b. Abd-Elrahman, Immunoinformatics approach for multiepitopes vaccine prediction against glycoprotein B of avian infectious laryngotracheitis virus. 2019. 2019.

  23. Li W et al. Peptide vaccine: progress and challenges. 2014;2(3):515-36.

  24. Hall T, Biosciences I, Carlsbad CJGBB. BioEdit: Important Softw Mol Biology. 2011;2(1):60–1.

    Google Scholar 

  25. Kumar S et al. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. 2016;33(7):1870-74.

  26. Doytchinova IA, D.R.J. .B.b. Flower, VaxiJen: a server for prediction of protective antigens, tumour antigens and subunit vaccines. 2007;8(1):1-7.

  27. Gasteiger E et al. ExPASy: the proteomics server for in-depth protein knowledge and analysis. 2003;31(13):3784-8.

  28. Buchan DW et al. Scalable web services for the PSIPRED Protein Analysis Workbench. 2013;41(W1):W349-W357.

  29. Ferrè F, Clote PJB. Disulfide Connectivity Prediction Using Secondary Struct Inform Diresidue Frequencies. 2005;21(10):2336–46.

    Google Scholar 

  30. Altschul SF et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. 1997;25(17):3389-402.

  31. Marchler-Bauer A et al. CDD: a conserved domain database for interactive domain family analysis. 2007;35(suppl_1):D237-D240.

  32. Eddy SRJB. Profile Hidden Markov Models. 1998;14(9):755–63.

    CAS  Google Scholar 

  33. El-Gebali S et al. The Pfam protein families database in 2019. 2019;47(D1):D427-D432.

  34. Vita R et al. The immune epitope database (IEDB) 3.0. 2015;43(D1):D405-D412.

  35. Larsen JEP, Lund O, Nielsen MJIr. Improved Method Predicting Linear B-cell Epitopes. 2006;2(1):1–7.

    Google Scholar 

  36. Emini EA et al. Induction of hepatitis a virus-neutralizing antibody by a virus-specific synthetic peptide. 1985;55(3):836–9.

  37. Kolaskar AS. And P.C.J.F.l. Tongaonkar, a semi-empirical method for prediction of antigenic determinants on protein antigens. 1990;276(1–2):172–4.

  38. Lundegaard C, Lund O, Nielsen MJB. Accurate approximation method for prediction of class I MHC affinities for peptides of length 8, 10 and 11 using prediction tools trained on 9mers. 2008;24(11):1397–1398.

  39. Bui H-H et al. Development of an epitope conservancy analysis tool to facilitate the design of epitope-based diagnostics and vaccines. 2007;8(1):1–6.

  40. Dimitrov I, et al. AllerTOP v 2—a Serv silico Prediction Allergens. 2014;20(6):1–6.

    CAS  Google Scholar 

  41. Gupta S, et al. silico Approach Predicting Toxic Peptides Proteins. 2013;8(9):e73957.

    CAS  Google Scholar 

  42. Kelley LA et al. The Phyre2 web portal for protein modeling, prediction and analysis. 2015;10(6):845–858.

  43. Källberg M et al. Template-based protein structure modeling using the RaptorX web server. 2012;7(8):1511–22.

  44. Peng J, Xu JJPS, Function, Bioinformatics. RaptorX: Exploiting Struct Inform Protein Alignment Stat Inference. 2011;79(S10):161–71.

    CAS  Google Scholar 

  45. Peng J, Xu JJPS, Function, Bioinformatics. A multiple-template approach to protein threading. 2011. 79(6): pp. 1930–9.

  46. Chan WM et al. User’s manual for Chimera grid tools, version 1.8. 2003.

  47. Yang Z, Bogdan P, Nazarian SJSr. An in silico deep learning approach to multi-epitope vaccine design: a SARS-CoV-2 case study. 2021;11(1):1–21.

  48. Khalid H, U.A.J.J.o.B I, Ashfaq. Exploring HCV genome to construct multi-epitope based subunit vaccine to battle HCV infection: Immunoinformatics based approach. 2020;108:103498.

  49. Verma SK, Yadav S. and A.J.A.b.r. Kumar, In silico prediction of B-and T-cell epitope on Lassa virus proteins for peptide based subunit vaccine design. 2015;4.

  50. Fridkis-Hareli M, Hoshino Y. Peptide-Based Immunotherapeutics and Vaccines. J Immunol Res. 2018;2018:9.

  51. Wang B et al. Peptide-Based Vaccination Therapy for Rheumatic Diseases. 2020;2020.

  52. Awadelkareem EA, S.A.J.T.M C, Ali. Vaccine design of coronavirus spike (S) glycoprotein in chicken: immunoinformatics and computational approaches. 2020;5(1):1–16.

  53. Yao B et al. Conformational B-cell epitope prediction on antigen protein structures: a review of current algorithms and comparison with common binding site prediction methods. 2013;8(4):e62249.

  54. Enayatkhani M et al. Reverse vaccinology approach to design a novel multi-epitope vaccine candidate against COVID-19: an in silico study. 2021;39(8):2857–72.

  55. Klein J, Sato A, Nikolaidis N. MHC, TSP, and the origin of species: from immunogenetics to evolutionary genetics. Annu Rev Genet. 2007;41:281–304. https://doi.org/10.1146/annurev.genet.41.110306.130137.

    Article  PubMed  CAS  Google Scholar 

  56. Nene V, Svitek N, Toye P et al. Designing bovine T cell vaccines via reverse immunology, Ticks and Tick-borne Diseases, 3;3:188–92, 2012.

  57. Sequencing BG, Consortium A, Elsik CG, Tellam RL, Worley andKC. The genome sequence of taurine cattle: a window to ruminant biology and evolution. Science. 2009;324:522–8.

    Article  Google Scholar 

  58. Zimin AV, Delcher AL, Florea L et al. A whole-genome assembly of the domestic cow, Bos taurus. Genome Biol, 10, no. 4, article R42, 2009.

Download references

Acknowledgements

Author want to acknowledgement Sudanese Veterinarian Association for capacity building and Training, Khartoum, Sudan for their support and engorgement. In addition, all classmate and staff of master program of Bioinformatics, Department of Molecular Biology and Bioinformatics, Collage of Veterinary Medicine, University of Bahri, Khartoum, Sudan.

Funding

Author want to acknowledgement Sudanese Veterinarian Association for capacity building and Training, Khartoum, Sudan for their support and engorgement.

Author information

Authors and Affiliations

Authors

Contributions

Eman, A. Awadelkareem and Hamdoun S designed this study, accomplished the data collection and analyze and results’ writing. Eman, A. Awadelkareem interpreted the data and wrote the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Eman A. Awadelkareem.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have 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/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Awadelkareem, E.A., Hamdoun, S. Insilco vaccine design of spike and hemagglutinin esterase proteins of bovine coronavirus. transl med commun 9, 28 (2024). https://doi.org/10.1186/s41231-024-00190-y

Download citation

  • Published:

  • DOI: https://doi.org/10.1186/s41231-024-00190-y

Keywords