Ribavirin, Remdesivir, Sofosbuvir, Galidesivir, and Tenofovir against SARS-CoV-2 RNA dependent RNA polymerase (RdRp): A molecular docking study

Abdo A Elfiky, Abdo A Elfiky

Abstract

Aims: A new human coronavirus (HCoV), which has been designated SARS-CoV-2, began spreading in December 2019 in Wuhan City, China causing pneumonia called COVID-19. The spread of SARS-CoV-2 has been faster than any other coronaviruses that have succeeded in crossing the animal-human barrier. There is concern that this new virus will spread around the world as did the previous two HCoVs-Severe Acute Respiratory Syndrome (SARS) and Middle East Respiratory Syndrome (MERS)-each of which caused approximately 800 deaths in the years 2002 and 2012, respectively. Thus far, 11,268 deaths have been reported from the 258,842 confirmed infections in 168 countries.

Main methods: In this study, the RNA-dependent RNA polymerase (RdRp) of the newly emerged coronavirus is modeled, validated, and then targeted using different anti-polymerase drugs currently on the market that have been approved for use against various viruses.

Key findings: The results suggest the effectiveness of Ribavirin, Remdesivir, Sofosbuvir, Galidesivir, and Tenofovir as potent drugs against SARS-CoV-2 since they tightly bind to its RdRp. In addition, the results suggest guanosine derivative (IDX-184), Setrobuvir, and YAK as top seeds for antiviral treatments with high potential to fight the SARS-CoV-2 strain specifically.

Significance: The availability of FDA-approved anti-RdRp drugs can help treat patients and reduce the danger of the mysterious new viral infection COVID-19. The drugs mentioned above can tightly bind to the RdRp of the SARS-CoV-2 strain and thus may be used to treat the disease. No toxicity measurements are required for these drugs since they were previously tested prior to their approval by the FDA.

Keywords: COVID-19; Drug repurposing; Molecular docking; RdRp; SARS-CoV-2; Structural bioinformatics.

Conflict of interest statement

Declaration of competing interest The author declares that there is no competing interest in this work.

Copyright © 2020 Elsevier Inc. All rights reserved.

Figures

Graphical abstract
Graphical abstract
Fig. 1
Fig. 1
(A) Phylogenetic analysis of the RdRps of 7 HCoVs (229E, NL63, HKU1, OC43, MERS, SARS, and Wuhan SARS-CoV-2) (cladogram); (B) Multiple sequence alignment of the HCoV RdRp sequences utilizing the Clustal Omega web server and represented by ESpript. Red highlights indicate identical residues while yellow-highlighted residues are less conserved. Secondary structures are described at the top of the MSA for SARS RdRp (PDB ID: 6NUR, chain A), while the surface accessibility is shown at the bottom (blue: highly accessible; white: buried). Active site aspartates (D255 and D256) are marked on the MSA. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.)
Fig. 2
Fig. 2
Bar graph representing the binding energies (in kcal/mol) calculated by AutoDock Vina software for GTP, UTP, CTP, ATP, Cinnamaldehyde, Thymoquinone, Galidesivir, Remdesivir, Tenofovir, Sofosbuvir, Ribavirin, Uprifosbuvir, Setrobuvir, Balaprevir, MK0608, R7128, IDX-184, 2′C-methylcytidine, BMS-986094, YAK, PSI-6130, PSI-6206, R1479, and Valopectibine against SARS-CoV-2 RdRp (blue), SARS HCoV RdRp (orange), and HCV NS5B RdRp (gray). The dashed-red rectangles mark the best three compounds. The dashed green box marks the five drugs (approved by FDA against HCV, EBOV, and HIV). (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.)
Fig. 3
Fig. 3
Interactions established after docking the (A) GTP, Ribavirin, Sofosbuvir, and Tenofovir and (B) IDX-184, Setrobuvir, and YAK against SARS-CoV-2 RdRp. Drugs are in orange, while cyan lines represent the active protein site pockets. Solid blue lines depict H-bonds, while hydrophobic interactions are gray dashed lines. Additionally, salt bridges, π-cation stacking, and halogen contacts are represented by yellow spheres connected by black dashed lines, orange dashed lines, and cyan lines, respectively. RdRp residues are labeled, and the docking scores are listed under each complex. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.)

