Potential therapeutic target identification in the novel 2019 coronavirus: insight from homology modeling and blind docking study

Olanrewaju Ayodeji Durojaye1, Talifhani Mushiana2, Henrietta Onyinye Uzoeto3, Samuel Cosmas4, Victor Malachy Udowo5, Abayomi Gaius Osotuyi6, Glory Omini Ibiang3, Miapeh Kous Gonlepa7
1School of Life Sciences, Department of Molecular and Cell Biology, University of Science and Technology of China, Hefei, China
2School of Chemistry and Material Sciences, Department of Chemistry, University of Science and Technology of China, Hefei, China
3Department of Biological Sciences, Coal City University, Emene, Enugu State, Nigeria
4Department of Biochemistry, University of Nigeria, Nsukka, Enugu State, Nigeria
5Institute of Metal Research, Chinese Academy of Sciences, Shenyang, China
6School of Earth and Space Sciences, University of Science and Technology of China, Hefei, China
7School of Public Affairs, Department of Public Administration, University of Science and Technology of China, Hefei, China

Tóm tắt

Abstract Background The 2019-nCoV which is regarded as a novel coronavirus is a positive-sense single-stranded RNA virus. It is infectious to humans and is the cause of the ongoing coronavirus outbreak which has elicited an emergency in public health and a call for immediate international concern has been linked to it. The coronavirus main proteinase which is also known as the 3C-like protease (3CLpro) is a very important protein in all coronaviruses for the role it plays in the replication of the virus and the proteolytic processing of the viral polyproteins. The resultant cytotoxic effect which is a product of consistent viral replication and proteolytic processing of polyproteins can be greatly reduced through the inhibition of the viral main proteinase activities. This makes the 3C-like protease of the coronavirus a potential and promising target for therapeutic agents against the viral infection. Results This study describes the detailed computational process by which the 2019-nCoV main proteinase coding sequence was mapped out from the viral full genome, translated and the resultant amino acid sequence used in modeling the protein 3D structure. Comparative physiochemical studies were carried out on the resultant target protein and its template while selected HIV protease inhibitors were docked against the protein binding sites which contained no co-crystallized ligand. Conclusion In line with results from this study which has shown great consistency with other scientific findings on coronaviruses, we recommend the administration of the selected HIV protease inhibitors as first-line therapeutic agents for the treatment of the current coronavirus epidemic.

Từ khóa


Tài liệu tham khảo

Huang C, Wang Y, Li X et al (2020) Clinical features of patients with 2019 novel coronavirus in Wuhan, China. Lancet. https://doi.org/10.1016/S0140-6736(20)30183-5

Lu R, Zhao X, Li J et al (2020) Genomic characterization and epidemiology of 2019 novel coronavirus: implications of virus origins and receptor binding. Lancet. https://doi.org/10.1016/S0140-6736(20)30251-8

Zhu N, Zhang D, Wang W et al (2020) A novel coronavirus from patients with pneumonia in China, 2019. N Engl J Med. https://doi.org/10.1056/NEJMoa2001017

Chan JF, Yuan S, Kok KH et al (2020) A familial cluster of pneumonia associated with the 2019 novel coronavirus indicating person-to-person transmission: a study of a family cluster. Lancet. https://doi.org/10.1016/S0140-6736(20)30154-9

WHO main website. https://www.who.int (accessed 5 Feb 5 2020)

Phan LT, Nguyen TV, Luong QC et al (2020) Importation and human-to-human transmission of a novel coronavirus in Vietnam. N Engl J Med. https://doi.org/10.1056/NEJMc2001272

National Health Commission of the People’s Republic of China. http://www.nhc.gov.cn (Assessed on 6 Feb 2020).

Rothe C, Schunk M, Sothmann P et al (2020) Transmission of 2019-nCoV infection from an asymptomatic contact in Germany. N Engl J Med. https://doi.org/10.1056/NEJMc2001468

