1
|
Das PK, Adil O, DeGregorio AP, Sumita M, Shamsi MH. Pseudouridine-modified RNA probe for label-free electrochemical detection of nucleic acids on 2D MoS 2 nanosheets. Analyst 2024; 149:1310-1317. [PMID: 38247383 DOI: 10.1039/d3an01832f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2024]
Abstract
RNA modification, particularly pseudouridine (Ψ), has played an important role in the development of the mRNA-based COVID-19 vaccine. This is because Ψ enhances RNA stability against nuclease activity and decreases the anti-RNA immune response. Ψ also provides structural flexibility to RNA by enhancing base stacking compared with canonical nucleobases. In this report, we demonstrate the first application of pseudouridine-modified RNA as a probe (Ψ-RNA) for label-free nucleic acid biosensing. It is known that MoS2 has a differential affinity for nucleic acids, which may be translated into a unique electronic signal. Herein, the Ψ-RNA probe interacts with the pristine MoS2 surface and causes a change in interfacial electrochemical charge transfer in the MoS2 nanosheets. Compared with an unmodified RNA probe, Ψ-RNA exhibited faster adsorption and higher affinity for MoS2. Moreover, Ψ-RNA could bind to complementary RNA and DNA targets with almost equal affinity when engaged with the MoS2 surface. Ψ-RNA maintained robust interactions with the MoS2 surface following the hybridization event, perhaps through its extra amino group. The detection sensitivity of the Ψ-RNA/MoS2 platform was as low as 500 attomoles, while the results also indicate that the probe can distinguish between complementary targets, single mismatches, and non-complementary nucleic acid sequences with statistical significance. This proof-of-concept study shows that the Ψ-RNA probe may solve numerous problems of adsorption-based biosensing platforms due to its stability and structural flexibility.
Collapse
Affiliation(s)
- Prabhangshu Kumer Das
- School of Chemical and Biomolecular Sciences, Southern Illinois University, Carbondale, IL 62901, USA.
| | - Omair Adil
- School of Chemical and Biomolecular Sciences, Southern Illinois University, Carbondale, IL 62901, USA.
| | - Anthony P DeGregorio
- Department of Chemistry, Southern Illinois University Edwardsville, IL, 62026, USA
| | - Minako Sumita
- Department of Chemistry, Southern Illinois University Edwardsville, IL, 62026, USA
| | - Mohtashim Hassan Shamsi
- School of Chemical and Biomolecular Sciences, Southern Illinois University, Carbondale, IL 62901, USA.
| |
Collapse
|
2
|
The Effects of Flexibility on dsDNA–dsDNA Interactions. Life (Basel) 2022; 12:life12050699. [PMID: 35629366 PMCID: PMC9147707 DOI: 10.3390/life12050699] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2022] [Revised: 05/02/2022] [Accepted: 05/06/2022] [Indexed: 11/29/2022] Open
Abstract
A detailed understanding of the physical mechanism of ion-mediated dsDNA interactions is important in biological functions such as DNA packaging and homologous pairing. We report the potential of mean force (PMF) or the effective solvent mediated interactions between two parallel identical dsDNAs as a function of interhelical separation in 0.15 M NaCl solution. Here, we study the influence of flexibility of dsDNAs on the effective interactions by comparing PMFs between rigid models and flexible ones. The role of flexibility of dsDNA pairs in their association is elucidated by studying the energetic properties of Na+ ions as well as the fluctuations of ions around dsDNAs. The introduction of flexibility of dsDNAs softens the vdW contact wall and induces more counterion fluctuations around dsDNAs. In addition, flexibility facilitates the Na+ ions dynamics affecting their distribution. The results quantify the extent of attraction influenced by dsDNA flexibility and further emphasize the importance of non-continuum solvation approaches.
Collapse
|
3
|
Wang Y, Guan X, Zhang S, Liu Y, Wang S, Fan P, Du X, Yan S, Zhang P, Chen HY, Li W, Zhang D, Huang S. Structural-profiling of low molecular weight RNAs by nanopore trapping/translocation using Mycobacterium smegmatis porin A. Nat Commun 2021; 12:3368. [PMID: 34099723 PMCID: PMC8185011 DOI: 10.1038/s41467-021-23764-y] [Citation(s) in RCA: 35] [Impact Index Per Article: 11.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2020] [Accepted: 05/06/2021] [Indexed: 12/15/2022] Open
Abstract
Folding of RNA can produce elaborate tertiary structures, corresponding to their diverse roles in the regulation of biological activities. Direct observation of RNA structures at high resolution in their native form however remains a challenge. The large vestibule and the narrow constriction of a Mycobacterium smegmatis porin A (MspA) suggests a sensing mode called nanopore trapping/translocation, which clearly distinguishes between microRNA, small interfering RNA (siRNA), transfer RNA (tRNA) and 5 S ribosomal RNA (rRNA). To further profit from the acquired event characteristics, a custom machine learning algorithm is developed. Events from measurements with a mixture of RNA analytes can be automatically classified, reporting a general accuracy of ~93.4%. tRNAs, which possess a unique tertiary structure, report a highly distinguishable sensing feature, different from all other RNA types tested in this study. With this strategy, tRNAs from different sources are measured and a high structural conservation across different species is observed in single molecule.
Collapse
MESH Headings
- Machine Learning
- MicroRNAs/chemistry
- MicroRNAs/genetics
- MicroRNAs/metabolism
- Molecular Dynamics Simulation
- Molecular Weight
- Mycobacterium smegmatis/genetics
- Mycobacterium smegmatis/metabolism
- Nanopores
- Nucleic Acid Conformation
- Porins/chemistry
- Porins/genetics
- Porins/metabolism
- RNA/chemistry
- RNA/genetics
- RNA/metabolism
- RNA Folding
- RNA Transport
- RNA, Ribosomal, 5S/chemistry
- RNA, Ribosomal, 5S/genetics
- RNA, Ribosomal, 5S/metabolism
- RNA, Small Interfering/chemistry
- RNA, Small Interfering/genetics
- RNA, Small Interfering/metabolism
- RNA, Transfer/chemistry
- RNA, Transfer/genetics
- RNA, Transfer/metabolism
Collapse
Affiliation(s)
- Yuqin Wang
- State Key Laboratory of Analytical Chemistry for Life Sciences, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China
- Chemistry and Biomedicine Innovation Center (ChemBIC), Nanjing University, Nanjing, China
| | - Xiaoyu Guan
- College of Computer Science and Technology, Nanjing University of Aeronautics and Astronautics, MIIT Key Laboratory of Pattern Analysis and Machine Intelligence, Nanjing, China
| | - Shanyu Zhang
- State Key Laboratory of Analytical Chemistry for Life Sciences, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China
- Chemistry and Biomedicine Innovation Center (ChemBIC), Nanjing University, Nanjing, China
| | - Yao Liu
- State Key Laboratory of Analytical Chemistry for Life Sciences, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China
- Chemistry and Biomedicine Innovation Center (ChemBIC), Nanjing University, Nanjing, China
| | - Sha Wang
- State Key Laboratory of Analytical Chemistry for Life Sciences, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China
- Chemistry and Biomedicine Innovation Center (ChemBIC), Nanjing University, Nanjing, China
| | - Pingping Fan
- State Key Laboratory of Analytical Chemistry for Life Sciences, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China
- Chemistry and Biomedicine Innovation Center (ChemBIC), Nanjing University, Nanjing, China
| | - Xiaoyu Du
- State Key Laboratory of Analytical Chemistry for Life Sciences, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China
- Chemistry and Biomedicine Innovation Center (ChemBIC), Nanjing University, Nanjing, China
| | - Shuanghong Yan
- State Key Laboratory of Analytical Chemistry for Life Sciences, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China
- Chemistry and Biomedicine Innovation Center (ChemBIC), Nanjing University, Nanjing, China
| | - Panke Zhang
- State Key Laboratory of Analytical Chemistry for Life Sciences, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China
| | - Hong-Yuan Chen
- State Key Laboratory of Analytical Chemistry for Life Sciences, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China
| | - Wenfei Li
- Collaborative Innovation Center of Advanced Microstructures, National Laboratory of Solid State Microstructure, Department of Physics, Nanjing University, Nanjing, China
| | - Daoqiang Zhang
- College of Computer Science and Technology, Nanjing University of Aeronautics and Astronautics, MIIT Key Laboratory of Pattern Analysis and Machine Intelligence, Nanjing, China
| | - Shuo Huang
- State Key Laboratory of Analytical Chemistry for Life Sciences, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China
- Chemistry and Biomedicine Innovation Center (ChemBIC), Nanjing University, Nanjing, China
| |
Collapse
|
4
|
Liu T, Yu T, Zhang S, Wang Y, Zhang W. Thermodynamic and kinetic properties of a single base pair in A-DNA and B-DNA. Phys Rev E 2021; 103:042409. [PMID: 34005973 DOI: 10.1103/physreve.103.042409] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2020] [Accepted: 03/27/2021] [Indexed: 11/07/2022]
Abstract
Double stranded DNA can adopt different forms, the so-called A-, B-, and Z-DNA, which play different biological roles. In this work, the thermodynamic and the kinetic parameters for the base-pair closing and opening in A-DNA and B-DNA were calculated by all-atom molecular dynamics simulations at different temperatures. The thermodynamic parameters of the base pair in B-DNA were in good agreement with the experimental results. The free energy barrier of breaking a single base stack results from the enthalpy increase ΔH caused by the disruption of hydrogen bonding and base-stacking interactions, as well as water and base interactions. The free energy barrier of base pair closing comes from the unfavorable entropy loss ΔS caused by the restriction of torsional angles and hydration. It was found that the enthalpy change ΔH and the entropy change ΔS for the base pair in A-DNA are much larger than those in B-DNA, and the transition rates between the opening and the closing state for the base pair in A-DNA are much slower than those in B-DNA. The large difference of the enthalpy and entropy change for forming the base pair in A-DNA and B-DNA results from different hydration in A-DNA and B-DNA. The hydration pattern observed around DNA is an accompanying process for forming the base pair, rather than a follow-up of the conformation.
Collapse
Affiliation(s)
- Taigang Liu
- Department of Physics Wuhan University, Wuhan 430072, China.,School of Medical Engineering, Xinxiang Medical University, Xinxiang 453003, China
| | - Ting Yu
- Department of Physics Wuhan University, Wuhan 430072, China
| | - Shuhao Zhang
- Department of Physics Wuhan University, Wuhan 430072, China
| | - Yujie Wang
- Department of Physics Wuhan University, Wuhan 430072, China.,Department of Physics and Telecommunication Engineering, Zhoukou Normal University, Zhoukou 466000, China
| | - Wenbing Zhang
- Department of Physics Wuhan University, Wuhan 430072, China
| |
Collapse
|
5
|
Wang Y, Wang Z, Wang Y, Liu T, Zhang W. The nearest neighbor and next nearest neighbor effects on the thermodynamic and kinetic properties of RNA base pair. J Chem Phys 2018; 148:045101. [PMID: 29390847 DOI: 10.1063/1.5013282] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The thermodynamic and kinetic parameters of an RNA base pair with different nearest and next nearest neighbors were obtained through long-time molecular dynamics simulation of the opening-closing switch process of the base pair near its melting temperature. The results indicate that thermodynamic parameters of GC base pair are dependent on the nearest neighbor base pair, and the next nearest neighbor base pair has little effect, which validated the nearest-neighbor model. The closing and opening rates of the GC base pair also showed nearest neighbor dependences. At certain temperature, the closing and opening rates of the GC pair with nearest neighbor AU is larger than that with the nearest neighbor GC, and the next nearest neighbor plays little role. The free energy landscape of the GC base pair with the nearest neighbor GC is rougher than that with nearest neighbor AU.
Collapse
Affiliation(s)
- Yujie Wang
- Department of Physics, Wuhan University, Wuhan, Hubei, People's Republic of China
| | - Zhen Wang
- Department of Physics, Wuhan University, Wuhan, Hubei, People's Republic of China
| | - Yanli Wang
- Department of Physics, Wuhan University, Wuhan, Hubei, People's Republic of China
| | - Taigang Liu
- Department of Physics, Wuhan University, Wuhan, Hubei, People's Republic of China
| | - Wenbing Zhang
- Department of Physics, Wuhan University, Wuhan, Hubei, People's Republic of China
| |
Collapse
|
6
|
Š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: 327] [Impact Index Per Article: 54.5] [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
|
7
|
Sun T, Mirzoev A, Korolev N, Lyubartsev AP, Nordenskiöld L. All-Atom MD Simulation of DNA Condensation Using Ab Initio Derived Force Field Parameters of Cobalt(III)-Hexammine. J Phys Chem B 2017; 121:7761-7770. [DOI: 10.1021/acs.jpcb.7b03793] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Affiliation(s)
- Tiedong Sun
- School
of Biological Sciences, Nanyang Technological University, Singapore 637551
| | - Alexander Mirzoev
- School
of Biological Sciences, Nanyang Technological University, Singapore 637551
| | - Nikolay Korolev
- School
of Biological Sciences, Nanyang Technological University, Singapore 637551
| | - Alexander P. Lyubartsev
- Department
of Materials and Environmental Chemistry, Stockholm University, 10691 Stockholm, Sweden
| | - Lars Nordenskiöld
- School
of Biological Sciences, Nanyang Technological University, Singapore 637551
| |
Collapse
|
8
|
Wang Y, Gong S, Wang Z, Zhang W. The thermodynamics and kinetics of a nucleotide base pair. J Chem Phys 2017; 144:115101. [PMID: 27004898 DOI: 10.1063/1.4944067] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The thermodynamic and kinetic parameters of an RNA base pair were obtained through a long-time molecular dynamics simulation of the opening-closing switch process of the base pair near its melting temperature. The thermodynamic parameters were in good agreement with the nearest-neighbor model. The opening rates showed strong temperature dependence, however, the closing rates showed only weak temperature dependence. The transition path time was weakly temperature dependent and was insensitive to the energy barrier. The diffusion constant exhibited super-Arrhenius behavior. The free energy barrier of breaking a single base stack results from the enthalpy increase, ΔH, caused by the disruption of hydrogen bonding and base-stacking interactions. The free energy barrier of base pair closing comes from the unfavorable entropy loss, ΔS, caused by the restriction of torsional angles. These results suggest that a one-dimensional free energy surface is sufficient to accurately describe the dynamics of base pair opening and closing, and the dynamics are Brownian.
Collapse
Affiliation(s)
- Yujie Wang
- Department of Physics, Wuhan University, Wuhan, Hubei, People's Republic of China
| | - Sha Gong
- Department of Physics, Wuhan University, Wuhan, Hubei, People's Republic of China
| | - Zhen Wang
- Department of Physics, Wuhan University, Wuhan, Hubei, People's Republic of China
| | - Wenbing Zhang
- Department of Physics, Wuhan University, Wuhan, Hubei, People's Republic of China
| |
Collapse
|
9
|
Vangaveti S, Ranganathan SV, Chen AA. Advances in RNA molecular dynamics: a simulator's guide to RNA force fields. WILEY INTERDISCIPLINARY REVIEWS-RNA 2016; 8. [DOI: 10.1002/wrna.1396] [Citation(s) in RCA: 35] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/18/2016] [Revised: 08/17/2016] [Accepted: 08/25/2016] [Indexed: 01/04/2023]
Affiliation(s)
- Sweta Vangaveti
- The RNA Institute; University at Albany State University of New York; Albany NY USA
| | | | - Alan A. Chen
- The RNA Institute; University at Albany State University of New York; Albany NY USA
- Department of Chemistry; University at Albany State University of New York; Albany NY USA
| |
Collapse
|
10
|
Lee TS, Radak BK, Harris ME, York DM. A Two-Metal-Ion-Mediated Conformational Switching Pathway for HDV Ribozyme Activation. ACS Catal 2016; 6:1853-1869. [PMID: 27774349 DOI: 10.1021/acscatal.5b02158] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
RNA enzymes serve as a potentially powerful platform from which to design catalysts and engineer new biotechnology. A fundamental understanding of these systems provides insight to guide design. The hepatitis delta virus ribozyme (HDVr) is a small, self-cleaving RNA motif widely distributed in nature, that has served as a paradigm for understanding basic principles of RNA catalysis. Nevertheless, questions remain regarding the precise roles of divalent metal ions and key nucleotides in catalysis. In an effort to establish a reaction mechanism model consistent with available experimental data, we utilize molecular dynamics simulations to explore different conformations and metal ion binding modes along the HDVr reaction path. Building upon recent crystallographic data, our results provide a dynamic model of the HDVr reaction mechanism involving a conformational switch between multiple non-canonical G25:U20 base pair conformations in the active site. These local nucleobase dynamics play an important role in catalysis by modulating the metal binding environments of two Mg2+ ions that support catalysis at different steps of the reaction pathway. The first ion plays a structural role by inducing a base pair flip necessary to obtain the catalytic fold in which C75 moves towards to the scissile phosphate in the active site. Ejection of this ion then permits a second ion to bind elsewhere in the active site and facilitate nucleophile activation. The simulations collectively describe a mechanistic scenario that is consistent with currently available experimental data from crystallography, phosphorothioate substitutions, and chemical probing studies. Avenues for further experimental verification are suggested.
Collapse
Affiliation(s)
- Tai-Sung Lee
- Center for Integrative Proteomics Research and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, United States
| | - Brian K. Radak
- Center for Integrative Proteomics Research and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, United States
- Argonne National Laboratory, Argonne, Illinois 60439, United State
| | - Michael E. Harris
- Department of Biochemistry, Case Western Reserve University School of Medicine, Cleveland, Ohio 44106, United States
| | - Darrin M. York
- Center for Integrative Proteomics Research and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, United States
| |
Collapse
|
11
|
Anisimov VM, Lamoureux G, Vorobyov IV, Huang N, Roux B, MacKerell AD. Determination of Electrostatic Parameters for a Polarizable Force Field Based on the Classical Drude Oscillator. J Chem Theory Comput 2015; 1:153-68. [PMID: 26641126 DOI: 10.1021/ct049930p] [Citation(s) in RCA: 235] [Impact Index Per Article: 26.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A procedure to determine the electrostatic parameters has been developed for a polarizable empirical force field based on the classical Drude oscillator model. Atomic charges and polarizabilities for a given molecule of interest were derived from restrained fitting to quantum-mechanical electrostatic potentials (ESP) calculated at the B3LYP/ cc-pVDZ or B3LYP/aug-cc-pVDZ levels on grid points located on concentric Connolly surfaces. The determination of the atomic polarizabilities requires a series of perturbed ESP maps, each one representing the electronic response of the molecule in the presence of a background charge placed on Connolly surfaces primarily along chemical bonds and lone pairs. Reference values for the partial atomic charges were taken from the CHARMM27 additive all-atom force field, and those for the polarizabilities were based on adjusted Miller's ahp atomic polarizability values. The fitted values of atomic polarizabilities were scaled to reflect the reduced polarization expected for the condensed media and/or to correct for the systematic underestimation of experimental molecular polarizabilities by B3LYP calculations. Following correction of the polarizabilities, the atomic charges were adjusted to reproduce gas-phase dipole moments. The developed scheme has been tested on a set of small molecules representing functional moieties of nucleic acids. The derived electrostatic parameters have been successfully applied in a preliminary polarizable molecular dynamics simulation of a DNA octamer in a box of water with sodium counterions. Thus, this study confirms the feasibility of the use of a polarizable force field based on a classical Drude model for simulations of biomolecules in the condensed phase.
Collapse
Affiliation(s)
- Victor M Anisimov
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, Maryland, Département de physique, Université de Montréal, C.P. 6128, succ. centre-ville, Montréal, Québec, Canada H3C 3J7, and Department of Biochemistry, Weill Medical College of Cornell University, New York, New York
| | - Guillaume Lamoureux
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, Maryland, Département de physique, Université de Montréal, C.P. 6128, succ. centre-ville, Montréal, Québec, Canada H3C 3J7, and Department of Biochemistry, Weill Medical College of Cornell University, New York, New York
| | - Igor V Vorobyov
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, Maryland, Département de physique, Université de Montréal, C.P. 6128, succ. centre-ville, Montréal, Québec, Canada H3C 3J7, and Department of Biochemistry, Weill Medical College of Cornell University, New York, New York
| | - Niu Huang
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, Maryland, Département de physique, Université de Montréal, C.P. 6128, succ. centre-ville, Montréal, Québec, Canada H3C 3J7, and Department of Biochemistry, Weill Medical College of Cornell University, New York, New York
| | - Benoît Roux
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, Maryland, Département de physique, Université de Montréal, C.P. 6128, succ. centre-ville, Montréal, Québec, Canada H3C 3J7, and Department of Biochemistry, Weill Medical College of Cornell University, New York, New York
| | - Alexander D MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, Maryland, Département de physique, Université de Montréal, C.P. 6128, succ. centre-ville, Montréal, Québec, Canada H3C 3J7, and Department of Biochemistry, Weill Medical College of Cornell University, New York, New York
| |
Collapse
|
12
|
Vanommeslaeghe K, MacKerell AD. CHARMM additive and polarizable force fields for biophysics and computer-aided drug design. BIOCHIMICA ET BIOPHYSICA ACTA 2015; 1850:861-871. [PMID: 25149274 PMCID: PMC4334745 DOI: 10.1016/j.bbagen.2014.08.004] [Citation(s) in RCA: 174] [Impact Index Per Article: 19.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/07/2014] [Revised: 08/10/2014] [Accepted: 08/12/2014] [Indexed: 11/29/2022]
Abstract
BACKGROUND Molecular Mechanics (MM) is the method of choice for computational studies of biomolecular systems owing to its modest computational cost, which makes it possible to routinely perform molecular dynamics (MD) simulations on chemical systems of biophysical and biomedical relevance. SCOPE OF REVIEW As one of the main factors limiting the accuracy of MD results is the empirical force field used, the present paper offers a review of recent developments in the CHARMM additive force field, one of the most popular biomolecular force fields. Additionally, we present a detailed discussion of the CHARMM Drude polarizable force field, anticipating a growth in the importance and utilization of polarizable force fields in the near future. Throughout the discussion emphasis is placed on the force fields' parametrization philosophy and methodology. MAJOR CONCLUSIONS Recent improvements in the CHARMM additive force field are mostly related to newly found weaknesses in the previous generation of additive force fields. Beyond the additive approximation is the newly available CHARMM Drude polarizable force field, which allows for MD simulations of up to 1μs on proteins, DNA, lipids and carbohydrates. GENERAL SIGNIFICANCE Addressing the limitations ensures the reliability of the new CHARMM36 additive force field for the types of calculations that are presently coming into routine computational reach while the availability of the Drude polarizable force fields offers an inherently more accurate model of the underlying physical forces driving macromolecular structures and dynamics. This article is part of a Special Issue entitled "Recent developments of molecular dynamics".
Collapse
Affiliation(s)
- K Vanommeslaeghe
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, MD 21201, USA
| | - A D MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, MD 21201, USA.
| |
Collapse
|
13
|
Galindo-Murillo R, Roe DR, Cheatham TE. Convergence and reproducibility in molecular dynamics simulations of the DNA duplex d(GCACGAACGAACGAACGC). Biochim Biophys Acta Gen Subj 2014; 1850:1041-1058. [PMID: 25219455 DOI: 10.1016/j.bbagen.2014.09.007] [Citation(s) in RCA: 119] [Impact Index Per Article: 11.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2014] [Revised: 09/01/2014] [Accepted: 09/03/2014] [Indexed: 12/12/2022]
Abstract
BACKGROUND The structure and dynamics of DNA are critically related to its function. Molecular dynamics simulations augment experiment by providing detailed information about the atomic motions. However, to date the simulations have not been long enough for convergence of the dynamics and structural properties of DNA. METHODS Molecular dynamics simulations performed with AMBER using the ff99SB force field with the parmbsc0 modifications, including ensembles of independent simulations, were compared to long timescale molecular dynamics performed with the specialized Anton MD engine on the B-DNA structure d(GCACGAACGAACGAACGC). To assess convergence, the decay of the average RMSD values over longer and longer time intervals was evaluated in addition to assessing convergence of the dynamics via the Kullback-Leibler divergence of principal component projection histograms. RESULTS These molecular dynamics simulations-including one of the longest simulations of DNA published to date at ~44μs-surprisingly suggest that the structure and dynamics of the DNA helix, neglecting the terminal base pairs, are essentially fully converged on the ~1-5μs timescale. CONCLUSIONS We can now reproducibly converge the structure and dynamics of B-DNA helices, omitting the terminal base pairs, on the μs time scale with both the AMBER and CHARMM C36 nucleic acid force fields. Results from independent ensembles of simulations starting from different initial conditions, when aggregated, match the results from long timescale simulations on the specialized Anton MD engine. GENERAL SIGNIFICANCE With access to large-scale GPU resources or the specialized MD engine "Anton" it is possible for a variety of molecular systems to reproducibly and reliably converge the conformational ensemble of sampled structures. This article is part of a Special Issue entitled: Recent developments of molecular dynamics.
Collapse
Affiliation(s)
- Rodrigo Galindo-Murillo
- Department of Medicinal Chemistry, L.S. Skaggs Pharmacy Institute, University of Utah, Salt Lake City, UT 84112, USA
| | - Daniel R Roe
- Department of Medicinal Chemistry, L.S. Skaggs Pharmacy Institute, University of Utah, Salt Lake City, UT 84112, USA
| | - Thomas E Cheatham
- Department of Medicinal Chemistry, L.S. Skaggs Pharmacy Institute, University of Utah, Salt Lake City, UT 84112, USA.
| |
Collapse
|
14
|
Mukherjee S, Bhattacharyya D. Influence of divalent magnesium ion on DNA: molecular dynamics simulation studies. J Biomol Struct Dyn 2013; 31:896-912. [DOI: 10.1080/07391102.2012.713780] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
|
15
|
Weber G. Mesoscopic model parametrization of hydrogen bonds and stacking interactions of RNA from melting temperatures. Nucleic Acids Res 2012; 41:e30. [PMID: 23087379 PMCID: PMC3592459 DOI: 10.1093/nar/gks964] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/05/2022] Open
Abstract
Information about molecular interactions in DNA can be obtained from experimental melting temperature data by using mesoscopic statistical physics models. Here, we extend the technique to RNA and show that the new parameters correctly reproduce known properties such as the stronger hydrogen bonds of AU base pairs. We also were able to calculate a complete set of elastic constants for all 10 irreducible combinations of nearest neighbours (NNs). We believe that this is particularly useful as experimentally derived information about RNA elasticity is relatively scarce. The melting temperature prediction using the present model improves over those from traditional NN model, providing thus an alternative way to calculate these temperatures for RNA. Additionally, we calculated the site-dependent base pair oscillation to explain why RNA shows larger oscillation amplitudes despite having stronger AU hydrogen bonds.
Collapse
Affiliation(s)
- Gerald Weber
- Departamento de Física, Universidade Federal de Minas Gerais, 31270-901 Belo Horizonte-MG, Brazil.
| |
Collapse
|
16
|
Halder S, Bhattacharyya D. Structural Variations of Single and Tandem Mismatches in RNA Duplexes: A Joint MD Simulation and Crystal Structure Database Analysis. J Phys Chem B 2012; 116:11845-56. [PMID: 22953716 DOI: 10.1021/jp305628v] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Sukanya Halder
- Biophysics
Division and ‡Computational Science Division, Saha Institute of Nuclear Physics, Kolkata, West Bengal, 700 064, India
| | - Dhananjay Bhattacharyya
- Biophysics
Division and ‡Computational Science Division, Saha Institute of Nuclear Physics, Kolkata, West Bengal, 700 064, India
| |
Collapse
|
17
|
Beššeová I, Banáš P, Kührová P, Košinová P, Otyepka M, Šponer J. Simulations of A-RNA Duplexes. The Effect of Sequence, Solute Force Field, Water Model, and Salt Concentration. J Phys Chem B 2012; 116:9899-916. [DOI: 10.1021/jp3014817] [Citation(s) in RCA: 59] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Ivana Beššeová
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska
135, 612 65 Brno, Czech Republic
| | - Pavel Banáš
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska
135, 612 65 Brno, Czech Republic
- Regional Centre of Advanced
Technologies and Materials, Department of Physical Chemistry, Faculty
of Science, Palacky University, tr. 17
listopadu 12, 771 46, Olomouc, Czech Republic
| | - Petra Kührová
- Regional Centre of Advanced
Technologies and Materials, Department of Physical Chemistry, Faculty
of Science, Palacky University, tr. 17
listopadu 12, 771 46, Olomouc, Czech Republic
| | - Pavlína Košinová
- Regional Centre of Advanced
Technologies and Materials, Department of Physical Chemistry, Faculty
of Science, Palacky University, tr. 17
listopadu 12, 771 46, Olomouc, Czech Republic
| | - Michal Otyepka
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska
135, 612 65 Brno, Czech Republic
- Regional Centre of Advanced
Technologies and Materials, Department of Physical Chemistry, Faculty
of Science, Palacky University, tr. 17
listopadu 12, 771 46, Olomouc, Czech Republic
| | - Jiří Šponer
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska
135, 612 65 Brno, Czech Republic
- CEITEC - Central European Institute of Technology, Campus Bohunice, Kamenice
5, 625 00 Brno, Czech Republic
| |
Collapse
|
18
|
Nakanishi K, Weinberg DE, Bartel DP, Patel DJ. Structure of yeast Argonaute with guide RNA. Nature 2012; 486:368-74. [PMID: 22722195 PMCID: PMC3853139 DOI: 10.1038/nature11211] [Citation(s) in RCA: 245] [Impact Index Per Article: 20.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2012] [Accepted: 05/11/2012] [Indexed: 01/07/2023]
Abstract
The RNA-induced silencing complex, comprising Argonaute and guide RNA, mediates RNA interference. Here we report the 3.2 Å crystal structure of Kluyveromyces polysporus Argonaute (KpAGO) fortuitously complexed with guide RNA originating from small-RNA duplexes autonomously loaded by recombinant KpAGO. Despite their diverse sequences, guide-RNA nucleotides 1-8 are positioned similarly, with sequence-independent contacts to bases, phosphates and 2'-hydroxyl groups pre-organizing the backbone of nucleotides 2-8 in a near-A-form conformation. Compared with prokaryotic Argonautes, KpAGO has numerous surface-exposed insertion segments, with a cluster of conserved insertions repositioning the N domain to enable full propagation of guide-target pairing. Compared with Argonautes in inactive conformations, KpAGO has a hydrogen-bond network that stabilizes an expanded and repositioned loop, which inserts an invariant glutamate into the catalytic pocket. Mutation analyses and analogies to ribonuclease H indicate that insertion of this glutamate finger completes a universally conserved catalytic tetrad, thereby activating Argonaute for RNA cleavage.
Collapse
Affiliation(s)
- Kotaro Nakanishi
- Structural Biology Program, Memorial Sloan-Kettering Cancer Center, New York, NY 10065, USA
| | - David E. Weinberg
- Whitehead Institute for Biomedical Research, 9 Cambridge Center, Cambridge, MA 02142, USA
- Howard Hughes Medical Institute and Department of Biology, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | - David P. Bartel
- Whitehead Institute for Biomedical Research, 9 Cambridge Center, Cambridge, MA 02142, USA
- Howard Hughes Medical Institute and Department of Biology, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | - Dinshaw J. Patel
- Structural Biology Program, Memorial Sloan-Kettering Cancer Center, New York, NY 10065, USA
| |
Collapse
|
19
|
Bereźniak T, Jäschke A, Smith JC, Imhof P. Stereoselection in the diels-alderase ribozyme: A molecular dynamics study. J Comput Chem 2012; 33:1603-14. [DOI: 10.1002/jcc.22993] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2011] [Revised: 03/05/2012] [Accepted: 03/18/2012] [Indexed: 01/03/2023]
|
20
|
Zhu X, Lopes PE, MacKerell AD. Recent Developments and Applications of the CHARMM force fields. WILEY INTERDISCIPLINARY REVIEWS. COMPUTATIONAL MOLECULAR SCIENCE 2012; 2:167-185. [PMID: 23066428 PMCID: PMC3468154 DOI: 10.1002/wcms.74] [Citation(s) in RCA: 127] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Empirical force fields commonly used to describe the condensed phase properties of complex systems such as biological macromolecules are continuously being updated. Improvements in quantum mechanical (QM) methods used to generate target data, availability of new experimental target data, incorporation of new classes of compounds and new theoretical developments (eg. polarizable methods) make force-field development a dynamic domain of research. Accordingly, a number of improvements and extensions of the CHARMM force fields have occurred over the years. The objective of the present review is to provide an up-to-date overview of the CHARMM force fields. A limited presentation on the historical aspects of force fields will be given, including underlying methodologies and principles, along with a brief description of the strategies used for parameter development. This is followed by information on the CHARMM additive and polarizable force fields, including examples of recent applications of those force fields.
Collapse
Affiliation(s)
- Xiao Zhu
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, MD 21201
| | - Pedro E.M. Lopes
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, MD 21201
| | - Alexander D. MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, MD 21201
| |
Collapse
|
21
|
Fulle S, Gohlke H. Flexibility analysis of biomacromolecules with application to computer-aided drug design. Methods Mol Biol 2012; 819:75-91. [PMID: 22183531 DOI: 10.1007/978-1-61779-465-0_6] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/31/2023]
Abstract
Flexibility characteristics of biomacromolecules can be efficiently determined down to the atomic level by a graph-theoretical technique as implemented in the FIRST (Floppy Inclusion and Rigid Substructure Topology) and ProFlex software packages. The method has been successfully applied to a series of protein and nucleic acid structures. Here, we describe practical guidelines for setting up and performing a flexibility analysis, discuss current bottlenecks of the approach, and provide sample applications as to how this technique can support computer-aided drug design approaches.
Collapse
Affiliation(s)
- Simone Fulle
- Institute of Pharmaceutical and Medicinal Chemistry, Department of Mathematics and Natural Sciences, Heinrich-Heine-University, Düsseldorf, Germany
| | | |
Collapse
|
22
|
Denning EJ, Priyakumar UD, Nilsson L, MacKerell AD. Impact of 2'-hydroxyl sampling on the conformational properties of RNA: update of the CHARMM all-atom additive force field for RNA. J Comput Chem 2011; 32:1929-43. [PMID: 21469161 PMCID: PMC3082605 DOI: 10.1002/jcc.21777] [Citation(s) in RCA: 288] [Impact Index Per Article: 22.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2010] [Revised: 01/24/2011] [Accepted: 01/30/2011] [Indexed: 01/02/2023]
Abstract
Here, we present an update of the CHARMM27 all-atom additive force field for nucleic acids that improves the treatment of RNA molecules. The original CHARMM27 force field parameters exhibit enhanced Watson-Crick base pair opening which is not consistent with experiment, whereas analysis of molecular dynamics (MD) simulations show the 2'-hydroxyl moiety to almost exclusively sample the O3' orientation. Quantum mechanical (QM) studies of RNA related model compounds indicate the energy minimum associated with the O3' orientation to be too favorable, consistent with the MD results. Optimization of the dihedral parameters dictating the energy of the 2'-hydroxyl proton targeting the QM data yielded several parameter sets, which sample both the base and O3' orientations of the 2'-hydroxyl to varying degrees. Selection of the final dihedral parameters was based on reproduction of hydration behavior as related to a survey of crystallographic data and better agreement with experimental NMR J-coupling values. Application of the model, designated CHARMM36, to a collection of canonical and noncanonical RNA molecules reveals overall improved agreement with a range of experimental observables as compared to CHARMM27. The results also indicate the sensitivity of the conformational heterogeneity of RNA to the orientation of the 2'-hydroxyl moiety and support a model whereby the 2'-hydroxyl can enhance the probability of conformational transitions in RNA.
Collapse
Affiliation(s)
- Elizabeth J. Denning
- Department of Pharmaceutical Sciences, School of Pharmacy, University
of Maryland, Baltimore, MD 21201
| | - U. Deva Priyakumar
- Department of Pharmaceutical Sciences, School of Pharmacy, University
of Maryland, Baltimore, MD 21201
| | - Lennart Nilsson
- Department of Pharmaceutical Sciences, School of Pharmacy, University
of Maryland, Baltimore, MD 21201
| | - Alexander D. MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University
of Maryland, Baltimore, MD 21201
| |
Collapse
|
23
|
Paciello G, Acquaviva A, Ficarra E, Deriu MA, Macii E. A molecular dynamics study of a miRNA:mRNA interaction. J Mol Model 2011; 17:2895-906. [DOI: 10.1007/s00894-011-0991-x] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2010] [Accepted: 01/24/2011] [Indexed: 12/25/2022]
|
24
|
Wienken CJ, Baaske P, Duhr S, Braun D. Thermophoretic melting curves quantify the conformation and stability of RNA and DNA. Nucleic Acids Res 2011; 39:e52. [PMID: 21297115 PMCID: PMC3082908 DOI: 10.1093/nar/gkr035] [Citation(s) in RCA: 48] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
Abstract
Measuring parameters such as stability and conformation of biomolecules, especially of nucleic acids, is important in the field of biology, medical diagnostics and biotechnology. We present a thermophoretic method to analyse the conformation and thermal stability of nucleic acids. It relies on the directed movement of molecules in a temperature gradient that depends on surface characteristics of the molecule, such as size, charge and hydrophobicity. By measuring thermophoresis of nucleic acids over temperature, we find clear melting transitions and resolve intermediate conformational states. These intermediate states are indicated by an additional peak in the thermophoretic signal preceding most melting transitions. We analysed single nucleotide polymorphisms, DNA modifications, conformational states of DNA hairpins and microRNA duplexes. The method is validated successfully against calculated melting temperatures and UV absorbance measurements. Interestingly, the methylation of DNA is detected by the thermophoretic amplitude even if it does not affect the melting temperature. In the described setup, thermophoresis is measured all-optical in a simple setup using a reproducible capillary format with only 250 nl probe consumption. The thermophoretic analysis of nucleic acids shows the technique’s versatility for the investigation of nucleic acids relevant in cellular processes like RNA interference or gene silencing.
Collapse
Affiliation(s)
- Christoph J Wienken
- Physics Department and Center for NanoScience, Ludwig-Maximilians-Universität München, Amalienstrasse 54, 80799 München, Germany
| | | | | | | |
Collapse
|
25
|
Halder S, Bhattacharyya D. Structural Stability of Tandemly Occurring Noncanonical Basepairs within Double Helical Fragments: Molecular Dynamics Studies of Functional RNA. J Phys Chem B 2010; 114:14028-40. [DOI: 10.1021/jp102835t] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Affiliation(s)
- Sukanya Halder
- Biophysics Division, Saha Institute of Nuclear Physics, 1/AF, Bidhannagar, Kolkata—700 064, India
| | - Dhananjay Bhattacharyya
- Biophysics Division, Saha Institute of Nuclear Physics, 1/AF, Bidhannagar, Kolkata—700 064, India
| |
Collapse
|
26
|
Taniguchi Y, Kurose Y, Nishioka T, Nagatsugi F, Sasaki S. The alkyl-connected 2-amino-6-vinylpurine (AVP) crosslinking agent for improved selectivity to the cytosine base in RNA. Bioorg Med Chem 2010; 18:2894-901. [PMID: 20346683 DOI: 10.1016/j.bmc.2010.03.008] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2010] [Revised: 03/03/2010] [Accepted: 03/04/2010] [Indexed: 11/25/2022]
Abstract
We have previously reported that the 2-amino-6-vinylpurine (AVP) nucleoside exhibits a highly efficient and selective crosslinking reaction toward cytosine and displayed an improved antisense inhibition in cultured cells. In this study, we further investigated the alkyl-connected AVP nucleoside analogs for more efficient crosslinking to the cytosine base (rC) of the target RNA. We synthesized three AVP analogs which connect the 2-amino-6-vinylpurine unit to the 2'-deoxyribose through a methylene, an ethylene, or a butylene linker. The ODN incorporating the AVP analog with the methylene or the butylene linker showed a slightly higher crosslinking to the target rC of RNA than the original AVP with no linker. In contrast, the AVP with the ethylene linker formed a selective and efficient crosslink to the rC of the target RNA.
Collapse
Affiliation(s)
- Yosuke Taniguchi
- Graduate School of Pharmaceutical Sciences, Kyushu University, 3-1-1 Maidashi, Higashi-ku, Fukuoka 812-8582, Japan
| | | | | | | | | |
Collapse
|
27
|
Wereszczynski J, Andricioaei I. Conformational and Solvent Entropy Contributions to the Thermal Response of Nucleic Acid-Based Nanothermometers. J Phys Chem B 2010; 114:2076-82. [DOI: 10.1021/jp911681n] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Jeff Wereszczynski
- Department of Chemistry, University of California, Irvine, California 92697
| | - Ioan Andricioaei
- Department of Chemistry, University of California, Irvine, California 92697
| |
Collapse
|
28
|
Banáš P, Jurečka P, Walter NG, Šponer J, Otyepka M. Theoretical studies of RNA catalysis: hybrid QM/MM methods and their comparison with MD and QM. Methods 2009; 49:202-16. [PMID: 19398008 PMCID: PMC2753711 DOI: 10.1016/j.ymeth.2009.04.007] [Citation(s) in RCA: 75] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2009] [Revised: 04/07/2009] [Accepted: 04/07/2009] [Indexed: 11/28/2022] Open
Abstract
Hybrid QM/MM methods combine the rigor of quantum mechanical (QM) calculations with the low computational cost of empirical molecular mechanical (MM) treatment allowing to capture dynamic properties to probe critical atomistic details of enzyme reactions. Catalysis by RNA enzymes (ribozymes) has only recently begun to be addressed with QM/MM approaches and is thus still a field under development. This review surveys methodology as well as recent advances in QM/MM applications to RNA mechanisms, including those of the HDV, hairpin, and hammerhead ribozymes, as well as the ribosome. We compare and correlate QM/MM results with those from QM and/or molecular dynamics (MD) simulations, and discuss scope and limitations with a critical eye on current shortcomings in available methodologies and computer resources. We thus hope to foster mutual appreciation and facilitate collaboration between experimentalists and theorists to jointly advance our understanding of RNA catalysis at an atomistic level.
Collapse
Affiliation(s)
- Pavel Banáš
- Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, tr. Svobody 26, 771 46 Olomouc, Czech Republic
| | - Petr Jurečka
- Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, tr. Svobody 26, 771 46 Olomouc, Czech Republic
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic
| | - Nils G. Walter
- Department of Chemistry, Single Molecule Analysis Group, University of Michigan, 930 N. University Ave., Ann Arbor, MI 48109-1055, USA
| | - Jiří Šponer
- Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, tr. Svobody 26, 771 46 Olomouc, Czech Republic
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic
| | - Michal Otyepka
- Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc, tr. Svobody 26, 771 46 Olomouc, Czech Republic
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic
| |
Collapse
|
29
|
Ouyang D, Zhang H, Herten DP, Parekh HS, Smith SC. Flexibility of Short-Strand RNA in Aqueous Solution as Revealed by Molecular Dynamics Simulation: Are A-RNA and A´-RNA Distinct Conformational Structures? Aust J Chem 2009. [DOI: 10.1071/ch09090] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
We use molecular dynamics simulations to compare the conformational structure and dynamics of a 21-base pair RNA sequence initially constructed according to the canonical A-RNA and A′-RNA forms in the presence of counterions and explicit water. Our study aims to add a dynamical perspective to the solid-state structural information that has been derived from X-ray data for these two characteristic forms of RNA. Analysis of the three main structural descriptors commonly used to differentiate between the two forms of RNA – namely major groove width, inclination and the number of base pairs in a helical twist – over a 30 ns simulation period reveals a flexible structure in aqueous solution with fluctuations in the values of these structural parameters encompassing the range between the two crystal forms and more. This provides evidence to suggest that the identification of distinct A-RNA and A′-RNA structures, while relevant in the crystalline form, may not be generally relevant in the context of RNA in the aqueous phase. The apparent structural flexibility observed in our simulations is likely to bear ramifications for the interactions of RNA with biological molecules (e.g. proteins) and non-biological molecules (e.g. non-viral gene delivery vectors).
Collapse
|
30
|
Bendic C, Enache M, Volanschi E. Analysis of actinomycin D-DNA model complexes using a quantum-chemical criterion: Mulliken overlap populations. J Mol Graph Model 2008; 24:10-6. [PMID: 15936964 DOI: 10.1016/j.jmgm.2005.03.004] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/10/2005] [Revised: 03/30/2005] [Accepted: 03/31/2005] [Indexed: 10/25/2022]
Abstract
The binding of the antitumoral drug actinomycin D to single- and double-stranded DNA was investigated using molecular modeling in the frame of MM+ molecular mechanics and AM1 semi-empirical method. Two other programs, especially conceived to analyze hydrogen-bonding patterns in biological macromolecules, HBexplore, based on geometrical criteria and SHB_interactions, based on quantum-chemical criteria (Mulliken overlap populations), were also used. The results account for the non-cooperative intercalative binding process previously investigated, and outline the contribution of specific hydrogen bonding as well as CH...O(N) and other atom-atom intermolecular interactions to the stabilization of the actinomycin D-DNA complexes. They also support the hemi-intercalation model proposed in literature for the actinomycin D-ssDNA complex.
Collapse
Affiliation(s)
- Cezar Bendic
- Department of Physical Chemistry, Faculty of Chemistry, University of Bucharest, Blvd. Elisabeta 4-12, Bucharest RO-70346, Romania
| | | | | |
Collapse
|
31
|
Abstract
RNA requires conformational dynamics to undergo its diverse functional roles. Here, a new topological network representation of RNA structures is presented that allows analyzing RNA flexibility/rigidity based on constraint counting. The method extends the FIRST approach, which identifies flexible and rigid regions in atomic detail in a single, static, three-dimensional molecular framework. Initially, the network rigidity of a canonical A-form RNA is analyzed by counting on constraints of network elements of increasing size. These considerations demonstrate that it is the inclusion of hydrophobic contacts into the RNA topological network that is crucial for an accurate flexibility prediction. The counting also explains why a protein-based parameterization results in overly rigid RNA structures. The new network representation is then validated on a tRNA(ASP) structure and all NMR-derived ensembles of RNA structures currently available in the Protein Data Bank (with chain length >/=40). The flexibility predictions demonstrate good agreement with experimental mobility data, and the results are superior compared to predictions based on two previously used network representations. Encouragingly, this holds for flexibility predictions as well as mobility predictions obtained by constrained geometric simulations on these networks. Potential applications of the approach to analyzing the flexibility of DNA and RNA/protein complexes are discussed.
Collapse
|
32
|
Priyakumar UD, Mackerell AD. Atomic detail investigation of the structure and dynamics of DNA.RNA hybrids: a molecular dynamics study. J Phys Chem B 2008; 112:1515-24. [PMID: 18197661 DOI: 10.1021/jp709827m] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
DNA.RNA hybrid duplexes are biologically important molecules and are shown to have potential therapeutic properties. To investigate the relationship between structures, energetics, solvation and RNase H activity of hybrid duplexes in comparison with pure DNA and RNA duplexes, a molecular dynamics study using the CHARMM27 force field was undertaken. The structural properties of all four nucleic acids considered are in very good agreement with the experimental data. The backbone dihedral angles and the puckering of the (deoxy)ribose indicate that the purine rich strands retain their A-/B-like properties but the pyrimidine rich DNA strand undergoes A-B conformational transitions. The minor groove widths of the hybrid structures are narrower than those in the RNA duplex, a requirement for RNase H binding. In addition, sampling of noncanonical phosphodiester backbone dihedrals by the DNA strands, differential solvation properties and helical properties, most notably rise, are suggested to contribute to hybrids being RNase H substrates. Differential RNase H activity toward hybrids containing purine versus pyrimidine rich RNA strands is suggested to be due to sampling of values of the phosphodiester backbone dihedrals in the DNA strands. Notably, the present results indicate that hybrids have decreased flexibility as compared to RNA, in contrast to previous reports.
Collapse
Affiliation(s)
- U Deva Priyakumar
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, Maryland 21201, USA
| | | |
Collapse
|
33
|
Kawai K, Osakada Y, Sugimoto A, Fujitsuka M, Majima T. Hole Transfer Rates in A-Form DNA/2′-OMeRNA Hybrid. Chemistry 2007; 13:2386-91. [PMID: 17163564 DOI: 10.1002/chem.200601210] [Citation(s) in RCA: 16] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
The hole transfer rates in the DNA/DNA B-form duplex and DNA/2'-OMeRNA A-form duplex were measured which occurred in the time range of approximately 100 micros. The hole transfer rates in the A-form duplexes were slower and more strongly dependent on the temperature compared to those in the B-form duplexes, suggesting that the A-form is more rigid than the B-form duplex in this time scale.
Collapse
Affiliation(s)
- Kiyohiko Kawai
- The Institute of Scientific and Industrial Research (SANKEN), Osaka University, Mihogaoka 8-1, Ibaraki, Osaka 567-0047, Japan.
| | | | | | | | | |
Collapse
|
34
|
McDowell SE, Špačková N, Šponer J, Walter NG. Molecular dynamics simulations of RNA: an in silico single molecule approach. Biopolymers 2007; 85:169-84. [PMID: 17080418 PMCID: PMC2018183 DOI: 10.1002/bip.20620] [Citation(s) in RCA: 123] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
RNA molecules are now known to be involved in the processing of genetic information at all levels, taking on a wide variety of central roles in the cell. Understanding how RNA molecules carry out their biological functions will require an understanding of structure and dynamics at the atomistic level, which can be significantly improved by combining computational simulation with experiment. This review provides a critical survey of the state of molecular dynamics (MD) simulations of RNA, including a discussion of important current limitations of the technique and examples of its successful application. Several types of simulations are discussed in detail, including those of structured RNA molecules and their interactions with the surrounding solvent and ions, catalytic RNAs, and RNA-small molecule and RNA-protein complexes. Increased cooperation between theorists and experimentalists will allow expanded judicious use of MD simulations to complement conceptually related single molecule experiments. Such cooperation will open the door to a fundamental understanding of the structure-function relationships in diverse and complex RNA molecules. .
Collapse
Affiliation(s)
- S. Elizabeth McDowell
- Biophysics Research Division, Single Molecule Analysis Group, University of Michigan, 930 N. University Avenue, Ann Arbor, MI 48109-1055
| | - Nad'a Špačková
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Prague
| | - Nils G. Walter
- Department of Chemistry, Single Molecule Analysis Group, University of Michigan, 930 N. University Avenue, Ann Arbor, MI 48109-1055
| |
Collapse
|
35
|
Réblová K, Lankas F, Rázga F, Krasovska MV, Koca J, Sponer J. Structure, dynamics, and elasticity of free 16s rRNA helix 44 studied by molecular dynamics simulations. Biopolymers 2006; 82:504-20. [PMID: 16538608 DOI: 10.1002/bip.20503] [Citation(s) in RCA: 67] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
Molecular dynamics (MD) simulations were employed to investigate the structure, dynamics, and local base-pair step deformability of the free 16S ribosomal helix 44 from Thermus thermophilus and of a canonical A-RNA double helix. While helix 44 is bent in the crystal structure of the small ribosomal subunit, the simulated helix 44 is intrinsically straight. It shows, however, substantial instantaneous bends that are isotropic. The spontaneous motions seen in simulations achieve large degrees of bending seen in the X-ray structure and would be entirely sufficient to allow the dynamics of the upper part of helix 44 evidenced by cryo-electron microscopic studies. Analysis of local base-pair step deformability reveals a patch of flexible steps in the upper part of helix 44 and in the area proximal to the bulge bases, suggesting that the upper part of helix 44 has enhanced flexibility. The simulations identify two conformational substates of the second bulge area (bottom part of the helix) with distinct base pairing. In agreement with nuclear magnetic resonance (NMR) and X-ray studies, a flipped out conformational substate of conserved 1492A is seen in the first bulge area. Molecular dynamics (MD) simulations reveal a number of reversible alpha-gamma backbone flips that correspond to transitions between two known A-RNA backbone families. The flipped substates do not cumulate along the trajectory and lead to a modest transient reduction of helical twist with no significant influence on the overall geometry of the duplexes. Despite their considerable flexibility, the simulated structures are very stable with no indication of substantial force field inaccuracies.
Collapse
Affiliation(s)
- Kamila Réblová
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 61265 Brno, Czech Republic
| | | | | | | | | | | |
Collapse
|
36
|
Salsbury FR, Clodfelter JE, Gentry MB, Hollis T, Scarpinato KD. The molecular mechanism of DNA damage recognition by MutS homologs and its consequences for cell death response. Nucleic Acids Res 2006; 34:2173-85. [PMID: 16648361 PMCID: PMC1450329 DOI: 10.1093/nar/gkl238] [Citation(s) in RCA: 48] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Abstract
We determined the molecular mechanism of cell death response by MutS homologs in distinction to the repair event. Key protein–DNA contacts differ in the interaction of MutS homologs with cisplatinated versus mismatched DNA. Mutational analyses of protein–DNA contacts, which were predicted by molecular dynamics (MD) simulations, were performed. Mutations in suggested interaction sites can affect repair and cell death response independently, and to different extents. A glutamate residue is identified as the key contact with cisplatin-DNA. Mutation of the residue increases cisplatin resistance due to increased non-specific DNA binding. In contrast, the conserved phenylalanine that is instrumental and indispensable for mismatch recognition during repair is not required for cisplatin cytotoxicity. These differences in protein–DNA interactions are translated into localized conformational changes that affect nucleotide requirements and inter-subunit interactions. Specifically, the ability for ATP binding/hydrolysis has little consequence for the MMR-dependent damage response. As a consequence, intersubunit contacts are altered that most likely affect the interaction with downstream proteins. We here describe the interaction of MutS homologs with DNA damage, as it differs from the interaction with a mismatch, and its structural translation into all other functional regions of the protein as a mechanism to initiate cell death response and concomitantly inhibit repair.
Collapse
Affiliation(s)
| | - Jill E. Clodfelter
- Department of Cancer Biology, Wake Forest University School of Medicine, Medical Center BoulevardWinston-Salem, NC 27157, USA
| | - Michael B. Gentry
- Department of Cancer Biology, Wake Forest University School of Medicine, Medical Center BoulevardWinston-Salem, NC 27157, USA
| | - Thomas Hollis
- Department of Biochemistry, Wake Forest University School of Medicine, Medical Center BoulevardWinston-Salem, NC 27157, USA
| | - Karin Drotschmann Scarpinato
- Department of Cancer Biology, Wake Forest University School of Medicine, Medical Center BoulevardWinston-Salem, NC 27157, USA
- To whom correspondence should be addressed. Tel: +1 336 713 4077; Fax: +1 336 716 0255;
| |
Collapse
|
37
|
Barthel A, Zacharias M. Conformational transitions in RNA single uridine and adenosine bulge structures: a molecular dynamics free energy simulation study. Biophys J 2006; 90:2450-62. [PMID: 16399833 PMCID: PMC1403169 DOI: 10.1529/biophysj.105.076158] [Citation(s) in RCA: 33] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Extra unmatched nucleotides (single base bulges) are common structural motifs in folded RNA molecules and can participate in RNA-ligand binding and RNA tertiary structure formation. Often these processes are associated with conformational transitions in the bulge region such as flipping out of the bulge base from an intrahelical stacked toward a looped out state. Knowledge of the flexibility of bulge structures and energetics of conformational transitions is an important prerequisite to better understand the function of this RNA motif. Molecular dynamics simulations were performed on single uridine and adenosine bulge nucleotides at the center of eight basepair RNA molecules and indicated larger flexibility of the bulge bases compared to basepaired regions. The umbrella sampling method was applied to study the bulge base looping out process and accompanying conformational and free energy changes. Looping out toward the major groove resulted in partial disruption of adjacent basepairs and was found to be less favorable compared to looping out toward the minor groove. For both uridine and adenosine bulges, a positive free energy change for full looping out was obtained which was approximately 1.5 kcal mol-1 higher in the case of the adenosine compared to the uridine bulge system. The simulations also indicated stable partially looped out states with the bulge bases located in the RNA minor groove and forming base triples with 5'-neighboring basepairs. In the case of the uridine bulge this state was more stable than the intrahelical stacked bulge structure. Induced looping out toward the minor groove involved crossing of an energy barrier of approximately 3.5 kcal mol-1 before reaching the base triple state. A continuum solvent analysis of intermediate bulge states indicated that electrostatic interactions stabilize looped out and base triple states, whereas van der Waals interactions and nonpolar contributions favor the stacked bulge conformation.
Collapse
Affiliation(s)
- André Barthel
- School of Engineering and Science, International University Bremen, D-28759 Bremen, Germany
| | | |
Collapse
|
38
|
Fohrer J, Hennig M, Carlomagno T. Influence of the 2'-hydroxyl group conformation on the stability of A-form helices in RNA. J Mol Biol 2005; 356:280-7. [PMID: 16376377 DOI: 10.1016/j.jmb.2005.11.043] [Citation(s) in RCA: 41] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2005] [Revised: 11/09/2005] [Accepted: 11/14/2005] [Indexed: 11/18/2022]
Abstract
The 2'-hydroxyl group plays fundamental roles in both the structure and the function of RNA, and is the major determinant of the conformational and thermodynamic differences between RNA and DNA. Here, we report a conformational analysis of 2'-OH groups of the HIV-2 TAR RNA by means of NMR scalar coupling measurements in solution. Our analysis supports the existence of a network of water molecules spanning the minor groove of an RNA A-form helix, as has been suggested on the basis of a high-resolution X-ray study of an RNA duplex. The 2'-OH protons of the lower stem nucleotides of the TAR RNA project either towards the O3' or towards the base, where the 2'-OH group can favorably participate in H-bonding interactions with a water molecule situated in the nucleotide base plane. We observe that the k(ex) rate of the 2'-OH proton with the bulk solvent anti-correlates with the base-pair stability, confirming the involvement of the 2'-OH group in a collective network of H-bonds, which requires the presence of canonical helical secondary structure. The methodology and conformational analysis presented here are broadly applicable and facilitate future studies aimed to correlate the conformation of the 2'-OH group with both the structure and the function of RNA and RNA-ligand complexes.
Collapse
Affiliation(s)
- Jörg Fohrer
- Department of NMR Based Structural Biology, Max Planck Institute for Biophysical Chemistry, Göttingen, Germany
| | | | | |
Collapse
|
39
|
Chatterjee S, Pathmasiri W, Chattopadhyaya J. The 5-Me of thyminyl (T) interaction with the neighboring nucleobases dictate the relative stability of isosequential DNA-RNA hybrid duplexes. Org Biomol Chem 2005; 3:3911-5. [PMID: 16240008 DOI: 10.1039/b511139k] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Relative energetic contributions from the base-pairing [DeltaG(degrees)(bp)] vis-a-vis stacking [DeltaG(degrees)(stacking)] to the total free-energy of stabilization [DeltaG(degrees)(37)] for 14 pairs of isosequential hybrid DNA-RNA duplexes (taken from E. A. Lesnik and S. M. Freier, Biochemistry ,1995, 34, 10807) have been dissected in order to understand the differences in the intrinsic nature of the electrostatic forces that are responsible for the self-assembly of the heteroduplexes compared to homoduplexes. The pKa differences between the monomeric nucleotide 3'-ethylphosphates [(d/rN)pEt] as well as nucleotide 3',5'-bis-ethylphosphates [Etp(d/rN)pEt] in both 2'-deoxy (dN) and ribo (rN) series (N = A/G/C/T/U), as the model donor and acceptor (in which stacking is completely eliminated) mimicking those of the internucleotide monomer building blocks of a duplex, can be qualitatively used (P. Acharya, P. Cheruku, S. Chatterjee, S. Acharya and J. Chattopadhyaya, J. Am. Chem. Soc., 2004, 126, 2862) to understand the strength of base-pairing energies in different DNA-RNA (DR), RNA-DNA (RD), DNA-DNA (DD), and RNA-RNA (RR) duplexes. The study has led us to show the following. (1) As the number of excess %T in DR duplexes compared to the isosequential RD duplexes increase the differences in their thermal stabilization [DeltaT(m)](DR-RD) increase and vice-versa(2) The total relative stabilizations, [DeltaDeltaG(degrees)(37)](DR-RD) among the 14 pairs of isosequential DR and RD duplexes (E. A. Lesnik and S. M. Freier, Biochemistry ,1995, 34, 10807) are wholly dependent on the differences in the number of 5-Me(T) stacking interactions with the nearest-neighbors in the D strands of DR duplexes compared to that of the RD duplexes (3) In the relative stabilization of the DR or RD duplexes differences in the free-energy of stackings [DeltaDeltaG(degrees)(stacking)](DR-RD) play a more significant role than the differences in the free-energy of base-pairing, [DeltaDeltaG(degrees)(bp)](DR-RD). In contradistinction, our experimental data shows that RNA-RNA duplexes are more stable than DNA-DNA duplexes because of larger energy gain from the base-pairing in the former compared to the latter (P. Acharya, P. Cheruku, S. Chatterjee, S. Acharya and J. Chattopadhyaya, J. Am. Chem. Soc., 2004, 126, 2862).
Collapse
Affiliation(s)
- Subhrangsu Chatterjee
- Department of Bioorganic Chemistry, Box 581, Biomedical Center, Uppsala University, SE-75123, Uppsala, Sweden
| | | | | |
Collapse
|
40
|
Dixit SB, Beveridge DL, Case DA, Cheatham TE, Giudice E, Lankas F, Lavery R, Maddocks JH, Osman R, Sklenar H, Thayer KM, Varnai P. Molecular dynamics simulations of the 136 unique tetranucleotide sequences of DNA oligonucleotides. II: sequence context effects on the dynamical structures of the 10 unique dinucleotide steps. Biophys J 2005; 89:3721-40. [PMID: 16169978 PMCID: PMC1366942 DOI: 10.1529/biophysj.105.067397] [Citation(s) in RCA: 189] [Impact Index Per Article: 9.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Molecular dynamics (MD) simulations including water and counterions on B-DNA oligomers containing all 136 unique tetranucleotide basepair steps are reported. The objective is to obtain the calculated dynamical structure for at least two copies of each case, use the results to examine issues with regard to convergence and dynamical stability of MD on DNA, and determine the significance of sequence context effects on all unique dinucleotide steps. This information is essential to understand sequence effects on DNA structure and has implications on diverse problems in the structural biology of DNA. Calculations were carried out on the 136 cases embedded in 39 DNA oligomers with repeating tetranucleotide sequences, capped on both ends by GC pairs and each having a total length of 15 nucleotide pairs. All simulations were carried out using a well-defined state-of-the-art MD protocol, the AMBER suite of programs, and the parm94 force field. In a previous article (Beveridge et al. 2004. Biophysical Journal. 87:3799-3813), the research design, details of the simulation protocol, and informatics issues were described. Preliminary results from 15 ns MD trajectories were presented for the d(CpG) step in all 10 unique sequence contexts. The results indicated the sequence context effects to be small for this step, but revealed that MD on DNA at this length of trajectory is subject to surprisingly persistent cooperative transitions of the sugar-phosphate backbone torsion angles alpha and gamma. In this article, we report detailed analysis of the entire trajectory database and occurrence of various conformational substates and its impact on studies of context effects. The analysis reveals a possible direct correspondence between the sequence-dependent dynamical tendencies of DNA structure and the tendency to undergo transitions that "trap" them in nonstandard conformational substates. The difference in mean of the observed basepair step helicoidal parameter distribution with different flanking sequence sometimes differs by as much as one standard deviation, indicating that the extent of sequence effects could be significant. The observations reveal that the impact of a flexible dinucleotide such as CpG could extend beyond the immediate basepair neighbors. The results in general provide new insight into MD on DNA and the sequence-dependent dynamical structural characteristics of DNA.
Collapse
Affiliation(s)
- Surjit B Dixit
- Chemistry Department and Molecular Biophysics Program, Wesleyan University, Middletown, Connecticut 06459, USA
| | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
41
|
Duchardt E, Schwalbe H. Residue specific ribose and nucleobase dynamics of the cUUCGg RNA tetraloop motif by MNMR 13C relaxation. JOURNAL OF BIOMOLECULAR NMR 2005; 32:295-308. [PMID: 16211483 DOI: 10.1007/s10858-005-0659-x] [Citation(s) in RCA: 39] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/04/2023]
Abstract
The dynamics of the nucleobase and the ribose moieties in a 14-nt RNA cUUCGg hairpin-loop uniformly labeled with 13C and 15N were studied by 13C spin relaxation experiments. R1, R1rho and the 13C-[1H] steady-state NOE of C6 and C1' in pyrimidine and C8 and C1' in purine residues were obtained at 298 K. The relaxation data were analyzed by the model-free formalism to yield dynamic information on timescales of pico-, nano- and milli-seconds. An axially symmetric diffusion tensor with an overall rotational correlation time tau(c) of 2.31 +/- 0.13 ns and an axial ratio of 1.35 +/- 0.02 were determined. Both findings are in agreement with hydrodynamic calculations. For the nucleobase carbons, the validity of different reported 13C chemical shift anisotropy values (Stueber, D. and Grant, D. M., 2002 J. Am. Chem. Soc. 124, 10539-10551; Fiala et al., 2000 J. Biomol. NMR 16, 291-302; Sitkoff, D. and Case, D. A., 1998 Prog. NMR Spectroscopy 32, 165-190) is discussed. The resulting dynamics are in agreement with the structural features of the cUUCGg motif in that all residues are mostly rigid (0.82 < S2 < 0.96) in both the nucleobase and the ribose moiety except for the nucleobase of U7, which is protruding into solution (S2 = 0.76). In general, ribose mobility follows nucleobase dynamics, but is less pronounced. Nucleobase dynamics resulting from the analysis of 13C relaxation rates were found to be in agreement with 15N relaxation data derived dynamic information (Akke et al., 1997 RNA 3, 702-709).
Collapse
Affiliation(s)
- Elke Duchardt
- Institute for Organic Chemistry and Chemical Biology, Center for Biomolecular Magnetic Resonance, Johann Wolfgang Goethe-University Frankfurt, Marie-Curie Str. 11, D-60439, Frankfurt/Main, Germany
| | | |
Collapse
|
42
|
Hart K, Nyström B, Ohman M, Nilsson L. Molecular dynamics simulations and free energy calculations of base flipping in dsRNA. RNA (NEW YORK, N.Y.) 2005; 11:609-618. [PMID: 15811914 PMCID: PMC1370749 DOI: 10.1261/rna.7147805] [Citation(s) in RCA: 38] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/05/2004] [Accepted: 01/26/2005] [Indexed: 05/24/2023]
Abstract
The family of adenosine deaminases acting on RNA (ADARs) targets adenosines in RNA that is mainly double stranded. Some substrates are promiscuously deaminated whereas others, such as the mammalian glutamate receptor B (gluR-B) pre-mRNA, are more selectively deaminated. Many DNA/RNA-base modification enzymes use a base flipping mechanism to be able to reach their target base and it is believed that ADARs function in a similar way. In this study we used molecular dynamics (MD) simulations to describe two sites on the gluR-B pre-mRNA, the selectively targeted R/G site and the nontargeted 46 site, in an attempt to explain the substrate specificity. We used regular MD and also a forced base flipping method with umbrella sampling to calculate the free energy of base opening. Spontaneous opening of the mismatched adenosine was observed for the R/G site but not for the 46 site.
Collapse
Affiliation(s)
- Katarina Hart
- Department of Biosciences at NOVUM, Center for Structural Biochemistry, Karolinska Institutet, SE-141 57 Huddinge, Sweden
| | | | | | | |
Collapse
|
43
|
Chapter 6 Molecular Modeling and Atomistic Simulation of Nucleic Acids. ACTA ACUST UNITED AC 2005. [DOI: 10.1016/s1574-1400(05)01006-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/05/2023]
|
44
|
Pérez A, Noy A, Lankas F, Luque FJ, Orozco M. The relative flexibility of B-DNA and A-RNA duplexes: database analysis. Nucleic Acids Res 2004; 32:6144-51. [PMID: 15562006 PMCID: PMC534633 DOI: 10.1093/nar/gkh954] [Citation(s) in RCA: 111] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
An extensive analysis of structural databases is carried out to investigate the relative flexibility of B-DNA and A-RNA duplexes in crystal form. Our results show that the general anisotropic concept of flexibility is not very useful to compare the deformability of B-DNA and A-RNA duplexes, since the flexibility patterns of B-DNA and A-RNA are quite different. In other words, 'flexibility' is a dangerous word for describing macromolecules, unless it is clearly defined. A few soft essential movements explain most of the natural flexibility of A-RNA, whereas many are necessary for B-DNA. Essential movements occurring in naked B-DNAs are identical to those necessary to deform DNA in DNA-protein complexes, which suggest that evolution has designed DNA-protein complexes so that B-DNA is deformed according to its natural tendency. DNA is generally more flexible, but for some distortions A-RNA is easier to deform. Local stiffness constants obtained for naked B-DNAs and DNA complexes are very close, demonstrating that global distortions in DNA necessary for binding to proteins are the result of the addition of small concerted deformations at the base-pair level. Finally, it is worth noting that in general the picture of the relative deformability of A-RNA and DNA derived from database analysis agrees very well with that derived from state of the art molecular dynamics (MD) simulations.
Collapse
Affiliation(s)
- Alberto Pérez
- Molecular Modelling and Bioinformatics Unit, Institut de Recerca Biomèdica, Parc Científic de Barcelona, Josep Samitier 1-5, Barcelona 08028, Spain
| | | | | | | | | |
Collapse
|
45
|
Noy A, Pérez A, Lankas F, Javier Luque F, Orozco M. Relative Flexibility of DNA and RNA: a Molecular Dynamics Study. J Mol Biol 2004; 343:627-38. [PMID: 15465050 DOI: 10.1016/j.jmb.2004.07.048] [Citation(s) in RCA: 75] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2004] [Revised: 06/30/2004] [Accepted: 07/06/2004] [Indexed: 10/26/2022]
Abstract
State of the art molecular dynamics simulations are used to study the structure, dynamics, molecular interaction properties and flexibility of DNA and RNA duplexes in aqueous solution. Special attention is paid to the deformability of both types of structures, revisiting concepts on the relative flexibility of DNA and RNA duplexes. Our simulations strongly suggest that the concepts of flexibility, rigidity and deformability are much more complex than usually believed, and that it is not always true that DNA is more flexible than RNA.
Collapse
Affiliation(s)
- Agnes Noy
- Molecular Modeling and Bioinformatics Unit, Institut de Recerca Biomèdica, Parc Científic de Barcelona, Josep Samitier 1-5, Barcelona 08028, Spain
| | | | | | | | | |
Collapse
|
46
|
Abstract
Empirical force field-based studies of biological macromolecules are becoming a common tool for investigating their structure-activity relationships at an atomic level of detail. Such studies facilitate interpretation of experimental data and allow for information not readily accessible to experimental methods to be obtained. A large part of the success of empirical force field-based methods is the quality of the force fields combined with the algorithmic advances that allow for more accurate reproduction of experimental observables. Presented is an overview of the issues associated with the development and application of empirical force fields to biomolecular systems. This is followed by a summary of the force fields commonly applied to the different classes of biomolecules; proteins, nucleic acids, lipids, and carbohydrates. In addition, issues associated with computational studies on "heterogeneous" biomolecular systems and the transferability of force fields to a wide range of organic molecules of pharmacological interest are discussed.
Collapse
Affiliation(s)
- Alexander D Mackerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, 20 Penn Street, Baltimore, Maryland 21201, USA.
| |
Collapse
|
47
|
Abstract
In moving towards the simulation of larger nucleic acid assemblies over longer timescales that include more accurate representations of the environment, we are nearing the end of an era characterized by single nanosecond molecular dynamics simulation of nucleic acids. We are excited by the promise and predictability of the modeling methods, yet remain prudently cautious of sampling and force field limitations. Highlights include the accurate representation of subtle drug-DNA interactions, the detailed study of modified and unusual nucleic acid structures, insight into the influence of dynamics on the structure of DNA, and exploration of the interaction of solvent and ions with nucleic acids.
Collapse
Affiliation(s)
- Thomas E Cheatham
- Department of Medicinal Chemistry, University of Utah, 2000 East, 30 South, Skaggs Hall 201, Salt Lake City, Utah 84112, USA.
| |
Collapse
|