References

    1. Bogoch I.I., Watts A., Thomas-Bachli A., Huber C., Kraemer M.U.G., Khan K. Pneumonia of unknown etiology in Wuhan, China: potential for international spread via commercial air travel. Journal of Travel Medicine. 2020;27(2)
    1. Hui D.S., I Azhar E., Madani T.A., Ntoumi F., Kock R., Dar O. The continuing 2019-nCoV epidemic threat of novel coronaviruses to global health—the latest 2019 novel coronavirus outbreak in Wuhan, China. Int. J. Infect. Dis. 2020;91:264–266.
    1. Rothan H.A., Byrareddy S.N. The epidemiology and pathogenesis of coronavirus disease (COVID-19) outbreak. J. Autoimmun. 2020 In press.
    1. Organization WH . World Health Organization; 2020. Surveillance Case Definitions for Human Infection With Novel Coronavirus (nCoV): Interim Guidance v1, January 2020.
    1. Organization WH . World Health Organization; 2020. Infection Prevention and Control During Health Care When Novel Coronavirus (nCoV) Infection Is Suspected: Interim Guidance, January 2020.
    1. Organization WH . World Health Organization; 2020. Laboratory Testing of Human Suspected Cases of Novel Coronavirus (nCoV) Infection: Interim Guidance, 10 January 2020.
    1. Parr J. Pneumonia in China: lack of information raises concerns among Hong Kong health workers. BMJ. 2020;368(56) doi: 10.1136/bmj.m56. Published 2020 Jan 8.
    1. Yang L. 2020. China Confirms Human-to-Human Transmission of Coronavirus.
    1. Chan J.F., Lau S.K., To K.K., Cheng V.C., Woo P.C., Yuen K.-Y. Middle East respiratory syndrome coronavirus: another zoonotic betacoronavirus causing SARS-like disease. Clin. Microbiol. Rev. 2015;28:465–522.
    1. Elfiky A.A., Mahdy S.M., Elshemey W.M. Quantitative structure-activity relationship and molecular docking revealed a potency of anti-hepatitis C virus drugs against human corona viruses. J. Med. Virol. 2017;89:1040–1047.
    1. Ibrahim I.M., Abdelmalek D.H., Elshahat M.E., Elfiky A.A. COVID-19 spike-host cell receptor GRP78 binding site prediction. J. Infect. 2020 In press.
    1. WHO . WHO; 2016. Middle East Respiratory Syndrome Coronavirus (MERS-CoV)
    1. Báez-Santos Y.M., Mielech A.M., Deng X., Baker S., Mesecar A.D. Catalytic function and substrate specificity of the papain-like protease domain of nsp3 from the Middle East respiratory syndrome coronavirus. J. Virol. 2014;88:12511–12527.
    1. Hemida M.G., Alnaeem A. Some one health based control strategies for the Middle East respiratory syndrome coronavirus. One Health. 2019;8
    1. Organization WH . World Health Organization; 2019. Clinical Management of Severe Acute Respiratory Infection When Middle East Respiratory Syndrome Coronavirus (MERS-CoV) Infection Is Suspected: Interim Guidance.
    1. Elfiky A.A. Zika viral polymerase inhibition using anti-HCV drugs both in market and under clinical trials. J. Med. Virol. 2016;88:2044–2051.
    1. Elfiky A.A. Zika virus: novel guanosine derivatives revealed strong binding and possible inhibition of the polymerase. Futur. Virol. 2017;12:721–728.
    1. Elfiky A.A. Novel guanosine derivatives as anti-HCV NS5b polymerase: a QSAR and molecular docking study. Med. Chem. 2019;15:130–137.
    1. Elfiky A.A., Elshemey W.M. IDX-184 is a superior HCV direct-acting antiviral drug: a QSAR study. Med. Chem. Res. 2016;25:1005–1008.
    1. Elfiky A.A., Elshemey W.M. Molecular dynamics simulation revealed binding of nucleotide inhibitors to ZIKV polymerase over 444 nanoseconds. J. Med. Virol. 2018;90:13–18.
    1. Elfiky A.A., Elshemey W.M., Gawad W.A., Desoky O.S. Molecular modeling comparison of the performance of NS5b polymerase inhibitor (PSI-7977) on prevalent HCV genotypes. Protein J. 2013;32:75–80.
    1. Elfiky A.A., Ismail A. Molecular dynamics and docking reveal the potency of novel GTP derivatives against RNA dependent RNA polymerase of genotype 4a HCV. Life Sci. 2019;238
    1. Elfiky A.A., Ismail A.M. Molecular modeling and docking revealed superiority of IDX-184 as HCV polymerase inhibitor. Futur. Virol. 2017;12:339–347.
    1. Ganesan A., Barakat K. Applications of computer-aided approaches in the development of hepatitis C antiviral agents. Expert Opin. Drug Discovery. 2017;12:407–425.
    1. Doublie S., Ellenberger T. The mechanism of action of T7 DNA polymerase. Curr. Opin. Struct. Biol. 1998;8:704–712.
    1. Elfiky A.A. Anti-HCV, nucleotide inhibitors, repurposing against COVID-19. Life Sci. 2020;248
    1. Elfiky A.A., Ismail A.M. Molecular docking revealed the binding of nucleotide/side inhibitors to Zika viral polymerase solved structures. SAR QSAR Environ. Res. 2018;29:409–418.
    1. Berman H., Henrick K., Nakamura H. Announcing the worldwide Protein Data Bank. Nat. Struct. Mol. Biol. 2003;10:980.
    1. Wu C., Liu Y., Yang Y., Zhang P., Zhong W., Wang Y. Analysis of therapeutic targets for SARS-CoV-2 and discovery of potential drugs by computational methods. Acta Pharm. Sin. B. 2020 In press.
    1. NCBI National center of biotechnology informatics (NCBI) database website. 2020.
    1. Biasini M., Bienert S., Waterhouse A., Arnold K., Studer G., Schmidt T. SWISS-MODEL: modelling protein tertiary and quaternary structure using evolutionary information. Nucleic Acids Res. 2014;42
    1. Kirchdoerfer R.N., Ward A.B. Structure of the SARS-CoV nsp12 polymerase bound to nsp7 and nsp8 co-factors. Nat. Commun. 2019;10:2342.
    1. SAVES . 2020. Structural Analysis and Verification Server Website.
    1. Williams C.J., Headd J.J., Moriarty N.W., Prisant M.G., Videau L.L., Deis L.N. MolProbity: more and better reference data for improved all-atom structure validation. Protein Sci. 2018;27:293–315.
    1. Laskowski R.A., Rullmann J.A.C., MacArthur M.W., Kaptein R., Thornton J.M. AQUA and PROCHECK-NMR: programs for checking the quality of protein structures solved by NMR. J. Biomol. NMR. 1996;8:477–486.
    1. Eisenberg D., Lüthy R., Bowie J.U. Methods in Enzymology. Elsevier; 1997. VERIFY3D: assessment of protein models with three-dimensional profiles; pp. 396–404.
    1. Joan Pontius JraSJW. Deviations from standard atomic volumes as a quality measure for protein crystal structures. J. Mol. Biol. 1996;264:121–136.
    1. Hooft R.W., Vriend G., Sander C., Abola E.E. Errors in protein structures. Nature. 1996;381:272.
    1. Elfiky AA. The antiviral Sofosbuvir against mucormycosis: an in silico perspective. Future Virology. 2019;14(11):739–744. 0:null.
    1. Summers K.L., Mahrok A.K., Dryden M.D., Stillman M.J. Structural properties of metal-free apometallothioneins. Biochem. Biophys. Res. Commun. 2012;425:485–492.
    1. Lii J.H., Allinger N.L. Molecular mechanics. The MM3 force field for hydrocarbons. 3. The van der Waals’ potentials and crystal data for aliphatic and aromatic hydrocarbons. J. Am. Chem. Soc. 1989;111:8576–8582.
    1. Trott O., Olson A.J. AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J. Comput. Chem. 2010;31:455–461.
    1. Salentin S., Schreiber S., Haupt V.J., Adasme M.F., Schroeder M. PLIP: fully automated protein–ligand interaction profiler. Nucleic Acids Res. 2015;43
    1. Becke A.D. Density-functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 1993;98:5648–5652.
    1. Leach A. 2nd edition. Prentice Hall; 2001. Molecular Modelling: Principles and Applications.
    1. Stewart J.J.P. Optimization of parameters for semiempirical methods V: modification of NDDO approximations and application to 70 elements. J. Mol. Model. 2007;13:1173–1213.
    1. Morris G.M., Huey R., Lindstrom W., Sanner M.F., Belew R.K., Goodsell D.S. AutoDock4 and AutoDockTools4: automated docking with selective receptor flexibility. J. Comput. Chem. 2009;30:2785–2791.

Source: PubMed

3
Tilaa