World Health Organization (August 14, 2003). Alert, verification and public health management of SARS in the post-outbreak period (http://www.who.int.csr/sars/postoutbreak/en/). Accessed 3 Mar 2020

Chan HL, Tsui SK, Sung JJ (2003) Coronavirus in severe acute respiratory syndrome (SARS). Trends Mol Med 9:323–325

Leng Q, Bentwich Z (2003) A novel coronavirus and SARS. N Engl J Med 349:709

Lee TW, Cherney MM, Huitema C, Liu J, James KE, Powers JC, Eltis LD, James MNG (2005) Crystal structures of the main peptidase from the SARS coronavirus inhibited by a substrate-like aza-peptide epoxide. J Mol Biol 353:1137–1151

Shi JH, Wei Z, Song JX (2004) Dissection study on the SARS 3C-like protease reveals the critical role of the extra domain in dimerization of the enzyme: defining the extra domain as a new target for design of highly-specific protease inhibitors. J Biol Chem 279:24765–24773

Huang C, Wei P, Fan K, Liu Y, Lai L (2004) 3C-like proteinase from SARS coronavirus catalyzes substrate hydrolysis by a general base mechanism. Biochemistry 43:4568–4574

Chen H, Wei P, Huang C, Tan L, Liu Y, Lai L (2006) Only one protomer is active in the dimer of SARS 3C-like proteinase. J Biol Chem 281:13894–13898

Fan K, Wei P, Feng Q, Chen S, Huang C, Ma L, Lai B, Pei J, Liu Y, Chen J et al (2004) Biosynthesis, purification, and substrate specificity of severe acute respiratory syndrome coronavirus 3C-like proteinase. J Biol Chem 279:1637–1642

Sirois S, Wei DQ, Du Q, Chou KC et al (2004) J Chem Inf Comput Sci 44:1111–1122

Vastag B et al (2003) JAMA 290:1695–1696

Cao B, Wang Y, Wen D et al (2020) A trial of lopinavir-ritonavir in adults hospitalized with severe covid-19. N Engl J Med. https://doi.org/10.1056/NEJMoa2001282 published Online First: 2020/03/19

Benson D, Karsch-Mizrachi I, Lipman DJ, Ostell J, Wheeler DL et al.. "GenBank". Nucleic Acids Research. 36 (Database): D25–D30. (2008).doi:https://doi.org/10.1093/nar/gkm929. PMC2238942. PMID18073190

Rice P, Longden I, Bleasby A (2000) EMBOSS: the European molecular biology open software suite. Trends Genet 16(6):276–277

Etzold T, Argos P (1993) SRS, an indexing and retrieval tool for flat file data libraries. CompApplBiosci. 9:49–57

Aniba MR, Poch O, Thompson JD (2010) Issues in bioinformatics benchmarking: the case study of multiple sequence alignment. Nucleic Acids Res 38:7353–7363

Remmert M, Biegert A, Hauser A, Söding J (2012) HHblits: lightning-fast iterative protein sequence searching by HMM-HMM alignment. Nat Methods 9:173–175

Boeckmann B, Bairoch A, Apweiler R, Blatter MC, Estreicher A, Gasteiger E, Martin MJ, Michoud K, O'Donovan C, Phan I et al (2003) The SWISS-PROT protein knowledgebase and its supplement TrEMBL in 2003. Nucleic Acids Res 31:365–370

Benkert P, Biasini M, Schwede T (2011) Toward the estimation of the absolute quality of individual protein structure models. Bioinformatics. 27:343–350

Williams et al (2018) MolProbity: more and better reference data for improved all-atom structure validation. Protein Sci 27:293–315

Berg JM, John LT, Lubert S. "Chapter 2: Protein Composition and Structure." Biochemistry. New York: W. H. Freeman, 2007. N. pag. Print.

Bertoni M, Kiefer F, Biasini M, Bordoli L, Schwede T. Modeling protein quaternary structure of homo- and hetero-oligomers beyond binary interactions by homology. Scientific Reports. 7 (2017).

Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, Ferrin TE (2004) UCSF chimera—a visualization system for exploratory research and analysis. J Comput Chem 25(13):1605–1612

Delano WL (2002) The PyMOL Molecular Graphics System. http://www.pymol.org. http://ci.nii.ac.jp/naid/10020095229/en/ (Accessed 5 Aug 2020)

Chou PY, Fasman GD (1974) Prediction of protein conformation. Biochemistry. 13(2):222–245

Gasteiger E, Hoogland C, Gattiker A, Duvaud S, Wilkins MR, Appel RD, Bairoch A. Protein Identification and Analysis Tools on the ExPASy Server. In: John M. Walker (ed): The Proteomics Protocols Handbook, Humana Press; 2005. p. 571-607.

Jones DT, Taylor WR, Thornton JM (1992) The rapid generation of mutation data matrices from protein sequences. Comput Appl Biosci 8:275–282

Kumar S, Stecher G, Tamura K. (2015). MEGA7: Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets.Molecular Biology and Evolution.

Bunin BA, Siesel B, Morales G, Bajorath, J.(2007). Chemoinformatics: Theory, Practice, & Products. Dordrecht, the Netherlands: Springer Netherlands. p. 87. ISBN 978-1-4020-5000-8.

Hwangseo P, Jinuk L,Sangyoub L. Critical assessment of the automated AutoDock as a new docking tool for virtual screening, Proteins: Structure, Function, and Bioinformatics, 2006; 65(3):549–554, doi:https://doi.org/10.1002/prot.21183, PMID16988956.

Moult J, Fidelis K, Rost B, Hubbard T, Tramontano A (2005) Critical assessment of methods of protein structure prediction (CASP) round 6. Proteins 61(Suppl 7):3–7

Word JM, Lovell SC, LaBean TH, Taylor HC, Zalis ME, Presley BK, Richardson JS, Richardson DC (1999) Visualizing and quantifying molecular goodness-of-fit: small-probe contact dots with explicit hydrogen atoms. J MolBiol 285(4):1711–1733

Arendall WB, Tempel W, Richardson JS, Zhou W, Wang S, Davis IW, Liu ZJ, Rose JP, Carson WM, Luo M et al (2005) A test of enhancing model accuracy in high-throughput crystallography. J StructFunct Genomics 6(1):1–11

Lovell SC, Word JM, Richardson JS, Richardson DC (2000) The penultimate rotamer library. Proteins 40(3):389–408

Berkholz DS, Krenesky PB, Davidson JR, Karplus PA. Protein geometry database: a flexible engine to explore backbone conformations and their relationships to covalent geometry. Nucleic Acids Res 2010; 38:320–325. [PubMed: 19906726].

Woody RW. Circular dichroism spectrum of peptides in the poly(pro)II conformation. J Am Chem Soc 2009; 131:8234–8245. [PubMed: 19462996].

Gill SC, von Hippel PH (1989) Calculation of protein extinction coefficients from amino acid sequence data. Anal Biochem 182:319–326

Gonda DK, Bachmair A, Wunning I, Tobias JW, Lane WS, Varshavsky AJ (1989) Universality and structure of the N-end rule. J BiolChem 264(16):700–712

Tobias JW, Shrader TE, Rocap G, Varshavsky A (1991) The N-end rule in bacteria. Science. 254:1374–1377

Guruprasad K, Reddy BVB, Pandit MW (1990) Correlation between stability of a protein and its dipeptide composition: a novel approach for predicting in vivo stability of a protein from its primary sequence. Protein Eng 4:155–161

Ikai AJ (1980) Thermostability and aliphatic index of globular proteins. J Biochem 88:1895–1898

Abrusan G, Marsh JA. "Alpha helices are more robust to mutations than beta strands". PLoS Comput Biol. 2016; 12(12):1-16. doi:https://doi.org/10.1371/journal.pcbi.1005242. PMC 5147804 . PMID 27935949.

Rocklin GJ, et al. (2017). "Global analysis of protein folding using massively parallel design, synthesis, and testing". Science. 2017; 357(6347):168-175. doi:https://doi.org/10.1126/science.aan0693. PMC 5568797 . PMID 28706065.

Hasegawa M, Horai S (1991) Time of the deepest root for polymorphism in human mitochondrial DNA. J Mol Evol 32:37–42

Horai S, Hayasaka K (1990) Intraspecific nucleotide sequence differences in the major noncoding region of human mitochondrial DNA. Am J Hum Genet 46:828–842

Jin L, Nei M (1990) Limitation of the evolutionary parsimony method of phylogenetic analysis. Mol Biol Evol 7:82–102

Hetényi DC, van der S. (2002) Efficient docking of peptides to proteins without prior knowledge of the binding site. Protein Sci 11:1729–1737

Brooijmans N, Kuntz ID (2003) Molecular recognition and docking algorithms. Annu Rev Biophys Biomol Struct 32:335–373