- 1Division of HIV/AIDS and Sex-transmitted Virus Vaccines, Institute for Biological Product Control, National Institutes for Food and Drug Control and WHO Collaborating Center for Standardization and Evaluation of Biologicals, Beijing, China
- 2Wuhan Institute of Biological Products Co., Ltd., Wuhan, China
- 3State Key Laboratory of Common Mechanism Research for Major Diseases, Suzhou Institute of Systems Medicine, Chinese Academy of Medical Sciences and Peking Union Medical College, Suzhou, Jiangsu, China
- 4Guangzhou National Laboratory, Guangzhou, China
- 5Research and Development Department, Beijing Yunling Biotechnology Co., Ltd., Beijing, China
- 6Changping Laboratory, Beijing, China
- 7Graduate School of Peking Union Medical College, Beijing, China
Introduction: Hepatitis E virus (HEV), with heightened virulence in immunocompromised individuals and pregnant women, is a pervasive threat in developing countries. A globaly available vaccine against HEV is currently lacking.
Methods: We designed a multi-epitope vaccine based on protein ORF2 and ORF3 of HEV using immunoinformatics.
Results: The vaccine comprised 23 nontoxic, nonallergenic, soluble peptides. The stability of the docked peptide vaccine-TLR3 complex was validated by molecular dynamic simulations. The induction of effective cellular and humoral immune responses by the multi-peptide vaccine was verified by simulated immunization.
Discussion: These findings provide a foundation for future HEV vaccine studies.
1 Introduction
Hepatitis E virus (HEV) causes acute viral hepatitis and is spread mainly by fecal–oral transmission (Dalton et al., 2008; Nan and Zhang, 2016). It is self-limiting or asymptomatic in healthy individuals, but causes high mortality in immunocompromised patients, such as pregnant women. The case fatality rate in pregnant women is up to 30% (Qian et al., 2023). There are eight genotypes of Orthohepevirus A species in the Hepeviridae family, of which four genotypes infect humans (Marion et al., 2020). Genotypes 1 and 2 exclusively infect humans, whereas genotypes 3 and 4 infect not only humans, but also animals, including swine. Therefore, they are considered zoonotic pathogens and have wide host tropism (Smith et al., 2014; Songtanin et al., 2023).
Hepatitis E virus (HEV) is a 7.2-kb positive-strand RNA virus belonging to the family Hepeviridae, and has three open reading frames (ORFs): ORF1, ORF2, and ORF3 (Smith et al., 2014). ORF1 encodes a large nonstructural polyprotein responsible for HEV replication and transcription. ORF2 is a capsid protein with approximately 660 amino acid residues, and most neutralizing sites are located on its surface. Generally, there are two forms of the ORF2 protein, namely intracellular ORF2 and secreted ORF2 (Montpellier et al., 2018; Yin et al., 2018). HEV is a quasi-enveloped hepatovirus, with the membrane supplied by the host ESCRT system. The membrane protects the virus from the neutralizing antibodies of its host (Feng et al., 2014). ORF3 conducts an ion channel function with a hydrophobic sequence, which is considered a viroporin. ORF3 encodes a protein containing 112 amino acids that is presented on the surface of enveloped HEV (Takahashi et al., 2008). The phosphorylated form of the ORF3 protein interacts with the non-glycosylated form of ORF2 (Tyagi et al., 2002). The a motif containing amino acid proline, serine, alanine, proline (PSAP) motif of ORF3 is necessary for the release of the virions from infected cells (Yamada et al., 2009; Ahmad et al., 2011; Ding et al., 2017). Previous studies have demonstrated that ORF3 not only enhances the production of interferon (IFN), but is also involved in the downregulation of the Toll-like receptor 3 (TLR3) and TLR7 downstream signaling pathways (He et al., 2016; Lei et al., 2018). ORF4, which was discovered in genotype 1, is responsible for the enhancement of viral replication of genotype 1 under endoplasmic reticulum (ER) stress, and this was verified by exogenous introduction of ORF4 into genotype 3 virus (Nair et al., 2016; Yadav et al., 2021).
Two existing hepatitis E vaccines have been evaluated in clinical trials. The Hecolin® vaccine antigen contains amino acids 368–606 of the ORF protein (pORF2) of HEV genotype 1, expressed in Escherichia coli. It can provide protection against HEV for approximately 4.5 years (Zhu et al., 2010; Zhang et al., 2015). The recombinant HEV vaccine (rHEV) produced by GlaxoSmithKline (GSK) contains amino acids 112–607 of the ORF2 protein of HEV genotype 1, expressed in baculovirus (Li et al., 2015).
Because the virus is quasi-enveloped, an antibody targeting ORF2 will not fully detect the enveloped virus in cell culture or serum, unless the virus is treated with sufficient detergent and protease (Takahashi et al., 2008). It has been reported that an ORF3 antibody can neutralize this enveloped virus to some extent, but not the virus in feces (Takahashi et al., 2008). A bacterially expressed ORF3 peptide effectively reduced the viral titer, reduced the duration of viremia and fecal shedding, and even partly prevented experimental hepatitis induced with two HEV genotypes: genotype 1 and 4 (Ma et al., 2009). Upon immunization with the bacterially expressed ORF3 protein of strains VaHEV from laying hens and YT-aHEV from broilers, the nucleic acid of avian HEV detected in cloacal swabs vanished within 7 days of infection. This suggests that the immunized antibodies inhibit the proliferation of the virus and provide protection against infection (Syed et al., 2017). HEV ORF3 equipped with the myotropic adeno-associated virus vector (AAVMYO3) induced dose-dependent ORF3 antibodies in mice and the serum immunized by ORF3 inhibited infection by enveloped HEV in vivo by 50% compared with that of the control (Maurer et al., 2022). Humans infected with HEV via blood transfusion have been reported in many countries (Takahashi et al., 2010).
In regard to the protective roles of ORF3 for enveloped HEV, we first introduce ORF3 as a candidate for multi-epitopic vaccine design.
Immunoinformatic methods have been successfully applied to multiple types of vaccines targeting different pathogens, including HPV, PEDV, Mycobacterium tuberculosis, SARS-CoV-2, and HEV (Behmard et al., 2020; Samavia et al., 2022; Arega et al., 2023; Hou et al., 2023; Ikram et al., 2023; Tîrziu et al., 2023). Multi-epitope vaccines designed with this method will stimulate the CD4+, CD8+ T-cell response and the B-cell response.
In this study, we screened peptides from the ORF2 and ORF3 proteins that stimulated T- and B-cell responses. We also assessed various aspects of the polypeptide vaccine, including its immunogenicity, allergenicity, and physiochemical properties. This study should provide a universal vaccine targeting ORF2 and ORF3 that will eliminate both enveloped and non-enveloped HEV.
2 Materials and methods
The flowchart for the design of the multi-epitope vaccine is shown in Figure 1.
2.1 HEV protein sequences
The amino acid sequences of different HEV proteins, including 961 ORF2 sequences and 746 ORF3 sequences of genotypes 1, 3, and 4, were retrieved from the National Center of Biotechnology Information protein database.
2.2 T-cell epitope prediction, including cytotoxic T-lymphocyte (CTL) epitopes and helper T-lymphocyte (HTL) epitopes
The Immune Epitope Database and Analysis Resource (IEDB)1 was used to predict the CD8 T-cell epitopes of the ORF2 and ORF3 proteins. NetMHCPan v4.1 was the prediction method used. Human was selected as the source species. A percentile rank of < 0.5 was considered the threshold (Reynisson et al., 2020).
Immune Epitope Database (IEDB) was also used to predict the MHC-II-binding (HTL) epitopes (15-mer) of the ORF2 and ORF3 proteins against human Human leukocyte Antigens (HLAs), such as HLA-DRB1*01:01, HLA-DRB1*03:01, HLA-DRB1*04:01, HLA-DRB1*04:05, HLA-DRB1*07:01, HLA-DRB1*08:02, HLA-DRB1*09:01, HLA-DRB1*11:01, HLA-DRB1*12:01, HLA-DRB1*13:02, and HLA-DRB1*15:01, with NN-align, SMM-align, CombLib, and Sturniolo. A percentile rank of < 2 was determined as the threshold.
2.3 B-cell epitope prediction
Three different methods were used to predict B-cell epitopes: the Emini surface accessibility scale, the Kolaskar and Tongaonkar (1990) antigenicity scale, and BepiPred v2.0 (Emini et al., 1985; Kolaskar and Tongaonkar, 1990; Closter et al., 2017). Duplicates in the epitopes predicted with the three methods were removed, and the remaining epitopes were used as candidate B-cell epitopes.
2.4 Assessment of screened epitopes for antigenicity, allergenicity, and toxicity
The antigenic potential of the T- and B-cell epitopes was predicted with VaxiJen v2.0, applying a threshold value of 0.5 (Doytchinova and Flower, 2007). The predicted T- and B-cell epitopes were then further evaluated in terms of their toxicity and allergenicity using the ToxinPred (Gupta et al., 2013) and AllergenFP v2.0 (Ivan et al., 2014) servers, respectively.
2.5 Construction of the polypeptide vaccine
Any multi-epitope vaccine constructed must be neither allergenic nor toxic. Good solubility of the vaccine when expressed at high levels must also be taken into account. Appropriate linkers, such as common rigid linker (EAAAK), fexible connecting peptide (GPGPG), The KK linker is a target sequence of the lysosomal protease histone B, which is one of the important proteases used for antigen processing in the context of MHC-II antigen presentation (KK), and AAY linkers act as cleavage sites for proteasomes in mammalian cells, contributing to the formation of natural epitopes and preventing the formation of “ligand epitopes”, thereby increasing the immunogenicity of epitope-based vaccines (AAY), were added to link the screened CTL, HTL and linear B-lymphocyte (LBL) epitopes. Human β-defensin 3 and LL37 were also added to the N-terminus of the vaccine sequence with the EAAAK linker, to increase the immunogenic capacity of the multi-epitope vaccine (Lande et al., 2021; Samad et al., 2022).
2.6 Immunogenic, allergenic, and physiochemical properties of the multi-peptide vaccine
The antigenicity of the multi-epitope vaccine polypeptide was predicted with the VaxiJen v2.0 tool, with a threshold value of 0.456. The allergenicity of the vaccine was analyzed with AllerTOP v.2.0 (Dimitrov et al., 2014). The ProtParam server was used to assess the physical and chemical properties of the construct (Majid and Andleeb, 2019), including the amino acid composition, molecular weight, theoretical isoelectric point (pI), grand average of hydropathicity (GRAVY), aliphatic and instability index, and in vitro and in vivo half-lives.
2.7 Vaccine structure modeling, refinement, and validation
The SOPMA server was used to predict the secondary structure of the polypeptide vaccine (Geourjon and Deléage, 1995). AlphaFold2 was used to model the tertiary structure of the vaccine. The structure was refined with the Galaxy server and validated with methods such as Ramachandran plots, and ProSA-web (Wiederstein and Sippl, 2007; Junsu et al., 2012; Wang et al., 2016).
2.8 Molecular docking and molecular dynamic (MD) simulation
The docking of the vaccine structure to TLR3 (PDB: 1ZIW) was performed by the HADDOCK server (van Zundert et al., 2016; Honorato et al., 2021). Molecular simulation was performed with GROMACS (version 2018.6) at 30 ns for 150,000 steps. Spc216.gro water was used as the solvent in which to simulate the proteins (Abraham et al., 2015). The addition of Na+ and Cl– ions with the GROMACS tool Genion was used to neutralize the overall charge (Ikram et al., 2023). The energy of the protein structure was minimized in 50,000 steps and the protein was equilibrated with NVT isothermal–isochoric and NPT ensembles in 50,000 steps. MD simulation at 1500,000 steps was used.
2.9 Trial immune simulation in silico
Immune simulation was performed with the C-ImmSim server, which simulates the immune response with both systems’ biology techniques and information supported by data-driven predictive methods (Nicolas et al., 2010). The simulation was performed in 105 steps, and the default option was used for all HLA types of the host. Three different sets of injections were designed: virus alone, vaccine alone, and vaccine followed by virus. The shared sequence of the ORF2 and ORF3, respectively, was used as a substitute for the viral antigen sequence.
3 Results
3.1 T-cell and B-cell epitopes of HEV
The final screened CTL antigenic (MHC-I) epitopes, HTL (MHC-II) epitopes, and LBL epitopes in the ORF2 and ORF3 proteins of HEV are listed in Tables 1, 2. After they were checked with a series of server tools, the antigenic, nontoxic, and nonallergenic T-cell and B-cell epitopes that induce the cytokines interferon-gamma (IFN-γ), interleukin 4 (IL4), and IL10 were selected for the construction of the multi-epitope vaccine.
Table 1. Cytotoxic T-lymphocyte (CTL) epitopes and helper T-lymphocyte (HTL) epitopes in the ORF2 and ORF3 proteins selected for the polypeptide vaccine.
Table 2. Linear B-lymphocyte (LBL) epitopes in the ORF2 and ORF3 proteins selected for the polypeptide vaccine.
3.2 Construction of the multi-epitope vaccine polypeptide
In total, nine CTL, eight HTL, and six LBL epitopic peptides were fused together with the GPGPG, KK, and AAY linkers, respectively, to create the multi-epitope vaccine construct. To boost the immunogenicity of the construct, human β-defensin 3 (GII NTLQKYYCRVRGGRCAVLSCLPKEEQIGKCSTRGRKCCRRKK) and LL37 (LLGDFFRKSKEKIGKEFKRIVQRIKDFLRNLVPRTES) were added as adjuvants to the amino terminus of the polypeptide, joined to the first CTL epitope with the EAAAK linker. An additional start codon was added at the start of the sequence. The primary structure of the multi-epitope subunit vaccine construct included a total of 433 amino acids. The sequence is as follows: MGIINTLQKYYCRVRGGRCAVLSCLPKEEQIGKCSTRGRKCCR RKKEAAAKLLGDFFRKSKEKIGKEFKRIVQRIKDFLRNLVPRTE SEAAAKAVPAVVSGVGPGPGRLHYRNQGWGPGPGTEASNYA QYGPGPGATEASNYAQYGPGPGIMATEASNYGPGPGTPYTG ALGLGPGPGRYRPLVPNAGPGPGGLLDFALELGPGPGVVRATI RYRGPGPGRGAILRRQYNLSTSPLKKATIRYRPLVPNAVKKAL GLLDFALELEKKAILRRQYNLSTSPLTKKATIRYRPLVPNAVGK KLRRQYNLSTSPLTSKKTIRYRPLVPNAVGGYKKYRVVRATIRY RPLKKPSPSPIAAYDYDNQHEQDRPTPSPAPSRPFAAYALGLLD FALELAAYRRQYNLSAAYENAQQDAAYDYDNQHEQDRPTPS PAPSR.
3.3 Immunogenicity, allergenicity, and physiochemical properties of the vaccine candidate
The assessment of the immunogenic, allergenic and solubility characteristics of the construct suggested that the multi-epitope vaccine designed here was a priority (Table 3). The calculated molecular weight was 47.4 kDa and the pI 10.20, indicating the basic nature of the vaccine. The instability index, 42.80, indicated its good stability. The GRAVY value was −0.636, indicating its hydrophilic characteristics, which allows it to connect with other proteins. The aliphatic index of the peptide was 71.76, indicating the high thermostability of the vaccine. The estimated half-life was 30 h in mammalian reticulocytes in vitro, and > 20 h in yeast in vivo, and > 10 h in E. coli in vivo.
3.4 Secondary structure and three-dimensional (3D) structure of multi-epitope vaccine
The secondary structure of the construct showed that the polypeptide contained 49.65% coil, 31.18% helix, 16.40% strand, and 2.77% beta turn (Figure 2). The sequence of the multi-epitope vaccine was submitted to AlphaFold2. AlphaFold2 uses the predicted local distance difference test loss (pLDDT) to directly derive the per-position confidence of predicted models (Figures 3A, B). The pLDDT score of this model was 75.45. The 3D structure was confirmed with Pro-SA with a Z-score of −4.39, indicating the good quality of the vaccine. On Ramachandran plots, the percentage of residues in the most-favored regions was 93.2%, and only 1.2% residues were in the region that was not allowed. Taken together, these parameters indicate the high quality of the multi-peptide vaccine (Figures 3C, D).
Figure 3. Predicted 3D structure of the multi-epitope vaccine and its validation with Pro-SA and Ramachandran plots. (A) Predicted structure of multi-epitope vaccine. (B) TLR3-docked structure with a HADDOCK score of –60.97 and a Z-score of –2.4. (C) Pro-SA with a Z-score of –4.39. (D) Ramachandran plots with 93.2% of residues in most-favored regions.
3.5 Docking the vaccine structure to TLR3 and molecular dynamic (MD) simulation
HADDOCK clustered 160 structures in 11 clusters, which represented 80% of the water-refined models generated by HADDOCK. The best structure of the top 10 clusters was selected according to its HADDOCK score of −60.97 and Z-score of −2.4 (Figure 3B). The more negative the Z-score, the more optimal the structure.
Molecular dynamic (MD) simulation of the complex formed between the polypeptide vaccine and TLR3 was performed. The stability of the docked complex was reflected in the root mean square deviation (RMSD), which had an average value 1.3 nm (Figure 4A). The root mean square fluctuation (RMSF) was between 0.2 and 2.0 nm, with a mean value 0.6 nm (Figure 4B), reflecting the good stability of the whole complex. The high-fluctuation residues (including R666, G730, G807, R931, D995, and Q1051) were distributed in the flexible loop area.
Figure 4. Molecular dynamic simulations. (A) RMSD, with an average value of 1.3 nm. (B) RMSF, with a mean value of 0.6 nm.
3.6 Simulated immunization with C-ImmSim
The assessment of antibody levels provides crucial insights into the efficacy of vaccines and the immune response to pathogens. The role of vaccines in influencing the immune response to infection can be clarified by examining the dynamics of CTL and HTL populations. Our results showed increased immunoglobulin M (IgM) and IgG after 20 days of stimulation with antigen, whereas increases of IgG2 were less obvious (Figure 5A). B memory (y2) and B-isotype IgM levels were elevated (Figure 5B). In general, the levels of TH memory (y2), TC, and TH cells were also increased (Figures 5C, D). The IFN-γ level peaked after 10 days of stimulation (Figure 5E).
Figure 5. Immunization simulated with C-ImmSim. (A) The virus, immunoglobulins, and immunocomplexes. IgM and IgG levels clearly increased, whereas IgG2 levels increased less obviously. (B) B memory (y2) and B-isotype IgM levels increased. (C,D) TH memory (y2), TC, and TH cells generally showed an increase. (E) IFN-γ levels increased notably after stimulation.
4 Discussion
Hepatitis E virus (HEV) is a zoonotic virus that is transmitted through fecal and oral transmission. There is currently no universal vaccine available for HEV (Dalton et al., 2008; Nan and Zhang, 2016). Traditional vaccines targeting ORF2 have been developed for many years (Zhu et al., 2010; Zhang et al., 2015). Multi-epitope vaccines are an emerging area of research interest. An advantage of multi-epitope vaccines is that they reduce adverse reactions (i.e., allergy) and the antigenic load, eliciting a more specific immune response to conserved epitopes (Oli et al., 2020). A combination of the dominant epitopes of different antigens is superior to a single antigen epitope at inducing an immune response (Samad et al., 2022). Multi-epitope vaccines have been used to inhibit a variety of pathogens. For example, multi-epitope vaccines targeting Epstein–Barr virus (EBV) and influenza virus successfully elicit robust neutralizing antibodies and T-cell responses. The EBV polyepitope vaccine, which includes 20 CD8+ T-cell epitopes, induces high frequencies of CD4+ T cells and CD8+ T cells, which are maintained for more than 7 months (Dasari et al., 2023). It has been reported that the norovirus (NoV) P particle, a 24-self-assembly of the protruding domain, displays H16, M2e, and NP9 on its surface, which elicit long-lasting hemagglutinin (HA)-specific neutralizing antibodies, M2e-CD4+ and CD8+ T-cell responses, and a nucleoprotein-specific CTL response, respectively (Nie et al., 2023).
The HEV exists in two states, quasi-enveloped virions, which have a surface membrane and non-enveloped naked virions. The ORF3 and ORF2 antigens are located on the surface of these two types of virions, respectively (Takahashi et al., 2008). Considering the reduction in antibody-mediated neutralization of HEV resulting in antibody evasion by the secreted form of ORF2 (Yin et al., 2018) and the partly protective ORF3 presented on the surface of quasi-enveloped HEV (Takahashi et al., 2008), we first induced ORF3 as an epitopic candidate combined with ORF2 as a multi-peptide vaccine. As previous studies have demonstrated, ORF3 can partially neutralize quasi-enveloped HEV (Ma et al., 2009; Syed et al., 2017; Maurer et al., 2022). The designed vaccine will therefore potentially possess efficacy against both non-enveloped and quasi-enveloped HEV.
In this study, we constructed peptide epitopic vaccines targeting HEV protein ORF2 and ORF3 through computational biology. We screened 9 CTL antigenic (MHC-I) epitopes, 8 HTL (MHC-II) epitopes, and 6 LBL epitopes that covered a wide population of up to 90.60%. On evaluation, these antigens showed good immunogenicity, low allergenicity and non-toxicity. The multi-epitope vaccine is connected through GPGPG, KK, AAY, and EAAAK connectors, and human β-defensin 3 and LL37 are used as adjuvants to enhance the immune response. As the 3D structure validated by Ramachandran plots, the percentage of residues in the most-favored regions was up to 93.2%, and only 1.2% of residues were disallowed. This is in accordance with the good quality model for Ramachandran plots, which is over 90% in the most favored regions (Hooft et al., 1997). Our multi-epitopic vaccine has potential efficacy and low allergenicity against HEV. The Z-score of the 3D structure docked with TLR3 was −2.4 by HADDOCK. The more negative the Z-score, the better the docked structure. This score suggests that the docked structure is of good quality. It has been reported that abundant TLR3 and IFNγ levels were helpful to inhibit HEV and promote recovery (Majumdar et al., 2015; Li et al., 2019), it appears likely that our design will lead to a promising result. Under simulated immunization with the peptide vaccine, both IgM and IgG levels were significantly elevated, as were the abundance of B cells, T cells, and IFN-γ. It suggests the multi-peptide vaccine will induce robust cellular and humoral immune responses effectively. These epitopes may provide a reference for the design of other ORF2 and ORF3 vaccines in the future.
The epitopes that cause excellent immune responses are mainly concentrated on ORF2 and ORF3 (Dawson et al., 1992; Khudyakov et al., 1993; Riddell et al., 2000; Srivastava et al., 2007; Coursaget et al., 2010; Favorov et al., 2010) and the existing HEV epitopic vaccines mainly target ORF2. In these studies, the parameter settings like the selected HLA range and the server threshold settings determine the final peptide selection. The selected peptides all differ between vaccine studies but effective immunoreactions have been induced using these different strategies (Samavia et al., 2022; Anwar et al., 2023; Ikram et al., 2023). The findings of our study are novel in that ORF3 was included as an additional target along with ORF2.
5 Conclusion
Hepatitis E is a self-limiting disease, but presents a serious threat to pregnant women and immunocompromised individuals (Qian et al., 2023). There is no global vaccine currently available against HEV. Using an immunoinformatic approach, we designed a polypeptide vaccine against the ORF2 and ORF3 proteins of HEV, which was capable of inducing strong T- and B-cell responses by simulated immunization. This is the first study to consider ORF3 as a multi-epitope candidate to combat HEV. Our design provides potential efficacy against both non-enveloped and quasi-enveloped HEV, however, the efficacy of this vaccine requires experimental validation.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
QL: Data curation, Formal analysis, Methodology, Software, Writing – original draft, Writing – review & editing. HW: Data curation, Methodology, Software, Writing – review & editing. JM: Formal analysis, Writing – review & editing. JWa: Formal analysis, Methodology, Software, Writing – original draft. JWu: Formal analysis, Writing – review & editing. SL: Formal analysis, Writing – review & editing. JT: Formal analysis, Writing – review & editing. JN: Project administration, Writing – review & editing. WH: Project administration, Writing – review & editing.
Funding
The authors declare that financial support was received for the research, authorship, and/or publication of this article. This research was funded by the National Key Research and Development Program of China (2023YFC2307900).
Conflict of interest
HW was employed by Wuhan Institute of Biological Products Co., Ltd. JWu was employed by Beijing Yunling Biotechnology Co., Ltd.
The remaining 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.
Footnotes
References
Abraham, M. J., Murtola, T., Schulz, R., Páll, S., Smith, J. C., Hess, B., et al. (2015). GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX 2, 19–25. doi: 10.1016/j.softx.2015.06.001
Ahmad, I., Holla, R. P., and Jameel, S. (2011). Molecular virology of hepatitis E virus. Virus Res. 161, 47–58. doi: 10.1016/j.virusres.2011.02.011
Anwar, T., Ismail, S., Parvaiz, F., Abbasi, S. W., Al-Abbasi, F., Alghamdi, A., et al. (2023). Computational design of experimentally validated multi-epitopes vaccine against hepatitis E virus: An immunological approach. PLoS One 18:e0294663. doi: 10.1371/journal.pone.0294663
Arega, A. M., Dhal, A. K., Pattanaik, K. P., Nayak, S., and Mahapatra, R. K. (2023). An immunoinformatics-based study of Mycobacterium tuberculosis region of difference-2 uncharacterized protein (Rv1987) as a potential subunit vaccine candidate for preliminary ex vivo analysis. Appl. Biochem. Biotechnol. doi: 10.1007/s12010-023-04658-9 [Epub ahead of print].
Behmard, E., Soleymani, B., Najafi, A., and Barzegari, E. (2020). Immunoinformatic design of a COVID-19 subunit vaccine using entire structural immunogenic epitopes of SARS-CoV-2. Sci. Rep. 10:20864. doi: 10.1038/s41598-020-77547-4
Closter, J. M., Bjoern, P., Morten, N., and Paolo, M. (2017). BepiPred-2.0: Improving sequence-based B-cell epitope prediction using conformational epitopes. Nucleic Acids Res. 45, W24–W29. doi: 10.1093/nar/gkx346
Coursaget, P., Buisson, Y., Depril, N., Cann, P., Chabaud, M., Molini, C., et al. (2010). Mapping of linear B cell epitopes on open reading frames 2- and 3-encoded proteins of hepatitis E virus using synthetic peptides. FEMS Microbiol. Lett. 109, 251–255. doi: 10.1111/j.1574-6968.1993.tb06176.x
Dalton, H. R., Bendall, R., Ijaz, S., and Banks, M. (2008). Hepatitis E: An emerging infection in developed countries. Lancet Infect. Dis. 8, 698–709. doi: 10.1016/S1473-3099(08)70255-X
Dasari, V., McNeil, L. K., Beckett, K., Solomon, M., Ambalathingal, G., Thuy, T. L., et al. (2023). Lymph node targeted multi-epitope subunit vaccine promotes effective immunity to EBV in HLA-expressing mice. Nat. Commun. 14:4371. doi: 10.1038/s41467-023-39770-1
Dawson, G., Chau, K., Cabal, C., Yarbough, P., Reyes, G., and Mushahwar, I. (1992). Solid-phase enzyme-linked immunosorbent assay for hepatitis E virus IgG and IgM antibodies utilizing recombinant antigens and synthetic peptides. J. Virol. Methods 38, 175–186. doi: 10.1016/0166-0934(92)90180-l
Dimitrov, I., Bangov, I., Flower, D. R., and Doytchinova, I. (2014). AllerTOP v.2–a server for in silico prediction of allergens. J. Mol. Model. 20:2278. doi: 10.1007/s00894-014-2278-5
Ding, Q., Heller, B., Capuccino, J. M., Song, B., Nimgaonkar, I., Hrebikova, G., et al. (2017). Hepatitis E virus ORF3 is a functional ion channel required for release of infectious particles. Proc. Natl. Acad. Sci. U.S.A. 114, 1147–1152. doi: 10.1073/pnas.1614955114
Doytchinova, I., and Flower, D. (2007). VaxiJen: A server for prediction of protective antigens, tumour antigens and subunit vaccines. BMC Bioinformatics 8:4. doi: 10.1186/1471-2105-8-4
Emini, E. A., Hughes, J. V., Perlow, D. S., and Boger, J. (1985). Induction of hepatitis A virus-neutralizing antibody by a virus-specific synthetic peptide. J. Virol. 55, 836–839. doi: 10.1128/JVI.55.3.836-839.1985
Favorov, M. O., Fields, H. A., Purdy, M. A., Yashina, T. L., and Margolis, H. S. (2010). Serologic identification of hepatitis E virus infections in epidemic and endemic settings. J. Med. Virol. 36, 246–250.
Feng, Z., Hirai-Yuki, A., McKnight, K. L., and Lemon, S. M. (2014). Naked viruses that aren’t always naked: Quasi-enveloped agents of acute hepatitis. Annu. Rev. Virol. 1, 539–560. doi: 10.1146/annurev-virology-031413-085359
Geourjon, C., and Deléage, G. (1995). SOPMA: Significant improvements in protein secondary structure prediction by consensus prediction from multiple alignments. Comput. Appl. Biosci. 11, 681–684. doi: 10.1093/bioinformatics/11.6.681
Gupta, S., Kapoor, P., Chaudhary, K., Gautam, A., Kumar, R., and Raghava, G. P. (2013). In silico approach for predicting toxicity of peptides and proteins. PLoS One 8:e73957. doi: 10.1371/journal.pone.0073957
He, M., Wang, M., Huang, Y., Peng, W., Zheng, Z., Xia, N., et al. (2016). The ORF3 protein of genotype 1 hepatitis E virus suppresses TLR3-induced NF-κB signaling via TRADD and RIP1. Sci. Rep. 6:27597. doi: 10.1038/srep27597
Honorato, R. V., Koukos, P. I., Jiménez-García, B., Tsaregorodtsev, A., Verlato, M., Giachetti, A., et al. (2021). Structural biology in the clouds: The WeNMR-EOSC ecosystem. Front. Mol. Biosci. 8:729513. doi: 10.3389/fmolb.2021.729513
Hooft, R. W., Sander, C., and Vriend, G. (1997). Objectively judging the quality of a protein structure from a Ramachandran plot. Comput. Appl. Biosci. 13, 425–430. doi: 10.1093/bioinformatics/13.4.425
Hou, W., Wu, H., Wang, S., Wang, W., Wang, B., and Wang, H. (2023). Designing a multi-epitope vaccine to control porcine epidemic diarrhea virus infection using immunoinformatics approaches. Front. Microbiol. 14:1264612. doi: 10.3389/fmicb.2023.1264612
Ikram, A., Alzahrani, B., Zaheer, T., Sattar, S., Rasheed, S., Aurangzeb, M., et al. (2023). An in silico deep learning approach to multi-epitope vaccine design: A hepatitis e virus case study. Vaccines (Basel) 11:710. doi: 10.3390/vaccines11030710
Ivan, D., Lyudmila, N., Irini, D., and Ivan, B. (2014). AllergenFP: Allergenicity prediction by descriptor fingerprints. Bioinformatics 6, 846–851. doi: 10.1093/bioinformatics/btt619
Junsu, K., Hahnbeom, P., Lim, H., and Chaok, S. (2012). GalaxyWEB server for protein structure prediction and refinement. Nucleic Acids Res. 40:W294.
Khudyakov, Y. E., Khudyakova, N. S., Fields, H. A., Jue, D., Starling, C., Favorov, M. O., et al. (1993). Epitope mapping in proteins of hepatitis E virus. Virology 194, 89–96.
Kolaskar, A. S., and Tongaonkar, P. C. (1990). A semi-empirical method for prediction of antigenic determinants on protein antigens. FEBS Lett. 276, 172–174.
Lande, R., Pietraforte, I., Mennella, A., Palazzo, R., and Frasca, L. (2021). Complementary effects of carbamylated and citrullinated LL37 in autoimmunity and inflammation in systemic lupus erythematosus. Int. J. Mol. Sci. 22:1650. doi: 10.3390/ijms22041650
Lei, Q., Li, L., Zhang, S., Li, T., Zhang, X., Ding, X., et al. (2018). HEV ORF3 downregulates TLR7 to inhibit the generation of type I interferon via impairment of multiple signaling pathways. Sci. Rep. 8:8585. doi: 10.1038/s41598-018-26975-4
Li, S. W., Zhao, Q., Wu, T., Chen, S., Zhang, J., and Xia, N. S. (2015). The development of a recombinant hepatitis E vaccine HEV 239. Hum. Vaccin. Immunother. 11, 908–914. doi: 10.1080/21645515.2015.1008870
Li, Y., Qu, C., Yu, P., Ou, X., Pan, Q., and Wang, W. (2019). The interplay between host innate immunity and hepatitis E virus. Viruses 11:541.
Ma, H., Song, X., Harrison, T. J., Li, R., Huang, G., Zhang, H., et al. (2009). Immunogenicity and efficacy of a bacterially expressed HEV ORF3 peptide, assessed by experimental infection of primates. Arch. Virol. 154, 1641–1648. doi: 10.1007/s00705-009-0496-4
Majid, M., and Andleeb, S. (2019). Designing a multi-epitopic vaccine against the enterotoxigenic Bacteroides fragilis based on immunoinformatics approach. Sci. Rep. 9:19780. doi: 10.1038/s41598-019-55613-w
Majumdar, M., Ratho, R. K., Chawla, Y., and Singh, M. P. (2015). Role of TLR gene expression and cytokine profiling in the immunopathogenesis of viral hepatitis E. J. Clin. Virol. 73, 8–13. doi: 10.1016/j.jcv.2015.09.011
Marion, O., Lhomme, S., Nayrac, M., Dubois, M., Pucelle, M., Requena, M., et al. (2020). Hepatitis E virus replication in human intestinal cells. Gut 69, 901–910. doi: 10.1136/gutjnl-2019-319004
Maurer, L., El Andari, J., Rapti, K., Spreyer, L., Steinmann, E., Grimm, D., et al. (2022). Induction of hepatitis E virus anti-ORF3 antibodies from systemic administration of a muscle-specific adeno-associated virus (AAV) vector. Viruses 14:266. doi: 10.3390/v14020266
Montpellier, C., Wychowski, C., Sayed, I. M., Meunier, J. C., Saliou, J. M., Ankavay, M., et al. (2018). Hepatitis E virus lifecycle and identification of 3 forms of the ORF2 capsid protein. Gastroenterology 154:211-223.e218. doi: 10.1053/j.gastro.2017.09.020
Nagashima, S., Takahashi, M., Kobayashi, T., Nishizawa, T., Nishiyama, T., Okamoto, H., et al. (2017). Characterization of the quasi-enveloped hepatitis E virus particles released by the cellular exosomal pathway. J. Virol. 91:e00822-17. doi: 10.1128/JVI.00822-17
Nair, V. P., Anang, S., Subramani, C., Madhvi, A., Bakshi, K., Srivastava, A., et al. (2016). Endoplasmic reticulum stress induced synthesis of a novel viral factor mediates efficient replication of genotype-1 hepatitis E virus. PLoS Pathog. 12:e1005521. doi: 10.1371/journal.ppat.1005521
Nan, Y., and Zhang, Y. J. (2016). Molecular biology and infection of hepatitis E virus. Front. Microbiol. 7:1419. doi: 10.3389/fmicb.2016.01419
Nicolas, R., Ole, L., Massimo, B., Filippo, C., and Vladimir, B. (2010). Computational immunology meets bioinformatics: The use of prediction tools for molecular binding in the simulation of the immune system. PLoS One 5:e9862. doi: 10.1371/journal.pone.0009862
Nie, J., Wang, Q., Jin, S., Yao, X., Xu, L., Chang, Y., et al. (2023). Self-assembled multiepitope nanovaccine based on NoV P particles induces effective and lasting protection against H3N2 influenza virus. Nano Res. 16, 7337–7346. doi: 10.1007/s12274-023-5395-6
Oli, A. N., Obialor, W. O., Ifeanyichukwu, M. O., Odimegwu, D. C., Okoyeh, J. N., Emechebe, G. O., et al. (2020). Immunoinformatics and vaccine development: An overview. Immunotargets Ther. 9, 13–30.
Qi, Y., Zhang, F., Zhang, L., Harrison, T. J., Huang, W., Zhao, C., et al. (2015). Hepatitis E virus produced from cell culture has a lipid envelope. PLoS One 10:e0132503. doi: 10.1371/journal.pone.0132503
Qian, Z., Cong, C., Li, Y., Bi, Y., He, Q., Li, T., et al. (2023). Quantification of host proteomic responses to genotype 4 hepatitis E virus replication facilitated by pregnancy serum. Virol. J. 20:111. doi: 10.1186/s12985-023-02080-5
Reynisson, B., Barra, C., Kaabinejadian, S., Hildebrand, W. H., Peters, B., and Nielsen, M. (2020). Improved prediction of MHC II antigen presentation through integration and motif deconvolution of mass spectrometry MHC eluted ligand data. J. Proteome Res. 19, 2304–2315. doi: 10.1021/acs.jproteome.9b00874
Riddell, M. A., Li, F., and Anderson, D. A. (2000). Identification of immunodominant and conformational epitopes in the capsid protein of hepatitis e virus by using monoclonal antibodies. J. Virol. 74:8011. doi: 10.1128/jvi.74.17.8011-8017.2000
Samad, A., Ahammad, F., Nain, Z., Alam, R., Imon, R. R., Hasan, M., et al. (2022). Designing a multi-epitope vaccine against SARS-CoV-2: An immunoinformatics approach. J. Biomol. Struct. Dyn. 40, 14–30. doi: 10.1080/07391102.2020.1792347
Samavia, N., Fahed, P., Yasir, W., Tasneem, A., and Syeda, N. (2022). Prediction of promiscuous epitopes in ORF2 of hepatitis E virus: An in-silico approach. Afr. Health Sci. 22, 626–639. doi: 10.4314/ahs.v22i3.67
Smith, D. B., Simmonds, P., et al. Members, Of The, International Committee, On The (2014). Consensus proposals for classification of the family Hepeviridae. J. Gen. Virol. 95, 2223–2232.
Songtanin, B., Molehin, A. J., Brittan, K., Manatsathit, W., and Nugent, K. (2023). Hepatitis E virus infections: Epidemiology, genetic diversity, and clinical considerations. Viruses 15:1389.
Srivastava, R., Aggarwal, R., Jameel, S., Puri, P., and Naik, S. (2007). Cellular immune responses in acute hepatitis E virus infection to the viral open reading frame 2 protein. Viral Immunol. 20, 56–65. doi: 10.1089/vim.2006.0053
Syed, S. F., Sun, Y., Du, T., Chen, Y., Liu, B., Wang, X., et al. (2017). Evaluation of recombinant Chinese avian hepatitis E virus (CaHEV) ORF2 and ORF3 proteins for protection of chickens against CaHEV infection. Vaccine 35, 3482–3489. doi: 10.1016/j.vaccine.2017.05.030
Takahashi, M., Tanaka, T., Takahashi, H., Hoshino, Y., Nagashima, S., and Jirintai, et al. (2010). Hepatitis E Virus (HEV) strains in serum samples can replicate efficiently in cultured cells despite the coexistence of HEV antibodies: Characterization of HEV virions in blood circulation. J. Clin. Microbiol. 48, 1112–1125. doi: 10.1128/jcm.02002-09
Takahashi, M., Yamada, K., Hoshino, Y., Takahashi, H., Ichiyama, K., Tanaka, T., et al. (2008). Monoclonal antibodies raised against the ORF3 protein of hepatitis E virus (HEV) can capture HEV particles in culture supernatant and serum but not those in feces. Arch. Virol. 153, 1703–1713. doi: 10.1007/s00705-008-0179-6
Tîrziu, A., Avram, S., Madă, L., Crişan-Vida, M., Popovici, C., Popovici, D., et al. (2023). Design of a synthetic long peptide vaccine targeting HPV-16 and -18 using immunoinformatic methods. Pharmaceutics 15:1798. doi: 10.3390/pharmaceutics15071798
Tyagi, S., Korkaya, H., Zafrullah, M., Jameel, S., and Lal, S. K. (2002). The phosphorylated form of the ORF3 protein of hepatitis E virus interacts with its non-glycosylated form of the major capsid protein, ORF2. J. Biol. Chem. 277, 22759–22767. doi: 10.1074/jbc.M200185200
van Zundert, G., Rodrigues, J., Trellet, M., Schmitz, C., Kastritis, P., Karaca, E., et al. (2016). The HADDOCK2.2 web server: User-friendly integrative modeling of biomolecular complexes. J. Mol. Biol. 428, 720–725. doi: 10.1016/j.jmb.2015.09.014
Wang, W., Xia, M., Chen, J., Deng, F., Yuan, R., Zhang, X., et al. (2016). Data set for phylogenetic tree and RAMPAGE Ramachandran plot analysis of SODs in Gossypium raimondii and G. arboreum. Data Brief 9, 345–348. doi: 10.1016/j.dib.2016.05.025
Wiederstein, M., and Sippl, M. J. (2007). ProSA-web: Interactive web service for the recognition of errors in three-dimensional structures of proteins. Nucleic Acids Res. 35, W407–W410. doi: 10.1093/nar/gkm290
Yadav, K. K., Boley, P. A., Fritts, Z., and Kenney, S. P. (2021). Ectopic expression of genotype 1 hepatitis E virus ORF4 increases genotype 3 HEV viral replication in cell culture. Viruses 13:75. doi: 10.3390/v13010075
Yamada, K., Takahashi, M., Hoshino, Y., Takahashi, H., Ichiyama, K., Nagashima, S., et al. (2009). ORF3 protein of hepatitis E virus is essential for virion release from infected cells. J. Gen. Virol. 90, 1880–1891. doi: 10.1099/vir.0.010561-0
Yin, X., Ying, D., Lhomme, S., Tang, Z., Walker, C. M., Xia, N., et al. (2018). Origin, antigenicity, and function of a secreted form of ORF2 in hepatitis E virus infection. Proc. Natl. Acad. Sci. U.S.A. 115, 4773–4778. doi: 10.1073/pnas.1721345115
Zhang, Z. J., Zhang, Z. X., Huang, H. S., Wu, T., Hu, H. Y., Wang, W. Z., et al. (2015). Long-term efficacy of a hepatitis E vaccine. N. Engl. J. Med. 372, 914–922.
Keywords: multi-epitope, vaccine, HEV, ORF2, ORF3
Citation: Lu Q, Wu H, Meng J, Wang J, Wu J, Liu S, Tong J, Nie J and Huang W (2024) Multi-epitope vaccine design for hepatitis E virus based on protein ORF2 and ORF3. Front. Microbiol. 15:1372069. doi: 10.3389/fmicb.2024.1372069
Received: 17 January 2024; Accepted: 07 March 2024;
Published: 21 March 2024.
Edited by:
Alexandro Guterres, Oswaldo Cruz Foundation (Fiocruz), BrazilReviewed by:
Kush Kumar Yadav, The Ohio State University, United StatesZoya Shafat, Institute of Management Technology-Dubai, United Arab Emirates
Copyright © 2024 Lu, Wu, Meng, Wang, Wu, Liu, Tong, Nie and Huang. 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: Jianhui Nie, niejianhui@nifdc.org.cn; Weijin Huang, huangweijin@nifdc.org.cn