Next Article in Journal
Does Avian Coronavirus Co-Circulate with Avian Paramyxovirus and Avian Influenza Virus in Wild Ducks in Siberia?
Next Article in Special Issue
Prediction of Antigenic Distance in Influenza A Using Attribute Network Embedding
Previous Article in Journal
Special Issue “HPV in the Head and Neck Region 2.0”
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 15
Issue 5
10.3390/v15051120
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Development of a Multi-Epitope Universal mRNA Vaccine Candidate for Monkeypox, Smallpox, and Vaccinia Viruses: Design and In Silico Analyses

by
Nino Rcheulishvili
1,†,
Jiawei Mao
1,†,
Dimitri Papukashvili
1,
Shunping Feng
1,
Cong Liu
1,
Xidan Yang
1,2,
Jihui Lin
1,2,
Yunjiao He
1,* and
Peng George Wang
1,*
1
Department of Pharmacology, School of Medicine, Southern University of Science and Technology, Shenzhen 518000, China
2
School of Nursing, Southwest Medical University, Luzhou 646000, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Viruses 2023, 15(5), 1120; https://doi.org/10.3390/v15051120
Submission received: 14 April 2023 / Revised: 3 May 2023 / Accepted: 4 May 2023 / Published: 7 May 2023
(This article belongs to the Special Issue Efficacy, Safety, and Immunogenicity of Vaccines against Viruses: From Network Medicine to Clinical Experimentation)
Browse Figures
Versions Notes

Abstract

:
Notwithstanding the presence of a smallpox vaccine that is effective against monkeypox (mpox), developing a universal vaccine candidate against monkeypox virus (MPXV) is highly required as the mpox multi-country outbreak has increased global concern. MPXV, along with variola virus (VARV) and vaccinia virus (VACV), belongs to the Orthopoxvirus genus. Due to the genetic similarity of antigens in this study, we have designed a potentially universal mRNA vaccine based on conserved epitopes that are specific to these three viruses. In order to design a potentially universal mRNA vaccine, antigens A29, A30, A35, B6, and M1 were selected. The conserved sequences among the three viral species—MPXV, VACV, and VARV—were detected, and B and T cell epitopes containing the conserved elements were used for the design of the multi-epitope mRNA construct. Immunoinformatics analyses demonstrated the stability of the vaccine construct and optimal binding to MHC molecules. Humoral and cellular immune responses were induced by immune simulation analyses. Eventually, based on in silico analysis, the universal mRNA multi-epitope vaccine candidate designed in this study may have a potential protection against MPXV, VARV, and VACV that will contribute to the advancement of prevention strategies for unpredictable pandemics.

1. Introduction

Unforeseeable outbreaks of infectious diseases are causing a rise in the worldwide risk to public health. Monkeypox virus (MPXV) has emerged in May 2022 and affected more than 86,000 people until now [1]. Monkeypox (mpox) is a zoonotic disease that has been largely neglected although there were cases in its endemic areas—West and Central Africa. MPXV, along with variola virus (VARV), vaccinia virus (VACV), and cowpox virus, belongs to the genus Orthopoxvirus, family Poxviridae. MPXV is a large, 200–300 nm, brick-shaped virus that exists in two different infectious forms: extracellular enveloped virions (EVs) and intracellular mature virions (MVs). EVs have an extra envelope compared to MVs. The genome of MPXV is about 197 kb linear double-stranded DNA (dsDNA) that encodes approximately 190 proteins [2,3,4]. Due to their intricate structure, many viral antigens and their functions still need to be studied.
Except for Africa, there was an outbreak in the United States (US) in 2003 when imported prairie dogs from Ghana spread the virus, and, consequently, there were 47 confirmed cases [5,6]. Interestingly, the prairie dogs were housed with African rodents, so the virus could have been transmitted via the rodents [7]. In addition, single cases of mpox were identified in different countries such as the US [8,9], Israel, the United Kingdom (UK), and Singapore and all were linked to travel to Nigeria [10]. What caused the MPXV virus outbreak this time? This and many other questions are unanswered until now. However, it is assumed that eradicating smallpox and ceasing vaccination in 1980 globally led to increased cases of mpox. This theory is quite convincing as more than 70% of people are unvaccinated against the smallpox virus today, and most cases of mpox take place in younger people. Interestingly, most of the cases occur in men who have sex with men (MSM), and close physical contact plays a key role in transmission [11,12,13]. Indeed, the sex-related infection rate was always observed even before the current outbreak—males predominated over females, while children accounted for the majority of cases [7]. The fatality rates range from 1% to 11% [14], and the disease is more severe in young children [15,16]. There are two genetic clades of MPXV: the Central African (Democratic Republic of the Congo (DRC), previously known as Zaire) and the West African clades. Out of these two clades, the Central African clade is more virulent and deadly. Unlike the coronavirus disease 2019 (COVID-19) and influenza viruses, the MPXV is a dsDNA [17] virus and is more stable. Generally, the mutation rate of DNA viruses is much lower than RNA viruses. Therefore, the assumption of an MPXV mutation as the only reason for the current outbreak would not be rational. However, a few mutations still take place, and it should not be completely neglected. Even though a smallpox vaccine is quite effective against the MPXV virus, there is still a necessity for developing a new, universal vaccine based on conserved elements. Moreover, apparently, the smallpox vaccine cannot completely protect from MPXV [18]. Additionally, in the worst-case scenario, if certain mutations take place in MPXV and the effectiveness of the smallpox vaccine to MPXV decreases, the situation will also worsen. Although the daily confirmed cases of mpox are significantly decreasing, it is globally spread, and the re-emergence of the same clade or even the spread of the Central African clade is anticipated. Remarkably, messenger RNA (mRNA)-based vaccines have revolutionized the field of vaccinology due to their favorable safety profile, low-cost manufacturing, high potency, and rapid development among other impressive advantages [19,20,21]. mRNA vaccine contains the antigen sequence that is translated into the corresponding protein after the introduction into the host body. After the antigen is released, it is recognized by antigen-presenting cells (APCs) and phagocytosed. The vaccine antigen is then processed into small peptides that are presented on the cell surface via major histocompatibility complex (MHC) I and II, and the cellular and humoral immune responses are induced [20].
In this research, five antigen proteins were selected from MPXV, VACV, and VARV. The corresponding sequences were retrieved and aligned. Conserved sequences were detected in each protein. Relevant T cell epitopes that were tested through experimentation were searched in the immune epitope database (IEDB), sorted according to the optimal results, and selected, while B cell epitopes were predicted via IEDB due to the lack of B cell epitope data in the database. The multi-epitope mRNA vaccine made up of the conserved epitopes of A29, A30, A35, B6, and M1 proteins (encoded by the genes A29L, A30L, A35R, B6R, and M1R, respectively) was constructed. Various properties of the vaccine were assessed in silico. Based on the obtained results, the mRNA construct proposed in this study has potentially high efficiency and elicits protection from MPXV, VARV, and VACV. The proposed construct represents a favorable candidate for developing the mRNA vaccine for global purposes. The contribution of the design of this potentially universal vaccine candidate is one step forward in the advancement of vaccine development and pandemic alertness.

2. Methods

2.1. Antigen Selection

Five antigens that are common for MPXV, VARV, and VACV were selected: A29 (A30 in VARV, A27 in VACV), A30 (A31 in VARV, A28 in VACV), A35 (A36 in VARV, A33 in VACV), M1 (M1 in VARV, L1 in VACV), and B6 (B7 in VARV, B5 in VACV). These proteins are highly conserved among the orthopoxviruses and serve important functions in immune response [3,22,23,24,25,26,27]. All five antigens were selected based on their functions in immune response and their successful application in previous vaccine studies [24,27,28,29]. The selected antigens along with their functions are given in Table 1.

2.2. Selection of Conserved Regions and Epitopes

The amino acid sequences of A29, A30, A35, B6, and M1 of MPXV and the corresponding proteins of VARV and VACV were retrieved from the National Center for Biotechnology Information (NCBI) database in FASTA format. In total, 372 sequences were downloaded for the antigen A29: 159 A29 (MPXV), 120 A30 (VARV), and 93 A27 (VACV). This protein is 110 amino acids in length, thus, the sequences longer or shorter than 110 amino acids were removed, and the rest (176 sequences) were aligned together. A total of 433 sequences were retrieved for the antigen A30: 363 A30 (MPXV), 8 A31 (VARV), and 62 A28 (VACV). As the full length of this protein is 146 amino acids, all the shorter or longer sequences were removed, and the remaining 158 sequences were aligned together. A total of 365 sequences of the antigen A35 were downloaded: 153 A35 (MPXV), 109 A36 (VARV), and 103 A33 (VACV). As the full-length sequence of A35 is 181 amino acids, A36 is 184, and A33 is 185, all the records shorter than 181 and longer than 185 were removed, and the remaining 194 sequences were aligned. In total 365 sequences of antigen B6 were retrieved: 180 B6 (MPXV), 6 B7 (VARV), and 179 B5 (VACV). The full length of this protein is 317 amino acids, thus, all the sequences shorter or longer than 317 amino acids were removed, and the rest of the sequences (196) were aligned. A total of 391 sequence records of antigen M1 were retrieved from NCBI: 82 M1 (MPXV), 6 M1 (VARV), and 303 L1 (VACV), which were aligned together. As the full length of this protein is 250 amino acids, all the sequences shorter or longer than 250 were removed, and the remaining 346 sequences were aligned together. The conserved sequences were identified via the bioinformatics software Jalview 2.11.1.4 [30]. The filter of conservation threshold was given 10 “below threshold”.
Experimentally tested T cell epitopes and MHC ligands of each antigen were found on IEDB [31] according to the species of orthopoxviruses (MPXV, VARV, and VACV) and host (human). IEDB is an excellent database supported by experimental data of humans and animals on antibody and T cell epitopes. IEDB also incorporates immunoinformatics tools for epitope prediction [31]. As there were no B cell epitopes found on IEDB, the Bepipred linear epitope prediction tool (v2.0) was used to predict B cell epitopes [32]. Ultimately, experimentally tested T cell epitopes were retrieved from IEDB; epitopes containing conserved sequences were ranked and the most studied ones were selected. On the other hand, linear B cell epitopes were predicted using the consensus sequence of each antigen.

