1
|
Boychuk BTA, Wetmore SD. Assessment of Density Functional Theory Methods for the Structural Prediction of Transition and Post-Transition Metal-Nucleic Acid Complexes. J Chem Theory Comput 2023. [PMID: 37399186 DOI: 10.1021/acs.jctc.3c00127] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/05/2023]
Abstract
Understanding the structure of metal-nucleic acid systems is important for many applications such as the design of new pharmaceuticals, metal detection platforms, and nanomaterials. Herein, we explore the ability of 20 density functional theory (DFT) functionals to reproduce the crystal structure geometry of transition and post-transition metal-nucleic acid complexes identified in the Protein Data Bank and Cambridge Structural Database. The environmental extremes of the gas phase and implicit water were considered, and analysis focused on the global and inner coordination geometry, including the coordination distances. Although gas-phase calculations were unable to describe the structure of 12 out of the 53 complexes in our test set regardless of the DFT functional considered, accounting for the broader environment through implicit solvation or constraining the model truncation points to crystallographic coordinates generally afforded agreement with the experimental structure, suggesting that functional performance for these systems is likely due to the models rather than the methods. For the remaining 41 complexes, our results show that the reliability of functionals depends on the metal identity, with the magnitude of error varying across the periodic table. Furthermore, minimal changes in the geometries of these metal-nucleic acid complexes occur upon use of the Stuttgart-Dresden effective core potential and/or inclusion of an implicit water environment. The overall top three performing functionals are ωB97X-V, ωB97X-D3(BJ), and MN15, which reliably describe the structure of a broad range of metal-nucleic acid systems. Other suitable functionals include MN15-L, which is a cheaper alternative to MN15, and PBEh-3c, which is commonly used in QM/MM calculations of biomolecules. In fact, these five methods were the only functionals tested to reproduce the coordination sphere of Cu2+-containing complexes. For metal-nucleic acid systems that do not contain Cu2+, ωB97X and ωB97X-D are also suitable choices. These top-performing methods can be utilized in future investigations of diverse metal-nucleic acid complexes of relevance to biology and material science.
Collapse
Affiliation(s)
- Briana T A Boychuk
- Department of Chemistry and Biochemistry, University of Lethbridge, 4401 University Drive West, Lethbridge, AB T1K 3M4, Canada
| | - Stacey D Wetmore
- Department of Chemistry and Biochemistry, University of Lethbridge, 4401 University Drive West, Lethbridge, AB T1K 3M4, Canada
| |
Collapse
|
2
|
Mráziková K, Kruse H, Mlýnský V, Auffinger P, Šponer J. Multiscale Modeling of Phosphate···π Contacts in RNA U-Turns Exposes Differences between Quantum-Chemical and AMBER Force Field Descriptions. J Chem Inf Model 2022; 62:6182-6200. [PMID: 36454943 DOI: 10.1021/acs.jcim.2c01064] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/03/2022]
Abstract
Phosphate···π, also called anion···π, contacts occur between nucleobases and anionic phosphate oxygens (OP2) in r(GNRA) and r(UNNN) U-turn motifs (N = A,G,C,U; R = A,G). These contacts were investigated using state-of-the-art quantum-chemical methods (QM) to characterize their physicochemical properties and to serve as a reference to evaluate AMBER force field (AFF) performance. We found that phosphate···π interaction energies calculated with the AFF for dimethyl phosphate···nucleobase model systems are less stabilizing in comparison with double-hybrid DFT and that minimum contact distances are larger for all nucleobases. These distance stretches are also observed in large-scale AFF vs QM/MM computations and classical molecular dynamics (MD) simulations on several r(gcGNRAgc) tetraloop hairpins when compared to experimental data extracted from X-ray/cryo-EM structures (res. ≤ 2.5 Å) using the WebFR3D bioinformatic tool. MD simulations further revealed shifted OP2/nucleobase positions. We propose that discrepancies between the QM and AFF result from a combination of missing polarization in the AFF combined with too large AFF Lennard-Jones (LJ) radii of nucleobase carbon atoms in addition to an exaggerated short-range repulsion of the r-12 LJ repulsive term. We compared these results with earlier data gathered on lone pair···π contacts in CpG Z-steps occurring in r(UNCG) tetraloops. In both instances, charge transfer calculations do not support any significant n → π* donation effects. We also investigated thiophosphate···π contacts that showed reduced stabilizing interaction energies when compared to phosphate···π contacts. Thus, we challenge suggestions that the experimentally observed enhanced thermodynamic stability of phosphorothioated r(GNRA) tetraloops can be explained by larger London dispersion.
Collapse
Affiliation(s)
- Klaudia Mráziková
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65Brno, Czech Republic.,National Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kamenice 5, 625 00Brno, Czech Republic
| | - Holger Kruse
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65Brno, Czech Republic
| | - Vojtěch Mlýnský
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65Brno, Czech Republic
| | - Pascal Auffinger
- Architecture et Réactivité de l'ARN, Université de Strasbourg, Institut de Biologie Moléculaire et Cellulaire du CNRS, Strasbourg67084, France
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65Brno, Czech Republic
| |
Collapse
|
3
|
Balasubramanian S, Senapati S. Dynamics and Barrier of Movements of Sodium and Potassium Ions Across the Oxytricha nova G-Quadruplex Core. J Phys Chem B 2020; 124:11055-11066. [PMID: 33238706 DOI: 10.1021/acs.jpcb.0c04826] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
G-quadruplexes (GQs) are highly stable noncanonical forms of nucleic acids that are present in important genomic regions. The central core of the GQ is lined up by four closely spaced carbonyl groups from the G-quartets, and the resulting electrostatic repulsion is neutralized by the coordinating cations. In spite of several reports on GQ structure and cation-GQ interactions, the atomic- to molecular-level understanding of the ion dynamics and ion exchange in the GQ core is quite poor. Here, we attempt to elucidate the mechanism of Na+ and K+ binding to the GQ core and trace the exchange of these ions with the ions in bulk by means of all-atomic molecular dynamics (MD) simulations. One of the most studied GQs, Oxytricha nova telomeric G-quadruplex (OxyGQ), is taken as the representative GQ. Subsequently, umbrella sampling MD simulations were performed to elucidate the energetics of ion translocation from one end to the other end of the GQ central core. Our study highlights the importance of ion hydration for the uptake and correct positioning of the cations in the core. The free-energy landscape of ion transport has shown favorable in-plane binding of Na+ ions with GQ quartets, which matches very well with the crystal structure. The binding of K+ ions, on the other hand, was out-of-plane and its translocation required a larger barrier to cross.
Collapse
Affiliation(s)
- Sangeetha Balasubramanian
- Department of Biotechnology, BJM School of Biosciences, Indian Institute of Technology Madras, Chennai 600 036, India
| | - Sanjib Senapati
- Department of Biotechnology, BJM School of Biosciences, Indian Institute of Technology Madras, Chennai 600 036, India
| |
Collapse
|
4
|
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
|
5
|
Ahmed A, Angell R, Oxenford S, Worthington J, Williams N, Barton N, Fowler TG, O’Flynn DE, Sunose M, McConville M, Vo T, Wilson WD, Karim SA, Morton JP, Neidle S. Asymmetrically Substituted Quadruplex-Binding Naphthalene Diimide Showing Potent Activity in Pancreatic Cancer Models. ACS Med Chem Lett 2020; 11:1634-1644. [PMID: 32832034 PMCID: PMC7429975 DOI: 10.1021/acsmedchemlett.0c00317] [Citation(s) in RCA: 28] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2020] [Accepted: 07/16/2020] [Indexed: 12/13/2022] Open
Abstract
Targeting of genomic quadruplexes is an approach to treating complex human cancers. We describe a series of tetra-substituted naphthalene diimide (ND) derivatives with a phenyl substituent directly attached to the ND core. The lead compound (SOP1812) has 10 times superior cellular and in vivo activity compared with previous ND compounds and nanomolar binding to human quadruplexes. The pharmacological properties of SOP1812 indicate good bioavailability, which is consistent with the in vivo activity in xenograft and genetic models for pancreatic cancer. Transcriptome analysis shows that it down-regulates several cancer gene pathways, including Wnt/β-catenin signaling.
Collapse
Affiliation(s)
- Ahmed
A. Ahmed
- School
of Pharmacy, University College London, London WC1N 1AX, United Kingdom
| | - Richard Angell
- School
of Pharmacy, University College London, London WC1N 1AX, United Kingdom
| | - Sally Oxenford
- School
of Pharmacy, University College London, London WC1N 1AX, United Kingdom
| | - Jenny Worthington
- Axis
Bio Discovery Services, Ltd., Coleraine, Northern Ireland BT51 3RP, United Kingdom
| | - Nicole Williams
- Axis
Bio Discovery Services, Ltd., Coleraine, Northern Ireland BT51 3RP, United Kingdom
| | - Naomi Barton
- Sygnature
Discovery Limited, Nottingham NG1 1GR, United Kingdom
| | | | | | - Mihiro Sunose
- Sygnature
Discovery Limited, Nottingham NG1 1GR, United Kingdom
| | | | - Tam Vo
- Department
of Chemistry and Center for Biotechnology and Drug Design, Georgia State University, Atlanta, Georgia 30303-3083, United States
| | - W. David Wilson
- Department
of Chemistry and Center for Biotechnology and Drug Design, Georgia State University, Atlanta, Georgia 30303-3083, United States
| | - Saadia A. Karim
- Cancer
Research UK Beatson Institute, Glasgow G61 1BD, United
Kingdom
| | - Jennifer P. Morton
- Cancer
Research UK Beatson Institute, Glasgow G61 1BD, United
Kingdom
- Institute
of Cancer Sciences, University of Glasgow, Glasgow G61 1QH, United Kingdom
| | - Stephen Neidle
- School
of Pharmacy, University College London, London WC1N 1AX, United Kingdom
| |
Collapse
|
6
|
Molecular dynamics simulations of G-quadruplexes: The basic principles and their application to folding and ligand binding. ANNUAL REPORTS IN MEDICINAL CHEMISTRY 2020. [DOI: 10.1016/bs.armc.2020.04.002] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
|
7
|
Kruse H, Šponer J. Revisiting the Potential Energy Surface of the Stacked Cytosine Dimer: FNO-CCSD(T) Interaction Energies, SAPT Decompositions, and Benchmarking. J Phys Chem A 2019; 123:9209-9222. [PMID: 31560201 DOI: 10.1021/acs.jpca.9b05940] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Nucleobase stacking interactions are crucial for the stability of nucleic acids. This study investigates base stacking energies of the cytosine homodimer in different configurations, including intermolecular separation plots, detailed twist dependence, and displaced structures. Highly accurate ab initio quantum chemical single point energies using an energy function based on MP2 complete basis set extrapolation ([6 → 7]ZaPa-NR) and a CCSD(T)/cc-pVTZ-F12 high-level correction are presented as new reference data, providing the most accurate stacking energies of nucleobase dimers currently available. Accurate SAPT2+(3)δMP2 energy decomposition is used to obtain detailed insights into the nature of base stacking interactions at varying vertical distances and twist values. The ab initio symmetry adapted perturbation theory (SAPT) energy decomposition suggests that the base stacking originates from an intricate interplay between dispersion attraction, short-range exchange-repulsion, and Coulomb interaction. The interpretation of the SAPT data is a complex issue as key energy terms vary substantially in the region of optimal (low energy) base stacking geometries. Thus, attempts to highlight one leading stabilizing SAPT base stacking term may be misleading and the outcome strongly depends on the used geometries within the range of geometries sampled in nucleic acids upon thermal fluctuations. Modern dispersion-corrected density functional theory (among them DSD-BLYP-D3, ωB97M-V, and ωB97M-D3BJ) is benchmarked and often reaches up to spectroscopic accuracy (below 1 kJ/mol). The classical AMBER force field is benchmarked with multiple different sets of point-charges (e.g. HF, DFT, and MP2-based) and is found to produce reasonable agreement with the benchmark data.
Collapse
Affiliation(s)
- Holger Kruse
- Institute of Biophysics of the Czech Academy of Sciences , Královopolská 135 , CZ-61265 Brno , Czech Republic
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences , Královopolská 135 , CZ-61265 Brno , Czech Republic.,Central European Institute of Technology , Masaryk University , Kamenice 753/5 , 62500 Brno , Czech Republic
| |
Collapse
|
8
|
Sutyak KB, Lee W, Zavalij PV, Gutierrez O, Davis JT. Templating and Catalyzing [2+2] Photocycloaddition in Solution Using a Dynamic G-Quadruplex. Angew Chem Int Ed Engl 2018; 57:17146-17150. [PMID: 30395701 DOI: 10.1002/anie.201811202] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2018] [Indexed: 11/08/2022]
Abstract
We describe a templating/covalent capture strategy that enables photochemical formation of 8 cyclobutanes in one noncovalent assembly. This process was characterized by experiment and quantum mechanical/molecular mechanics (ONIOM) calculations. Thus, KI and 16 units of 5'-cinnamate guanosine form a G-quadruplex where C=C π bonds in neighboring G4 -quartets are separated by 3.3 Å, enabling [2+2] photocycloaddition in solution. This reaction is high-yielding (>90 %), regio- and diastereoselective. Since all components are in dynamic equilibrium this photocycloaddition is catalytic in K+ .
Collapse
Affiliation(s)
- Keith B Sutyak
- Department of Chemistry and Biotechnology, University of Maryland, College Park, MD, 20742, USA
| | - Wes Lee
- Department of Chemistry and Biotechnology, University of Maryland, College Park, MD, 20742, USA
| | - Peter V Zavalij
- Department of Chemistry and Biotechnology, University of Maryland, College Park, MD, 20742, USA
| | - Osvaldo Gutierrez
- Department of Chemistry and Biotechnology, University of Maryland, College Park, MD, 20742, USA
| | - Jeffery T Davis
- Department of Chemistry and Biotechnology, University of Maryland, College Park, MD, 20742, USA
| |
Collapse
|
9
|
Sutyak KB, Lee W, Zavalij PV, Gutierrez O, Davis JT. Templating and Catalyzing [2+2] Photocycloaddition in Solution Using a Dynamic G‐Quadruplex. Angew Chem Int Ed Engl 2018. [DOI: 10.1002/ange.201811202] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Affiliation(s)
- Keith B. Sutyak
- Department of Chemistry and Biotechnology University of Maryland College Park MD 20742 USA
| | - Wes Lee
- Department of Chemistry and Biotechnology University of Maryland College Park MD 20742 USA
| | - Peter V. Zavalij
- Department of Chemistry and Biotechnology University of Maryland College Park MD 20742 USA
| | - Osvaldo Gutierrez
- Department of Chemistry and Biotechnology University of Maryland College Park MD 20742 USA
| | - Jeffery T. Davis
- Department of Chemistry and Biotechnology University of Maryland College Park MD 20742 USA
| |
Collapse
|
10
|
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
|
11
|
The estimation of H-bond and metal ion-ligand interaction energies in the G-Quadruplex ⋯ Mn+ complexes. J Mol Struct 2018. [DOI: 10.1016/j.molstruc.2018.02.023] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
|
12
|
Š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
|
13
|
Harikrishna S, Kotaru S, Pradeepkumar PI. Ligand-induced conformational preorganization of loops of c-MYC G-quadruplex DNA and its implications in structure-specific drug design. MOLECULAR BIOSYSTEMS 2018. [PMID: 28650023 DOI: 10.1039/c7mb00175d] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
Stabilization of a G-quadruplex (G4) DNA structure in the proto-oncogene c-MYC using small molecule ligands has emerged as an attractive strategy for the development of anticancer therapeutics. To understand the subtle structural changes in the G4 structure upon ligand binding, molecular dynamics (MD) simulations of c-MYC G4 DNA were carried out in a complex with six different potent ligands: 3AQN, 6AQN, 3APN, 360A, Nap-Et, and Nap-Pr. The results show that the ligands 3AQN, 6AQN, 3APN, and 360A stabilize the G4 structure by making stacking interactions with the top quartet. On the other hand, Nap-Et and Nap-Pr bind at the groove of the G4 structure. These groove binding ligands make crucial H-bond contacts with the guanines and electrostatic interactions with the phosphate backbone. Two-dimensional dynamic correlation maps unraveled the ligand-induced correlated motions between the guanines in the quartet and a di-nucleotide present in the propeller loop-2 of the G4 structure. Cluster analysis and ONIOM calculations revealed the structural dynamics in the loop of the quadruplex upon ligand binding. Overall, the results from the present study suggest that engineering specific contacts with the propeller loop can be an efficient way to design c-MYC G4-specific ligands.
Collapse
Affiliation(s)
- S Harikrishna
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai - 400076, India.
| | - Saikiran Kotaru
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai - 400076, India.
| | - P I Pradeepkumar
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai - 400076, India.
| |
Collapse
|
14
|
Šponer JE, Szabla R, Góra RW, Saitta AM, Pietrucci F, Saija F, Di Mauro E, Saladino R, Ferus M, Civiš S, Šponer J. Prebiotic synthesis of nucleic acids and their building blocks at the atomic level - merging models and mechanisms from advanced computations and experiments. Phys Chem Chem Phys 2018; 18:20047-66. [PMID: 27136968 DOI: 10.1039/c6cp00670a] [Citation(s) in RCA: 42] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
The origin of life on Earth is one of the most fascinating questions of contemporary science. Extensive research in the past decades furnished diverse experimental proposals for the emergence of first informational polymers that could form the basis of the early terrestrial life. Side by side with the experiments, the fast development of modern computational chemistry methods during the last 20 years facilitated the use of in silico modelling tools to complement the experiments. Modern computations can provide unique atomic-level insights into the structural and electronic aspects as well as the energetics of key prebiotic chemical reactions. Many of these insights are not directly obtainable from the experimental techniques and the computations are thus becoming indispensable for proper interpretation of many experiments and for qualified predictions. This review illustrates the synergy between experiment and theory in the origin of life research focusing on the prebiotic synthesis of various nucleic acid building blocks and on the self-assembly of nucleotides leading to the first functional oligonucleotides.
Collapse
Affiliation(s)
- Judit E Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, CZ-612 65 Brno, Czech Republic. and CEITEC - Central European Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, CZ-62500 Brno, Czech Republic
| | - Rafał Szabla
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, CZ-612 65 Brno, Czech Republic.
| | - Robert W Góra
- Theoretical Chemistry Group, Institute of Physical and Theoretical Chemistry, Wrocław University of Technology, Wybrzeże Wyspiańskiego 27, 50-370 Wrocław, Poland
| | - A Marco Saitta
- Sorbonne Universités, Université Pierre et Marie Curie Paris 6, CNRS, Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie, Muséum National d'Histoire Naturelle, Institut de Recherche pour le Développement, UMR 7590, F-75005 Paris, France
| | - Fabio Pietrucci
- Sorbonne Universités, Université Pierre et Marie Curie Paris 6, CNRS, Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie, Muséum National d'Histoire Naturelle, Institut de Recherche pour le Développement, UMR 7590, F-75005 Paris, France
| | - Franz Saija
- CNR-IPCF, Viale Ferdinando Stagno d'Alcontres 37, 98158 Messina, Italy
| | - Ernesto Di Mauro
- Dipartimento di Biologia e Biotecnologie "Charles Darwin", "Sapienza" Università di Roma, Piazzale Aldo Moro 5, Rome 00185, Italy
| | - Raffaele Saladino
- Dipartimento di Scienze Ecologiche e Biologiche Università della Tuscia, Via San Camillo De Lellis, 01100 Viterbo, Italy
| | - Martin Ferus
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, Dolejškova 3, CZ-182 23 Prague 8, Czech Republic
| | - Svatopluk Civiš
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, Dolejškova 3, CZ-182 23 Prague 8, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, CZ-612 65 Brno, Czech Republic. and CEITEC - Central European Institute of Technology, Masaryk University, Campus Bohunice, Kamenice 5, CZ-62500 Brno, Czech Republic
| |
Collapse
|
15
|
Mlýnský V, Kührová P, Jurečka P, Šponer J, Otyepka M, Banáš P. Mapping the Chemical Space of the RNA Cleavage and Its Implications for Ribozyme Catalysis. J Phys Chem B 2017; 121:10828-10840. [DOI: 10.1021/acs.jpcb.7b09129] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Vojtěch Mlýnský
- Regional Centre
of Advanced Technologies and Materials, Department of Physical Chemistry,
Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
- Scuola Internazionale Superiore di Studi Avanzati (SISSA), via
Bonomea 265, 34136 Trieste, Italy
| | - Petra Kührová
- Regional Centre
of Advanced Technologies and Materials, Department of Physical Chemistry,
Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Petr Jurečka
- Regional Centre
of Advanced Technologies and Materials, Department of Physical Chemistry,
Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Jiří Šponer
- Regional Centre
of Advanced Technologies and Materials, Department of Physical Chemistry,
Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
- Institute of Biophysics of the Czech Academy of Sciences, Kralovopolská 135, 612 65 Brno, Czech Republic
| | - Michal Otyepka
- Regional Centre
of Advanced Technologies and Materials, Department of Physical Chemistry,
Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Pavel Banáš
- Regional Centre
of Advanced Technologies and Materials, Department of Physical Chemistry,
Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
- Institute of Biophysics of the Czech Academy of Sciences, Kralovopolská 135, 612 65 Brno, Czech Republic
| |
Collapse
|
16
|
Alcolea Palafox M, Rastogi V, Singh S. Effect of the sulphur atom on geometry and spectra of the biomolecule 2-thiouracil and in the WC base pair 2-thiouridine-adenosine. Influence of water in the first hydration shell. J Biomol Struct Dyn 2017; 36:1225-1254. [DOI: 10.1080/07391102.2017.1318304] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Affiliation(s)
- M. Alcolea Palafox
- Facultad de Ciencias Químicas, Departamento de Química-Fisica1, Universidad Complutense, Madrid 28040, Spain
| | - V.K. Rastogi
- R.D. Foundation Group of Institutions, NH-58, Kadrabad, Modinagar, Ghaziabad, India
- Indian Spectroscopy Society, KC 68/1, Old Kavinagar, Ghaziabad 201 002, India
| | - S.P. Singh
- Department of Physics, Dr B R Ambedkar Govt Degree College, Mainpuri, India
| |
Collapse
|
17
|
Yurenko YP, Novotný J, Marek R. Weak Supramolecular Interactions Governing Parallel and Antiparallel DNA Quadruplexes: Insights from Large-Scale Quantum Mechanics Analysis of Experimentally Derived Models. Chemistry 2017; 23:5573-5584. [PMID: 28225208 DOI: 10.1002/chem.201700236] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2017] [Indexed: 12/30/2022]
Abstract
The topology and energetics of guanine (G) quadruplexes is governed by supramolecular interactions within their strands. In this work, an extensive quantum mechanical (QM) study has been performed to analyze supramolecular interactions that shape the stems of (4+0) parallel (P) and (2+2) antiparallel (AP) quadruplex systems. The large-scale (≈400 atoms) models of P and AP were constructed from high-quality experimental structures. The results provide evidence that each of the P and AP structures is shaped by a distinct network of supramolecular interactions. Analysis of electron topological characteristics of hydrogen bonds in P and AP systems indicates that the P model benefits from stronger intratetrad hydrogen bonding. For intertetrad stacking interactions, both noncovalent interaction plot and energy decomposition analysis approaches suggest that the stem of the P quadruplex benefits more from stacking than that of the AP stem; the difference in energetic stabilization for the two topologies is about 10 %. Stronger hydrogen-bonding and stacking interactions in the stem of the P quadruplex, relative to those in the AP system, can be an important indicator to explain the experimental observations that guanine-rich oligonucleotides tend to form all-parallel stems with an all-anti orientation of nucleobases. However, in addition to intrinsic stabilization, partial desolvation effects, which affect the energetics and dynamics of the G-quadruplex folding process, call for further investigations.
Collapse
Affiliation(s)
- Yevgen P Yurenko
- CEITEC - Central European Institute of Technology, Masaryk University, Kamenice 5/A4, 625 00, Brno, Czech Republic
| | - Jan Novotný
- CEITEC - Central European Institute of Technology, Masaryk University, Kamenice 5/A4, 625 00, Brno, Czech Republic
| | - Radek Marek
- CEITEC - Central European Institute of Technology, Masaryk University, Kamenice 5/A4, 625 00, Brno, Czech Republic
| |
Collapse
|
18
|
Cerón-Carrasco JP, Jacquemin D. Exposing the G-quadruplex to electric fields: the role played by telomeres in the propagation of DNA errors. Phys Chem Chem Phys 2017; 19:9358-9365. [DOI: 10.1039/c7cp01034f] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
We use quantum calculations to assess the impact of external electric fields on the stability of G-quadruplex, a key structure in telomere functionality.
Collapse
Affiliation(s)
- José Pedro Cerón-Carrasco
- Bioinformatics and High Performance Computing Research Group (BIO-HPC)
- Universidad Católica San Antonio de Murcia (UCAM)
- Murcia
- Spain
| | - Denis Jacquemin
- CEISAM
- UMR CNRS 6230
- BP 92208
- Université de Nantes
- 44322 Nantes Cedex 3
| |
Collapse
|
19
|
Szabla R, Havrila M, Kruse H, Šponer J. Comparative Assessment of Different RNA Tetranucleotides from the DFT-D3 and Force Field Perspective. J Phys Chem B 2016; 120:10635-10648. [PMID: 27681853 DOI: 10.1021/acs.jpcb.6b07551] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Classical force field (FF) molecular dynamics (MD) simulations of RNA tetranucleotides have substantial problems in reproducing conformer populations indicated by NMR experiments. To provide more information about the possible sources of errors, we performed quantum mechanical (QM, TPSS-D3/def2-TZVP) and molecular mechanics (MM, AMBER parm99bsc0+χOL3) calculations of different r(CCCC), r(GACC), and r(UUUU) conformers obtained from explicit solvent MD simulations. Solvent effects in the static QM and MM calculations were mimicked using implicit solvent models (COSMO and Poisson-Boltzmann, respectively). The comparison of QM and MM geometries and energies revealed that the two methodologies provide qualitatively consistent results in most of the cases. Even though we found some differences, these were insufficient to indicate any systematic corrections of the RNA FF terms that could improve the performance of classical MD in simulating tetranucleotides. On the basis of these findings, we inferred that the overpopulation of intercalated conformers in the MD simulations of RNA tetramers, which were not observed experimentally, might be predominantly caused by imbalanced water-solvent and water-water interactions. Apart from the large-scale QM calculations performed to assess the performance of the AMBER FF, a representative spectrum of faster QM methods was tested.
Collapse
Affiliation(s)
- Rafał Szabla
- Institute of Biophysics, Academy of Sciences of the Czech Republic , Královopolská 135, CZ-61265 Brno, Czech Republic
| | - Marek Havrila
- Institute of Biophysics, Academy of Sciences of the Czech Republic , Královopolská 135, CZ-61265 Brno, Czech Republic.,CEITEC - Central European Institute of Technology, Masaryk University , Campus Bohunice, Kamenice 5, CZ-62500 Brno, Czech Republic
| | - Holger Kruse
- Institute of Biophysics, Academy of Sciences of the Czech Republic , Královopolská 135, CZ-61265 Brno, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic , Královopolská 135, CZ-61265 Brno, Czech Republic.,CEITEC - Central European Institute of Technology, Masaryk University , Campus Bohunice, Kamenice 5, CZ-62500 Brno, Czech Republic
| |
Collapse
|