1
|
Zubova EA, Strelnikov IA. Experimental detection of conformational transitions between forms of DNA: problems and prospects. Biophys Rev 2023; 15:1053-1078. [PMID: 37974981 PMCID: PMC10643659 DOI: 10.1007/s12551-023-01143-9] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2023] [Accepted: 09/06/2023] [Indexed: 11/19/2023] Open
Abstract
Under different conditions, the DNA double helix can take different geometric forms. Of the large number of its conformations, in addition to the "canonical" B form, the A, C, and Z forms are widely known, and the D, Hoogsteen, and X forms are less known. DNA locally takes the A, C, and Z forms in the cell, in complexes with proteins. We compare different methods for detecting non-canonical DNA conformations: X-ray, IR, and Raman spectroscopy, linear and circular dichroism in both the infrared and ultraviolet regions, as well as NMR (measurement of chemical shifts and their anisotropy, scalar and residual dipolar couplings and inter-proton distances from NOESY (nuclear Overhauser effect spectroscopy) data). We discuss the difficulties in applying these methods, the problems of theoretical interpretation of the experimental results, and the prospects for reliable identification of non-canonical DNA conformations.
Collapse
Affiliation(s)
- Elena A. Zubova
- N.N. Semenov Federal Research Center for Chemical Physics, Russian Academy of Sciences, 4 Kosygin St., Moscow, 119991 Russia
| | - Ivan A. Strelnikov
- N.N. Semenov Federal Research Center for Chemical Physics, Russian Academy of Sciences, 4 Kosygin St., Moscow, 119991 Russia
| |
Collapse
|
2
|
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
|
3
|
Westhof E, Watson ZL, Zirbel CL, Cate JHD. Anionic G•U pairs in bacterial ribosomal rRNAs. RNA (NEW YORK, N.Y.) 2023; 29:1069-1076. [PMID: 37068913 DOI: 10.1261/rna.079583.123] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/08/2023] [Accepted: 04/05/2023] [Indexed: 06/18/2023]
Abstract
Wobble GU pairs (or G•U) occur frequently within double-stranded RNA helices interspersed between standard G=C and A-U Watson-Crick pairs. Another type of G•U pair interacting via their Watson-Crick edges has been observed in the A site of ribosome structures between a modified U34 in the tRNA anticodon triplet and G + 3 in the mRNA. In such pairs, the electronic structure of the U is changed with a negative charge on N3(U), resulting in two H-bonds between N1(G)…O4(U) and N2(G)…N3(U). Here, we report that such pairs occur in other highly conserved positions in ribosomal RNAs of bacteria in the absence of U modification. An anionic cis Watson-Crick G•G pair is also observed and well conserved in the small subunit. These pairs are observed in tightly folded regions.
Collapse
Affiliation(s)
- Eric Westhof
- Architecture et Réactivité de l'ARN, Université de Strasbourg, Institut de biologie moléculaire et cellulaire du CNRS, F-67084 Strasbourg, France
| | - Zoe L Watson
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- California Institute for Quantitative Biosciences (QB3), University of California, Berkeley, California 94720, USA
| | - Craig L Zirbel
- Department of Mathematics and Statistics, Bowling Green State University, Bowling Green, Ohio 43403, USA
| | - Jamie H D Cate
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Department of Molecular and Cell Biology, University of California, Berkeley, California 94720, USA
- Molecular Biophysics and Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
- Innovative Genomics Institute, University of California, Berkeley, California 94720, USA
| |
Collapse
|
4
|
Abstract
There is growing interest in therapeutic intervention that targets disease-relevant RNAs using small molecules. While there have been some successes in RNA-targeted small-molecule discovery, a deeper understanding of structure-activity relationships in pursuing these targets has remained elusive. One of the best-studied tertiary-structured RNAs is the theophylline aptamer, which binds theophylline with high affinity and selectivity. Although not a drug target, this aptamer has had many applications, especially pertaining to genetic control circuits. Heretofore, no compound has been shown to bind the theophylline aptamer with greater affinity than theophylline itself. However, by carrying out a high-throughput screen of low-molecular-weight compounds, several unique hits were identified that are chemically distinct from theophylline and bind with up to 340-fold greater affinity. Multiple atomic-resolution X-ray crystal structures were determined to investigate the binding mode of theophylline and four of the best hits. These structures reveal both the rigidity of the theophylline aptamer binding pocket and the opportunity for other ligands to bind more tightly in this pocket by forming additional hydrogen-bonding interactions. These results give encouragement that the same approaches to drug discovery that have been applied so successfully to proteins can also be applied to RNAs.
Collapse
|
5
|
Kiliszek A, Pluta M, Bejger M, Rypniewski W. Structure and thermodynamics of a UGG motif interacting with Ba2+ and other metal ions: accommodating changes in the RNA structure and the presence of a G(syn)-G(syn) pair. RNA (NEW YORK, N.Y.) 2022; 29:rna.079414.122. [PMID: 36319090 PMCID: PMC9808570 DOI: 10.1261/rna.079414.122] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/11/2022] [Accepted: 10/11/2022] [Indexed: 06/16/2023]
Abstract
The self-complementary triplet 5'UGG3'/5'UGG3' is a particular structural motif containing noncanonical G-G pair and two U·G wobble pairs. It constitutes a specific structural and electrostatic environment attracting metal ions, particularly Ba2+ ions. Crystallographic research has shown that two Ba2+ cations are located in the major groove of the helix and interact directly with the UGG triplet. A comparison with the unliganded structure has revealed global changes in the RNA structure in the presence of metal ions, whereas thermodynamic measurements have shown increased stability. Moreover, in the structure with Ba2+, an unusual noncanonical G(syn)-G(syn) pair is observed instead of the common G(anti)-G(syn). We further elucidate the metal binding properties of the UGG/UGG triplet by performing crystallographic and thermodynamic studies using DSC and UV melting with other metal ions. The results explain the preferences of the UGG sequence for Ba2+ cations and point to possible applications of this metal-binding propensity.
Collapse
|
6
|
Ramachandran V, Mainan A, Roy S. Dynamic effects of the spine of hydrated magnesium on viral RNA pseudoknot structure. Phys Chem Chem Phys 2022; 24:24570-24581. [PMID: 36193826 DOI: 10.1039/d2cp01075e] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
In the cellular environment, a viral RNA Pseudoknot (PK) structure is responsive to its prevailing ion atmosphere created by a mixture of monovalent and divalent cations. We investigate the influence of such a mixed-salt environment on RNA-PK structure at an atomic resolution through three sets of 1.5 μs explicit solvent molecular dynamics (MD) simulations and also by building a dynamic counterion-condensation (DCC) model at varying divalent Mg2+ concentrations. The DCC model includes explicit interaction of the ligand and adjacent chelated Mg2+ by extending the recently developed generalized Manning condensation model. Its implementation within an all-atom structure-based molecular dynamics framework bolsters its opportunity to explore large-length scale and long-timescale phenomena associated with complex RNA systems immersed in its dynamic ion environment. In the present case of RNA-PK, both explicit MD and DCC simulations reveal a spine of hydrated-Mg2+ to induce stem-I and stem-II closure where the minor groove between these stems is akin to breathing. Mg2+ mediated minor-groove narrowing is coupled with local base-flipping dynamics of a base triple and quadruple, changing the stem structure of such RNA. Contrary to the conversational view of the indispensable need for Mg2+ for the tertiary structure of RNA, the study warns about the plausible detrimental effect of specific Mg2+-phosphate interactions on the RNA-PK structure beyond a certain concentration of Mg2+.
Collapse
Affiliation(s)
- Vysakh Ramachandran
- Department of Chemical Sciences, Indian Institute of Science Education and Research Kolkata, Mohanpur, West Bengal 741246, India.
| | - Avijit Mainan
- Department of Chemical Sciences, Indian Institute of Science Education and Research Kolkata, Mohanpur, West Bengal 741246, India.
| | - Susmita Roy
- Department of Chemical Sciences, Indian Institute of Science Education and Research Kolkata, Mohanpur, West Bengal 741246, India.
| |
Collapse
|
7
|
Kolev SK, St. Petkov P, Milenov TI, Vayssilov GN. Sodium and Magnesium Ion Location at the Backbone and at the Nucleobase of RNA: Ab Initio Molecular Dynamics in Water Solution. ACS OMEGA 2022; 7:23234-23244. [PMID: 35847262 PMCID: PMC9280761 DOI: 10.1021/acsomega.2c01327] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
The interactions between Na+ or Mg2+ ions with different parts of single-stranded RNA molecules, namely, the oxygen atoms from the phosphate groups or the guanine base, in water solution have been studied using first-principles molecular dynamics. Sodium ions were found to be much more mobile than Mg2+ ions and readily underwent transitions between a state directly bonded to RNA oxygen atoms and a completely solvated state. The inner solvation shell of Na+ ions fluctuated stochastically at a femtosecond timescale coordinating on average 5 oxygen atoms for bonded Na+ ions and 5.5 oxygen atoms for solvated Na+ ions. In contrast, the inner solvation shell of Mg2+ ions was stable in both RNA-bonded and completely solvated states. In both cases, Mg2+ ions coordinated 6 oxygen atoms from the inner solvation shell. Consistent with their stable solvation shells, Mg2+ ions were more effective than Na+ ions in stabilizing the RNA backbone conformation. The exclusion zones between the first and second solvation shells, solvation shell widths, and angles for binding to carbonyl oxygen of guanine for solvated Na+ or Mg2+ ions exhibited a number of quantitative differences when compared with RNA crystallographic data. The presented results support the distinct capacity of Mg2+ ions to support the RNA structure not only in the crystal phase but also in the dynamic water environment both on the side of the phosphate moiety and on the side of the nucleobase.
Collapse
Affiliation(s)
- Stefan K. Kolev
- Institute
of Electronics, Bulgarian Academy of Sciences, 72 Tzarigradsko Chaussee Blvd., Sofia 1784, Bulgaria
| | - Petko St. Petkov
- Faculty
of Chemistry and Pharmacy, University of
Sofia, Boulevard James
Bouchier 1, Sofia 1126, Bulgaria
| | - Teodor I. Milenov
- Institute
of Electronics, Bulgarian Academy of Sciences, 72 Tzarigradsko Chaussee Blvd., Sofia 1784, Bulgaria
| | - Georgi N. Vayssilov
- Faculty
of Chemistry and Pharmacy, University of
Sofia, Boulevard James
Bouchier 1, Sofia 1126, Bulgaria
| |
Collapse
|
8
|
Hydrated magnesium (II) with H-bonded hexamethylenetetramine and co-ligands: synthesis, structural characterization, thermal decomposition, and hirshfeld surface analysis. CHEMICAL PAPERS 2022. [DOI: 10.1007/s11696-022-02322-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
9
|
Mair S, Erharter K, Renard E, Brillet K, Brunner M, Lusser A, Kreutz C, Ennifar E, Micura R. Towards a comprehensive understanding of RNA deamination: synthesis and properties of xanthosine-modified RNA. Nucleic Acids Res 2022; 50:6038-6051. [PMID: 35687141 PMCID: PMC9226506 DOI: 10.1093/nar/gkac477] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2022] [Revised: 04/28/2022] [Accepted: 05/20/2022] [Indexed: 12/25/2022] Open
Abstract
Nucleobase deamination, such as A-to-I editing, represents an important posttranscriptional modification of RNA. When deamination affects guanosines, a xanthosine (X) containing RNA is generated. However, the biological significance and chemical consequences on RNA are poorly understood. We present a comprehensive study on the preparation and biophysical properties of X-modified RNA. Thermodynamic analyses revealed that base pairing strength is reduced to a level similar to that observed for a G•U replacement. Applying NMR spectroscopy and X-ray crystallography, we demonstrate that X can form distinct wobble geometries with uridine depending on the sequence context. In contrast, X pairing with cytidine occurs either through wobble geometry involving protonated C or in Watson–Crick-like arrangement. This indicates that the different pairing modes are of comparable stability separated by low energetic barriers for switching. Furthermore, we demonstrate that the flexible pairing properties directly affect the recognition of X-modified RNA by reverse transcription enzymes. Primer extension assays and PCR-based sequencing analysis reveal that X is preferentially read as G or A and that the ratio depends on the type of reverse transcriptase. Taken together, our results elucidate important properties of X-modified RNA paving the way for future studies on its biological significance.
Collapse
Affiliation(s)
- Stefan Mair
- Institute of Organic Chemistry, Center for Molecular Biosciences Innsbruck, University of Innsbruck, Innsbruck 6020, Austria
| | - Kevin Erharter
- Institute of Organic Chemistry, Center for Molecular Biosciences Innsbruck, University of Innsbruck, Innsbruck 6020, Austria
| | - Eva Renard
- Architecture et Réactivité de l'ARN - CNRS UPR 9002, Institut de Biologie Moléculaire et Cellulaire, Université de Strasbourg, 67000 Strasbourg, France
| | - Karl Brillet
- Architecture et Réactivité de l'ARN - CNRS UPR 9002, Institut de Biologie Moléculaire et Cellulaire, Université de Strasbourg, 67000 Strasbourg, France
| | - Melanie Brunner
- Institute of Molecular Biology, Biocenter, Medical University of Innsbruck, Innsbruck 6020, Austria
| | - Alexandra Lusser
- Institute of Molecular Biology, Biocenter, Medical University of Innsbruck, Innsbruck 6020, Austria
| | - Christoph Kreutz
- Institute of Organic Chemistry, Center for Molecular Biosciences Innsbruck, University of Innsbruck, Innsbruck 6020, Austria
| | - Eric Ennifar
- Architecture et Réactivité de l'ARN - CNRS UPR 9002, Institut de Biologie Moléculaire et Cellulaire, Université de Strasbourg, 67000 Strasbourg, France
| | - Ronald Micura
- Institute of Organic Chemistry, Center for Molecular Biosciences Innsbruck, University of Innsbruck, Innsbruck 6020, Austria
| |
Collapse
|
10
|
Vázquez-Durán A, Téllez-Isaías G, Hernández-Rodríguez M, Ruvalcaba RM, Martínez J, Nicolás-Vázquez MI, Aceves-Hernández JM, Méndez-Albores A. The Ability of Chlorophyll to Trap Carcinogen Aflatoxin B 1: A Theoretical Approach. Int J Mol Sci 2022; 23:6068. [PMID: 35682746 PMCID: PMC9181093 DOI: 10.3390/ijms23116068] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2022] [Revised: 05/19/2022] [Accepted: 05/26/2022] [Indexed: 01/27/2023] Open
Abstract
The coordination of one and two aflatoxin B1 (AFB1, a potent carcinogen) molecules with chlorophyll a (chl a) was studied at a theoretical level. Calculations were performed using the M06-2X method in conjunction with the 6-311G(d,p) basis set, in both gas and water phases. The molecular electrostatic potential map shows the chemical activity of various sites of the AFB1 and chl a molecules. The energy difference between molecular orbitals of AFB1 and chl a allowed for the establishment of an intermolecular interaction. A charge transfer from AFB1 to the central cation of chl a was shown. The energies of the optimized structures for chl a show two configurations, unfolded and folded, with a difference of 15.41 kcal/mol. Chl a appeared axially coordinated to the plane (α-down or β-up) of the porphyrin moiety, either with the oxygen atom of the ketonic group, or with the oxygen atom of the lactone moiety of AFB1. The complexes of maximum stability were chl a 1-α-E-AFB1 and chl a 2-β-E-AFB1, at -36.4 and -39.2 kcal/mol, respectively. Additionally, with two AFB1 molecules were chl a 1-D-2AFB1 and chl a 2-E-2AFB1, at -60.0 and -64.8 kcal/mol, respectively. Finally, biosorbents containing chlorophyll could improve AFB1 adsorption.
Collapse
Affiliation(s)
- Alma Vázquez-Durán
- Unidad de Investigación Multidisciplinaria L14 (Alimentos, Micotoxinas, y Micotoxicosis), Facultad de Estudios Superiores Cuautitlán, Universidad Nacional Autónoma de México, Cuautitlán Izcalli, Estado de Mexico 54714, Mexico; (A.V.-D.); (J.M.A.-H.); (A.M.-A.)
| | | | - Maricarmen Hernández-Rodríguez
- Laboratorio de Cultivo Celular, Sección de Posgrado e Investigación, Escuela Superior de Medicina, Instituto Politécnico Nacional, Ciudad de Mexico 11340, Mexico;
| | - René Miranda Ruvalcaba
- Departamento de Ciencias Químicas, Facultad de Estudios Superiores Cuautitlán Campo 1, Universidad Nacional Autónoma de México, Avenida 1o de Mayo s/n, Colonia Santa María las Torres, Cuautitlán Izcalli, Estado de Mexico 54740, Mexico;
| | - Joel Martínez
- Facultad de Ciencias Químicas, Universidad Autónoma de San Luis Potosí, San Luis Potosi 78210, Mexico
| | - María Inés Nicolás-Vázquez
- Departamento de Ciencias Químicas, Facultad de Estudios Superiores Cuautitlán Campo 1, Universidad Nacional Autónoma de México, Avenida 1o de Mayo s/n, Colonia Santa María las Torres, Cuautitlán Izcalli, Estado de Mexico 54740, Mexico;
| | - Juan Manuel Aceves-Hernández
- Unidad de Investigación Multidisciplinaria L14 (Alimentos, Micotoxinas, y Micotoxicosis), Facultad de Estudios Superiores Cuautitlán, Universidad Nacional Autónoma de México, Cuautitlán Izcalli, Estado de Mexico 54714, Mexico; (A.V.-D.); (J.M.A.-H.); (A.M.-A.)
| | - Abraham Méndez-Albores
- Unidad de Investigación Multidisciplinaria L14 (Alimentos, Micotoxinas, y Micotoxicosis), Facultad de Estudios Superiores Cuautitlán, Universidad Nacional Autónoma de México, Cuautitlán Izcalli, Estado de Mexico 54714, Mexico; (A.V.-D.); (J.M.A.-H.); (A.M.-A.)
| |
Collapse
|
11
|
Neumann J, Schwierz N. Artificial Intelligence Resolves Kinetic Pathways of Magnesium Binding to RNA. J Chem Theory Comput 2022; 18:1202-1212. [PMID: 35084846 PMCID: PMC8830046 DOI: 10.1021/acs.jctc.1c00752] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
Magnesium is an indispensable cofactor in countless vital processes. In order to understand its functional role, the characterization of the binding pathways to biomolecules such as RNA is crucial. Despite the importance, a molecular description is still lacking since the transition from the water-mediated outer-sphere to the direct inner-sphere coordination is on the millisecond time scale and therefore out of reach for conventional simulation techniques. To fill this gap, we use transition path sampling to resolve the binding pathways and to elucidate the role of the solvent in the binding process. The results reveal that the molecular void provoked by the leaving phosphate oxygen of the RNA is immediately filled by an entering water molecule. In addition, water molecules from the first and second hydration shell couple to the concerted exchange. To capture the intimate solute-solvent coupling, we perform a committor analysis as the basis for a machine learning algorithm that derives the optimal deep learning model from thousands of scanned architectures using hyperparameter tuning. The results reveal that the properly optimized deep network architecture recognizes the important solvent structures, extracts the relevant information, and predicts the commitment probability with high accuracy. Our results provide detailed insights into the solute-solvent coupling which is ubiquitous for kosmotropic ions and governs a large variety of biochemical reactions in aqueous solutions.
Collapse
Affiliation(s)
- Jan Neumann
- Allianz Global Investors GmbH, Bockenheimer Landstrasse 42, 60323 Frankfurt am Main, Germany
| | - Nadine Schwierz
- Department of Theoretical Biophysics, Max-Planck-Institute of Biophysics, 60438 Frankfurt am Main, Germany
| |
Collapse
|
12
|
Tirumalai MR, Rivas M, Tran Q, Fox GE. The Peptidyl Transferase Center: a Window to the Past. Microbiol Mol Biol Rev 2021; 85:e0010421. [PMID: 34756086 PMCID: PMC8579967 DOI: 10.1128/mmbr.00104-21] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023] Open
Abstract
In his 2001 article, "Translation: in retrospect and prospect," the late Carl Woese made a prescient observation that there was a need for the then-current view of translation to be "reformulated to become an all-embracing perspective about which 21st century Biology can develop" (RNA 7:1055-1067, 2001, https://doi.org/10.1017/s1355838201010615). The quest to decipher the origins of life and the road to the genetic code are both inextricably linked with the history of the ribosome. After over 60 years of research, significant progress in our understanding of how ribosomes work has been made. Particularly attractive is a model in which the ribosome may facilitate an ∼180° rotation of the CCA end of the tRNA from the A-site to the P-site while the acceptor stem of the tRNA would then undergo a translation from the A-site to the P-site. However, the central question of how the ribosome originated remains unresolved. Along the path from a primitive RNA world or an RNA-peptide world to a proto-ribosome world, the advent of the peptidyl transferase activity would have been a seminal event. This functionality is now housed within a local region of the large-subunit (LSU) rRNA, namely, the peptidyl transferase center (PTC). The PTC is responsible for peptide bond formation during protein synthesis and is usually considered to be the oldest part of the modern ribosome. What is frequently overlooked is that by examining the origins of the PTC itself, one is likely going back even further in time. In this regard, it has been proposed that the modern PTC originated from the association of two smaller RNAs that were once independent and now comprise a pseudosymmetric region in the modern PTC. Could such an association have survived? Recent studies have shown that the extant PTC is largely depleted of ribosomal protein interactions. It is other elements like metallic ion coordination and nonstandard base/base interactions that would have had to stabilize the association of RNAs. Here, we present a detailed review of the literature focused on the nature of the extant PTC and its proposed ancestor, the proto-ribosome.
Collapse
Affiliation(s)
- Madhan R. Tirumalai
- Department of Biology and Biochemistry, University of Houston, Houston, Texas, USA
| | - Mario Rivas
- Department of Biology and Biochemistry, University of Houston, Houston, Texas, USA
| | - Quyen Tran
- Department of Biology and Biochemistry, University of Houston, Houston, Texas, USA
| | - George E. Fox
- Department of Biology and Biochemistry, University of Houston, Houston, Texas, USA
| |
Collapse
|
13
|
Largy E, König A, Ghosh A, Ghosh D, Benabou S, Rosu F, Gabelica V. Mass Spectrometry of Nucleic Acid Noncovalent Complexes. Chem Rev 2021; 122:7720-7839. [PMID: 34587741 DOI: 10.1021/acs.chemrev.1c00386] [Citation(s) in RCA: 39] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Nucleic acids have been among the first targets for antitumor drugs and antibiotics. With the unveiling of new biological roles in regulation of gene expression, specific DNA and RNA structures have become very attractive targets, especially when the corresponding proteins are undruggable. Biophysical assays to assess target structure as well as ligand binding stoichiometry, affinity, specificity, and binding modes are part of the drug development process. Mass spectrometry offers unique advantages as a biophysical method owing to its ability to distinguish each stoichiometry present in a mixture. In addition, advanced mass spectrometry approaches (reactive probing, fragmentation techniques, ion mobility spectrometry, ion spectroscopy) provide more detailed information on the complexes. Here, we review the fundamentals of mass spectrometry and all its particularities when studying noncovalent nucleic acid structures, and then review what has been learned thanks to mass spectrometry on nucleic acid structures, self-assemblies (e.g., duplexes or G-quadruplexes), and their complexes with ligands.
Collapse
Affiliation(s)
- Eric Largy
- Univ. Bordeaux, CNRS, INSERM, ARNA, UMR 5320, U1212, IECB, F-33600 Pessac, France
| | - Alexander König
- Univ. Bordeaux, CNRS, INSERM, ARNA, UMR 5320, U1212, IECB, F-33600 Pessac, France
| | - Anirban Ghosh
- Univ. Bordeaux, CNRS, INSERM, ARNA, UMR 5320, U1212, IECB, F-33600 Pessac, France
| | - Debasmita Ghosh
- Univ. Bordeaux, CNRS, INSERM, ARNA, UMR 5320, U1212, IECB, F-33600 Pessac, France
| | - Sanae Benabou
- Univ. Bordeaux, CNRS, INSERM, ARNA, UMR 5320, U1212, IECB, F-33600 Pessac, France
| | - Frédéric Rosu
- Univ. Bordeaux, CNRS, INSERM, IECB, UMS 3033, F-33600 Pessac, France
| | - Valérie Gabelica
- Univ. Bordeaux, CNRS, INSERM, ARNA, UMR 5320, U1212, IECB, F-33600 Pessac, France
| |
Collapse
|
14
|
Cruz-León S, Grotz KK, Schwierz N. Extended magnesium and calcium force field parameters for accurate ion-nucleic acid interactions in biomolecular simulations. J Chem Phys 2021; 154:171102. [PMID: 34241062 DOI: 10.1063/5.0048113] [Citation(s) in RCA: 18] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Abstract
Magnesium and calcium play an essential role in the folding and function of nucleic acids. To correctly describe their interactions with DNA and RNA in biomolecular simulations, an accurate parameterization is crucial. In most cases, the ion parameters are optimized based on a set of experimental solution properties such as solvation free energies, radial distribution functions, water exchange rates, and activity coefficient derivatives. However, the transferability of such bulk-optimized ion parameters to quantitatively describe biomolecular systems is limited. Here, we extend the applicability of our previous bulk-optimized parameters by including experimental binding affinities toward the phosphate oxygen on nucleic acids. In particular, we systematically adjust the combination rules that are an integral part of the pairwise interaction potentials of classical force fields. This allows us to quantitatively describe specific ion binding to nucleic acids without changing the solution properties in the most simple and efficient way. We show the advancement of the optimized Lorentz combination rule for two representative nucleic acid systems. For double-stranded DNA, the optimized combination rule for Ca2+ significantly improves the agreement with experiments, while the standard combination rule leads to unrealistically distorted DNA structures. For the add A-riboswitch, the optimized combination rule for Mg2+ improves the structure of two specifically bound Mg2+ ions as judged by the experimental distance to the binding site. Including experimental binding affinities toward specific ion binding sites on biomolecules, therefore, provides a promising perspective to develop a more accurate description of metal cations for biomolecular simulations.
Collapse
Affiliation(s)
- Sergio Cruz-León
- Department of Theoretical Biophysics, Max Planck Institute of Biophysics, Max-von-Laue-Str. 3, 60438 Frankfurt am Main, Germany
| | - Kara K Grotz
- Department of Theoretical Biophysics, Max Planck Institute of Biophysics, Max-von-Laue-Str. 3, 60438 Frankfurt am Main, Germany
| | - Nadine Schwierz
- Department of Theoretical Biophysics, Max Planck Institute of Biophysics, Max-von-Laue-Str. 3, 60438 Frankfurt am Main, Germany
| |
Collapse
|
15
|
Grotz KK, Cruz-León S, Schwierz N. Optimized Magnesium Force Field Parameters for Biomolecular Simulations with Accurate Solvation, Ion-Binding, and Water-Exchange Properties. J Chem Theory Comput 2021; 17:2530-2540. [PMID: 33720710 PMCID: PMC8047801 DOI: 10.1021/acs.jctc.0c01281] [Citation(s) in RCA: 30] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2020] [Indexed: 12/31/2022]
Abstract
Magnesium ions play an essential role in many vital processes. To correctly describe their interactions in molecular dynamics simulations, an accurate parametrization is crucial. Despite the importance and considerable scientific effort, current force fields based on the commonly used 12-6 Lennard-Jones interaction potential fail to reproduce a variety of experimental solution properties. In particular, no parametrization exists so far that simultaneously reproduces the solvation free energy and the distance to the water oxygens in the first hydration shell. Moreover, current Mg2+ force fields significantly underestimate the rate of water exchange leading to unrealistically slow exchange kinetics. In order to make progress in the development of improved models, we systematically optimize the Mg2+ parameters in combination with the TIP3P water model in a much larger parameter space than previously done. The results show that a long-ranged interaction potential and modified Lorentz-Berthelot combination rules allow us to accurately reproduce multiple experimental properties including the solvation free energy, the distances to the oxygens of the first hydration shell, the hydration number, the activity coefficient derivative in MgCl2 solutions, the self-diffusion coefficient, and the binding affinity to the phosphate oxygen of RNA. Matching this broad range of thermodynamic properties, we present two sets of optimal parameters: MicroMg yields water exchange on the microsecond timescale in agreement with experiments. NanoMg yields water exchange on the nanosecond timescale facilitating the direct observation of ion-binding events. As shown for the example of the add A-riboswitch, the optimized parameters correctly reproduce the structure of specifically bound ions and permit the de novo prediction of Mg2+-binding sites in biomolecular simulations.
Collapse
Affiliation(s)
- Kara K. Grotz
- Department of Theoretical Biophysics, Max-Planck-Institute of Biophysics, Frankfurt am Main 60438, Germany
| | - Sergio Cruz-León
- Department of Theoretical Biophysics, Max-Planck-Institute of Biophysics, Frankfurt am Main 60438, Germany
| | - Nadine Schwierz
- Department of Theoretical Biophysics, Max-Planck-Institute of Biophysics, Frankfurt am Main 60438, Germany
| |
Collapse
|
16
|
Jones RA, Steckelberg AL, Vicens Q, Szucs MJ, Akiyama BM, Kieft JS. Different tertiary interactions create the same important 3D features in a distinct flavivirus xrRNA. RNA (NEW YORK, N.Y.) 2021; 27:54-65. [PMID: 33004436 PMCID: PMC7749634 DOI: 10.1261/rna.077065.120] [Citation(s) in RCA: 23] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/01/2020] [Accepted: 09/27/2020] [Indexed: 05/12/2023]
Abstract
During infection by a flavivirus (FV), cells accumulate noncoding subgenomic flavivirus RNAs (sfRNAs) that interfere with several antiviral pathways. These sfRNAs are formed by structured RNA elements in the 3' untranslated region (UTR) of the viral genomic RNA, which block the progression of host cell exoribonucleases that have targeted the viral RNA. Previous work on these exoribonuclease-resistant RNAs (xrRNAs) from mosquito-borne FVs revealed a specific three-dimensional fold with a unique topology in which a ring-like structure protectively encircles the 5' end of the xrRNA. Conserved nucleotides make specific tertiary interactions that support this fold. Examination of more divergent FVs reveals differences in their 3' UTR sequences, raising the question of whether they contain xrRNAs and if so, how they fold. To answer this, we demonstrated the presence of an authentic xrRNA in the 3' UTR of the Tamana bat virus (TABV) and solved its structure by X-ray crystallography. The structure reveals conserved features from previously characterized xrRNAs, but in the TABV version these features are created through a novel set of tertiary interactions not previously seen in xrRNAs. This includes two important A-C interactions, four distinct backbone kinks, several ordered Mg2+ ions, and a C+-G-C base triple. The discovery that the same overall architecture can be achieved by very different sequences and interactions in distantly related flaviviruses provides insight into the diversity of this type of RNA and will inform searches for undiscovered xrRNAs in viruses and beyond.
Collapse
MESH Headings
- 3' Untranslated Regions
- Animals
- Base Pairing
- Base Sequence
- Cations, Divalent
- Crystallography, X-Ray
- Encephalitis Virus, Murray Valley/genetics
- Encephalitis Virus, Murray Valley/metabolism
- Encephalitis Virus, Murray Valley/ultrastructure
- Exoribonucleases/chemistry
- Exoribonucleases/metabolism
- Flaviviridae/genetics
- Flaviviridae/metabolism
- Flaviviridae/ultrastructure
- Host-Pathogen Interactions/genetics
- Magnesium/chemistry
- Magnesium/metabolism
- RNA Folding
- RNA, Untranslated/chemistry
- RNA, Untranslated/genetics
- RNA, Untranslated/metabolism
- RNA, Viral/chemistry
- RNA, Viral/genetics
- RNA, Viral/metabolism
- Viruses, Unclassified/genetics
- Viruses, Unclassified/metabolism
- Viruses, Unclassified/ultrastructure
- Zika Virus/genetics
- Zika Virus/metabolism
- Zika Virus/ultrastructure
Collapse
Affiliation(s)
- Rachel A Jones
- Department of Biochemistry and Molecular Genetics, University of Colorado Denver School of Medicine, Aurora, Colorado 80045, USA
| | - Anna-Lena Steckelberg
- Department of Biochemistry and Molecular Genetics, University of Colorado Denver School of Medicine, Aurora, Colorado 80045, USA
| | - Quentin Vicens
- Department of Biochemistry and Molecular Genetics, University of Colorado Denver School of Medicine, Aurora, Colorado 80045, USA
| | - Matthew J Szucs
- Department of Biochemistry and Molecular Genetics, University of Colorado Denver School of Medicine, Aurora, Colorado 80045, USA
| | - Benjamin M Akiyama
- Department of Biochemistry and Molecular Genetics, University of Colorado Denver School of Medicine, Aurora, Colorado 80045, USA
| | - Jeffrey S Kieft
- Department of Biochemistry and Molecular Genetics, University of Colorado Denver School of Medicine, Aurora, Colorado 80045, USA
- RNA BioScience Initiative, University of Colorado Denver School of Medicine, Aurora, Colorado 80045, USA
| |
Collapse
|
17
|
Auffinger P, Ennifar E, D'Ascenzo L. Deflating the RNA Mg 2+ bubble. Stereochemistry to the rescue! RNA (NEW YORK, N.Y.) 2020; 27:rna.076067.120. [PMID: 33268500 PMCID: PMC7901845 DOI: 10.1261/rna.076067.120] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/25/2020] [Accepted: 11/20/2020] [Indexed: 05/03/2023]
Abstract
Proper evaluation of the ionic structure of biomolecular systems through X ray and cryo-EM techniques remains challenging but is essential for advancing our understanding of the underlying structure/activity/solvent relationships. However, numerous studies overestimate the number of Mg2+ in deposited structures due to assignment errors finding their origin in improper consideration of stereochemical rules. Herein, to tackle such issues, we re-evaluate the PDBid 6QNR and 6SJ6 models of the ribosome ionic structure. We establish that stereochemical principles need to be carefully pondered when evaluating ion binding features, even when K+ anomalous signals are available as it is the case for the 6QNR PDB entry. For ribosomes, assignment errors can result in misleading conceptions of their solvent structure. For instance, present stereochemical analysis result in a significant decrease of the number of assigned Mg2+ in 6QNR, suggesting that K+ and not Mg2+ is the prevalent ion in the ribosome 1st solvation shell. We stress that the use of proper stereochemical guidelines in combination or not with other identification techniques, such as those pertaining to the detection of transition metals, of some anions and of K+ anomalous signals, is critical for deflating the current Mg2+ bubble witnessed in many ribosome and other RNA structures. We also stress that for the identification of lighter ions such as Mg2+, Na+, …, for which no anomalous signals can be detected, stereochemistry coupled with high resolution structures (<2.4 Å) remain the best currently available option.
Collapse
|
18
|
Kognole AA, MacKerell AD. Contributions and competition of Mg 2+ and K + in folding and stabilization of the Twister ribozyme. RNA (NEW YORK, N.Y.) 2020; 26:1704-1715. [PMID: 32769092 PMCID: PMC7566569 DOI: 10.1261/rna.076851.120] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/15/2020] [Accepted: 08/02/2020] [Indexed: 06/11/2023]
Abstract
Native folded and compact intermediate states of RNA typically involve tertiary structures in the presence of divalent ions such as Mg2+ in a background of monovalent ions. In a recent study, we have shown how the presence of Mg2+ impacts the transition from partially unfolded to folded states through a "push-pull" mechanism where the ion both favors and disfavors the sampling of specific phosphate-phosphate interactions. To further understand the ion atmosphere of RNA in folded and partially folded states results from atomistic umbrella sampling and oscillating chemical potential grand canonical Monte Carlo/molecular dynamics (GCMC/MD) simulations are used to obtain atomic-level details of the distributions of Mg2+ and K+ ions around Twister RNA. Results show the presence of 100 mM Mg2+ to lead to increased charge neutralization over that predicted by counterion condensation theory. Upon going from partially unfolded to folded states, overall charge neutralization increases at all studied ion concentrations that, while associated with an increase in the number of direct ion-phosphate interactions, is fully accounted for by the monovalent K+ ions. Furthermore, K+ preferentially interacts with purine N7 atoms of helical regions in partially unfolded states, thereby potentially stabilizing the helical regions. Thus, both secondary helical structures and formation of tertiary structures leads to increased counterion condensation, thereby stabilizing those structural features of Twister. Notably, it is shown that K+ can act as a surrogate for Mg2+ by participating in specific interactions with nonsequential phosphate pairs that occur in the folded state, explaining the ability of Twister to self-cleave at submillimolar Mg2+ concentrations.
Collapse
Affiliation(s)
- Abhishek A Kognole
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, Maryland 21201, USA
| | - Alexander D MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, Maryland 21201, USA
| |
Collapse
|
19
|
Cruz-León S, Schwierz N. Hofmeister Series for Metal-Cation-RNA Interactions: The Interplay of Binding Affinity and Exchange Kinetics. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2020; 36:5979-5989. [PMID: 32366101 PMCID: PMC7304902 DOI: 10.1021/acs.langmuir.0c00851] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/04/2023]
Abstract
A large variety of physicochemical properties involving RNA depends on the type of metal cation present in solution. In order to gain microscopic insight into the origin of these ion specific effects, we apply molecular dynamics simulations to describe the interactions of metal cations and RNA. For the three most common ion binding sites on RNA, we calculate the binding affinities and exchange rates of eight different mono- and divalent metal cations. Our results reveal that binding sites involving phosphate groups preferentially bind metal cations with high charge density (such as Mg2+) in inner-sphere conformations while binding sites involving N7 or O6 atoms preferentially bind cations with low charge density (such as K+). The binding affinity therefore follows a direct Hofmeister series at the backbone but is reversed at the nucleobases leading to a high selectivity of ion binding sites on RNA. In addition, the exchange rates for cation binding cover almost 5 orders of magnitude, leading to a vastly different time scale for the lifetimes of contact pairs. Taken together, the site-specific binding affinities and the specific lifetime of contact pairs provide the microscopic explanation of ion specific effects observed in a wide variety of macroscopic RNA properties. Finally, combining the results from atomistic simulations with extended Poisson-Boltzmann theory allows us to predict the distribution of metal cations around double-stranded RNA at finite concentrations and to reproduce the results of ion counting experiments with good accuracy.
Collapse
|
20
|
|
21
|
Importance of potassium ions for ribosome structure and function revealed by long-wavelength X-ray diffraction. Nat Commun 2019; 10:2519. [PMID: 31175275 PMCID: PMC6555806 DOI: 10.1038/s41467-019-10409-4] [Citation(s) in RCA: 103] [Impact Index Per Article: 20.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2018] [Accepted: 05/06/2019] [Indexed: 11/08/2022] Open
Abstract
The ribosome, the largest RNA-containing macromolecular machinery in cells, requires metal ions not only to maintain its three-dimensional fold but also to perform protein synthesis. Despite the vast biochemical data regarding the importance of metal ions for efficient protein synthesis and the increasing number of ribosome structures solved by X-ray crystallography or cryo-electron microscopy, the assignment of metal ions within the ribosome remains elusive due to methodological limitations. Here we present extensive experimental data on the potassium composition and environment in two structures of functional ribosome complexes obtained by measurement of the potassium anomalous signal at the K-edge, derived from long-wavelength X-ray diffraction data. We elucidate the role of potassium ions in protein synthesis at the three-dimensional level, most notably, in the environment of the ribosome functional decoding and peptidyl transferase centers. Our data expand the fundamental knowledge of the mechanism of ribosome function and structural integrity.
Collapse
|
22
|
Westhof E, Yusupov M, Yusupova G. The multiple flavors of GoU pairs in RNA. J Mol Recognit 2019; 32:e2782. [PMID: 31033092 PMCID: PMC6617799 DOI: 10.1002/jmr.2782] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2019] [Revised: 03/02/2019] [Accepted: 03/14/2019] [Indexed: 11/10/2022]
Abstract
Wobble GU pairs (or GoU) occur frequently within double‐stranded RNA helices interspersed within the standard G═C and A─U Watson‐Crick pairs. However, other types of GoU pairs interacting on their Watson‐Crick edges have been observed. The structural and functional roles of such alternative GoU pairs are surprisingly diverse and reflect the various pairings G and U can form by exploiting all the subtleties of their electronic configurations. Here, the structural characteristics of the GoU pairs are updated following the recent crystallographic structures of functional ribosomal complexes and the development in our understanding of ribosomal translation.
Collapse
Affiliation(s)
- Eric Westhof
- Architecture et Réactivité de l'ARN, Institut de biologie moléculaire et cellulaire du CNRS, Université de Strasbourg, Strasbourg, France
| | - Marat Yusupov
- Department of Integrated Structural Biology, Institute of Genetics and Molecular and Cellular Biology, INSERM, U964, CNRS, UMR7104, Université de Strasbourg, Illkirch, France
| | - Gulnara Yusupova
- Department of Integrated Structural Biology, Institute of Genetics and Molecular and Cellular Biology, INSERM, U964, CNRS, UMR7104, Université de Strasbourg, Illkirch, France
| |
Collapse
|
23
|
Cassone G, Kruse H, Sponer J. Interactions between cyclic nucleotides and common cations: an ab initio molecular dynamics study. Phys Chem Chem Phys 2019; 21:8121-8132. [PMID: 30932112 DOI: 10.1039/c8cp07492e] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
We present the first, to the best of our knowledge, ab initio molecular dynamics (AIMD) investigation on three aqueous solutions where an abasic cyclic nucleotide model is solvated in the presence of distinct cations (i.e., Na+, K+ and Mg2+). We elucidate the typical modalities of interaction between those ionic species and the nucleotide moiety by first-principles numerical simulations, starting from an inner-shell binding configuration on a time scale of 100 ps (total simulation time of ∼600 ps). Whereas the strong "structure-maker" Mg2+ is permanently bound to one of the two oxygen atoms of the phosphate group of the nucleotide model, Na+ and K+ show binding times τb of 65 ps and 10-15 ps, respectively, thus reflecting their chemical nature in aqueous solutions. Furthermore, we qualitatively relate these findings to approximate free-energy barriers of the cations' unbinding obtained by means of exploratory well-tempered metadynamics. With the aim of shedding light on the features of commonly employed force-fields (FFs), classical MD simulations (almost 200 trajectories with a total simulation time of ∼18 μs) using the biomolecular AMBER FF are also reported. By choosing several combinations of the parametrization for the water environment (i.e., TIP3P, SPC/E and OPC) and cations (i.e., Joung-Cheatham, Li-Merz 12-6 and Li-Merz 12-6-4), we found significant differences in the radial distribution functions and residence times compared to the ab initio results. The Na+ and K+ ions wrongly show quasi-identical radial distribution functions and the Li & Merz 12-6-4 Lennard-Jones parameters for Mg2+ were found to be essential in quickly reaching the binding state consistent with AIMD.
Collapse
Affiliation(s)
- Giuseppe Cassone
- Institute of Biophysics of the Czech Academy of Sciences, v.v.i, Královopolská 135, 61265 Brno, Czech Republic.
| | | | | |
Collapse
|