Engineering and Development of a Novel Chimera of HBHA-OMP25 as a Potent Recombinant Vaccine Against Brucellosis: A Bioinformatics Investigation

authors:

avatar Maryam Faroukhpoor 1 , avatar Amin Jaydari ORCID 1 , * , avatar Ali Forouharmehr 2 , avatar Nemat Shams 1 , avatar Narges Nazifi 3

Department of Pathobiology, Faculty of Veterinary Medicine, Lorestan University, Khorramabad, Iran
Department of Animal Science, Faculty of Agriculture, Lorestan University, Khorramabad, Iran
Department of Basic Science, Faculty of Veterinary Medicine, Lorestan University, Khorramabad, Iran
Journal of Inflammatory Diseases: Vol.27, issue 2; e150899
published online: August 4, 2024
article type: Research Article
received: July 7, 2024
accepted: August 4, 2024
DOI: https://doi.org/10.5812/jid-150899

how to cite: Faroukhpoor M, Jaydari A, Forouharmehr A, Shams N, Nazifi N. Engineering and Development of a Novel Chimera of HBHA-OMP25 as a Potent Recombinant Vaccine Against Brucellosis: A Bioinformatics Investigation. J Inflamm Dis. 2023;27(2):e150899. https://doi.org/10.5812/jid-150899.

Abstract

Background:

Brucellosis is a zoonotic disease widely reported in most developing countries. Despite vaccination programs for domestic animals, this disease continues to cause economic losses and human infections.

Objectives:

The current project aimed to design a recombinant construct using a bioinformatics approach against brucellosis.

Methods:

To design the recombinant construct, amino acid sequences of heparin-binding hemagglutinin adhesin (HBHA) and outer membrane protein 25 (OMP25) were collected from databases, and then the HBHA-OMP25 chimera was assembled using an appropriate linker. The most important physical and chemical features, antigenicity, secondary structure, three-dimensional structure, and structural refinement of the designed chimera were evaluated using the best online tools. Moreover, the interaction between the HBHA moiety of the designed chimera (as a ligand) and toll-like receptor 4 (TLR4)/MD2 (as a receptor) was assessed using protein-protein docking. Finally, the codons of the coding DNA sequence of the designed construct were adapted using the JCat tool and cloned into pET21a (+).

Results:

The findings confirmed that the designed construct, with a protein weight of 45.172 kDa and an antigenicity index of 0.7184, was an antigenic protein capable of triggering immune responses. Additionally, our results demonstrated that the HBHA moiety of the designed chimera could interact with the TLR4/MD2 receptor using nine hydrogen bonds. Finally, our results showed that the designed chimera, with a codon adaptation index (CAI) of 0.97, was successfully cloned into pET21a (+) between the BamHI and NotI restriction enzymes.

Conclusions:

It seems that our designed chimera could be an appropriate candidate for combating brucellosis.

1. Background

Brucellosis is a zoonotic disease that infects about 500,000 people annually (1). This disease has been observed in many geographical regions, including the Mediterranean, Middle East, Central and Southern America, Africa, and Asia (2). Brucellosis can be transmitted through the secretions of infected animals, but the consumption of unpasteurized dairy products, such as goat milk and fresh cheese, is the most common cause of infection. For humans, antibiotic therapy is recommended to treat brucellosis. In contrast, infected animals are often slaughtered due to a lack of drug adaptation and affordability (3).

To prevent brucellosis in animals, various strategies are employed, including inoculation, adherence to sanitary principles, and detection and elimination of infected animals. In general, attenuated vaccines are used to control brucellosis in humans. Although these vaccines have had some success, their application has been limited due to low efficiency and the issue of antibiotic resistance. Therefore, given the sanitary and economic importance of brucellosis, along with serious problems in its treatment and control, novel and precise methods should be investigated.

An efficient vaccine used to control brucellosis should have several advantages, including stability in the host body, the absence of disease triggering in the host, prevention of abortion, and the capacity for large-scale production (4). Currently, recombinant vaccines (RVs) employing different antigenic proteins are recommended as an intriguing method to control various infectious diseases (5). These vaccines not only have all the advantages of an efficient vaccine but can also be flexibly designed for specific applications. Despite the advantages of RVs, low antigenicity is known as their most significant disadvantage; consequently, a molecular adjuvant should be conjugated with RVs to enhance their immune responses (6).

