ABSTRACT
Ebola virus disease (EVD) causes outbreaks and epidemics in West Africa that persist until today. The envelope glycoprotein of Ebola virus (GP) consists of two subunits, GP1 and GP2, and plays a key role in anchoring or fusing the virus to the host cell in its active form on the virion surface. Toremifene (TOR) is a ligand that mainly acts as an estrogen receptor antagonist; however, a recent study showed a strong and efficient interaction with GP. In this context, we aimed to evaluate the energetic affinity features involved in the interaction between GP and toremifene by computer simulation techniques using the Molecular Fractionation Method with Conjugate Caps (MFCC) scheme and quantum-mechanical (QM) calculations, as well as missense mutations to assess protein stability. We identified ASP522, GLU100, TYR517, THR519, LEU186, LEU515 as the most attractive residues in the EBOV glycoprotein structure that form the binding pocket. We divided toremifene into three regions and evaluated that region i was more important than region iii and region ii for the formation of the TOR-GP1/GP2 complex, which might control the molecular remodeling process of TOR. The mutations that caused more destabilization were ARG134, LEU515, TYR517 and ARG559, while those that caused stabilization were GLU523 and ASP522. TYR517 is a critical residue for the binding of TOR, and is highly conserved among EBOV species. Our results may help to elucidate the mechanism of drug action on the GP protein of the Ebola virus and subsequently develop new pharmacological approaches against EVD.Communicated by Ramaswamy H. Sarma.
ABSTRACT
The emergent Zika virus (ZIKV) infection has become a threat to global health due to its association with severe neurological abnormalities, namely Guillain-Barré Syndrome (GBS) in adults and Congenital Zika virus Syndrome (CZS) in neonates. Many studies are nowadays being conducted to find an effective antiviral drug against ZIKV. In particular, NS2B-NS3 protease is an attractive drug target due to its essential function in viral replication, although a drug is not yet commercially available. In this context, we present here a comparative structural study, based on quantum chemistry calculations, to analyze the intermolecular binding energies between the crystallographic structure of NS2B-NS3 protease and dipeptide boronic acid (cn-716) and aldehyde (acyl-KR-aldehyde) peptidomimetic inhibitors, by using the molecular fractionation with conjugate caps (MFCC) scheme within the density functional theory (DFT) formalism. Most intermolecular interactions in cn-716/NS2B-NS3 (acyl-KR-aldehyde/NS2B-NS3) are due to the amino acid residues Asp83*, His51, Asp129, Ser81*, Gly133, Ala132, Tyr161, Asn152 and Asp75 (Asp83*, Asp129, His51, Asn152, Tyr161, Tyr130, Gly153, Gly151, Asp75, Pro131, and Gly82). Additionally, we have considered missense mutation analysis of these residues to evaluate the destabilization and the increase of the flexibility of the protease, showing that mutation of the residues Tyr161 and Tyr130 causes more impact. Our simulations are a valuable tool for a better understanding of the binding mechanism of recognized inhibitors of NS2B-NS3 protease, and can lead to the rational design and development of novel anti-Zika drugs with improved efficiency.
Subject(s)
Aldehydes/pharmacology , Antiviral Agents/pharmacology , Boronic Acids/pharmacology , Dipeptides/pharmacology , Protease Inhibitors/pharmacology , Zika Virus/drug effects , Aldehydes/chemistry , Antiviral Agents/chemistry , Boronic Acids/chemistry , Density Functional Theory , Dipeptides/chemistry , Microbial Sensitivity Tests , Molecular Conformation , Peptidomimetics/antagonists & inhibitors , Peptidomimetics/metabolism , Protease Inhibitors/chemistry , RNA Helicases/antagonists & inhibitors , RNA Helicases/metabolism , Serine Endopeptidases/metabolism , Viral Nonstructural Proteins/antagonists & inhibitors , Viral Nonstructural Proteins/metabolism , Zika Virus/metabolismABSTRACT
We intend to investigate the drug-binding energy of each nucleotide inside the aminoglycoside hygromycin B (hygB) binding site of 30S ribosomal RNA (rRNA) subunit by using the molecular fractionation with conjugate caps (MFCC) strategy based on the density functional theory (DFT), considering the functional LDA/PWC, OBS, and the dielectric constant parametrization. Aminoglycosides are bactericidal antibiotics that have high affinity to the prokaryotic rRNA, inhibiting the synthesis of proteins by acting on the main stages of the translation mechanism, whereas binding to rRNA 16S, a component of the 30S ribosomal subunit in prokaryotes. The identification of the nucleotides presenting the most negative binding energies allows us to stabilize hygB in a suitable binding pocket of the 30S ribosomal subunit. In addition, it should be highlighted that mutations in these residues may probably lead to resistance to ribosome-targeting antibiotics. Quantum calculations of aminoglycoside hygromycin B-ribosome complex might contribute to further quantum studies with antibiotics like macrolides and other aminoglycosides.
