Abstract
Monkeypox (mpox) is a zoonotic infectious disease caused by the mpox virus. Mpox symptoms are similar to smallpox with less severity and lower mortality. As yet mpox virus is not characterized by as high transmissibility as some severe acute respiratory syndrome 2 (SARS-CoV-2) variants, still, it is spreading, especially among men who have sex with men (MSM). Thus, taking preventive measures, such as vaccination, is highly recommended. While the smallpox vaccine has demonstrated considerable efficacy against the mpox virus due to the antigenic similarities, the development of a universal anti-mpox vaccine remains a necessary pursuit. Recently, nucleic acid vaccines have garnered special attention owing to their numerous advantages compared to traditional vaccines. Importantly, DNA vaccines have certain advantages over mRNA vaccines. In this study, a potentially universal DNA vaccine candidate against mpox based on conserved epitopes was designed and its efficacy was evaluated via an immunoinformatics approach. The vaccine candidate demonstrated potent humoral and cellular immune responses in silico, indicating the potential efficacy in vivo and the need for further research.
1. Introduction
The history of pandemics has shown us that not a single virus outbreak should be neglected as sooner or later the re-emergence of the same or even more virulent strain may occur with severe outcomes. The spread of the monkeypox (mpox) virus was considered limited to Central and West African countries until May 2022 when it crossed the continental borders and gained global health emergency status by World Health Organization (WHO) on 23 July 2022 (Papukashvili et al., 2022b; World Health Organization, 2022b). Mpox is a zoonotic disease with similar clinical manifestation as the smallpox virus which also belongs to the genus Orthopoxvirus with an additional symptom of lymph node enlargement. Mpox disease is characterized by a rash appearing 1–3 days following the onset of fever and lymphadenopathy. The rash is usually distributed on the whole body, especially on extremities (Sklenovská and Van Ranst, 2018), genitals, and anus (Kmiec and Kirchhoff, 2022; Mileto et al., 2022; Moschese et al., 2022; Sah et al., 2022; World Health Organization, 2022a). The mpox virus is a large, ~280 nm X ~220, brick/oval-shaped virus with a large linear double-stranded ~197 kb long DNA genome (Papukashvili et al., 2022b) with more than 190 open reading frames (ORFs) (Shchelkunov et al., 2002; Kugelman et al., 2014; Kmiec and Kirchhoff, 2022; Vandenbogaert et al., 2022; Zhu et al., 2022). The mpox virus produces two infectious forms of virions from the infected cells—an intracellular mature virion (IMV) and an extracellular enveloped virion (EEV) with an extra envelope. The structure of the mpox virus is complex, and many of the viral antigens are not well-studied. Infection by EEV is more efficient than by IMV (Locker et al., 2000).
One of the main reasons for the recent mpox multi-country outbreaks is the cessation of smallpox vaccination in 1980 which seems reasonable as most of the current mpox cases are detected in smallpox-unvaccinated people (Papukashvili et al., 2022b). Although the smallpox vaccine which is based on the vaccinia virus, another representative of Orthopoxvirus, is effective against the mpox virus (Zandi et al., 2023), there is no specific mpox vaccine available till now. The modified vaccinia virus Ankara-Bavarian Nordic (MVA-BN), also known as Imvamune, Jynneos, or Imvanex, is a third-generation authorized smallpox vaccine and is also used as the mpox vaccine (Zaeck et al., 2022). Although the currently available vaccine is effective for the mpox virus (CDC, 2022; Kandeel et al., 2023), Zaeck et al. have demonstrated that immunization series with the two-shot MVA-BN vaccine, in non-primed individuals, yields relatively low levels of mpox virus-neutralizing antibodies (Zaeck et al., 2022). The line of preclinical studies has made efforts in the development of the mpox vaccine. Hooper et al. have shown that the DNA vaccine encoding the antigens L1, A27, A33, and B5R protected rhesus macaques from the lethal challenge of the mpox virus (Hooper et al., 2004). Other studies have provided evidence of protection from the mpox virus challenge in non-human primates after immunization with A27, A33, B5, and L1 proteins (Buchman et al., 2010; Heraud et al., 2022). However, these vaccines were strain-specific. Importantly, recently worldwide spread mpox virus is of the West African clade with a mortality rate of approximately 1–3.6% (Kmiec and Kirchhoff, 2022; Kozlov, 2022; Yang, 2022), while the mpox of Central African clade is deadlier with a mortality rate of up to 11% (Jezek et al., 1987; Shafaati and Zandi, 2022; Yang, 2022). Considering the possible severe outcomes of whether the mpox virus of the Central African clade is spread, developing a new universal vaccine against mpox is urgently needed. Nucleic acid vaccines have gained enormous attention in recent years due to their efficacy, safety, cost-effectiveness, and time-saving features compared with conventional approaches. Apparently, the available mRNA vaccines for the current coronavirus disease 2019 (COVID-19) have high efficacy and safety profiles (Pardi et al., 2018; Corbett et al., 2020; Walsh et al., 2020; Haas et al., 2021; Kowalzik et al., 2021; Liu et al., 2022b). On the other hand, studies on DNA vaccines demonstrated no less advantageous features than mRNA vaccines including the requirement of even less production time and high stability (Flingai et al., 2013; Williams, 2013; Chavda et al., 2021b; Xia et al., 2021; Wang et al., 2022). More importantly, there is already an available DNA vaccine for COVID-19 approved in India (Sheridan, 2021). According to the current situation of available vaccines and the general epidemiological picture, developing a universal next-generation effective anti-mpox vaccine is necessary. To design a potentially universal DNA vaccine contender against mpox virus, in the present research, according to their functions and immunogenicity (Shchelkunov et al., 2002; Hooper et al., 2004; Buchman et al., 2010; Hirao et al., 2011; Meng et al., 2018; Zhang et al., 2023), sequences of four antigen proteins (A5L, A15L, A35R, and B6R) were retrieved and aligned, and the conserved sequences were selected to predict B-cell, MHC-I, and MHC-II binding epitopes. Except for the importance of immune response and protection against viral challenges, another rationale behind choosing these representative antigens was to assess the effectiveness of the novel antigen combination, aiming to propose another candidate for an anti-mpox vaccine. After the design of the vaccine, the final construct was optimized, and the structure and various characteristics of the vaccine were predicted. Immune simulation analyses demonstrated the strong humoral and cellular responses that warrant the contribution of the vaccine candidate in potentially universal multi-epitope DNA vaccine development.
2. Materials and methods
2.1. Selection of antigens, collection of data, and selection of conserved sequences
Four antigen proteins, namely, A5L, A15L, A35R, and B6R, were chosen according to their function in the life cycle and immunogenicity of the mpox virus. This selection also aimed to assess the immune responses generated by the vaccine against this combination of antigens. They are the orthologs of vaccinia virus glycoproteins—A4L, A14L, A33R, and B5R, respectively (Shchelkunov et al., 2002). These antigens are highly conserved among the mpox, vaccinia, and variola viruses. The protein sequences of mpox virus antigens reported on the National Center for Biotechnology Information (NCBI) database were retrieved. The number of downloaded sequences for A5L, A15L, A35R, and B6R was 1819, 355, 349, and 376, respectively. All the sequences of A5L were used for further analyses as they met the selection criteria such as the known length of the protein sequence. For the rest of the three proteins, the redundant sequences were removed according to their length and quality−179 sequences of A15L, 178 sequences of A35R, and 178 sequences of B6R were used for further analyses. The software Jalview 2.11.1.4 was used for the alignment and conservancy analysis (Waterhouse et al., 2009). The threshold filter of conservation was adjusted to 10 below the threshold, and the consensus sequence was used for B-cell and T-cell (MHC-I and MHC-II binding) epitope prediction.
2.2. Epitope prediction and selection
NetMHCpan EL 4.1 server was applied for the prediction of cytotoxic T lymphocyte (CTL) epitopes (https://services.healthtech.dtu.dk/service.php?NetMHCpan-4.1) (Reynisson et al., 2020). This server predicts the binding of peptides to any MHC-I molecule of the known sequence using artificial neural networks (ANNs) that are a collection of simple interconnected algorithms processing information according to the external input. The globally most prevalent alleles HLA-A*01:01, HLA-B*07:02, and HLA-B*08:01 were selected from the “Allele Frequency Net Database” and used for the prediction of the MHC-I binding epitope. The default parameters of the server were used: the threshold for strong binders was set to 0.5 (% rank), while for weak binders it was 2%. The predicted peptides were identified as strong binders when the % rank was below the set threshold (0.5%), while the weak binders were considered when the peptide rank was above 0.5% and below 2%. For the prediction of helper T lymphocyte (HTL) epitopes, the NetMHCIIpan-4.0 server which performs the prediction of peptide binding to any MHC-II molecule of the known sequence using ANNs was used. The worldwide most prevalent MHC-II alleles DRB1*07:01 and DRB1*15:01 were selected for the prediction of MHC-II binding peptides of the mpox virus conserved sequences. According to the default parameters, the predicted peptides were identified as strong binders when the % rank was below the set threshold (1%), while the weak binders were considered when the peptide rank was above 1% and below 2%. The epitopes with the strongest binding capacity to MHC-I and MHC-II alleles were selected. Both NetMHCpan-4.1 and NetMHCIIpan-4.0 use the NNAlign_MA machine learning framework, which has ANNs at its core. NetMHCpan-4.1, a computational tool, has been trained using a comprehensive dataset comprising over 850,000 peptides with quantitative binding affinity (BA) and mass spectrometry-eluted ligand (EL) information. The BA data utilized in the training of NetMHCpan-4.1 encompass 170 MHC molecules from various species—humans (HLA-A, HLA-B, HLA-C, and HLA-E), mice (H-2), cattle (BoLA), primates (Patr, Mamu, and Gogo), equine (Eqca), and swine (SLA). The EL data used in the analysis include 177 MHC molecules from a range of species, such as humans (HLA-A, B, C, E), mice (H-2), cattle (BoLA), primates (Patr, Mamu, and Gogo), swine (SLA), equine (Eqca), and dogs (DLA). NetMHCIIpan-4.0 has been trained using a comprehensive dataset comprising more than 500,000 measurements of BA and EL. This dataset covers the three human MHC class II isotypes, namely HLA-DR, HLA-DQ, and HLA-DP, along with the mouse molecules (H-2). The inclusion of EL data expands the coverage of MHC-II molecules, as the BA data cover 59 molecules, while the EL data cover 74 molecules. Both NetMHCpan-4.1 and NetMHCIIpan-4.0 are the currently recommended algorithms in Immune Epitope Database (IEDB). They demonstrate cutting-edge performance and surpass other algorithms in benchmark experiments (Reynisson et al., 2020).
For predicting the linear B-cell epitopes, the IEDB with the method of Bepipred linear epitope prediction tool (v2.0) was utilized (Jespersen et al., 2017). The BepiPred-2.0 server performs the prediction of B-cell epitopes from a protein sequence via a random forest algorithm trained on epitopes and non-epitope amino acids determined from crystal structures. The threshold was set at 0.6 meaning that the residues scoring above the threshold were predicted to be part of the epitopes. The B-cell epitopes were then selected.
2.3. Multi-epitope DNA vaccine design
After the prediction of B-cell, MHC-I, and MHC-II epitopes, the epitopes with the most optimal features were selected. Multi-peptide vaccine construct was designed and epitopes of A5L, A15L, A35R, and B6R were fused via the appropriate linkers. The linkers play a crucial role in multi-epitope vaccine development as their reasonable selection allows to enhance structural stability, flexibility, and proper folding, as well as to increase the immune response by allowing to include multiple epitopes in a single vaccine construct. Thus, the well-described linkers KK, GPGPG, and EAAAK were used to connect different peptide components. EAAAK is a rigid linker that provides an alpha helix-forming structure between the domains that enhances the stability and maintains the constant spacing between the domains of the fusion protein as well as preserves their functions (Chen et al., 2013). B-cell epitopes were linked with a flexible KK linker (Rahmani et al., 2021; Tarrahimofrad et al., 2021). The KK linker allows the presentation of each linked epitope to antibodies while avoiding antibody induction against the whole joined sequences resulting in the antibody reactivity to each of the B-cell epitopes (Yano et al., 2005). MHC-I and MHC-II epitopes were conjugated together with another flexible linker GPGPG which is a glycine-rich linker, in addition enhancing the solubility of vaccine construct, and also provides high accessibility and flexibility for the adjacent epitopes (Dong et al., 2020; Khan et al., 2021; Martinelli, 2022). Except for the epitopes, to rationally optimize the fused multi-peptide construct for the optimal immunogenic outcome, adjuvants were also fused. Particularly, 45 amino acid long human β-defensin 3 (hBD3) with immunomodulatory effect (Hoover et al., 2003; Ali et al., 2017; Qamar et al., 2020) was adjoined to the N-terminal of the vaccine after the leading sequence of tissue plasminogen activator which also facilitates the antigen presentation (Ahammad and Lira, 2020). hBD3 is one of the top five adjuvants utilized for COVID-19 subunit vaccines (Mekonnen et al., 2022). Furthermore, defensins play a crucial role in defending against pathogen infections as they effectively bridge the innate and adaptive immune responses through leukocytes such as dendritic cells (DCs) and T cells (Bellamkonda et al., 2022). The rigid linker EAAAK (Nezafat et al., 2014; Tarrahimofrad et al., 2021) was used to conjugate the hBD3 with the pan-HLA-DR-binding epitope (PADRE) which triggers antigen-specific CD4+ T cells (Alexander et al., 2000). EAAAK was used to connect the PADRE sequence to the rest of the construct as well. For future in vitro validation of the expression of the multi-epitope construct, the C-terminal polyhistidine tag-−6xHis—was also adjoined (Boyoglu-Barnum et al., 2021). The untranslated regions (UTRs) were included for ensuring the proper stability and translation of the multi-epitope construct when it is transcribed into mRNA and for allowing the translation efficiency and stability of mRNA while validating the protein expression in vitro and in vivo. The following UTRs were used: cytomegalovirus (HCMV) immediate early (IE) partial sequence as 5′UTR, and human growth hormone partial sequence as a 3′UTR as they have demonstrated the effective stimulation of the protein production in vivo (Rybakova et al., 2019).
2.4. Vaccine properties—In silico assessment
The vaccine should have wide human population coverage. Thus, the IEDB population coverage analysis tool was employed to evaluate the worldwide human population coverage of the designed DNA vaccine (Bui et al., 2006). The allergenicity of the designed vaccine was predicted via the AllerCatPro server which predicts the allergenic proteins based on the similarity of both their amino acid sequences and 3D structures compared with the comprehensive dataset of already known allergens from the WHO/International Union of Immunological Societies, Comprehensive Protein Allergen Resource, Food Allergy Research and Resource Program, UniProtKB, and Allergome (Maurer-Stroh et al., 2019; Nguyen et al., 2022). VaxiJen v2.0 was applied to predict the antigenicity of the multi-epitope protein construct (Doytchinova and Flower, 2007a,b). The threshold was set as 0.4, and default parameters were employed (Ghafouri et al., 2022).
The multi-epitope vaccine physico-chemical properties were assessed by ProtParam (Gasteiger et al., 2005). Molecular weight (MW), the composition of atoms and amino acids, instability, estimated half-life, theoretical isoelectric point (pI), aliphatic indexes, and grand average of hydropathicity (GRAVY) of the vaccine construct were evaluated. The MW was determined by summing amino acids' average isotopic masses and the average isotopic mass of a single water molecule. pI was calculated by amino acids' pKa value. Half-life estimation describes the necessary time from protein synthesis until the disappearance of its half-amount from the cell. The instability index was used to determine how stable is the protein in a test tube. The value <40 estimates the protein as stable. GRAVY calculates the hydropathicity index of all amino acids divided by the length of the sequence. The larger number denotes more hydrophobicity of amino acids (Gasteiger et al., 2005). An aliphatic index is used to characterize the relative volume of protein which is occupied by aliphatic side chains and is considered a positive factor in increasing thermostability. The solubility of the vaccine protein was predicted via DeepSoluE. This server is a deep-learning predictor, and its prediction algorithm outperforms the other servers for the prediction of protein solubility (Wang and Zou, 2023).
2.5. Prediction and quality assessment of the vaccine's 3D structure
The 3D structure of the designed vaccine was constructed via RoseTTAFold (Baek et al., 2021). This tool uses a three-track network: integrating protein sequence patterns, amino acid interactions, and 3D structure to accurately predict protein structure and interaction (Baek et al., 2021). After the 3D structure of the vaccine protein was predicted, it was optimized with GalaxyRefine and FG-MD which is a molecular dynamics (MD)-based algorithm for protein structure refinement at the atomic level (Zhang et al., 2011; Heo et al., 2013). This immunoinformatics approach uses molecular dynamics simulation and relaxes and refines the protein structure (Heo et al., 2013). Predicting and knowing the potential structure is important as it determines the function of the protein that allows its application. As the predicted 3D structure provides the understanding of its interactions with the immune system, it is then used for further analyses, e.g., prediction of B-cell discontinuous epitopes and molecular docking. The tertiary structure quality of the 3D protein was then verified with the Ramachandran plot and ERRAT. The Ramachandran plot illustrates energetically permissible regions for the backbone dihedral angles ψ and ϕ of amino acid residues within a protein structure allowing the assessment of 3D protein structure quality (Nelson et al., 2021). ERRAT also represents the quality but for non-bonded interactions, and a higher score of ERRAT indicates a higher quality of the tertiary structure of the protein.
2.6. Presence of conformational B-cell epitopes in the vaccine
To predict the discontinuous B-cell epitopes, ElliPro was employed. This tool is the recourse of the immune epitope database (IEDB). It allows for predicting the conformational B-cell epitopes according to the antigen's tertiary structure. ElliPro associates predicted discontinuous epitopes with protrusion index (PI) which is averaged over epitopes' residues. The approximation of the antigen's tertiary structure is achieved via the number of ellipsoids. The PI score is defined based on the residue's mass center which is lying outside the largest ellipsoid. A higher value suggests higher solvent accessibility which, on the other hand, is crucial in protein stability and folding (Ponomarenko et al., 2008).
2.7. In silico immunization
The C-ImmSim server was utilized to analyze the simulated humoral and cellular immune responses elicited against the multi-epitope vaccine protein. After submitting the sequence of vaccines in the FASTA format, the server predicts the immune responses via a position-specific scoring matrix and machine learning algorithms (Rapin et al., 2011). Two HLA-A (01:01 and 02:01), two HLA-B (07:02 and 08:01), and two DRB (07:01 and 15:01) were selected through Allele Frequency Net Database (AFND). The prime, second, and booster doses of the vaccine were injected with 4 weeks gap. The volume and steps of the simulation were 10 and 270, respectively. The vaccine did not contain LPS. The random seed was 12,345, and the injection time periods were 1, 85, and 169 (Rapin et al., 2010).
2.8. Docking of vaccine antigen with host receptors
To evaluate the binding capacity of the vaccine antigen to its recognition receptors of the host, molecular docking analysis was conducted by the server of ClusPro. ClusPro determines the molecular docking of two proteins providing an outcome via presenting putative protein complexes in a ranked list. Ligand conformation, orientation, and position along with the assessment of binding affinity are the major properties that determine molecular docking. Ultimately, the electrostatic-favored protein–protein docked complexes with favorable desolvation-free energies are picked (Kozakov et al., 2017). Analysis of vaccine antigen's docking with MHC-I (HLA-A*02:01) (6TDS) and MHC-II (HLA-DRB1*01:01) (1AQD) host receptors was performed. The MHC molecules were used to dock the following epitopes in the vaccine construct: MHC-I ligands—ETNDLVTNVY (A5L), SPTRTWKVL (A15L), DSDVSQEVRKY (A35R), and VSDYVSELY (B6R); MHC-II ligands—APILLPSSTAPVLKP (A5L) and HSDYKSFEDAKANCA (A35R). The tertiary structure of these epitopes was predicted with AlphaFold2 (Jumper et al., 2021; Mirdita et al., 2022). The protein databank (PDB) files 6TDS (MHC-I) and 1AQD (MHC-II) were edited and cleaned to remove heteroatoms, bound peptides, and water molecules. The server of ClusPro performs the docking of rigid body via sampling billions of conformations, clusters 1000 lowest energy structures to find the most favorable models, and refines selected structures via minimization of energy (Kozakov et al., 2017).
3. Results
3.1. Design of the final construct
All the predicted B-cell, MHC-I, and MHC-II epitopes are given in Supplementary materials 1–12. In total, 30 epitopes were used for the final vaccine construct design. The open reading frame (ORF) of the vaccine construct is shown in Supplementary material. The selected four antigens of the mpox virus are shown in Figure 1. The location, functions, and immune properties of the selected antigens are listed in Table 1. The strategy of in silico design and evaluation of a potentially universal DNA vaccine for mpox are shown in Figure 2. All the epitopes included in the proposed multi-epitope vaccine construct along with the DNA vaccine are given in Figure 3.
Figure 1
Table 1
| Antigen name | Location | Function/properties | Ref. |
|---|---|---|---|
| A5L | IMV | Immunodominant virion core protein (281 aa) needed for the assembly and disassembly of virion; A highly antigenic protein of the viral core; Enhances the effect of cytotoxic T lymphocytes in the mpox challenge model; Similar to VACV A4L protein. | (Maa and Esteban, 1987; Shchelkunov et al., 2002; Hirao et al., 2011) |
| A15L | IMV | Essential inner membrane protein (90 aa) of IMV; Immunodominant antigen in smallpox vaccine; Similar to VACV A14L protein. | (Shchelkunov et al., 2002; Meng et al., 2018) |
| A35R | EEV | Envelope glycoprotein (181 aa); Cell-to-cell spread of virion; Neutralizing antibody target; Similar to VACV A33R protein | (Shchelkunov et al., 2002; Hirao et al., 2011) |
| B6R | EEV | Palmitoylated glycoprotein (317 aa); Required for efficient cell spread; Complement control; Similar to VACV B5R protein | (Shchelkunov et al., 2002; Hirao et al., 2011) |
Selected antigens of the mpox virus and their functions.
Mpox, monkeypox; IMV, intracellular mature virus; EEV, extracellular enveloped virus.
Figure 2
Figure 3
The population coverage analyses demonstrated that the 5 MHC-I alleles (HLA-A*01:01, HLA-A*02:01, HLA-B*07:02, HLA-B*08:01, and HLA-A*02:12) with IC50 binding affinity <500 nM with the MHC-I epitopes and the 13 MHC-II alleles (HLA-DRB1*07:01, HLA-DRB1*01:01, HLA-DRB1*09:01, HLA-DRB1*10:01, HLA-DRB1*04:01, HLA-DRB1*15:01, HLA-DRB1*14:01, HLA-DRB1*13:01, HLA-DRB1*08:01, HLA-DRB1*03:01, HLA-DRB1*11:01, HLA-DRB1*12:01, and HLA-DRB1*16:01) with IC50 binding affinity <500 nM with MHC-II epitopes cover 95.21% of the human population, which means that these epitopes could cover approximately 95% of the population from different regions of the world (Supplementary Figure 1 and Supplementary Table 1).
3.2. Analyses of protein structure and assessment of physico-chemical characteristics
The predicted 3D structure of the vaccine construct along with the quality evaluation is given in Figure 4. The distribution of amino acid residues in the Ramachandran plot was as follows: refined structure had 279 residues in the most favored region, 93 residues in the additionally allowed region, 13 in the generously allowed region, and only 7 residues in the disallowed region. The Ramachandran plot allows us to visualize the energetically allowed and disallowed regions for the dihedral angles. For example, in poor quality homology models, many dihedral angles can be found in the forbidden regions of the plot which indicates the problems with the structure. The favored regions in the Ramachandran plot correspond to the regular secondary structures. The ERRAT score was 88.7218 meaning that all these parameters indicate the favorable overall quality of the vaccine antigen protein.
Figure 4
The instability index of the proposed vaccine construct was 29.12, which classifies the vaccine as stable. The aliphatic index was 60.22 which indicates the high thermostability of the vaccine protein. GRAVY was −0.678, stipulating hydrophilicity. The half-life in vitro was estimated to be 30 h, while the in vivo in the yeast and E. coli were estimated to be over 20 h and over 10 h, respectively. The number of amino acids was 542, MW was 56.55076 kDa, and the theoretical pI was 9.28. The vaccine construct contains 46 negatively charged amino acid residues (Asp+Glu) and 64 positively charged amino acid residues (Arg+Lys). The probability of the vaccine solubility was 0.7963, indicating the solubility of the vaccine protein (Wang and Zou, 2023).
3.3. Presence of discontinuous B-cell epitopes
All three conformational epitopes along with the scores are listed in Table 2. The locations of each discontinuous epitope in the tertiary structure of the antigen protein are illustrated in Figure 5.
Table 2
| # | Epitopes | Residue number | Score |
|---|---|---|---|
| 1 | A:G1, A:I2, A:I3, A:N4, A:T5, A:L6, A:Q7, A:K8, A:Y9, A:Y10, A:C11, A:V13, A:R14, A:G16, A:R17, A:C18, A:A19, A:V20, A:L21, A:S22, A:C23, A:L24, A:E28, A:Q29, A:I30, A:G31, A:K32, A:C33, A:S34, A:T35, A:R36, A:G37, A:R38, A:K39, A:C40, A:C41, A:R42, A:K109, A:V110, A:V111, A:S112, A:S113, A:T114, A:T115, A:Q116, A:Y117, A:D118, A:H119, A:K120, A:K121, A:K122, A:G123, A:S124, A:D125, A:G126, A:N127, A:P128, A:I129, A:T130, A:K131, A:T132, A:T133, A:S134, A:D135, A:Y136, A:Q137, A:D138, A:S139, A:D140, A:V141, A:S142, A:Q143, A:E144, A:V145, A:R146, A:K147, A:K148, A:S149, A:E150, A:L151, A:Y152, A:D153, A:K154, A:P155, A:L156, A:K157, A:K158, A:P159, A:D160, A:D161, A:E162, A:T163, A:D164, A:L165, A:K167, A:L168, A:K170, A:G171, A:P172, A:G173, A:P174, A:G175, A:D179, A:P203, A:G204, A:Y472, A:H474, A:K475, A:E476, A:G477, A:Y485, A:S487, A:F488, A:G499, A:P500, A:G501, A:D502, A:S503, A:G504, A:Y505, A:H506, A:S507, A:L508, A:D509, A:P510, A:N511, A:A512, A:V513, A:C514, A:E515, A:T516, A:G517, A:P518, A:G519, A:P520, A:G521, A:M522, A:C523, A:T524, A:V525, A:S526, A:D527, A:Y528, A:V529, A:S530, A:E531, A:L532, A:Y533, A:D534, A:K535, A:P536, A:H537, A:H538, A:H539, A:H540, A:H541, A:H542 | 157 | 0.746 |
| 2 | A:A67, A:K68, A:Q69, A:S70, A:D71, A:G72, A:S73, A:I74, A:S75, A:C76, A:K78, A:T79, A:T80, A:A81, A:Q82, A:Q83, A:N84, A:P85, A:N86, A:P87, A:G88, A:A89, A:Q92, A:G272, A:P273, A:G274, A:S275, A:P276, A:T277, A:R278, A:T279, A:W280, A:K281, A:V282, A:G283, A:P284, A:G285, A:P286, A:G287, A:D288, A:S289, A:D290, A:V291, A:S292, A:Q293, A:E294, A:V295, A:R296, A:K297, A:Y298, A:G299, A:P300, A:G301, A:P302, A:G303, A:Y304, A:Q305, A:D306, A:S307, A:D308, A:V309, A:S310, A:Q311, A:E312, A:V313, A:G314, A:P315, A:G316, A:P317, A:G318, A:L319, A:P320, A:N321, A:K322, A:S323, A:D324, A:V325, A:L326, A:G327, A:P328, A:G329, A:P330, A:G331, A:N332, A:K333, A:R334, A:K335, A:R336, A:I338, A:G339, A:L340, A:G341, A:Y354, A:G355, A:P356, A:G357, A:V368, A:G369, A:P370, A:G371, A:P372, A:G373, A:K374, A:P375, A:G384, A:P385, A:G386, A:P387, A:G388, A:H394, A:K395, A:L396, A:G397, A:P398, A:G399, A:P400, A:G401, A:P403, A:D422, A:K429, A:E432, A:R433, A:Q435, A:K436, A:G437, A:P438, A:G439, A:P440, A:G446, A:G456, A:G457, A:P458, A:G459, A:P460 | 134 | 0.663 |
| 3 | A:P416, A:G417, A:P418 | 3 | 0.504 |
Predicted discontinuous B-cell epitopes.
Figure 5
3.4. Assessment of allergenicity, antigenicity, and molecular docking
A vaccine protein construct did not show evidence of allergenicity allowing us to presume that the multi-epitope construct can be considered with low allergenic potential. On the other hand, the multi-peptide was demonstrated to be a probable antigen with a score of 0.4248.
The molecular docking with MHC molecules demonstrated stable binding. ClusPro identifies the most probable complex models by determining the largest clusters available. The interaction positions of the docked complex of vaccine epitopes and MHC molecules are shown in Figure 6. The interaction energy is given in Table 3 where the cluster size is the size of the largest cluster found, while the lowest energy is the Gibbs free energy of the best complex conformation (Kozakov et al., 2017).
Figure 6
Table 3
| Peptide | Cluster size | Center weighted score | Lowest energy |
|---|---|---|---|
| Molecular docking of epitopes with 6TDS (MHC-I) | |||
| ETNDLVTNVY | 252 | −909.7 | −1,089.0 |
| SPTRTWKVL | 396 | −1,083.7 | −1,239.3 |
| DSDVSQEVRKY | 217 | −647.4 | −794.9 |
| VSDYVSELY | 344 | −1,194.4 | −1,194.4 |
| Molecular docking of epitopes with 1AQD (MHC-II) | |||
| APILLPSSTAPVLKP | 190 | −1,296.4 | −1,296.4 |
| HSDYKSFEDAKANCA | 244 | −845.1 | −1,023.0 |
Molecular docking of selected vaccine epitopes and MHC molecules.
3.5. Evaluation of immune response induced by immune simulation
In silico analyses demonstrated that the immune response triggered by the proposed vaccine was compatible with the responses generally induced via in vivo immunization. The immune responses after the additional two booster doses of vaccine were stronger compared to the prime immunization. The immunization lowered the antigen level, while high levels of antibodies (IgG and IgM) were produced. Immunoglobulin and antigen levels varied over time. The antigen abundance peaked at each injection time point (Figure 7A). The humoral response was induced after each shot, and the antibody levels remained elevated during the weeks after the last vaccine injection. The immunization stimulated the B cells (Figures 7B, C). The ignorable number of anergic cells and activation of CTLs was observed (Figure 7D). The counts of HTLs demonstrate that the duplication phase starts immediately after each injection (Figure 7E). Days after each injection and upon successful interaction with antigen-presenting cells, HTLs start to duplicate and differentiate into memory cells. HTLs also foster the secretion of cytokines and humoral response. High levels of IFN-γ, TGF-β, and IL-10 were induced upon each injection (Figure 7F). The other parameters of immune responses are shown in Supplementary Figure 2.
Figure 7
4. Discussion
Since the mpox virus crossed African borders and caused a multi-country outbreak (Supplementary Figure 3) with increased cases of human-to-human transmission, the global concern has increased; thus, there is a strong need for the development of specific drugs and vaccines. Although the available smallpox vaccine is effective for mpox prevention (Shafaati and Zandi, 2023), except for the painful immunization procedure, adverse effects were also reported (Maurer et al., 2003). Certain groups of the population are vulnerable to vaccination as serious side effects including pericarditis and myocarditis have been noted. A risk of recombination of genes between the mpox virus and attenuated or live poxvirus-based vaccines exists (Lum et al., 2022). In addition, after immunization with the smallpox vaccine, a manifestation of mpox disease still has been observed (Meyer et al., 2002). All the abovementioned disadvantages of the smallpox vaccine can be potentially overcome by the new DNA-based universal mpox vaccine candidate as, on the one hand, DNA vaccines are safe and simple to manufacture compared to their conventional counterparts. On the other hand, the mpox vaccine designed in this study has the potential to be effective against various strains of mpox viruses and smallpox virus since these two species share a high level of similarity. Noteworthily, the majority of the cases occur in men who have sex with men (MSM). Close physical contact is indeed crucial in transmission (Heskin et al., 2022; Martínez et al., 2022; Orviz et al., 2022). The sex-related infection tendency was always reported before the current multi-country outbreak (Sklenovská and Van Ranst, 2018).
The report of human-to-dog transmission of the mpox virus also complicates combating the outbreak (Seang et al., 2022). Additionally, the re-emergence risk of smallpox (mortality rate 10–75%) (Javier et al., 2022) or risk of spreading to the Central African clade of the mpox virus, which is much deadlier compared with the currently spread West African clade (Bunge et al., 2022), is very real and has a huge potential to cause severe outcomes. The recent multi-country outbreak also demonstrated that the mpox virus which belongs to the West African clade can be life-threatening (Meyer et al., 2002). All the abovementioned conditions indicate the urgency of designing a new, specific next-generation mpox vaccine that will be effective, safe, and easy to develop.
In general, nucleic acid-based approaches find ways through various applications such as vaccines against bacterial (Wang et al., 2023) and viral (Hobernik and Bros, 2018) infections, protein replacement (Papukashvili et al., 2022a; Vavilis et al., 2023), and cancer treatment (Hobernik and Bros, 2018; Liu et al., 2022a). In silico-designed and evaluated multi-peptide vaccines are thermodynamically stable, effective, specific, and easily and inexpensively developed compared with conventional vaccines.
In this study, four antigens of the mpox virus that are the target of immunity were selected according to their functions and role in immune response (Table 1). A5L is an immunodominant core protein that is necessary for the assembly and disassembly of the virion (Maa and Esteban, 1987; Shchelkunov et al., 2002). A4L, which is the ortholog of the A5L protein in the vaccinia virus, has been used as a part of the DNA vaccine along with the orthologs of A35R and B6R, among the other antigens and induced strong immune response in non-human primates, providing protection from the mpox virus (Hirao et al., 2011). A15L (similar to A14L in the vaccinia virus) is an inner membrane immunodominant protein which is a target of antibodies (Shchelkunov et al., 2002; Meng et al., 2018). In addition to the aforementioned rationale for selecting these antigens, another reason was to make a vaccine by combining these four antigens to assess the immune response. This approach aimed to generate an additional vaccine candidate as having multiple vaccine candidates would greatly enhance the effectiveness of pandemic prevention efforts. After the conservancy analyses of each protein sequence downloaded from the NCBI database, B- and T-cell conserved epitopes were predicted to design a DNA vaccine. DNA vaccines are a type of nucleic acid vaccine with a number of advantages over conventional immunization strategies (Supplementary Figure 4).
DNA vaccines instruct the host cells to produce antigens similar to the virus and the body becomes ready for the future possible infection, which is capable to stop the spreading of the virus and, thus, a manifestation of clinical symptoms (Chavda et al., 2021a; Qin et al., 2021). After the delivery into the host cell, the DNA is transported to the nucleus where transcription takes place. Subsequently, the mRNA is conveyed to the cytoplasm where it undergoes translation to produce the antigens for the vaccine. Figure 8 depicts the schematic representation of the mechanism of action for DNA vaccines upon delivery into the host cell.
Figure 8
Remarkably, although DNA vaccine enters the host nucleus, which is the main concern related to DNA vaccines, the likelihood of its integration into the host genome is extremely low (Williams, 2013). Indeed, there is a number of clinical trials on the application of DNA vaccines for the prevention of various infectious diseases: NCT04591184, NCT01498718, NCT01487876, and NCT04445389. Furthermore, India has designated the emergency use of the first DNA vaccine against COVID-19 in 2021 (Sheridan, 2021) and there are a number of other COVID-19 DNA vaccine candidates in clinical (Sheridan, 2021; Silveira et al., 2021) and preclinical development (Shafaati et al., 2022; Wang et al., 2022). More importantly, DNA vaccines have shown protective efficacy against the mpox virus in non-human primates (Hooper et al., 2004; Hirao et al., 2011). The steps of developing a DNA-based vaccine from the antigen selection to the commercial availability are vaccine design, optimization with proper linkers and addition of immunomodulator sequences, synthesis, transformation into competent bacterial cells for its amplification, extraction, purification, animal immunization, assessment of immune responses, clinical trials, approval, and manufacturing (Papukashvili et al., 2022b; Rcheulishvili et al., 2022). Steps that are needed for developing a potentially universal DNA-based anti-mpox multi-epitope vaccine are given in Figure 9.
Figure 9
Designing a vaccine employing immunoinformatics approaches significantly shortens the time of vaccine development and lays the groundwork for the rational design of the effective vaccine, while in the case of conventional techniques, the development of a vaccine may take decades of laboratory work (Rappuoli and Aderem, 2011). Leveraging immunoinformatics enables us to predict the potential outcome of the vaccine. Upon achieving favorable results with the designed vaccine, it becomes more rational to proceed with in vitro and in vivo screening of the multi-epitope construct.
Leveraging immunoinformatics tools enables us to forecast the potential efficacy of the vaccine. Upon achieving favorable outcomes with the designed vaccine, it becomes more logical to proceed with in vitro and in vivo screening of the multi-epitope construct.
The links of the employed servers and tools are given in Supplementary Table 2. Apparently, there are a number of in silico studies which designed the vaccine and computationally evaluated its efficacy, stability, etc., e.g., vaccines for COVID-19 (Dong et al., 2020; Oliveira et al., 2020), influenza (Behbahani et al., 2021; Sharma et al., 2021; Rcheulishvili et al., 2023), mpox (Hirao et al., 2011; Akhtar et al., 2022), and other viruses (Ali et al., 2019; Mahmudul et al., 2020; Ros-Lucas et al., 2020).
In this study, the vaccine construct designed here was found to have favorable physico-chemical properties and induce strong cellular and humoral immune responses. The immune simulation analysis shows that immunization with the multi-epitope vaccine candidate induces the production of HTLs and CTLs and stimulates the B-cell population, antibodies, and cytokines. The computationally designed vaccine in this study requires to be validated with in vitro and in vivo studies to confirm the outcomes obtained in this study.
5. Conclusion
In summary, the development of new, next-generation, specific vaccine candidates against the mpox virus is undoubtedly essential. The available immunoinformatics approaches allow for the rational design of the potentially protective vaccine and facilitate the process of vaccine development. The results obtained in this study demonstrate that the multi-epitope vaccine designed here may be suggested as an auspicious vaccine candidate which has the prospect of eliciting strong immune responses and providing protection against the mpox virus. In addition, the strategy of developing a universal multi-epitope DNA vaccine used in this study may have a positive impact on the development of a potentially universal vaccine against the mpox virus and other viruses and, thus, will aid in averting future outbreaks or pandemics.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding authors.
Author contributions
NR conceived and designed the study and wrote the original draft. NR and JM contributed to the methodology and investigation. NR, DP, and SF contributed to the visualization. CL and XW revised the manuscript. YH supervised the study. PW supervised and administered the project. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by the Shenzhen Science and Technology Innovation Program (Grant No. KQTD20200909113758004).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2023.1203355/full#supplementary-material
References
1
AhammadI.LiraS. S. (2020). Designing a novel mRNA vaccine against SARS-CoV-2: an immunoinformatics approach. Int. J. Biol. Macromol.162, 820–837. 10.1016/j.ijbiomac.2020.06.213
2
AkhtarN.KaushikV.GrewalR. K.WaniA. K. (2022). Immunoinformatics-Aided Design of a Peptide Based Multiepitope Vaccine Targeting Glycoproteins and Membrane Proteins against Monkeypox Virus. Viruses14, 2374. 10.3390/v14112374
3
AlexanderJ.del GuercioM.-F.MaewalA.QiaoL.FikesJ.ChesnutR. W.et al. (2000). Linear PADRE T helper epitope and carbohydrate B cell epitope conjugates induce specific high titer IgG antibody responses. J. Immunol.164, 1625–1633. 10.4049/jimmunol.164.3.1625
4
AliA.KhanA.KaushikA. C.WangY.AliS. S.JunaidM.et al. (2019). Immunoinformatic and systems biology approaches to predict and validate peptide vaccines against Epstein–Barr virus (EBV). Sci. Rep. 9, 720. 10.1038/s41598-018-37070-z
5
AliM.PandeyR. K.KhatoonN.NarulaA.MishraA.PrajapatiV. K. (2017). Exploring dengue genome to construct a multi-epitope based subunit vaccine by utilizing immunoinformatics approach to battle against dengue infection. Sci. Rep.7, 9232. 10.1038/s41598-017-09199-w
6
BaekM.DimaioF.AnishchenkoI.DauparasJ.OvchinnikovS.LeeG. R.et al. (2021). Accurate prediction of protein structures and interactions using a 3-track neural network. Science373, 871–876. 10.1126/science.abj8754
7
BehbahaniM.MoradiM.MohabatkarH. (2021). In silico design of a multi-epitope peptide construct as a potential vaccine candidate for Influenza A based on neuraminidase protein. Silico Pharmacol.9, 36. 10.1007/s40203-021-00095-w
8
BellamkondaN.LambeU. P.SawantS.NandiS. S.ChakrabortyC.ShuklaD. (2022). Immune response to SARS-CoV-2 vaccines. Biomedicines10, 1464. 10.3390/biomedicines10071464
9
Boyoglu-BarnumS.EllisD.GillespieR. A.HutchinsonG. B.ParkY.MoinS. M.et al. (2021). Quadrivalent influenza nanoparticle vaccines induce broad protection. Nature592, 623–628. 10.1038/s41586-021-03365-x
10
BuchmanG. W.CohenM. E.XiaoY.Richardson-HarmanN.SilveraP.DeTollaL. J.et al. (2010). A protein-based smallpox vaccine protects non-human primates from a lethal monkeypox virus challenge. Vaccine28, 6627–6636. 10.1016/j.vaccine.2010.07.030
11
BuiH. H.SidneyJ.DinhK.SouthwoodS.NewmanM. J.SetteA. (2006). Predicting population coverage of T-cell epitope-based diagnostics and vaccines. BMC Bioinformatics7, 153. 10.1186/1471-2105-7-153
12
BungeE. M.HoetB.ChenL.LienertF.WeidenthalerH.BaerL. R.et al. (2022). The changing epidemiology of human monkeypox—a potential threat? A systematic review. PLoS Negl. Trop. Dis.16, e0010141. 10.1371/journal.pntd.0010141
13
CDC (2022). JYNNEOS Vaccine. Available online at: https://www.cdc.gov/poxvirus/mpox/interim-considerations/jynneos-vaccine.html (accessed June 30, 2023).
14
ChavdaV. P.HossainK.BeladiyaJ.ApostolopoulosV. (2021a). Nucleic acid vaccines for COVID-19: a paradigm shift in the vaccine development arena. Biologics1, 337–356. 10.3390/biologics1030020
15
ChavdaV. P.PandyaR.ApostolopoulosV. (2021b). DNA vaccines for SARS-CoV-2: toward third-generation vaccination era. Expert Rev. Vaccines20, 1549–1560. 10.1080/14760584.2021.1987223
16
ChenX.ZaroJ. L.ShenW. C. (2013). Fusion protein linkers: property, design and functionality. Adv. Drug Deliv. Rev.65, 1357–1369. 10.1016/j.addr.2012.09.039
17
CorbettK. S.FlynnB.FouldsK. E.FrancicaJ. R.Boyoglu-BarnumS.WernerA. P.et al. (2020). Evaluation of the mRNA-1273 vaccine against SARS-CoV-2 in nonhuman primates. N. Engl. J. Med.383, 1544–1555. 10.1056/NEJMoa2024671
18
DongR.ChuZ.YuF.ZhaY. (2020). Contriving multi-epitope subunit of vaccine for COVID-19: immunoinformatics approaches. Front. Immunol.11, 1784. 10.3389/fimmu.2020.01784
19
DoytchinovaI. A.FlowerD. R. (2007a). Identifying candidate subunit vaccines using an alignment-independent method based on principal amino acid properties. Vaccine25, 856–866. 10.1016/j.vaccine.2006.09.032
20
DoytchinovaI. A.FlowerD. R. (2007b). VaxiJen: a server for prediction of protective antigens, tumour antigens and subunit vaccines. BMC Bioinform.8, 4. 10.1186/1471-2105-8-4
21
FlingaiS.CzerwonkoM.GoodmanJ.KudchodkarS. B.MuthumaniK.WeinerD. B. (2013). Synthetic DNA vaccines: improved vaccine potency by electroporation and co-delivered genetic adjuvants. Front. Immunol.4, 354. 10.3389/fimmu.2013.00354
22
GasteigerE.HooglandC.GattikerA.DuvaudS.WilkinsM. R.AppelR. D.et al. (2005). Protein identification and analysis tools on the ExPASy server; the proteomics protocols handbook. Humana Press112, 531–552. 10.1385/1-59259-890-0:571
23
GhafouriF.CohanR. A.SamimiH.HosseiniA.NaderiM.NoorbakhshF.et al. (2022). Development of a multiepitope vaccine against SARS-CoV-2: immunoinformatics study. JMIR Bioinforma. Biotechnol.3, e36100. 10.2196/36100
24
HaasE. J.AnguloF. J.McLaughlinJ. M.AnisE.SingerS. R.KhanF.et al. (2021). Impact and effectiveness of mRNA BNT162b2 vaccine against SARS-CoV-2 infections and COVID-19 cases, hospitalisations, and deaths following a nationwide vaccination campaign in Israel: an observational study using national surveillance data. Lancet397, 1819–1829. 10.1016/S0140-6736(21)00947-8
25
HeoL.ParkH.SeokC. (2013). GalaxyRefine: protein structure refinement driven by side-chain repacking. Nucleic Acids Res.41, 384–388. 10.1093/nar/gkt458
26
HeraudJ.Edghill-smithY.AyalaV.KaliszI.ParrinoJ.KalyanaramanV. S.et al. (2022). Subunit recombinant vaccine protects against monkeypox. J. Immunol.177, 2552–2564. 10.4049/jimmunol.177.4.2552
27
HeskinJ.BelfieldA.MilneC.BrownN.WaltersY.ScottC.et al. (2022). Transmission of monkeypox virus through sexual contact – a novel route of infection. J. Infect.85, 334–363. 10.1016/j.jinf.2022.05.028
28
HiraoL. A.Draghia-akliR.PriggeJ. T.YangM.SatishchandranA.WuL.et al. (2011). Multivalent smallpox DNA vaccine delivered by intradermal electroporation drives protective immunity in nonhuman primates against lethal monkeypox challenge. J. Infect. Dis.203, 95–102. 10.1093/infdis/jiq017
29
HobernikD.BrosM. (2018). DNA vaccines—how far from clinical use?Int. J. Mol. Sci.19, 3605. 10.3390/ijms19113605
30
HooperJ. W.ThompsonE.WilhelmsenC.ZimmermanM.IchouM. A.SteffenS. E.et al. (2004). Smallpox DNA vaccine protects nonhuman primates against lethal monkeypox. J. Virol.78, 4433–4443. 10.1128/JVI.78.9.4433-4443.2004
31
HooverD. M.WuZ.TuckerK.LuW.LubkowskiJ. (2003). Antimicrobial characterization of human β-defensin 3 derivatives. Antimicrob. Agents Chemother.47, 2804–2809. 10.1128/AAC.47.9.2804-2809.2003
32
JavierF.SánchezM.Martínez-sellésM.Ruiz-galianaJ.RamosP. D. L. (2022). Monkeypox in humans: a new outbreak. Off. J. Spanish Soc. Chemother.6, 2022.
33
JespersenM. C.PetersB.NielsenM.MarcatiliP. (2017). BepiPred-2.0: improving sequence-based B-cell epitope prediction using conformational epitopes. Nucleic Acids Res.45, W24–W29. 10.1093/nar/gkx346
34
JezekZ.SzczeniowskiM.PalukuK. M. (1987). Human monkeypox: clinical features of 282 patients. J. Infect. Dis.156, 293–298. 10.1093/infdis/156.2.293
35
JumperJ.EvansR.PritzelA.GreenT.FigurnovM.RonnebergerO.et al. (2021). Highly accurate protein structure prediction with AlphaFold. Nature596, 583–589. 10.1038/s41586-021-03819-2
36
KandeelM.MorsyM. A.El-lateefH. M. A.MarzokM.El-BeltagiH. S.KhodairK. M.et al. (2023). Efficacy of the modified vaccinia Ankara virus vaccine and the replication-competent vaccine ACAM2000 in monkeypox prevention Mahmoud. Int. Immunopharmacol.119, 110206. 10.1016/j.intimp.2023.110206
37
KhanM. T.IslamR.JerinT. J.MahmudA.KhatunS.KobirA.et al. (2021). Immunoinformatics and molecular dynamics approaches: next generation vaccine design against West Nile virus. PLoS ONE16, e0253393. 10.1371/journal.pone.0253393
38
KmiecD.KirchhoffF. (2022). Monkeypox: a new threat?Int. J. Mol. Sci.23, 7866. 10.3390/ijms23147866
39
KowalzikF.SchreinerD.JensenC.TeschnerD.GehringS.ZeppF. (2021). mRNA-based vaccines. Vaccines9, 390. 10.3390/vaccines9040390
40
KozakovD.HallD. R.XiaB.PorterK. A.PadhornyD.YuehC.et al. (2017). The ClusPro web server for protein-protein docking. Nat. Protoc.12, 255–278. 10.1038/nprot.2016.169
41
KozlovM. (2022). Monkeypox outbreaks: 4 key questions researchers have. Nature606, 238–239. 10.1038/d41586-022-01493-6
42
KugelmanJ. R.JohnstonS. C.MulembakaniP. M.KisaluN.LeeM. S.KorolevaG.et al. (2014). Genomic variability of monkeypox virus among humans, democratic Republic of the Congo. Emerg. Infect. Dis.20, 232–239. 10.3201/eid2002.130118
43
LiuC.PapukashviliD.DongY.WangX.HuX.YangN.et al. (2022a). Identification of tumor antigens and design of mRNA vaccine for colorectal cancer based on the immune subtype. Front. Cell Dev. Biol.9, 783527. 10.3389/fcell.2021.783527
44
LiuC.RcheulishviliN.ShenZ.PapukashviliD.XieF.WangZ.et al. (2022b). Development of an LNP-Encapsulated mRNA-RBD vaccine against SARS-CoV-2 and its variants. Pharmaceutics14, 1101. 10.3390/pharmaceutics14051101
45
LockerJ. K.KuehnA.SchleichS.RutterG.HohenbergH.WepfR.et al. (2000). Entry of the two infectious forms of vaccinia virus at the plasma membane is signaling-dependent for the IMV but not the EEV. Mol. Biol. Cell11, 2497–2511. 10.1091/mbc.11.7.2497
46
LumF.-M.RuestaA. T.TayM. Z.LinR. T. P.LyeD. C.RéniaL.et al. (2022). Monkeypox: disease epidemiology, host immunity and clinical interventions. Nat. Rev. Immunol.22, 597–613. 10.1038/s41577-022-00775-4
47
MaaJ. S.EstebanM. (1987). Structural and functional studies of a 39,000-Mr immunodominant protein of vaccinia virus. J. Virol.61, 3910–3919. 10.1128/jvi.61.12.3910-3919.1987
48
MahmudulH.ShifulI.SouravC.AbuH. M.KaziF. A.ZiaulF. J.et al. (2020). Contriving a chimeric polyvalent vaccine to prevent infections caused by herpes simplex virus (Type-1 and Type-2): an exploratory immunoinformatic approach. J. Biomol. Struct. Dyn.38, 2898–2915. 10.1080/07391102.2019.1647286
49
MartinelliD. D. (2022). In silico vaccine design: a tutorial in immunoinformatics. Healthc. Anal.2, 100044. 10.1016/j.health.2022.100044
50
MartínezJ. I.MontalbánE. G.BuenoS. J.MartínezF. M.Juli,áA. N.DíazS.et al. (2022). Monkeypox outbreak predominantly affecting men who have sex with men, Madrid, Spain, 26 April to 16 June 2022. Eurosurveillance27, 2200471. 10.2807/1560-7917.ES.2022.27.27.2200471
51
MaurerD. M.HarringtonB.ArmyD.HospitalC.HoodF.LaneJ. M. (2003). Smallpox vaccine: contraindications, administration, and adverse reactions. Am. Fam. Physician68, 889–896.
52
Maurer-StrohS.KrutzN. L.KernP. S.GunalanV.NguyenM. N.LimviphuvadhV.et al. (2019). AllerCatPro-prediction of protein allergenicity potential from the protein sequence. Bioinformatics35, 3020–3027. 10.1093/bioinformatics/btz029
53
MekonnenD.MengistH. M.JinT. (2022). SARS-CoV-2 subunit vaccine adjuvants and their signaling pathways. Expert Rev. Vaccines21, 69–81. 10.1080/14760584.2021.1991794
54
MengX.KaeverT.YanB.TraktmanP.ZajoncD. M.PetersB.et al. (2018). Characterization of murine antibody responses to vaccinia virus envelope protein A14 reveals an immunodominant antigen lacking of effective neutralization targets. Virology518, 284–292. 10.1016/j.virol.2018.03.005
55
MeyerH.PerrichotM.StemmlerM.EmmerichP.SchmitzH.VaraineF.et al. (2002). Outbreaks of disease suspected of being due to human monkeypox virus infection in the Democratic Republic of Congo in 2001. J. Clin. Microbiol.40, 2919–2921. 10.1128/JCM.40.8.2919-2921.2002
56
MiletoD.RivaA.CutreraM.MoscheseD.ManconA.MeroniL.et al. (2022). New challenges in human monkeypox outside Africa: a review and case report from Italy. Travel Med. Infect. Dis.49, 102386. 10.1016/j.tmaid.2022.102386
57
MirditaM.SchützeK.MoriwakiY.HeoL. (2022). ColabFold: making protein folding accessible to all. Nat. Methods19, 679–682. 10.1038/s41592-022-01488-1
58
MoscheseD.GiacomelliA.BeltramiM.PozzaG.MiletoD.ReatoS.et al. (2022). Hospitalisation for monkeypox in Milan, Italy. Travel Med. Infect. Dis.49, 102417. 10.1016/j.tmaid.2022.102417
59
NelsonD. L.CoxM. M.HoskinsA. A. (2021). Lehninger Principles of Biochemistry. 8th ed. Austin: Macmillan Learning.
60
NezafatN.GhasemiY.JavadiG.KhoshnoudM. J.OmidiniaE. (2014). A novel multi-epitope peptide vaccine against cancer: an in silico approach. J. Theor. Biol.349, 121–134. 10.1016/j.jtbi.2014.01.018
61
NguyenM. N.KrutzN. L.LimviphuvadhV.LopataA. L.GerberickG. F.Maurer-StrohS. (2022). AllerCatPro 2.0: a web server for predicting protein allergenicity potential. Nucleic Acids Res.50, W36–W43. 10.1093/nar/gkac446
62
OliveiraS. C.MagalhãesM. T. Q.De HomanE. J. (2020). Immunoinformatic analysis of SARS-CoV-2 nucleocapsid protein and identification of COVID-19 vaccine targets. Front. Immunol.11, 587615. 10.3389/fimmu.2020.587615
63
OrvizE.NegredoA.AyerdiO.VázquezA.Muñoz-GomezA.MonzónS.et al. (2022). Monkeypox outbreak in Madrid (Spain): clinical and virological aspects. J. Infect.85, 412–417. 10.1016/j.jinf.2022.07.005
64
PapukashviliD.RcheulishviliN.LiuC.JiY.HeY.WangP. G. (2022a). Self-amplifying RNA approach for protein replacement therapy. Int. J. Mol. Sci.23, 12884. 10.3390/ijms232112884
65
PapukashviliD.RcheulishviliN.LiuC.WangX.HeY.WangP. G. (2022b). Strategy of developing nucleic acid-based universal monkeypox vaccine candidates. Front. Immunol.13, 1050309. 10.3389/fimmu.2022.1050309
66
PardiN.HoganM. J.PorterF. W.WeissmanD. (2018). mRNA vaccines-a new era in vaccinology. Nat. Rev. Drug Discov.17, 261–279. 10.1038/nrd.2017.243
67
PonomarenkoJ.BuiH. H.LiW.FussederN.BourneP. E.SetteA.et al. (2008). ElliPro: a new structure-based tool for the prediction of antibody epitopes. BMC Bioinform.9, 514. 10.1186/1471-2105-9-514
68
QamarM. T. U.ShokatZ.MuneerI.AshfaqU. A.JavedH.AnwarF.et al. (2020). Multiepitope-based subunit vaccine design and evaluation against respiratory syncytial virus using reverse vaccinology approach. Vaccines8, 288. 10.3390/vaccines8020288
69
QinF.XiaF.ChenH.CuiB.FengY.ZhangP.et al. (2021). A guide to nucleic acid vaccines in the prevention and treatment of infectious diseases and cancers: from basic principles to current applications. Front. Cell Dev. Biol.9, 633776. 10.3389/fcell.2021.633776
70
RahmaniA.BaeeM.SalekiK.MoradiS.NouriH. R. (2021). Applying high throughput and comprehensive immunoinformatics approaches to design a trivalent subunit vaccine for induction of immune response against emerging human coronaviruses SARS-CoV, MERS-CoV and SARS-CoV-2. J. Biomol. Struct. Dyn.40, 6097–6113. 10.21203/rs.3.rs-92515/v1
71
RapinN.LundO.BernaschiM.CastiglioneF. (2010). Computational immunology meets bioinformatics: the use of prediction tools for molecular binding in the simulation of the immune system. PLoS ONE5, e9862. 10.1371/journal.pone.0009862
72
RapinN.LundO.CastiglioneF. (2011). Immune system simulation online. Bioinformatics27, 2013–2014. 10.1093/bioinformatics/btr335
73
RappuoliR.AderemA. (2011). A 2020 vision for vaccines against HIV, tuberculosis and malaria. Nature473, 463–469. 10.1038/nature10124
74
RcheulishviliN.MaoJ.PapukashviliD.LiuC.WangZ.ZhaoJ.et al. (2023). Designing multi-epitope mRNA construct as a universal influenza vaccine candidate for future epidemic/pandemic preparedness. Int. J. Biol. Macromol.226, 885–899. 10.1016/j.ijbiomac.2022.12.066
75
RcheulishviliN.PapukashviliD.LiuC.JiY.HeY.WangP. G. (2022). Promising strategy for developing mRNA-based universal influenza virus vaccine for human population, poultry, and pigs– focus on the bigger picture. Front. Immunol.13, 1025884. 10.3389/fimmu.2022.1025884
76
ReynissonB.AlvarezB.PaulS.PetersB.NielsenM. (2020). NetMHCpan-4.1 and NetMHCIIpan-4.0: improved predictions of MHC antigen presentation by concurrent motif deconvolution and integration of MS MHC eluted ligand data. Nucleic Acids Res.48, W449–W454. 10.1093/nar/gkaa379
77
Ros-LucasA.Correa-FizF.Bosch-CamósL.RodriguezF.Alonso-PadillaJ. (2020). Computational analysis of african swine fever virus protein space for the design of an epitope-based vaccine ensemble. Pathogens. 9, 1078.
78
RybakovaY.KowalskiP. S.HuangY.GonzalezJ. T.HeartleinM. W.DerosaF.et al. (2019). mRNA delivery for therapeutic Anti-HER2 antibody expression in vivo. Mol. Ther.27, 1415–1423. 10.1016/j.ymthe.2019.05.012
79
SahR.AbdelaalA.RedaA.KatameshB. E.ManirambonaE. (2022). Monkeypox and its possible sexual transmission: where are we now with its evidence?Pathogens11, 924. 10.3390/pathogens11080924
80
SeangS.BurrelS.TodescoE.LeducqV.MonselG.PluartD.et al. (2022). Evidence of human-to-dog transmission of monkeypox virus. Lancet400, 658–659. 10.1016/S0140-6736(22)01487-8
81
ShafaatiM.MirzaeiR.SaidijamM.AmirheidariB.SoleimaniM.MahakiH.et al. (2022). A brief review on DNA vaccines in the era of COVID-19. Future Virol.17, 49–66. 10.2217/fvl-2021-0170
82
ShafaatiM.ZandiM. (2022). State-of-the-art on monkeypox virus: an emerging zoonotic disease. Infection50, 1425–1430. 10.1007/s15010-022-01935-3
83
ShafaatiM.ZandiM. (2023). Human monkeypox (hMPXV) re-emergence: host immunity status and current vaccines landscape. J. Med. Virol.95, e28251. 10.1002/jmv.28251
84
SharmaS.KumariV.VasantB.MukherjeeA.PandeyR.KondabagilK. (2021). Immunoinformatics approach for a novel multi-epitope subunit vaccine design against various subtypes of Influenza A virus. Immunobiology226, 152053. 10.1016/j.imbio.2021.152053
85
ShchelkunovS. N.TotmeninA. V.SafronovP. F.MikheevM. V.GutorovV. V.RyazankinaO. I.et al. (2002). Analysis of the monkeypox virus genome. Virology197, 172–194. 10.1006/viro.2002.1446
86
SheridanC. (2021). First COVID-19 DNA vaccine approved, others in hot pursuit. Nat. Biotechnol.39, 1479–1482. 10.1038/d41587-021-00023-5
87
SilveiraM. M.MoreiraG. M. S. G.MendonçaM. (2021). DNA vaccines against COVID-19: Perspectives and challenges. Life Sci.267, 118919. 10.1016/j.lfs.2020.118919
88
SklenovskáN.Van RanstM. (2018). Emergence of monkeypox as the most important orthopoxvirus infection in humans. Front. Public Heal.6, 241. 10.3389/fpubh.2018.00241
89
TarrahimofradH.RahimnahalS.ZamaniJ.JahangirianE.AminzadehS. (2021). Designing a multi-epitope vaccine to provoke the robust immune response against influenza A H7N9. Sci. Rep.11, 24485. 10.1038/s41598-021-03932-2
90
VandenbogaertM.KwasiborskiA.GonofioE.DescorpsS.SelekonB.AndyA.et al. (2022). Nanopore sequencing of a monkeypox virus strain isolated from a pustular lesion in the Central African Republic. Sci. Rep.12, 10768. 10.1038/s41598-022-15073-1
91
VavilisT.StamoulaE.AinatzoglouA.SachinidisA.LamprinouM.DardalasI.et al. (2023). mRNA in the context of protein replacement therapy. Pharmaceutics15, 166. 10.3390/pharmaceutics15010166
92
WalshE. E.FrenckR. W.FalseyA. R.KitchinN.AbsalonJ.GurtmanA.et al. (2020). Safety and immunogenicity of Two RNA-based Covid-19 vaccine candidates. N. Engl. J. Med.383, 2439–2450. 10.1056/NEJMoa2027906
93
WangC.ZouQ. (2023). Prediction of protein solubility based on sequence physicochemical patterns and distributed representation information with DeepSoluE. BMC Biol.21, 12. 10.1186/s12915-023-01510-8
94
WangX.LiuC.RcheulishviliN.PapukashviliD.XieF.ZhaoJ.et al. (2023). Strong immune responses and protection of PcrV and OprF-I mRNA vaccine candidates against Pseudomonas aeruginosa. npj Vaccines8, 76. 10.1038/s41541-023-00672-4
95
WangX.RcheulishviliN.CaiJ.LiuC.XieF.HuX.et al. (2022). Development of DNA vaccine candidate against SARS-CoV-2. Viruses14, 1049. 10.3390/v14051049
96
WaterhouseA. M.ProcterJ. B.MartinD. M. A.ClampM.BartonG. J. (2009). Jalview Version 2-a multiple sequence alignment editor and analysis workbench. Bioinformatics25, 1189–1191. 10.1093/bioinformatics/btp033
97
WilliamsJ. (2013). Vector design for improved DNA vaccine efficacy, safety and production. Vaccines1, 225–249. 10.3390/vaccines1030225
98
World Health Organization (2022a). Monkeypox. Available online at: https://www.who.int/news-room/fact-sheets/detail/monkeypox (accessed September 16, 2022).
99
World Health Organization (2022b). WHO Director-General Declares the Ongoing Monkeypox Outbreak a Public Health Emergency of International Concern. Geneva: WHO.
100
XiaD.JinR.ByagathvalliG.YuH.YeL.LuC. Y.et al. (2021). An ultra-low-cost electroporator with microneedle electrodes (ePatch) for SARS-CoV-2 vaccination. Proc. Natl. Acad. Sci. U. S. A.118, e2110817118. 10.1073/pnas.2110817118
101
YangZ. (2022). Monkeypox: a potential global threat?J. Med. Virol.94, 4034–4036. 10.1002/jmv.27884
102
YanoA.OnozukaA.Asahi-OzakiY.ImaiS.HanadaN.MiwaY.et al. (2005). An ingenious design for peptide vaccines. Vaccine23, 2322–2326. 10.1016/j.vaccine.2005.01.031
103
ZaeckL. M.LamersM. M.VerstrepenB. E.BestebroerT. M.RoyenM. E.van GötzH.et al. (2022). Low levels of monkeypox virus neutralizing antibodies after MVA-BN vaccination in healthy individuals. Nat. Med.29, 270–278. 10.1038/s41591-022-02090-w
104
ZandiM.ShafaatiM.HosseiniF. (2023). Mechanisms of immune evasion of monkeypox virus. Front. Microbiol.14, 1106247. 10.3389/fmicb.2023.1106247
105
ZhangJ.LiangY.ZhangY. (2011). Atomic-level protein structure refinement using fragment guided molecular dynamics conformation sampling. Structure19, 1784–1795. 10.1016/j.str.2011.09.022
106
ZhangR. R.WangZ. J.ZhuY. L.TangW.ZhouC.ZhaoS. Q.et al. (2023). Rational development of multicomponent mRNA vaccine candidates against mpox. Emerg. Microbes Infect.12, 2192815. 10.1080/22221751.2023.2192815
107
ZhuM.JiJ.ShiD.LuX.WangB.WuN.et al. (2022). Unusual global outbreak of monkeypox: what should we do?Front. Med.16, 507–517. 10.1007/s11684-022-0952-z
Summary
Keywords
monkeypox, mpox virus, universal vaccine, DNA vaccine, immunoinformatics, in silico, epidemics, pandemics
Citation
Rcheulishvili N, Mao J, Papukashvili D, Feng S, Liu C, Wang X, He Y and Wang PG (2023) Design, evaluation, and immune simulation of potentially universal multi-epitope mpox vaccine candidate: focus on DNA vaccine. Front. Microbiol. 14:1203355. doi: 10.3389/fmicb.2023.1203355
Received
10 April 2023
Accepted
03 July 2023
Published
21 July 2023
Volume
14 - 2023
Edited by
Sinosh Skariyachan, St. Pius X College, India
Reviewed by
Anand Anbarasu, VIT University, India; Milad Zandi, Tehran University of Medical Sciences, Iran
Updates
Copyright
© 2023 Rcheulishvili, Mao, Papukashvili, Feng, Liu, Wang, He and Wang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yunjiao He heyj@sustech.edu.cnPeng George Wang wangp6@sustech.edu.cn
†These authors have contributed equally to this work
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.