- 1Department of Biotechnology, School of Biological Sciences, Science Campus, Alagappa University, Karaikudi, India
- 2Laboratory of Bioinformatics and Computational Biology, Department of Genetics and Molecular Biology, Federal University of Rio de Janeiro State (UNIRIO), Rio de Janeiro, Brazil
The novel coronavirus SARS-CoV-2 disease “COVID-19” emerged in China and rapidly spread to other countries; due to its rapid worldwide spread, the WHO has declared this as a global emergency. As there is no specific treatment prescribed to treat COVID-19, the seeking of suitable therapeutics among existing drugs seems valuable. The structure availability of coronavirus macromolecules has encouraged the finding of conceivable anti-SARS-CoV-2 therapeutics through in silico analysis. The results reveal that quinoline,1,2,3,4-tetrahydro-1-[(2-phenylcyclopropyl)sulfonyl]-trans-(8CI) and saquinavir strongly interact with the active site (Cys-His catalytic dyad), thereby are predicted to hinder the activity of SARS-CoV-2 3CLpro. Out of 113 quinoline-drugs, elvitegravir and oxolinic acid are able to interact with the NTP entry-channel and thus interfere with the RNA-directed 5′-3′ polymerase activity of SARS-CoV-2 RdRp. The bioactivity-prediction results also validate the outcome of the docking study. Moreover, as SARS-CoV-2 Spike-glycoprotein uses human ACE2-receptor for viral entry, targeting the Spike-RBD-ACE2 has been viewed as a promising strategy to control the infection. The result shows rilapladib is the only quinoline that can interrupt the Spike-RBD-ACE2 complex. In conclusion, owing to their ability to target functional macromolecules of SARS-CoV-2, along with positive ADMET properties, quinoline,1,2,3,4-tetrahydro-1-[(2-phenylcyclopropyl)sulfonyl]-trans-(8CI), saquinavir, elvitegravir, oxolinic acid, and rilapladib are suggested for the treatment of COVID-19.
Introduction
The current outbreak of coronavirus disease 2019 (COVID-19), caused by the severe acute syndrome coronavirus-2 (SARS-CoV-2), has been considered as a major anxiety of the twenty first century (Dey et al., 2020). As of May 07, 2020, WHO states that over 3,672,238 cases have been authoritatively affirmed, including 254,045 deaths around the globe. The pathognomonic symptoms of COVID-19 are fever, dry cough, shortness of breath, and dyspnea (Wu et al., 2020). In severe conditions, it causes hypercytokinemia, lymphopenia, disseminated intravascular coagulation, severe acute respiratory syndrome, kidney failure, and eventually death (Thomas-Rüddel et al., 2020). In general, SARS-CoV-2 is a positive sense, long (30,000 bp), single-strand RNA coronavirus that belongs to the family Coronaviridae and genus Betacoronavirus, which is highly similar to SARS-CoV (Zou et al., 2020). No specific medication for COVID-19 is accessible at the present time. Thus, researchers are seriously searching for suitable vaccines and therapeutic-drugs against COVID-19 (Yao et al., 2020). The fact is that discovery, as well as marketing, of new drugs frequently takes months to years (Stebbing et al., 2020), and hence looking for appropriate therapeutics among existing-drugs seems to be a promising strategy to control the current pandemic of COVID-19 in this critical time.
In the SARS-CoV-2 macromolecules, the large-polyproteins-encoded cysteine protease, called 3-chymotrypsin like protease [3CLpro or main protease (Mpro)], are essential for the viral life-cycle of novel coronavirus (Zhang H. et al., 2020). This enzyme plays a crucial role in the processing of viral polyproteins, which are indispensable for viral maturation and their infectivity (Khan et al., 2020). Subsequently, RNA-dependent RNA polymerase (RdRp) is a key enzyme essential for the viral replication of SARS-CoV-2 (Gao et al., 2020). Due to their crucial roles, these viral proteins are considered as imperative targets for developing antiviral compounds against COVID-19 (Wu et al., 2020). Recently, Choy et al. (2020) reported that the combination of HIV-protease inhibitors such as lopinavir/ritonavir effectively kills SARS-CoV-2 at the cellular level. Similarly, Wang M. et al. (2020) reported that the nucleotide analog RdRp-inhibitor, remdesivir, successfully inhibited SARS-CoV-2 in vitro. Hall and Ji (2020) reported zanamivir, indinavir, saquinavir, and remdesivir as SARS-CoV-2 3CLpro inhibitors using in silico analysis. Further, Elfiky (2020) also suggested ribavirin, remdesivir, sofosbuvir, galidesivir, and tenofovir as potent drugs against SARS-CoV-2 through docking analysis.
On the other hand, human angiotensin-converting enzyme 2 (ACE-2), a type-I integral membrane protein, has been considered to be the specific and functional receptor for the spike glycoprotein of SARS-CoV-2 (Patel et al., 2014). It is also well-known to play the main role in the rennin-angiotensin system (RAS), which is associated with the regulation of heart function and blood pressure hemeostasis (Oudit et al., 2003). The coronavirus entry into host cells is mediated by the spike glycoprotein, which is a surface transmembrane protein in SARS-CoV-2 (Zhao et al., 2020). The analysis of the receptor-binding motif (RDM) in the Spike glycoprotein revealed that most of the aminoacid residues essential for receptor-binding with ACE-2 were conserved between SARS-CoV-2 and SARS-CoV, demonstrating that these viruses use the same host receptor for cell entry (Yan et al., 2020). Hoffmann et al. (2020) proved that anti-human ACE-2 antibody (R&D Systems, Catalog #AF933) can inhibit the Spike protein-associated entry into cultured cells in vitro. Accordingly, human ACE2 is considered as a host target for the treatment of COVID-19 to avoid SARS-CoV-2 from entering host cells (Zhang L. et al., 2020).
The existing quinoline-based antimalarial drugs, hydroxychloroquine and chloroquine, have shown their potential in the treatment of COVID-19 (Kaur et al., 2010), which inspired us to identify the quinoline-based potent inhibitors against the therapeutic targets of SARS-CoV-2 using an in silico approach. Due to the drug-like properties and therapeutic potential, quinoline-derived compounds have sustained attention for developing novel drugs in future medicine (O'donnell et al., 2010). Quinolines are nitrogen-containing heterocyclic aromatic compounds, known to be versatile compounds because of their extensive uses in medicine, organic chemistry, and industrial chemistry (Prajapati et al., 2014). They are frequently found in several medicinal plants and are known to have antimalarial, anticancer, antibacterial, anti-fungal, anticonvulsant, anti-inflammatory, anthelminitc, cardiotonic, and analgesic activity (Hussaini, 2016). Some of the compounds with quinoline core are the preferred choice for the treatment of diverse ailments, especially cancer and malaria (Touret and de Lamballerie, 2020).
Materials and Methods
Ligand Preparation
Numerous medicinal plants and their phytocompounds have demonstrated their antiviral properties against a large group of viruses. Consequently, the phytocompounds of Diplocyclos palmatus leaf extract were subjected to docking analysis in the current study. In previous studies, the tropical medicinal plant of D. palmatus has been reported for its anti-biofilm, anti-infection, and anti-photoaging activity using Caenorhabditis elegans model (Alexpandi et al., 2019). The list of quinoline-drugs (total 113) was retrieved from DrugBank database (https://www.drugbank.ca/categories/DBCAT000788). The canonical SMILES of the compounds was retrieved from the PubChem database. The canonical SMILES of quinoline,1,2,3,4-tetrahydro-1-[(2-phenylcyclopropyl)sulfonyl]-trans-(8CI) was retrieved from the Guidechem database (https://www.guidechem.com/reference/dic-395649.html). Then, the PDB-format 3D-structure of compounds was downloaded from the Openbabel online server http://www.cheminfo.org/Chemistry/Cheminformatics/FormatConverter/index.html.
Protein Preparation
The 3D crystal protein-structures of SARS-CoV-2 3CLpro (PDB ID: 6LU7) (Hall and Ji, 2020), SARS-CoV-2 spike protein-ACE-2 receptor-binding domain (RBD) (PDB ID: 6M17) (Wu et al., 2020), and human ACE2 (PDB ID: 1R4L) (Joshi et al., 2020) were obtained from the RCSB PDB database (http://www.rcsb.org/pdb). The 3D crystal structures of SARS-CoV-2 RdRp generated through homology modeling using ICM 3.7.3 modeling software was gifted by Prof. Hua Li, Hubei Key Laboratory of Natural Medicinal Chemistry and Resource Evaluation, School of Pharmacy, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China (Wu et al., 2020). The energy minimization of targeted protein structures was performed using the YASARA server. The protein preparation was done with AutoDock Tools Version 1.5.6.
Molecular Docking
The virtual screening of best scoring compounds was performed using the iGEMDOCK with blind-mode docking. The iGEMDOCK tool is a graphical-automatic drug design system mainly used for structure-based virtual screening of drug molecules (Hsu et al., 2011). After the selection of best binding compounds, the interaction on the active domains of therapeutic targets (3CLpro, RdRp, and Spike-ACE2 complex) of selected compounds were analyzed using the AutoDock Vina tool. It is an open-source docking software, which extensively improves the average accuracy of the binding mode predictions of compounds better than other docking tools (Trott and Olson, 2010). It implements a competent optimization algorithm for estimating the affinity of protein-ligand interactions and predicting the plausible binding modes of compounds (Goodsell et al., 1996). Then, the ligand-protein interactions were visualized by Maestro 10 (Schrödinger) (Balasubramaniam et al., 2019). In the present study, to compare the selected quinolines with already reported anti-SARS-CoV-2 drugs, lopinavir (Hall and Ji, 2020) and remdesivir triphosphate (Gordon et al., 2020) were selected as positive inhibitors of SARS-CoV-2 for in silico analysis in the present study. Further, the selected quinolines and remdesivir triphosphate were compared with the parental nucleotides (NTPs) of SARS-CoV-2 RdRp for understanding the inhibition mode of viral replication.
In silico Drug-Likeliness and Bioactivity Prediction
The drug likeliness and bioactivity of quinolines were analyzed using the Molinspiration server (http://www.molinspiration.com). Molinspiration tool is a cheminformatics software that provides molecular properties as well as bioactivity prediction of compounds (Mabkhot et al., 2016). In the Molinspiration-based drug-likeness analysis, there are two important factors, including the lipophilicity level (log P) and polar surface area (PSA) directly associated with the pharmacokinetic properties (PK) of the compounds (Beetge et al., 2000). In the Molinspiration-based bioactivity analysis, the calculation of the bioactivity score of compounds toward GPCR ligands, ion channel modulators, kinase inhibitors, nuclear receptor ligands, protease inhibitors, and other enzyme targets were analyzed by sophisticated Bayesian statistics (Mabkhot et al., 2016). This was done as the protein families, such as G protein-coupled receptors (GPCR), ion channels, kinases, nuclear hormone receptors, proteases, and other enzymes (RdRp), are the major drug targets of most of the drugs (Hauser et al., 2017).
In silico ADMET Analysis
The PK properties, such as Absorption, Distribution, Metabolism, Excretion, and Toxicity (ADMET), of quinolines were predicted using the admerSAR v2.0 server (http://lmmd.ecust.edu.cn/admetsar2/). The admerSAR server is an open-source computational tool for prediction of ADMET properties of compounds, which makes it a practical platform for drug discovery and other pharmacological research (Guan et al., 2019). In the ADMET analysis, the absorption (A) of good drugs depends on factors such as membrane permeability [designated by colon cancer cell line (Caco-2)], human intestinal absorption (HIA), and the status of either P-glycoprotein substrate or inhibitor. The distribution (D) of drugs mainly depends on the ability to cross the blood-brain barrier (BBB). The metabolism (M) of drugs is calculated by the CYP, MATE1, and OATP1B1-OATP1B3 models. Excretion (E) of the drugs is estimated based on the renal OCT substrate. Then, the toxicity (T) of the drugs is predicted on the Human Ether-a-go-go- related gene inhibition, carcinogenic status, mutagenic status, and acute oral toxicity (Shen et al., 2012).
Results
Screening of Potent SARS-CoV-2 3CLpro Inhibitors
In Table S1, the iGEMDOCK-based virtual screening result reveals that among the 17 phytocompounds, the novel quinoline (NQ) identified from D. palmatus, quinoline,1,2,3,4-tetrahydro-1-[(2-phenylcyclopropyl)sulfonyl]-trans-(8CI), has a lower binding energy (-6.8 Kcal/mol) with SARS-CoV-2 3CLpro, compared to other phytocompounds. As shown in Figures 1D,I, the NQ builds five hydrophobic interactions (Cys145, Met49, Met165, Phe140, and Leu141), 8 polar interactions (His41, His164, Gln192, Gln189, His63, His172, Ser144, and Asn142), two negative-charged interactions (Asp187 and Glu166), and one unspecified residue interaction (Arg188). As shown in Figure 1G, the NQ forms a hydrogen bonding interaction with Glu166 (2.044 Å distance). Importantly, the NQ has effectively interacted with the active site of catalytic-dyad (Cys145 and His41) of the SARS-CoV-2 3CLpro (Figures 1B,F), and hence we assumed that the NQ is able to hinder the protease activity of SARS-CoV-2 like lopinavir (−6.6 Kcal/mol) (Figure S1), which is a protease-inhibitor based anti-SARS-CoV-2 drug reported for COVID-19 (Choy et al., 2020).
Figure 1. (A) The 3CLpro monomer has three domains, namely domain I (residues 8–101), domain II (residues 102–184), and domain III (residues 201–303). The red arrow mark indicates the region of the active site of 3CLpro, which is located in the gap between domains I and II, and has a CysHis catalytic dyad (Cys145 and His41). (B,C) Depicts the binding region of the quinoline,1,2,3,4-tetrahydro-1-[(2-phenylcyclopropyl)sulfonyl]-trans-(8CI) (−6.8 Kcal/mol) and saquinavir (−8.5 Kcal/mol) with SARS-CoV-2 3CLpro, respectively. (D,E) Shows the super-position view of the binding sites and their interacting amino acids of quinoline,1,2,3,4-tetrahydro-1-[(2-phenylcyclopropyl)sulfonyl]-trans-(8CI) and saquinavir with SARS-CoV-2 3CLpro, respectively. (F,H) Indicates the interaction with the active catalytic dyad of SARS-CoV-2 3CLpro (Cys145 and His41). (G) Hydrogen bond formation of quinoline,1,2,3,4-tetrahydro-1-[(2-phenylcyclopropyl)sulfonyl]-trans-(8CI) with Glu166 residues of 3CLpro at 2.044 Å distance. (I,J) Reveal the interacted aminoacid residues of SARS-CoV-2 3CLpro with quinoline,1,2,3,4-tetrahydro-1-[(2-phenylcyclopropyl)sulfonyl]-trans-(8CI) and saquinavir.
Based on the virtual screening results of SARS-CoV-2 3CLpro, rilapladib, saquinavir, oxolinic acid, elvitegravir, batefenterol, sitafloxacin, CP-609754, GSK-256066, alatrofloxacin, and quarfloxin were predicted to be the best compounds from the commercial quinoline-based drugs (Figure 2 and Table S2). Out of these compounds, only saquinavir (−8.5 Kcal/mol) was found to interact with the active catalytic-domain (Cys145 and His41) of the SARS-CoV-2 2CLpro (Figures 1H,C). As shown in Figures 1E,J, saquinavir forms seven hydrophobic interactions (Cys145, Met165, Met49, Leu27, Leu167, Leu141, and Phe140), 12 polar interactions (His41, Asn142, His163, Ser144, Ser46, Thr25, Thr26, Thr24, Gln192, Gln189, Thr190, and His164), one negative-charged interaction (Glu166), and one glycine interaction (Gly143) with SARS-CoV-2 3CLpro. Therefore, for the development of strong interactions at the CysHis catalytic dyad, saquinavir was predicted to obstruct the 3CLpro activity.
Figure 2. The binding patterns and aminoacid interactions of the best scoring quinolines with SARS-CoV-2 3CLpro. The red circle-marking indicates the interaction with the active site of SARS-CoV-2 3CLpro (CysHis catalytic dyad).
Screening of Potent SARS-CoV-2 RdRp Inhibitors
As shown in Figure 3 and Table S3, amongst the tested-quinolines, elvitegravir, oxolinic acid, saquinavir, garenoxacin, rilapladib, pelitinib, difloxacin, batefenterol, danofloxacin, and LGD-2226 were shown to be the best-docked compounds against SARS-CoV-2 RdRp. Figures 4A,C reveals that elvitegravir and oxolinic acid have the same binding energy (−7.1 Kcal/mol) and same binding site at SARS-CoV-2 RdRp, which is almost similar to natural nucleotides such as ATP (−7.6 Kcal/mol), UTP (−7.1 Kcal/mol), GTP (−7.7 Kcal/mol), and CTP (−7.1 Kcal/mol) (Figure S2). As shown in Figures 4B,I, elvitegravir forms five hydrophobic interactions (Tyr455, Ala554, Tyr619, Pro620, and Cys622), seven polar interactions (Thr556, Thr680, Ser681, Ser682, Thr687, Asn691, and Ser759), four negative-charged interactions (Asp452, Asp618, Asp623, and Asp760), and five unspecified residue interactions (Lys545, Arg553, Arg555, Lys621, and Lys798) with SARS-CoV-2 RdRp. Notably, elvitegravir builds a hydrogen bond interaction with Lys621 residue (1.931 Å distance), as shown in Figure 4F. Similarly, it formed five hydrophobic interactions (Tyr456, Met542, Ala554, Val557, and Ala558), six polar interactions (Thr556, Thr680, Ser681, Ser682, Thr687, and Asn691), two negative-charged interactions (Asp452 and Asp623), and five unspecified residue interactions (Lys545, Arg553, Arg555, Arg624, and Lys676) (Figures 4D,J). Subsequently, oxolinic acid forms three hydrogen bonding interactions with Thy456 (1.224Å), Ser682 (1.815Å), and Arg624 (2.904) residues of the SARS-CoV-2 RdRp (Figure 4H). As shown in Figures 4E,G, both elvitegravir and oxolinic acid have the ability to bind with the NTP binding channel (a set of hydrophilic residues such as Lys545, Arg553, and Arg555) of the SARS-CoV-2 RdRp, similar to remdesivir triphosphate (−7.8 Kcal/mol) (Figure S3), which is a RdRp-inhibitor based anti-SARS-CoV-2 agent reported for COVID-19 (Gao et al., 2020; Gordon et al., 2020). As shown in Figure S4, due to the same binding sites, we expect these two quinolines (elvitegravir and oxolinic acid) can readily interact with the NTP binding channel of SARS-CoV-2 RdRp, more quickly than the parental-nucleotides (NTPs) such as ATP, UTP, GTP, and CTP. Hence, these two quinoline drugs are anticipated to arrest the viral replication of SARS-CoV-2, as is seen with remdesivir triphosphate.
Figure 3. The binding patterns and aminoacid interactions of the best scoring quinolines with SARS-CoV-2 RdRp. The red circle-marking indicates the interaction with the NTP entry channel of SARS-CoV-2.
Figure 4. (A,B) Shows the binding region and close-up view of the interacting amino acids of elvitegravir (−7.1 Kcal/mol) with SARS-CoV-2 RdRp, respectively. (C,D) Individually depicts the binding region and super-position view of the interacting amino acids of oxolinic acid (−7.1 Kcal/mol) with SARS-CoV-2 RdRp, respectively. (E,G) Visualizes the interaction with the NTP entry channel of SARS-CoV-2 RdRp (a set of hydrophilic residues such as Lys545, Arg553, and Arg555) of elvitegravir and oxolinic acid. (F,H) Shows the hydrogen bond formation of elvitegravir (Lys621 with a distance of 1.931 Å) and oxolinic acid (Thy456, Ser682, and Arg624 with a distance of 1.224, 1.815, and 2.904 Å, respectively). (I,J) Illustrate the interacted aminoacid residues of SARS-CoV-2 3CLpro with elvitegravir and oxolinic acid.
Screening of Potent Inhibitor for Spike Protein-RBD-ACE2 Interaction
In Figure S5 and Figure 6A, the crystal structure of spike protein-ACE2 (PDB: 6M17) revealed that the aminoacid residues of ACE2, including Gln24, Thr27, Asp30, Lys31, His34, Glu35, Glu37, Asp38, Tyr41, Gln42, Met82, Lys353, Gly354, Asp355, and Arg357, were recognized as the Spike protein receptor-binding domain (RBD) or entry receptor site to invade the target cells (Benítez-Cardoza and Vique-Sánchez, 2020). Based on the virtual screening results of ACE2 protein, CP-609754, saquinavir, rilapladib, quarfloxin, batefenterol, oxolinic acid, alatrofloxacin, dovitinib, GSK-256066, and rebamipide were predicted as the best compounds, exhibiting a high binding affinity to ACE2 with low energy (Figure 5 and Table S4). However, these quinolines were not predicted to interact with the RBD of the Spike-ACE2 complex (Figure 5). The only compound that could target the RBD interface between Spike and ACE2 was rilapladib, as shown in Figure 6B. Rilapladib was predicted to be positioned on the central shallow pit of the RBD of the Spike-ACE2 complex with strong interactions (Figure 6D). The aminoacid residues of the Spike-RBD-ACE2 complex that interact with rilapladib were His34, Glu35, Glu37, Asp38, and Leu39, as shown in Figure 6C. As a result of superimposing the Spike-RBD-ACE2 complex to the rilapladib-RBD-ACE2 complex, an individual overlap of rilapladib with the interface of ACE2 was observed, signifying that rilapladib may possibly interrupt the interaction of the Spike-RBD-ACE2 complex.
Figure 5. The binding patterns and aminoacid interactions of the top ten scoring quinolines with the RBD interface of the Spike-ACE2 complex. The red circle-marking indicates the interaction with the RBD interface of Spike-ACE2 complex.
Figure 6. (A) Shows the receptor binding domain (RBD) of SARS-CoV-2 Spike protein with human ACE-2 receptor for viral entry into host cells. (B) The interrupting mode of rilapladib (−9.7 Kcal/mol) on the interface of RBD of the Spike protein-ACE2 complex. (C) Super-position view of the binding interations of rilapladib on the interface of RBD of the Spike protein-ACE2 complex. (D) Shows the interacted aminoacid residues of rilapladib on the interface of RBD of the Spike protein-ACE2 complex.
In silico Prediction of Drug-Likeness Property and Bioactivity Score
The adequacy of therapeutic drugs mainly depends on the molecular property and bioactivity of the compounds (Shen et al., 2012). To predict the drug-likeness and bioactivity of the selected quinolines, the in silico molecular property assessment was performed using the Molinspiration tool. This tool measures the milogP value (Octanol-water partition coefficient logP) and TPSA (Topological polar surface area) values of the compounds using Bayesian statistics. The result shows that the milogP value of quinoline,1,2,3,4-tetrahydro-1-[(2-phenylcyclopropyl)sulfonyl]-trans-(8CI) (3.63), saquinavir (4.26), elvitegravir (3.58), and oxolinic acid (0.68) were predicted as having ideal lipophilicity (logP < 5) (Han et al., 2019); rilapladib (7.33) was predicted as having poor lipophilicity (logP > 5) in the aspect of absorption and permeation (Figure S6). The TPSA of the quinoline,1,2,3,4-tetrahydro-1-[(2-phenylcyclopropyl)sulfonyl]-trans-(8CI) (37.38), elvitegravir (88.77), oxolinic acid (77.77), and rilapladib (54.79) were < 100, showing that these compounds had superior oral-absorption or membrane permeability than saquinavir (166.75), lopinavir (119.99), and remdesivir triphosphate (289.53) (Bakht et al., 2010). On the other hand, the majority of drug targets of existing drugs are in one of the following protein families: G protein-coupled receptors (GPCR), ion channels, kinases, nuclear hormone receptors, proteases, and other enzymes. In Figure S7, the in silico bioactivity prediction analysis also divulges that the NQ and saquinavir were predicted as protease inhibitors. Subsequently, elvitegravir and oxolinic acid were predicted as enzyme inhibitors, which means that these are able to inhibit other enzymes, including RdRp enyzme, except G protein-coupled receptors (GPCR), ion channels, kinases, nuclear hormone receptors, and proteases. This is because the protein families are the major drug targets of most of the drugs (Hauser et al., 2017). These data support the outcome of the predicted in-silico activity of selected quinolines against SARS-CoV-2.
The half-life (t1/2) of drugs is a valuable pharmacokinetic factor as it provides an exact indication of the duration of time that the effect of the drug continues in an individual. The period of action of a drug is called its half-life. The distribution half-life (t1/2a) is the time required to divide the plasma concentration by two after reaching pseudo-equilibrium, and not the time needed to eliminate half the administered dose. The elimination half-life (t1/2b) of drugs is defined as the time required for the concentration of the drug to reach half of its original value in the plasma or the total amount in the body. As mentioned in Table 1, the elimination half-life (t1/2b) of saquinavir (670.84 g/mol) is 7–12 h and their elimination half-life is 9–15 h (Taylor et al., 2001). The distribution half-life (t1/2a) of elvitegravir is calculated at 8.7 h, and the elimination half-life (t1/2b) as 12.9 h (Cada et al., 2013). Oxolinic acid (261.23 g/mol) has 1.3 and 84 h as its distribution half-life and elimination half-life, respectively (Samuelsen et al., 2003). Unfortunately, the half-life of rilapladib (735.81 g/mol) as well as quinoline,1,2,3,4-tetrahydro-1-[(2-phenylcyclopropyl)sulfonyl]-, trans-(8CI) (313.41 g/mol) are not available in the database. On the other hand, the elimination half-life of SARS-CoV-2 protease inhibitor lopinavir (628.8 g/mol) measured at 6.9 ± 2.2 h. Similarly, the reported anti-COVID-19 drug remdesivir has a very short half-life (0.39 h) and hence the human esterases hastily converted the remdesivir into nucleoside triphosphate metabolite (remdesivir triphosphate). However, the produced remdesivir triphosphate has a longer half-life of ~20 h in humans, in which the metabolite acts as an NTP analog and slows down the viral replication of SARS-CoV-2 (Eastman et al., 2020). Due to the long half-life of the NTP analog of remdesivir triphosphate, only one dose is required for daily administration for treating COVID-19. The data obtained revealed that the selected quinolines are significantly long in terms of drug half-life, in which the quinolines treatment will be successful in the drug repurposing against COVID-19.
Table 1. Half-life, molecular weight, and pharmacological functions of the suggested quinoline drugs for drug-repurposing against COVID-19.
In silico ADMET Analysis
In Table 2, the absorption (A) analysis reveals that the NQ and oxolinic acid were predicted to have high Caco-2 permeability. A HIA of <30% is classified to be poorly absorbed. The result reveals that all the quinolines were predicted to be highly absorbed by the human intestine. P-glycoprotein, a member of ATP-binding trans-membrane glycoprotein, excretes incoming drugs or other chemicals from the cells (Äänismaa and Seelig, 2007). The results revealed that the NQ and oxolinic acid are non-substrates or non-inhibitors for P-glycoprotein and other compounds are substrates/inhibitors for P-glycoprotein. In the distribution (D) analysis, the BBB permeability, logBBB >0.3, is thought to cross the BBB easily (Han et al., 2019). Here, all the quinolines were predicted to be BBB+ and they can easily cross the BBB. In the metabolism (M) analysis, the two main sub-types of cytochrome P450s (CYP) are CYP2D6 and CYP3A4, which are essential enzyme-systems for drug metabolism in the liver (Sams et al., 2004). The results showed that all the quinolines were non-inhibitors to CYP2D6; only NQ and saquinavir were substrates for CYP2D6 and CYP3A4. This data suggested that these quinolines may possibly metabolize in the liver. OATP1B1 and OATP1B3 are transporters expressed on the sinusoidal-membrane of hepatocytes, which interact with therapeutic-drugs as their substrates or inhibitors, and cause clinically relevant drug-drug interactions (Shitara, 2010). The result found that all the selected quinolines were predicted as inhibitors for OATP1B1-OATP1B3. The excretion (E) of drugs is associated with their hydrophilicity and molecular weight. In the kidney, organic cation transporters (OCTs) and multidrug and toxin extrusion proteins (MATEs) are the foremost transporters for the clearance of cationic drugs into the urine (Motohashi and Inui, 2013). The results showed that the quinolines were predicted as non-inhibitors to MATE-1, OCT-1, and OCT-2, which indicates the safety elimination profile of the quinolines. In the toxicity (T) analysis, the result showed that all quinolines were non-carcinogen and non-eye corrosive. In conclusion, the predicted results indicate that the ADMET properties of the quinolines are almost similar to the reported anti-SARS-CoV-2 drugs (lopinavir and remdesivir triphosphate), which promotes repurposing of these quinoline drugs for the treatment of COVID-19.
Discussion
As significant functional biological-macromolecules of coronavirus, the viral protease (3CLpro) and RNA-dependent RNA-polymerase (RdRp) are indispensable for proteolytic processing of the polyproteins as well as viral replication and have been considered as promising drug targets in the treatment of viral diseases (Zumla et al., 2016). Several drugs, including hydroxychloroquine, chloroquine, arbidol, remdesivir, favipiravir, lopinavir/ritonavir, interferon-α, and ribavirin are undergoing clinical trials to assess their anti-viral efficacy and safety level in the treatment of COVID-19 (Dong et al., 2020). Most of the reported anti-SARS-CoV-2 drugs are protease-inhibitors or RdRp-inhibitors (Elfiky, 2020; Wu et al., 2020).
One of the most-characterized therapeutic targets among coronaviruses is inhibiting the 3CLpro activity since this enzyme is crucial for processing the polyproteins that are translated from the RNA molecules (Ghosh et al., 2007; Khan et al., 2020). The 3CLpro, also called Nsp5 (non-structural protein 5), is first routinely cleaved from polyproteins to produce mature enzymes, and subsequently further cleaves downstream Nsps at 11 cleavage sites to release Nsp4-Nsp16 (Wu et al., 2020). 3CLpro directly mediates the maturation of Nsps, which is fundamental in the life-cycle of SARS-CoV-2 (Zhang H. et al., 2020). The 3CLpro monomer has three domains, namely domain I (residues 8–101), domain II (residues 102–184), and domain III (residues 201–303), and a long loop (residues 185–200) links domains II and III, as shown in Figure 1A. The active site of 3CLpro is located in the gap between domains I and II, and has a CysHis catalytic dyad (Cys145 and His41) (Jin et al., 2019; Wu et al., 2020). The active site of SARS-CoV-2 3CLpro is located in the gap between domains I and II, and has a CysHis catalytic dyad (Cys145 and His41) (Muralidharan et al., 2020). The cleavage by 3CLpro arises at the glutamine residue in the P1 position of the substrate through the CysHis catalytic dyad, wherein cysteine thiol functions as the nucleophile in the proteolytic process (Chen et al., 2005). Hence, inhibiting the activity of this enzyme would arrest the viral replication of SARS-CoV-2.
We saw that quinoline,1,2,3,4-tetrahydro-1-[(2-phenylcyclopropyl)sulfonyl]-trans-(8CI) (NQ) and saquinavir can target main proteases through authoritative interaction to the catalytic dyad (Cys145 and His41) of SARS-CoV-2 3CLpro (Figures 1F,H), and along these lines are believed to hinder the protease activity, as has been reported with anti-SARS-CoV-2 agent (lopinavir). As referenced before, we identified the NQ from the methanolic leaf extract of D. palmatus using GC-MS analysis (Alexpandi et al., 2019). Unfortunately, the NQ is not available commercially. Henceforth, the present study essentially evaluated their drug-likeness, bioactivity, and ADMET properties through an in silico approach. The outcomes demonstrated that the NQ has low milogP (3.63) and low TPSA value (37.38), which authenticates the ideal lipophilicity (milogP < 5) nature and higher oral-absorption or membrane permeability (TPSA < 100) than other quinolines. These are the physicochemical properties that play a fundamental role in deciding the ADMET properties of compounds (Shen et al., 2012; Han et al., 2019). In ADMET analysis, the results reveal that the NQ was predicted as Caco-2+, HIA+, a non-inhibitor to P-glycoprotein, BBB+, a non-inhibitor to CYP2D6, a non-substrate to CYP3A4, and a non-inhibitor to excretion-related receptors such as OCT-1, OCT-2, and MATE-1. The NQ showing some toxic impact on the human either-a-go-go related genes, has even been predicted as a non-carcinogenic, non-mutagenic, non-eye irritant, and non-eye corrosive. These results unequivocally suggest that the novel phytocompound, quinoline,1,2,3,4-tetrahydro-1-[(2-phenylcyclopropyl)sulfonyl]-trans-(8CI), can be used as a protease-inhibitor drug for the treatment of COVID-19, but the compound needs to be synthesized.
Saquinavir, the first FDA-approved HIV-1 drug, has the ability to cleave between Tyr-Pro or Phe-Pro of the HIV polyproteins, which is rare in mammalian systems (Noble and Faulds, 1996). Hence, saquinavir does not interfere with mammalian proteases, symbolizing its safety level for humans (Ganguly et al., 2011). Saquinavir is a protease inhibitor that binds to the active site of the viral protease and thereby blocks cleavage of viral polyproteins and maturation of the HIV-1 and HIV-2 (Kim et al., 1998). Furthermore, (Tan et al., 2004) showed the in vitro antiviral activity of saquinavir against SARS-CoV-1. The present study revealed that saquinavir is able to bind with the catalytic dyad, and is thereby anticipated to interrupt 3CLpro activity (Figure 1H). So, we recommend saquinavir as the potent protease-inhibitor for drug-repurposing against COVID-19.
In the research of anti-SARS-CoV-2 drug designing, RdRp has been well-thought-out as an incredibly potent drug target due to its central role in RNA-synthesis from RNA-templates (Gao et al., 2020). In addition, RdRp-inhibitors do not show considerable toxicity or side effects on host cells (Dong et al., 2020). The active site of the SARS-CoV-2 RdRp domain is formed by the conserved polymerase motifs A-G, within 549th to 776th aminoacid residues, which are essential for the RNA-directed 5'-3' polymerase activity (Shannon et al., 2020). As in other RNA-polymerases, the template/primer entry (known as nucleoside triphosphate (NTP) entry) and budding strand are congregating in a central cavity where the RdRp motifs intercede RNA-template mediated RNA synthesis in SARS-CoV-2. The NTP entry channel of SARS-CoV-2 is placed in the set of hydrophilic residues, including Lys545, Arg553, and Arg555 in motif-F (Gao et al., 2020). Therefore, the nucleotide-analog antiviral inhibitors such as remdesivir and favipiravir also showed their antiviral potential against SARS-CoV-2 (Gordon et al., 2020). Our results revealed that both elvitegravir and oxolinic acid have the ability to bind with the NTP binding channel of SARS-CoV-2 with a low binding energy (−7.1 Kcal/mol) (Figures 4E,G), similar to the parental nucleotides such as ATP (−7.6 Kcal/mol), UTP (−7.1 Kcal/mol), GTP (−7.7 Kcal/mol), and CTP (−7.1 Kcal/mol). We believed that these two quinolines can more readily interact with the NTP binding channel than parental-nucleotides (especially than UTP and CTP), and thereby, possibly block the de novo addition of NTP to the 3'-OH strand, which leads to the arrest of viral replication (Gao et al., 2020), as illustrated in Figure S4. Elvitegravir is an integrase inhibitor used for the anti-retroviral treatment of HIV-1 (Shimura et al., 2008; Ramanathan et al., 2011). Oxolinic acid is a synthetic quinoline-derived antibiotic used to treat bacteria causing urinary tract infections (Sato et al., 2006; Irgi et al., 2015). Owing to the nucleotide-antagonistic behavior, we suggest that elvitegravir and oxolinic acid might be the potent RdRp inhibitors for SARS-CoV-2 in the treatment of COVID-19.
On the other hand, ACE2 is a type I transmembrane metallocarboxypeptidase, an enzyme that plays a crucial role in the rennin-angiotensin (RAS) system and is considered as a target for the treatment of hypertension (Burrell et al., 2004). ACE2 is widely distributed in the human body and has been associated with the protective function in the cardiovascular system and other organs (Yagil et al., 2003). In contrast, human ACE2 is the recognized functional receptor for the Spike glycoprotein of SARS-CoV-2 that initiates cell entry into host cells and viral replication in the target cells (Lan et al., 2020). The previous report proved that ACE2 knockout significantly reduces the viral load in mice after the experimental SARS-CoV infection (Hoffmann et al., 2020). As shown in Figure 3, though several quinolines could bind with ACE2, none was found to bind with the RBD of the ACE2-Spike complex. Moreover, these kinds of ACE2 inhibitors may not be appropriate for treating COVID-19 because these drugs can inhibit ACE2 enzyme activities, and cause lung injury and heart failure (Velkoska et al., 2016). To predict the RBD interface binding compound, we found one quinoline-drug, rilapladib, was targeting the RBD of the Spike-ACE2 complex, as shown in Figure 6B. Rilapladib, a hydroquinoline-based small molecule drug developed by GlaxoSmithKline was used as a lipoprotein-associated phospholipase A2 (Lp-PLA2) inhibitor or 1-alkyl-2-acetylglycerophosphocholine esterase inhibitor for treating atherosclerotic plaques and Alzheimer's disease (Shaddinger et al., 2014). It's worth mentioning that rilapladib was well-fitted into the interface of the RBD of the Spike-ACE2 complex (Figure 6B), in which lots of interactions with His34, Glu35, Glu37, Asp38, Leu39, Lys68, Ala71, Phe72, Glu75, Gln76, Leu79, and Lys353 create a strong binding with the interface of RBD and block the Spike-ACE2 interactions. Owing to the formation of possible interactions at the RBD interface of the Spike-ACE2 complex, the present study suggests rilapladib prevents ACE2-mediated viral entry of SARS-CoV-2.
Moreover, the in silico ADMET results demonstrated that these quinolines were non-toxic, non-carcinogenic, absorb in the human intestine, have Caco-2 permeability, do not inhibit CYP enzymes, are non-inhibitors for RCT, and non-inhibitors for Human Ether-a-go-go related genes, which suggested their significant pharmacokinetic properties. Further, the drug half-life of selected quinoline drugs are significantly long, in which these quinolines were expected to offer an efficient drug distribution against COVID-19. Overall, we believe that these quinolines may be efficient drug candidates for the development of efficient therapeutics against COVID-19 in this pandemic period.
Conclusion
In conclusion, the current investigation recommends that the novel phyto-quinoline (quinoline,1,2,3,4-tetrahydro-1-[(2-phenylcyclopropyl)sulfonyl]-trans-(8CI)) and the existing quinoline-based drugs saquinavir, elvitegravir, and oxolinic acid could be used as potent inhibitors for SARS-CoV-2. Subsequently, rilapladib is suggested for anti-Spike-RBD-ACE2 therapy to avoid ACE2-mediated viral entry of SARS-CoV-2 into the host cells. Notwithstanding, further in vitro and in vivo experiments are needed to transform these potential inhibitors into clinical drugs. We anticipate that this new finding could significantly impact the development of therapeutic agents for COVID-19 in the future.
Data Availability Statement
All datasets presented in this study are included in the article/Supplementary Material.
Author Contributions
RA: conceptualization, performed the experiments, data analysis, and writing-original draft. JD: data analysis and reviewing-original draft. SP: supervision and reviewing-original draft. AR: conceptualization, supervision, and reviewing-original draft. All authors contributed to the article and approved the submitted version.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
RA sincerely acknowledges Rashtriya Uchchatar Shiksha Abhiyan (RUSA), MHRD, Government of India [F. 24-51/2014-U, Policy (TN Multi-Gen), Dept. of Edn, GOI] for providing RUSA-2.0 Ph.D. Fellowship. We sincerely thank Prof. Hua Li, Huazhong University of Science and Technology, Wuhan, China for providing the modeled-structure of SARS-CoV-2 RdRp. RA thanks Mr. Ragupathi Thennarasu, DBT-BIF Research Assistant, ALU for assistance during data analysis. We also acknowledge the BIF (funded by DBT, File No. BT/BI/25/012/2012), DST-FIST [Grant No. SR/FST/LSI-639/2015(C)], UGC-SAP [Grant No. F.5-1/2018/DRS-II (SAP-II)], DST-PURSE [Grant No. SR/PURSE Phase 2/38 (G)], and RUSA- Phase 2.0, Government of India [F. 24-51/2014-U, Policy (TN Multi-Gen)] for providing instrumentation facilities.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2020.01796/full#supplementary-material
References
Äänismaa, P., and Seelig, A. (2007). P-Glycoprotein kinetics measured in plasma membrane vesicles and living cells. Biochemistry 46, 3394–3404. doi: 10.1021/bi0619526
Agarwal, S., Pal, D., and Mitra, A. K. (2007). Both P-gp and MRP2 mediate transport of Lopinavir, a protease inhibitor. Int. J. Pharm. 339, 139–147. doi: 10.1016/j.ijpharm.2007.02.036
Alexpandi, R., Prasanth, M. I., Ravi, A. V., Balamurugan, K., Durgadevi, R., Srinivasan, R., et al. (2019). Protective effect of neglected plant Diplocyclos palmatus on quorum sensing mediated infection of Serratia marcescens and UV-A induced photoaging in model Caenorhabditis elegans. J. Photoch. Photobio. B 201:111637. doi: 10.1016/j.jphotobiol.2019.111637
Bakht, M. A., Yar, M. S., Abdel-Hamid, S. G., Al Qasoumi, S. I., and Samad, A. (2010). Molecular properties prediction, synthesis and antimicrobial activity of some newer oxadiazole derivatives. Eur. J. Med. Chem. 45, 5862–5869. doi: 10.1016/j.ejmech.2010.07.069
Balasubramaniam, B., Alexpandi, R., and Darjily, D. R. (2019). Exploration of the optimized parameters for bioactive prodigiosin mass production and its biomedical applications in vitro as well as in silico. Biocatal. Agric. Biotechnol. 22:101385. doi: 10.1016/j.bcab.2019.101385
Barry, A. L., Jones, R. N., Thornsberry, C., Ayers, L. W., Gerlach, E. H., and Sommers, H. M. (1984). Antibacterial activities of ciprofloxacin, norfloxacin, oxolinic acid, cinoxacin, and nalidixic acid. Antimicrob. Agents. Ch. 25, 633–637.
Beetge, E., du Plessis, J., Müller, D. G., Goosen, C., and van Rensburg, F. J. (2000). The influence of the physicochemical characteristics and pharmacokinetic properties of selected NSAID's on their transdermal absorption. Int. J. Pharm. 193, 261–264. doi: 10.1016/S0378-5173(99)00340-3
Benítez-Cardoza, C. G., and Vique-Sánchez, J. L. (2020). Potential inhibitors of the interaction between ACE2 and SARS-CoV-2 (RBD), to develop a drug. Life. Sci. 256:117970. doi: 10.1016/j.lfs.2020.117970
Burrell, L. M., Johnston, C. I., Tikellis, C., and Cooper, M. E. (2004). ACE2, a new regulator of the renin–angiotensin system. Trends. Endocrin. Met. 15, 166–169. doi: 10.1016/j.tem.2004.03.001
Cada, D. J., Torres, S., Levien, T. L., and Baker, D. E. (2013). Elvitegravir/Cobicistat/Emtricitabine/Tenofovir disoproxil fumarate tablets. Hosp. Pharm. 48, 48–56. doi: 10.1310/hpj4801-48
Chen, L. R., Wang, Y. C., Lin, Y. W., Chou, S. Y., Chen, S. F., Liu, L. T., et al. (2005). Synthesis and evaluation of isatin derivatives as effective SARS coronavirus 3CL protease inhibitors. Bioorg. Med. Chem. Lett. 15, 3058–3062. doi: 10.1016/j.bmcl.2005.04.027
Choy, K. T., Wong, A. Y. L., Kaewpreedee, P., Sia, S. F., Chen, D., Hui, K. P. Y., et al. (2020). Remdesivir, lopinavir, emetine, and homoharringtonine inhibit SARS-CoV-2 replication in vitro. Antivir. Res. 178:104786. doi: 10.1016/j.antiviral.2020.104786
Dey, S. K., Rahman, M. M., Siddiqi, U. R., and Howlader, A. (2020). Analyzing the epidemiological outbreak of COVID-19: a visual exploratory data analysis approach. J. Med. Virol. 92, 632–638. doi: 10.1002/jmv.25743
Dong, L., Hu, S., and Gao, J. (2020). Discovering drugs to treat coronavirus disease 2019 (COVID-19). Drug. Discov. Ther. 14, 58–60. doi: 10.5582/ddt.2020.01012
Eastman, R. T., Roth, J. S., Brimacombe, K. R., Simeonov, A., Shen, M., Patnaik, S., et al. (2020). Remdesivir: a review of its discovery and development leading to emergency use authorization for treatment of COVID-19. ACS Cent. Sci. 6, 672–683. doi: 10.1021/acscentsci.0c00489
Elfiky, A. A. (2020). Ribavirin, Remdesivir, Sofosbuvir, Galidesivir, and Tenofovir against SARS-CoV-2 RNA dependent RNA polymerase (RdRp): a molecular docking study. Life. Sci. 253:117592. doi: 10.1016/j.lfs.2020.117592
Ganguly, A. K., Alluri, S. S., Caroccia, D., Biswas, D., Wang, C. H., Kang, E., et al. (2011). Design, synthesis, and X-ray crystallographic analysis of a novel class of HIV-1 protease inhibitors. J. Med. Chem. 54, 7176–7183. doi: 10.1021/jm200778q
Gao, Y., Yan, L., Huang, Y., Liu, F., Zhao, Y., Cao, L., et al. (2020). Structure of the RNA-dependent RNA polymerase from COVID-19 virus. Science 368, 779–782. doi: 10.1126/science.abb7498
Ghosh, A. K., Xi, K., Grum-Tokars, V., Xu, X., Ratia, K., Fu, W., et al. (2007). Structure-based design, synthesis, and biological evaluation of peptidomimetic SARS-CoV 3CLpro inhibitors. Bioorg. Med. Chem. Lett. 17, 5876–5880. doi: 10.1016/j.bmcl.2007.08.031
Goodsell, D. S., Morris, G. M., and Olson, A. J. (1996). Automated docking of flexible ligands: applications of AutoDock. J. Mol. Recognit. 9, 1–5.
Gordon, C. J., Tchesnokov, E. P., Woolner, E., Perry, J. K., Feng, J. Y., Porter, D. P., et al. (2020). Remdesivir is a direct-acting antiviral that inhibits RNA-dependent RNA polymerase from severe acute respiratory syndrome coronavirus 2 with high potency. J. Biol. Chem. 295:jbc.RA120.013679. doi: 10.1074/jbc.RA120.013679
Guan, L., Yang, H., Cai, Y., Sun, L., Di, P., Li, W., et al. (2019). ADMET-score–a comprehensive scoring function for evaluation of chemical drug-likeness. Medchemcomm 10, 148–157. doi: 10.1039/C8MD00472B
Hall, D. C. Jr, and Ji, H. F. (2020). A search for medications to treat COVID-19 via in silico molecular docking models of the SARS-CoV-2 spike glycoprotein and 3CL protease. Travel. Med. Infect. Dis. 35:101646. doi: 10.1016/j.tmaid.2020.101646
Han, Y., Zhang, X., Zhang, J., and Hu, C. Q. (2019). In silico ADME and toxicity prediction of ceftazidime and its impurities. Front. Pharmacol. 10:434. doi: 10.3389/fphar.2019.00434
Hauser, A. S., Attwood, M. M., Rask-Andersen, M., Schiöth, H. B., and Gloriam, D. E. (2017). Trends in GPCR drug discovery: new agents, targets and indications. Nat. Rev. Drug Discov. 16, 829–842. doi: 10.1038/nrd.2017.178
Hoffmann, M., Kleine-Weber, H., Schroeder, S., Krüger, N., Herrler, T., Erichsen, S., et al. (2020). SARS-CoV-2 cell entry depends on ACE2 and TMPRSS2 and is blocked by a clinically proven protease inhibitor. Cell 181, 271–280.e8. doi: 10.1016/j.cell.2020.02.052
Hsu, K. C., Chen, Y. F., Lin, S. R., and Yang, J. M. (2011). iGEMDOCK: a graphical environment of enhancing GEMDOCK using pharmacological interactions and post-screening analysis. BMC bioinformatics 12:S33. doi: 10.1186/1471-2105-12-S1-S33
Husna Ibrahim, N., Yahaya, M. F., Mohamed, W., Teoh, S. L., Hui, C. K., and Kumar, J. (2020). Pharmacotherapy of Alzheimer's Disease: seeking clarity in a time of uncertainty. Front. Pharmacol. 11:261. doi: 10.3389/fphar.2020.00261
Hussaini, S. M. A. (2016). Therapeutic significance of quinolines: a patent review (2013-2015). Expert. Opin. Ther. Pat. 26, 1201–1221. doi: 10.1080/13543776.2016.1216545
Irgi, E. P., Geromichalos, G. D., Balala, S., Kljun, J., Kalogiannis, S., Papadopoulos, A., et al. (2015). Cobalt (II) complexes with the quinolone antimicrobial drug oxolinic acid: structure and biological perspectives. RSC Adv. 5, 36353–36367. doi: 10.1039/C5RA05308K
Jin, Z., Du, X., Xu, Y., Deng, Y., Liu, M., Zhao, Y., et al. (2019). Structure of M pro from SARS-CoV-2 and discovery of its inhibitors. Nature 582, 289–293. doi: 10.1038/s41586-020-2223-y
Joshi, T., Joshi, T., Sharma, P., Mathpal, S., Pundir, H., Bhatt, V., et al. (2020). In silico screening of natural compounds against COVID-19 by targeting Mpro and ACE2 using molecular docking. Eur. Rev. Med. Pharmacol. Sci. 24, 4529–4536. doi: 10.26355/eurrev_202004_21036
Kaur, K., Jain, M., Reddy, R. P., and Jain, R. (2010). Quinolines and structurally related heterocycles as antimalarials. Eur. J. Med. Chem. 45, 3245–3264. doi: 10.1016/j.ejmech.2010.04.011
Khan, R. J., Jha, R. K., Amera, G., Jain, M., Singh, E., Pathak, A., et al. (2020). Targeting SARS-Cov-2: a systematic drug repurposing approach to identify promising inhibitors against 3C-like Proteinase and 2'-O-RiboseMethyltransferase. J. Biomol. Struct. Dyn. 1–40. doi: 10.1080/07391102.2020.1753577. [Epub ahead of print].
Kim, A. E., Dintaman, J. M., Waddell, D. S., and Silverman, J. A. (1998). Saquinavir, an HIV protease inhibitor, is transported by P-glycoprotein. J. Pharmacol. Exp. Ther. 286, 1439–1445.
Lampiris, H. W. (2012). Elvitegravir: a once-daily, boosted, HIV-1 integrase inhibitor. Expert. Rev. Anti-infe. 10, 13–20. doi: 10.1586/eri.11.157
Lan, J., Ge, J., Yu, J., Shan, S., Zhou, H., Fan, S., et al. (2020). Structure of the SARS-CoV-2 spike receptor-binding domain bound to the ACE2 receptor. Nature 581, 215–220. doi: 10.1038/s41586-020-2180-5
Mabkhot, Y. N., Alatibi, F., El-Sayed, N. N., Al-Showiman, S., Kheder, N. A., Wadood, A., et al. (2016). Antimicrobial activity of some novel armed thiophene derivatives and petra/osiris/molinspiration (POM) analyses. Molecules 21:222. doi: 10.3390/molecules21020222
Motohashi, H., and Inui, K. I. (2013). Organic cation transporter OCTs (SLC22) and MATEs (SLC47) in the human kidney. AAPS J. 15, 581–588. doi: 10.1208/s12248-013-9465-7
Muralidharan, N., Sakthivel, R., Velmurugan, D., and Gromiha, M. M. (2020). Computational studies of drug repurposing and synergism of lopinavir, oseltamivir and ritonavir binding with SARS-CoV-2 Protease against COVID-19. J. Biomol. Struct. Dyn. 1–6. doi: 10.1080/07391102.2020.1752802. [Epub ahead of print].
Noble, S., and Faulds, D. (1996). Saquinavir. Drugs. 52, 93–112. doi: 10.2165/00003495-199652010-00007
O'donnell, F., Smyth, T. J. P., Ramachandran, V. N., and Smyth, W. F. (2010). A study of the antimicrobial activity of selected synthetic and naturally occurring quinolines. Int. J. Antimicrob. Agents. 35, 30–38. doi: 10.1016/j.ijantimicag.2009.06.031
Oudit, G. Y., Crackower, M. A., Backx, P. H., and Penninger, J. M. (2003). The role of ACE2 in cardiovascular physiology. Trends. Cardiovas. Med. 13, 93–101. doi: 10.1016/S1050-1738(02)00233-5
Patel, V. B., Clarke, N., Wang, Z., Fan, D., Parajuli, N., Basu, R., et al. (2014). Angiotensin II induced proteolytic cleavage of myocardial ACE2 is mediated by TACE/ADAM-17: a positive feedback mechanism in the RAS. J. Mol. Cell. Cardiol. 66, 167–176. doi: 10.1016/j.yjmcc.2013.11.017
Prajapati, S. M., Patel, K. D., Vekariya, R. H., Panchal, S. N., and Patel, H. D. (2014). Recent advances in the synthesis of quinolines: a review. Rsc. Adv. 4, 24463–24476. doi: 10.1039/C4RA01814A
Ramanathan, S., Mathias, A. A., German, P., and Kearney, B. P. (2011). Clinical pharmacokinetic and pharmacodynamic profile of the HIV integrase inhibitor elvitegravir. Clin. Pharmacokinet. 50, 229–244. doi: 10.1128/AAC.25.5.633
Sams, C., Cocker, J., and Lennard, M. S. (2004). Biotransformation of chlorpyrifos and diazinon by human liver microsomes and recombinant human cytochrome P450s (CYP). Xenobiotica 34, 861–873. doi: 10.1080/00498250400017273
Samuelsen, O. B., Bergh, O., and Ervik, A. (2003). A single-dose pharmacokinetic study of oxolinic acid and vetoquinol, an oxolinic acid ester, in cod, Gadus morhua L., held in sea water at 8 degrees C and in vitro antibacterial activity of oxolinic acid against Vibrio anguillarum strains isolated from diseased cod. J. Fish Dis. 26, 339–347. doi: 10.1046/j.1365-2761.2003.00466.x
Sato, M., Motomura, T., Aramaki, H., Matsuda, T., Yamashita, M., Ito, Y., et al. (2006). Novel HIV-1 integrase inhibitors derived from quinolone antibiotics. J. Med. Chem. 49, 1506–1508. doi: 10.1021/jm0600139
Shaddinger, B. C., Xu, Y., Roger, J. H., Macphee, C. H., Handel, M., Baidoo, C. A., et al. (2014). Platelet aggregation unchanged by lipoprotein-associated phospholipase A2 inhibition: results from an in vitro study and two randomized phase I trials. PloS ONE 9:e83094. doi: 10.1371/journal.pone.0083094
Shannon, A., Le, N. T. T., Selisko, B., Eydoux, C., Alvarez, K., Guillemot, J. C., et al. (2020). Remdesivir and SARS-CoV-2: structural requirements at both nsp12 RdRp and nsp14 Exonuclease active-sites. Antivir. Res. 178:104793. doi: 10.1016/j.antiviral.2020.104793
Shen, M., Tian, S., Li, Y., Li, Q., Xu, X., Wang, J., et al. (2012). Drug-likeness analysis of traditional Chinese medicines: 1. Property distributions of drug-like compounds, non-drug-like compounds and natural compounds from traditional Chinese medicines. J. Cheminformat. 4:31. doi: 10.1186/1758-2946-4-31
Shimura, K., Kodama, E., Sakagami, Y., Matsuzaki, Y., Watanabe, W., Yamataka, K., et al. (2008). Broad antiretroviral activity and resistance profile of the novel human immunodeficiency virus integrase inhibitor elvitegravir (JTK-303/GS-9137). J. Virol. 82, 764–774. doi: 10.1128/JVI.01534-07
Shitara, Y. (2010). Clinical importance of OATP1B1 and OATP1B3 in drug-drug interactions. Drug. Metab. Pharmacok. 26, 220–227. doi: 10.2133/dmpk.DMPK-10-RV-094
Stebbing, J., Phelan, A., Griffin, I., Tucker, C., Oechsle, O., Smith, D., et al. (2020). COVID-19: combining antiviral and anti-inflammatory treatments. Lancet. Infect. Dis. 20, 400–402. doi: 10.1016/S1473-3099(20)30132-8
Tan, E. L., Ooi, E. E., Lin, C. Y., Tan, H. C., Ling, A. E., Lim, B., et al. (2004). Inhibition of SARS coronavirus infection in vitro with clinically approved antiviral drugs. Emerg. Infect. Dis. 10:581. doi: 10.3201/eid1004.030458
Taylor, S., Back, D. J., Drake, S. M., Workman, J., Reynolds, H., Gibbons, S. E., et al. (2001). Antiretroviral drug concentrations in semen of HIV-infected men: differential penetration of indinavir, ritonavir and saquinavir. J. Antimicrob. Chemother. 48, 351–354. doi: 10.1093/jac/48.3.351
Thomas-Rüddel, D., Winning, J., Dickmann, P., Ouart, D., Kortgen, A., Janssens, U., et al. (2020). Coronavirus disease 2019 “(COVID-19): update für Anästhesisten und Intensivmediziner März 2020. Der. Anaesthesist. 69, 225–235. doi: 10.1007/s00101-020-00758-x
Touret, F., and de Lamballerie, X. (2020). Of chloroquine and COVID-19. Antivir. Res. 177:104762. doi: 10.1016/j.antiviral.2020.104762
Trott, O., and Olson, A. J. (2010). AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J. Comput. Chem. 31, 455–461. doi: 10.1002/jcc.21334
Velkoska, E., Patel, S. K., and Burrell, L. M. (2016). Angiotensin converting enzyme 2 and diminazene: role in cardiovascular and blood pressure regulation. Curr. Opin. Nephrol. Hypertens. 25, 384–395. doi: 10.1097/MNH.0000000000000254
Walmsley, S., Bernstein, B., King, M., Arribas, J., Beall, G., Ruane, P., et al. (2002). Lopinavir–ritonavir versus nelfinavir for the initial treatment of HIV infection. N. Engl. J. Med. 346, 2039–2046. doi: 10.1056/NEJMoa012354
Wang, M., Cao, R., Zhang, L., Yang, X., Liu, J., Xu, M., et al. (2020). Remdesivir and chloroquine effectively inhibit the recently emerged novel coronavirus (2019-nCoV) in vitro. Cell. Res. 30, 269–271. doi: 10.1038/s41422-020-0282-0
Wang, Y., Zhang, D., Du, G., Du, R., Zhao, J., Jin, Y., et al. (2020). Remdesivir in adults with severe COVID-19: a randomised, double-blind, placebo-controlled, multicentre trial. The Lancet 395, 1569–1578. doi: 10.1016/S0140-6736(20)31022-9
Wright, H. T., Nurse, K. C., and Goldstein, D. J. (1981). Nalidixic acid, oxolinic acid, and novobiocin inhibit yeast glycyl-and leucyl-transfer RNA synthetases. Science 213, 455–456. doi: 10.1126/science.7017932
Wu, C., Liu, Y., Yang, Y., Zhang, P., Zhong, W., Wang, Y., et al. (2020). Analysis of therapeutic targets for SARS-CoV-2 and discovery of potential drugs by computational methods. Acta. Pharma. Sin. B 10, 766–788. doi: 10.1016/j.apsb.2020.02.008
Yagil, Y., Yagil, C., Crackower, M., and Penninger, J. (2003). OR-51: the protective role of ACE2 in hypertension. Am. J. Hypertens. 16:23A. doi: 10.1016/S0895-7061(03)00135-3
Yan, R., Zhang, Y., Li, Y., Xia, L., Guo, Y., and Zhou, Q. (2020). Structural basis for the recognition of SARS-CoV-2 by full-length human ACE2. Science 367, 1444–1448. doi: 10.1126/science.abb2762
Yao, X., Ye, F., Zhang, M., Cui, C., Huang, B., Niu, P., et al. (2020). In vitro antiviral activity and projection of optimized dosing design of hydroxychloroquine for the treatment of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Clin. Infect. Dis. doi: 10.1093/cid/ciaa237
Zhang, H., Penninger, J. M., Li, Y., Zhong, N., and Slutsky, A. S. (2020). Angiotensin-converting enzyme 2 (ACE2) as a SARS-CoV-2 receptor: molecular mechanisms and potential therapeutic target. Intensive. Care. Med. 46, 586–590. doi: 10.1007/s00134-020-05985-9
Zhang, L., Lin, D., Sun, X., Curth, U., Drosten, C., Sauerhering, L., et al. (2020). Crystal structure of SARS-CoV-2 main protease provides a basis for design of improved α-ketoamide inhibitors. Science 368, 409–412. doi: 10.1126/science.abb3405
Zhao, Y., Zhao, Z., Wang, Y., Zhou, Y., Ma, Y., and Zuo, W. (2020). Single-cell RNA expression profiling of ACE2, the receptor of SARS-CoV-2. Biorxiv [Preprint]. doi: 10.1101/2020.01.26.919985
Zou, L., Ruan, F., Huang, M., Liang, L., Huang, H., Hong, Z., et al. (2020). SARS-CoV-2 viral load in upper respiratory specimens of infected patients. N. Engl. J. Med. 382, 1177–1179. doi: 10.1056/NEJMc2001737
Keywords: drug repurposing, SARS-CoV-2, main-protease (3CLpro), RNA-dependent RNA-polymerase, spike-ACE2 complex, quinoline based-drugs
Citation: Alexpandi R, De Mesquita JF, Pandian SK and Ravi AV (2020) Quinolines-Based SARS-CoV-2 3CLpro and RdRp Inhibitors and Spike-RBD-ACE2 Inhibitor for Drug-Repurposing Against COVID-19: An in silico Analysis. Front. Microbiol. 11:1796. doi: 10.3389/fmicb.2020.01796
Received: 12 May 2020; Accepted: 09 July 2020;
Published: 23 July 2020.
Edited by:
Nejat Duzgunes, University of the Pacific, United StatesReviewed by:
Fouzia Sadiq, Shifa Tameer-e-Millat University, PakistanSiew Pheng Lim, Denka Life Innovation Research (DLIR), Singapore
Copyright © 2020 Alexpandi, De Mesquita, Pandian and Ravi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Arumugam Veera Ravi, YXZlZXJhcmF2aSYjeDAwMDQwO3JlZGlmZm1haWwuY29t