key: cord-0872664-oz99upzz
authors: Ahmad, Jamshaid; Ikram, Saima; Ahmad, Fawad; Rehman, Irshad Ur; Mushtaq, Maryam
title: SARS-CoV-2 RNA Dependent RNA Polymerase (RdRp) – A drug repurposing study
date: 2020-07-23
journal: Heliyon
DOI: 10.1016/j.heliyon.2020.e04502
sha: 6cccfdba0fa25743534071476f1251e9da3bfbe8
doc_id: 872664
cord_uid: oz99upzz

The outbreak of SARS-CoV-2 in December 2019 in China subsequently lead to a pandemic. Lack of vaccine and specific anti-viral drugs started a global health disaster. For a sustained control and protection, development of potential anti-viral drugs is one of the targeted approach. Although, designing and developing a panel of new drugs molecules are always encouraged. However, in the current emergency, drug repurposing study is one of the most effective and fast track option. The crystal structure of a SARS-CoV-2 (Severe acute respiratory syndrome coronavirus 2) RNA Dependent RNA Polymerase (RdRp) has recently been deciphered through X-ray crystallography. The single-chain of core RNA Dependent RNA Polymerase relies on virus-encoded cofactors nsp7 and two units of nsp8 for its optimum function. This study explored the FDA approved database of 7922 molecules and screened against the core polymerase along with cofactors. Here we report a panel of FDA approved drugs that show substantial interactions with key amino acid residues of the active site. Interestingly, some of the identified drugs (Ornipressin, Lypressin, Examorelin, Polymyxin B1) bind strongly within the binding pockets of both forms of RdRp. Besides, we found strong candidates for the complex form as well which include Nacortocin, Cistinexine, Cisatracurium (among others). These drugs have the potential to be considered while contriving therapeutic options.

COVID-19, an infectious disease caused by a novel strain of coronavirus SARS-CoV-2 (Severe Acute Respiratory Syndrome). It caught the attention of the world community after its outbreak in December 2019 from the city of Wuhan, China. Since then it has spread globally, WHO declared COVID-19 as pandemic on 11 March 2020 [1, 2] . The emergence of this new strain of coronavirus (SARS-CoV-2) which is highly contagious and caused the deaths of 444,813 (as of 17 th , June 2020, WHO) individuals worldwide, making it as one of the biggest challenge for the scientific community worldwide. In this global health emergency, drug repurposing (or repositioning) is one of the fast track option that involves screening of existing FDA approved drugs for the identification of potential molecules that can disrupt the function of key proteins of the SARS-CoV-2 and can be used for treatment against COVID-19. Coronaviruses are enveloped viruses belonging to family Coronoviridae, which is further subdivided into four genera i.e. alpha, beta, gamma, and delta coronavirus [3] . The members of the alpha and beta coronavirus causes infections in mammals however gamma and delta coronavirus causes disease in birds [4, 5] . SARS-CoV-2 belongs to genus beta-coronavirus. Other members of the same genus that are known to infect humans include HCoV-OC43, HCoV-HKU1, SARS-CoV, and MERS-CoV [4, 5] . SARS-CoV-2 genome is a single strand positive-sense RNA. The size of the genome is 29.8kb and has 14 Open Reading Frames (ORFs) which encode information for 27 structural and non-structural proteins [4] . At the 5`-end of the genome, there are two ORFs (ORF-1a and ORF-1ab) encoding two long stretches of poly-proteins, pp-1a and pp-1ab. These encode information for 15 non-structural proteins (nsp 1-10 and nsp12-16) [4] . The vital nonstructural proteins are; nsp3 (multi-domain protein including PL-pro domain), nsp5 (3CL chymotrypsin-like), nsp9 (helicase, may participate in viral replication), nsp12 (RNA dependant RNA polymerase) and nsp13 (helicase). The 3`-end of the genome encodes information for four structural and eight accessory proteins. The structural proteins are; Spike surface glycoproteins (S), Envelope (E), Matrix (M), and Nucleocapsid (N) proteins. The accessory proteins are; 3a, 3b, p6, 7a, 7b, 8b, 9b, and orf14 [4] .

