key: cord-1022630-ojuo3nfo
authors: Barros, Romulo O.; Junior, Fabio L. C. C.; Pereira, Wildrimak S.; Oliveira, Neiva M. N.; Ramos, Ricardo M.
title: Interaction of Drug Candidates with Various SARS-CoV-2 Receptors: An in Silico Study to Combat COVID-19
date: 2020-07-27
journal: J Proteome Res
DOI: 10.1021/acs.jproteome.0c00327
sha: c9271f04add2128f7e5c1356854a18a0c0373c17
doc_id: 1022630
cord_uid: ojuo3nfo

[Image: see text] The world is currently facing the COVID-19 pandemic caused by the SARS-CoV-2 virus. The pandemic is causing the death of people around the world, and public and social health measures to slow or prevent the spread of COVID-19 are being implemented with the involvement of all members of society. Research institutions are accelerating the discovery of vaccines and therapies for COVID-19. In this work, molecular docking was used to study (in silico) the interaction of 24 ligands, divided into four groups, with four SARS-CoV-2 receptors, Nsp9 replicase, main protease (Mpro), NSP15 endoribonuclease, and spike protein (S-protein) interacting with human ACE2. The results showed that the antimalarial drug Metaquine and anti-HIV antiretroviral Saquinavir interacted with all the studied receptors, indicating that they are potential candidates for multitarget drugs for COVID-19.

