Keywords
COVID-19, SARS-CoV-2, Docking, In silico, RNA-dependent RNA polymerase, Molecular dynamics simulation, Drug repurposing
This article is included in the Manipal Academy of Higher Education gateway.
This article is included in the Emerging Diseases and Outbreaks gateway.
This article is included in the Coronavirus collection.
COVID-19, SARS-CoV-2, Docking, In silico, RNA-dependent RNA polymerase, Molecular dynamics simulation, Drug repurposing
A novel disease was first discovered in late 2019 in a 41-year-old patient in China who was admitted to a health facility for severe respiratory syndrome and fever, dizziness, and cough. A new strain of RNA virus of the Coronaviridae family was isolated from the patient's bronchoalveolar lavage fluid and was originally unique initially referred to as 'WH-Human 1' coronavirus (and later known as '2019-nCoV') was. The novel virus was with 60–140 nm diameter with 9–12 nm long spikes, with the virion comparable to the corona of the sun. Viral genome analysis revealed that it is made of ~30k nucleotides, and has 89.1% nucleotide similarity to SARS-like CoV (genus Betacoronavirus, subgenus Sarbecovirus)1. The human CoV B814-strain was isolated in 1965 from the nasal discharge of a patient with a common cold. Since then, 30 specific strains were identified. Seven hCoVs have been reported previously causes disease, including HCoV-229E, HCoV-OC43, HCoV-NL63, HCoV-HKU1, SARS-CoV, MERS-CoV, and 2019-nCoV (now established as SARS-CoV-2)2. SARS-CoV emerged during the period 2000–2004, and it was discovered that the infection was descended from animals and bats as intermediate hosts3. MERS-CoV was diagnosed for the first time in Saudi Arabia4. The World Health Organization (WHO) named the disease coronavirus disease 2019 (COVID-19) on 12 January 2020. The same outbreak was declared a public health emergency of international concern by the WHO on 30-January-2020 due to rapid transmission mortality rate5. The International Committee on the Taxonomy of Viruses Coronavirus Study Group (CSG) proposed the name of nCoV as SARS-CoV-26.
Currently, clinical trials have approved the use of remdesivir, favipiravir, and dexamethasone, but the infection could not be cured or prevented. The remdesivir and favipiravir are comparatively expensive than dexamethasone and not affordable for people at developing countries. Vaccine development, convalescent plasma, interferon-based therapies, small molecule medications, cell-based therapy, and monoclonal antibodies (mAbs) are the various likely future medications and pathways being investigated7. The development of a drug is expensive and time-consuming with a high attrition rate, which is unacceptable in the current global emergency. Hence, there is a great interest in repurposing the existing drugs and speed-up to develop antiviral therapies8. The national health commission (NHC) of the People's Republic of China recommends the use of α-interferon (IFN-α), lopinavir/ritonavir, ribavirin, chloroquine phosphate, and arbidol as antiviral therapy against SARS-CoV-2. Other agents such as inhibition of fusion/entry (camostat mesylate, baricitinib, arbidol and chloroquine phosphate), the agents that disrupt SARS-CoV-2 replication (remdesivir, favipiravir, lopinavir/ritonavir), the agents that suppress excessive inflammatory response (corticosteroid), convalescent plasma, vaccines (inactivated vaccine, recombinant subunits vaccine, nucleic acid-based vaccine, adenoviral vector vaccine, recombinant influenza viral vector vaccine), a combination of traditional Chinese and Western medicine (ShuFeng JieDu capsules and Lianhua Qingwen capsules) are also gaining therapeutic interest9.
The SARS-CoV-2 genetic material is translated inside the host cell to two distinct groups of proteins; structural proteins, such as Spike (S), Nucleocapsid (N), Matrix (M) and Envelope (E), and non-structural proteins (nsp) such as RNA dependent RNA polymerase (nsp12), helicase (nsp13), papain-like protease (PLpro or nsp3) and main protease Mpro (also known nsp5 or 3CLpro)10–12. SARS-CoV-2 replication is facilitated through a multi-subunit replication/transcription complex of non-structural viral proteins. The key aspect of the nsp complex is an RNA-dependent RNA polymerase (RdRp; nsp12). The central RdRp domain is divided into three subdomains, the thumb, palm, and right-handed cup-like fingers, which is less than 500 long amino acids length11. A high level of homology has been found between SARS-CoV-2, MERS-CoV, and SARS-CoV RdRp13. The RdRp requires extra factors such as nsp7 and nsp8, for its activity14. Remdesivir, which was recently approved by US-FDA for COVID-19, inhibits RdRp15. Furthermore, researchers have explored RdRp-target for in silico repurposing of drugs. Hence, this proves that Nsp12 or RdRp is an attractive target for inhibiting viral replication. In the current study, we assessed USFDA approved drugs to impede the RdRp activity of SARS-CoV-2 through an in silico drug-repurposing strategy.
All computational simulation experiments were conducted on the Schrödinger Suite's (version 2020-2) Maestro graphical user interface (www.schrodinger.com; v12.3) on an Ubuntu software desktop workstation, Intel ® Xenon ® Gold 6130 CPU @ 2.10 GHz × 64 processors, Quadro P620 / PCle / SSE2 graphics card and 134.8 GB RAM. Alternatively, free software including Autodock 4.2.6, Gromacs 5.1, ArgusLab 4.0 and PDB2PQR 3.0 could be used for such study.
Molecules licensed by US-FDA were downloaded from https://www.drugbank.ca (date of access: 15 Feb 2020). The molecules were configured using LigPrep, wherein the molecules generated the 3D coordinates16. Using the Epik module, the suitable ionization state was predicted at pH 8.0, types of tautomers were produced, and proper chirality was estimated. Finally, the molecules' structure was minimized with energy using the OPLS3e force field17. Alternatively, the Zinc database can be used for FDA approved drug and ArgusLab 4.0, PDB2PQR can be used for ligand configuration and stabilising the ligand state at certain pH respectively.
The electron microscopical RdRp coordinate structure of SARS-CoV-2 linked to template-primer RNA which was 50- base pair long and remdesivir triphosphate was downloaded from Protein Data Bank (PDB), accession code 7BV2, at a resolution of 2.5 Å. The RdRp-NS7-NS8 complex was linked to three Mg2+ ions (coordinated by Asp760, Asp761, Asp618, and pyrophosphate), double-stranded RNA (14-base template strand and 11-base primary strand) and monophosphate remdesivir (RMP). The Protein Prep Wizard was used for the optimization of protein structure18. ArgusLAb 4.0 software, which is free tool for protein preparation for docking studies can be used for the same. In this process, the missing hydrogens, side chains, and other organic solvent molecules and water residues were eliminated. The proper ionization state was then created at pH 8.0 using PROPKA3, regenerating the hydrogen bond network, and finally minimizing the protein structure using a restrained minimization technique19,20.
Ligand docking of molecules licensed by US-FDA (~2800 Nos) was performed using the Schrodinger Glide v8.7 module21. Alternatively, AutoDock 4.2.6, an open access package for ligand docking, can be used for the same function. The binding site of the ligands on the protein was assessed using the Receptor Grid Generation tool using the centroid formed around the bound ligand. The docking was initially performed using high-throughput virtual screening (HTVS). Later, top molecules with strong binding to the active site of RdRp and docking scores were moved ahead for standard precision (SP)-docking. The single best pose of each molecule was saved during ligand docking22.
The MD-simulation was run on Schrodinger's Desmond module23. Alternatively, the free software NAMD with VMD molecular simulation could be used. The solvated water-soaked system was generated using the Desmond System Builder tool. The solvating system used in the experiment was the TIP3P model. An orthorhombic simulation was a box that generated at least 10 Å from the outer surface of the protein with periodic boundary conditions. The system was neutralized by adding a reasonable amount of counterions. By adding 0.15 M NaCl into the simulation panel, the isosmotic condition was conserved. A pre-defined equilibration procedure was performed before the simulation. The MD simulation was performed at an ambient pressure of 1.013 bar, and a temperature of 300°K, with 1000 frames saved to the trajectory, for 50 nsec period. Later the simulation was analyzed using a simulation interaction diagram, wherein protein-ligand root mean square deviation (RMSD), protein root mean square fluctuations (RMSF), ligand RMSF, protein-ligand contacts, ligand-protein contacts, ligand torsion profile. were analyzed.
SARS-CoV-2 RdRp is a central polymerase involved in the replication of RNA, is a big enzyme consisting of 932 amino acids. The amino acid framework of SARS-CoV RdRp linked to the Nsp7 and Nsp8 cofactors24. Structurally, the RdRp protein is divided into the N-terminal and polymerase domains. The N-terminal domain extends from amino acid residues 1 to 397. The polymerase domain extends from amino acids 398 to 919 and is similar to the previously reported SARS-CoV nsp12 (PDB ID: 6NUR). The RdRp polymerase domain is subdivided into three structurally different subunits, namely the finger subunit (398–581 and 628–687 amino acids), palm subunit (582–627 and 688–815 amino acids), and thumb subunit (816–919 amino acids), which embraces a "right-handed" shaped conformation15. The active site is located at the hinge interface between the finger and thumb subunits of the SARS-CoV-2 RdRp, which is an evolutionary conserved in both the coronaviruses reported earlier, ie. MERS and SARS-nsp12. The active site of nsp12 is situated in the middle of the substrate domain where the synthesis of RNA takes place as an RNA template is accessed from the template input channel and nucleoside triphosphate (NTP) from NTP input channel25. Structural assessment of SARS-CoV-2 defined in RdRp displayed the presence of seven retained motifs (A-G motifs), which play a key role in the association of substrate and nucleotides. They also play a key role in catalysis of the attachment of the NTP. Antiviral drugs function by binding to first replicated base pairs and thereby stopping the chain elongation cycle. The sensitivity of the RdRp binding to the RNA prototype has been enhanced by the interaction of the nsp7 and nsp8 co-factors with RdRp. The RNA polymerization activity appeared to be impaired in the presence of remdesivir active triphosphate (RTP), which is being investigated in clinical trials. Remdesivir blocked the polymerization process of RdRp at 1 mM and 10 mM ATP completely. However, remdesivir alone on its own or in its monophosphate form did not demonstrate activity at 5 mM concentration. Remdesivir monophosphate was found to be covalently linked to the primary strand in protein structure to the pyrophosphate moiety and three magnesium as catalytic ions. The ligand-protein interaction of remdesivir also has two hydrogen bonds, with the uridine base, which leads to a stable complex formation. Predominant interactions include those with K545 and R55 side chains, the interaction of magnesium ions with the backbone of the phosphate diester, as well as with the primary strand RNA. Seven conserved motifs from A to G represent the nsp12 RdRp catalytic active site, among which the palm subdomain with an SDD sequence (residues 759-761) in motif C, including amino acid residues D760 and D761 in the coordination of the two magnesium ions. The template-primer RNA, which oriented in the same fashion as template-primer RNA in the poliovirus RdRp elongation complex and residues involved in binding as well as interactions to 2′-OH groups, was found to be highly conserved15.
The docking provides details about binding affinity and orientation of the interactions between ligand and protein. HTVS was initially performed using ~2800 drug molecules. The HTVS docking is used to screen large numbers of ligands rapidly. However, HTVS seems to have more limited conformational sampling than docking with SP, and hence HTVS can not be used with ranking in place. Hence, molecules were advanced to the higher precision of docking known as SP. The SP docking was performed on 500 molecules, which exhibited good interactions with the RdRp and had topped in the docking score. Among the 500 molecules, 14 displaying top SP-docking scores and the conserved interactions with the active site as that of remdesivir were shortlisted based on visual inspection (Table 1). These molecules were able to interact with RdRp amino acid residues via a predominant metal coordination bond and hydrogen bonding with the active site. The interactions exhibited by rosoxacin, kappadione, tedizolid phosphate, levomefolic acid, etodolac, pitavastatin, montelukast, norfloxacin, cinoxacin, mefenamic acid, diflunisal, fospropofol, fluvastatin, and ridogrel was found to be prominent. Figure 1 depicts the binding mode of Pitavastatin, Ridogrel and Rosoxacin, which emerged as the best molecule by MD simulations.
The mode of interaction of selected molecules with the RdRp active site in the SP-docking is depicted in Table 1. The shortlisted molecules interact with conserved residue ARG555 and divalent magnesium in the manner Remdesivir interacted with RdRp. Drugs including rosoxacin, levomefolic acid, etodolac, kappadione, pitavastatin, montelukast, norfloxacin, cinoxacin, diflunisal, fospropofol and ridogrel have demonstrated π-cation interactions with ARG555. Norfloxacin and cinoxacin have demonstrated π-π stacking interactions with 20 unspecified residues: (U P: 20). In SP-docking, the drugs shortlisted made metal coordination and salt bridges with MG A:1004 and MG A:1005 (divalent magnesium ions) in the same way remdesivir was bound (Table 1). Rosoxacin, pitavastatin and ridogrel had strong interactions across metal coordination bonds among the molecules shortlisted by SP-docking, and π-cation with RdRp can be further tested for their in vitro function.
The RdRp-protein and ligand binding stability were checked based on the RMSD fluctuations during MD-simulations. In the MD simulation, the RMSD fluctuation was measured individually for the protein and ligand structures, and it falls within 3 Å. Hence the complex was considered to be stable. Besides, the persistence of the intermolecular interactions between ligand and RdRp was also monitored during the simulation. The analysis RMSD plot for the structures saved in the trajectory generated during the MD-simulation exhibited stable bindings for 14 drug molecules. For rosaxocin, pitavastatin and ridogrel in complex with the RdRp-protein, the RMSD fluctuations for the ligand and protein remained within 2.0 Å. Figure 2 depicts the ligand and protein RMSD plot for the shortlisted 14 drugs.
During the initial phase of MD simulations, the RdRp-bound pitavastatin complex exhibited higher fluctuations due to the equilibration. For the initial 10 nsec of the simulation, the protein backbone exhibited higher RMSD fluctuation. The latter 40 nsec the protein backbone fluctuations remained within the range of 0.6 Å, indicating stabilization of the protein structure. Similarly, the ligand showed consistently stable interactions throughout the simulation. The ligand RMSD fluctuations remained within the range of 0.5 Å for the remaining 40 nsec of the simulation indicated a stable ligand in the complex (Figure 2). The interaction with residues ASP618, ASP760, ASP761, and GLU811 exhibited maximum contacts with the ligand during MD simulation. The bridged hydrogen bonding was exhibited with the ligand was amino acid residue TYR619 over more than 30.0% of the simulation time in the 50 nsec trajectory. A significant metal coordination bond was formed with the divalent magnesium and residues ASP618, ASP760, ASP761, and GLU811. The majority of the interactions observed between the protein and ligand during the docking were consistently retained during the MD simulation (Figure 3). Pitavastatin is used as a primary hyperlipidemic agent and in mixed dyslipidemia. It lowers the elevated total cholesterol, low-density lipoprotein, apolipoprotein-B, triglycerides, and increases high-density lipoprotein. Statins are commonly prescribed in cardiovascular patients. Statins offer protection to vasculatures and are known to have pleiotropic effects in the body, especially in modulating the inflammatory process at the cellular level26. One of the reports on in silico studies analysed the statins interaction with Mpro, and they have also reported the pitavastatin as a hit-drug27. These combined Mpro and RdRp inhibitory properties of statins might help to tackle COVID-19.
The RdRp bound ridogrel complex exhibited a mixture of hydrophilic and ionic interactions during the MD simulation. For the initial ten nsec of the simulation, the protein backbone exhibited higher RMSD fluctuation. The latter 40 nsec the protein backbone fluctuations remained within the range of 0.6 Å, indicating stabilization of the protein structure (Figure 2). Similarly, the ligand showed a higher fluctuation for the initial 15 nsec of the simulation. The ligand RMSD fluctuations remained within the range of 1.5 Å for the remaining 35 nsec of the simulation, indicating a relatively stable ligand-protein complex. After the initial fluctuation for a period of 10 nsec due to the equilibration, the RMSD for the protein structures remained between 2.6-4.4 Å until the end of the simulation. Amino acid residues ASP760 followed by CYS618, ASP622 and ASP623 interacted to the greatest extent with the ligand. The residues CYS622 and ASP623 made bridged hydrogen bonding interaction with the ligand which accounted to have 75 and 74% interaction with the protein residues, respectively. The amino acid residues ASP618 and ASP760 made metal coordination bonds with divalent magnesium. The majority of the interactions observed between the protein and ligand during the docking were consistently retained during the MD simulation (Figure 4). Based on the RdRp interaction, we propose that ridogrel can be repurposed for COVID-19 as an antiviral agent. Ridogrel is a thromboxane synthase inhibitor and thromboxane or prostaglandin/endoperoxide receptor blocker. It is used for the prevention of systemic thromboembolism in acute myocardial infarction28. Redogrel is been replaced with aspirin in the therapy as aspirin was found to be clinical superior. The RdRp inhibitory properties of ridogrel might make it a useful drug in COVID-19. Hence, it can be recommended for further screening and optimize it by in vitro assays.
RdRp bound with rosaxocin exhibited a mixture of hydrophobic interactions and hydrophilic interactions during the MD-simulation. For the initial 20 nsec of the simulation, the protein backbone exhibited higher RMSD fluctuation. The latter 30 nsec the protein backbone fluctuations remained within the range of 0.8 Å, indicating stabilization of the protein structure (Figure 2). Similarly, the ligand showed a higher fluctuation for the initial 40 nsec of the simulation. The ligand RMSD fluctuations remained within the range of 1.4 Å for the remaining 30 nsec of the simulation indicated a stable ligand-protein complex. The residues ASP760, ARG555 and ASP623 made maximum contact with the ligand. A predominant π-cation interaction was exhibited with ARG555 residue (54%). A direct salt bridge interaction was observed with ARG555 with the carbonyl carbon, and a bridged hydrogen bonding interactions, which accounted for 92% was observed. Similarly, bridged hydrogen bonding interactions were observed with ARG553 and carboxyl electronegative oxygen (accounted for 64%). The amino acid residues LYS621, CYS622 made direct hydrogen-bonding interactions with carboxyl electronegative oxygen (accounted for 53%, 60%) and ASP623 made direct hydrogen bonding with unsaturated oxygen located in the ligand (58%)respectively over more than 30.0% of the simulation time in the 50 nsec simulation trajectory. Additionally, metal coordination interaction was observed with ASP760 with the divalent magnesium. The majority of the interactions observed between the protein and ligand during the docking were consistently retained during the MD simulation (Figure 5). Rosoxacin is a broad-spectrum quinolone antibiotic for the treatment of urinary tract infections (UTIs) and sexually transmitted diseases29. The popular anti-COVID-19 drugs chloroquine and hydroxychloroquine have the quinoline moiety30. Considering the PK/PD and the safety profiles of rosoxacin, it can be taken forward in clinical trials as a repurposed candidate for COVID-19. This study is the first report of a broad-spectrum antibiotic with a high affinity for RdRp identified by in silico tools.
The other molecules short listed are from the diverse pharmacological group can be repurposed by studying in vitro. These molecules include etodolac, montelukast, fluvastatin, kappadione, levomefolic acid, norfloxacin, cinoxacin, tedizolid phosphate, mefenamic acid, diflunisal and fospropofol. The details of these MD-simulation data and molecular interaction are represented in Figure 6–Fiugre 16.
Etodolac is an NSAID with anti-inflammatory, analgesic and antipyretic properties. It inhibits the cyclooxygenase (COX) and prevents the formation of peripheral prostaglandins which causes inflammation. It is officially prescribed for rheumatoid arthritis and osteoarthritis. Etodolac is 50-times more selective for COX-2 than COX-1, and the antipyresis is attributed to hypothalamic actions, cutaneous blood flow, and subsequent heat losses31. Hence, if this drug is repurposed for COVID-19, it will have a multidirectional approach because the drugs will be useful in bringing down the elevated body temperature and anti-inflammatory properties. Montelukast, a leukotriene receptor antagonist, used to treat asthma, exercise-induced bronchoconstriction, and allergic rhinitis32. The anti-inflammatory properties, as well as binding of montelukast to RdRp, might help in tackling COVID-19. Hence, montelukast can be tried for drug repurposing with detailed studies. Fluvastatin is a hypolipidemic agent belongs to the class statins33. Because of their pleiotropic role, statins can be explored further to repurpose to COVID-19. Kappadione or menadiol sodium diphosphate is a highly water-soluble vitamin K analogue approved by the FDA and marketed by Lilly. It is indicated in anticoagulant-induced prothrombin deficiency, and therapy of hypoprothrombinemia due to antibacterial therapy, hypoprothrombinemia secondary to factors limiting absorption or synthesis of vitamin K. Recently, a computational study reported that kappadione would tightly bind to PAX2 transcription factor which regulates the ABC-transporter. PAX2 binding by kappadione inhibits the PAX2-DNA interaction, and this would be beneficial for combating chemoresistance in pancreatic ductal adenocarcinoma34.
Levomefolic acid, also known as 5-methyltetrahydrofolate (5-MTHF) is an active form of folate found in food, which, when administered, increases folate (vitamin B9) concentration in blood. Levomefolic acid is preferred over folic acid supplement due to lower potential for the induction of vitamin B12 deficiency symptoms35. Hence, this drug is safest can be easily repurposed in COVID-19 patients. Norfloxacin is a broad-spectrum fluoroquinolone antibiotic. Norfloxacin acts by blocking DNA gyrase and used to treat UTI36. Similar to quinolones, fluoroquinolones also found to modulate the cytokines, and this property is an added advantage in repurposing norfloxacin for COVID-19. Cinoxacin is currently used to treat UTIs and as an alternative antimicrobial to oxolinic acid and nalidixic acid37. Tedizolid is an oxazolidinone class of antibiotics, a congener of linezolid, which is effective against multidrug-resistant Gram-positive bacteria. Tedizolid is useful in acute bacterial skin infections and was approved by the FDA in 2014. Currently available as both an oral tablet and as a powder for intravenous injection38. Mefenamic acid belongs to the class of aminobenzoic acids effective for the treatment of dysmenorrhea, rheumatoid arthritis, osteoarthritis, mild to moderate pain, inflammation, and fever39. Using this drug will have advantages in COVID-19, as it also tackles the inflammation and has the antipyretic effect. Diflunisal, a salicylate derivative has anti-inflammatory, analgesic, and antipyretic function as NSAID40. It modulates HMGB1 through Toll-like receptor 4 (TLR4). Diflunisal does not inhibit responses that rely on TLR4. The TLR modulation is one among the targets for COVID-19, diflunisal has advantages. Fospropopol is a hypnotic/sedative/anaesthetic short-acting drug prodrug of propofol, which is an anaesthetic drug. It has a high lipid solubility can be tried in COVID-19 based on its interactions with RdRp41.
Molecular docking of FDA-approved drugs for RdRp led us to shortlist 14 molecules. There were many proposals of RdRp inhibitors for repurposing, especially by in silico computational simulations. One previous report used the ZINC drug database to screen SARS-CoV-2 RdRp inhibitors13. They studied 78 commonly used antiviral drugs and drugs which are currently in clinical trials for COVID-19 for molecular simulations and reported the possible in silico repurposing of drugs. They have listed the top 20 drugs from the ZINC-drug database that included the currently used antiviral drug valganciclovir, anti-bacterial drugs such as chlorhexidine, ceftibuten, cefuroxime, novobiocin, the anti-asthmatic drugs fenoterol and cromolyn, the antitumour drugs such as fludarabine, idarubicin, the antimalarial drug atovaquone, the muscle relaxant pancuronium,the anti-allergic cortisone, the contraceptive tibolone, the hepatoprotective silybin, the anti-diarrheal drug diphenoxylate, the anti-amenorrheal bromocriptine, and chenodeoxycholic acid, which is used to dissolve gallstones. In the same study, the investigators shortlisted the top-20 natural compounds as RdRp inhibitors. Another group of researchers studied the feasibility of known RNA-polymerase inhibitors to repurpose as anti-SARS-CoV-2 drugs42. The drugs they studied are broad-spectrum antiviral agents such as remdesivir, 5-fluorouracil, ribavirin and favipiravir. However, they noted that the SARS-CoV-2 might evolve to acquire drug resistance mutations against these nucleoside inhibitors. The virtually combined deep learning and molecular docking simulations approach identified potentially useful 49 most promising FDA-approved drugs for repurposing in COVID-1943. The selected drugs in this study were anidulafungin, velpatasvir, glecaprevir, rifabutin, procaine penicillin G, tadalafil, riboflavin 5'-monophosphate, flavin adenine dinucleotide, terlipressin, desmopressin, elbasvir, oxatomide, enasidenib, edoxaban and selinexor. Another group of computational chemists studied the FDA-approved drugs by virtual screening for the three proteases of SARS-CoV-2, Mpro, PLpro and RdRp44. They have shortlisted promising antiviral drugs such as simeprevir, ledipasvir, idarubicin, saquinavir, ledipasivir, partitaprevir, glecaprevir, and velpatasvir. They also shortlisted several novel drugs such as antiviral raltegravir; antidiabetic amaryl, antibiotics retapamulin, rifimixin, and rifabutin; antiemetic fosaprepitant and netupitant. Another group of researchers used a multi-targetted approach and shortlisted molecules such as δ-viniferin, myricitrin, taiwanhomoflavone A, lactucopicrin 15-oxalate, nympholide A, afzelin, biorobin, hesperidin and phyllaemblicin B. These molecules have the affinity for SARS-CoV-2 Mpro, RdRp and hACE245. Furthermore, a group of researcher have proposed well-tolerated, cost-effective drugs to combat COVID-19 which has the affinity for both Mpro and RdRp. The compounds they have shortlisted are ergotamine, dihydroergotamine, conivaptan, paliperidone, and tipranavir46.
In the pathogenesis and pathophysiology of COVID-19, both innate and adaptive immunity are activated, and many inflammatory mediators are released47. Mitigating inflammation will reduce the lethality because death due to COVID-19 is partly due to cytokine storm and end-organ failure48–50. The immuno-inflammatory modulators have their role in the treatment of COVID-19, and the drug shortlisted in this study have those potentials have potentials to mitigate inflammation. The shortlisted 14 novel molecules for SARS-CoV-2 RdRp inhibitors from different pharmacological classes. These drugs include antimicrobials such as rosoxacin, norfloxacin, cinoxacin, tedizolid phosphate, vitamin analogues such as kappadione, levomefolic acid, anti-inflammatory agents such as etodolac, mefenamic acid, diflunisal, montelukast, statins such as pitavastatin, fluvastatin, a sedative fospropofol, and a thromboxane inhibitor ridogrel. These drugs are commonly used, well-tolerated and economical. Among the 14-drugs, pitavastatin, is a lipid lowering agent and has pleotropic effects. Ridogrel rosoxacin, a quinolone antibiotic indicated for the treatment of urinary tract infections and certain sexually transmitted diseases, Testing in vitro for RdRp inhibitory activity, and antiviral activity could confirm these molecules for the repurposing or repositioning for COVID-19 treatment.
In the absence of approved therapies for treatment or prevention, drug repurposing has provided valuable insight into the treatment of COVID-19. The antiviral drugs such as lopinavir, ritonavir, favipiravir and remdesivir, immunosuppressants such as sirolimus, anthelmintics like ivermectin, corticosteroids such as methylprednisolone, dexamethasone are currently under clinical trials for COVID-19. In this molecular docking study using structure-based virtual screening, we identified favourable drugs, namely rosoxacin, levomefolic acid, etodolac, kappadione, pitavastatin, montelukast, fluvastatin, norfloxacin, cinoxacin, tedizolid phosphate, mefenamic acid, diflunisal, fospropofol, and ridogrel. The MD simulation studies revealed pitavastatin, ridogrel and rosoxacin to be superior compared to other shortlisted drugs. Further, the therapeutic properties of these shortlisted molecules make it suitable for testing RdRp-inhibitory activity and antiviral activity in vitro and in vivo could confirm these molecules for repurposing in COVID-19.
All data underlying the results are available as part of the article and no additional source data are required.
Authors are thankful to the Department of Pharmaceutics, Manipal College of Pharmaceutical Sciences for facilities of Schrodinger’s Suit computer simulations.
Views | Downloads | |
---|---|---|
F1000Research | - | - |
PubMed Central
Data from PMC are received and updated monthly.
|
- | - |
Is the work clearly and accurately presented and does it cite the current literature?
Yes
Is the study design appropriate and is the work technically sound?
Yes
Are sufficient details of methods and analysis provided to allow replication by others?
Yes
If applicable, is the statistical analysis and its interpretation appropriate?
Not applicable
Are all the source data underlying the results available to ensure full reproducibility?
Yes
Are the conclusions drawn adequately supported by the results?
Yes
Competing Interests: No competing interests were disclosed.
Reviewer Expertise: Inflammation, computational biology, neurodegeneration, aging and Alzheimer's disease.
Is the work clearly and accurately presented and does it cite the current literature?
Yes
Is the study design appropriate and is the work technically sound?
Yes
Are sufficient details of methods and analysis provided to allow replication by others?
Yes
If applicable, is the statistical analysis and its interpretation appropriate?
I cannot comment. A qualified statistician is required.
Are all the source data underlying the results available to ensure full reproducibility?
Yes
Are the conclusions drawn adequately supported by the results?
Yes
Competing Interests: No competing interests were disclosed.
Reviewer Expertise: Synthetic Organic Chemistry, Medicinal chemistry, Natural Product Chemistry.
Alongside their report, reviewers assign a status to the article:
Invited Reviewers | ||
---|---|---|
1 | 2 | |
Version 1 23 Sep 20 |
read | read |
Provide sufficient details of any financial or non-financial competing interests to enable users to assess whether your comments might lead a reasonable person to question your impartiality. Consider the following examples, but note that this is not an exhaustive list:
Sign up for content alerts and receive a weekly or monthly email with all newly published articles
Already registered? Sign in
The email address should be the one you originally registered with F1000.
You registered with F1000 via Google, so we cannot reset your password.
To sign in, please click here.
If you still need help with your Google account password, please click here.
You registered with F1000 via Facebook, so we cannot reset your password.
To sign in, please click here.
If you still need help with your Facebook account password, please click here.
If your email address is registered with us, we will email you instructions to reset your password.
If you think you should have received this email but it has not arrived, please check your spam filters and/or contact for further assistance.
Comments on this article Comments (0)