RNA dependent RNA polymerase (RdRp) is involved in the replication and transcription of the SARS-CoV-2 genome. It is the cleavage product of the polyproteins 1a and 1ab from ORF1a and ORF1ab [6] . There is a high degree of conservation among RNA dependent RNA polymerases of different RNA viruses. The Core protein is a single chain of approximately 900 amino acids, showing minimal activity. However, enhanced activity is attained with the attachment of other key subunits [6] [7] [8] . The minimum complex for the proper activity of SARS-CoV-2 requires an attachment of nsp7/nsp8 and an additional nsp8 protein to the core protein [6] . The attachment site for the second additional nsp8 is different. Also, other nonstructural proteins (nsp (s)) are involved in this replication and transcription cascade [6] . The core protein looks like a cupped right hand which is further sub-divided into subdomain including finger domain (amino acid residues 398-581, 628-687), palm domain (582-627, 688-815), and thumb domain (816-919) [9] (Figure 1 ). Two additional Zn ions are also required for the structural stability of the RdRp. One of the Zn ion is attached to four amino acids residue (His295, Cys301, Cys306, and Cys310) in the N-terminal domain while the second Zn ion is attached to four amino acid residues (Cys487, His642, Cys645, Cys646) located in finger domain [6] . The presence of Zn in this location shows that it plays an essential role in stabilizing the overall three-dimensional structure of the protein. It has no direct role in the activity of the polymerase as both are quite distal to the active site. All the important sites including, template entry and binding, polymerase activity reaction site followed by the exit through the tunnel (thumb) remain highly conserved among coronaviruses including SARS-CoV-2 [6] (Figure 1 ). To the extent charge potential is concerned, the palm (RNA template and NTP binding site) has positive electrostatic potential. However, the other sites including attachment of nsp7/8 complex, nsp8 site, and the template exit site (thumb) is neutral [6] . The key residues that are involved in the interaction include; Tyr618, Cys622, Asn691, Asn695, Met755, Ile756, Leu757, Leu758, Ser759, Asp760, Asp761, Ala762, Val763, Glu811, Phe812, Cys813 and Ser814 (numbering about the recently solved structure PDB ID 6M71) [10, 11] . The active site key residues are adjacent aspartates i.e. Asp761 and Asp762, which are involved in the actual reaction of the RdRp enzyme [10, 12] .

The panel of different antiviral molecules has already been identified by targeting key proteins involved in different stages of the SARS-CoV-2 life cycle. Some of them are in a clinical trial stage. Through repurposing study, the existing antiviral drugs (Ribavirin, Sofosbuvir, Remdesiver, Tenofovir) have already been screened [10, 12] . Those in a clinical trial stage include; ASC09 and Ritonavir (both are HIV inhibitor), GS-5734 (a broadspectrum investigational antiviral), and Kaletra (inhibitor of 3-CL of SARS and MERS) [13] . Similarly, other compounds in a clinical trial (Favipiravir, used against influenza) didn't show strong activity in clinical isolates [13] , while some others with positive results need further investigations [14, 15] . Owing to its significance, the current study encompasses screening and identifying novel FDA approved drugs against the recently elucidated structure of SARS-CoV RNA Dependent RNA Polymerase. The possible lead molecules can further be explored to combat this disease.

The crystal structure of RdRp was downloaded from PDB (ID 6M71) with a reported resolution of 2.90Å. RNA Dependent RNA polymerase is a multimeric protein. The minimum complex required for its proper functioning is completed by attachment of three additional protein peptides (nsp7-nsp8, and one additional nsp8) to the core polymerase which is chain A and contains 851 amino acids residues. The current structure is unraveled in complex with all the required cofactors. To investigate the binding affinity of the drugs, we prepared both the structures of RdRp i.e. with and without cofactors. Protein preparation wizard of Maestro was used for protein preparation. Hydrogen atoms, missing residues, and loops were added. The side chains were fixed during protein preparation of both forms of RdRp (with and without cofactors). Besides, disulfide bonds were also created. Protonation states of amino acids were generated at pH 7.4, by using PROPKA to simulate physiological conditions [16, 17] . For the protein structures minimization, the OPLS-2005 force field was used [18] . Grid box was generated by selecting the active site residues as mentioned [10] . Active site residues Asp760 and 761 (DD) were treated as flexible during docking. Grid box was generated as 30, 30, 30 centred (x, y, z) of (114. 52, 114.11, 122 .91) Å as reported [10] . Receptor grid was generated using Maestro tools and the flexible rotation of certain amino acids residues (Tyr455, Tyr456, Thr462, Cys482, Tyr483, Ser549, Thr586, Thr591, Thr604, Tyr619, Cys622, Thr680, Ser681, Ser682, Thr686, Thr687, Tyr689, Ser692, Cys697, Thr701, Ser754, Ser759, Cys765, Ser778, Tyr788, Ser795, Cys799, Thr801, Cys813, Ser814, Thr817) were allowed. FDA approved drug database (~7922) was downloaded from the NIH chemical genomics Centre (NCGC) pharmaceutical collection (NPC) database. The database was prepared using the lig-prep module of Maestro [18] . During protein preparation, different combinations of enantiomers and tautomers were generated, which ultimately enhanced the total number of drug molecules. For the protonation states of the ligand at pH 7.4, Epik was used [19] . To check the binding efficacy of the FDA approved drugs, accurate and fast Glide SP-protocol was used for the molecular docking simulations [20, 21] . Glide Induced Fit Docking (IFD) was used to re-dock the top ten molecules in the binding pocket of RdRp [20] . IFD was performed in various steps. Initially, ligands were docked to rigid protein employing potential softened-potential glide docking with the van dar Waals radii scaling of 0.50/0.50, for the receptor and ligand respectively. In the following stage receptor, sampling and refinement were carried out. Residues with minimum one atom inside 5 Å of any of the 10 ligand poses were selected to conformational search and minimization whereas residues lying outside this range were fixed. In this manner, the protein flexibility was considered in the process. Additional redocking of the ligands were performed. For all docking calculations, Glide-XP was analyzed. 3D structure of the protein model, 2D structures of molecules and interaction maps were generated using Maestro (Schrodinger). For other illustration, Adobe Illustrator was used.

For the screening of novel potential inhibitors against RdRp, we selected FDA approved drugs database. A similar type of repurposing study has already been performed, in which the homology model of RdRp was generated, validated and various existing antiviral drugs were investigated for their binding potential within the active site of RdRp without cofactors [10, 12] . Furthermore, other proteins of the SARS-CoV-2 for which X-ray crystal structures were available have also been used as a potential drug target [22] . The current study is the first of its type in which recently determined X-Ray crystal structure of RNA Dependent RNA polymerase was used for the screening of potential inhibitors from the FDA approved drugs database. In the homology model, two aspartates at position 255 and 256 were reported as active site residues [12] . However, in the crystal structure, these similar active site aspartates are represented as Asp760 and Asp761. We have used both forms of the polymerase, a single-core chain, and along with its cofactors (holoenzyme, replication/transcription complex). Initially, a single chain of the RdRp polymerase was selected for the screening, followed by further screening for the holoenzyme.