2.3. Vaccine Design

For the construction of the multi-epitope mRNA vaccine, the epitopes were arranged orderly. B cell epitopes were joined with a flexible linker KK [33,34], while a GGGS linker was used to connect the T cell epitopes [35]. The leading sequence of tissue plasminogen activator (tPA) that augments the antigen presentation [36] was used as a signal sequence (MDAMKRGLCCVLLLCGAVFVSPS). The incorporation of the amino acid sequence of interleukin 6 (IL-6) was done to enhance the immunogenicity of the vaccine [37]. The pan-HLA DR binding epitope (PADRE sequence) was connected to the IL-6 sequence via a GGGS linker and to an epitope of the B cell via EAAAK [34,38]. A polyhistidine (6x) tag was placed on the C-terminal of the vaccine sequence for facilitating the fusion protein detection [39]. The complete open reading frame of the proposed vaccine is given in the Supplementary Materials.

2.4. Prediction of Vaccine Properties

To check whether the antigen could provoke an allergic response, the AlgPred server, a highly accurate tool, was used to predict its allergenicity [40]. As the prediction tool, a hybrid approach—Support Vector Machine (SVMc) algorithm + IgE epitope + ARPs BLAST + MAST—was selected [40,41]. This method detects the allergenicity of the protein based on the composition of amino acids and dipeptides using SVM, which is a motif-based technique using the software MAST. Ultimately, the tool specifies the antigen as an allergen if it contains a segment identical to the known IgE epitopes or similar to allergen-representative proteins [40,41]. For the prediction of antigenicity, server Vaxijen v2.0, which is based on an alignment-independent prediction of protective antigens, was used [42]. Vaxijen categorizes viral, bacterial, and tumor antigens according to their physicochemical characteristics [43].
ProtParam was used to calculate the physicochemical properties [44]. The following parameters were characterized: the molecular weight (MW) of the multi-peptide, the atomic and amino acid composition, the theoretical isoelectric point (pI), the instability, the estimated half-life, the grand average of hydropathicity (GRAVY), and the aliphatic indexes. MW was calculated by adding the average isotopic masses of amino acids and the average isotopic mass of one water molecule. In addition to the MW, the pI, which plays a crucial role in pH-dependent properties, was computed based on the pKa value of amino acids, which plays an essential role in characterizing the pH-dependency of the protein. The estimated half-life denotes the duration required from the synthesis of the protein until its decay and reduction to half of its original amount within the cell. The instability index is used to evaluate the stability of the protein in a test tube. An instability index of less than 40 is considered an indication of stability. When this measure is >40, the protein of interest is considered unstable. To compute GRAVY, the hydropathy values of each amino acid are summed up and then divided by the total residue number in the protein. A higher number denotes that the amino acids are more hydrophobic [44]. An aliphatic index is used to characterize the protein relative volume that is taken up by amino acids that have aliphatic side chains, and it is considered a positive factor in augmenting its thermostability.

2.5. Tertiary Structure Prediction and Evaluation of Its Quality; Discontinuous B Cell Epitope Prediction

RoseTTAFold was used to generate the tertiary structure of the vaccine protein [45]. By combining the sequence data of a one-dimensional protein with the two-dimensional data of the distances between the amino acids, as well as the prediction of the three-dimensional atomic structure, this tool is able to forecast the configuration and interplay of the protein [45]. After preparing the structure of the vaccine protein, it was refined with GalaxyRefine [46]. Through the use of molecular dynamics simulation, this immunoinformatics approach reconstructs and reorganizes the amino acid side chains, while also relaxing the overall protein structure [46]. Subsequently, the Ramachandran plot and ERRAT were used to verify the quality of the protein tertiary structure. The Ramachandran plot visualizes the energetically allowed regions for backbone dihedral angles ψ against φ of amino acid residues in protein structures, which allows the testing of the quality of protein structures [47]. The ERRAT score serves as a quality indicator for the non-bonded interactions, where a larger score suggests a superior quality of the protein tertiary structure. ElliPro was employed for the prediction of B cell conformational epitopes based on the 3D structure of the protein [48].

2.6. Immune Simulation

