Using in-silica Analysis and Reverse Vaccinology Approach for COVID-19 Vaccine Development

Ajay Agarwal


Background: The recent pandemic of COVID19 that has struck the world is yet to be battled by a potential cure. Countless lives have been claimed due to the existing pandemic and the societal normalcy has been damaged permanently. As a result, it becomes crucial for academic researchers in the field of bioinformatics to combat the existing pandemic. Materials and Methods: The study involved collecting the virulent strain sequence of SARS-nCoV19 for the country USA against human host through publically available bioinformatics databases. Using in-silica analysis, reverse vaccinology, and 3-D modelling, two leader proteins were identified to be potential vaccine candidates for development of a multi-epitope drug. Results: It was revealed that the two leader proteins ORF1ab MT326102 and MT326715 had the highest extinction coefficient and the lowest score on the GRAVY. Along with the given parameters, these leader proteins were highly stable and were also antigenic in nature. The two selected epitopes were then docked against their respective alleles to obtain the global energy scores, which was the lowest of all possible pairs. Conclusion: The epitopes which displayed the lowest global energy score on docking with the alleles were selected and proposed as successful and potential vaccine candidates for multi-epitope vaccine development.


Doi: 10.28991/SciMedJ-2020-02-SI-9

Full Text: PDF


In-silica Analysis; SARS-CoV2; COVID-19; Vaccine; Bioinformatics; Reverse Vaccinology; Strain.


Q&A on coronaviruses (COVID-19), Available online: (accessed on 17 April 2020).

Coronavirus disease (COVID-19) pandemic, Available online: (accessed on 31 March 2020).

R&D Blueprint and COVID-19, Available online: (accessed on 30 March 2020).

Pickett, B. E., Sadat, E. L., Zhang, Y., Noronha, J. M., Squires, R. B., Hunt, V., … Scheuermann, R. H. (2011). ViPR: an open bioinformatics database and analysis resource for virology research. Nucleic Acids Research, 40(D1), D593–D598. doi:10.1093/nar/gkr859.

Charif, D., & Lobry, J. R. (2007). SeqinR 1.0-2: A Contributed Package to the R Project for Statistical Computing Devoted to Biological Sequences Retrieval and Analysis. Biological and Medical Physics, Biomedical Engineering, 207–232. doi:10.1007/978-3-540-35306-5_10.

Pagès, H., Aboyoun, P., Gentleman, R., & DebRoy, S. (2017). Biostrings: Efficient manipulation of biological strings. R package version, 2(0).

Osorio, D., Rondón-Villarrea, P. and Torres, R., 2015. Peptides: a package for data mining of antimicrobial peptides. R Journal, 7(1).

Doytchinova, I. A., & Flower, D. R. (2007). VaxiJen: a server for prediction of protective antigens, tumour antigens and subunit vaccines. BMC Bioinformatics, 8(1). doi:10.1186/1471-2105-8-4.

Meunier, M., Guyard-Nicodème, M., Hirchaud, E., Parra, A., Chemaly, M., & Dory, D. (2016). Identification of Novel Vaccine Candidates againstCampylobacterthrough Reverse Vaccinology. Journal of Immunology Research, 2016, 1–9. doi:10.1155/2016/5715790.

Vita, R., Overton, J.A., Greenbaum, J.A., Ponomarenko, J., Clark, J.D., Cantrell, J.R., Wheeler, D.K., Gabbard, J.L., Hix, D., Sette, A. and Peters, B., (2015). The immune epitope database (IEDB) 3.0. Nucleic acids research, 43(D1), pp. D405-D412.

Ullah, M. A., Sarkar, B., & Islam, S. S. (2020). Exploiting the reverse vaccinology approach to design novel subunit vaccines against Ebola virus. Immunobiology, 225(3), 151949. doi:10.1016/j.imbio.2020.151949.

Thevenet, P., Shen, Y., Maupetit, J., Guyon, F., Derreumaux, P., & Tuffery, P. (2012). PEP-FOLD: an updated de novo structure prediction server for both linear and disulfide bonded cyclic peptides. Nucleic Acids Research, 40(W1), W288–W293. doi:10.1093/nar/gks419.

Shen, Y., Maupetit, J., Derreumaux, P., & Tufféry, P. (2014). Improved PEP-FOLD Approach for Peptide and Miniprotein Structure Prediction. Journal of Chemical Theory and Computation, 10(10), 4745–4758. doi:10.1021/ct500592m.

Lamiable, A., Thévenet, P., Rey, J., Vavrusa, M., Derreumaux, P., & Tufféry, P. (2016). PEP-FOLD3: fasterde novostructure prediction for linear peptides in solution and in complex. Nucleic Acids Research, 44(W1), W449–W454. doi:10.1093/nar/gkw329.

Zhou, P., Jin, B., Li, H., & Huang, S.-Y. (2018). HPEPDOCK: a web server for blind peptide–protein docking based on a hierarchical algorithm. Nucleic Acids Research, 46(W1), W443–W450. doi:10.1093/nar/gky357.

Full Text: PDF

DOI: 10.28991/SciMedJ-2020-02-SI-9


  • There are currently no refbacks.

Copyright (c) 2020 Ajay Agarwal