key: cord-0988223-jxbzts7v
authors: Singh, Sanjay Kumar; Upadhyay, Atul Kumar; Reddy, M. Sudhakara
title: Screening of potent drug inhibitors against SARS-CoV-2 RNA polymerase: an in silico approach
date: 2021-01-24
journal: 3 Biotech
DOI: 10.1007/s13205-020-02610-w
sha: 942f8cd68a5ce7dca818bf4ea054c06a7601be93
doc_id: 988223
cord_uid: jxbzts7v

COVID-19 has emerged as a rapidly escalating serious global health issue, affecting every section of population in a detrimental way. Present situation invigorated researchers to look for potent targets, development as well as repurposing of conventional therapeutic drugs. NSP12, a RNA polymerase, is key player in viral RNA replication and, hence, viral multiplication. In our study, we have screened a battery of FDA-approved drugs against SARS-CoV-2 RNA polymerase using in silico molecular docking approach. Identification of potent inhibitors against SARS-CoV-2 NSP12 (RNA polymerase) were screeened from FDA approved drugs by virtual screening for therapeutic applications in treatment of COVID-19. In this study, virtual screening of 1749 antiviral drugs was executed using AutoDock Vina in PyRx software. Binding affinities between NSP12 and drug molecules were determined using Ligplot(+) and PyMOL was used for visualization of docking between interacting residues. Screening of 1749 compounds resulted in 14 compounds that rendered high binding affinity for NSP12 target molecule. Out of 14 compounds, 5 compounds which include 3a (Paritaprevir), 3d (Glecaprevir), 3h (Velpatasvir), 3j (Remdesivir) and 3l (Ribavirin) had a binding affinity of − 10.2 kcal/mol, −9.6 kcal/mol, − 8.5 kcal/mol, − 8.0 kcal/mol and − 6.8 kcal/mol, respectively. Moreover, a number of hydrophobic interactions and hydrogen bonding between these 5 compounds and NSP12 active site were observed. Further, 3l (Ribavirin) was docked with 6M71 and molecular dynamic simulation of the complex was also performed to check the stability of the conformation. In silico analysis postulated the potential of conventional antiviral drugs in treatment of COVID-19. However, these finding may be further supported by experimental data for its possible clinical application in present scenario.

COVID-19 (coronavirus disease) emerged as grave health threat which started from Wuhan (China) in early December 2019 and rapidly spread across the world in few weeks (Loeffelholz and Tang 2020; Bogoch et al. 2020) . As of, 19 January 2021 over 93 million cases and 2 million deaths have been reported globally, with over 4.7 million new cases and nearly 93,000 new deaths reported over the past week. (WHO 2021) . SARS-CoV-2 belongs to Coronaviridae family, which causes infection in birds, mammals and humans (Rodriguez-Morales et al. 2020 ). In the last two decades, Coronaviruses were responsible for three huge outbreaks, which includes the SARS epidemic (2002-03) (Drosten et al. 2003) , the MERS epidemic (2012) (Zaki et al. 2012) and the present COVID-19 outbreak (Chan et al. 2020) .

SARS-CoV-2 is highly similar to other coronaviruses in terms of their genomic organization, which consist of single-stranded RNA as their genetic material (Mousavizadeh and Ghasemi 2020) . Proteins of Human coronaviruses are characterized into structural proteins and non-structural proteins. Structural proteins are further classified as S (spike), E (envelope), M (membrane), and N (nucleocapsid) proteins; whereas, non-structural proteins are categorized as proteases (NSP3 and NSP5) and RNA-dependent RNA polymerase (RdRp, also called as NSP12) (Elfiky 2017; Wu et al. 2020 ).

(2021) 11:93 1 3 93 Page 2 of 13 NSP12 catalyzes the viral RNA synthesis and, thus, plays a central role in viral replication as well as multiplication process, along with NSP7 and NSP8 proteins as co-factors (Gao et al. 2020) . NSP12 plays an important role in viral replication and has been targeted for treatment of several viral infections, including Hepatitis C, Zika and other Coronaviruses (Ganesan and Barakat 2017; Elfiky 2019 Elfiky , 2020 Elfiky and Elshemey 2018; Hui et al. 2020) . Till date, no specific treatment against SARS-CoV-2 has been reported despite of ardent efforts. Newer strategies for development or repurposing of conventional drugs are dire need for the present grave situation. Number of FDA approved antiviral drugs could have potential therapeutic options for SARS-CoV-2 infection. Also, a number of FDAapproved drugs are under investigation in randomized controlled clinical trials, including remdesivir (GS-5734) Gordon et al. 2020) . Nowadays, computational biology and molecular simulation studies open new possibilities of drug screening and may aid in the development of COVID-19 treatment (Dong et al. 2020) .

In this study, virtual screening of 1749 antiviral drugs was performed with NSP12 protein template by molecular docking and molecular dynamics simulation approach. The interactions of drug and protein moiety were also studied to understand the molecular interaction between drugs and proteins. The results have been found to be promising and indicate potential antiviral activity of selected drugs against SARS-CoV-2, which can be further validated through experimental settings.

The sequences corresponding to SARS-CoV nsp12 and SARS-CoV-2 nsp12 were retrieved from the NCBI Gen-eBank database. The related structural information of the SARS-CoV (PDB ID: 6NUR) and SARS-CoV-2 (PDB ID: 6M71) were collected from RCSB database (https ://www. rcsb.org/). Pairwise sequence comparisons of SARS-CoV-2 nsp12 and SARS-CoVnsp12 were performed using Clustal Omega webserver. The nsp12 protein from SARS-CoV exhibits 97% identity with SARS-CoV-2 nsp12.

For determining the specific functional part of protein, active site prediction is necessary. Potential active site of the target protein SARS-Cov-2 RNA-dependent RNA polymerase (6M71) was determined using PrankWeb (Jendele et al. 2019 ).

FDA-approved drugs dataset from Zinc15 library (Sterling and Irwin 2015) and antiviral agents (Accession Number: DBCAT000066) from DrugBank (Wishart et al., 2018) were retrieved for virtual screening. These molecules were downloaded in Structure Data File (SDF) format. Open Babel (http://openb abel.org) tool was used to convert various file formats. PyRx (Dallakyan and Olson 2015) was initially used to minimize compounds energy and convert all molecules to AutoDock Ligand (PDBQT) format.

Virtual screening was executed using the published crystal structure of SARS-CoV-2 nsp12 (PDB ID: 6M71). Auto-Dock Vina (Trott and Olson 2010) in PyRx 0.8 was used to perform virtual screening. The compounds without any specified binding sites were docked against whole surface of the protein. The value of the grid box was set to center_x = 121.8, center_y = 123.5, center_z = 127.0 while size_x = 74.8, size_y = 84.5, and size_z = 106.0 with the default exhaustiveness value of 8. The outcomes of docking results were reported in the form of binding energy. LigPlot + program was used for the analysis of post-docking results (Laskowski and Swindells 2011) . Using PyMOL, the docked complexes with lowest binding affinity values were further analyzed for hydrogen and hydrophobic bond interaction analysis. The resulting ligands have high potential to be used as drug candidates.

The ADME and drug-likeness predictions of compounds were carried out using SwissADME (http://www.swiss adme.ch/) (Daina et al. 2017) . The SMILES of compounds have been used in SwissADME web tool as input. Further, ADMET and the pharmacokinetic properties were evaluated using admetSAR (http://lmmd.ecust .edu.cn/admet sar2) web server to ensure the druggability potential of compounds (Yang et al. 2019) .

The structure of SARS-Cov-2 RNA-dependent RNA polymerase (PDB_ID: 6M71) was prepared by automated protein preparation protocol in Discovery Studio 4.1 (DS). Pre-processing of the selected structure in terms of fixing the missing atoms in the incomplete residues, deletion of co-crystallized H2O molecules, modeling of loops and protonation of selected residues were performed by applying CHARM force-filed. Side-by-side, the ligand Ribavirin (DB00811), a synthetic guanosine analog broad-spectrum antiviral effective against a number of RNA viruses, was processed using 'Prepare Ligands' protocol in DS 4.1 to generate protonation state, and have molecular geometry for docking and simulation studies. Molecular docking and molecular dynamics (MD) analysis of these pre-processed molecules were performed by 'define and edit binding site' tool to specify the active site. Ligand molecule was allowed to dock into the active site of 6M71 structure using CDOCKER tool of Discovery studio. Selection of the poses was done on the basis of interaction energy of these conformations. To check the stability of the docked complex, molecular dynamics simulation was performed for the period of 100 ns using academic version of "Desmond" program. An orthorhombic water box shape of 10 × 10 × 10 Å 3 dimensions was used as the solvent system using TIP3P. Neutralization of the system was done by adding Na + ions and 0.15 M salt concentration. Molecular dynamics simulations using NPT ensemble was performed with the temperature at 310 K and pressure at 1 bar using Nose-Hoover Chain as thermostat and Martyna-Tobias-Klein as barostat, respectively. The simulations were run for a period of 100 ns with time step of 1.0 femtoseconds (fs).

In Coronaviridae, the protein subunit containing the domain RdRp is recognized as non-structural protein NSP12. NSP12 is one of the most powerful enzymes essential to genome replication and transcription success of RNA viruses (Venkataraman et al. 2018 ). SARS-CoV-2 nsp12 is nearly similar to SARS-CoV (97% identity and 98% similarity) ( Fig. 1) . Most of the residues which were not identical have similar properties (strongly or weakly) and also overall architecture of nsp12-nsp7-nsp8 complex in both are almost similar.

The SARS-CoV-2 nsp12 structure includes a right-hand RdRp domain, and an N-terminal extension domain that is unique to nidovirus (Lehmann et al. 2015) . Architecture of polymerase domain is conserved in viral polymerase family (McDonald 2013) and it contains three subdomains; a finger (L366-A581 and K621-G679), a palm (residues T582-P620 and T680-Q815) and a thumb sub-domain (H816-E920) (Fig. 2) . The active site is established in the palm domain by conserved polymerase motifs (A-G) (Gao et al. 2020) .

Virtual screening (VS) has undeniably changed and enhanced the process of drug discovery, and has been considered as one of the most promising drug design techniques for identifying hit molecules as starting points in medicinal chemistry (Lavecchia et al. 2013; Lionta et al. 2014) . Recently in various research, FDA-approved drugs were repurposed against SARS-CoV-2 (Chandel et al. 2020a (Chandel et al. , 2020b . Compounds were downloaded in Structure Data File (SDF) format either from ZINC database or DrugBank. Energy of each compound was minimized and all converted to (PDBQT) format for virtual screening. Virtual screening was executed using the published electron microscopic structure of SARS-CoV-2 NSP12 PDB ID: 6M71 (Gao et al. 2020 ) with a resolution 2.90 Å. Table 1 describes the docking results which have been shown to be strong binding affinities. We studied the interactions between SARS-CoV-2 NSP12 and these strongest inhibitors. The graphical representation of the 14 best screened compounds is depicted in the pictorial form in Fig. 3 (Fig. 3) .

Compound 3a (CID 45110509), compound 3d (CID 66828839), compound 3 h (CID 67683363), compound 3j (CID 121304016) and compound 3 l (CID 37542) exhibit hydrogen bond interactions in the palm sub-domain which contain active site residues shown in Table 2 . Compound 3b, compound 3c, compound 3e, compound 3f, compound 3 g, compound 3i and compound 3 k were showing good binding affinity, hydrogen bond and hydrophobic interaction but not in the active site containing palm sub-domain.

The drug-likeness predictions of 14 compounds were carried out using SwissADME. Lipinski 's rule of 5 states that a molecule should meet the following conditions for any ligand to be considered drug-like: molecular weight < 500 Dalton, number of donors of H-bonds < 5, number of acceptors of H-bonds < 10 and LogP < 5 (Lipinski 2004) . Out of 14 FDA-approved compounds, only 3 (3f, 3 l and 3n) followed the Lipinski's rule of five and are provided in Table 3 .

Further, all these 14 best docked ligands were further subjected to admetSAR for evaluation of ADMET properties (Table 4 ). Blood-Brain Barrier penetration was shown by all except 3c. Hepatotoxicity was only shown by 3e. The acute oral toxicity for the all compounds was estimated as class III. All compounds demonstrated no Caco-2 permeability. Similarly, intestinal absorption (human) was observed in all compounds. Honey bee toxicity was detected in 3b, 3c, 3e and 3 h compounds. Most of the ligands were CYP3A4 substrate except 3f, 3 k, 3 l, 3 m and 3n compounds, which were non-substrates. Only four ligands (3i, 3 l, 3 m and 3n) were non-inhibitors of CYP3A4, while rest were inhibitors. The water solubility of most of the docked compounds is greater than −3.

Compound 3a (Paritaprevir) is an antiviral drug used to treat chronic hepatitis C as part of the combination treatment (Wilkins et al. 2015) . This drug was also studied against nsp 15 protein by Khan et al. (2020) and against nsp13 helicase and nsp14 of SARS-CoV-2 by Gurung (2020) . The compound 3a showed the highest binding affinity (− 10.2 kcal/ mol) and had hydrogen bonds with amino acids Arg553, Asp623, Thr687, Ser759; whereas hydrophobic interaction with Tyr455, Tyr619, Pro620, Lys621, Cys622, Ser681, Ser682, Ala688, Asn691 and Asp760 residues (Fig. 4a, b) .

Compound 3d (Glecaprevir) has been reported to have the binding affinity (− 10.1 kcal/mol) to the receptor via hydrogen bond interaction with Arg553, Arg555, Tyr619, Cys622; while, the hydrophobic interaction was with the residues Thr556, Asp618, Pro620, Lys621, Asp623, Ser682, Thr687, Asp760 and Asp761 (Fig. 5a, b) . Glecaprevir is an antiviral agent and NS3/4A protease inhibitor of the Hepatitis C virus (HCV) that inhibits the replication of viral RNA (Zeuzem et al. 2018) . Glecaprevir is a useful therapy in conjunction with Pibrentasvir for patients who have encountered clinical failure from other inhibitors of the NS3/4A protease (Zeuzem et al. 2018) . Glecaprevir was also investigated for the treatment of COVID-19 infection (Mahdian et al. 2020; Chtita et al. 2020) .

Compound 3 h (Velpatasvir), a NS5A inhibitor is a part of combination therapy to treat chronic hepatitis C (Heo and Deeks 2018) . With limited side effects, usually fatigue and headache, Velpatasvir may be attractive as a therapeutic to fight the new coronavirus (Chen et al. 2020) . It was reported that velpatasvir has the binding affinity (− 8.5 kcal/mol) and bound to the receptor through hydrogen bond interaction with Arg33, Tyr129, Thr141, Lys780; whereas, the hydrophobic interaction was with Val31, Tyr32, Lys47, Tyr122, Asn138, Asp140, Ser709, Thr710, Gly774 and Asn781 residues (Fig. 6a, b) . CID 45110509 − 10.2 − 10 − 9.9 − 9.7 − 9.7 − 9.6 − 9.6 − 9.5 − 9.5 CID 71171 − 10.1 − 9.6 − 9.5 − 9.5 − 9.4 − 9.2 − 9.2 − 9.1 − 9 CID 8223 − 10. Compound 3j (Remdesivir) is a promising candidate based on existing evidence, and clinical trials of remdesivir are now underway among adults with COVID-19 in hospitals (NIH 2020). Remdesivir is a direct-acting antiviral agent used for the treatment of serious 2019 coronavirus disease (COVID-19) patients Baby et al. 2020) . It showed a mixed outcome with a reasonable side effect in COVID-19 patients . It was found that Remdesivir has the binding affinity (− 8.0 kcal/mol), and residues Lys47, Tyr129, Ser709, Asp711 and Asn781 were connected with hydrogen bonds; while, residues Tyr32, Ala46, His133, Asp135, Thr710, III   III   III   III   III   II   III   III   III   III   III   III   III Gln773, Gly774, Lys780 and Ser784 had hydrophobic interactions (Fig. 7a, b) . Compound 3 l (Ribavirin) is a synthetic guanosine analog broad-spectrum antiviral effective against a number of RNA viruses (Galli et al. 2018) . Our results of molecular docking and molecular dynamic simulation have confirmed that the Ribavirin have best interactions with 6M71 and also showed energetically stable conformations along with best druggability properties (Fig. 8) . In silico screening, molecular mechanics, and molecular dynamics simulation (MDS) research indicate that at the receptor-binding domain that recognizes hACE-2 (RBD-hACE2) interface, ribavirin has strong interaction (Tiwari et al. 2020) . It has the binding affinity (− 6.8 kcal/mol) and formed 5 hydrogen bonds (Trp617, Asp761, Lys798, His810 and Glu811) and 3 hydrophobic contacts (Trp800, Phe812, Cys813) (Fig. 9a, b) .

Many studies are underway in establishing conventional antiviral drugs to counter SARS-CoV-2 infection. In a short time, a number of potential drug inhibitors have been screened and identified through docking-based virtual screening using various in silico approaches. In this study, docking studies were performed for NSP12 protein and 1615 FDA-approved antiviral compounds from ZINC15 database and 149 molecules of known antiviral activities. Out of 1764 screened compounds, 14 compounds exhibited strong binding affinities with NSP12 protein. Based on binding affinity, 5 compounds which include 3a (Paritaprevir), 3d (Glecaprevir), 3 h (Velpatasvir), 3j (Remdesivir) and 3 l (Ribavirin) were selected for further analysis. These five drugs showed high binding affinity, and also exhibited hydrophobic interactions and hydrogen bonding in the palm sub-domain of NSP12. 3 l (Ribavirin) was docked with 6M71 and molecular dynamic simulation of the complex was also performed to check the stability of the conformation. This study established the potential of conventional antiviral drugs in treatment of SARS-CoV-2 infection. The resulting compounds have high potential to be used as drug candidates. However, more extensive screening along with experimental results may also support the results of present finding and can open new future possibilities for drug discovery against SARS-CoV-2 infection.

Targeting SARS-CoV-2 main protease: a computational drug repurposing study

Pneumonia of unknown etiology in Wuhan, China: potential for international spread via commercial air travel

Genomic characterization of the 2019 novel human-pathogenic coronavirus isolated from a patient with atypical pneumonia after visiting Wuhan

Structure-based drug repurposing for targeting Nsp9 replicase and spike proteins of severe acute respiratory syndrome coronavirus 2

In silico identification of potent FDA approved drugs against Coronavirus COVID-19 main protease: a drug repurposing approach

Prediction of the SARS-CoV-2 (2019-nCoV) 3C-like protease (3CL pro ) structure: virtual screening reveals velpatasvir, ledipasvir, and other drug repurposing candidates

Discovery of potent SARS-CoV-2 inhibitors from approved antiviral drugs via docking screening

SwissADME: a free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules

Small-molecule library screening by docking with PyRx

A guideline for homology modeling of the proteins from newly discovered betacoronavirus, 2019 novel coronavirus (2019-nCoV)

Identification of a novel coronavirus in patients with severe acute respiratory syndrome

Zika virus: novel guanosine derivatives revealed strong binding and possible inhibition of the polymerase

Novel guanosine derivatives as anti-HCV NS5b polymerase: a QSAR and molecular docking study

Anti-HCV, nucleotide inhibitors, repurposing against COVID-19

Molecular dynamics simulation revealed binding of nucleotide inhibitors to ZIKV polymerase over 444 nanoseconds

Antiviral effect of ribavirin against HCV associated with increased frequency of G-to-A and C-to-U transitions in infectious cell culture model

Applications of computer-aided approaches in the development of hepatitis C antiviral agents

Structure of the RNA-dependent RNA polymerase from COVID-19 virus

Remdesivir is a direct-acting antiviral that inhibits RNA-dependent RNA polymerase from severe acute respiratory syndrome coronavirus 2 with high potency

In silico structure modelling of SARS-CoV-2 Nsp13 helicase and Nsp14 and repurposing of FDA approved antiviral drugs as dual inhibitors

Sofosbuvir/Velpatasvir/Voxilaprevir: a review in chronic hepatitis C

The continuing 2019-nCoV epidemic threat of novel coronaviruses to global health: the latest 2019 novel coronavirus outbreak in Wuhan, China

Prank-Web: a web server for ligand binding site prediction and visualization

Identification of promising antiviral drug candidates against non-structural protein

-2: an in silico assisted drug-repurposing study

LigPlot+: multiple ligand-protein interaction diagrams for drug discovery

Virtual screening strategies in drug discovery: a critical review

Discovery of an essential nucleotidylating activity associated with a newly delineated conserved domain in the RNA polymerase-containing protein of all nidoviruses

Updated approaches against SARS-CoV-2

Structure-based virtual screening for drug discovery: principles, applications and recent advances

Lead-and drug-like compounds: the rule-of-five revolution

Laboratory diagnosis of emerging human coronavirus infections-the state of the art

Drug repurposing using computational methods to identify therapeutic options for COVID-19

RNA synthetic mechanisms employed by diverse families of RNA viruses

Genotype and phenotype of COVID-19: Their roles in pathogenesis

National Institutes of Health, NIH Clinical Trial of Remdesivir to Treat COVID-19 Begins

History is repeating itself: probable zoonotic spillover as the cause of the 2019 novel coronavirus epidemic

Remdesivir in COVID-19: a critical review of pharmacology, pre-clinical and clinical studies

ZINC 15-ligand discovery for everyone

Template-dependent inhibition of coronavirus RNA-dependent RNA polymerase by remdesivir reveals a second mechanism of action

Denovo designing, retrosynthetic analysis, and combinatorial synthesis of a hybrid antiviral (VTAR-01) to inhibit the interaction of SARS-CoV2 spike glycoprotein with human angiotensin-converting enzyme 2

AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading

RNA dependent RNA polymerases: insights from structure, function and evolution

Diagnosis and management of hepatitis C

DrugBank 5.0: a major update to the DrugBank database for 2018

Analysis of therapeutic targets for SARS-CoV-2 and discovery of potential drugs by computational methods

admetSAR 2.0: web-service for prediction and optimization of chemical ADMET properties

Isolation of a novel coronavirus from a man with pneumonia in Saudi Arabia

Glecaprevir-pibrentasvir for 8 or 12 weeks in HCV genotype 1 or 3 infection

Acknowledgements Authors acknowledge Department of Biotechnol-