Humoral and cellular immune responses driven against the proposed multi-epitope vaccine protein were analyzed via the C-ImmSim server (https://kraken.iac.rm.cnr.it/C-IMMSIM/; accessed on 31 August 2022). This immunoinformatics approach provides the service of immune simulation and characterizes the immune system responses, both humoral and cellular, towards vaccination. For predicting immune epitopes and analyzing immune interactions, C-ImmSim applies a position-specific scoring matrix in combination with machine learning techniques. After uploading the sequence of vaccine antigen in FASTA format, the server predicted the immune responses [49]. The “Allele Frequency Net Database” identified the most widespread HLA-A, HLA-B, and DRB alleles globally. The outcome indicated HLA-A*02:01, HLA-A*01:01, HLA-B*07:02, HLA-B*08:01, DRB1*07:01, and DRB1*15:01. The vaccine was administered thrice with an interval of four weeks between each injection. The simulation volume was adjusted at 10, and the simulation progressed through 270 steps. The vaccine used during the simulation did not comprise LPS. The random seed was 12345, and time periods were set at 1, 85, and 169 [50].

2.7. In Silico Validation of Vaccine Protein Binding to the Host Receptors

In order to assess the binding capacity of vaccine protein to antigen recognition receptors, molecular docking was carried out with a ClusPro server. ClusPro calculates the docking interaction between two protein structures and provides a list of potential complexes in order of priority. It is performed based on a predicted conformation of ligand, orientation, and position as well as binding affinity analysis. Ultimately, the complexes that exhibit good electrostatic and desolvation energies are chosen [51]. A molecular docking of vaccine construct with MHC-I (HLA-A*02:01) (6TDS) and MHC-II (HLA-DRB1*01:01) (1AQD) host receptors was carried out. The following epitopes of the vaccine construct were docked with the MHC molecules: TLFPGDDDL (A29/A30/A27), NTLSERISSK (M1/M1/L1)—MHC-I ligands; FFIVVATAAVCLLFI (A30/A31/A28), LSMITMSAFLIVRLN (A35/A36/A33), ASYISCTANSWNVIP (B6/B7/B5), and KIQNVIIDECY (M1/M1/L1)—MHC-II ligands. The tertiary structure of these epitopes of the vaccine construct was predicted with AlphaFold2 [52,53]. The PDB files 6TDS and 1AQD were edited and cleaned to remove heteroatoms, bound peptides, and water molecules.

3. Results

3.1. Selecting Conserved Epitopes

The selected antigens in the virus are illustrated in Figure 1. All the conserved sequences and epitopes of each antigen are given in the Supplementary Materials (Tables S1–S10). The strategy for designing a potentially universal MPXV mRNA vaccine based on the conserved epitopes is given in Figure 2. The ultimately selected epitopes that were used for the vaccine design along with the MPXV multi-epitope universal mRNA vaccine construct and the plasmid vector are given in Figure 3.

3.2. Assessment of Structure

Using RoseTTAFold, the tertiary structure of the protein of the vaccine was predicted and subsequently refined by GalaxyRefine (Figure 4). The tool generated five models out of which the one with the better characteristics was selected for the docking analysis. Of the amino acids, 93.0% were located in the most favored region of the Ramachandran, 6.3% were in the allowed region, 0.4% were in regions that were generously allowed, and 0.4% of the amino acid residues were in the outlier region. These parameters, along with the ERRAT with an overall quality factor of 92.9825, indicate the desired level of the quality of the vaccine product (Figure 5).

3.3. Physicochemical Properties, Allergenicity, and Antigenicity Analyses

According to the physicochemical properties predicted by ProtParam, the vaccine sequence contains 647 amino acid residues, and its MW is 68.97777 kDa. The theoretical pI was computed to be 9.31. The vaccine contains 55 negatively charged residues (Asp + Glu). The instability index was 39.82, which classifies the vaccine as stable. The estimated half-life in vitro was 100 h, while the estimated half-lives in vivo in yeast and E. coli were 20 h and 10 h, respectively. The aliphatic index was 82.32, which implies that the vaccine construct has a high level of thermostability. The vaccine’s GRAVY score was −0.233, which suggests that the protein is hydrophilic.
The allergenicity prediction via the AlgPred server showed that the multi-epitope vaccine was non-allergen, while the server VaxiJen v2.0 with the default parameters and the threshold 0.4 [54] demonstrated its probable antigenicity with the antigenic score 0.4656.

3.4. Conformational B Cell Epitopes

In total, four discontinuous B cell epitopes were predicted with the scores 0.69–0.859. The sizes varied between 44 and 109 residues. The presence of the conformational epitopes in the vaccine product indicates its potential capacity to induce a humoral immune response when the vaccine is administered in vivo. All the predicted discontinuous epitopes are given in Table 2 and illustrated in Figure 6.

3.5. Molecular Docking and Immune Responses following Immune Simulation

The molecular docking analysis of epitopes with the MHC-I and MHC-II molecules revealed the stable binding (Figure 7). The molecular docking results are given in the Supplementary Materials (Tables S11 and S12). In total, three doses were administered four weeks apart for the immune simulation—day 0, day 28, and day 56 with an 8 h offset. As expected, the immune response was also compatible with the immune responses that are elicited in general after in vivo immunization. The immune responses after receiving additional and booster doses were more robust compared to the first shot. The immune responses following the second and third shots were stronger in comparison to the prime immunization. The antigen level was decreased while elevated levels of antibodies, including IgG1, IgG1 + IgG2, IgM, and IgG + IgM, were detected. Antigen abundance peaked at each vaccine injection (Figure 8A). The humoral response surged after each immunization, and the antibody levels remained high for weeks after the last vaccine injection. Following immunization, the activation of the B cell population, resulting in an increased number of B cells producing antibodies, was observed (Figure 8B). The counts of CD8 T-cytotoxic lymphocytes show the ignorable number of anergic cells and activation of T-cytotoxic cells (Figure 8C). The counts of CD4 T-helper lymphocytes demonstrate that the duplication phase starts immediately after each injection (Figure 8D). The cytokine production was also manifested upon immunization; particularly, high levels of IFN-γ, TGF-β, and IL-10 were induced by the infection (Figure 8E). The other parameters of the immune response elicited by the vaccine are presented in the Supplementary Materials (Figure S1).

4. Discussion

Mpox was known to be a rare zoonotic disease caused by the MPXV that usually emerged in rural rainforest regions of African countries [7]. The genus Orthopoxvirus contains four species that are pathogenic to humans: VARV (the causative agent of smallpox), MPXV, VACV, and cowpox. Smallpox and mpox are often life-threatening diseases, whereas vaccinia and cowpox are generally associated with local lesions. MPXV received this name as it was first isolated from laboratory macaque monkeys imported from Singapore to Copenhagen, Denmark, in 1958 [55]; however, it is considered that the host and the main source of MPXV are rodents such as squirrels, Gambian rats, etc. [56]. The first human case was documented in the Democratic Republic of the Congo (DRC) in 1970 while smallpox surveillance was taking place [7]. Particularly, a nine-month-old boy was admitted to the hospital with suspected smallpox. The specimens were sent to the Smallpox Reference Centre of the World Health Organization (WHO) where MPXV was isolated [57]. The symptoms of mpox include fever, malaise, respiratory symptoms, fatigue, lymphadenopathy, headache, and muscle aches, along with the main symptom—a rash similar to pimples and blisters that are itchy and painful and eventually crust over and are healed. The manifestation of mpox is similar to other rash diseases, hence, misdiagnoses often take place. Instead of mpox, mostly, chickenpox caused by Varicella zoster virus is mistakenly diagnosed [15,58,59]. The rush is usually concentrated on the face and extremities including palms and soles of the feet; however, according to the evidence during the current mpox outbreak, it is also localized in the perineal/perianal area as well as on the genitals [60,61]. Although mpox usually resolves by itself, the recent 2022 outbreak demonstrated that it is a life-threatening disease [62]. Indeed, after all, MPXV is categorized as a high-threat virus that belongs to biosafety level 3, according to EU regulations [63]. The possible reasons why MPXV is being spread more than usual include the increased exotic animal trade and international travel [58], the introduction from a single origin with super-spreader events [64], the long period of MPXV cryptic dissemination in humans as well as in animals in non-endemic countries [64], the affected human immunity levels due to the COVID-19 global pandemic while gaining the adaptability by the MPXV, and the waning of the immunity against smallpox in the world population which is caused by terminating the vaccination since 1980 when the disease was announced eradicated by WHO [7,55,56,65]. The latter one seems more reasonable as today only approximately 30% of the world population is vaccinated against smallpox [58] and the majority of infected people are younger than 50 and have never been vaccinated against smallpox. Vaccination to smallpox induced coincident immunity to MPXV, but smallpox eradication and halted vaccination as well as lack of vaccination triggered the MPXV to acquire clinical relevance [66]. Although the smallpox vaccine protects from MPXV [67], there are cases of MPXV infection in patients who are vaccinated against smallpox [18]. Mpox has the potential to grow as a global threat which necessitates the development of a specific vaccine. Moreover, if MPXV is spread easily due to the waned immunity against smallpox, smallpox itself might also re-emerge anytime soon, and that could be a huge plague. In addition, the “forgotten” viruses should get more attention as from time to time they manage to emerge, and in the case of the highly deadly virus, it would cause an unprecedentedly adverse outcome for the human population, e.g., Crimean-Congo hemorrhagic fever (CCHF) virus, Zika, Ebola, etc. Although they are considered to be limited geographically, the possibility of their spreading wider is high. MPXV was also reckoned as geographically limited and did not get much attention; however, in 2003, the outbreak took place outside of Africa for the first time which led to the increased attention to this virus. Moreover, there is a number of cases that are not reported due to the lack of surveillance systems in endemic areas, which most possibly resulted in the expansion of geographic areas [68]. According to the history of epidemics and pandemics, although the number of cases is not increasing recently, the new wave of MPXV can pose a higher risk to the health of the world population. Hence, it is important to design a more specific and universal vaccine against MPXV that will have fewer side effects. As MPXV, VARV, and VACV are from the same genus and their genomes are conserved, it is prudent to develop a vaccine based on the conserved elements. In this study, five antigens—A29, A30, A35, M1, and B6—were selected according to their importance for the viral life cycle. Importantly, these antigens have shown favorable outcomes when they were used in the line of studies and elicited protection from the MPXV challenge. Buchman et al., as well as Heraud et al., have challenged the non-human primates intradermally with the lethal dose of MPXV after the vaccination with a subunit vaccine consisting of VACV membrane proteins A33, B5, L1, and A27 [28,29] (see the names of the corresponding proteins of MPXV and VARV in Table 1), and all the animals survived. Hooper et al. used a DNA vaccine comprising the same antigens of VACV and demonstrated that the immunization with DNA vaccine encoding these four antigens protected rhesus macaques from a severe disease after the lethal MPXV challenge [27]. On the other hand, Hirao et al. used a DNA vaccine encoding eight antigens of VACV—A4, A56, F9, H3, A27, A33, B5, and L1—to immunize cynomolgus macaques and showed that the vaccination elicited protective immunity in the animals [24]. As to the fifth selected antigen in the current study, A28, it is one of at least eight transmembrane proteins in the VACV entry/fusion complex subunit and is evidenced to be a target of neutralizing and protective antibodies in rabbits [69].
The vaccine used for the eradication of smallpox was the live VACV-based vaccine that necessitated special precautions to protect from spreading the VACV from their vaccination spot [70]. Although the current VACV-based vaccine that is currently used against MPXV is a highly effective attenuated Modified Vaccinia Ankara-Bavarian Nordic (MVA-BN) virus [71,72], the development of a potentially universal next-generation vaccine against MPXV, VARV, and similar viruses from the same genus is reasonable. Nucleic acid vaccines, specifically mRNA-based vaccines, have recently garnered significant interest due to their remarkable advantages over other types of nucleic acid and conventional vaccines [73]. mRNA vaccines are highly potent, safe, cost-effective, and rapidly manufactured [19,20,21]. Indeed, the current COVID-19 pandemic has proven the superiority of mRNA vaccines over any other vaccines [74,75,76]. The advantages of mRNA-based immunization include efficacy, safety (no risk of genome integration as it is directly translated into the cytoplasm unlike DNA vaccines), cell-free manufacturing, and fast production. Hence, as mRNA technology represents a promising strategy, designing an mRNA vaccine against MPXV, VARV, and VACV is definitely reasonable [77,78].
Smallpox is the only infectious disease that is considered to be globally eradicated since 1980 [79]. However, the course of pandemic events in world history has revealed that a new pandemic can emerge unexpectedly despite the rapid advancement of science with numerous valuable studies. Thus, even though smallpox has been eliminated, there is a risk of re-emergence. Moreover, in the past, it killed millions of people until the effective vaccine was developed [80]. The current outbreak of mpox is an example of it. In addition, although VACV is not regarded as dangerous, evidently, it causes infection [81], and there is a possibility that it will become more virulent in the future. Accordingly, designing a potentially universal vaccine that can be protective to MPXV as a current emergency, VARV as a potential re-emergency, and VACV as its effective application as a vaccine seems prudent.
Noteworthily, a number of excellent immunoinformatics tools are available for finding the conserved sequences of the selected antigens [82,83,84], for finding experimentally tested epitopes or predicting the epitopes on the IEDB [31], and for the prediction of certain properties and immune response induced by the designed vaccine [30,40,41,42,43,44,45,46,47,49,51,85,86,87]. These in silico analyses save time and allow us to predict the potential outcome of the designed vaccine. Moreover, rather than commencing the time-consuming and expensive in vitro/in vivo experiments directly, using the immunoinformatics analysis that is favorable to formulate the universal MPXV mRNA vaccine and screen the multi-epitope construct containing conserved elements using in silico approaches seems very reasonable. Indeed, a line of research has been dedicated to developing vaccines for COVID-19 [41,88], influenza [89,90,91], and other viruses [92,93,94] using bioinformatics approaches. Until now, there is no study focusing on designing a multi-epitope universal mRNA vaccine against MPXV, VARV, and VACV together. Hence, in the current study, epitopes of T cells that have already undergone experimental testing were chosen for their potential to enhance the favorable results of the proposed mRNA vaccine. On the other hand, B cell epitopes were predicted via IEDB prediction software [32] as there were no experimentally tested B cell epitopes found in the database. As the optimization of the final construct is of great importance and linkers that connect the epitopes to each other play a crucial role, flexible KK [41] and GGGS [35] linkers along with a rigid EAAAK linker [34,38] were used. The flexible linkers used here improved folding and stability in fusion proteins, while the EAAAK attached the peptides to each other to achieve an immunologically active multi-peptide [33].
In silico immunization showed that the vaccine designed for MPXV, VARV, and VACV in this study elicits optimal humoral and cellular immune responses following repeated immunization. In silico immunization with the proposed vaccine stimulated the B cell population, and increased the production of immunoglobulins, CD8 T-cytotoxic and CD4 T-helper cells, memory cells, as well as cytokines. In addition, the molecular docking showed favorable binding to the MHC molecules and other physicochemical parameters of the vaccine demonstrated stability and optimal characteristics of the vaccine.
The in vivo validation of the efficacy and safety of the potentially universal mRNA vaccine, which was designed and tested in silico in this study, is unquestionably necessary. For this, the multi-epitope antigen should be subcloned into an expression vector, transformed into DH5α competent cells, amplified, extracted and purified, linearized, in vitro transcribed, and encapsulated with lipid nanoparticles (LNPs). Then, the animal immunization experiment should be conducted for the assessment of immune responses. The schematic illustration of the steps needed for the development of a potentially universal mRNA vaccine candidate for MPXV, VARV, and VACV is shown in Figure 9.

5. Conclusions

Taken together, it is conceivable that mpox outbreaks will take place more frequently in the future [95]. Thus, more research is urgently needed for developing a specific and universal MPXV vaccine and therapeutics for the prevention and treatment of the disease, respectively. To the best of our knowledge, this is the first study focusing on designing of multi-epitope universal mRNA MPXV vaccine that will be potentially effective against VARV and VACV as well. Indeed, immunoinformatics presents a promising strategy for contributing to the rapid development of mRNA vaccines, which could be achieved through subsequent pre-clinical and clinical investigations. This would enable the prevention, mitigation, and management of potential outbreaks before they escalate into deadly pandemics. Moreover, given that predicting the timing and location of the next pandemic is unfeasible, it is essential to have readily developed vaccine candidates that can be used swiftly to contain outbreaks before they escalate into worldwide pandemics. In summary, the current study demonstrated that the multi-epitope construct for mpox proposed here seems to be an auspicious candidate owing to the potential to induce an immune response against the two clades of MPXV and same family viruses—VACV and VARV—that might lay the groundwork to the alertness for epidemics and pandemics.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/v15051120/s1, Figure S1: Immune simulation results; Table S1: List of T cell epitopes (antigen A29L-A30L-A27L) containing conserved amino acid sequences (origin VACV); Table S2: List of B cell epitopes (antigen A29L-A30L-A27L); Table S3: List of T cell epitopes (antigen A30L/A31L/A28L) containing conserved amino acid sequences (origin VACV); Table S4: List of B cell epitopes (antigen A30L/A31L/A28L); Table S5: List of T cell epitopes (antigen A35R-A36R-A33R) containing conserved amino acid sequences (origin VACV); Table S6: List of B cell epitopes (antigen A35R-A36R-A33R); Table S7: List of T cell epitopes (antigen B6R-B7R-B5R) containing conserved amino acid sequences (origin VACV and VARV); Table S8: List of B cell epitopes (antigen B6R-B7R-B5R); Table S9: List of T cell epitopes (antigen M1R-M1R-L1R) containing conserved amino acid sequences (origin VACV); Table S10: List of B cell epitopes (antigen M1R-M1R-L1R); Table S11: Molecular docking results of the final vaccine constructs with MHC-I; Table S12: Molecular docking results of the final vaccine constructs with MHC-II.

Author Contributions

Conceptualization, N.R.; methodology, N.R. and J.M.; software, N.R. and J.M.; writing—original draft preparation, N.R.; writing—review and editing, D.P. and C.L.; visualization, D.P., S.F., X.Y. and J.L.; supervision, Y.H.; project administration, P.G.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by SHENZHEN SCIENCE AND TECHNOLOGY INNOVATION PROGRAM, grant number KQTD20200909113758004.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets analyzed for this study can be found in the NCBI database [https://www.ncbi.nlm.nih.gov/; accessed on 19 August 2022] and the Immune Epitope Database (IEDB) [https://www.iedb.org/; accessed on 19 August 2022]. The data that support the findings of this study are available from the corresponding authors upon reasonable request. All data created during this study are included in this article. Further inquiries can be directed to the corresponding authors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Centers for Disease Control and Prevention (CDC). 2022 Monkeypox Outbreak Global Map. Available online: https://www.cdc.gov/poxvirus/monkeypox/response/2022/world-map.html (accessed on 6 March 2023).
  2. Vandenbogaert, M.; Kwasiborski, A.; Gonofio, E.; Descorps, S.; Selekon, B.; Andy, A.; Meyong, N.; Ouilibona, R.S.; Gessain, A.; Manuguerra, J.C.; et al. Nanopore Sequencing of a Monkeypox Virus Strain Isolated from a Pustular Lesion in the Central African Republic. Sci. Rep. 2022, 12, 10768. [Google Scholar] [CrossRef] [PubMed]
  3. Shchelkunov, S.N.; Totmenin, A.V.; Safronov, P.F.; Mikheev, M.V.; Gutorov, V.V.; Ryazankina, O.I.; Petrov, N.A.; Babkin, I.V.; Uvarova, E.A.; Sandakhchiev, L.S.; et al. Analysis of the Monkeypox Virus Genome. Virology 2002, 197, 172–194. [Google Scholar] [CrossRef] [PubMed]
  4. Zhu, M.; Ji, J.; Shi, D.; Lu, X.; Wang, B.; Wu, N.; Wu, J. Unusual Global Outbreak of Monkeypox: What Should We Do? Front. Med. 2022, 16, 507–517. [Google Scholar] [CrossRef] [PubMed]
  5. Centers for Disease Control and Prevention (CDC). Past U.S. Cases and Outbreaks. Available online: https://www.cdc.gov/poxvirus/monkeypox/outbreak/us-outbreaks.html (accessed on 16 September 2022).
  6. Gong, Q.; Wang, C.; Chuai, X.; Chiu, S. Monkeypox Virus: A Re-Emergent Threat to Humans. Virol. Sin. 2022, 37, 477–482. [Google Scholar] [CrossRef]
  7. Sklenovská, N.; Ranst, M. Van Emergence of Monkeypox as the Most Important Orthopoxvirus Infection in Humans. Front. Public Health 2018, 6, 241. [Google Scholar] [CrossRef]
  8. Rao, A.K.; Schulte, J.; Chen, T.-H.; Hughes, C.M.; Davidson, W.; Neff, J.M.; Markarian, M.; Delea, K.C.; Wada, S.; Liddell, A.; et al. Monkeypox in a Traveler Returning from Nigeria—Dallas, Texas, July 2021. MMWR. Morb. Mortal. Wkly. Rep. 2022, 71, 509–516. [Google Scholar] [CrossRef]
  9. Costello, V.; Sowash, M.; Gaur, A.; Cardis, M.; Pasieka, H.; Wortmann, G.; Ramdeen, S. Imported Monkeypox from International Traveler, Maryland, USA, 2021. Emerg. Infect. Dis. 2022, 28, 1002–1005. [Google Scholar] [CrossRef]
  10. Haider, N.; Guitian, J.; Simons, D.; Asogun, D.; Ansumana, R.; Honeyborne, I.; Velavan, T.P.; Ntoumi, F.; Valdoleiros, S.R.; Petersen, E.; et al. Increased Outbreaks of Monkeypox Highlight Gaps in Actual Disease Burden in Sub-Saharan Africa and in Animal Reservoirs. Int. J. Infect. Dis. 2022, 122, 107–111. [Google Scholar] [CrossRef]
  11. Martínez, J.I.; Montalbán, E.G.; Bueno, S.J.; Martínez, F.M.; Juliá, A.N.; Díaz, S.; Marín, N.G.; Deorador, E.C.; Forte, A.N.; García, M.A.; et al. Monkeypox Outbreak Predominantly Affecting Men Who Have Sex with Men, Madrid, Spain, 26 April to 16 June 2022. Eurosurveillance 2022, 27, 2200471. [Google Scholar] [CrossRef]
  12. Heskin, J.; Belfield, A.; Milne, C.; Brown, N.; Walters, Y.; Scott, C.; Bracchi, M.; Moore, L.S.; Mughal, N.; Rampling, T.; et al. Transmission of Monkeypox Virus through Sexual Contact—A Novel Route of Infection. J. Infect. 2022, 85, 334–363. [Google Scholar] [CrossRef]
  13. Orviz, E.; Negredo, A.; Ayerdi, O.; Vázquez, A.; Muñoz-Gomez, A.; Monzón, S.; Clavo, P.; Zaballos, A.; Vera, M.; Sánchez, P.; et al. Monkeypox Outbreak in Madrid (Spain): Clinical and Virological Aspects. J. Infect. 2022, 85, 412–417. [Google Scholar] [CrossRef]
  14. Beer, E.M.; Id, V.B.R. A Systematic Review of the Epidemiology of Human Monkeypox Outbreaks and Implications for Outbreak Strategy. PLoS ONE 2019, 13, e0007791. [Google Scholar] [CrossRef] [PubMed]
  15. Jezek, Z.; Szczeniowski, M.; Paluku, K.M. Human Monkeypox: Clinical Features of 282 Patients. J. Infect. Dis. 1987, 156, 293–298. [Google Scholar] [CrossRef] [PubMed]
  16. Anderson, M.G.; Frenkel, L.D.; Homann, S.; Guffey, J. A Case of Severe Monkeypox Virus Disease in an American Child: Emerging Infections and Changing Professional Values. Pediatr. Infect. Dis. J. 2003, 22, 1093–1096. [Google Scholar] [CrossRef]
  17. Shchelkunov, S.N.; Totmenin, A.V.; Babkin, I.V.; Safronov, P.F.; Ryazankina, O.I.; Petrov, N.A.; Gutorov, V.V.; Uvarova, E.A.; Mikheev, M.V.; Sisler, J.R.; et al. Human Monkeypox and Smallpox Viruses: Genomic Comparison. FEBS Lett. 2001, 509, 66–70. [Google Scholar] [CrossRef] [PubMed]
  18. Huhn, G.D.; Bauer, A.M.; Yorita, K.; Graham, M.B.; Sejvar, J.; Likos, A.; Damon, I.K.; Reynolds, M.G.; Kuehnert, M.J. Clinical Characteristics of Human Monkeypox, and Risk Factors for Severe Disease. Clin. Infect. Dis. 2005, 41, 1742–1751. [Google Scholar] [CrossRef] [PubMed]
  19. Liu, C.; Rcheulishvili, N.; Shen, Z.; Papukashvili, D.; Xie, F.; Wang, Z.; Wang, X.; He, Y.; Wang, P.G. Development of an LNP-Encapsulated MRNA-RBD Vaccine against SARS-CoV-2 and Its Variants. Pharmaceutics 2022, 14, 1101. [Google Scholar] [CrossRef]
  20. Pardi, N.; Hogan, M.J.; Porter, F.W.; Weissman, D. MRNA Vaccines-a New Era in Vaccinology. Nat. Rev. Drug Discov. 2018, 17, 261–279. [Google Scholar] [CrossRef]
  21. Kowalzik, F.; Schreiner, D.; Jensen, C.; Teschner, D.; Gehring, S.; Zepp, F. MRNA-Based Vaccines. Vaccines 2021, 9, 390. [Google Scholar] [CrossRef]
  22. Pütz, M.M.; Midgley, C.M.; Law, M.; Smith, G.L. Quantification of Antibody Responses against Multiple Antigens of the Two Infectious Forms of Vaccinia Virus Provides a Benchmark for Smallpox Vaccination. Nat. Med. 2006, 12, 1310–1315. [Google Scholar] [CrossRef]
  23. Wolffe, E.J.; Vijaya, S.; Moss, B. A Myristylated Membrane Protein Encoded by the Vaccinia Virus L1R Open Reading Frame Is the Target of Potent Neutralizing Monoclonal Antibodies. Virology 1995, 211, 53–63. [Google Scholar] [CrossRef] [PubMed]
  24. Hirao, L.A.; Draghia-akli, R.; Prigge, J.T.; Yang, M.; Satishchandran, A.; Wu, L.; Hammarlund, E.; Khan, A.S.; Babas, T.; Rhodes, L.; et al. Multivalent Smallpox DNA Vaccine Delivered by Intradermal Electroporation Drives Protective Immunity in Nonhuman Primates Against Lethal Monkeypox Challenge. J. Infect. Dis. 2011, 203, 95–102. [Google Scholar] [CrossRef]
  25. Gilchuk, I.; Gilchuk, P.; Sapparapu, G.; Lampley, R.; Singh, V.; Kose, N.; Blum, D.L.; Hughes, L.J.; Satheshkumar, P.S.; Townsend, M.B.; et al. Cross-Neutralizing and Protective Human Antibody Specificities to Poxvirus Infections. Cell 2016, 167, 684–694.e9. [Google Scholar] [CrossRef] [PubMed]
  26. Moss, B. Smallpox Vaccines: Targets of Protective Immunity. Immunol. Rev. 2011, 239, 8–26. [Google Scholar] [CrossRef] [PubMed]
  27. Hooper, J.W.; Thompson, E.; Wilhelmsen, C.; Zimmerman, M.; Ichou, M.A.; Steffen, S.E.; Schmaljohn, C.S.; Schmaljohn, A.L.; Jahrling, P.B. Smallpox DNA Vaccine Protects Nonhuman Primates against Lethal Monkeypox. J. Virol. 2004, 78, 4433–4443. [Google Scholar] [CrossRef] [PubMed]
  28. Buchman, G.W.; Cohen, M.E.; Xiao, Y.; Richardson-Harman, N.; Silvera, P.; DeTolla, L.J.; Davis, H.L.; Eisenberg, R.J.; Cohen, G.H.; Isaacs, S.N. A Protein-Based Smallpox Vaccine Protects Non-Human Primates from a Lethal Monkeypox Virus Challenge. Vaccine 2010, 28, 6627–6636. [Google Scholar] [CrossRef] [PubMed]
  29. Heraud, J.; Edghill-smith, Y.; Ayala, V.; Kalisz, I.; Parrino, J.; Kalyanaraman, V.S.; Manischewitz, J.; King, L.R.; Hryniewicz, A.; Trindade, C.J.; et al. Subunit Recombinant Vaccine Protects against Monkeypox. J. Immunol. 2022, 177, 2552–2564. [Google Scholar] [CrossRef]
  30. Waterhouse, A.M.; Procter, J.B.; Martin, D.M.A.; Clamp, M.; Barton, G.J. Jalview Version 2-a Multiple Sequence Alignment Editor and Analysis Workbench. Bioinformatics 2009, 25, 1189–1191. [Google Scholar] [CrossRef]
  31. Vita, R.; Mahajan, S.; Overton, J.A.; Dhanda, S.K.; Martini, S.; Cantrell, J.R.; Wheeler, D.K.; Sette, A.; Peters, B.; Immune, T.; et al. The Immune Epitope Database (IEDB): 2018 Update. Nucleic Acids Res. 2019, 47, 339–343. [Google Scholar] [CrossRef]
  32. Jespersen, M.C.; Peters, B.; Nielsen, M.; Marcatili, P. BepiPred-2.0: Improving Sequence-Based B-Cell Epitope Prediction Using Conformational Epitopes. Nucleic Acids Res. 2017, 45, W24–W29. [Google Scholar] [CrossRef]
  33. Rahmani, A.; Baee, M.; Saleki, K.; Moradi, S.; Nouri, H.R. 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. 2021, 40, 6097–6113. [Google Scholar] [CrossRef] [PubMed]
  34. Tarrahimofrad, H.; Rahimnahal, S.; Zamani, J.; Jahangirian, E.; Aminzadeh, S. Designing a Multi-epitope Vaccine to Provoke the Robust Immune Response against Influenza A H7N9. Sci. Rep. 2021, 11, 24485. [Google Scholar] [CrossRef] [PubMed]
  35. Sanches, R.C.O.; Tiwari, S.; Ferreira, L.C.G.; Oliveira, F.M.; Lopes, M.D.; Passos, M.J.F.; Maia, E.H.B.; Taranto, A.G.; Kato, R.; Azevedo, V.A.C.; et al. Immunoinformatics Design of Multi-Epitope Peptide-Based Vaccine Against Schistosoma Mansoni Using Transmembrane Proteins as a Target. Front. Immunol. 2021, 12, 621706. [Google Scholar] [CrossRef] [PubMed]
  36. Ahammad, I.; Lira, S.S. Designing a Novel MRNA Vaccine against SARS-CoV-2: An Immunoinformatics Approach. Int. J. Biol. Macromol. 2020, 162, 820–837. [Google Scholar] [CrossRef]
  37. Yi, L.; Lee, Y.; Izzard, L.; Hurt, A.C. A Review of DNA Vaccines Against Influenza. Front. Immunol. 2018, 9, 1568. [Google Scholar] [CrossRef]
  38. Nezafat, N.; Ghasemi, Y.; Javadi, G.; Khoshnoud, M.J.; Omidinia, E. A Novel Multi-Epitope Peptide Vaccine against Cancer: An in Silico Approach. J. Theor. Biol. 2014, 349, 121–134. [Google Scholar] [CrossRef]
  39. Boyoglu-Barnum, S.; Ellis, D.; Gillespie, R.A.; Hutchinson, G.B.; Park, Y.; Moin, S.M.; Acton, O.J.; Ravichandran, R.; Murphy, M.; Pettie, D.; et al. Quadrivalent Influenza Nanoparticle Vaccines Induce Broad Protection. Nature 2021, 592, 623–628. [Google Scholar] [CrossRef]
  40. Saha, S.; Raghava, G.P.S. AlgPred: Prediction of Allergenic Proteins and Mapping of IgE Epitopes. Nucleic Acids Res. 2006, 34, 202–209. [Google Scholar] [CrossRef]
  41. Dong, R.; Chu, Z.; Yu, F.; Zha, Y. Contriving Multi-Epitope Subunit of Vaccine for COVID-19: Immunoinformatics Approaches. Front. Immunol. 2020, 11, 1784. [Google Scholar] [CrossRef]
  42. Doytchinova, I.A.; Flower, D.R. Identifying Candidate Subunit Vaccines Using an Alignment-Independent Method Based on Principal Amino Acid Properties. Vaccine 2007, 25, 856–866. [Google Scholar] [CrossRef]
  43. Doytchinova, I.A.; Flower, D.R. VaxiJen: A Server for Prediction of Protective Antigens, Tumour Antigens and Subunit Vaccines. BMC Bioinform. 2007, 8, 4. [Google Scholar] [CrossRef] [PubMed]
  44. Gasteiger, E.; Hoogland, C.; Gattiker, A.; Duvaud, S.; Wilkins, M.R.; Appel, R.D.; Bairoch, A. Protein Identification and Analysis Tools on the ExPASy Server; The Proteomics Protocols Handbook. Hum. Press 2005, 112, 531–552. [Google Scholar] [CrossRef]
  45. Baek, M.; Dimaio, F.; Anishchenko, I.; Dauparas, J.; Ovchinnikov, S.; Lee, G.R.; Wang, J.; Cong, Q.; Kinch, L.N.; Schaeffer, R.D.; et al. Accurate Prediction of Protein Structures and Interactions Using a 3-Track Neural Network. Science 2021, 373, 871–876. [Google Scholar] [CrossRef] [PubMed]
  46. Heo, L.; Park, H.; Seok, C. GalaxyRefine: Protein Structure Refinement Driven by Side-Chain Repacking. Nucleic Acids Res. 2013, 41, 384–388. [Google Scholar] [CrossRef] [PubMed]
  47. Nelson, D.L.; Cox, M.M.; Hoskins, A.A. Lehninger Principles of Biochemistry, 8th ed.; Macmillan Learning: Austin, TX, USA, 2021; ISBN 9781429234146. [Google Scholar]
  48. Ponomarenko, J.; Bui, H.H.; Li, W.; Fusseder, N.; Bourne, P.E.; Sette, A.; Peters, B. ElliPro: A New Structure-Based Tool for the Prediction of Antibody Epitopes. BMC Bioinform. 2008, 9, 514. [Google Scholar] [CrossRef] [PubMed]
  49. Rapin, N.; Lund, O.; Castiglione, F. Immune System Simulation Online. Bioinformatics 2011, 27, 2013–2014. [Google Scholar] [CrossRef] [PubMed]
  50. Rapin, N.; Lund, O.; Bernaschi, M.; Castiglione, F. Computational Immunology Meets Bioinformatics: The Use of Prediction Tools for Molecular Binding in the Simulation of the Immune System. PLoS ONE 2010, 5, e9862. [Google Scholar] [CrossRef]
  51. Kozakov, D.; Hall, D.R.; Xia, B.; Porter, K.A.; Padhorny, D.; Yueh, C.; Beglov, D.; Vajda, S.; Biology, Q. The ClusPro Web Server for Protein-Protein Docking. Nat. Protoc. 2017, 12, 255–278. [Google Scholar] [CrossRef]
  52. Mirdita, M.; Schütze, K.; Moriwaki, Y.; Heo, L. ColabFold: Making Protein Folding Accessible to All. Nat. Methods 2022, 19, 679–682. [Google Scholar] [CrossRef]
  53. Jumper, J.; Evans, R.; Pritzel, A.; Green, T.; Figurnov, M.; Ronneberger, O.; Tunyasuvunakool, K.; Bates, R.; Žídek, A.; Potapenko, A.; et al. Highly Accurate Protein Structure Prediction with AlphaFold. Nature 2021, 596, 583–589. [Google Scholar] [CrossRef]
  54. Ghafouri, F.; Cohan, R.A.; Samimi, H.; Hosseini, A.; Naderi, M.; Noorbakhsh, F.; Haghpanah, V. Development of a Multiepitope Vaccine Against SARS-CoV-2: Immunoinformatics Study. JMIR Bioinforma. Biotechnol. 2022, 3, e36100. [Google Scholar] [CrossRef] [PubMed]
  55. Xiang, Y.; White, A. Monkeypox Virus Emerges from The Shadow of Its More Infamous Cousin: Family Biology Matters. Emerg. Microbes Infect. 2022, 11, 1768–1777. [Google Scholar] [CrossRef] [PubMed]
  56. Durski, K.N.; Mccollum, A.M.; Nakazawa, Y.; Petersen, B.W.; Reynolds, M.G.; Briand, S.; Djingarey, M.H.; Olson, V.; Damon, I.K.; Khalakdina, A. Emergence of Monkeypox—West and Central Africa, 1970–2017. MMWR. Morb. Mortal. Wkly. Rep. 2018, 67, 306–310. [Google Scholar] [CrossRef] [PubMed]
  57. Ladnyj, D.; Ziegler, P.; Kima, E. A Human Infection Caused by Monkeypox Virus in Basankusu Territory, Democratic Republic of the Congo. Bull. World Health Organ. 1972, 46, 593–597. [Google Scholar]
  58. Kmiec, D.; Kirchhoff, F. Monkeypox: A New Threat? Int. J. Mol. Sci. 2022, 23, 7866. [Google Scholar] [CrossRef] [PubMed]
  59. Osadebe, L.; Hughes, C.M.; Lushima, R.S.; Kabamba, J.; Nguete, B.; Malekani, J.; Pukuta, E.; Karhemere, S.; Tamfum, J.M.; Okitolonda, E.W.; et al. Enhancing Case Definitions for Surveillance of Human Monkeypox in the Democratic Republic of Congo. PLoS Negl. Trop. Dis. 2017, 11, e0005857. [Google Scholar] [CrossRef] [PubMed]
  60. Pfäfflin, F.; Wendisch, D.; Scherer, R.; Jürgens, L.; Godzick, G.; Eva, N.; Tober, P.; Miriam, L.; Stegemann, S.; Max, V.; et al. Monkeypox In-Patients with Severe Anal Pain. Infection 2022, 51, 483–487. [Google Scholar] [CrossRef]
  61. Jang, Y.R.; Lee, M.; Shin, H.; Kim, J.; Choi, M.; Kim, Y.M.; Lee, M.J.; Kim, J.; Na, H.K.; Kim, J.Y.; et al. The First Case of Monkeypox in the Republic of Korea. Infect. Dis. Microbiol. Parasitol. 2022, 37, e224. [Google Scholar] [CrossRef]
  62. Meyer, H.; Perrichot, M.; Stemmler, M.; Emmerich, P.; Schmitz, H.; Varaine, F.; Shungu, R.; Tshioko, F.; Formenty, P. Outbreaks of Disease Suspected of Being Due to Human Monkeypox Virus Infection in the Democratic Republic of Congo in 2001. J. Clin. Microbiol. 2002, 40, 2919–2921. [Google Scholar] [CrossRef]
  63. Tian, D.; Zheng, T. Comparison and Analysis of Biological Agent Category Lists Based on Biosafety and Biodefense. PLoS ONE 2014, 9, e101163. [Google Scholar] [CrossRef]
  64. Isidro, J.; Borges, V.; Pinto, M.; Sobral, D.; Santos, J.D.; Nunes, A.; Mixão, V.; Ferreira, R.; Santos, D.; Duarte, S.; et al. Phylogenomic Characterization and Signs of Microevolution in the 2022 Multi-Country Outbreak of Monkeypox Virus. Nat. Med. 2022, 28, 1569–1572. [Google Scholar] [CrossRef] [PubMed]
  65. Weaver, J.R.; Isaacs, S.N. Monkeypox Virus and Insights into Its Immunomodulatory Proteins. Immunol. Rev. 2009, 225, 96–113. [Google Scholar] [CrossRef] [PubMed]
  66. Nguyen, P.; Ajisegiri, W.S.; Costantino, V.; Chughtai, A.A.; Macintyre, C.R. Reemergence of Human Monkeypox and Declining Population Immunity in the Context of Urbanization, Nigeria, 2017–2020. Emerg. Infect. Dis. 2021, 27, 1007–1014. [Google Scholar] [CrossRef]
  67. Burki, T. Investigating Monkeypox. Lancet 2022, 399, 2254–2255. [Google Scholar] [CrossRef] [PubMed]
  68. Yang, Z. Monkeypox: A Potential Global Threat? J. Med. Virol. 2022, 94, 4034–4036. [Google Scholar] [CrossRef] [PubMed]
  69. Nelson, G.E.; Sisler, J.R.; Chandran, D.; Moss, B. Vaccinia Virus Entry/Fusion Complex Subunit A28 Is a Target of Neutralizing and Protective Antibodies. Virology 2008, 380, 394–401. [Google Scholar] [CrossRef]
  70. Rizk, J.G.; Lippi, G.; Henry, B.M.; Forthal, D.N.; Rizk, Y. Prevention and Treatment of Monkeypox. Drugs 2022, 82, 957–963. [Google Scholar] [CrossRef] [PubMed]
  71. Kandeel, M.; Morsy, M.A.; El-lateef, H.M.A.; Marzok, M.; El-Beltagi, H.S.; Khodair, K.M.A.; Albokhadaim, I.; Venugopala, K.N. Efficacy of the Modified Vaccinia Ankara Virus Vaccine and the Replication-Competent Vaccine ACAM2000 in Monkeypox Prevention Mahmoud. Int. Immunopharmacol. 2023, 119, 110206. [Google Scholar] [CrossRef]
  72. Food and Drug Administration (FDA). Use of JYNNEOS (Smallpox and Monkeypox Vaccine, Live, Nonreplicating) for Preexposure Vaccination of Persons at Risk for Occupational Exposure to Orthopoxviruses: Recommendations of the Advisory Committee on Immunization Practices—United States. 2022. Available online: https://www.fda.gov/media/131078/download#:~:text=JYNNEOS is an attenuated%2C live,prevention of smallpox and monkeypox (accessed on 4 May 2023).
  73. Turner, J.S.; Halloran, J.A.O.; Kalaidina, E.; Kim, W.; Schmitz, A.J.; Zhou, J.Q.; Lei, T.; Thapa, M.; Chen, R.E.; Case, J.B.; et al. SARS-CoV-2 MRNA Vaccines Induce Persistent Human Germinal Centre Responses. Nature 2021, 596, 109–113. [Google Scholar] [CrossRef]
  74. Walsh, E.E.; Frenck, R.W.; Falsey, A.R.; Kitchin, N.; Absalon, J.; Gurtman, A.; Lockhart, S.; Neuzil, K.; Mulligan, M.J.; Bailey, R.; et al. Safety and Immunogenicity of Two RNA-Based COVID-19 Vaccine Candidates. N. Engl. J. Med. 2020, 383, 2439–2450. [Google Scholar] [CrossRef]
  75. Haas, E.J.; Angulo, F.J.; McLaughlin, J.M.; Anis, E.; Singer, S.R.; Khan, F.; Brooks, N.; Smaja, M.; Mircus, G.; Pan, K.; et al. 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. Lancet 2021, 397, 1819–1829. [Google Scholar] [CrossRef] [PubMed]
  76. Corbett, K.S.; Flynn, B.; Foulds, K.E.; Francica, J.R.; Boyoglu-Barnum, S.; Werner, A.P.; Flach, B.; O’Connell, S.; Bock, K.W.; Minai, M.; et al. Evaluation of the MRNA-1273 Vaccine against SARS-CoV-2 in Nonhuman Primates. N. Engl. J. Med. 2020, 383, 1544–1555. [Google Scholar] [CrossRef]
  77. Chaudhary, N.; Weissman, D.; Whitehead, K.A. MRNA Vaccines for Infectious Diseases: Principles, Delivery and Clinical Translation. Nat. Rev. Drug Discov. 2021, 20, 817–838. [Google Scholar] [CrossRef] [PubMed]
  78. Damase, T.R.; Sukhovershin, R.; Boada, C.; Taraballi, F.; Pettigrew, R.I.; Cooke, J.P. The Limitless Future of RNA Therapeutics. Front. Bioeng. Biotechnol. 2021, 9, 628137. [Google Scholar] [CrossRef] [PubMed]
  79. Osterholm, M.T.; Gellin, B. Confronting 21st-Century Monkeypox. Science 2022, 377, 349. [Google Scholar] [CrossRef]
  80. Belongia, E.A.; Clinic, M. Smallpox Vaccine: The Good, the Bad, and the Ugly. Clin. Med. Res. 2003, 1, 87–92. [Google Scholar] [CrossRef] [PubMed]
  81. Lu, B.; Cui, L.; Gu, M.; Shi, C.; Sun, C.; Zhao, K. Outbreak of Vaccinia Virus Infection. Emerg. Infect. Dis. 2019, 25, 2017–2020. [Google Scholar]
  82. Sievers, F.; Wilm, A.; Dineen, D.; Gibson, T.J.; Karplus, K.; Li, W.; Lopez, R.; Mcwilliam, H.; Remmert, M.; Söding, J.; et al. Fast, Scalable Generation of High-Quality Protein Multiple Sequence Alignments Using Clustal Omega. Mol. Syst. Biol. 2011, 7, 539. [Google Scholar] [CrossRef] [PubMed]
  83. Sievers, F.; Higgins, D.G. Clustal Omega for Making Accurate Alignments of Many Protein Sequences. Protein Sci. 2018, 27, 135–145. [Google Scholar] [CrossRef]
  84. Sievers, F.; Barton, G.J.; Higgins, D.G. Multiple Sequence Alignments. Bioinformatics, 4th ed.; Baxevanis, A.D., Bader, G.D., Wishart, D.S., Eds.; Wiley: Hoboken, NJ, USA, 2020; Volume 227, ISBN 978-1-119-33558-0. [Google Scholar]
  85. González-Pech, R.A.; Stephens, T.G.; Chan, C.X. Commonly Misunderstood Parameters of NCBI BLAST and Important Considerations for Users. Bioinformatics 2019, 35, 2697–2698. [Google Scholar] [CrossRef]
  86. Zhang, J.; Liang, Y.; Zhang, Y. Atomic-Level Protein Structure Refinement Using Fragment Guided Molecular Dynamics Conformation Sampling. Structure 2011, 19, 1784–1795. [Google Scholar] [CrossRef] [PubMed]
  87. Kozakov, D.; Brenke, R.; Comeau, S.R.; Vajda, S. PIPER: An FFT-Based Protein Docking Program with Pairwise Potentials. Proteins 2006, 65, 392–406. [Google Scholar] [CrossRef] [PubMed]
  88. Oliveira, S.C.; Magalhães, M.T.Q. De; Homan, E.J. Immunoinformatic Analysis of SARS-CoV-2 Nucleocapsid Protein and Identi Fi Cation of COVID-19 Vaccine Targets. Front. Immunol. 2020, 11, 587615. [Google Scholar] [CrossRef]
  89. Sharma, S.; Kumari, V.; Vasant, B.; Mukherjee, A.; Pandey, R.; Kondabagil, K. Immunoinformatics Approach for a Novel Multi-Epitope Subunit Vaccine Design against Various Subtypes of Influenza A Virus. Immunobiology 2021, 226, 152053. [Google Scholar] [CrossRef] [PubMed]
  90. Behbahani, M.; Moradi, M.; Mohabatkar, H. In Silico Design of a Multi-epitope Peptide Construct as a Potential Vaccine Candidate for Influenza A Based on Neuraminidase Protein. In Silico Pharmacol. 2021, 9, 36. [Google Scholar] [CrossRef] [PubMed]
  91. Rcheulishvili, N.; Mao, J.; Papukashvili, D.; Liu, C.; Wang, Z.; Zhao, J.; Xie, F.; Pan, X.; Ji, Y.; He, Y.; et al. Designing multi-epitope mRNA construct as a universal influenza vaccine candidate for future epidemic/pandemic preparedness. Int J Biol Macromol. 2023, 226, 885–899. [Google Scholar] [CrossRef] [PubMed]
  92. Ros-Lucas, A.; Correa-Fiz, F.; Bosch-Camós, L.; Rodriguez, F.; Alonso-Padilla, J. Computational Analysis of African Swine Fever Virus Protein Space for the Design of an Epitope-Based Vaccine Ensemble. Pathogens 2020, 9, 1078. [Google Scholar] [CrossRef]
  93. Mahmudul, H.; Shiful, I.; Sourav, C.; Abu, H.M.; Kazi, F.A.; Ziaul, F.J.; Nazmul, H.; Shakhawat, H.F.; Nazmul, H. 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. 2020, 38, 2898–2915. [Google Scholar] [CrossRef]
  94. Ali, A.; Khan, A.; Kaushik, A.C.; Wang, Y.; Ali, S.S.; Junaid, M.; Saleem, S.; Cho, W.C.S.; Mao, X.; Wei, D. Immunoinformatic and Systems Biology Approaches to Predict and Validate Peptide Vaccines against Epstein–Barr Virus (EBV). Sci. Rep. 2019, 9, 720. [Google Scholar] [CrossRef]
  95. Papukashvili, D.; Rcheulishvili, N.; Liu, C.; Wang, X.; He, Y.; Wang, P.G. Strategy of developing nucleic acid-based universal monkeypox vaccine candidates. Front Immunol. 2022, 13, 1050309. [Google Scholar] [CrossRef]
Viruses 15 01120 g001 550
Figure 1. Structure of MPXV and selected viral antigens. Red font indicates the names of the antigens in MPXV, VARV, and VACV. Underline indicates the special emphasis on the name of antigen in MPXV. MPXV, monkeypox virus; VARV, variola virus; VACV, vaccinia virus.
Figure 1. Structure of MPXV and selected viral antigens. Red font indicates the names of the antigens in MPXV, VARV, and VACV. Underline indicates the special emphasis on the name of antigen in MPXV. MPXV, monkeypox virus; VARV, variola virus; VACV, vaccinia virus.
Viruses 15 01120 g001
Viruses 15 01120 g002 550
Figure 2. Strategy of designing a potentially universal MPXV mRNA vaccine. MPXV, monkeypox virus. MPXV, monkeypox virus; VARV, variola virus; VACV, vaccinia virus.
Figure 2. Strategy of designing a potentially universal MPXV mRNA vaccine. MPXV, monkeypox virus. MPXV, monkeypox virus; VARV, variola virus; VACV, vaccinia virus.
Viruses 15 01120 g002
Viruses 15 01120 g003 550
Figure 3. Schematic illustration of the designed mRNA and plasmid. (A) Table of selected conserved epitopes for the mRNA vaccine design. (B) Plasmid scheme. (C) Final design of the multi-epitope mRNA vaccine construct.
Figure 3. Schematic illustration of the designed mRNA and plasmid. (A) Table of selected conserved epitopes for the mRNA vaccine design. (B) Plasmid scheme. (C) Final design of the multi-epitope mRNA vaccine construct.
Viruses 15 01120 g003
Viruses 15 01120 g004 550
Figure 4. Tertiary structure of the protein comprising multiple epitopes.
Figure 4. Tertiary structure of the protein comprising multiple epitopes.
Viruses 15 01120 g004
Viruses 15 01120 g005 550
Figure 5. Ramachandran plot assessing vaccine structure quality. Each dot shows the amino acids. The positioning of these dots reflects the backbone dihedral angles ψ versus φ of the amino acids present in the vaccine product. Triangles indicate glycine, squares indicate prolines, while circles denote non-glycine and non-proline residues.
Figure 5. Ramachandran plot assessing vaccine structure quality. Each dot shows the amino acids. The positioning of these dots reflects the backbone dihedral angles ψ versus φ of the amino acids present in the vaccine product. Triangles indicate glycine, squares indicate prolines, while circles denote non-glycine and non-proline residues.
Viruses 15 01120 g005
Viruses 15 01120 g006 550
Figure 6. Predicted discontinuous epitopes of B cell epitopes in the vaccine product. The subfigure numbers (1–4) correspond to the conformational B cell epitope numbers listed in Table 2.
Figure 6. Predicted discontinuous epitopes of B cell epitopes in the vaccine product. The subfigure numbers (1–4) correspond to the conformational B cell epitope numbers listed in Table 2.
Viruses 15 01120 g006
Viruses 15 01120 g007 550
Figure 7. Representation of predicted molecular docking between epitopes of a multi-epitope vaccine construct and host receptors. The MHC molecules are displayed in cyan color and the epitope is shown in magenta. (A) Docked complex of epitope TLFPGDDDL (A29/A30/A27) and MHC-I. (B) Docked complex of epitope NTLSERISSK (M1/M1/L1) and MHC-I. (C) Docked complex of epitope LSMITMSAFLIVRLN (A35/A36/A33) and MHC-II. (D) Docked complex of epitope ASYISCTANSWNVIP (B6/B7/B5) and MHC-II. (E) Docked complex of epitope KIQNVIIDECY (M1/M1/L1) and MHC-II. (F) Docked complex of epitope FFIVVATAAVCLLFI (A30/A31/A28) and MHC-II.
Figure 7. Representation of predicted molecular docking between epitopes of a multi-epitope vaccine construct and host receptors. The MHC molecules are displayed in cyan color and the epitope is shown in magenta. (A) Docked complex of epitope TLFPGDDDL (A29/A30/A27) and MHC-I. (B) Docked complex of epitope NTLSERISSK (M1/M1/L1) and MHC-I. (C) Docked complex of epitope LSMITMSAFLIVRLN (A35/A36/A33) and MHC-II. (D) Docked complex of epitope ASYISCTANSWNVIP (B6/B7/B5) and MHC-II. (E) Docked complex of epitope KIQNVIIDECY (M1/M1/L1) and MHC-II. (F) Docked complex of epitope FFIVVATAAVCLLFI (A30/A31/A28) and MHC-II.
Viruses 15 01120 g007
Viruses 15 01120 g008 550
Figure 8. Immune responses elicited after in silico immunization. (A) Antibodies produced in response to immune simulation. Immunoglobulin subclasses are given in different colors. (B) Representation of evolution of B cell populations following the immunization with three doses. (C) Generation of CD8 T-cytotoxic lymphocytes induced by antigen exposure. The resting state denotes T-cytotoxic cells not actively engaged in fighting the infection, while the anergic state denotes the unresponsiveness of the CD8 T-cytotoxic cells due to repeated and prolonged exposure to the antigen. (D) Production of CD4 T-helper cells induced by exposure to the vaccine antigen. (E) The cytokine profile in response to the vaccine injections. The inserted graph shows the IL-2 level with the Simpson Index. D indicates a danger signal. (F) Schedule of in silico immunization.
Figure 8. Immune responses elicited after in silico immunization. (A) Antibodies produced in response to immune simulation. Immunoglobulin subclasses are given in different colors. (B) Representation of evolution of B cell populations following the immunization with three doses. (C) Generation of CD8 T-cytotoxic lymphocytes induced by antigen exposure. The resting state denotes T-cytotoxic cells not actively engaged in fighting the infection, while the anergic state denotes the unresponsiveness of the CD8 T-cytotoxic cells due to repeated and prolonged exposure to the antigen. (D) Production of CD4 T-helper cells induced by exposure to the vaccine antigen. (E) The cytokine profile in response to the vaccine injections. The inserted graph shows the IL-2 level with the Simpson Index. D indicates a danger signal. (F) Schedule of in silico immunization.
Viruses 15 01120 g008
Viruses 15 01120 g009 550
Figure 9. A roadmap for developing a potentially universal multi-epitope mRNA vaccine candidate against MPXV. MPXV, monkeypox virus.
Figure 9. A roadmap for developing a potentially universal multi-epitope mRNA vaccine candidate against MPXV. MPXV, monkeypox virus.
Viruses 15 01120 g009
Table
Table 1. Selected antigens among MPXV, VARV, and VACV, and their location in the virus, function, and characteristics.
Table 1. Selected antigens among MPXV, VARV, and VACV, and their location in the virus, function, and characteristics.
Name (MPXV)Name (VARV)Name (VACV)LocationFunction and Characteristics
A29A30A27MVSurface membrane fusion protein;
Binds to cell surface heparan;
Neutralizing antibody target
A30A31A28MVEnvelope protein;
Virus entry into a host;
Cell–cell fusion (syncytial formation);
Neutralizing antibody target
A35A36A33EVEnvelope glycoprotein;
Formation of actin-containing microvilli and cell-to-cell spread of virion;
Neutralizing antibody target;
Target of complement-mediated cytolysis
B6B7B5EVPalmitylated glycoprotein;
Required for efficient cell spread;
Complement control
M1M1L1MVMyristylated surface membrane protein;
Virus entry into a host;
Neutralizing antibody target
Notes: MPXV, monkeypox virus; VARV, variola virus; VACV, vaccinia virus.
Table
Table 2. Discontinuous B cell epitopes predicted by ElliPro.
Table 2. Discontinuous B cell epitopes predicted by ElliPro.
#EpitopesResidue NumberScore
1A:K287, A:E289, A:Q290, A:T291, A:S292, A:K293, A:K294, A:V295, A:S296, A:E297, A:L298, A:Y299, A:N300, A:K301, A:P302, A:L303, A:Y304, A:K305, A:K306, A:E307, A:E308, A:K309, A:N310, A:G311, A:N312, A:T313, A:S314, A:W315, A:N316, A:D317, A:T318, A:V319, A:K320, A:K321, A:P322, A:D323, A:D324, A:E325, A:T326, A:D327, A:L328, A:S329, A:K330, A:L331440.859
2A:I511, A:V512, A:R513, A:L514, A:N515, A:Q516, A:C517, A:M518, A:S519, A:A520, A:N521, A:G522, A:G523, A:G524, A:S525, A:A526, A:S527, A:Y528, A:I529, A:S530, A:C531, A:T532, A:A533, A:N534, A:S535, A:N537, A:I539, A:G546, A:V547, A:I548, A:H549, A:L550, A:S551, A:C552, A:K553, A:S554, A:G555, A:F556, A:I557, A:L558, A:T559, A:G560, A:G561, A:G562, A:S563, A:M564, A:K565, A:T566, A:I567, A:S568, A:V569, A:V570, A:S589, A:V592, A:L593, A:V594, A:C595, A:S596, A:C597, A:N598, A:G599, A:G600, A:G601, A:S602, A:K603, A:I604, A:Q605, A:N606, A:V607, A:I608, A:I609, A:D610, A:E611, A:C612, A:Y613, A:G614, A:G615, A:G616, A:S617, A:A618, A:A619, A:L620, A:F621, A:M622, A:Y623, A:Y624, A:A625, A:K626, A:R627, A:G628, A:G629, A:G630, A:S631, A:N632, A:T633, A:L634, A:S635, A:E636, A:R637, A:I638, A:S639, A:S640, A:K641, A:H642, A:H643, A:H644, A:H645, A:H646, A:H6471090.761
3A:K212, A:D216, A:N217, A:K218, A:K221, A:G224, A:G225, A:T226, A:P227, A:A228, A:K229, A:K230, A:Q231, A:D232, A:V233, A:N234, A:D235, A:T236, A:I237, A:S238, A:D239, A:K240, A:K241, A:G242, A:P243, A:N244, A:N245, A:T246, A:R247, A:K248, A:K249, A:S250, A:T251, A:H252, A:R253, A:K254, A:V255, A:F401, A:I402, A:V403, A:V404, A:A405, A:T406, A:A407, A:A408, A:V409, A:C410, A:L411, A:L412, A:F413, A:I414, A:G415, A:G416, A:G417, A:S418, A:M419, A:N420, A:S421, A:L422, A:S423, A:I424, A:F425, A:F426, A:I427, A:V428, A:V429, A:A430, A:T431, A:A432, A:A433700.733
4A:D9, A:V10, A:A11, A:A12, A:P13, A:H14, A:R15, A:Q16, A:P17, A:L18, A:T19, A:S20, A:S21, A:E22, A:R23, A:I24, A:D25, A:M66, A:A67, A:E68, A:K69, A:D70, A:G71, A:C72, A:F73, A:Q74, A:S75, A:G76, A:F77, A:N78, A:E79, A:E80, A:T81, A:C82, A:L83, A:V84, A:K85, A:I86, A:I87, A:T88, A:L91, A:V120, A:L121, A:I122, A:Q123, A:F124, A:L125, A:Q126, A:K128, A:A129, A:K130, A:N131, A:L132, A:D133, A:A134, A:I135, A:T136, A:T137, A:P138, A:D139, A:P140, A:T141, A:T142, A:A144, A:S145, A:T148, A:K149, A:S176, A:L177, A:R178, A:A179, A:L180, A:R181, A:Q182, A:M183, A:G184, A:G185, A:G186, A:S187, A:A188, A:K189, A:F190, A:V191, A:A192, A:A193, A:W194, A:T195, A:K197, A:A198, A:E201900.69
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Rcheulishvili, N.; Mao, J.; Papukashvili, D.; Feng, S.; Liu, C.; Yang, X.; Lin, J.; He, Y.; Wang, P.G. Development of a Multi-Epitope Universal mRNA Vaccine Candidate for Monkeypox, Smallpox, and Vaccinia Viruses: Design and In Silico Analyses. Viruses 2023, 15, 1120. https://doi.org/10.3390/v15051120

AMA Style

Rcheulishvili N, Mao J, Papukashvili D, Feng S, Liu C, Yang X, Lin J, He Y, Wang PG. Development of a Multi-Epitope Universal mRNA Vaccine Candidate for Monkeypox, Smallpox, and Vaccinia Viruses: Design and In Silico Analyses. Viruses. 2023; 15(5):1120. https://doi.org/10.3390/v15051120

Chicago/Turabian Style

Rcheulishvili, Nino, Jiawei Mao, Dimitri Papukashvili, Shunping Feng, Cong Liu, Xidan Yang, Jihui Lin, Yunjiao He, and Peng George Wang. 2023. "Development of a Multi-Epitope Universal mRNA Vaccine Candidate for Monkeypox, Smallpox, and Vaccinia Viruses: Design and In Silico Analyses" Viruses 15, no. 5: 1120. https://doi.org/10.3390/v15051120

APA Style

Rcheulishvili, N., Mao, J., Papukashvili, D., Feng, S., Liu, C., Yang, X., Lin, J., He, Y., & Wang, P. G. (2023). Development of a Multi-Epitope Universal mRNA Vaccine Candidate for Monkeypox, Smallpox, and Vaccinia Viruses: Design and In Silico Analyses. Viruses, 15(5), 1120. https://doi.org/10.3390/v15051120

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop