51
|
Maurina SF, O'Sullivan JP, Sharma G, Pineda Rodriguez DC, MacFadden A, Cendali F, Henen MA, Vögeli B, Kieft JS, Glasgow A, Steckelberg AL. An Evolutionarily Conserved Strategy for Ribosome Binding and Host Translation Inhibition by β-coronavirus Non-structural Protein 1. J Mol Biol 2023; 435:168259. [PMID: 37660941 PMCID: PMC10543557 DOI: 10.1016/j.jmb.2023.168259] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2023] [Revised: 08/15/2023] [Accepted: 08/28/2023] [Indexed: 09/05/2023]
Abstract
An important pathogenicity factor of SARS-CoV-2 and related coronaviruses is Non-structural protein 1 (Nsp1), which suppresses host gene expression and stunts antiviral signaling. SARS-CoV-2 Nsp1 binds the ribosome to inhibit translation through mRNA displacement and induces degradation of host mRNAs. Here we show that Nsp1-dependent host shutoff is conserved in diverse coronaviruses, but only Nsp1 from β-Coronaviruses (β-CoV) inhibits translation through ribosome binding. The C-terminal domain (CTD) of all β-CoV Nsp1s confers high-affinity ribosome binding despite low sequence conservation. Modeling of interactions of four Nsp1s with the ribosome identified the few absolutely conserved amino acids that, together with an overall conservation in surface charge, form the β-CoV Nsp1 ribosome-binding domain. Contrary to previous models, the Nsp1 ribosome-binding domain is an inefficient translation inhibitor. Instead, the Nsp1-CTD likely functions by recruiting Nsp1's N-terminal "effector" domain. Finally, we show that a cis-acting viral RNA element has co-evolved to fine-tune SARS-CoV-2 Nsp1 function, but does not provide similar protection against Nsp1 from related viruses. Together, our work provides new insight into the diversity and conservation of ribosome-dependent host-shutoff functions of Nsp1, knowledge that could aid future efforts in pharmacological targeting of Nsp1 from SARS-CoV-2 and related human-pathogenic β-CoVs. Our study also exemplifies how comparing highly divergent Nsp1 variants can help to dissect the different modalities of this multi-functional viral protein.
Collapse
|
52
|
Xie L, Lin F, Dong P, Li Y. MAb Targeting a Link Between ExoN and MTase of TGEV NSP14. Monoclon Antib Immunodiagn Immunother 2023; 42:178-181. [PMID: 37855908 PMCID: PMC10621669 DOI: 10.1089/mab.2023.0010] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2023] [Accepted: 09/29/2023] [Indexed: 10/20/2023] Open
Abstract
Porcine transmissible gastroenteritis virus (TGEV) infection results in severe gastrointestinal disease manifesting vomiting, diarrhea in neonatal porcine, with extremely high mortality. Monoclonal antibody (MAb) specific to TGEV nonstructural protein (NSP)14 that contains two functional domains, exonuclease (ExoN) and methyltransferase (MTase) domains, may help elucidate the role of NSP14 in the viral life-cycle. In this study, we developed a murine MAb, designated 12F1, against TGEV NSP14 using traditional cell-fusion technique. It was shown the MAb can exclusively bind to viral NSP14, as evidenced by the results of indirect fluorescent assay and western blotting. Intriguingly, epitope screening assay shown that 12F1 targets a hinge region connecting ExoN and N7-MTase of NSP14.
Collapse
|
53
|
Andrade MA, Mottin M, Sousa BKDP, Barbosa JARG, Dos Santos Azevedo C, Lasse Silva C, Gonçalves de Andrade M, Motta FN, Maulay-Bailly C, Amand S, Santana JMD, Horta Andrade C, Grellier P, Bastos IMD. Identification of novel Zika virus NS3 protease inhibitors with different inhibition modes by integrative experimental and computational approaches. Biochimie 2023; 212:143-152. [PMID: 37088408 DOI: 10.1016/j.biochi.2023.04.004] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2022] [Revised: 02/14/2023] [Accepted: 04/07/2023] [Indexed: 04/25/2023]
Abstract
Zika virus (ZIKV) infection is associated with severe neurological disorders and congenital malformation. Despite efforts to eradicate the disease, there is still neither vaccine nor approved drugs to treat ZIKV infection. The NS2B-NS3 protease is a validated drug target since it is essential to polyprotein virus maturation. In the present study, we describe an experimental screening of 2,320 compounds from the chemical library of the Muséum National d'Histoire Naturelle of Paris on ZIKV NS2B-NS3 protease. A total of 96 hits were identified with 90% or more of inhibitory activity at 10 μM. Amongst the most active compounds, five were analyzed for their inhibitory mechanisms by kinetics assays and computational approaches such as molecular docking. 2-(3-methoxyphenoxy) benzoic acid (compound 945) show characteristics of a competitive inhibition (Ki = 0.49 μM) that was corroborated by its molecular docking at the active site of the NS2B-NS3 protease. Taxifolin (compound 2292) behaves as an allosteric inhibitor whereas 3,8,9-trihydroxy-2-methyl-1H-phenalen-1-one (compound 128), harmol (compound 368) and anthrapurpurin (compound 1499) show uncompetitive inhibitions. These new NS2B-NS3 protease inhibitors are valuable hits to further hit-to-lead optimization.
Collapse
|
54
|
Abaeva IS, Arhab Y, Miścicka A, Hellen CUT, Pestova TV. In vitro reconstitution of SARS-CoV-2 Nsp1-induced mRNA cleavage reveals the key roles of the N-terminal domain of Nsp1 and the RRM domain of eIF3g. Genes Dev 2023; 37:844-860. [PMID: 37821106 PMCID: PMC10620056 DOI: 10.1101/gad.350829.123] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2023] [Accepted: 09/19/2023] [Indexed: 10/13/2023]
Abstract
SARS CoV-2 nonstructural protein 1 (Nsp1) is the major pathogenesis factor that inhibits host translation using a dual strategy of impairing initiation and inducing endonucleolytic cleavage of cellular mRNAs. To investigate the mechanism of cleavage, we reconstituted it in vitro on β-globin, EMCV IRES, and CrPV IRES mRNAs that use unrelated initiation mechanisms. In all instances, cleavage required Nsp1 and only canonical translational components (40S subunits and initiation factors), arguing against involvement of a putative cellular RNA endonuclease. Requirements for initiation factors differed for these mRNAs, reflecting their requirements for ribosomal attachment. Cleavage of CrPV IRES mRNA was supported by a minimal set of components consisting of 40S subunits and eIF3g's RRM domain. The cleavage site was located in the coding region 18 nt downstream from the mRNA entrance, indicating that cleavage occurs on the solvent side of the 40S subunit. Mutational analysis identified a positively charged surface on Nsp1's N-terminal domain (NTD) and a surface above the mRNA-binding channel on eIF3g's RRM domain that contain residues essential for cleavage. These residues were required for cleavage on all three mRNAs, highlighting general roles of the Nsp1 NTD and eIF3g's RRM domain in cleavage per se, irrespective of the mode of ribosomal attachment.
Collapse
|
55
|
Fink EA, Bardine C, Gahbauer S, Singh I, Detomasi TC, White K, Gu S, Wan X, Chen J, Ary B, Glenn I, O'Connell J, O'Donnell H, Fajtová P, Lyu J, Vigneron S, Young NJ, Kondratov IS, Alisoltani A, Simons LM, Lorenzo‐Redondo R, Ozer EA, Hultquist JF, O'Donoghue AJ, Moroz YS, Taunton J, Renslo AR, Irwin JJ, García‐Sastre A, Shoichet BK, Craik CS. Large library docking for novel SARS-CoV-2 main protease non-covalent and covalent inhibitors. Protein Sci 2023; 32:e4712. [PMID: 37354015 PMCID: PMC10364469 DOI: 10.1002/pro.4712] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2023] [Revised: 05/29/2023] [Accepted: 06/21/2023] [Indexed: 06/25/2023]
Abstract
Antiviral therapeutics to treat SARS-CoV-2 are needed to diminish the morbidity of the ongoing COVID-19 pandemic. A well-precedented drug target is the main viral protease (MPro ), which is targeted by an approved drug and by several investigational drugs. Emerging viral resistance has made new inhibitor chemotypes more pressing. Adopting a structure-based approach, we docked 1.2 billion non-covalent lead-like molecules and a new library of 6.5 million electrophiles against the enzyme structure. From these, 29 non-covalent and 11 covalent inhibitors were identified in 37 series, the most potent having an IC50 of 29 and 20 μM, respectively. Several series were optimized, resulting in low micromolar inhibitors. Subsequent crystallography confirmed the docking predicted binding modes and may template further optimization. While the new chemotypes may aid further optimization of MPro inhibitors for SARS-CoV-2, the modest success rate also reveals weaknesses in our approach for challenging targets like MPro versus other targets where it has been more successful, and versus other structure-based techniques against MPro itself.
Collapse
|
56
|
Godoy AS, Nakamura AM, Douangamath A, Song Y, Noske GD, Gawriljuk VO, Fernandes RS, Pereira H, Oliveira K, Fearon D, Dias A, Krojer T, Fairhead M, Powell A, Dunnet L, Brandao-Neto J, Skyner R, Chalk R, Bajusz D, Bege M, Borbás A, Keserű GM, von Delft F, Oliva G. Allosteric regulation and crystallographic fragment screening of SARS-CoV-2 NSP15 endoribonuclease. Nucleic Acids Res 2023; 51:5255-5270. [PMID: 37115000 PMCID: PMC10250223 DOI: 10.1093/nar/gkad314] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2022] [Revised: 04/11/2023] [Accepted: 04/13/2023] [Indexed: 04/29/2023] Open
Abstract
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of coronavirus disease 2019 (COVID-19). The NSP15 endoribonuclease enzyme, known as NendoU, is highly conserved and plays a critical role in the ability of the virus to evade the immune system. NendoU is a promising target for the development of new antiviral drugs. However, the complexity of the enzyme's structure and kinetics, along with the broad range of recognition sequences and lack of structural complexes, hampers the development of inhibitors. Here, we performed enzymatic characterization of NendoU in its monomeric and hexameric form, showing that hexamers are allosteric enzymes with a positive cooperative index, and with no influence of manganese on enzymatic activity. Through combining cryo-electron microscopy at different pHs, X-ray crystallography and biochemical and structural analysis, we showed that NendoU can shift between open and closed forms, which probably correspond to active and inactive states, respectively. We also explored the possibility of NendoU assembling into larger supramolecular structures and proposed a mechanism for allosteric regulation. In addition, we conducted a large fragment screening campaign against NendoU and identified several new allosteric sites that could be targeted for the development of new inhibitors. Overall, our findings provide insights into the complex structure and function of NendoU and offer new opportunities for the development of inhibitors.
Collapse
|
57
|
Ghavami G, Adibzadeh S, Amiri S, Sardari S. Combined in silico strategy for repurposing DrugBank entries towards introducing potential anti-SARS-CoV-2 drugs. Can J Physiol Pharmacol 2023; 101:268-285. [PMID: 36848647 DOI: 10.1139/cjpp-2022-0309] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/01/2023]
Abstract
The emergence of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) from China in December 2019 led to the coronavirus disorder 2019 pandemic, which has affected tens of millions of humans worldwide. Various in silico research via bio-cheminformatics methods were performed to examine the efficiency of a range of repurposed approved drugs with a new role as anti-SARS-CoV-2 drugs. The current study has been performed to screen the approved drugs in the DrugBank database based on a novel bioinformatics/cheminformatics strategy to repurpose available approved drugs towards introducing them as a possible anti-SARS-CoV-2 drug. As a result, 96 approved drugs with the best docking scores passed through several relevant filters were presented as the candidate drugs with potential novel antiviral activities against the SARS-CoV-2 virus.
Collapse
|
58
|
Roy KK, Jyothi D, Paul U, Sukla S. Identification and validation of novel non-nucleoside class of molecules inhibiting the dengue virus replication. J Biomol Struct Dyn 2023; 41:13993-14002. [PMID: 36970872 DOI: 10.1080/07391102.2023.2192796] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2022] [Accepted: 01/29/2023] [Indexed: 03/29/2023]
Abstract
There is currently no drug available to treat mosquito-borne dengue. The C-terminal RNA-dependent RNA polymerase (RdRp) domain in the non-structural type 5 (NS5) protein of the dengue virus (DENV) is essential for viral RNA synthesis and replication, and therefore, it is an attractive target for the anti-dengue drug development. We report herein the discovery and validation of two novel non-nucleoside classes of small molecules as DENV RdRp inhibitors. Firstly, using the refined X-ray structure of the DENV NS5 RdRp domain (PDB-ID: 4V0R), we conducted docking, binding free-energy studies, and short-scale molecular dynamics simulation to investigate the binding sites of known small molecules that led to the optimized protein-ligand complex. Subsequently, protein structure-based screening of a commercial database (∼500,000 synthetic compounds), pre-filtered for the drug-likeness, led to the top-ranked 171 molecules, which was then subjected to structural diversity analysis and clustering. This process led to six structurally distinct and best-scored compounds that were procured from the commercial vendor, and then subjected to the in vitro testing in the MTT and dengue infection assays. It revealed two unique and structurally unique compounds, KKR-D-02 and KKR-D-03, exhibiting 84 and 81% reductions, respectively, in DENV copy number in repeated assays in comparison to the virus-infected cell controls. These active compounds represent novel scaffolds for further structure-based discovery of novel candidate molecules for the intervention of dengue.Communicated by Ramaswamy H. Sarma.
Collapse
|
59
|
Yang KS, Blankenship LR, Kuo STA, Sheng YJ, Li P, Fierke CA, Russell DH, Yan X, Xu S, Liu WR. A Novel Y-Shaped, S-O-N-O-S-Bridged Cross-Link between Three Residues C22, C44, and K61 Is Frequently Observed in the SARS-CoV-2 Main Protease. ACS Chem Biol 2023; 18:449-455. [PMID: 36629751 PMCID: PMC10023456 DOI: 10.1021/acschembio.2c00695] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
As the COVID-19 pathogen, SARS-CoV-2 relies on its main protease (MPro) for pathogenesis and replication. During crystallographic analyses of MPro crystals that were exposed to the air, a uniquely Y-shaped, S-O-N-O-S-bridged post-translational cross-link that connects three residues C22, C44, and K61 at their side chains was frequently observed. As a novel covalent modification, this cross-link serves potentially as a redox switch to regulate the catalytic activity of MPro, a demonstrated drug target of COVID-19. The formation of this linkage leads to a much more open active site that can potentially be targeted for the development of novel SARS-CoV-2 antivirals. The structural rearrangement of MPro by this cross-link indicates that small molecules that lock MPro in the cross-linked form can potentially be used with other active-site-targeting molecules such as paxlovid for synergistic effects in inhibiting SARS-CoV-2 viral replication.
Collapse
|
60
|
Cheng J, Hao Y, Shi Q, Hou G, Wang Y, Wang Y, Xiao W, Othman J, Qi J, Wang Y, Chen Y, Yu G. Discovery of Novel Chinese Medicine Compounds Targeting 3CL Protease by Virtual Screening and Molecular Dynamics Simulation. Molecules 2023; 28:molecules28030937. [PMID: 36770604 PMCID: PMC9921503 DOI: 10.3390/molecules28030937] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2022] [Revised: 12/23/2022] [Accepted: 01/13/2023] [Indexed: 01/19/2023] Open
Abstract
The transmission and infectivity of COVID-19 have caused a pandemic that has lasted for several years. This is due to the constantly changing variants and subvariants that have evolved rapidly from SARS-CoV-2. To discover drugs with therapeutic potential for COVID-19, we focused on the 3CL protease (3CLpro) of SARS-CoV-2, which has been proven to be an important target for COVID-19 infection. Computational prediction techniques are quick and accurate enough to facilitate the discovery of drugs against the 3CLpro of SARS-CoV-2. In this paper, we used both ligand-based virtual screening and structure-based virtual screening to screen the traditional Chinese medicine small molecules that have the potential to target the 3CLpro of SARS-CoV-2. MD simulations were used to confirm these results for future in vitro testing. MCCS was then used to calculate the normalized free energy of each ligand and the residue energy contribution. As a result, we found ZINC15676170, ZINC09033700, and ZINC12530139 to be the most promising antiviral therapies against the 3CLpro of SARS-CoV-2.
Collapse
|
61
|
Ren J, Vaid TM, Lee H, Ojeda I, Johnson ME. Evaluation of interactions between the hepatitis C virus NS3/4A and sulfonamidobenzamide based molecules using molecular docking, molecular dynamics simulations and binding free energy calculations. J Comput Aided Mol Des 2023; 37:53-65. [PMID: 36427108 PMCID: PMC9839505 DOI: 10.1007/s10822-022-00490-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/20/2022] [Accepted: 11/17/2022] [Indexed: 11/26/2022]
Abstract
The Hepatitis C Virus (HCV) NS3/4A is an attractive target for the treatment of Hepatitis C infection. Herein, we present an investigation of HCV NS3/4A inhibitors based on a sulfonamidobenzamide scaffold. Inhibitor interactions with HCV NS3/4A were explored by molecular docking, molecular dynamics simulations, and MM/PBSA binding free energy calculations. All of the inhibitors adopt similar molecular docking poses in the catalytic site of the protease that are stabilized by hydrogen bond interactions with G137 and the catalytic S139, which are known to be important for potency and binding stability. The quantitative assessments of binding free energies from MM/PBSA correlate well with the experimental results, with a high coefficient of determination, R2 of 0.92. Binding free energy decomposition analyses elucidate the different contributions of Q41, F43, H57, R109, K136, G137, S138, S139, A156, M485, and Q526 in binding different inhibitors. The importance of these sidechain contributions was further confirmed by computational alanine scanning mutagenesis. In addition, the sidechains of K136 and S139 show crucial but distinct contributions to inhibitor binding with HCV NS3/4A. The structural basis of the potency has been elucidated, demonstrating the importance of the R155 sidechain conformation. This extensive exploration of binding energies and interactions between these compounds and HCV NS3/4A at the atomic level should benefit future antiviral drug design.
Collapse
|
62
|
Naceri S, Marc D, Blot R, Flatters D, Camproux AC. Druggable Pockets at the RNA Interface Region of Influenza A Virus NS1 Protein Are Conserved across Sequence Variants from Distinct Subtypes. Biomolecules 2022; 13:biom13010064. [PMID: 36671449 PMCID: PMC9855689 DOI: 10.3390/biom13010064] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2022] [Revised: 12/24/2022] [Accepted: 12/25/2022] [Indexed: 12/31/2022] Open
Abstract
Influenza A viruses still represent a major health issue, for both humans and animals. One of the main viral proteins of interest to target is the NS1 protein, which counters the host immune response and promotes viral replication. NS1 is a homodimer composed of a dimeric RNA-binding domain (RBD), which is structurally stable and conserved in sequence, and two effector domains that are tethered to the RBD by linker regions. This linker flexibility leads to NS1 polymorphism and can therefore exhibit different forms. Previously, we identified a putative drug-binding site, located in the RBD interface in a crystal structure of NS1. This pocket could be targeted to block RNA binding and inhibit NS1 activities. The objective of the present study is to confirm the presence of this druggable site, whatever the sequence variants, in order to develop a universal therapeutic compound that is insensitive to sequence variations and structural flexibility. Using a set of four NS1 full-length structures, we combined different bioinformatics approaches such as pocket tracking along molecular dynamics simulations, druggability prediction and classification. This protocol successfully confirmed a frequent large binding-site that is highly druggable and shared by different NS1 forms, which is promising for developing a robust NS1-targeted therapy.
Collapse
|
63
|
Mottin M, de Paula Sousa BK, de Moraes Roso Mesquita NC, de Oliveira KIZ, Noske GD, Sartori GR, de Oliveira Albuquerque A, Urbina F, Puhl AC, Moreira-Filho JT, Souza GE, Guido RV, Muratov E, Neves BJ, da Silva JHM, Clark AE, Siqueira-Neto JL, Perryman AL, Oliva G, Ekins S, Andrade CH. Discovery of New Zika Protease and Polymerase Inhibitors through the Open Science Collaboration Project OpenZika. J Chem Inf Model 2022; 62:6825-6843. [PMID: 36239304 PMCID: PMC9923514 DOI: 10.1021/acs.jcim.2c00596] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
The Zika virus (ZIKV) is a neurotropic arbovirus considered a global threat to public health. Although there have been several efforts in drug discovery projects for ZIKV in recent years, there are still no antiviral drugs approved to date. Here, we describe the results of a global collaborative crowdsourced open science project, the OpenZika project, from IBM's World Community Grid (WCG), which integrates different computational and experimental strategies for advancing a drug candidate for ZIKV. Initially, molecular docking protocols were developed to identify potential inhibitors of ZIKV NS5 RNA-dependent RNA polymerase (NS5 RdRp), NS3 protease (NS2B-NS3pro), and NS3 helicase (NS3hel). Then, a machine learning (ML) model was built to distinguish active vs inactive compounds for the cytoprotective effect against ZIKV infection. We performed three independent target-based virtual screening campaigns (NS5 RdRp, NS2B-NS3pro, and NS3hel), followed by predictions by the ML model and other filters, and prioritized a total of 61 compounds for further testing in enzymatic and phenotypic assays. This yielded five non-nucleoside compounds which showed inhibitory activity against ZIKV NS5 RdRp in enzymatic assays (IC50 range from 0.61 to 17 μM). Two compounds thermally destabilized NS3hel and showed binding affinity in the micromolar range (Kd range from 9 to 35 μM). Moreover, the compounds LabMol-301 inhibited both NS5 RdRp and NS2B-NS3pro (IC50 of 0.8 and 7.4 μM, respectively) and LabMol-212 thermally destabilized the ZIKV NS3hel (Kd of 35 μM). Both also protected cells from death induced by ZIKV infection in in vitro cell-based assays. However, while eight compounds (including LabMol-301 and LabMol-212) showed a cytoprotective effect and prevented ZIKV-induced cell death, agreeing with our ML model for prediction of this cytoprotective effect, no compound showed a direct antiviral effect against ZIKV. Thus, the new scaffolds discovered here are promising hits for future structural optimization and for advancing the discovery of further drug candidates for ZIKV. Furthermore, this work has demonstrated the importance of the integration of computational and experimental approaches, as well as the potential of large-scale collaborative networks to advance drug discovery projects for neglected diseases and emerging viruses, despite the lack of available direct antiviral activity and cytoprotective effect data, that reflects on the assertiveness of the computational predictions. The importance of these efforts rests with the need to be prepared for future viral epidemic and pandemic outbreaks.
Collapse
|
64
|
Shi L, Wen Z, Song Y, Wang J, Yu D. Computational investigation of potent inhibitors against SARS-CoV-2 2'-O-methyltransferase (nsp16): Structure-based pharmacophore modeling, molecular docking, molecular dynamics simulations and binding free energy calculations. J Mol Graph Model 2022; 117:108306. [PMID: 36063745 PMCID: PMC9385381 DOI: 10.1016/j.jmgm.2022.108306] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/24/2021] [Revised: 08/06/2022] [Accepted: 08/08/2022] [Indexed: 01/14/2023]
Abstract
The Coronavirus Disease 2019 (COVID-19), caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, has created unprecedented public health and economic crises around the world. SARS-CoV-2 2'-O-methyltransferase (nsp16) adds a "cap" to viral RNA to maintain the stability of viral RNA, and inhibition of nsp16 activity may reduce viral proliferation, making this protein an attractive drug target. Here, we report the identification of several small molecule inhibitors of nsp16 by virtual screening. First, the nsp16-sinefungin complex (PDB ID: 6WKQ) was selected from the protein data bank. Asp6912, Cys6913, Asp6897 and Asp6928 were determined to be the key amino acids for sinefungin binding in the crystal structure of nsp16-sinefungin complex by molecular dynamics simulation. The complex structures in the stable binding trajectory of nsp16-sinefungin were than clustered through molecular dynamics RMSD analysis. Six clusters were generated, and six representative structures were selected to construct the pharmacophore based on the structure. These six pharmacophores were superimposed on the binding pocket to simplify and pick the common characteristics. The compounds obtained by the pharmacophore screening from Bionet and Chembiv databases were docked into the nsp16 active pocket. The candidate compounds were selected according to the molecular docking score and then screened by MM/GBSA. Finally, four candidate compounds were obtained. Four sets of 150ns molecular dynamics simulations were performed to determine whether candidate compounds could maintain stable interactions with key amino acids. The results of MD and MM/PBSA energy decomposition indicated that C1 and C2 could form a stable complex system with nsp16, and could form strong hydrogen bonds and salt bridges with the key amino acid Asp6897 and Asp6928. This study thus identifies and attempts to validate for the first time the potential inhibitory activities of C1 and C2 against nsp16, allowing the development of potent anti-COVID-19 drugs and unique treatment strategies.
Collapse
|
65
|
Di Matteo F, Frumenzio G, Chandramouli B, Grottesi A, Emerson A, Musiani F. Computational Study of Helicase from SARS-CoV-2 in RNA-Free and Engaged Form. Int J Mol Sci 2022; 23:ijms232314721. [PMID: 36499049 PMCID: PMC9738952 DOI: 10.3390/ijms232314721] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2022] [Revised: 11/16/2022] [Accepted: 11/22/2022] [Indexed: 11/27/2022] Open
Abstract
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the pandemic that broke out in 2020 and continues to be the cause of massive global upheaval. Coronaviruses are positive-strand RNA viruses with a genome of ~30 kb. The genome is replicated and transcribed by RNA-dependent RNA polymerase together with accessory factors. One of the latter is the protein helicase (NSP13), which is essential for viral replication. The recently solved helicase structure revealed a tertiary structure composed of five domains. Here, we investigated NSP13 from a structural point of view, comparing its RNA-free form with the RNA-engaged form by using atomistic molecular dynamics (MD) simulations at the microsecond timescale. Structural analyses revealed conformational changes that provide insights into the contribution of the different domains, identifying the residues responsible for domain-domain interactions in both observed forms. The RNA-free system appears to be more flexible than the RNA-engaged form. This result underlies the stabilizing role of the nucleic acid and the functional core role of these domains.
Collapse
|
66
|
Mondal S, Chen Y, Lockbaum GJ, Sen S, Chaudhuri S, Reyes AC, Lee JM, Kaur AN, Sultana N, Cameron MD, Shaffer SA, Schiffer CA, Fitzgerald KA, Thompson PR. Dual Inhibitors of Main Protease (M Pro) and Cathepsin L as Potent Antivirals against SARS-CoV2. J Am Chem Soc 2022; 144:21035-21045. [PMID: 36356199 PMCID: PMC9662648 DOI: 10.1021/jacs.2c04626] [Citation(s) in RCA: 18] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2022] [Indexed: 11/12/2022]
Abstract
Given the current impact of SARS-CoV2 and COVID-19 on human health and the global economy, the development of direct acting antivirals is of paramount importance. Main protease (MPro), a cysteine protease that cleaves the viral polyprotein, is essential for viral replication. Therefore, MPro is a novel therapeutic target. We identified two novel MPro inhibitors, D-FFRCMKyne and D-FFCitCMKyne, that covalently modify the active site cysteine (C145) and determined cocrystal structures. Medicinal chemistry efforts led to SM141 and SM142, which adopt a unique binding mode within the MPro active site. Notably, these inhibitors do not inhibit the other cysteine protease, papain-like protease (PLPro), involved in the life cycle of SARS-CoV2. SM141 and SM142 block SARS-CoV2 replication in hACE2 expressing A549 cells with IC50 values of 8.2 and 14.7 nM. Detailed studies indicate that these compounds also inhibit cathepsin L (CatL), which cleaves the viral S protein to promote viral entry into host cells. Detailed biochemical, proteomic, and knockdown studies indicate that the antiviral activity of SM141 and SM142 results from the dual inhibition of MPro and CatL. Notably, intranasal and intraperitoneal administration of SM141 and SM142 lead to reduced viral replication, viral loads in the lung, and enhanced survival in SARS-CoV2 infected K18-ACE2 transgenic mice. In total, these data indicate that SM141 and SM142 represent promising scaffolds on which to develop antiviral drugs against SARS-CoV2.
Collapse
|
67
|
Ornelas MY, Thomas AY, Johnson Rosas LI, Scoville RO, Mehta AP. Synthetic Platforms for Characterizing and Targeting of SARS-CoV-2 Genome Capping Enzymes. ACS Synth Biol 2022; 11:3759-3771. [PMID: 36331143 PMCID: PMC9662071 DOI: 10.1021/acssynbio.2c00359] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2022] [Indexed: 11/06/2022]
Abstract
Essential viral enzymes have been successfully targeted to combat the diseases caused by emerging pathogenic RNA viruses (e.g., viral RNA-dependent RNA polymerase). Because of the conserved nature of such viral enzymes, therapeutics targeting these enzymes have the potential to be repurposed to combat emerging diseases, e.g., remdesivir, which was initially developed as a potential Ebola treatment, then was repurposed for COVID-19. Our efforts described in this study target another essential and highly conserved, but relatively less explored, step in RNA virus translation and replication, i.e., capping of the viral RNA genome. The viral genome cap structure disguises the genome of most RNA viruses to resemble the mRNA cap structure of their host and is essential for viral translation, propagation, and immune evasion. Here, we developed a synthetic, phenotypic yeast-based complementation platform (YeRC0M) for molecular characterization and targeting of SARS-CoV-2 genome-encoded RNA cap-0 (guanine-N7)-methyltransferase (N7-MTase) enzyme (nsp14). In YeRC0M, the lack of yeast mRNA capping N7-MTase in yeast, which is an essential gene in yeast, is complemented by the expression of functional viral N7-MTase or its variants. Using YeRC0M, we first identified important protein domains and amino acid residues that are essential for SARS-CoV-2 nsp14 N7-MTase activity. We also expanded YeRC0M to include key nsp14 variants observed in emerging variants of SARS-CoV-2 (e.g., delta variant of SARS-CoV-2 encodes nsp14 A394V and nsp14 P46L). We also combined YeRC0M with directed evolution to identify attenuation mutations in SARS-CoV-2 nsp14. Because of the high sequence similarity of nsp14 in emerging coronaviruses, these observations could have implications on live attenuated vaccine development strategies. These data taken together reveal key domains in SARS-CoV-2 nsp14 that can be targeted for therapeutic strategies. We also anticipate that these readily tractable phenotypic platforms can also be used for the identification of inhibitors of viral RNA capping enzymes as antivirals.
Collapse
|
68
|
Padhi AK, Tripathi T. Hotspot residues and resistance mutations in the nirmatrelvir-binding site of SARS-CoV-2 main protease: Design, identification, and correlation with globally circulating viral genomes. Biochem Biophys Res Commun 2022; 629:54-60. [PMID: 36113178 PMCID: PMC9450486 DOI: 10.1016/j.bbrc.2022.09.010] [Citation(s) in RCA: 18] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2022] [Accepted: 09/02/2022] [Indexed: 12/15/2022]
Abstract
Shortly after the onset of the COVID-19 pandemic, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has acquired numerous variations in its intracellular proteins to adapt quickly, become more infectious, and ultimately develop drug resistance by mutating certain hotspot residues. To keep the emerging variants at bay, including Omicron and subvariants, FDA has approved the antiviral nirmatrelvir for mild-to-moderate and high-risk COVID-19 cases. Like other viruses, SARS-CoV-2 could acquire mutations in its main protease (Mpro) to adapt and develop resistance against nirmatrelvir. Employing a unique high-throughput protein design technique, the hotspot residues, and signatures of adaptation of Mpro having the highest probability of mutating and rendering nirmatrelvir ineffective were identified. Our results show that ∼40% of the designed mutations in Mpro already exist in the globally circulating SARS-CoV-2 lineages and several predicted mutations. Moreover, several high-frequency, designed mutations were found to be in corroboration with the experimentally reported nirmatrelvir-resistant mutants and are naturally occurring. Our work on the targeted design of the nirmatrelvir-binding site offers a comprehensive picture of potential hotspot sites and resistance mutations in Mpro and is thus crucial in comprehending viral adaptation, robust antiviral design, and surveillance of evolving Mpro variations.
Collapse
|
69
|
Shu B, Ooi JSG, Tan AWK, Ng TS, Dejnirattisai W, Mongkolsapaya J, Fibriansah G, Shi J, Kostyuchenko VA, Screaton GR, Lok SM. CryoEM structures of the multimeric secreted NS1, a major factor for dengue hemorrhagic fever. Nat Commun 2022; 13:6756. [PMID: 36347841 PMCID: PMC9643530 DOI: 10.1038/s41467-022-34415-1] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2022] [Accepted: 10/25/2022] [Indexed: 11/11/2022] Open
Abstract
Dengue virus infection can cause dengue hemorrhagic fever (DHF). Dengue NS1 is multifunctional. The intracellular dimeric NS1 (iNS1) forms part of the viral replication complex. Previous studies suggest the extracellular secreted NS1 (sNS1), which is a major factor contributing to DHF, exists as hexamers. The structure of the iNS1 is well-characterised but not that of sNS1. Here we show by cryoEM that the recombinant sNS1 exists in multiple oligomeric states: the tetrameric (stable and loose conformation) and hexameric structures. Stability of the stable and loose tetramers is determined by the conformation of their N-terminal domain - elongated β-sheet or β-roll. Binding of an anti-NS1 Fab breaks the loose tetrameric and hexameric sNS1 into dimers, whereas the stable tetramer remains largely unbound. Our results show detailed quaternary organization of different oligomeric states of sNS1 and will contribute towards the design of dengue therapeutics.
Collapse
|
70
|
Wilson IM, Frazier MN, Li JL, Randall TA, Stanley RE. Biochemical Characterization of Emerging SARS-CoV-2 Nsp15 Endoribonuclease Variants. J Mol Biol 2022; 434:167796. [PMID: 35995266 PMCID: PMC9389836 DOI: 10.1016/j.jmb.2022.167796] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2022] [Revised: 08/08/2022] [Accepted: 08/15/2022] [Indexed: 12/02/2022]
Abstract
Global sequencing efforts from the ongoing COVID-19 pandemic, caused by the novel coronavirus SARS-CoV-2, continue to provide insight into the evolution of the viral genome. Coronaviruses encode 16 nonstructural proteins, within the first two-thirds of their genome, that facilitate viral replication and transcription as well as evasion of the host immune response. However, many of these viral proteins remain understudied. Nsp15 is a uridine-specific endoribonuclease conserved across all coronaviruses. The nuclease activity of Nsp15 helps the virus evade triggering an innate immune response. Understanding how Nsp15 has changed over the course of the pandemic, and how mutations affect its RNA processing function, will provide insight into the evolution of an oligomerization-dependent endoribonuclease and inform drug design. In combination with previous structural data, bioinformatics analyses of 1.9 + million SARS-CoV-2 sequences revealed mutations across Nsp15's three structured domains (N-terminal, Middle, EndoU). Selected Nsp15 variants were characterized biochemically and compared to wild type Nsp15. We found that mutations to important catalytic residues decreased cleavage activity but increased the hexamer/monomer ratio of the recombinant protein. Many of the highly prevalent variants we analyzed led to decreased nuclease activity as well as an increase in the inactive, monomeric form. Overall, our work establishes how Nsp15 variants seen in patient samples affect nuclease activity and oligomerization, providing insight into the effect of these variants in vivo.
Collapse
|
71
|
Pourfarjam Y, Ma Z, Kim IK. ATP enhances the error-prone ribonucleotide incorporation by the SARS-CoV-2 RNA polymerase. Biochem Biophys Res Commun 2022; 625:53-59. [PMID: 35947915 PMCID: PMC9344795 DOI: 10.1016/j.bbrc.2022.07.087] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2022] [Revised: 07/22/2022] [Accepted: 07/22/2022] [Indexed: 01/18/2023]
Abstract
The novel Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2 or COVID-19) has caused a global pandemic. The SARS-CoV-2 RNA genome is replicated by a conserved "core" replication-transcription complex (RTC) containing an error-prone RNA-dependent RNA polymerase holoenzyme (holo-RdRp, nsp12-nsp7-nsp8) and a RNA proofreading nuclease (nsp14-nsp10). Although structures and functions of SARS-CoV-2 holo-RdRp have been extensively studied and ribonucleotide-analog inhibitors, such as Remdesivir, have been treated for COVID-19 patients, the substrate and nucleotide specificity of SARS-CoV-2 holo-RdRp remain unknown. Here, our biochemical analysis of SARS-CoV-2 holo-RdRp reveals that it has a robust DNA-dependent RNA polymerase activity, in addition to its intrinsic RNA-dependent RNA polymerase activity. Strikingly, SARS-CoV-2 holo-RdRp fully extends RNAs with a low-fidelity even when only ATP and pyrimidine nucleotides, in particular CTP, are provided. This ATP-dependent error-prone ribonucleotide incorporation by SARS-CoV-2 holo-RdRp resists excision by the RNA proofreading nuclease in vitro. Our collective results suggest that a physiological concentration of ATP likely contributes to promoting the error-prone incorporation of ribonucleotides and ribonucleotide-analogs by SARS-CoV-2 holo-RdRp and provide a useful foundation to develop ribonucleotide analogs as an effective therapeutic strategy to combat coronavirus-mediated outbreak.
Collapse
|
72
|
Madushanka A, Verma N, Freindorf M, Kraka E. Papaya Leaf Extracts as Potential Dengue Treatment: An In-Silico Study. Int J Mol Sci 2022; 23:ijms232012310. [PMID: 36293162 PMCID: PMC9610845 DOI: 10.3390/ijms232012310] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2022] [Revised: 10/04/2022] [Accepted: 10/06/2022] [Indexed: 11/05/2022] Open
Abstract
Dengue fever (DF), dengue hemorrhagic fever (DHF), and dengue shock syndrome (DSS) cause serious public health problems, with nearly 390 million people affected and 20,000 deaths per year in tropical and subtropical countries. Despite numerous attempts, no antiviral drug or vaccine is currently available to combat the manifestation. The challenge of discovering an efficient vaccine is enhanced by the surplus presence of efficient vectors and drug resistance from the virus. For centuries, papaya (Carica papaya) extracts have been traditionally used to treat DF, DHF, and DSS. In the present study, we systematically investigated seven compounds isolated from papaya leaf extract with regard to their potential as inhibitors for non-structural (NS) proteins, NS3 and NS5, which play a crucial role in viral RNA replication. The computational tools applied stretched across classical molecular docking, molecular dynamics (MD) simulations and SwissADME used to calculate binding affinities; binding free energies; Absorption, Distribution, Metabolism, and Excretion (ADME); and drug-likeness properties, thus, identifying Kaempferol, Chlorogenic acid, and Quercetin as potential candidates, with Kaempferol and Quercetin scoring best. Therefore, for the Kaempferol and Quercetin complexes, hybrid quantum mechanical/molecular mechanical (QM/MM) geometry and frequency calculations were performed, followed by the local mode analysis developed in our group to quantify Kaempferol-NS and Quercetin-NS hydrogen bonding. Given the non-toxic nature and the wide availability of the Kaempferol and Quercetin papaya extract in almost all of the susceptible regions, and our results showing high NS3 and NS5 binding affinities and energies, strong hydrogen bonding with both NS3 and NS5, and excellent ADME properties, we suggest Kaempferol and Quercetin as a strong NS3 and NS5 inhibitor to be further investigated in vitro.
Collapse
|
73
|
Sajid Jamal QM, Alharbi AH, Ahmad V. Identification of doxorubicin as a potential therapeutic against SARS-CoV-2 (COVID-19) protease: a molecular docking and dynamics simulation studies. J Biomol Struct Dyn 2022; 40:7960-7974. [PMID: 33826483 PMCID: PMC8043163 DOI: 10.1080/07391102.2021.1905551] [Citation(s) in RCA: 18] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2020] [Accepted: 03/15/2021] [Indexed: 02/08/2023]
Abstract
After one year, the COVID-19 pandemic caused by SARS-CoV-2 is still the largest concern for the scientific community. Of the many recognized drug targets of SARS-CoV-2, the main protease is one of the most important target due to its function in viral replication. We conducted an in silico study with repurposing drugs of antibiotics class against virus protease and peptidase using AutoDock tool. The following significant binding energy interaction was observed with protease (PDB: 6LU7) like piperacillin -7.25; tobramycin -9.20 and doxorubicin (Doxo) -10.04 kcal/mol and with peptidase (PDB: 2GTB) piperacillin -7.08; tobramycin -8.54 and Doxo -9.89 kcal/mol. Furthermore, the interaction and stability behavior of the Doxo-protease and Doxo-peptidase complexes were analyzed for a 100-nanosecond (ns) time. Calculated RMSD values observed using molecular dynamics simulation (MDS) were found to be 0.15-0.25 nm, RMSF calculation per residues showed a value near 0.2 nm and Rg values remained approximately 2.25 nm. MM-PBSA analysis of total binding energy calculation of Doxo-protease and Doxo-peptidase complexes are found to be -148.692 and -105.367 kJ/mol, respectively. Moreover, amino acid residue ASP-197 showed the lowest contribution binding energy i.e. -18.1185 kJ/mol, and amino acid residue ASP-187 showed -17.0267 kJ/mol contribution energy. Thus, significant docking interaction and stable dynamicity of Doxo-protease complex with time was suggested that Doxo could be a choice to inhibit potentially the viral proteases that could prevent the entry inside the host cell to control the COVID-19 disease. Communicated by Ramaswamy H. Sarma.
Collapse
|
74
|
Frazier MN, Riccio AA, Wilson IM, Copeland WC, Stanley RE. Recent insights into the structure and function of coronavirus ribonucleases. FEBS Open Bio 2022; 12:1567-1583. [PMID: 35445579 PMCID: PMC9110870 DOI: 10.1002/2211-5463.13414] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2022] [Revised: 04/07/2022] [Accepted: 04/19/2022] [Indexed: 11/10/2022] Open
Abstract
Coronaviruses use approximately two-thirds of their 30-kb genomes to encode nonstructural proteins (nsps) with diverse functions that assist in viral replication and transcription, and evasion of the host immune response. The SARS-CoV-2 pandemic has led to renewed interest in the molecular mechanisms used by coronaviruses to infect cells and replicate. Among the 16 Nsps involved in replication and transcription, coronaviruses encode two ribonucleases that process the viral RNA-an exonuclease (Nsp14) and an endonuclease (Nsp15). In this review, we discuss recent structural and biochemical studies of these nucleases and the implications for drug discovery.
Collapse
|
75
|
Clemente CM, Freiberger MI, Ravetti S, Beltramo DM, Garro AG. An in silico analysis of Ibuprofen enantiomers in high concentrations of sodium chloride with SARS-CoV-2 main protease. J Biomol Struct Dyn 2022; 40:5653-5664. [PMID: 33459192 PMCID: PMC7832455 DOI: 10.1080/07391102.2021.1872420] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2020] [Accepted: 01/02/2021] [Indexed: 12/01/2022]
Abstract
2020 will be remembered worldwide for the outbreak of Coronavirus disease (COVID-19), which quickly spread until it was declared as a global pandemic. The main protease (Mpro) of SARS-CoV-2, a key enzyme in coronavirus, represents an attractive pharmacological target for inhibition of SARS-CoV-2 replication. Here, we evaluated whether the anti-inflammatory drug Ibuprofen, may act as a potential SARS-CoV-2 Mpro inhibitor, using an in silico study. From molecular dynamics (MD) simulations, we also evaluated the influence of ionic strength on the affinity and stability of the Ibuprofen-Mpro complexes. The docking analysis shows that R(-)Ibuprofen and S(+)Ibuprofen isomers can interact with multiple key residues of the main protease, through hydrophobic interactions and hydrogen bonds, with favourable binding energies (-6.2 and -5.7 kcal/mol, respectively). MM-GBSA and MM-PBSA calculations confirm the affinity of these complexes, in terms of binding energies. It also demonstrates that the ionic strength modifies significantly their binding affinities. Different structural parameters calculated from the MD simulations (120 ns) reveal that these complexes are conformational stable in the different conditions analysed. In this context, the results suggest that the condition 2 (0.25 NaCl) bind more tightly the Ibuprofen to Mpro than the others conditions. From the frustration analysis, we could characterize two important regions (Cys44-Pro52 and Linker loop) of this protein involved in the interaction with Ibuprofen. In conclusion, our findings allow us to propose that racemic mixtures of the Ibuprofen enantiomers might be a potential treatment option against SARS-CoV-2 Mpro. However, further research is necessary to determinate their possible medicinal use.Communicated by Ramaswamy H. Sarma.
Collapse
|