Based on docking scores (Glide SP, Glide IFD) and interactions with crucial residues and ligand efficiencies, the top-ranked molecules were selected and analyzed. Table 1 explains the review of the top ten screened molecules and the existing usage of these drugs. Antibiotic used against multidrugresistant gramnegative bacterial infection [39, 40] Based on the maximum occupancy of compounds bound to the reported binding tunnel of RdRp, top-ranked conformations were selected considering the hydrogen bonds and other non-covalent interactions. Ornipressin binds into the RdRp binding tunnel with a docking score of -10.37kcal/mol and the ligand efficiency of -0.14. Ligand interaction analysis of the Ornipressin/RdRp complex shows that ligand mostly made H-bonds and salt bridges. Asp865 is involved in H-bond and salt bridge while Asp760, Thr591, Gln815, Ser814, Cys813, Glu811, Tyr619 formed H-bonds. Ornipressin also depicted some interactions with the hydrophobic and polar residues (Figure 2 ). Figure S2) . The binding potential of Argiprespocin with RdRp was mainly governed by H-bonds. Argiprespocin made six H-bonds with Lys798, Glu811, Leu758, and Ser814. While two H-bonds and one salt bridge were observed for Asp761, one π-cation and H-bond was observed for Trp617 with a docking score of -9.87kcal/mol and ligand efficiency of -0.137 (Supplementary Figure S2) . Whereas, Demoxytocin showed ten H-bonds with both active site Asp760 and Asp761 and other key residues e.g. Trp617, Tyr619, Lys621, Ser682, Glu811, Lys621, Tyr619, Trp617, Ser682 and Glu811 with dock score -9.68kcal/mol and ligand efficiency of -0.142 (Supplementary Figure S3) . Carbetocin, a synthetic analogue of oxytocin, also has the potential to inhibit RdRp (single chain) and showed stable interaction with a docking score of -9.479 kcal/mol and ligand efficiency of -0.137. Eight H-bonds were observed with Lys551, Arg553, Thr556, Trp617, Tyr619, Asp623, and Asp760 (Supplementary Figure S3) .

Another screened drug Lypressin showed a docking score of -9.413 kcal/mol and ligand efficiency of -0.129. Ligand interactions showed salt bridge and H-bonding with Asp761 and only H-bonding with other residues including Arg553, Arg555, Thr556, Trp617, Tyr619, Arg624, Asp760, Asp761, Ala762 and Glu811 (Supplementary Figure S4) . Examorelin, a screened inhibitor of this study showed most of H-bonds with Asp760, Cys813, Ser814, Gln815, and Glu811. Salt bridge was also observed with Glu811. Arg555 showed π-cation interaction with a docking score of -9.30kcal/mol and ligand efficiency of -0.143 (Supplementary Figure S4) . Colistin (polymyxin E, polypeptide antibiotics) showed most of the H-bonding with Lys551, Trp617, Tyr619, Asp618, Ser682, Asp684, Asn691, and both catalytic residues i.e. Asp760, Asp761, with a docking score of -9.24kcal/mol and ligand efficiency of -0.113 (Supplementary Figure S5) . Polymyxin B1 also bound tightly with RdRp core protein through salt bridges and H-bonding. Catalytic Asp760 and Asp761 made four salt bridges and three H-bonds with RdRp which shows its strong binding affinity. Moreover, H-bonds were formed with Arg-555, Trp-617, Asp681, Asn691, Cys813, and Ser814 with a docking score of -9.22kcal/mol and ligand efficiency of -0.109 (Supplementary Figure S5) .