The end of 2019 and the beginning of 2020 were marked by the pandemic caused by the severe acute respiratory syndrome coronavirus 2. The first cases of infection by a novel coronavirus (SARS-CoV-2) emerged in December 2019 and were related to exposure to the Huanan Seafood Wholesale Market in the city of Wuhan, Hubei, China. 1, 2 The COVID-19 then spread worldwide, and by July 14, 2020, there were 12 964 809 confirmed cases and 570 288 deaths globally. 3 Like SARS-CoV, SARS-CoV-2 seems to use the ACE2 receptor to enter the target cells. 4−8 Alveolar epithelial type II cells are 83% of ACE2-expressing cells, 5 which may explain the higher damage of the infection caused to the lungs, where COVID-19 showed high expression among alveolar epithelial in immunostaining study. 9 There is also evidence that a furinlike cleavage site in the S-protein of SARS-CoV-2 may be implicated in viral invasion of the cells. 10 In most severe cases of COVID-19, patients require endotracheal intubation and mechanical ventilation, as they may develop severe pneumonia, pulmonary edema, ARDS, sepsis shock, or multiple organ failure. 2−11 The disease affects most severely males more than females, and patients with chronic cardiovascular and cerebrovascular diseases, diabetes, obesity, smoking history, and other comorbidities, resulting in lower immunity, have the worst prognostics. 12−14 Research institutions are accelerating the discovery of vaccines and therapies for COVID-19. In this work, molecular docking was used to study (in silico) the interaction of 24 ligands with four SARS-CoV-2 receptors. Crystallographic structures of SARS-CoV-2 receptors, Nsp9 replicase protein, main protease (Mpro), and NSP15 endoribonuclease were used. One of the receptors was the model built by homology modeling of the spike protein (S-protein) and the human ACE2 receptor, and obtained in a study by Smith and Smith, 2020 . 15 The objective of this work was to study the interactions performed by the ligands on these receptors and to indicate drug candidates for COVID-19. and Smith, 2020. 15 Group 3 is made up of five anti-HIV antiretrovirals, and Group 04 is made up of ten ligands from the database of our research laboratory used in previous studies ( Table 2) . Table 1 shows that Receptors 1, 2, and 3 were obtained from the Protein Data Bank (PDB) (https://www.rcsb.org/ ). 16 Receptor 4 was built by homology modeling and obtained from ref 15.

Molecular docking calculations were performed using AutoDock Vina software (https://vina.scripps.edu/). 17 The ligands and proteins were prepared for the calculations with AutoDock Tools (ADT) 1.5.6. 18 Hydrogens were added to both receptors and ligands individually. Gasteiger charges were then calculated by ADT, and nonpolar hydrogens were merged. The size of the grid box was set to 22.5 Å for each axis. The grid boxes were centered on the coordinates of the atoms of residues located in the region of the active site and interface region, as shown in Table 1 . The number of modes was set to 50, and the exhaustiveness was set to 24. A total of 96 molecular docking calculations were performed. Binding energy (ΔG bind ) values better or equal to −7 kcal/mol were established as a criterion for the efficiency of the interaction. Table 2 shows the results of the molecular docking calculations and references related to the ligands. The ligands in Groups 1 and 3 were selected for a drug repurposing study for COVID-19, the ligands in Group 2 were selected to analyze the interaction in the receptors of this study, and the ligands in Group 4 are chemical structures of natural compounds.

Regarding the ligands of Group 1, Table 2 shows that hydroxychloroquine and azithromycin did not interact with the receptors, which may indicate that these drugs act in different steps in the viral cycle other than the ones that involve the proteins in this study. The literature cites three probable mechanisms of action for hydroxychloroquine: conformational modification of angiotensin-converting enzyme type 2 (ACE 2), alteration of the pH of the endosome, and reduction and inhibition in the release of proinflammatory cytokines (TNF-α and IL-6). Mefloquine showed interaction at Receptors 2 (−7.6 kcal/mol) and 3 (−7.3 kcal/mol). Nelfinavir and PL-69 showed interaction in three of the four studied receptors. Metaquine (N,N′-Bis(7chloro-4-quinolyl)-m-phenylenediamine), mentioned in the literature as an antimalarial substance, 26 presented relevant results of interaction in all receptors, mainly in Receptor 3, where the binding energy was equal to −8.6 kcal/mol (see Figure 3 ). The ligands in Group 02 showed relevant results in receptor 4 according to the study the of Smith and Smith, 2020. 15 In this study, Pemirolast and Benserazide ligands did not show any relevant interaction with the other three receptors. Luteolin showed an interaction with a binding energy value of −7.3 kcal/mol with Receptors 2 and 3. The results obtained for the three ligands at Receptor 4 (−7.4 kcal/mol) are in accordance with the study of Smith and Smith, 2020. 15 Among Group 3 ligands, Lopinavir and Ritonavir showed interaction at two of the four receptors, Piperaquine and Dolutegravir had better results and showed interactions with three receptors. Saquinavir showed interaction with all receptors, with highlights to binding energy value of −8.6 kcal/mol at Receptor 4 (see Figure 4) .

Regarding the ligands in Group 4, these did not present relevant results. What stood out the most in this group was Fulvic Acid, with a result of binding energy equal to −7 kcal/ mol with Receptors 2 and 4. Epicatechin showed the same result with Receptor 2. The results show that Saquinavir (Group 1) and Metaquine (Group 3) showed interaction with all the receptors in this study, indicating that these two drugs can be repurposed for the treatment of COVID-19. Nelfinavir and PL-69 (Group 1), Luteolin (Group 2), and Piperaquine and Dolutegravir (Group 3) interacted with three receptors and cannot be excluded as multitarget drug candidates for COVID-19. In silico studies can be carried out by changing the geometry of these ligands in order to improve the binding energy in the studied receptors; besides, in vitro and in vivo studies must be carried out to make promising drug candidates.

Metaquine is a potent, soluble, and bioavailable antimalarial substance. Metaquin shows activity in vivo (oral ID50 of 25 μmol/kg) against Plasmodium berghei and in vitro (0.17 μM) against Plasmodium falciparum K1 multidrug-resistant. Metaquin shows strong affinity for the putative heme antimalarial receptor. 26 Fielding et al. 27 showed that the inclusion of trifluoromethyl (CF3) in the chemical structure of Metaquine had a great impact on the drug interaction with the heme antimalarial receptor.

Saquinavir is a specific inhibitor drug for HIV protease. Jayaswal et al. 34 carried out an in silico study with modified Saquinavir structures to test the interaction with the HIV protease active site, with the aim of finding structures that present better binding energy than Saquinavir. Khadim et al. 33 conducted an in silico study of the interaction of SARS-Cov-2 main protease (Mpro) with anti-HIV drugs. The results obtained are in agreement with our study indicating that Saquinavir presented interactions with Mpro. Figures 1−4 show information about the complexes with the best binding energy for each receptor. Regarding the LigPlot+ diagrams, Figures 1C, 2C , 3C, and 4C, the hydrogen bonds are represented by green dotted lines and the hydrophobic contacts are represented by red semicircles.

The value that appears on the hydrogen bonds corresponds to the distance, in angstroms, between the amino acid residue of the protein and the ligand. Figures S1, S2 , and S3 show interaction of Saquinavir at Receptors 1, 2, and 3, respectively, and Figure S4 shows the interaction of Metaquine at Receptor 4. Figure 1C shows that Metaquine performed two interactions by hydrogen bonds with Receptor 1 (residues Asn28 and Thr78) and 11 interactions by hydrophobic contacts (residues Phe76, Val111, Leu46, Leu30, Asp27, Ala29, Asp48, Thr80, Ser47, Lys87, and Leu107). Figure 2A shows the conformation of Metaquine in a cavity of Receptor 2 near the Asn142 residue (in red) in the region of the active site (see Table 1 ). In Figure 2C , it is shown that Metaquine interacted by hydrophobic contacts with 15 residues (Asn142−active site, Glu166, Phe140, Leu141, Leu27, His163, Gly143, Cys145, Thr25, Thr45, His41, Met49, Ser46, Gln189, and Cys44). Figure 3A shows the conformation of the Metaquine in a cavity of Receptor 3 near His250 (in red) in the region of the active site.

In Receptor 3, Metaquine performed an interaction by hydrogen bond (Ser294) and nine interactions by hydrophobic contacts (Val292, Lys290, His250−active site, Thr341, Trp333, His235, Tyr343, Cys293, and Leu346). Figure 4 shows a part of the model created for human ACE2 receptor in green and part of the model created for spike protein (Sprotein) in blue.

These parts of the two models were selected for molecular docking studies (see ref 15 ). As it can be seen in Figure 4A , Journal of Proteome Research pubs.acs.org/jpr Article Saquinavir is positioned in the expected region, between the ACE2 receptor and the S-protein. Figure 4C shows that Saquinavir performed four interactions by hydrogen bond, three of them with the ACE2 receptor (Asn33, Glu37, and Phe390) and one with the S-protein (Lys417). In addition, it performed 18 interactions by hydrophobic contact with both models.

The results confirmed the interaction of Metaquine (N,N′-Bis(7-chloro-4-quinolyl)-m-phenylenediamine) and Saquinavir in all receptors studied, indicating them as candidates for multitarget drugs that can be repurposed to the treatment of COVID-19. In vitro, in vivo, and clinical tests should be performed to confirm the effectiveness of these drugs in the treatment of COVID-19.

■ ASSOCIATED CONTENT

The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.jproteome.0c00327.

Figures S1−S4: Global structure of the conformations resulted from docking (PDF)

Early Transmission Dynamics in Wuhan, China, of Novel Coronavirus−Infected Pneumonia

Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: a descriptive study

Situation Report

Angiotensin-converting enzyme 2 (ACE2) as a SARS-CoV-2 receptor: molecular mechanisms and potential therapeutic target

Singlecell RNA expression profiling of ACE2, the putative receptor of Wuhan 2019-nCov. bioRxiv

A pneumonia outbreak associated with a new coronavirus of probable bat origin

Evolution of the novel coronavirus from the ongoing Wuhan outbreak and modeling of its spike protein for risk of human transmission

SARS-CoV-2 Cell Entry Depends on ACE2 and TMPRSS2 and Is Blocked by a Clinically Proven Protease Inhibitor

The spike glycoprotein of the new coronavirus 2019-nCoV contains a furin-like cleavage site absent in CoV of the same clade

Review of the Clinical Characteristics of Coronavirus Disease 2019 (COVID-19)

Novel Coronavirus (COVID-19) based on current evidence

Coronavirus Disease (COVID-19): A primer for emergency physicians

Prevalence of comorbidities in the novel Wuhan coronavirus (COVID-19) infection: a systematic review and meta-analysis

Repurposing Therapeutics for COVID-19: Supercomputer-Based Docking to the SARS-CoV-2 Viral Spike Protein and Viral Spike Protein-Human ACE2 Interface

The Protein Data

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

AutoDock4 and AutoDockTools4: Automated Docking with Selective Receptor Flexibility

Hydroxychloroquine and azithromycin as a treatment of COVID-19: results of an open-label non-randomized clinical trial

Chloroquine and hydroxychloroquine in the treatment of COVID-19 with or without diabetes: A systematic search and a narrative review with a special reference to India and other developing countries

Azithromycin treats diffuse panbronchiolitis by Journal of Proteome Research pubs

targeting T cells via inhibition of mTOR pathway

Adjunctive therapy with azithromycin for moderate and severe acute respiratory distress syndrome: a retrospective, propensity score-matching analysis of prospectively collected data at a single center

Spiramycin and azithromycin, safe for administration to children, exert antiviral activity against enterovirus A71 in vitro and in vivo

Ontological and bioinformatic analysis of anti-coronavirus drugs and their implication for drug repurposing against COVID-19. Preprints

Mapping Antimalarial Pharmacophores as a Useful Tool for the Rapid Discovery of Drugs Effective in Vivo: Design, Construction, Characterization, and Pharmacology of Metaquine

Modulation of Antimalarial Activity at a Putative Bisquinoline Receptor In Vivo Using Fluorinated Bisquinolines

Virtual Screening based prediction of potential drugs for COVID-19. Preprints

Nelfinavir was predicted to be a potential inhibitor of 2019-nCov main protease by an integrative approach combining homology modelling, molecular docking and binding free energy calculation. bioRxiv

Synthesis, structure-activity relationship, and mode-of-action studies of antimalarial reversed chloroquine compounds

Incorporation of graphene oxide−NiO nanocomposite and n-hexyl-3-methylimidazolium hexafluoro phosphate into carbon paste electrode: application as an electrochemical sensor for simultaneous determination of benserazide, levodopa and tryptophan

Luteolin sensitizes the antitumor effect of cisplatin in drug-resistant ovarian cancer via induction of apoptosis and inhibition of cell migration and invasion

pH Dependent Differential Binding Behavior of Prtotease Inhibitor Molecular Drugs for SARS-COV-2. ChemRxiv

Examining pharmacodynamic and pharmacokinetic properties of eleven analogues of saquinavir for HIV protease inhibition

Recent Progress in understanding 2019 Novel Coronavirus associated with Human Respiratory Disease: Detection, Mechanism and Treatment

Patients of COVID-19 may benefit from sustained lopinavir-combined regimen and the increase of eosinophil may predict the outcome of COVID-19 progression

Case of the Index Patient Who Caused Tertiary Transmission of Coronavirus Disease 2019 in Korea: the Application of Lopinavir/Ritonavir for the Treatment of COVID-19 Pneumonia Monitored by Quantitative RT-PCR

Intermittent screening and treatment or intermittent western Kenya: an open-label, three-group, randomised controlled superiority trial

Dolutegravir plus Abacavir−Lamivudine for the Treatment of HIV-1 Infection

Carbohydrate-derived fulvic acid is a highly promising topical agent to enhance healing of wounds infected with drug-resistant pathogens

Effects of macromolecular humic/fulvic acid on Cd (II) adsorption onto reed-derived biochar as compared with tannic acid

Characterization of the dissolution behavior of piperine/cyclodextrins inclusion complexes

De novo transcriptome analysis of rose-scented geranium provides insights into the metabolic specificity of terpene and tartaric acid biosynthesis

Linalool-and α-terpineol-induced programmed cell death in Chlamydomonas reinhardtii

Factors affecting levels of ferulic acid, ethyl ferulate and taste-active pyroglutamyl peptides in sake

Epicatechin alleviates inflammation in lipopolysaccharide-induced acute lung injury in mice by inhibiting the p38 MAPK signaling pathway

Epicatechin protective effects on bleomycin-induced pulmonary oxidative stress and fibrosis in mice

Pires de Moura do Amaral, F.; et al. Phytochemical study, molecular docking, genotoxicity and therapeutic efficacy of the aqueous extract of the stem bark of Ximenia americana L. in the treatment of experimental COPD in rats

The capsaicin receptor TRPV1 is the first line defense protecting Journal of Proteome Research pubs

from acute non damaging heat: a translational approach

Ceśar, I. C. Liquid chromatographytandem mass spectrometry bioanalytical method for the determination of kavain in mice plasma: Application to a pharmacokinetic study

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

UCSF Chimeraa visualization system for exploratory research and analysis