1
|
Frezza E, Laage D, Duboué-Dijon E. Molecular Origin of Distinct Hydration Dynamics in Double Helical DNA and RNA Sequences. J Phys Chem Lett 2024; 15:4351-4358. [PMID: 38619551 DOI: 10.1021/acs.jpclett.4c00629] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/16/2024]
Abstract
Water molecules are essential to determine the structure of nucleic acids and mediate their interactions with other biomolecules. Here, we characterize the hydration dynamics of analogous DNA and RNA double helices with unprecedented resolution and elucidate the molecular origin of their differences: first, the localization of the slowest hydration water molecules─in the minor groove in DNA, next to phosphates in RNA─and second, the markedly distinct hydration dynamics of the two phosphate oxygen atoms OR and OS in RNA. Using our Extended Jump Model for water reorientation, we assess the relative importance of previously proposed factors, including the local topography, water bridges, and the presence of ions. We show that the slow hydration dynamics at RNA OR sites is not due to bridging water molecules but is caused by both the larger excluded volume and the stronger initial H-bond next to OR, due to the different phosphate orientations in A-form double helical RNA.
Collapse
Affiliation(s)
- Elisa Frezza
- Université Paris Cité, CNRS, CiTCoM, Paris 75006, France
| | - Damien Laage
- PASTEUR, Department of Chemistry, École Normale Supérieure-PSL, Sorbonne Université, CNRS, Paris 75005, France
| | - Elise Duboué-Dijon
- Université Paris Cité, CNRS, Laboratoire de Biochimie Théorique, 13 rue Pierre et Marie Curie, Paris 75005, France
| |
Collapse
|
2
|
Drici N. The influence of the hydrogen-bond network on the structure and dynamics of the RAPRKKG heptapeptide and its mutants. J Mol Graph Model 2023; 125:108598. [PMID: 37586130 DOI: 10.1016/j.jmgm.2023.108598] [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: 05/27/2023] [Revised: 08/02/2023] [Accepted: 08/08/2023] [Indexed: 08/18/2023]
Abstract
The structural behaviour of the RAPRKKG heptapeptide after individual or multiple mutations was inspected through molecular dynamics simulation. The nature of the mutations provided information on the flexibility of the heptapeptide and on how water molecules establish hydrogen bonds with it. The structural behaviour of the wild-type and the mutated structures were measured through the analysis of protein‒protein and protein‒solvent hydrogen bonds. The conformational behaviours of the different structures were analysed through free energy landscape analysis. The flexibility characteristics of the mutants seem to depend on the reorganization of water molecules and their static or dynamic behaviour around amino acid side chains.
Collapse
Affiliation(s)
- Nedjoua Drici
- University of Mostaganem, Abdelhamid Ibn Badis, Faculty of Exact Sciences and Informatics, Chemin des cretes ex INES, Mostaganem, 27000, Algeria; Laboratoire de Chimie Physique Macromoleculaire LCPM, University of Oran1 Ahmed benbella, Oran, 31000, Algeria.
| |
Collapse
|
3
|
Gabas F, Conte R, Ceotto M. Quantum Vibrational Spectroscopy of Explicitly Solvated Thymidine in Semiclassical Approximation. J Phys Chem Lett 2022; 13:1350-1355. [PMID: 35109652 PMCID: PMC8842300 DOI: 10.1021/acs.jpclett.1c04087] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2021] [Accepted: 01/28/2022] [Indexed: 06/14/2023]
Abstract
In this paper, we demonstrate the possibility to perform spectroscopy simulations of solvated biological species taking into consideration quantum effects and explicit solvation. We achieve this goal by interfacing our recently developed divide-and-conquer approach for semiclassical initial value representation molecular dynamics with the polarizable AMOEBABIO18 force field. The method is applied to the study of solvation of the thymidine nucleoside in two different polar solvents, water and N,N-dimethylformamide. Such systems are made of up to 2476 atoms. Experimental evidence concerning the different behavior of thymidine in the two solvents is well reproduced by our study, even though quantitative estimates are hampered by the limited accuracy of the classical force field employed. Overall, this study shows that semiclassically approximate quantum dynamical studies of explicitly solvated biological systems are both computationally affordable and insightful.
Collapse
|
4
|
Yuan Y, Fu S, Huo D, Su W, Zhang R, Wei J. Multipolar electrostatics for hairpin and pseudoknots in RNA: Improving the accuracy of force field potential energy function. J Comput Chem 2021; 42:771-786. [PMID: 33586809 DOI: 10.1002/jcc.26497] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2020] [Revised: 01/17/2021] [Accepted: 01/24/2021] [Indexed: 12/19/2022]
Abstract
Molecular dynamics (MD) simulations that rely on force field methods has been widely used to explore the structure and function of RNAs. However, the current commonly used force fields are limited by the electrostatic description offered by atomic charge, dipole and at most quadrupole moments, failing to capture the anisotropic picture of electronic features. Actually, the distribution of electrons around atomic nuclei is not spherically symmetric but is geometry dependent. A multipolar electrostatic model based on high rank multipole moments is described in this work, which allows us to combine polarizability and anisotropy of electron density. RNA secondary structure was taken as a research system, and its substructures including stem, loops (hairpin loop, bulge loop, internal loop, and multi-branch loop), and pseudoknots (H-type and K-type) were investigated, respectively, as well as the hairpin. First, the atom-atom electrostatic properties derived from one chain of a duplex RNA 2MVY in our previous work (Ref. 58) were measured by the pilot RNA systems of hairpin, hairpin loop, stem, and H-type pseudoknot, respectively. The prediction results were not satisfactory. Consequently, to obtain a general set of electrostatic parameters for RNA force fields, the convergence behavior of the atom-atom electrostatic interactions in the pilot RNA systems was explored using high rank atomic multipole moments. The pilot RNA systems were cut into four types of different-sized molecular fragments, and the single nucleotide fragment and nucleotide-paired fragment proved to be the most reasonable systems for base-unpairing regions and base-pairing regions to investigate the convergence behavior of all types of atom-atom electrostatic interactions, respectively. Transferability of the electrostatic properties drawn from the pilot RNA systems to the corresponding test systems was also investigated. Furthermore, the convergence behavior of atomic electrostatic interactions in other substructures including bulge loop, internal loop, multi-branch loop, and K-type pseudoknot was expected to be modeled via the hairpin.
Collapse
Affiliation(s)
- Yongna Yuan
- School of Information Science & Engineering, Lanzhou University, Lanzhou, Gansu, China
| | - Shaowei Fu
- School of Information Science & Engineering, Lanzhou University, Lanzhou, Gansu, China
| | - Dongxu Huo
- School of Information Science & Engineering, Lanzhou University, Lanzhou, Gansu, China
| | - Wei Su
- School of Information Science & Engineering, Lanzhou University, Lanzhou, Gansu, China
| | - Ruisheng Zhang
- School of Information Science & Engineering, Lanzhou University, Lanzhou, Gansu, China
| | - Jiaxuan Wei
- School of Information Science & Engineering, Lanzhou University, Lanzhou, Gansu, China
| |
Collapse
|
5
|
Kalra K, Gorle S, Cavallo L, Oliva R, Chawla M. Occurrence and stability of lone pair-π and OH-π interactions between water and nucleobases in functional RNAs. Nucleic Acids Res 2020; 48:5825-5838. [PMID: 32392301 PMCID: PMC7293021 DOI: 10.1093/nar/gkaa345] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2019] [Revised: 04/21/2020] [Accepted: 05/07/2020] [Indexed: 01/15/2023] Open
Abstract
We identified over 1000 instances of water-nucleobase stacking contacts in a variety of RNA molecules from a non-redundant set of crystal structures with resolution ≤3.0 Å. Such contacts may be of either the lone pair-π (lp-π) or the OH-π type, in nature. The distribution of the distances of the water oxygen from the nucleobase plane peaks at 3.5 Å for A, G and C, and approximately at 3.1-3.2 Å for U. Quantum mechanics (QM) calculations confirm, as expected, that the optimal energy is reached at a shorter distance for the lp-π interaction as compared to the OH-π one (3.0 versus 3.5 Å). The preference of each nucleobase for either type of interaction closely correlates with its electrostatic potential map. Furthermore, QM calculations show that for all the nucleobases a favorable interaction, of either the lp-π or the OH-π type, can be established at virtually any position of the water molecule above the nucleobase skeleton, which is consistent with the uniform projection of the OW atoms over the nucleobases ring we observed in the experimental occurrences. Finally, molecular dynamics simulations of a model system for the characterization of water-nucleobase stacking contacts confirm the stability of these interactions also under dynamic conditions.
Collapse
Affiliation(s)
- Kanav Kalra
- King Abdullah University of Science and Technology (KAUST), Physical Sciences and Engineering Division, Kaust Catalysis Center, Thuwal 23955-6900, Saudi Arabia
| | - Suresh Gorle
- Department of Biochemistry and Molecular Biology, University of Texas Medical Branch, Galveston, TX 77555, USA
| | - Luigi Cavallo
- King Abdullah University of Science and Technology (KAUST), Physical Sciences and Engineering Division, Kaust Catalysis Center, Thuwal 23955-6900, Saudi Arabia
| | - Romina Oliva
- Department of Sciences and Technologies, University Parthenope of Naples, Centro Direzionale Isola C4, I-80143 Naples, Italy
| | - Mohit Chawla
- King Abdullah University of Science and Technology (KAUST), Physical Sciences and Engineering Division, Kaust Catalysis Center, Thuwal 23955-6900, Saudi Arabia
| |
Collapse
|
6
|
Kumar N, Marx D. How do ribozymes accommodate additional water molecules upon hydrostatic compression deep into the kilobar pressure regime? Biophys Chem 2019; 252:106192. [DOI: 10.1016/j.bpc.2019.106192] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2019] [Revised: 05/23/2019] [Accepted: 05/23/2019] [Indexed: 12/19/2022]
|
7
|
Galindo-Murillo R, Cheatham TE. Lessons learned in atomistic simulation of double-stranded DNA: Solvation and salt concerns [Article v1.0]. ACTA ACUST UNITED AC 2019; 1. [PMID: 33073182 DOI: 10.33011/livecoms.1.2.9974] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Nucleic acids are highly charged macromolecules sensitive to their surroundings of water, salt, and other biomolecules. Molecular dynamics simulations with accurate biomolecular force fields provide a detailed atomistic view into DNA and RNA that has been useful to study the structure and dynamics of these molecules and their biological relevance. In this work we study the Drew-Dickerson dodecamer duplex with the sequence d(GCGCAATTGCGC)2 in three different salt concentrations and using different monvalent salt types to detect possible structural influence. Overall, the DNA shows no major structural changes regardless of amount or type of monovalent ions used. Our results show that only at very high salt conditions (5M) is a small structural effect observed in the DNA duplex, which mainly consist of narrowing of the grooves due to increased residence of ions. We also present the importance of sampling time to achieve a converged ensemble, which is of major relevance in any simulation to avoid biased or non-meaningful results.
Collapse
Affiliation(s)
- Rodrigo Galindo-Murillo
- Department of Medicinal Chemistry, L. S. Skaggs Pharmacy Institute, University of Utah, Salt Lake City, UT 84112
| | - Thomas E Cheatham
- Department of Medicinal Chemistry, L. S. Skaggs Pharmacy Institute, University of Utah, Salt Lake City, UT 84112
| |
Collapse
|
8
|
Wei W, Luo J, Waldispühl J, Moitessier N. Predicting Positions of Bridging Water Molecules in Nucleic Acid-Ligand Complexes. J Chem Inf Model 2019; 59:2941-2951. [PMID: 30998377 DOI: 10.1021/acs.jcim.9b00163] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
Over the past two decades, interests in DNA and RNA as drug targets have been growing rapidly. Following the trends observed with protein drug targets, computational approaches for drug design have been developed for this new class of molecules. Our efforts toward the development of a universal docking program, Fitted, led us to focus on nucleic acids. Throughout the development of this docking program, efforts were directed toward displaceable water molecules which must be accurately located for optimal docking-based drug discovery. However, although there is a plethora of methods to place water molecules in and around protein structures, there is, to the best of our knowledge, no such fully automated method for nucleic acids, which are significantly more polar and solvated than proteins. We report herein a new method, Splash'Em (Solvation Potential Laid around Statistical Hydration on Entire Macromolecules) developed to place water molecules within the binding cavity of nucleic acids. This fast method was shown to have high agreement with water positions in crystal structures and will therefore provide essential information to medicinal chemists.
Collapse
|
9
|
Deb I, Pal R, Sarzynska J, Lahiri A. Reparameterizations of theχTorsion and Lennard-JonesσParameters Improve the Conformational Characteristics of Modified Uridines. J Comput Chem 2016; 37:1576-88. [DOI: 10.1002/jcc.24374] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2016] [Accepted: 03/05/2016] [Indexed: 01/12/2023]
Affiliation(s)
- Indrajit Deb
- Department of Biophysics, Molecular Biology and Bioinformatics; University of Calcutta; 92 APC Road Kolkata West Bengal 700009 India
- Institute of Bioorganic Chemistry, Polish Academy of Sciences; Noskowskiego 12/14 Poznan 61-704 Poland
| | - Rupak Pal
- Department of Biophysics, Molecular Biology and Bioinformatics; University of Calcutta; 92 APC Road Kolkata West Bengal 700009 India
| | - Joanna Sarzynska
- Institute of Bioorganic Chemistry, Polish Academy of Sciences; Noskowskiego 12/14 Poznan 61-704 Poland
| | - Ansuman Lahiri
- Department of Biophysics, Molecular Biology and Bioinformatics; University of Calcutta; 92 APC Road Kolkata West Bengal 700009 India
| |
Collapse
|
10
|
Kührová P, Otyepka M, Šponer J, Banáš P. Are Waters around RNA More than Just a Solvent? - An Insight from Molecular Dynamics Simulations. J Chem Theory Comput 2015; 10:401-11. [PMID: 26579919 DOI: 10.1021/ct400663s] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Abstract
Hydrating water molecules are believed to be an inherent part of the RNA structure and have a considerable impact on RNA conformation. However, the magnitude and mechanism of the interplay between water molecules and the RNA structure are still poorly understood. In principle, such hydration effects can be studied by molecular dynamics (MD) simulations. In our recent MD studies, we observed that the choice of water model has a visible impact on the predicted structure and structural dynamics of RNA and, in particular, has a larger effect than type, parametrization, and concentration of the ions. Furthermore, the water model effect is sequence dependent and modulates the sequence dependence of A-RNA helical parameters. Clearly, the sensitivity of A-RNA structural dynamics to the water model parametrization is a rather spurious effect that complicates MD studies of RNA molecules. These results nevertheless suggest that the sequence dependence of the A-RNA structure, usually attributed to base stacking, might be driven by the structural dynamics of specific hydration. Here, we present a systematic MD study that aimed to (i) clarify the atomistic mechanism of the water model sensitivity and (ii) discover whether and to what extent specific hydration modulates the A-RNA structural variability. We carried out an extended set of MD simulations of canonical A-RNA duplexes with TIP3P, TIP4P/2005, TIP5P, and SPC/E explicit water models and found that different water models provided a different extent of water bridging between 2'-OH groups across the minor groove, which in turn influences their distance and consequently also inclination, roll, and slide parameters. Minor groove hydration is also responsible for the sequence dependence of these helical parameters. Our simulations suggest that TIP5P is not optimal for RNA simulations.
Collapse
Affiliation(s)
- Petra Kührová
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , tr. 17. Listopadu 12, 771 46, Olomouc, Czech Republic
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , tr. 17. Listopadu 12, 771 46, Olomouc, Czech Republic.,Institute of Biophysics, Academy of Sciences of the Czech Republic , Kralovopolska 135, 612 65 Brno, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic , Kralovopolska 135, 612 65 Brno, Czech Republic.,CEITEC - Central European Institute of Technology , Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , tr. 17. Listopadu 12, 771 46, Olomouc, Czech Republic.,Institute of Biophysics, Academy of Sciences of the Czech Republic , Kralovopolska 135, 612 65 Brno, Czech Republic
| |
Collapse
|
11
|
Aduri R, Psciuk BT, Saro P, Taniga H, Schlegel HB, SantaLucia J. AMBER Force Field Parameters for the Naturally Occurring Modified Nucleosides in RNA. J Chem Theory Comput 2015; 3:1464-75. [PMID: 26633217 DOI: 10.1021/ct600329w] [Citation(s) in RCA: 145] [Impact Index Per Article: 16.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Classical molecular dynamics (MD) simulations are useful for characterizing the structure and dynamics of biological macromolecules, ultimately, resulting in elucidation of biological function. The AMBER force field is widely used and has well-defined bond length, bond angle, partial charge, and van der Waals parameters for all the common amino acids and nucleotides, but it lacks parameters for many of the modifications found in nucleic acids and proteins. Presently there are 107 known naturally occurring modifications that play important roles in RNA stability, folding, and other functions. Modified nucleotides are found in almost all transfer RNAs, ribosomal RNAs of both the small and large subunits, and in many other functional RNAs. We developed force field parameters for the 107 modified nucleotides currently known to be present in RNA. The methodology used for deriving the modified nucleotide parameters is consistent with the methods used to develop the Cornell et al. force field. These parameters will improve the functionality of AMBER so that simulations can now be readily performed on diverse RNAs having post-transcriptional modifications.
Collapse
Affiliation(s)
- Raviprasad Aduri
- Department of Chemistry, Wayne State University, Detroit, Michigan 48202
| | - Brian T Psciuk
- Department of Chemistry, Wayne State University, Detroit, Michigan 48202
| | - Pirro Saro
- Department of Chemistry, Wayne State University, Detroit, Michigan 48202
| | - Hariprakash Taniga
- Department of Chemistry, Wayne State University, Detroit, Michigan 48202
| | | | - John SantaLucia
- Department of Chemistry, Wayne State University, Detroit, Michigan 48202
| |
Collapse
|
12
|
Basu I, Manna M, Mukhopadhyay C. Insights into the behavioral difference of water in the presence of GM1. FEBS Lett 2015; 589:3887-92. [DOI: 10.1016/j.febslet.2015.11.014] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2015] [Revised: 11/16/2015] [Accepted: 11/16/2015] [Indexed: 11/26/2022]
|
13
|
Deb I, Sarzynska J, Nilsson L, Lahiri A. Conformational preferences of modified uridines: comparison of AMBER derived force fields. J Chem Inf Model 2014; 54:1129-42. [PMID: 24697757 DOI: 10.1021/ci400582a] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The widespread occurrence of modified residues in RNA sequences necessitates development of accurate parameters for these modifications for reliable modeling of RNA structure and dynamics. A comprehensive set of parameters for the 107 naturally occurring RNA modifications was proposed by Aduri et al. (J. Chem. Theory Comput. 2007, 3, 1464-1475) for the AMBER FF99 force field. In this work, we tested these parameters on a set of modified uridine residues, namely, dihydrouridine, 2-thiouridine, 4-thiouridine, pseudouridine, and uridine-5-oxyacetic acid, by performing molecular dynamics and replica exchange molecular dynamics simulations of these nucleosides. Although our simulations using the FF99 force field did not, in general, reproduce the experimentally observed conformational characteristics well, combination of the parameter set with recent revisions of the FF99 force field for RNA showed noticeable improvement for some of the nucleosides.
Collapse
Affiliation(s)
- Indrajit Deb
- Department of Biophysics, Molecular Biology & Bioinformatics, University of Calcutta , Kolkata 700009, West Bengal, India
| | | | | | | |
Collapse
|
14
|
Ren P, Chun J, Thomas DG, Schnieders MJ, Marucho M, Zhang J, Baker NA. Biomolecular electrostatics and solvation: a computational perspective. Q Rev Biophys 2012; 45:427-91. [PMID: 23217364 PMCID: PMC3533255 DOI: 10.1017/s003358351200011x] [Citation(s) in RCA: 135] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
An understanding of molecular interactions is essential for insight into biological systems at the molecular scale. Among the various components of molecular interactions, electrostatics are of special importance because of their long-range nature and their influence on polar or charged molecules, including water, aqueous ions, proteins, nucleic acids, carbohydrates, and membrane lipids. In particular, robust models of electrostatic interactions are essential for understanding the solvation properties of biomolecules and the effects of solvation upon biomolecular folding, binding, enzyme catalysis, and dynamics. Electrostatics, therefore, are of central importance to understanding biomolecular structure and modeling interactions within and among biological molecules. This review discusses the solvation of biomolecules with a computational biophysics view toward describing the phenomenon. While our main focus lies on the computational aspect of the models, we provide an overview of the basic elements of biomolecular solvation (e.g. solvent structure, polarization, ion binding, and non-polar behavior) in order to provide a background to understand the different types of solvation models.
Collapse
Affiliation(s)
- Pengyu Ren
- Department of Biomedical Engineering, The University of Texas at Austin
| | | | | | | | - Marcelo Marucho
- Department of Physics and Astronomy, The University of Texas at San Antonio
| | - Jiajing Zhang
- Department of Biomedical Engineering, The University of Texas at Austin
| | - Nathan A. Baker
- To whom correspondence should be addressed. Pacific Northwest National Laboratory, PO Box 999, MSID K7-29, Richland, WA 99352. Phone: +1-509-375-3997,
| |
Collapse
|
15
|
Spasic A, Serafini J, Mathews DH. The Amber ff99 Force Field Predicts Relative Free Energy Changes for RNA Helix Formation. J Chem Theory Comput 2012; 8:2497-2505. [PMID: 23112748 PMCID: PMC3482406 DOI: 10.1021/ct300240k] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
The ability of the Amber ff99 force field to predict relative free energies of RNA helix formation was investigated. The test systems were three hexaloop RNA hairpins with identical loops and varying stems. The potential of mean force of stretching the hairpins from the native state to an extended conformation was calculated with umbrella sampling. Because the hairpins have identical loop sequence, the differences in free energy changes are only from the stem composition. The Amber ff99 force field was able to correctly predict the order of stabilities of the hairpins, although the magnitude of the free energy change is larger than that determined by optical melting experiments. The two measurements cannot be compared directly because the unfolded state in the optical melting experiments is a random coil, while the end state in the umbrella sampling simulations was an elongated chain. The calculations can be compared to reference data by using a thermodynamic cycle. By applying the thermodynamic cycle to the transitions between the hairpins using simulations and nearest neighbor data, agreement was found to be within the sampling error of simulations, thus demonstrating that ff99 force field is able to accurately predict relative free energies of RNA helix formation.
Collapse
Affiliation(s)
- Aleksandar Spasic
- Department of Biochemistry & Biophysics and Center for RNA Biology, University of Rochester Medical Center, Rochester, New York
| | - John Serafini
- Department of Biochemistry & Biophysics and Center for RNA Biology, University of Rochester Medical Center, Rochester, New York
| | - David H. Mathews
- Department of Biochemistry & Biophysics and Center for RNA Biology, University of Rochester Medical Center, Rochester, New York
- Department of Biostatistics & Computational Biology, University of Rochester Medical Center, Rochester, New York
| |
Collapse
|
16
|
Guga P, Koziołkiewicz M. Phosphorothioate nucleotides and oligonucleotides - recent progress in synthesis and application. Chem Biodivers 2012; 8:1642-81. [PMID: 21922655 DOI: 10.1002/cbdv.201100130] [Citation(s) in RCA: 62] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Affiliation(s)
- Piotr Guga
- Centre of Molecular and Macromolecular Studies, Polish Academy of Sciences, Department of Bioorganic Chemistry, Sienkiewicza 112, PL-90-363 Łódź.
| | | |
Collapse
|
17
|
Giel-Pietraszuk M, Fedoruk-Wyszomirska A, Barciszewski J. Effect of high hydrostatic pressure on hydration and activity of ribozymes. Mol Biol Rep 2010; 37:3713-9. [PMID: 20204525 DOI: 10.1007/s11033-010-0024-3] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2010] [Accepted: 02/17/2010] [Indexed: 10/19/2022]
Abstract
Formation and stabilization of RNA structure in the cell depends on its interaction with solvent and metal ions. High hydrostatic pressure (HHP) is a convenient tool in an analysis of the role of small molecules in the structure stabilization of biological macromolecules. Analysis of HHP effect and various concentrations of ions showed that water induce formation of the active ribozyme structure. So, it is clear that water is the driving force of conformational changes of nucleic acid.
Collapse
Affiliation(s)
- Małgorzata Giel-Pietraszuk
- Institute of Bioorganic Chemistry of the Polish Academy of Sciences, Noskowskiego 12/14, 61-704, Poznan, Poland
| | | | | |
Collapse
|
18
|
Michael LA, Chenault JA, Miller BR, Knolhoff AM, Nagan MC. Water, Shape Recognition, Salt Bridges, and Cation–Pi Interactions Differentiate Peptide Recognition of the HIV Rev-Responsive Element. J Mol Biol 2009; 392:774-86. [DOI: 10.1016/j.jmb.2009.07.047] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2009] [Revised: 07/06/2009] [Accepted: 07/16/2009] [Indexed: 10/20/2022]
|
19
|
A short guide for molecular dynamics simulations of RNA systems. Methods 2009; 47:187-97. [DOI: 10.1016/j.ymeth.2008.09.020] [Citation(s) in RCA: 64] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2008] [Accepted: 09/19/2008] [Indexed: 01/13/2023] Open
|
20
|
Dong F, Wagoner JA, Baker NA. Assessing the performance of implicit solvation models at a nucleic acid surface. Phys Chem Chem Phys 2008; 10:4889-902. [PMID: 18688533 PMCID: PMC2538626 DOI: 10.1039/b807384h] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
Implicit solvation models are popular alternatives to explicit solvent methods due to their ability to "pre-average" solvent behavior and thus reduce the need for computationally-expensive sampling. Previously, we have demonstrated that Poisson-Boltzmann models for polar solvation and integral-based models for nonpolar solvation can reproduce explicit solvation forces in a low-charge density protein system. In the present work, we examine the ability of these continuum models to describe solvation forces at the surface of a RNA hairpin. While these models do not completely describe all of the details of solvent behavior at this highly-charged biomolecular interface, they do provide a reasonable description of average solvation forces and therefore show significant promise for developing more robust implicit descriptions of solvent around nucleic acid systems for use in biomolecular simulation and modeling. Additionally, we observe fairly good transferability in the nonpolar model parameters optimized for protein systems, suggesting its robustness for modeling general nonpolar solvation phenomena in biomolecular systems.
Collapse
Affiliation(s)
- Feng Dong
- Merck & Co., Inc., 770 Sumneytown Pike, P.O. Box 4, WP42-330, West Point, PA 19486, USA. E-mail:
| | - Jason A. Wagoner
- Department of Chemistry, Stanford University, 333 Campus Drive #121, Mailbox 13, Stanford, CA 94305-5080, USA. E-mail:
| | - Nathan A. Baker
- To whom correspondence should be addressed. Department of Biochemistry and Molecular Biophysics, Center for Computational Biology, Washington University in St. Louis, 700. S. Euclid Ave., St. Louis, MO 63110, USA. E-mail:
| |
Collapse
|
21
|
Abstract
Enzymatic catalysis by RNA was discovered 25 years ago, yet mechanistic insights are emerging only slowly. Thought to be metalloenzymes at first, some ribozymes proved more versatile than anticipated when shown to utilize their own functional groups for catalysis. Recent evidence suggests that some may also judiciously place structural water molecules to shuttle protons in acid-base catalyzed reactions.
Collapse
Affiliation(s)
- Nils G Walter
- Department of Chemistry, Single Molecule Analysis Group, University of Michigan, 930 North University Avenue, Ann Arbor, MI 48019-1055, USA.
| |
Collapse
|
22
|
Auffinger P, Hashem Y. Nucleic acid solvation: from outside to insight. Curr Opin Struct Biol 2007; 17:325-33. [PMID: 17574833 DOI: 10.1016/j.sbi.2007.05.008] [Citation(s) in RCA: 112] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2007] [Revised: 03/28/2007] [Accepted: 05/31/2007] [Indexed: 11/18/2022]
Abstract
Nucleic acids are polyanionic molecules that were historically considered to be solely surrounded by a shell of water molecules and a neutralizing cloud of monovalent and divalent cations. In this respect, recent experimental and theoretical reports demonstrate that water molecules within complex nucleic acid structures can display very long residency times, and assist drug binding and catalytic reactions. Finally, anions can also bind to these polyanionic systems. Many of these recent insights are provided by state-of-the-art molecular dynamics simulations of nucleic acid systems, which will be described together with relevant methodological issues.
Collapse
Affiliation(s)
- Pascal Auffinger
- Architecture et réactivité de l'ARN, Université Louis Pasteur de Strasbourg, CNRS, IBMC, 15 rue René Descartes, 67084 Strasbourg, France.
| | | |
Collapse
|
23
|
Piton N, Mu Y, Stock G, Prisner TF, Schiemann O, Engels JW. Base-specific spin-labeling of RNA for structure determination. Nucleic Acids Res 2007; 35:3128-43. [PMID: 17452362 PMCID: PMC1891445 DOI: 10.1093/nar/gkm169] [Citation(s) in RCA: 123] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023] Open
Abstract
To facilitate the measurement of intramolecular distances in solvated RNA systems, a combination of spin-labeling, electron paramagnetic resonance (EPR), and molecular dynamics (MD) simulation is presented. The fairly rigid spin label 2,2,5,5-tetramethyl-pyrrolin-1-yloxyl-3-acetylene (TPA) was base and site specifically introduced into RNA through a Sonogashira palladium catalyzed cross-coupling on column. For this purpose 5-iodo-uridine, 5-iodo-cytidine and 2-iodo-adenosine phosphoramidites were synthesized and incorporated into RNA-sequences. Application of the recently developed ACE chemistry presented the main advantage to limit the reduction of the nitroxide to an amine during the oligonucleotide automated synthesis and thus to increase substantially the reliability of the synthesis and the yield of labeled oligonucleotides. 4-Pulse Electron Double Resonance (PELDOR) was then successfully used to measure the intramolecular spin-spin distances in six doubly labeled RNA-duplexes. Comparison of these results with our previous work on DNA showed that A- and B-Form can be differentiated. Using an all-atom force field with explicit solvent, MD simulations gave results in good agreement with the measured distances and indicated that the RNA A-Form was conserved despite a local destabilization effect of the nitroxide label. The applicability of the method to more complex biological systems is discussed.
Collapse
Affiliation(s)
- Nelly Piton
- Institute of Organic Chemistry and Chemical Biology, J. W. Goethe-University, Max-von-Laue Strasse 7, 60438 Frankfurt am Main, Germany, Institute of Physical and Theoretical Chemistry, J. W. Goethe-University, Max-von-Laue Strasse 7, 60438 Frankfurt am Main, Germany and Center of Biological Magnetic Resonance, Nanyang Technological University, 60 Nanyang Drive, Singapore 637551, Singapore
| | - Yuguang Mu
- Institute of Organic Chemistry and Chemical Biology, J. W. Goethe-University, Max-von-Laue Strasse 7, 60438 Frankfurt am Main, Germany, Institute of Physical and Theoretical Chemistry, J. W. Goethe-University, Max-von-Laue Strasse 7, 60438 Frankfurt am Main, Germany and Center of Biological Magnetic Resonance, Nanyang Technological University, 60 Nanyang Drive, Singapore 637551, Singapore
| | - Gerhard Stock
- Institute of Organic Chemistry and Chemical Biology, J. W. Goethe-University, Max-von-Laue Strasse 7, 60438 Frankfurt am Main, Germany, Institute of Physical and Theoretical Chemistry, J. W. Goethe-University, Max-von-Laue Strasse 7, 60438 Frankfurt am Main, Germany and Center of Biological Magnetic Resonance, Nanyang Technological University, 60 Nanyang Drive, Singapore 637551, Singapore
| | - Thomas F. Prisner
- Institute of Organic Chemistry and Chemical Biology, J. W. Goethe-University, Max-von-Laue Strasse 7, 60438 Frankfurt am Main, Germany, Institute of Physical and Theoretical Chemistry, J. W. Goethe-University, Max-von-Laue Strasse 7, 60438 Frankfurt am Main, Germany and Center of Biological Magnetic Resonance, Nanyang Technological University, 60 Nanyang Drive, Singapore 637551, Singapore
| | - Olav Schiemann
- Institute of Organic Chemistry and Chemical Biology, J. W. Goethe-University, Max-von-Laue Strasse 7, 60438 Frankfurt am Main, Germany, Institute of Physical and Theoretical Chemistry, J. W. Goethe-University, Max-von-Laue Strasse 7, 60438 Frankfurt am Main, Germany and Center of Biological Magnetic Resonance, Nanyang Technological University, 60 Nanyang Drive, Singapore 637551, Singapore
| | - Joachim W. Engels
- Institute of Organic Chemistry and Chemical Biology, J. W. Goethe-University, Max-von-Laue Strasse 7, 60438 Frankfurt am Main, Germany, Institute of Physical and Theoretical Chemistry, J. W. Goethe-University, Max-von-Laue Strasse 7, 60438 Frankfurt am Main, Germany and Center of Biological Magnetic Resonance, Nanyang Technological University, 60 Nanyang Drive, Singapore 637551, Singapore
- *To whom correspondence should be addressed.+49-69-798-29150+49-69-798-29148 Correspondence may also be addressed to Olav Schiemann. +49-69-798-29786+49-69-798-29404
| |
Collapse
|
24
|
Abstract
UNLABELLED SwS, based on a statistical analysis of crystallographic structures deposited in the NDB, is designed to provide an exhaustive overview of the solvation of nucleic acid structural elements through the generation of 3D solvent density maps. A first version (v1.0) of this web service focuses on the interaction of DNA, RNA and hybrid base pairs linked by two or three hydrogen bonds with water, cations and/or anions. Data provided by SwS are updated on a weekly basis and can be used by: (i) those involved in molecular dynamics simulation studies for validation purposes; (ii) crystallographers for help in the interpretation of solvent density maps; and all those involved in (iii) drug design and, more generally, in (iv) nucleic acid structural studies. SwS provides also statistical data related to the frequency of occurrence of different types of base pairs in crystallographic structures and the conformation of the involved nucleotides. This web service has been designed to allow a maximum of flexibility in terms of queries and has also been developed with didactic considerations in mind. AVAILABILITY http://www-ibmc.u-strasbg.fr/arn/sws.html
Collapse
Affiliation(s)
- Pascal Auffinger
- A.R.N., Université Louis Pasteur, IBMC-CNRS, Strasbourg, France.
| | | |
Collapse
|
25
|
Fiorucci S, Golebiowski J, Cabrol-Bass D, Antonczak S. Molecular simulations reveal a new entry site in quercetin 2,3-dioxygenase. A pathway for dioxygen? Proteins 2006; 64:845-50. [PMID: 16786599 DOI: 10.1002/prot.21042] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
Molecular dynamics simulations performed on quercetin 2,3-dioxygenase have shown the existence of a channel linking the bulk solvent and the cavity of the enzyme. Although much is known about the the oxygenolysis reaction catalyzed by this enzyme, the way dioxygen enters the active site has not been firmly established. The size, orientation and hydrophobic character of this channel suggests that it could provide an entrance for molecular dioxygen into the cavity. Free energy calculations show that such a process is likely to occur.
Collapse
Affiliation(s)
- Sébastien Fiorucci
- LCMBA, UMR-CNRS 6001, Faculté des Sciences, Université de Nice-Sophia Antipolis, Nice, France
| | | | | | | |
Collapse
|
26
|
Rhodes MM, Réblová K, Šponer J, Walter NG. Trapped water molecules are essential to structural dynamics and function of a ribozyme. Proc Natl Acad Sci U S A 2006; 103:13380-5. [PMID: 16938834 PMCID: PMC1569172 DOI: 10.1073/pnas.0605090103] [Citation(s) in RCA: 80] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Ribozymes are catalytically competent examples of highly structured noncoding RNAs, which are ubiquitous in the processing and regulation of genetic information. Combining explicit-solvent molecular dynamics simulation and single molecule fluorescence spectroscopy approaches, we find that a ribozyme from a subviral plant pathogen exhibits a coupled hydrogen bonding network that communicates dynamic structural rearrangements throughout the catalytic core in response to site-specific chemical modification. Trapped long-residency water molecules are critical for this network and only occasionally exchange with bulk solvent as they pass through a breathing interdomain base stack. These highly structured water molecules line up in a string that may potentially also be involved in specific base catalysis. Our observations suggest important, still underappreciated roles for specifically bound water molecules in the structural dynamics and function of noncoding RNAs.
Collapse
MESH Headings
- Base Pairing
- Base Sequence
- Catalysis
- Catalytic Domain
- Cations, Divalent/chemistry
- Computer Simulation
- Fluorescence Resonance Energy Transfer
- Hydrogen Bonding
- Kinetics
- Magnesium/chemistry
- Models, Molecular
- Models, Theoretical
- Mutation
- Nucleic Acid Conformation
- Protons
- RNA Precursors/chemistry
- RNA Precursors/genetics
- RNA Precursors/metabolism
- RNA, Catalytic/chemistry
- RNA, Catalytic/genetics
- RNA, Catalytic/metabolism
- RNA, Viral/chemistry
- RNA, Viral/genetics
- RNA, Viral/metabolism
- Solvents/chemistry
- Spectrometry, Fluorescence
- Static Electricity
- Water/chemistry
Collapse
Affiliation(s)
- Maria M. Rhodes
- *Department of Chemistry, Single Molecule Analysis Group, University of Michigan, 930 North University Avenue, Ann Arbor, MI 48109-1055
| | - Kamila Réblová
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic; and
| | - Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic; and
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, 166 10 Prague 6, Czech Republic
- To whom correspondence may be addressed. E-mail:
or
| | - Nils G. Walter
- *Department of Chemistry, Single Molecule Analysis Group, University of Michigan, 930 North University Avenue, Ann Arbor, MI 48109-1055
- To whom correspondence may be addressed. E-mail:
or
| |
Collapse
|
27
|
Rázga F, Zacharias M, Réblová K, Koca J, Sponer J. RNA kink-turns as molecular elbows: hydration, cation binding, and large-scale dynamics. Structure 2006; 14:825-35. [PMID: 16698544 DOI: 10.1016/j.str.2006.02.012] [Citation(s) in RCA: 43] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2005] [Revised: 02/02/2006] [Accepted: 02/21/2006] [Indexed: 10/24/2022]
Abstract
The presence of Kink-turns (Kt) at key functional sites in the ribosome (e.g., A-site finger and L7/L12 stalk) suggests that some Kink-turns can confer flexibility on RNA protuberances that regulate the traversal of tRNAs during translocation. Explicit solvent molecular dynamics demonstrates that Kink-turns can act as flexible molecular elbows. Kink-turns are associated with a unique network of long-residency static and dynamical hydration sites that is intimately involved in modulating their conformational dynamics. An implicit solvent conformational search confirms the flexibility of Kink-turns around their X-ray geometries and identifies a second low-energy region with open structures that could correspond to Kink-turn geometries seen in solution experiments. An extended simulation of Kt-42 with the factor binding site (helices 43 and 44) shows that the local Kt-42 elbow-like motion fully propagates beyond the Kink-turn, and that there is no other comparably flexible site in this rRNA region. Kink-turns could mediate large-scale adjustments of distant RNA segments.
Collapse
Affiliation(s)
- Filip Rázga
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 61265 Brno, Czech Republic
| | | | | | | | | |
Collapse
|
28
|
Krasovska MV, Sefcikova J, Réblová K, Schneider B, Walter NG, Sponer J. Cations and hydration in catalytic RNA: molecular dynamics of the hepatitis delta virus ribozyme. Biophys J 2006; 91:626-38. [PMID: 16617077 PMCID: PMC1483112 DOI: 10.1529/biophysj.105.079368] [Citation(s) in RCA: 106] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
The hepatitis delta virus (HDV) ribozyme is an RNA enzyme from the human pathogenic HDV. Cations play a crucial role in self-cleavage of the HDV ribozyme, by promoting both folding and chemistry. Experimental studies have revealed limited but intriguing details on the location and structural and catalytic functions of metal ions. Here, we analyze a total of approximately 200 ns of explicit-solvent molecular dynamics simulations to provide a complementary atomistic view of the binding of monovalent and divalent cations as well as water molecules to reaction precursor and product forms of the HDV ribozyme. Our simulations find that an Mg2+ cation binds stably, by both inner- and outer-sphere contacts, to the electronegative catalytic pocket of the reaction precursor, in a position to potentially support chemistry. In contrast, protonation of the catalytically involved C75 in the precursor or artificial placement of this Mg2+ into the product structure result in its swift expulsion from the active site. These findings are consistent with a concerted reaction mechanism in which C75 and hydrated Mg2+ act as general base and acid, respectively. Monovalent cations bind to the active site and elsewhere assisted by structurally bridging long-residency water molecules, but are generally delocalized.
Collapse
Affiliation(s)
- Maryna V Krasovska
- Institute of Biophysics, Academy of Sciences of the Czech Republic, 61265 Brno, Czech Republic
| | | | | | | | | | | |
Collapse
|
29
|
Abstract
Explicit solvent molecular dynamics (MD) simulations were carried out for sarcin-ricin domain (SRD) motifs from 23S (Escherichia coli) and 28S (rat) rRNAs. The SRD motif consists of GAGA tetraloop, G-bulged cross-strand A-stack, flexible region and duplex part. Detailed analysis of the overall dynamics, base pairing, hydration, cation binding and other SRD features is presented. The SRD is surprisingly static in multiple 25 ns long simulations and lacks any non-local motions, with root mean square deviation (r.m.s.d.) values between averaged MD and high-resolution X-ray structures of 1-1.4 A. Modest dynamics is observed in the tetraloop, namely, rotation of adenine in its apex and subtle reversible shift of the tetraloop with respect to the adjacent base pair. The deformed flexible region in low-resolution rat X-ray structure is repaired by simulations. The simulations reveal few backbone flips, which do not affect positions of bases and do not indicate a force field imbalance. Non-Watson-Crick base pairs are rigid and mediated by long-residency water molecules while there are several modest cation-binding sites around SRD. In summary, SRD is an unusually stiff rRNA building block. Its intrinsic structural and dynamical signatures seen in simulations are strikingly distinct from other rRNA motifs such as Loop E and Kink-turns.
Collapse
MESH Headings
- Animals
- Base Pairing
- Binding Sites
- Carbohydrates/chemistry
- Cations/chemistry
- Computer Simulation
- Crystallography, X-Ray
- Endoribonucleases/metabolism
- Escherichia coli/genetics
- Fungal Proteins/metabolism
- Hydrogen Bonding
- Models, Molecular
- Nucleic Acid Conformation
- RNA, Ribosomal, 23S/chemistry
- RNA, Ribosomal, 23S/metabolism
- RNA, Ribosomal, 28S/chemistry
- RNA, Ribosomal, 28S/metabolism
- Rats
- Ricin/metabolism
- Water/chemistry
Collapse
Affiliation(s)
- Nad'a Špačková
- Institute of Biophysics, Academy of Sciences of the Czech RepublicKrálovopolská 135, 612 65 Brno, Czech Republic
- To whom correspondence should be addressed. Tel: +420 541 517 109; Fax: +420 541 212 179;
| | - Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech RepublicKrálovopolská 135, 612 65 Brno, Czech Republic
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech RepublicFlemingovo náměstí 2, 166 10 Prague 6, Czech Republic
| |
Collapse
|
30
|
Bhide SY, Berkowitz ML. Structure and dynamics of water at the interface with phospholipid bilayers. J Chem Phys 2005; 123:224702. [PMID: 16375490 DOI: 10.1063/1.2132277] [Citation(s) in RCA: 98] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We have performed two molecular-dynamics simulations to study the structural and dynamical properties of water at the interface with phospholipid bilayers. In one of the simulations the bilayer contained neutral phospholipid molecules, dioleoylphosphatidylcholine (DOPC); in the second simulation the bilayer contained charged lipid molecules, dioleoylphosphatidylserine (DOPS). From the density profile of water we observe that water next to the DOPS bilayer is more perturbed as compared to water near the DOPC bilayer. Using an energetic criterion for the determination of hydrogen bonding we find that water molecules create strong hydrogen bonds with the headgroups of the phospholipid molecules. Due to the presence of these bonds and also due to the confinement of water, the translational and orientational dynamics of water at the interface are slowed down. The degree of slowing down of the dynamics depends upon the location of water molecules near a lipid headgroup.
Collapse
Affiliation(s)
- Shreyas Y Bhide
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA.
| | | |
Collapse
|
31
|
Koplin J, Mu Y, Richter C, Schwalbe H, Stock G. Structure and dynamics of an RNA tetraloop: a joint molecular dynamics and NMR study. Structure 2005; 13:1255-67. [PMID: 16154083 DOI: 10.1016/j.str.2005.05.015] [Citation(s) in RCA: 27] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2005] [Revised: 03/31/2005] [Accepted: 05/13/2005] [Indexed: 10/25/2022]
Abstract
Molecular dynamics simulations of the RNA tetraloop 5'-CGCUUUUGCG-3' with high melting temperature and significant conformational heterogeneity in explicit water solvent are presented and compared to NMR studies. The NMR data allow for a detailed test of the theoretical model, including the quality of the force field and the conformational sampling. Due to the conformational heterogeneity of the tetraloop, high temperature (350 K) and locally enhanced sampling simulations need to be invoked. The Amber98 force field leads to a good overall agreement with experimental data. Based on NMR data and a principal component analysis of the 350 K trajectory, the dynamic structure of the tetraloop is revealed. The principal component free energy surface exhibits four minima, which correspond to well-defined conformational structures that differ mainly by their base stacking in the loop region. No correlation between the motion of the sugar rings and the stacking dynamics of the loop bases is found.
Collapse
Affiliation(s)
- Jessica Koplin
- Institute of Physical and Theoretical Chemistry, Center for Biomolecular Magnetic Resonance, Johann Wolfgang Goethe-University, Frankfurt, Marie-Curie-Str 11, D-60439 Frankfurt/Main, Germany
| | | | | | | | | |
Collapse
|
32
|
Mu Y, Stock G. Conformational dynamics of RNA-peptide binding: a molecular dynamics simulation study. Biophys J 2005; 90:391-9. [PMID: 16239331 PMCID: PMC1367046 DOI: 10.1529/biophysj.105.069559] [Citation(s) in RCA: 31] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Molecular dynamics simulations of the binding of the heterochiral tripeptide KkN to the transactivation responsive (TAR) RNA of HIV-1 is presented, using an all-atom force field with explicit water. To obtain starting structures for the TAR-KkN complex, semirigid docking calculations were performed that employ an NMR structure of free TAR RNA. The molecular dynamics simulations show that the starting structures in which KkN binds to the major groove of TAR (as it is the case for the Tat-TAR complex of HIV-1) are unstable. On the other hand, the minor-groove starting structures are found to lead to several binding modes, which are stabilized by a complex interplay of stacking, hydrogen bonding, and electrostatic interactions. Although the ligand does not occupy the binding position of Tat protein, it is shown to hinder the interhelical motion of free TAR RNA. The latter is presumably necessary to achieve the conformational change of TAR RNA to bind Tat protein. Considering the time evolution of the trajectories, the binding process is found to be ligand-induced and cooperative. That is, the conformational rearrangement only occurs in the presence of the ligand and the concerted motion of the ligand and a large part of the RNA binding site is necessary to achieve the final low-energy binding state.
Collapse
Affiliation(s)
- Yuguang Mu
- School of Biological Sciences, Nanyang Technological University, Singapore and School of Physics and Microelectronics, Shandong University, Jinan, China
| | | |
Collapse
|
33
|
Cojocaru V, Klement R, Jovin TM. Loss of G-A base pairs is insufficient for achieving a large opening of U4 snRNA K-turn motif. Nucleic Acids Res 2005; 33:3435-46. [PMID: 15956103 PMCID: PMC1150281 DOI: 10.1093/nar/gki664] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Upon binding to the 15.5K protein, two tandem-sheared G–A base pairs are formed in the internal loop of the kink-turn motif of U4 snRNA (Kt-U4). We have reported that the folding of Kt-U4 is assisted by protein binding. Unstable interactions that contribute to a large opening of the free RNA (‘k–e motion’) were identified using locally enhanced sampling molecular dynamics simulations, results that agree with experiments. A detailed analysis of the simulations reveals that the k–e motion in Kt-U4 is triggered both by loss of G–A base pairs in the internal loop and backbone flexibility in the stems. Essential dynamics show that the loss of G–A base pairs is correlated along the first mode but anti-correlated along the third mode with the k–e motion. Moreover, when enhanced sampling was confined to the internal loop, the RNA adopted an alternative conformation characterized by a sharper kink, opening of G–A base pairs and modified stacking interactions. Thus, loss of G–A base pairs is insufficient for achieving a large opening of the free RNA. These findings, supported by previously published RNA structure probing experiments, suggest that G–A base pair formation occurs upon protein binding, thereby stabilizing a selective orientation of the stems.
Collapse
Affiliation(s)
| | | | - Thomas M. Jovin
- To whom correspondence should be addressed. Tel: +49 551 2011382; Fax: +49 551 2011467;
| |
Collapse
|
34
|
Rázga F, Koca J, Sponer J, Leontis NB. Hinge-like motions in RNA kink-turns: the role of the second a-minor motif and nominally unpaired bases. Biophys J 2005; 88:3466-85. [PMID: 15722438 PMCID: PMC1305493 DOI: 10.1529/biophysj.104.054916] [Citation(s) in RCA: 74] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Kink-turn (K-turn) motifs are asymmetric internal loops found at conserved positions in diverse RNAs, with sharp bends in phosphodiester backbones producing V-shaped structures. Explicit-solvent molecular dynamics simulations were carried out for three K-turns from 23S rRNA, i.e., Kt-38 located at the base of the A-site finger, Kt-42 located at the base of the L7/L12 stalk, and Kt-58 located in domain III, and for the K-turn of human U4 snRNA. The simulations reveal hinge-like K-turn motions on the nanosecond timescale. The first conserved A-minor interaction between the K-turn stems is entirely stable in all simulations. The angle between the helical arms of Kt-38 and Kt-42 is regulated by local variations of the second A-minor (type I) interaction between the stems. Its variability ranges from closed geometries to open ones stabilized by insertion of long-residency waters between adenine and cytosine. The simulated A-minor geometries fully agree with x-ray data. Kt-58 and Kt-U4 exhibit similar elbow-like motions caused by conformational change of the adenosine from the nominally unpaired region. Despite the observed substantial dynamics of K-turns, key tertiary interactions are stable and no sign of unfolding is seen. We suggest that some K-turns are flexible elements mediating large-scale ribosomal motions during the protein synthesis cycle.
Collapse
Affiliation(s)
- Filip Rázga
- National Centre for Biomolecular Research, Brno, Czech Republic
| | | | | | | |
Collapse
|
35
|
Chapter 6 Molecular Modeling and Atomistic Simulation of Nucleic Acids. ACTA ACUST UNITED AC 2005. [DOI: 10.1016/s1574-1400(05)01006-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/05/2023]
|
36
|
Auffinger P, Bielecki L, Westhof E. Anion binding to nucleic acids. Structure 2004; 12:379-88. [PMID: 15016354 DOI: 10.1016/j.str.2004.02.015] [Citation(s) in RCA: 102] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2003] [Revised: 12/01/2003] [Accepted: 12/07/2003] [Indexed: 11/20/2022]
Abstract
Nucleic acids are generally considered as efficient cation binders. Therefore, the likelihood that negatively charged ions might intrude their first hydration shell is rarely considered. Here, we show on the basis of (i) a survey of the Nucleic Acid Database, (ii) several structures extracted from the Cambridge Structural Database, and (iii) molecular dynamics simulations, that the nucleotide electropositive edges involving mainly amino, imino, and hydroxyl groups can cast specific anion binding sites. These binding sites constitute also good locations for the binding of the negatively charged groups of the Asp and Glu residues or the nucleic acid phosphate groups. Furthermore, it is observed in several instances that anions, like water molecules and cations, do mediate protein/nucleic acid interactions. Thus, anions as well as negatively charged groups are directly involved in specific recognition and folding phenomena involving polyanionic nucleic acids.
Collapse
Affiliation(s)
- Pascal Auffinger
- Institut de Biologie Moléculaire et Cellulaire du CNRS, Modélisations et Simulations des Acides Nucléiques, UPR 9002, 15, rue René Descartes, 67084 Strasbourg Cedex, France.
| | | | | |
Collapse
|
37
|
Golebiowski J, Antonczak S, Fernandez-Carmona J, Condom R, Cabrol-Bass D. Closing loop base pairs in RNA loop-loop complexes: structural behavior, interaction energy and solvation analysis through molecular dynamics simulations. J Mol Model 2004; 10:408-17. [PMID: 15597210 DOI: 10.1007/s00894-004-0216-7] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2004] [Accepted: 09/10/2004] [Indexed: 10/26/2022]
Abstract
Nanosecond molecular dynamics using the Ewald summation method have been performed to elucidate the structural and energetic role of the closing base pair in loop-loop RNA duplexes neutralized by Mg2+ counterions in aqueous phases. Mismatches GA, CU and Watson-Crick GC base pairs have been considered for closing the loop of an RNA in complementary interaction with HIV-1 TAR. The simulations reveal that the mismatch GA base, mediated by a water molecule, leads to a complex that presents the best compromise between flexibility and energetic contributions. The mismatch CU base pair, in spite of the presence of an inserted water molecule, is too short to achieve a tight interaction at the closing-loop junction and seems to force TAR to reorganize upon binding. An energetic analysis has allowed us to quantify the strength of the interactions of the closing and the loop-loop pairs throughout the simulations. Although the water-mediated GA closing base pair presents an interaction energy similar to that found on fully geometry-optimized structure, the water-mediated CU closing base pair energy interaction reaches less than half the optimal value.
Collapse
Affiliation(s)
- Jérôme Golebiowski
- Laboratoire Arômes, Synthèses, Interactions, Faculté des sciences, Université de Nice-Sophia Antipolis, Nice Cedex 2, 06108, France.
| | | | | | | | | |
Collapse
|
38
|
Réblová K, Spacková N, Sponer JE, Koca J, Sponer J. Molecular dynamics simulations of RNA kissing-loop motifs reveal structural dynamics and formation of cation-binding pockets. Nucleic Acids Res 2004; 31:6942-52. [PMID: 14627827 PMCID: PMC290250 DOI: 10.1093/nar/gkg880] [Citation(s) in RCA: 56] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Explicit solvent molecular dynamics (MD) simulations were carried out for three RNA kissing-loop complexes. The theoretical structure of two base pairs (2 bp) complex of H3 stem-loop of Moloney murine leukemia virus agrees with the NMR structure with modest violations of few NMR restraints comparable to violations present in the NMR structure. In contrast to the NMR structure, however, MD shows relaxed intermolecular G-C base pairs. The core region of the kissing complex forms a cation-binding pocket with highly negative electrostatic potential. The pocket shows nanosecond-scale breathing motions coupled with oscillations of the whole molecule. Additional simulations were carried out for 6 bp kissing complexes of the DIS HIV-1 subtypes A and B. The simulated structures agree well with the X-ray data. The subtype B forms a novel four-base stack of bulged-out adenines. Both 6 bp kissing complexes have extended cation-binding pockets in their central parts. While the pocket of subtype A interacts with two hexacoordinated Mg2+ ions and one sodium ion, pocket of subtype B is filled with a string of three delocalized Na+ ions with residency times of individual cations 1-2 ns. The 6 bp complexes show breathing motions of the cation-binding pockets and loop major grooves.
Collapse
Affiliation(s)
- Kamila Réblová
- National Center for Biomolecular Research, Masaryk University, Kotlárská 2, 611 37 Brno, Czech Republic
| | | | | | | | | |
Collapse
|
39
|
Golebiowski J, Antonczak S, Di-Giorgio A, Condom R, Cabrol-Bass D. Molecular dynamics simulation of hepatitis C virus IRES IIId domain: structural behavior, electrostatic and energetic analysis. J Mol Model 2003; 10:60-8. [PMID: 14691674 DOI: 10.1007/s00894-003-0170-9] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2003] [Accepted: 10/31/2003] [Indexed: 11/26/2022]
Abstract
The dynamic behavior of the HCV IRES IIId domain is analyzed by means of a 2.6-ns molecular dynamics simulation, starting from an NMR structure. The simulation is carried out in explicit water with Na+ counterions, and particle-mesh Ewald summation is used for the electrostatic interactions. In this work, we analyze selected patterns of the helix that are crucial for IRES activity and that could be considered as targets for the intervention of inhibitors, such as the hexanucleotide terminal loop (more particularly its three consecutive guanines) and the loop-E motif. The simulation has allowed us to analyze the dynamics of the loop substructure and has revealed a behavior among the guanine bases that might explain the different role of the third guanine of the GGG triplet upon molecular recognition. The accessibility of the loop-E motif and the loop major and minor groove is also examined, as well as the effect of Na+ or Mg2+ counterion within the simulation. The electrostatic analysis reveals several ion pockets, not discussed in the experimental structure. The positions of these ions are useful for locating specific electrostatic recognition sites for potential inhibitor binding.
Collapse
Affiliation(s)
- Jérôme Golebiowski
- Laboratoire Arômes Synthèses Interactions, Faculté des Sciences de Nice Sophia-Antipolis, 06108 cedex 2, Nice, France.
| | | | | | | | | |
Collapse
|
40
|
Vicens Q, Westhof E. Molecular recognition of aminoglycoside antibiotics by ribosomal RNA and resistance enzymes: an analysis of x-ray crystal structures. Biopolymers 2003; 70:42-57. [PMID: 12925992 DOI: 10.1002/bip.10414] [Citation(s) in RCA: 123] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
The potential of RNA molecules to be used as therapeutic targets by small inhibitors is now well established. In this fascinating wide-open field, aminoglycoside antibiotics constitute the most studied family of RNA binding drugs. Within the last three years, several x-ray crystal structures were solved for aminoglycosides complexed to one of their main natural targets in the bacterial cell, the decoding aminoacyl-tRNA site (A site). Other crystallographic structures have revealed the binding modes of aminoglycosides to the three existing types of resistance-associated enzymes. The present review summarizes the various aspects of the molecular recognition of aminoglycosides by these natural RNA or protein receptors. The analysis and the comparisons of the detailed interactions offer insights that are helpful in designing new generations of antibiotics.
Collapse
Affiliation(s)
- Quentin Vicens
- Institut de Biologie Moléculaire et Cellulaire du CNRS, Modélisation et simulations des Acides Nucléiques, UPR 9002, Université Louis Pasteur, 15 rue René Descartes, 67084 Strasbourg Cedex, France
| | | |
Collapse
|
41
|
Réblová K, Spacková N, Stefl R, Csaszar K, Koca J, Leontis NB, Sponer J. Non-Watson-Crick basepairing and hydration in RNA motifs: molecular dynamics of 5S rRNA loop E. Biophys J 2003; 84:3564-82. [PMID: 12770867 PMCID: PMC1302943 DOI: 10.1016/s0006-3495(03)75089-9] [Citation(s) in RCA: 93] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022] Open
Abstract
Explicit solvent and counterion molecular dynamics simulations have been carried out for a total of >80 ns on the bacterial and spinach chloroplast 5S rRNA Loop E motifs. The Loop E sequences form unique duplex architectures composed of seven consecutive non-Watson-Crick basepairs. The starting structure of spinach chloroplast Loop E was modeled using isostericity principles, and the simulations refined the geometries of the three non-Watson-Crick basepairs that differ from the consensus bacterial sequence. The deep groove of Loop E motifs provides unique sites for cation binding. Binding of Mg(2+) rigidifies Loop E and stabilizes its major groove at an intermediate width. In the absence of Mg(2+), the Loop E motifs show an unprecedented degree of inner-shell binding of monovalent cations that, in contrast to Mg(2+), penetrate into the most negative regions inside the deep groove. The spinach chloroplast Loop E shows a marked tendency to compress its deep groove compared with the bacterial consensus. Structures with a narrow deep groove essentially collapse around a string of Na(+) cations with long coordination times. The Loop E non-Watson-Crick basepairing is complemented by highly specific hydration sites ranging from water bridges to hydration pockets hosting 2 to 3 long-residing waters. The ordered hydration is intimately connected with RNA local conformational variations.
Collapse
Affiliation(s)
- Kamila Réblová
- National Center for Biomolecular Research, Brno, Czech Republic
| | | | | | | | | | | | | |
Collapse
|
42
|
Ennifar E, Walter P, Dumas P. A crystallographic study of the binding of 13 metal ions to two related RNA duplexes. Nucleic Acids Res 2003; 31:2671-82. [PMID: 12736317 PMCID: PMC156032 DOI: 10.1093/nar/gkg350] [Citation(s) in RCA: 122] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023] Open
Abstract
Metal ions, and magnesium in particular, are known to be involved in RNA folding by stabilizing secondary and tertiary structures, and, as cofactors, in RNA enzymatic activity. We have conducted a systematic crystallographic analysis of cation binding to the duplex form of the HIV-1 RNA dimerization initiation site for the subtype-A and -B natural sequences. Eleven ions (K+, Pb2+, Mn2+, Ba2+, Ca2+, Cd2+, Sr2+, Zn2+, Co2+, Au3+ and Pt4+) and two hexammines [Co (NH3)6]3+ and [Ru (NH3)6]3+ were found to bind to the DIS duplex structure. Although the two sequences are very similar, strong differences were found in their cation binding properties. Divalent cations bind almost exclusively, as Mg2+, at 'Hoogsteen' sites of guanine residues, with a cation-dependent affinity for each site. Notably, a given cation can have very different affinities for a priori equivalent sites within the same molecule. Surprisingly, none of the two hexammines used were able to efficiently replace hexahydrated magnesium. Instead, [Co (NH3)4]3+ was seen bound by inner-sphere coordination to the RNA. This raises some questions about the practical use of [Co (NH3)6]3+ as a [Mg (H2O)6]2+ mimetic. Also very unexpected was the binding of the small Au3+ cation exactly between the Watson-Crick sites of a G-C base pair after an obligatory deprotonation of N1 of the guanine base. This extensive study of metal ion binding using X-ray crystallography significantly enriches our knowledge on the binding of middleweight or heavy metal ions to RNA, particularly compared with magnesium.
Collapse
Affiliation(s)
- Eric Ennifar
- Institut de Biologie Moléculaire et Cellulaire, CNRS-UPR 9002, 15 rue René Descartes, 67084 Strasbourg cedex, France
| | | | | |
Collapse
|