In general, outer membrane proteins (OMPs) are known as immunogenic proteins in Brucella that can strongly stimulate immune responses (1). Different members of OMPs have been evaluated for Brucella; among these proteins, outer membrane protein 25 (OMP25) has been strongly recommended for RV design due to its protective properties (7). Mycobacterium tuberculosis heparin-binding hemagglutinin adhesin (HBHA) has been extensively applied as a molecular adjuvant for designing RVs (6, 8). This protein, as a novel antagonist of the toll-like receptor 4 (TLR4) receptor without toxicity in the host body, can play a role in the maturation of dendritic cells and enhancement of immune responses (9).

2. Objectives

The present study was conducted to design and engineer an RV based on the OMP25 protein of Brucella (as an antigenic moiety) and the HBHA protein of M. tuberculosis (as a molecular adjuvant moiety).

3. Methods

3.1. Protein Achievement

To carry out the current project, the amino acid sequences of OMP25 and HBHA proteins, with accession numbers JF918759.1 and NC-000962.39 (respectively), were obtained from the National Center for Biotechnology Information (https://www.ncbi.nlm.nih.gov/). Additionally, the PDB file of the TLR4/MD2 receptor, with an accession number of 2z64, was obtained from the Protein Data Bank (https://www.rcsb.org/).

3.2 Engineering a Chimera of HBHA-OMP25

To engineer a chimera of HBHA-OMP25, CLC Main WorkBench software version 5.2 was employed. In this process, the first signal peptide sequences of both proteins (HBHA and OMP25) were removed, and the correction of the reading frames was evaluated. Finally, the amino acid sequence of the HBHA protein (at the N terminal) was linked to the amino acid sequence of the OMP25 protein (at the C terminal) using an appropriate rigid linker.

3.3 Investigation of Important Features

To investigate the physical and chemical features (e.g., protein weight, GRAVY, pI, stability indicator, and aliphatic indicator) of the HBHA-OMP25 chimera, the ProtParam tool was used. Additionally, the VaxiJen tool with a cutoff of 0.5 was employed to evaluate the antigenicity of the designed construct, while its allergenicity was investigated using the AllerTOP tool.

3.4. Investigation of Different Structures

To investigate the secondary structure of the chimera, the SOPMA tool was applied. In this case, the amino acid sequence of the chimera was uploaded to the tool. Additionally, to predict the three-dimensional (3D) structure of the chimera, the I-TASSER tool was used. Refinement of the 3D structure was conducted using the GalaxyRefine tool, and the best-refined model was chosen based on the outputs from the VADAR tool.

3.5. Protein-Protein Docking

To investigate the docking between the HBHA moiety of the chimera and the TLR4/MD2 receptor, a protein-protein docking analysis was performed using the ClusPro tool. For this purpose, the PDB format of the best-refined model of the chimera and the PDB file of the TLR4/MD2 receptor were uploaded to the tool. Additionally, the hydrogen bonds of the complex were analyzed using LigPlot+ software.

3.6. Codon Adaptation and Computer-Aided Cloning

In this project, codon adaptation of the chimera nucleotide sequence for expression in a prokaryotic system was performed using the JCat tool (http://www.jcat.de/). The JCat tool employs a precise algorithm to predict the codon adaptation index (CAI), as presented by Carbone et al. (10). This algorithm can prepare a reference set of genes that are highly representative of the bias. Finally, the cloning potential of the adapted sequence into pET21a (+) was verified. In this process, BamHI and NotI restriction enzyme sites were added to the 5' and 3' regions of the adapted sequence.

4. Results

4.1. Engineering a Chimera of HBHA-OMP25

The engineered chimera, from the N-terminal to the C-terminal, contained HBHA and OMP25 moieties, respectively (Figure 1). Additionally, an EAAAK linker, consisting of five amino acid residues, was embedded between the two moieties (Figure 1). It must be noted that the entire designed chimera contained 417 amino acid residues.

Final construct of the engineered chimera which was designed by CLC WorkBench software. Molecular adjuvant [heparin-binding hemagglutinin adhesin (HBHA)] moiety, rigid linker (EAAAK) and antigenic moiety [outer membrane protein 25 (OMP25)] have been illustrated by green, red and blue, respectively.
Final construct of the engineered chimera which was designed by CLC WorkBench software. Molecular adjuvant [heparin-binding hemagglutinin adhesin (HBHA)] moiety, rigid linker (EAAAK) and antigenic moiety [outer membrane protein 25 (OMP25)] have been illustrated by green, red and blue, respectively.
Figure 1.

Final construct of the engineered chimera which was designed by CLC WorkBench software. Molecular adjuvant [heparin-binding hemagglutinin adhesin (HBHA)] moiety, rigid linker (EAAAK) and antigenic moiety [outer membrane protein 25 (OMP25)] have been illustrated by green, red and blue, respectively.

4.2. Investigation of Important Features

Based on the physical and chemical results, the protein weight, pI, aliphatic index, instability index, and GRAVY of the chimera were 45.172 kDa, 9.03, 80.46, 34.95, and -0.450, respectively. Additionally, the half-life of the engineered chimera in mammals, yeast, and E. coli was 30 hours, > 20 hours, and > 10 hours, respectively. According to the results, the antigenicity score of the engineered chimera was 0.7184, classifying it as an antigen, while the allergenicity results confirmed that the engineered chimera was a non-allergenic protein.

4.3. Investigation of Different Structures

As mentioned, to predict the secondary structure of the engineered chimera, the SOPMA tool was used. The results demonstrated that the engineered chimera contained 52.30% alpha helix, 17.97% extended strand, 4.61% beta turn, and 25.12% random coil (Figure 2). The 3D structure of the engineered chimera was predicted by the I-TASSER tool with a C-score of -1.32 (Figure 3). Furthermore, the Ramachandran analysis confirmed that 89% of the amino acid residues of the most refined model were in the core region, while only 1% were in the disallowed region (Figure 4).

Secondary structure of the engineered chimera. The structure contains different states including extended strand, random coil, beta turn and alpha helix.
Secondary structure of the engineered chimera. The structure contains different states including extended strand, random coil, beta turn and alpha helix.
Figure 2.

Secondary structure of the engineered chimera. The structure contains different states including extended strand, random coil, beta turn and alpha helix.

3D structure of the engineered chimera predicted by I-TASSER tool. Heparin-binding hemagglutinin adhesin (HBHA) moiety, EAAAK linker and outer membrane protein 25 (OMP25) moiety have been shown by green, red and blue, respectively.
3D structure of the engineered chimera predicted by I-TASSER tool. Heparin-binding hemagglutinin adhesin (HBHA) moiety, EAAAK linker and outer membrane protein 25 (OMP25) moiety have been shown by green, red and blue, respectively.
Figure 3.

3D structure of the engineered chimera predicted by I-TASSER tool. Heparin-binding hemagglutinin adhesin (HBHA) moiety, EAAAK linker and outer membrane protein 25 (OMP25) moiety have been shown by green, red and blue, respectively.

Ramachandran plot of the most excellent- refined model. As shown in this plot, 89% of the amino acid were in the core region (red area).
Ramachandran plot of the most excellent- refined model. As shown in this plot, 89% of the amino acid were in the core region (red area).
Figure 4.

Ramachandran plot of the most excellent- refined model. As shown in this plot, 89% of the amino acid were in the core region (red area).

4.4. Protein-Protein Docking

As mentioned, in the current project, a protein-protein docking strategy was employed to investigate the interaction between the HBHA moiety of the engineered chimera and the TLR4/MD2 receptor. The results confirmed that the HBHA moiety could successfully dock with the TLR4/MD2 receptor with a lowest energy of -0.777 kcal/mol (Figure 5A). Additionally, the results from LigPlot+ software revealed that nine hydrogen bonds were formed in this interaction (Figure 5B).

Interaction between the chimera of heparin-binding hemagglutinin adhesin (HBHA)- outer membrane protein 25 (OMP25) and toll-like receptor 4 (TLR4)/MD2 was visualized by PyMol software. A, protein-protein docking between the HBHA moiety (yellow) and the TLR4/MD2 receptor (green); B, amino acids involved in the interaction (green) between the HBHA moiety and TLR4/MD2 receptor; C, hydrogen bonds (green lines) formed between the HBHA moiety and TLR4/MD2 receptor.
Interaction between the chimera of heparin-binding hemagglutinin adhesin (HBHA)- outer membrane protein 25 (OMP25) and toll-like receptor 4 (TLR4)/MD2 was visualized by PyMol software. A, protein-protein docking between the HBHA moiety (yellow) and the TLR4/MD2 receptor (green); B, amino acids involved in the interaction (green) between the HBHA moiety and TLR4/MD2 receptor; C, hydrogen bonds (green lines) formed between the HBHA moiety and TLR4/MD2 receptor.
Figure 5.

Interaction between the chimera of heparin-binding hemagglutinin adhesin (HBHA)- outer membrane protein 25 (OMP25) and toll-like receptor 4 (TLR4)/MD2 was visualized by PyMol software. A, protein-protein docking between the HBHA moiety (yellow) and the TLR4/MD2 receptor (green); B, amino acids involved in the interaction (green) between the HBHA moiety and TLR4/MD2 receptor; C, hydrogen bonds (green lines) formed between the HBHA moiety and TLR4/MD2 receptor.

4.5. Codon Adaptation and Computer-Aided Cloning

To perform codon adaptation of the nucleotide sequence of the engineered chimera, the JCat tool was used. The results of this analysis demonstrated that the CAI of the sequence before and after adaptation were 0.375 and 0.970, respectively, whereas the GC contents of the sequence before and after adaptation were 61.9% and 51.8%, respectively (Figure 6). Finally, the results of computer-aided cloning revealed that the nucleotide sequence of the engineered chimera could be correctly cloned between the BamHI and NotI enzymes of the multiple cloning sites of the pET21a (+) vector (Figure 7).

Codon adaptation graph of the engineered chimera nucleotide sequence with codon adaptation index (CAI) = 0.97
Codon adaptation graph of the engineered chimera nucleotide sequence with codon adaptation index (CAI) = 0.97
Figure 6.

Codon adaptation graph of the engineered chimera nucleotide sequence with codon adaptation index (CAI) = 0.97

Computer-aided cloning of the engineered chimera nucleotide sequence. The inserted gene (red region) with 1258 nucleotides in length was cloned between BamHI and NotI enzymes at 5' and 3' regions, respectively. The most important regions of the vector including T7 promoter and His tag have been shown with green and purple, respectively.
Computer-aided cloning of the engineered chimera nucleotide sequence. The inserted gene (red region) with 1258 nucleotides in length was cloned between BamHI and NotI enzymes at 5' and 3' regions, respectively. The most important regions of the vector including T7 promoter and His tag have been shown with green and purple, respectively.
Figure 7.

Computer-aided cloning of the engineered chimera nucleotide sequence. The inserted gene (red region) with 1258 nucleotides in length was cloned between BamHI and NotI enzymes at 5' and 3' regions, respectively. The most important regions of the vector including T7 promoter and His tag have been shown with green and purple, respectively.

5. Discussion

Brucellosis is a zoonotic disease prevalent in most developing countries (11). Despite livestock vaccination efforts, brucellosis continues to cause significant economic damages and human infections (12-14). Currently, antibiotic therapy is used to treat brucellosis; however, this method is not entirely efficient (13). In contrast, a vaccination strategy proves to be more effective and affordable (13, 14). Therefore, vaccine development is strongly recommended to eradicate brucellosis. Since the beginning of the 20th century, numerous studies have been conducted to develop an ideal vaccine for brucellosis (15). So far, different types of vaccines, such as live and inactivated vaccines, have been used to control brucellosis. Despite their success, these vaccines have several disadvantages (16).

Currently, a third generation of vaccines, which combines bioinformatics and genetics, is being developed. This generation leads to the production of safe vaccines that only use the immunogenic components of a pathogen (17, 18). The current project aimed to engineer and develop an HBHA-OMP25 vaccine against brucellosis using bioinformatics. In this case, the OMP25 protein of Brucella melitensis was considered an immunogenic protein. Outer membrane protein 25 is an outer membrane protein with a molecular weight of 25 kDa found in all species of Brucella (4). It has been reported that OMP25 plays an important role in triggering the pathogenicity and viability of B. melitensis in its host (19). Bowden et al. (20) and Commander et al. (21) reported that using OMP25 as a recombinant protein could lead to resistance against brucellosis in mice. Moreover, Yousefi et al. (19) confirmed that OMP25 can be used as an effective vaccine candidate against brucellosis.

In the current study, besides OMP25, the HBHA protein was used as an adjuvant to enhance the potency of the engineered chimera. The HBHA protein is a TLR4 receptor antagonist known to strengthen immune responses, making it useful as a molecular adjuvant (22). As mentioned in this study, the EAAAK linker was used to assemble the HBHA-OMP25 chimera. Fusion protein linkers play a key role in maintaining the functions of chimeric proteins. EAAAK is a rigid linker that preserves the functionality of different moieties within a chimeric protein (23).

In continuation, the physical and chemical features of the engineered chimera were investigated. Our analysis revealed that the molecular weight of the protein was 45.172 kDa. It has been confirmed that proteins with a molecular weight greater than 10 kDa can evade the renal system of the kidney, consequently increasing their half-life in the host (17). The instability index is an important parameter related to the stability of a protein. In general, a protein is considered stable when the instability index is less than 40, whereas proteins with an index greater than 40 are considered unstable. Our results showed that the instability index of our engineered chimera was 34.95, indicating it is a stable protein (24, 25).

Additionally, antigenicity results showed that our engineered chimera, with an antigenicity score of 0.7184, could stimulate the immune system. The 3D structure analysis of our engineered chimera revealed that 89% of amino acid residues were located in the core region, indicating appropriate 3D structure modeling. Our protein-protein docking results demonstrated that the HBHA moiety of the engineered chimera could interact with the TLR4/MD2 receptor with the lowest energy of -0.777 kcal/mol and nine hydrogen bonds, as shown in Table 1. To express our engineered chimera in a prokaryotic expression system, the codons of its nucleotide sequence were adapted for this system. Our results showed that the CAI (the most important indicator) of the adapted sequence improved from 0.375 to 0.970. It must be noted that the maximum value for CAI is 1; therefore, when CAI reaches 1, all codons are fully adapted (10).

Table 1.

List of Amino Acid Residues Involved in Hydrogen Bonds of Heparin-Binding Hemagglutinin Adhesin and Toll-Like Receptor/MD2 Interaction

HBHA MoietyTLR4/MD2 ReceptorLength of the Bond
Isoleucine 273Serine 1162.84
Phenylalanine 417Serine 1162.73
Lysine 274Glutamic acid 642.54
Lysine 274Glutamic acid 642.58
Lysine 274Asparagine 1142.67
Serine 272Histidine 1642.70
Aspartic acid 9Tryptophan 232.94
Alanine 20Glutamic acid 493.02
Alanine 16Phenylalanine 243.31

Abbreviations: HBHA, heparin-binding hemagglutinin adhesin; TLR4, toll-like receptor 4.

5.1. Conclusions

In this study, based on the OMP25 of B. melitensis and the HBHA protein (as a molecular adjuvant), a novel chimera of HBHA-OMP25 was designed. Our various results confirmed that the designed chimera could not only be an appropriate RV against brucellosis but also be easily produced using a prokaryotic expression system.

Acknowledgements

References

  • 1.

    Yousefi S, Tahmoorespur M, Sekhavati MH. B and T‐Cell Epitope Prediction of the OMP25 Antigen for Developing Brucella melitensis Vaccines for Sheep. Iran J Appl Anim Sci. 2015;5(3):639-8.

  • 2.

    Luna-Martinez JE, Mejia-Teran C. Brucellosis in Mexico: current status and trends. Vet Microbiol. 2002;90(1-4):19-30. [PubMed ID: 12414130]. https://doi.org/10.1016/s0378-1135(02)00241-9.

  • 3.

    Zowghi E, Ebadi A. Typing of brucella strains isolated in Iran. Arch Inst Razi. 1982;33(1):109-14. https://doi.org/10.22092/ari.1982.108881.

  • 4.

    Golshani M, Rafati S, Dashti A, Gholami E, Siadat SD, Oloomi M, et al. Vaccination with recombinant L7/L12-truncated Omp31 protein induces protection against Brucella infection in BALB/c mice. Mol Immunol. 2015;65(2):287-92. [PubMed ID: 25723468]. https://doi.org/10.1016/j.molimm.2015.01.009.

  • 5.

    Rashidian E, Forouharmehr A, Nazifi N, Jaydari A, Shams N. Computer-Aided Design of a Novel Poly-Epitope Protein in Fusion with an Adjuvant as a Vaccine Candidate Against Leptospirosis. Curr Proteomics. 2021;18(2):113-23. https://doi.org/10.2174/1570164617666200319144331.

  • 6.

    Bazgir M, Shams N, Forouharmehr A, Jaydari A, Nazifi N. [Engineering a Recombinant Chimeric Protein Consisting of OmpL Antigen Conjugated to HBHA Molecular Adjuvant (HBHA-OmpL) with the Aim of Developing a Vaccine Against Salmonella Typhimurium]. Vet Res Biol Products. 2023;36(2):100-11. Persian. https://doi.org/10.22092/vj.2022.359453.1995.

  • 7.

    Cloeckaert A, Vizcaino N, Paquet JY, Bowden RA, Elzer PH. Major outer membrane proteins of Brucella spp.: past, present and future. Vet Microbiol. 2002;90(1-4):229-47. [PubMed ID: 12414146]. https://doi.org/10.1016/s0378-1135(02)00211-0.

  • 8.

    Rashidian E, Shakarami Gandabeh Z, Forouharmehr A, Nazifi N, Shams N, Jaydari A. Immunoinformatics Approach to Engineer a Potent Poly-epitope Fusion Protein Vaccine Against Coxiella burnetii. Int J Pept Res Ther. 2020;26(4):2191-201. https://doi.org/10.1007/s10989-019-10013-6.

  • 9.

    Lei Y, Shao J, Ma F, Lei C, Chang H, Zhang Y. Enhanced efficacy of a multi-epitope vaccine for type A and O foot‑and-mouth disease virus by fusing multiple epitopes with Mycobacterium tuberculosis heparin-binding hemagglutinin (HBHA), a novel TLR4 agonist. Mol Immunol. 2020;121:118-26. [PubMed ID: 32199211]. https://doi.org/10.1016/j.molimm.2020.02.018.

  • 10.

    Carbone A, Zinovyev A, Kepes F. Codon adaptation index as a measure of dominating codon bias. Bioinformatics. 2003;19(16):2005-15. [PubMed ID: 14594704]. https://doi.org/10.1093/bioinformatics/btg272.

  • 11.

    Pappas G, Papadimitriou P, Akritidis N, Christou L, Tsianos EV. The new global map of human brucellosis. Lancet Infect Dis. 2006;6(2):91-9. [PubMed ID: 16439329]. https://doi.org/10.1016/S1473-3099(06)70382-6.

  • 12.

    D'Anastasio R, Zipfel B, Moggi-Cecchi J, Stanyon R, Capasso L. Possible brucellosis in an early hominin skeleton from sterkfontein, South Africa. PLoS One. 2009;4(7). e6439. [PubMed ID: 19649274]. [PubMed Central ID: PMC2713413]. https://doi.org/10.1371/journal.pone.0006439.

  • 13.

    Estein SM, Cassataro J, Vizcaino N, Zygmunt MS, Cloeckaert A, Bowden RA. The recombinant Omp31 from Brucella melitensis alone or associated with rough lipopolysaccharide induces protection against Brucella ovis infection in BALB/c mice. Microbes Infect. 2003;5(2):85-93. [PubMed ID: 12650766]. https://doi.org/10.1016/s1286-4579(02)00075-8.

  • 14.

    Olsen SC, Stoffregen WS. Essential role of vaccines in brucellosis control and eradication programs for livestock. Expert Rev Vaccines. 2005;4(6):915-28. [PubMed ID: 16372886]. https://doi.org/10.1586/14760584.4.6.915.

  • 15.

    Nicoletti P. Vaccination against Brucella. Adv Biotechnol Processes. 1990;13:147-68. [PubMed ID: 2185782].

  • 16.

    Avila-Calderon ED, Lopez-Merino A, Sriranganathan N, Boyle SM, Contreras-Rodriguez A. A history of the development of Brucella vaccines. Biomed Res Int. 2013;2013:743509. [PubMed ID: 23862154]. [PubMed Central ID: PMC3686056]. https://doi.org/10.1155/2013/743509.

  • 17.

    Jaydari A, Nazifi N, Forouharmehr A. Computational design of a novel multi-epitope vaccine against Coxiella burnetii. Hum Immunol. 2020;81(10-11):596-605. [PubMed ID: 32718721]. https://doi.org/10.1016/j.humimm.2020.05.010.

  • 18.

    Nazifi N, Tahmoorespur M, Sekhavati MH, Haghparast A, Behroozikhah MA. Engineering, Cloning and Expression of DNA Sequence Coding of OMP31 Epitope of Brucella melitensis linked to IL-2 in Escherichia coli. Int J Infect. 2018;5(3). e68974. https://doi.org/10.5812/iji.68974.

  • 19.

    Yousefi S, Tahmoorespur M, Sekhavati MH. Cloning, expression and molecular analysis of Iranian Brucella melitensis Omp25 gene for designing a subunit vaccine. Res Pharm Sci. 2016;11(5):412-8. [PubMed ID: 27920824]. [PubMed Central ID: PMC5122831]. https://doi.org/10.4103/1735-5362.192493.

  • 20.

    Bowden RA, Cloeckaert A, Zygmunt MS, Dubray G. Evaluation of immunogenicity and protective activity in BALB/c mice of the 25-kDa major outer-membrane protein of Brucella melitensis (Omp25) expressed in Escherichia coli. J Med Microbiol. 1998;47(1):39-48. [PubMed ID: 9449948]. https://doi.org/10.1099/00222615-47-1-39.

  • 21.

    Commander NJ, Spencer SA, Wren BW, MacMillan AP. The identification of two protective DNA vaccines from a panel of five plasmid constructs encoding Brucella melitensis 16M genes. Vaccine. 2007;25(1):43-54. [PubMed ID: 17049676]. https://doi.org/10.1016/j.vaccine.2006.07.046.

  • 22.

    Forouharmehr A. Engineering an efficient poly-epitope vaccine against Toxoplasma gondii infection: A computational vaccinology study. Microb Pathog. 2021;152:104646. [PubMed ID: 33242641]. https://doi.org/10.1016/j.micpath.2020.104646.

  • 23.

    Chen X, Zaro JL, Shen WC. Fusion protein linkers: property, design and functionality. Adv Drug Deliv Rev. 2013;65(10):1357-69. [PubMed ID: 23026637]. [PubMed Central ID: PMC3726540]. https://doi.org/10.1016/j.addr.2012.09.039.

  • 24.

    Gamage DG, Gunaratne A, Periyannan GR, Russell TG. Applicability of Instability Index for In vitro Protein Stability Prediction. Protein Pept Lett. 2019;26(5):339-47. [PubMed ID: 30816075]. https://doi.org/10.2174/0929866526666190228144219.

  • 25.

    Guruprasad K, Reddy BV, Pandit MW. Correlation between stability of a protein and its dipeptide composition: a novel approach for predicting in vivo stability of a protein from its primary sequence. Protein Eng. 1990;4(2):155-61. [PubMed ID: 2075190]. https://doi.org/10.1093/protein/4.2.155.