To investigate the difference in binding energies, molecules were re-docked in the binding pocket of RdRp by using the Glide-Induced Fit Docking protocol. No significant changes were observed for drugs (except for Colistin which didn't show any pose in IFD) in both the protocols. A similar type of interaction was observed in Glide-IFD as found in Glide-SP (Table 1) . Glide-IFD docking poses of selected top 10 molecules are given in Supplementary figures (Figure S6 and S7) . All complexes showed promising binding and maximum occupancy of crucial binding pocket residues. It has been observed that the top ten screened drugs binding affinity was stronger with RdRp core protein.

As previously reported, Sofosbuvir, Ribavirin, Tenofovir, Setrobuvir, Guanosine derivative (IDX-184) were shown with a docking score of -7.5, -7.8, -6.9, -9.3, -9.0 kcal/mol respectively against RdRp single-chain [12] . In current study, we have also docked Ribavirin, Remdesivir, Sofosbuvir, Galidesivir, Tenofovir, Guanosine derivative (IDX-184) and Setrobuvir with a docking score of -5.94, -5.57, -4.94, -7.20, -4.52, -8.57, -4.45 kcal/mol respectively (Supplementry Figure S8) .The difference in the docking score may be due to the use of different docking algorithms. These antivirals were analysed on the basis of interactions with crucial residues and compared the interactions with screened molecules. Our identified molecules showed more stable H-bonds, salt bridges, and ionic interactions as compared to reported antivirals. For a potent drug, it is essential to interact with key amino acid residues in the active site pocket of the polymerase. Once the drug binds to the active site of RdRP, the access of positive sense RNA of the virus which serves as a template for anti-genome is blocked. No new viral genomes are synthesized and mature virus particles cannot be assembled and propagated further. In this study, all the identified FDA approved drugs have shown stable interaction with key residues, along with lowest binding energies which indicates the potential for these drugs to inhibit the activity of RdRp, therefore contributing to stop viral proliferation. The URL link of additional top 90 molecules including 2D structures, docking scores, ligand efficiencies and glide energies are given in the supplementary data.

The minimal complex of RdRp along with essential cofactors (nsp7-nsp8, nsp8) required for enhanced activity is used for this drug repurposing study. We found a panel of molecules that showed significant binding with the RdRp complex. Molecules were analyzed based on docking energies, ligand efficiencies, and other protein-ligand interactions. Some of the screened molecules remained similar as found for single chain of RdRP, they therefore can be taken as common inhibitors for both single chain and in complex form of RdRp. The results of top 14 selected molecules from Glide/SP along with their 2D structures, ligand efficiency, interacting residues and existing usage of the drugs are shown in table 2. Top 14 molecules including common molecules among both structures were redocked through Glide IFD. Docking score of Gide/IFD are shown in table 2. For a complex, increased docking score was observed in Glide/SP. This phenomenon might be due to the allosteric behaviour of the polymerase, as it is a common observation that upon binding of the molecules on an allosteric site, it brings conformational changes in the active site. This then facilitates the efficient binding of the substrate within the active site. The attachment of cofactors (nsp8/nsp7, nsp8) with a core protein (nsp12) might have similar effect on the activity of RdRp. Additionally, the binding of some of the molecules in both the forms depend on the chemical structures of the screened molecules, changed interacting residues, steric hindrance, and stereo chemistry of the active site. Ventilatorassociated pneumonia [37] , nosocomial pneumonia [38] , antibiotic against Nacartocin is one of the lead molecules with a docking score of -13.943 kcal/mol and ligand efficiency of -0.202. Nacartocin showed H-bonds and non-covalent interactions with key residues and surface accessible residues (Figure 3 ). π -cation interaction with Arg-553 and non-covalent interactions with crucial residues were observed for Cistinexine, with a docking score of -13.687 kcal/mol and a ligand efficiency of -0.201. Whereas a neurotransmitter blocking agent Cisatracurium showed two H-bonds with Arg553 and one salt bridge with Glu811 with a docking score of -13.064, and a ligand efficiency of -0.195. Other non-covalent interactions were also observed with polar, nonpolar, and hydrophobic residues (Supplementary Figure S9) . Pegamotecan is another lead molecule that showed a strong binding affinity towards RdRp with their essential cofactors. H-bond with Lys621 and other non-covalent interactions with key residues were observed with a docking score of -12.782 kcal/mol and ligand efficiency of -0.178 respectively. Among the top 14 screened drugs Polymyxin B1 is one of the common drugs, identified in both forms of RdRp. Almost similar types of interactions were observed for Polymyxin B1 in the catalytic site of the RdRP complex as found during core protein interactions. Polymyxin B1 showed two salt bridges and one H-bond with Asp760 and Asp761 and two H-bonds with Arg-555 with a docking score of -12.646 kcal/mol and ligand efficiency of -0.149 (Supplementary Figure S9) . Non-covalent interactions and one H-bond with Trp800 were observed for Ediratide with a docking score of -12.682 kcal/mol and with ligand efficiency of -0.18. Asp760, Ser759, and Glu811 made H-bonds with Sulfomyxin through docking score of -12.468 kcal/mol and ligand efficiency of -0.152 (Supplementary Figure S10) . π-cation with Arg553 and H-bonds with Tyr619, Arg86, and Asp833 were observed for Diagastrin, which is a peptide analogue of gastrin showed a docking score of -12.313 kcal/mol and ligand efficiency of -0.173. Two salt bridges with Asp760 and Asp623 were observed for Ditercalinium chloride with a docking score of -12.308 kcal/mol and ligand efficiency of -0.228 (Supplementary Figure  S10 ).

Flavonoid drug, Benzquercin showed two π-cation interactions with Arg553 and Lys798 and other non-covalent interactions with polar, non-polar and hydrophobic residues with docking score of -12.174 kcal/mol and ligand efficiency of -0.214 (Supplementary Figure S11) .

Examorelin, Lypressin, Ornipressin, and Colistin are also common drugs in both form of RdRp. Only one H-bond with His810 and other non-covalent interactions were observed for Examorelin showed a docking score of -12.139 kcal/mol and ligand efficiency of -0.187. Lypressin made three H-bonds with Arg555, Lys621, and Cys622 with a docking score of -11.923 kcal/mol and ligand efficiency of -0.163 followed by Ornipressin, showed a docking score of -11.717 kcal/mol (ligand efficiency of -0.163). Salt bridge with Asp623 and H-bond with Lys621 is observed during interactions however, Colistin showed two H-bonds with Ile548 and Ser814 with a docking score of -11.077 kcal/mol (-0.135). Colistin is also known as Polymyxin E (Supplementary Figure S11) . IFD poses of top 14 molecules are given in supplementary figures S12A and S12B. Pegamotecan, Cisatracurium, and colistin didn't show any pose during IFD docking simulation. Colistin showed similar results in both structures with and without co-factors. Based on docking analysis and its reconfirmation, we didn't find any significant difference in the docking energies, therefore our identified molecules have the potential to strongly bind in the binding tunnel of both single chain and complex from of the RdRp. The URL link of 3D structures, docking scores, ligand efficiencies, and Glide energies of additional top 100 molecules against RdRp complex are given in supplementary data. In a drug repurposing study, the aim is to screen existing approved drugs and suggest their reuse for new medical complications, being the most feasible approach as compared to the de-novo drug designing and development. The benefits associated with repurposed drugs include; safety for human use, no escalating cost, and reduced timeline for its development. For instance, Zidovudine, which was originally used for the treatment in cancer, was repurposed as the first anti-HIV drug, approved by the FDA in 1987 [50, 51] . Similarly, Rituximab, originally used in the treatment of various cancers was approved by the FDA in 2006 as a repurposed drug against rheumatoid arthritis [52] . Raloxifene, another drug used for osteoporosis, has also been approved for the treatment of breast cancer [52] . Aspirin, which was primarily used as an analgesic is now recommended and approved by the FDA in 2015 against colorectal cancer and in cardiovascular diseases [52] . In the last few years, about one-third of the FDA approved drugs relate to repurposed drugs [52, 53] . All of the screened molecules identified in this study are FDA approved drugs. We used fast and accurate docking algorithms Glide-SP to screen potential hits against RdRp (core and complex). Glide-IFD was used to compare and confirm the difference between the binding energies of top screened molecules. We have identified molecules for two forms of the RdRp individually. Some of the molecules have shown promising interactions with both forms (core and holoenzyme) of the RdRp, e.g. Polymyxin B1, Ornipressin, Lypressin, and Exemorelin. Although all the top ten screened molecules for both the forms showed strong binding. Polymyxin B1 shows strong binding with key amino acid residues in core as well as in the holoenzyme form. Polymyxin B1 is currently used as an antibiotic for the treatment of respiratory tract infections i.e. pneumonia (nosocomial, healthcare, or ventilator-associated), due to the multidrug-resistant gram-negative bacteria. For the treatment of respiratory tract infection, it is generally used as oral inhalation via nebulization. The general mechanism of action of Polymyxin involves the destabilization of the integrity of the bacterial membrane, thereby killing gram-negative bacterial cells through membrane lysis. It interacts with the lipopolysaccharides of the Gram-negative bacteria and disrupts the outer membrane. The other known mechanism is by vesicle-vesicle contact pathway, in which Polymyxin induces an osmotic imbalance, that ultimately leads to membrane lysis. Based on the findings of the current study, Polymyxin B1 can be used for the treatment of COVID-19 by blocking the active site of RdRp.

Due to the current COVID-19 pandemic, to date, 444,813 people have lost their lives. SARS-CoV-2A is spreading worldwide at an alarming rate. To overcome this infection, the current study was designed to screen all those drugs which are available in the market or currently in a clinical trial. For a potent drug, it is essential to interact with key amino acid residues in the active site pocket of the polymerase. Once the drug binds to the active sites of RdRp, the access of positive sense RNA of the virus which serves as a template for anti-genome is blocked. No new viral genomes are synthesized and mature virus particles cannot be assembled and propagated further. In this study, all the identified FDA approved drugs have shown stable interaction with key residues, along with lowest binding energies which indicates the potential for these drugs to inhibit the activity of RdRp, therefore contributing to stop viral proliferation. The top candidate drugs we found that interact with a single chain core RdRp include; Ornipressin, Otosiban, Lanreotide, Argiprestocin, Demoxytocin, Carbetocin, Lypressin, Examorelin, Colistin, and Polymyxin B1. Additionally, we also identified compounds that interact with the complex form of RdRp (holoenzyme). It includes; Nacartocin, Cistinexine, Cisatracurium, Pegamopecan, Ebiratide, Sulfomyxine, Diagastrin, Ditercalinium chloride, and Benzquercin. Among these lead molecules, Ornipressin, Lypressin, Examorelin, Polymyxin B1 showed strong binding efficacy with both core and holoenzyme. Although all of the screened molecules showed strong binding affinity towards RdRp (core and holoenzyme), however Polymyxin B1 and Colistin are already in use against respiratory tract infections. Based on our findings, we, therefore, suggest that these screened drugs can be clinically tested to handle this intricate infection.

The SARS-CoV-2 outbreak: what we know

Clinical Features of 69 Cases with Coronavirus Disease

Structure, Function, and Evolution of Coronavirus Spike Proteins

Genome Composition and Divergence of the Novel Coronavirus (2019-nCoV) Originating in China

Inferring the hosts of coronavirus using dual statistical models based on nucleotide composition

Structure of the SARS-CoV nsp12 polymerase bound to nsp7 and nsp8 co-factors

Biochemical characterization of a recombinant SARS coronavirus nsp12 RNA-dependent RNA polymerase capable of copying viral RNA templates

One severe acute respiratory syndrome coronavirus protein complex integrates processive RNA polymerase and exonuclease activities

RNA synthetic mechanisms employed by diverse families of RNA viruses

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

The mechanism of action of T7 DNA polymerase

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

Coronavirus puts drug repurposing on the fast track

Of chloroquine and COVID-19

Use of Hydroxychloroquine and Chloroquine During the COVID-19 Pandemic: What Every Clinician Should Know

Protein and ligand preparation: parameters, protocols, and influence on virtual screening enrichments

Very fast prediction and rationalization of pKa values for protein-ligand complexes

The OPLS [optimized potentials for liquid simulations] potential functions for proteins, energy minimizations for crystals of cyclic peptides and crambin

Epik: a software program for pK( a ) prediction and protonation state generation for drug-like molecules

Glide: a new approach for rapid, accurate docking and scoring. 1. Method and assessment of docking accuracy

Glide: a new approach for rapid, accurate docking and scoring. 2. Enrichment factors in database screening

Identification of Chymotrypsin-like Protease Inhibitors of SARS-CoV-2 Via Integrated Computational Approach

The pharmacology of ornipressin (POR-8): a local vasoconstrictor used in surgery

Ornipressin in the treatment of functional renal failure in decompensated liver cirrhosis: effects on renal hemodynamics and atrial natriuretic factor

Management of preterm labor: atosiban or nifedipine? International journal of women's health

The safety of lanreotide for neuroendocrine tumor

Vasopressin and its analogues in shock states: a review

Use of oral oxytocics for stimulation of labor in cases of premature rupture of the membranes at term: a randomized comparative study of prostaglandin E2 tablets and demoxytocin resoriblets

Carbetocin versus oxytocin in caesarean section with high risk of post-partum haemorrhage

Usefulness in patients who manifest allergies to other antidiuretic hormone preparations

Growth hormone-releasing activity of hexarelin in humans. A dose-response study

Effects of GHRP-2 and hexarelin, two synthetic GH-releasing peptides, on GH, prolactin, ACTH and cortisol levels in man. Comparison with the effects of GHRH, TRH and hCRH

The clinical use of colistin in patients with cystic fibrosis

Effect of aerosolized colistin on multidrug-resistant Pseudomonas aeruginosa in bronchial secretions of patients without cystic fibrosis

Combination therapy with intravenous colistin for management of infections due to multidrug-resistant Gram-negative bacteria in patients without cystic fibrosis

Treatment of multidrug-resistant Acinetobacter baumannii ventilator-associated pneumonia (VAP) with intravenous colistin: a comparison with imipenem-susceptible VAP

Aerosolized colistin for the treatment of nosocomial pneumonia due to multidrug-resistant Gram-negative bacteria in patients without cystic fibrosis

Colistin: the revival of polymyxins for the management of multidrug-resistant gram-negative bacterial infections

Polymyxin B for the treatment of multidrug-resistant pathogens: a critical review

Nacartocin--analogue of oxytocin with enhanced natriuretic properties: natriuretic and hemodynamic characteristics

The effect of a new expectorant drug on mucus transport in chronic bronchitis

Use of cisatracurium in critical care: a review of the literature

A phase II study of pegylated-camptothecin (pegamotecan) in the treatment of locally advanced and metastatic gastric and gastro-oesophageal junction adenocarcinoma

The neurotrophic effects of ebiratide, an analog of ACTH4-9, on cultured septal cells and aged rats

Sulfomycins, a series of new sulfur-containing antibiotics. I. Isolation, purification and properties

Effect of somatostatin analogs on gastric acid secretion in dogs and rats

Ditercalinium chloride, a pro-anticancer drug, intimately associates with mammalian mitochondrial DNA and inhibits its replication

In-silico identification and evaluation of plant flavonoids as dengue NS2B/NS3 protease inhibitors using molecular docking and simulation approach

Azido-3'-deoxythymidine (BW A509U): an antiviral agent that inhibits the infectivity and cytopathic effect of human T-lymphotropic virus type III/lymphadenopathy-associated virus in vitro

Administration of 3'-azido-3'-deoxythymidine, an inhibitor of HTLV-III/LAV replication, to patients with AIDS or AIDS-related complex

Drug repurposing: progress, challenges and recommendations

Challenges and opportunities with drug repurposing: finding strategies to find alternative uses of therapeutics

We are thankful to Shelvia Malik and other Schrodinger team members for facilitating us with the licence to use Schrodinger Suite which enabled us to successfully complete this study.