1. Introduction {#s0005}
===============

Dengue virus or as commonly called DENV is a single stranded RNA virus that infects approximately 390 million people each year ([@bb0005]), putting more than two-fifth of the world\'s population under the threat of this efficacious virus. The dengue fever has, thus become one of the most widespread disease ([@bb0300]). The virus belongs to the family *Flaviviridae* and genus *Flavivirus* ([@bb0330]; [@bb0025]). DENV is an arbovirus, having two known mosquito vectors *Aedes aegypti* ([@bb0110]) and *Aedes albopictus* ([@bb0185]). The positive stranded RNA genome of dengue virus is of 10.7 Kb size and composed of three structural proteins (Envelope, Capsid, Membrane) and seven non-structural proteins (NS1, NS2A, NS2B,NS3, NS4A, NS4B, NS5) ([@bb0150]; [@bb0125]; [@bb0300]). There are atleast four serotypes and they show 65% similarity in the genome structure ([@bb0020]; [@bb0270]).

The dengue infection is caused by one of the four serotypes of DENV that are spread by *Aedes* mosquito ([@bb0155]). During primary infection, the body develops immune responses in the form of antibodies against the particular serotype attacked ([@bb0285]). But the main complexity of DENV arises during the secondary infection with another serotype, leading to serious version of dengue infection like Dengue Haemorrhagic fever (DHF) and Dengue Shock Syndrome (DSS) ([@bb0055]; [@bb0210]; [@bb0285]). This is caused due to the antibodies produced during primary attack which complicate the secondary DENV infection by a phenomenon known as Antibody Dependent Enhancement (ADE) ([@bb0075]; [@bb0105]). During ADE, there is a cross reaction between the antibodies of the primary infection and virus of secondary infection such that there is an increased infection in macrophages and monocytes ([@bb0335]; [@bb0080]). These challenges bring the importance of an archetypal dengue vaccine which can provide life time immunity against all the serotypes ([@bb0310]; [@bb0280]).

Currently, the vaccine candidates that are under various stages of clinical trial are the live attenuated viruses ([@bb0315]), chimeric vaccine ([@bb0120]), recombinant vaccine with adjuvants ([@bb0145]), reverse vaccinology ([@bb0035]), purified and inactivated virions ([@bb0095]), subunit proteins and plasmid DNA ([@bb0305]). Among these, live attenuated DENGVAXIA or CYD-TDV, a tetravalent chimeric dengue vaccine ([@bb0290]), developed by Sanofi Pasteur in December 2015, is the first licensed vaccine in some Asian and Latin American countries ([@bb0260]; [@bb0100]). These clinical manifestations caused by the vaccine are ascribed to inefficiency of the vaccine in producing competent T- cells that protect against DENV disease ([@bb0170]). Moreover, the vaccine does not encode any non- structural proteins which are required by the virus to evade immune response of the host ([@bb0130]; [@bb0230]). All these studies imply that a vaccine that is tetravalent and simultaneously prevents antibody- dependent enhancement (ADE) needs to be designed urgently. These concerns led to the need for a relatively new technique of vaccine development i.e. Epitope or synthetic peptide based vaccines. As DENV has both structural and non-structural proteins for its viral activity ([@bb0245]), conserved epitopes may prove to be useful in designing synthetic peptide based vaccine. This can be easily initiated in today\'s time, as there is no dearth of information about genome sequences in the databases ([@bb0140]; [@bb0295]).

Our study provides answers to all above mentioned concerns as we could successfully find conserved epitopes in structural as well as non- structural proteins. These conserved epitopes can be used for further analysis and exploited to develop the vaccine against the deadly dengue virus as well as in rapid diagnostics.

2. Materials and methods {#s0010}
========================

2.1. Retrieval of sequences and alignment {#s0015}
-----------------------------------------

A sum total of 23,622 partial sequences of structural and non-structural proteins of dengue virus were retrieved from NCBI ([www.ncbi.nih.gov](http://www.ncbi.nih.gov){#ir0005}) ([Table 1](#t0005){ref-type="table"} ). After retrieving, the sequences were aligned using multiple sequence alignment program CLUSTAL_X ([@bb0320]). These aligned sequence files were then used for further analysis.Table 1Number of sequences for each protein of different serotypes of DENV analyzed for this study.Table 1ProteinsDENV-1DENV-2DENV-3DENV-4Envelope (E)500500500500Membrane (Prm)100010001000250Capsid (C)10001000857285NS110001000797201NS2A100100250100NS2B1000100500100NS31000500500100NS4A50080010001000NS4B1000500500100NS51000500882100

2.2. Search for conserved sequences {#s0020}
-----------------------------------

Sequences were used for detecting the species-specific signature sequences or **motifs**. Motifs were obtained for each serotype of DENV individually using an online tool multiple em for motif elicitation or MEME *Suite* ([@bb0030]). Motifs common for all serotypes were also obtained using the same method. In order to get a maximum number of motifs, the default setting was adjusted from 3 motifs to 10 motifs.

2.3. Linear B-cell epitope prediction {#s0025}
-------------------------------------

The motifs were then analyzed for the presence of B cell epitopes. The linear B cell epitopes were found using BCPRED and BEPIPRED tools of immune epitope database with default settings. Bepipred used Hidden Markov model for the prediction of B cell epitopes ([@bb0190]; [@bb0265]).

2.4. IEDB analysis {#s0030}
------------------

The immunogenicity of each epitope was checked using Kolaskar and Tangaonkar antigenicity method ([@bb0175]) with a default threshold value 0.9. Hydrophilicity of the antigenic epitopes, required to check the accessibility were found using Parker Hydrophilicity method ([@bb0250]) at a threshold value of 3.448 ([@bb0295]). Epitopes were checked for their surface accessibility using Emini surface accessibility method ([@bb0090]) with a threshold value of 1.00. Flexibility and Beta turns were checked using Karplus and Schulz Flexibility ([@bb0160]) and Chou and Fasman Beta-turn methods ([@bb0050]) respectively, with a threshold of 1.00 for both ([@bb0295]). Conservancy of epitopes was checked using Epitope Conservancy Analysis tool ([@bb0045]).

2.5. 3-D structure prediction {#s0035}
-----------------------------

Models of the 3D structures of DENV proteins were downloaded from RCSB PDB server (<https://www.rcsb.org>/) for mapping the epitopes. Then the location of predicted epitopes in the 3-D model was found using CHIMERA visualization tool ([@bb0255]). For DENV proteins where no 3D strcture was available, I-TASSER server was used to predict the 3-D structures ([@bb0350]; [@bb0275]). The predicted models were saved in pdb format files, which were later used to generate Ramachandran plot using PDBsum --PROCHECK software ([@bb0195]). This was done to verify the models generated by I-TASSER.

3. Results {#s0040}
==========

In the present study, comparative genomics was performed on DENV sequences to find novel vaccine targets. For this purpose, all the ten proteins (3 structural and 7 non-structural) of DENV-1, DENV-2, DENV-3, DENV-4 were analyzed individually by retrieving their partial protein sequences and motif analysis was performed by using MEME SUITE. These conserved motif sequences were then used for B-cell epitope analysis ([Table 2](#t0010){ref-type="table"} ). For the matter of space and clarity, we are reporting the figures for one of the conserved epitopes in the main text, while the others are presented in the supplementary materials ([Fig. 1](#f0005){ref-type="fig"}, [Fig. 2](#f0010){ref-type="fig"}, [Fig. 3](#f0015){ref-type="fig"} , Supplementary data).Table 2Epitopes found within the conserved motif region.Table 2[^1]Fig. 1Motif 'YQNKVVRVQRPAKNGTVMDVISSRDQRGSGQVGTYGLNTFTNMEAQLIRQ' containing the epitope 'QRGSGQV' as 100% conserved region in all sequences of protein NS5.Fig. 1Fig. 2(A) Linear B-cell epitope prediction, (B) Antigenicity prediction, (C) Surface accessibility prediction, (D) Hydrophilicity prediction, (E) flexibility prediction graphs of epitope 'QRGSGQV' of NS5. The graphs were plotted between the score for predicted epitope and the position of epitope in protein sequence. The red line represents the threshold value used for epitope prediction. The yellow region indicates the possible region of B-cell epitope in the protein sequence. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.)Fig. 2Fig. 3Three-dimensional visualization of predicted epitope 'QRGSGQV' of NS5 protein. PDB protein structure of dengue virus NS5 protein (4V0Q) was used to map 100% conserved epitope 'QRGSGQV' using the Chimera Visualization tool and the epitope was found to be present on beta turn and surface accessible.Fig. 3

In total 13 epitopes were found which were further checked for antigenicity, hydrophilicity and surface accessibility. The B-cell epitopes on native proteins are generally composed of hydrophillic amino acids on the protein surface that are topographically accessible to membrane-bound or free antibody. These epitopes tend to be located in flexible regions of an immunogen and display site mobility of epitopes which maximizes complementarity with the antibody\'s binding site, permitting an antibody to bind with an epitope that it might bind ineffectively if it were rigid. Surface accessibility and hydrophilicity of these predicted epitopes was therefore determined ([Table 3](#t0015){ref-type="table"} ). The flexibility of these epitopes and the presence of beta turn in their structure were also evaluated. Predicted epitopes of Envelope, NS1, NS3, NS4A, NS4B and NS5 proteins were found out to be antigenic as their score was found out to be higher than the threshold value whereas, few antigens were found to be non-antigenic ([Table 3](#t0015){ref-type="table"}). All the predicted epitopes were found to be flexible in nature ([Table 3](#t0015){ref-type="table"}).Table 3B-cell epitope analysis.Table 3ProteinPeptide sequenceAntigenicity prediction scoreHydrophillicity prediction scoreSurface accessibility prediction scoreFlexibility prediction scoreBeta turn prediction scoreInferenceConservancyTHRSHOLD SCORE→0.93.4481.0001.0001.000ETPQAPTS1.005\
Antigenic4.171\
Hydrophillic2.599\
Surface accessible1.0611.147EpitopeConserved only in D1RCPTQGE1.00\
Antigenic4.629\
Hydrophillic1.0\
Surface accessible1.0851.129Epitope100% ConservedCapsidIPPTAGIL1.066\
Antigenic−1.0\
Non-hydrophillic0.408\
Surface accessible1.0180.969Non-epitopeIPPTAGVL in D3PrMEHRRDKRS0.923\
Antigenic5.588\
Hydrophillic7.594\
Surface accessible1.0351.055EpitopeEHRREKRS in D2\
ERRREKRS in D4NS1GYHTQTA1.005\
Antigenic3.486\
Hydrophillic1.824\
Surface accessible1.0381.03EpitopeGYHTQIA in D2\
GYFTQTA in D1\
GYATQTA in D4PNTPECPNA1.02\
Antigenic4.657\
Hydrophillic1.033\
Surface accessible1.0531.293EpitopePNTPECPSA in D3\
PDTSECPNE in D4\
PETAECPNT in D2NS3SAAQRRGR0.968\
Antigenic4.400\
Hydrophillic2.371\
Surface accessible1.0931.156Epitope100% conservedPGTSGSPI0.973\
Antigenic5.343\
Hydrophillic1.467\
Surface accessible1.1301.426Epitope100% conservedNS4AQRTPQDNQL0.992\
Antigenic5.643\
Hydrophillic3.674\
Surface accessible1.0691.080Epitope100% conservedNS4BLQAKATREAQKRA1.037\
Antigenic5.029\
Hydrophillic4.181\
Surface accessible1.0441.126Epitope100% conservedNS5WHYDQDHPYKT1.034\
Antigenic4.857\
Hydrophillic3.436\
Surface accessible1.0261.236Epitope100% conserved in D2 onlySHWVPTSRTT1.008\
Antigenic3.529\
Hydrophillic2.635\
Surface accessible1.0811.04EpitopeSHWFPTSRTT in D4QRGSGQV1.007\
Antigenic4.343\
Hydrophillic1.395\
Surface accessible1.1131.137Epitope100% conserved

3.1. Conserved epitope in RNA-dependent RNA polymerase (RdRp) of DENV {#s0045}
---------------------------------------------------------------------

Among all the proteins of dengue virus, RNA-dependent RNA polymerase (RdRp) is encoded by non-structural gene NS5 and consists of 20 important enzymatic activities needed for viral propagation. Studies have been conducted on RdRp for designing drug against dengue virus as it is known to be the largest and highly conserved protein ([@bb0235]; [@bb0240]; [@bb0085]; [@bb0200]). Conserved motif 'YQNKVVRVQRPAKNGTVMDVISSRDQRGSGQVGTYGLNTFTNMEAQLIRQ' was found in NS5 ([Fig. 1](#f0005){ref-type="fig"}), which when inspected for presence of potential epitope resulted in 'QRGSGQV' ([Table 2](#t0010){ref-type="table"}). Linear B-cell epitope prediction, Antigenicity prediction, Surface accessibility prediction, Hydrophilicity prediction and flexibility prediction of this region ([Fig. 2](#f0010){ref-type="fig"}A--E) showed 'QRGSGQV' as the most desirable epitope. The epitope was further checked for conservancy using IEDB Conservancy analysis and was found to be 100% conserved in all serotypes of DENV ([Table 3](#t0015){ref-type="table"}), thus making it one of the most efficient epitope for vaccine development based on in silico analysis. PDB protein structure of dengue virus NS5 protein (4V0Q) ([@bb0355]) was used to map this epitope and the epitope was found to be located in beta turn region of the RNA-directed RNA polymerase of DENV structure ([Fig. 3](#f0015){ref-type="fig"}), thus making it accessible for interaction.

3.2. Highly conserved epitope 'RCPTQGE' on structural protein - Envelope {#s0050}
------------------------------------------------------------------------

As stated above, there are three structural and 7 non-structural proteins in DENV polyprotein. Envelope protein is highly variable and is notably displayed on the surface of DENV ([@bb0220]), thus making it the crucial candidate for vaccine development ([@bb0165]). In our study, we could find two epitopes in Envelope protein ([Table 2](#t0010){ref-type="table"}), but among them, epitope 'RCPTQGE' was found not only antigenic and surface accessible, but was also 100% conserved in all serotypes of DENV ([Table 3](#t0015){ref-type="table"}, Supplementary Figs. S1--S3). On mapping this epitope on 3D structure of Envelope protein CryoEM structure of Dengue virus envelope protein heterotetramer (pdb id [3J2P](pdb:3J2P){#ir0015}), it was found to be present in beta turn region and therefore, is accessible for interaction, making it yet another highly suitable candidate for vaccine target. As the structural proteins are highly variable, therefore it was difficult to find any conserved epitope for capsid and PrM proteins.

3.3. Other highly conserved epitopes of DENV serotypes {#s0055}
------------------------------------------------------

To our surprise, four more predicted epitopes (other than the conserved epitopes of NS5 and envelope proteins) were found to be 100% conserved i.e., NS3 (**SAAQRRGR**, **PGTSGSPI**); NS4A (**QRTPQDNQL**); NS4B (**LQAKATREAQKRA**) and in all DENV serotypes ([Table 3](#t0015){ref-type="table"}).

The three dimensional visualization of predicted epitopes of NS4A and NS4B protein was performed using I-TASSER server (Supplementary Figs. S12, S15) as the pdb structure was not available for these proteins. The predicted structure of protein sequences was validated by plotting Ramachandran graph using PROCheck software (data not shown). Remaining epitopes were found to be serotype-specific ([Table 3](#t0015){ref-type="table"}) and hence 3D structure prediction and validation was not performed for these epitopes.

3.4. Serotype-specific conserved epitopes of DENV {#s0060}
-------------------------------------------------

For developing a successful vaccine candidate, epitope common in all the serotypes is a perquisite. But for developing serotype-specific drugs and in diagnostics of serotype specific strains, serotype-specific conserved epitopes may prove to be useful. Hence, epitopes like 'TPQAPTS' of Envelope protein (specific for DENV-1 serotype); Capsid protein epitopes 'EHRREKRS' for DENV-2 and 'ERRREKRS' for DENV-4; NS1 epitope 'GYH/F/ATQT/IA' and 'PN/D/ETP/S/AECPN/SA/E/T'; NS5 epitope 'WHYDQDHPYKT' in DENV-2 and 'SHWV/FPTSRTT' epitope may be exploited for developing diagnostic kits for early detection. Among the above reported serotype specific epitopes, both the epitopes of NS1 are potential candidates for distinguishing one serotype from other.

4. Discussion {#s0065}
=============

Dengue is one of the most rapidly spreading mosquito-borne viral disease ([@bb0340]) and has emerged as one of the biggest threats to public health ([@bb0115]; [@bb0305]; [@bb0010]). There are four serotypes of DENV which show 65% similarity at genomic level and share same epidemics. Yet the difference in their interactions with antibodies is enough to make the development of a common vaccine, a mammoth task ([@bb0180]). The urgent need of dengue vaccine development is presented by the alarming rise in the dengue endemics ([@bb0065]; [@bb0345]). Although various vaccines are undergoing clinical trials worldwide but presently, there is no specific dengue therapeutics or vaccine available which provides overall protection against this viral infection. Recently, the only licensed DENV vaccine Dengvaxia, has been reported as a failure as it lacks DENV non-structural protein antigens, and is found to be incompetent to raise antibodies against NS1 ([@bb0130]).

In the present study, Comparative Genomics and immunoinformatics was used as a tool to explore the potential candidates of dengue virus to find novel drug and vaccine targets. The advantage of development of an epitope based vaccine over other vaccines is its ability to induce specific immune response without undesirable effects. Also, both time and expenditure needed to screen a large number of epitopes can be saved using such computational approaches ([@bb0060]). Similar studies have been conducted on Human *Coronavirus* ([@bb0295]), *Saint Louis encephalitis virus* ([@bb0140]), Rotaviruses ([@bb0225]), H1N1 influenza A virus strains ([@bb0040]) and Zika virus ([@bb0070]), however, the analysis has been conducted for only some selective proteins. Therefore our study is the first one to deal with all the proteins of any viral genome for designing an efficient vaccine.

To achieve our objective, conserved epitopes have been obtained from both structural and non-structural proteins. We could successfully found six highly conserved epitopes, one from structural protein and other five from non-structural proteins. Among structural proteins, only Envelope protein showed 100% conserved epitope (**RCPTQGE**), indicating this as a good target for vaccine designing. Other structural proteins like PrM and Capsid did not show any conservancy among the predicted epitopes. This was expected as the structural proteins are known to be highly variable even within serotypes ([@bb0220]). The predicted epitopes of Envelope protein (**RCPTQGE**); NS3 (**SAAQRRGR**, **PGTSGSPI**); NS4A (**QRTPQDNQL**); NS4B (**LQAKATREAQKRA**) and NS5 proteins (**QRGSGQV**) were found out to be antigenic, hydrophillic, surface accessible, flexible and consists of beta turn. Thus, these epitopes fulfill all the pre-requisite conditions of a desired epitope. Moreover, the functions of these proteins can be used as an asset for controlling the effect of the virus on the host. For instance, NS1 is an important cofactor for the formation of replication complex, and has also been exploited for the development of vaccine ([@bb0205]; [@bb0145]) whereas NS4B regulates the helicase activity of NS3 and NS4A-NS4B complex is required for genome replication ([@bb0215]; [@bb0325]). NS5 is the most conserved protein of DENV and its C-terminal containing RdRp (RNA dependent RNA polymerase) domain is crucial for genomic replication ([@bb0135]: [@bb0015]). Therefore, we can say that each of these proteins are essential for replication of the virus inside the host and therefore, they can be used as common vaccine target for all the four serotypes of Dengue virus. Serotype specific conserved epitopes can also be exploited for detection of Dengue at early stages. During the Phase IIb and III trials, DENGVAXIA\'s specific response against DENV2 serotype largely became negligible. So, a vaccine developed using these epitopes could overcome this problem as these predicted epitopes are found to be conserved in all four serotypes of DENV.

5. Conclusion {#s0070}
=============

Dengue Virus (DENV) has emerged as a potential threat to human health worldwide. The therapeutics against DENV are either not available in major parts of the world or if available in some endemic countries, are found to be inefficient. In this study we used computational approaches to find novel vaccine targets. It focuses to explore B-cell epitopes for all serotypes for each protein of dengue virus. We found six highly conserved epitopes in all serotypes of DENV \[Envelope protein (**RCPTQGE**); NS3 (**SAAQRRGR, PGTSGSPI**); NS4A (**QRTPQDNQL**); NS4B (**LQAKATREAQKRA**) and NS5 proteins (**QRGSGQV**)\]. Thus, our results suggest any of these proteins can be targeted to stimulate a specific immune response against all serotypes of DENV. These predicted epitopes would be the candidate target for the universal multi-subunit vaccine.

Nevertheless, further studies are needed to confirm the utility of these epitopes.

MV and RL conceptualized the work; MV, SB, KK, NM, SS,SG curated data; formal analysis was done by MV, PSD; MV, PSD were involved in fund acquisition and project administration; MV supervised the work; MV. SB, KK were involved in writing original draft; Review and editing was done by MV, PSD, RL.

Funding {#s0075}
=======

The funding to conduct this work was provided by Delhi University Innovative Scheme.

Ethical approval {#s0080}
================

Not applicable.

Conflict of interest {#s0085}
====================

The authors declare that they have no conflict of interest.

Authors\' information {#s0090}
=====================

This work is performed by the undergraduate students at Sri Venkateswara College, University of Delhi.

Declaration of interests {#s0095}
========================

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Appendix A. Supplementary data {#s0100}
==============================

Supplementary figuresImage 1

Authors acknowledge Dr. P. Hemalatha Reddy, Principal, Sri Venkateswara College, University of Delhi for constantly encouraging Undergraduate students for research.

Supplementary data to this article can be found online at <https://doi.org/10.1016/j.gene.2019.02.001>.

[^1]: Epitopes (marked in blue) found within the conserved motif region.
