1
|
Fereydooni Naghani F, Emamian S, Zare K. A comprehensive theoretical analysis on the intermolecular hydrogen bond interactions with the Lewis bases having multiple hydrogen bonding ability. J PHYS ORG CHEM 2022. [DOI: 10.1002/poc.4357] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
| | - Saeedreza Emamian
- Chemistry Department, Shahrood Branch Islamic Azad University Shahrood Iran
| | - Karim Zare
- Department of Chemistry, Science and Research Branch Islamic Azad University Tehran Iran
| |
Collapse
|
2
|
Gervasoni S, Spencer J, Hinchliffe P, Pedretti A, Vairoletti F, Mahler G, Mulholland AJ. A multiscale approach to predict the binding mode of metallo beta-lactamase inhibitors. Proteins 2022; 90:372-384. [PMID: 34455628 PMCID: PMC8944931 DOI: 10.1002/prot.26227] [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] [Received: 03/23/2021] [Revised: 06/09/2021] [Accepted: 08/18/2021] [Indexed: 02/03/2023]
Abstract
Antibiotic resistance is a major threat to global public health. β-lactamases, which catalyze breakdown of β-lactam antibiotics, are a principal cause. Metallo β-lactamases (MBLs) represent a particular challenge because they hydrolyze almost all β-lactams and to date no MBL inhibitor has been approved for clinical use. Molecular simulations can aid drug discovery, for example, predicting inhibitor complexes, but empirical molecular mechanics (MM) methods often perform poorly for metalloproteins. Here we present a multiscale approach to model thiol inhibitor binding to IMP-1, a clinically important MBL containing two catalytic zinc ions, and predict the binding mode of a 2-mercaptomethyl thiazolidine (MMTZ) inhibitor. Inhibitors were first docked into the IMP-1 active site, testing different docking programs and scoring functions on multiple crystal structures. Complexes were then subjected to molecular dynamics (MD) simulations and subsequently refined through QM/MM optimization with a density functional theory (DFT) method, B3LYP/6-31G(d), increasing the accuracy of the method with successive steps. This workflow was tested on two IMP-1:MMTZ complexes, for which it reproduced crystallographically observed binding, and applied to predict the binding mode of a third MMTZ inhibitor for which a complex structure was crystallographically intractable. We also tested a 12-6-4 nonbonded interaction model in MD simulations and optimization with a SCC-DFTB QM/MM approach. The results show the limitations of empirical models for treating these systems and indicate the need for higher level calculations, for example, DFT/MM, for reliable structural predictions. This study demonstrates a reliable computational pipeline that can be applied to inhibitor design for MBLs and other zinc-metalloenzyme systems.
Collapse
Affiliation(s)
- Silvia Gervasoni
- Department of Pharmaceutical Sciences, University of Milan, Milan, Italy
| | - James Spencer
- School of Cellular and Molecular Medicine, University of Bristol, Bristol, UK
| | - Philip Hinchliffe
- School of Cellular and Molecular Medicine, University of Bristol, Bristol, UK
| | | | - Franco Vairoletti
- Laboratorio de Química Farmacéutica, Departamento de Química Orgánica, Facultad de Química, Universidad de la República (UdelaR), Avda. General Flores 2124, Montevideo, Uruguay
| | - Graciela Mahler
- Laboratorio de Química Farmacéutica, Departamento de Química Orgánica, Facultad de Química, Universidad de la República (UdelaR), Avda. General Flores 2124, Montevideo, Uruguay
| | | |
Collapse
|
3
|
Non-covalent interactions abetted supramolecular arrangements of N-Substituted benzylidene acetohydrazide to direct its solid-state network. J Mol Struct 2021. [DOI: 10.1016/j.molstruc.2020.129827] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
|
4
|
Khalid M, Ali A, Abid S, Tahir MN, Khan MU, Ashfaq M, Imran M, Ahmad A. Facile Ultrasound‐Based Synthesis, SC‐XRD, DFT Exploration of the Substituted Acyl‐Hydrazones: An Experimental and Theoretical Slant towards Supramolecular Chemistry. ChemistrySelect 2020. [DOI: 10.1002/slct.202003589] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Muhammad Khalid
- Department of Chemistry Khwaja Fareed University of Engineering & Information Technology Rahim Yar Khan 64200 Pakistan
| | - Akbar Ali
- Department of Chemistry University of Malakand Chakdara 18800, Lower Dir Khyber Pakhtunkhwa Pakistan
- Department of Chemistry University of Sargodha, 40100 Pakistan
| | - Saba Abid
- Department of Chemistry Khwaja Fareed University of Engineering & Information Technology Rahim Yar Khan 64200 Pakistan
| | | | - Muhammad Usman Khan
- Department of Applied Chemistry Government College University Faisalabad 38000 Pakistan
| | - Muhammad Ashfaq
- Department of Physics University of Sargodha Sargodha Pakistan
| | - Muhammad Imran
- Department of Chemistry Faculty of Science King Khalid University, P.O. Box 9004 Abha 61413 Saudi Arabia
| | - Anees Ahmad
- Instituto de Química de Sao Carlos Universidade de Sao Paulo CEP 13560-970 Sao Carlos SP Brazil
| |
Collapse
|
5
|
Mráziková K, Mlýnský V, Kührová P, Pokorná P, Kruse H, Krepl M, Otyepka M, Banáš P, Šponer J. UUCG RNA Tetraloop as a Formidable Force-Field Challenge for MD Simulations. J Chem Theory Comput 2020; 16:7601-7617. [PMID: 33215915 DOI: 10.1021/acs.jctc.0c00801] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
Explicit solvent atomistic molecular dynamics (MD) simulations represent an established technique to study structural dynamics of RNA molecules and an important complement for diverse experimental methods. However, performance of molecular mechanical (MM) force fields (ff's) remains far from satisfactory even after decades of development, as apparent from a problematic structural description of some important RNA motifs. Actually, some of the smallest RNA molecules belong to the most challenging systems for MD simulations and, among them, the UUCG tetraloop is saliently difficult. We report a detailed analysis of UUCG MD simulations, depicting the sequence of events leading to the loss of the UUCG native state during MD simulations. The total amount of MD simulation data analyzed in this work is close to 1.3 ms. We identify molecular interactions, backbone conformations, and substates that are involved in the process. Then, we unravel specific ff deficiencies using diverse quantum mechanical/molecular mechanical (QM/MM) and QM calculations. Comparison between the MM and QM methods shows discrepancies in the description of the 5'-flanking phosphate moiety and both signature sugar-base interactions. Our work indicates that poor behavior of the UUCG tetraloop in simulations is a complex issue that cannot be attributed to one dominant and straightforwardly correctable factor. Instead, there is a concerted effect of multiple ff inaccuracies that are coupled and amplifying each other. We attempted to improve the simulation behavior by some carefully tailored interventions, but the results were still far from satisfactory, underlying the difficulties in development of accurate nucleic acid ff's.
Collapse
Affiliation(s)
- Klaudia Mráziková
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic.,National Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic
| | - Vojtěch Mlýnský
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Petra Kührová
- Regional Centre of Advanced Technologies and Materials, Faculty of Science, Palacký University, Šlechtitelů 27, 783 71 Olomouc, Czech Republic
| | - Pavlína Pokorná
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic.,National Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic
| | - Holger Kruse
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Miroslav Krepl
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Faculty of Science, Palacký University, Šlechtitelů 27, 783 71 Olomouc, Czech Republic
| | - Michal Otyepka
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Faculty of Science, Palacký University, Šlechtitelů 27, 783 71 Olomouc, Czech Republic
| | - Pavel Banáš
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Faculty of Science, Palacký University, Šlechtitelů 27, 783 71 Olomouc, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| |
Collapse
|
6
|
Khalid M, Ali A, Tariq J, Tahir MN, Aliabad HAR, Hussain I, Ashfaq M, Khan MU. Stabilization of Supramolecular Assembly of N‐Substituted Benzylidene Acetohydrazide Analogs by Non‐Covalent Interactions: A Concise Experimental and Theoretical Approach. ChemistrySelect 2020. [DOI: 10.1002/slct.202002653] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Affiliation(s)
- Muhammad Khalid
- Department of Chemistry Khwaja Fareed University of Engineering & Information Technology Rahim Yar Khan 64200 Pakistan
| | - Akbar Ali
- Department of Chemistry University of Malakand Chakdara 18800 Lower Dir Khyber Pakhtunkhwa Pakistan
- Department of Chemistry University of Sargodha 40100 Pakistan
| | - Jahrukh Tariq
- Department of Chemistry Khwaja Fareed University of Engineering & Information Technology Rahim Yar Khan 64200 Pakistan
| | | | | | - Ishtiaq Hussain
- Department of Pharmaceutical Sciences Abbottabad University of Science and Technology KPK Pakistan
| | - Muhammad Ashfaq
- Department of Physics University of Sargodha Sargodha Pakistan
| | | |
Collapse
|
7
|
Green Synthesis, SC-XRD, Non-Covalent Interactive Potential and Electronic Communication via DFT Exploration of Pyridine-Based Hydrazone. CRYSTALS 2020. [DOI: 10.3390/cryst10090778] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
Ultrasound-based synthesis at room temperature produces valuable compounds greener and safer than most other methods. This study presents the sonochemical fabrication and characterization of a pyridine-based halogenated hydrazone, (E)-2-((6-chloropyridin-2-yl)oxy)-N′-(2-hydroxybenzylidene) acetohydrazide (HBPAH). The NMR spectroscopic technique was used to determine the structure, while SC-XRD confirmed its crystalline nature. Our structural studies revealed that strong, inter-molecular attractive forces stabilize this crystalline organic compound. Moreover, the compound was optimized at the B3LYP/6-311G(d,p) level using the Crystallographic Information File (CIF). Natural bonding orbital (NBO) and natural population analysis (NPA) were performed at the same level using optimized geometry. Time-dependent density functional theory (DFT) was performed at the B3LYP/6-311G (d,p) method to calculate the frontier molecular orbitals (FMOs) and molecular electrostatic potential (MEP). The global reactivity descriptors were determined using HOMO and LUMO energy gaps. Theoretical calculations based on the Quantum Theory of Atoms in Molecules (QT-AIM) and Hirshfeld analyses identified the non-covalent and covalent interactions of the HBPAH compound. Consequently, QT-AIM and Hirshfeld analyses agree with experimental results.
Collapse
|
8
|
Emamian S, Lu T, Kruse H, Emamian H. Exploring Nature and Predicting Strength of Hydrogen Bonds: A Correlation Analysis Between Atoms-in-Molecules Descriptors, Binding Energies, and Energy Components of Symmetry-Adapted Perturbation Theory. J Comput Chem 2019; 40:2868-2881. [PMID: 31518004 DOI: 10.1002/jcc.26068] [Citation(s) in RCA: 425] [Impact Index Per Article: 85.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2019] [Revised: 08/06/2019] [Accepted: 08/22/2019] [Indexed: 01/08/2023]
Abstract
This work studies the underlying nature of H-bonds (HBs) of different types and strengths and tries to predict binding energies (BEs) based on the properties derived from wave function analysis. A total of 42 HB complexes constructed from 28 neutral and 14 charged monomers were considered. This set was designed to sample a wide range of HB strengths to obtain a complete view about HBs. BEs were derived with the accurate coupled cluster singles and doubles with perturbative triples correction (CCSD(T))(T) method and the physical components of the BE were investigated by symmetry-adapted perturbation theory (SAPT). Quantum theory of atoms-in-molecules (QTAIM) descriptors and other HB indices were calculated based on high-quality density functional theory wave functions. We propose a new and rigorous classification of H-bonds (HBs) based on the SAPT decomposition. Neutral complexes are either classified as "very weak" HBs with a BE ≥ -2.5 kcal/mol that are mainly dominated by both dispersion and electrostatic interactions or as "weak-to-medium" HBs with a BE varying between -2.5 and -14.0 kcal/mol that are only dominated by electrostatic interactions. On the other hand, charged complexes are divided into "medium" HBs with a BE in the range of -11.0 to -15.0 kcal/mol, which are mainly dominated by electrostatic interactions, or into "strong" HBs whose BE is more negative than -15.0 kcal/mol, which are mainly dominated by electrostatic together with induction interactions. Among various explored correlations between BEs and wave function-based HB descriptors, a fairly satisfactory correlation was found for the electron density at the bond critical point (BCP; ρBCP ) of HBs. The fitted equation for neutral complexes is BE/kcal/mol = - 223.08 × ρBCP /a. u. + 0.7423 with a mean absolute percentage error (MAPE) of 14.7%, while that for charged complexes is BE/kcal/mol = - 332.34 × ρBCP /a. u. - 1.0661 with a MAPE of 10.0%. In practice, these equations may be used for a quick estimation of HB BEs, for example, for intramolecular HBs or large HB networks in biomolecules. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Saeedreza Emamian
- Chemistry Department, Shahrood Branch, Islamic Azad University, Shahrood, Iran
| | - Tian Lu
- Beijing Kein Research Center for Natural Sciences, Beijing, 100022, People's Republic of China
| | - Holger Kruse
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65, Brno, Czech Republic
| | - Hamidreza Emamian
- Materials Science & Metallurgy Department, Technology & Processing of New Materials Research Center, Neyshabur Branch, Islamic Azad University, Neyshabur, Iran
| |
Collapse
|
9
|
Górecka KM, Krepl M, Szlachcic A, Poznański J, Šponer J, Nowotny M. RuvC uses dynamic probing of the Holliday junction to achieve sequence specificity and efficient resolution. Nat Commun 2019; 10:4102. [PMID: 31506434 PMCID: PMC6736871 DOI: 10.1038/s41467-019-11900-8] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2019] [Accepted: 08/09/2019] [Indexed: 12/31/2022] Open
Abstract
Holliday junctions (HJs) are four-way DNA structures that occur in DNA repair by homologous recombination. Specialized nucleases, termed resolvases, remove (i.e., resolve) HJs. The bacterial protein RuvC is a canonical resolvase that introduces two symmetric cuts into the HJ. For complete resolution of the HJ, the two cuts need to be tightly coordinated. They are also specific for cognate DNA sequences. Using a combination of structural biology, biochemistry, and a computational approach, here we show that correct positioning of the substrate for cleavage requires conformational changes within the bound DNA. These changes involve rare high-energy states with protein-assisted base flipping that are readily accessible for the cognate DNA sequence but not for non-cognate sequences. These conformational changes and the relief of protein-induced structural tension of the DNA facilitate coordination between the two cuts. The unique DNA cleavage mechanism of RuvC demonstrates the importance of high-energy conformational states in nucleic acid readouts.
Collapse
Affiliation(s)
- Karolina Maria Górecka
- Laboratory of Protein Structure, International Institute of Molecular and Cell Biology, 4 Trojdena St., 02-109, Warsaw, Poland
| | - Miroslav Krepl
- Institute of Biophysics of the Czech Academy of Sciences, Kralovopolska 135, 612 65, Brno, Czech Republic.
| | - Aleksandra Szlachcic
- Laboratory of Protein Structure, International Institute of Molecular and Cell Biology, 4 Trojdena St., 02-109, Warsaw, Poland
| | - Jarosław Poznański
- Institute of Biochemistry and Biophysics Polish Academy of Sciences, 5a Pawinskiego St., 02-106, Warsaw, Poland
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences, Kralovopolska 135, 612 65, Brno, Czech Republic
- Regional Centre of Advanced Technologies and Materials, Faculty of Science, Palacky University Olomouc, Slechtitelu 27, 771 46, Olomouc, Czech Republic
| | - Marcin Nowotny
- Laboratory of Protein Structure, International Institute of Molecular and Cell Biology, 4 Trojdena St., 02-109, Warsaw, Poland.
| |
Collapse
|
10
|
Pokorná P, Krepl M, Bártová E, Šponer J. Role of Fine Structural Dynamics in Recognition of Histone H3 by HP1γ(CSD) Dimer and Ability of Force Fields to Describe Their Interaction Network. J Chem Theory Comput 2019; 15:5659-5673. [DOI: 10.1021/acs.jctc.9b00434] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Pavlína Pokorná
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
- National Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic
| | - Miroslav Krepl
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Eva Bártová
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| |
Collapse
|
11
|
Litman J, Thiel AC, Schnieders MJ. Scalable Indirect Free Energy Method Applied to Divalent Cation-Metalloprotein Binding. J Chem Theory Comput 2019; 15:4602-4614. [DOI: 10.1021/acs.jctc.9b00147] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Jacob Litman
- University of Iowa, Department of Biochemistry, 51 Newton Road, 4-403 Bowen Science Building, Iowa City, Iowa 52242, United States
| | - Andrew C. Thiel
- University of Iowa, Department of Biomedical Engineering, 103 South Capitol Street, 5601 Seamans Center for the Engineering Arts and Sciences, Iowa City, Iowa 52242, United States
| | - Michael J. Schnieders
- University of Iowa, Department of Biochemistry, 51 Newton Road, 4-403 Bowen Science Building, Iowa City, Iowa 52242, United States
- University of Iowa, Department of Biomedical Engineering, 103 South Capitol Street, 5601 Seamans Center for the Engineering Arts and Sciences, Iowa City, Iowa 52242, United States
| |
Collapse
|
12
|
Pokorná P, Kruse H, Krepl M, Šponer J. QM/MM Calculations on Protein-RNA Complexes: Understanding Limitations of Classical MD Simulations and Search for Reliable Cost-Effective QM Methods. J Chem Theory Comput 2018; 14:5419-5433. [PMID: 30199638 DOI: 10.1021/acs.jctc.8b00670] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Although atomistic explicit-solvent Molecular Dynamics (MD) is a popular tool to study protein-RNA recognition, satisfactory MD description of protein-RNA complexes is not always achieved. Unfortunately, it is often difficult to separate MD simulation instabilities primarily caused by the simple point-charge molecular mechanics (MM) force fields from problems related to the notorious uncertainties in the starting structures. Herein, we report a series of large-scale QM/MM calculations on the U1A protein-RNA complex. This experimentally well-characterized system has an intricate protein-RNA interface, which is very unstable in MD simulations. The QM/MM calculations identify several H-bonds poorly described by the MM method and thus indicate the sources of instabilities of the U1A interface in MD simulations. The results suggest that advanced QM/MM computations could be used to indirectly rationalize problems seen in MM-based MD simulations of protein-RNA complexes. As the most accurate QM method, we employ the computationally demanding meta-GGA density functional TPSS-D3(BJ)/def2-TZVP level of theory. Because considerably faster methods would be needed to extend sampling and to study even larger protein-RNA interfaces, a set of low-cost QM/MM methods is compared to the TPSS-D3(BJ)/def2-TZVP data. The PBEh-3c and B97-3c density functional composite methods appear to be suitable for protein-RNA interfaces. In contrast, HF-3c and the tight-binding Hamiltonians DFTB3-D3 and GFN-xTB perform unsatisfactorily and do not provide any advantage over the MM description. These conclusions are supported also by similar analysis of a simple HutP protein-RNA interface, which is well-described by MD with the exception of just one H-bond. Some other methodological aspects of QM/MM calculations on protein-RNA interfaces are discussed.
Collapse
Affiliation(s)
- Pavlína Pokorná
- Institute of Biophysics of the Czech Academy of Sciences , Královopolská 135 , Brno 612 65 , Czech Republic
| | - Holger Kruse
- Institute of Biophysics of the Czech Academy of Sciences , Královopolská 135 , Brno 612 65 , Czech Republic
| | - Miroslav Krepl
- Institute of Biophysics of the Czech Academy of Sciences , Královopolská 135 , Brno 612 65 , Czech Republic
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences , Královopolská 135 , Brno 612 65 , Czech Republic.,CEITEC - Central European Institute of Technology, Masaryk University , Campus Bohunice, Kamenice 5 , Brno 625 00 , Czech Republic
| |
Collapse
|
13
|
Krepl M, Vögele J, Kruse H, Duchardt-Ferner E, Wöhnert J, Sponer J. An intricate balance of hydrogen bonding, ion atmosphere and dynamics facilitates a seamless uracil to cytosine substitution in the U-turn of the neomycin-sensing riboswitch. Nucleic Acids Res 2018; 46:6528-6543. [PMID: 29893898 PMCID: PMC6061696 DOI: 10.1093/nar/gky490] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2018] [Accepted: 05/23/2018] [Indexed: 01/04/2023] Open
Abstract
The neomycin sensing riboswitch is the smallest biologically functional RNA riboswitch, forming a hairpin capped with a U-turn loop-a well-known RNA motif containing a conserved uracil. It was shown previously that a U→C substitution of the eponymous conserved uracil does not alter the riboswitch structure due to C protonation at N3. Furthermore, cytosine is evolutionary permitted to replace uracil in other U-turns. Here, we use molecular dynamics simulations to study the molecular basis of this substitution in the neomycin sensing riboswitch and show that a structure-stabilizing monovalent cation-binding site in the wild-type RNA is the main reason for its negligible structural effect. We then use NMR spectroscopy to confirm the existence of this cation-binding site and to demonstrate its effects on RNA stability. Lastly, using quantum chemical calculations, we show that the cation-binding site is altering the electronic environment of the wild-type U-turn so that it is more similar to the cytosine mutant. The study reveals an amazingly complex and delicate interplay between various energy contributions shaping up the 3D structure and evolution of nucleic acids.
Collapse
Affiliation(s)
- Miroslav Krepl
- Institute of Biophysics of the Czech Academy of Sciences, Kralovopolska 135, 612 65 Brno, Czech Republic
| | - Jennifer Vögele
- Institute of Molecular Biosciences and Center for Biomolecular Magnetic Resonance (BMRZ), Goethe-University Frankfurt, Max-von-Laue-Str. 9, 60438 Frankfurt, Germany
| | - Holger Kruse
- Institute of Biophysics of the Czech Academy of Sciences, Kralovopolska 135, 612 65 Brno, Czech Republic
| | - Elke Duchardt-Ferner
- Institute of Molecular Biosciences and Center for Biomolecular Magnetic Resonance (BMRZ), Goethe-University Frankfurt, Max-von-Laue-Str. 9, 60438 Frankfurt, Germany
| | - Jens Wöhnert
- Institute of Molecular Biosciences and Center for Biomolecular Magnetic Resonance (BMRZ), Goethe-University Frankfurt, Max-von-Laue-Str. 9, 60438 Frankfurt, Germany
| | - Jiri Sponer
- Institute of Biophysics of the Czech Academy of Sciences, Kralovopolska 135, 612 65 Brno, Czech Republic
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, 17. listopadu 12, 771 46 Olomouc, Czech Republic
| |
Collapse
|
14
|
Šponer J, Bussi G, Krepl M, Banáš P, Bottaro S, Cunha RA, Gil-Ley A, Pinamonti G, Poblete S, Jurečka P, Walter NG, Otyepka M. RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview. Chem Rev 2018; 118:4177-4338. [PMID: 29297679 PMCID: PMC5920944 DOI: 10.1021/acs.chemrev.7b00427] [Citation(s) in RCA: 336] [Impact Index Per Article: 56.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2017] [Indexed: 12/14/2022]
Abstract
With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.
Collapse
Affiliation(s)
- Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences , Kralovopolska 135 , Brno 612 65 , Czech Republic
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Miroslav Krepl
- Institute of Biophysics of the Czech Academy of Sciences , Kralovopolska 135 , Brno 612 65 , Czech Republic
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| | - Sandro Bottaro
- Structural Biology and NMR Laboratory, Department of Biology , University of Copenhagen , Copenhagen 2200 , Denmark
| | - Richard A Cunha
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Alejandro Gil-Ley
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Giovanni Pinamonti
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Simón Poblete
- Scuola Internazionale Superiore di Studi Avanzati , Via Bonomea 265 , Trieste 34136 , Italy
| | - Petr Jurečka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| | - Nils G Walter
- Single Molecule Analysis Group and Center for RNA Biomedicine, Department of Chemistry , University of Michigan , Ann Arbor , Michigan 48109 , United States
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science , Palacky University Olomouc , 17. listopadu 12 , Olomouc 771 46 , Czech Republic
| |
Collapse
|