1
|
Zhang Z, Šponer J, Bussi G, Mlýnský V, Šulc P, Simmons CR, Stephanopoulos N, Krepl M. Atomistic Picture of Opening-Closing Dynamics of DNA Holliday Junction Obtained by Molecular Simulations. J Chem Inf Model 2023; 63:2794-2809. [PMID: 37126365 PMCID: PMC10170514 DOI: 10.1021/acs.jcim.3c00358] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/02/2023]
Abstract
Holliday junction (HJ) is a noncanonical four-way DNA structure with a prominent role in DNA repair, recombination, and DNA nanotechnology. By rearranging its four arms, HJ can adopt either closed or open state. With enzymes typically recognizing only a single state, acquiring detailed knowledge of the rearrangement process is an important step toward fully understanding the biological function of HJs. Here, we carried out standard all-atom molecular dynamics (MD) simulations of the spontaneous opening-closing transitions, which revealed complex conformational transitions of HJs with an involvement of previously unconsidered "half-closed" intermediates. Detailed free-energy landscapes of the transitions were obtained by sophisticated enhanced sampling simulations. Because the force field overstabilizes the closed conformation of HJs, we developed a system-specific modification which for the first time allows the observation of spontaneous opening-closing HJ transitions in unbiased MD simulations and opens the possibilities for more accurate HJ computational studies of biological processes and nanomaterials.
Collapse
Affiliation(s)
- Zhengyue Zhang
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 00 Brno, Czech Republic
- CEITEC─Central European Institute of Technology, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic
- National Center for Biomolecular Research, Faculty of Science, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 00 Brno, Czech Republic
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati (SISSA), via Bonomea 265, 34136 Trieste, Italy
| | - Vojtěch Mlýnský
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 00 Brno, Czech Republic
| | - Petr Šulc
- Biodesign Center for Molecular Design and Biomimetics, Arizona State University, 1001 S. McAllister Ave, Tempe, 85287 Arizona, United States
| | - Chad R Simmons
- Biodesign Center for Molecular Design and Biomimetics, Arizona State University, 1001 S. McAllister Ave, Tempe, 85287 Arizona, United States
| | - Nicholas Stephanopoulos
- Biodesign Center for Molecular Design and Biomimetics, Arizona State University, 1001 S. McAllister Ave, Tempe, 85287 Arizona, United States
| | - Miroslav Krepl
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 00 Brno, Czech Republic
- Regional Centre of Advanced Technologies and Materials, Czech Advanced Technology and Research Institute (CATRIN), Palacky University Olomouc, Slechtitelu 241/27, 783 71 Olomouc, Czech Republic
| |
Collapse
|
2
|
Ashwood B, Jones MS, Radakovic A, Khanna S, Lee Y, Sachleben JR, Szostak JW, Ferguson AL, Tokmakoff A. Direct monitoring of the thermodynamics and kinetics of DNA and RNA dinucleotide dehybridization from gaps and overhangs. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.04.10.536266. [PMID: 37090657 PMCID: PMC10120721 DOI: 10.1101/2023.04.10.536266] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/25/2023]
Abstract
Hybridization of short nucleic acid segments (<4 nucleotides) to single-strand templates occurs as a critical intermediate in processes such as non-enzymatic nucleic acid replication and toehold-mediated strand displacement. These templates often contain adjacent duplex segments that stabilize base pairing with single-strand gaps or overhangs, but the thermodynamics and kinetics of hybridization in such contexts are poorly understood due to experimental challenges of probing weak binding and rapid structural dynamics. Here we develop an approach to directly measure the thermodynamics and kinetics of DNA and RNA dinucleotide dehybridization using steady-state and temperature-jump infrared spectroscopy. Our results suggest that dinucleotide binding is stabilized through coaxial stacking interactions with the adjacent duplex segments as well as from potential non-canonical base pairing configurations and structural dynamics of gap and overhang templates revealed using molecular dynamics simulations. We measure timescales for dissociation ranging from 0.2 to 40 µs depending on the template and temperature. Dinucleotide hybridization and dehybridization involves a significant free energy barrier with characteristics resembling that of canonical oligonucleotides. Together, our work provides an initial step for predicting the stability and kinetics of hybridization between short nucleic acid segments and various templates.
Collapse
Affiliation(s)
- Brennan Ashwood
- Department of Chemistry, The University of Chicago, 5735 S. Ellis Avenue, Chicago, IL 60637
- The James Franck Institute and Institute for Biophysical Dynamics, The University of Chicago, 929 East 57 Street, Chicago, Illinois 60637, United States
| | - Michael S Jones
- Pritzker School of Molecular Engineering, The University of Chicago, 5640 South Ellis Avenue, Chicago, Illinois 60637, United States
| | | | - Smayan Khanna
- Pritzker School of Molecular Engineering, The University of Chicago, 5640 South Ellis Avenue, Chicago, Illinois 60637, United States
| | - Yumin Lee
- Department of Chemistry, The University of Chicago, 5735 S. Ellis Avenue, Chicago, IL 60637
| | - Joseph R Sachleben
- Biomolecular NMR Core Facility, Biological Sciences Division, The University of Chicago, Chicago, IL 60637, United States
| | - Jack W Szostak
- Department of Chemistry, The University of Chicago, 5735 S. Ellis Avenue, Chicago, IL 60637
| | - Andrew L Ferguson
- Pritzker School of Molecular Engineering, The University of Chicago, 5640 South Ellis Avenue, Chicago, Illinois 60637, United States
| | - Andrei Tokmakoff
- Department of Chemistry, The University of Chicago, 5735 S. Ellis Avenue, Chicago, IL 60637
- The James Franck Institute and Institute for Biophysical Dynamics, The University of Chicago, 929 East 57 Street, Chicago, Illinois 60637, United States
| |
Collapse
|
3
|
Liebl K, Zacharias M. Toward Force Fields with Improved Base Stacking Descriptions. J Chem Theory Comput 2023; 19:1529-1536. [PMID: 36795949 DOI: 10.1021/acs.jctc.2c01121] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/18/2023]
Abstract
Recent DNA force fields indicate good performance in describing flexibility and structural stability of double-stranded B-DNA. However, it is not clear how accurately base stacking interactions are represented that are critical for simulating structure formation processes and conformational changes. Based on the equilibrium nucleoside association and base pair nicking, we find that the recent Tumuc1 force field improves the description of base stacking compared to previous state-of-the-art force fields. Nevertheless, base pair stacking is still overstabilized compared to experiment. We propose a rapid method to reweight calculated free energies of stacking upon force field modifications in order to generate improved parameters. A decrease of the Lennard-Jones attraction between nucleo-bases alone appears insufficient; however, adjustments in the partial charge distribution on base atoms could help to further improve the force field description of base stacking.
Collapse
Affiliation(s)
- Korbinian Liebl
- Department of Chemistry, Chicago Center for Theoretical Chemistry, Institute for Biophysical Dynamics, and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, United States
| | - Martin Zacharias
- Physics Department and Center of Protein Assemblies, Technical University of Munich, Garching 85748, Germany
| |
Collapse
|
4
|
Gallardo A, Bogart BM, Dutagaci B. Protein-Nucleic Acid Interactions for RNA Polymerase II Elongation Factors by Molecular Dynamics Simulations. J Chem Inf Model 2022; 62:3079-3089. [PMID: 35686985 DOI: 10.1021/acs.jcim.2c00121] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
RNA polymerase II (Pol II) forms a complex with elongation factors to proceed to the elongation stage of the transcription process. In this work, we studied the elongation factor SPT5 and explored the protein-nucleic acid interactions for the isolated systems of KOW1 and KOW4 domains of SPT5 with DNA and RNA, respectively. We performed molecular dynamics (MD) simulations using three commonly used force fields that are CHARMM c36m, AMBER ff14sb, and ff19sb. Simulations showed strong protein-nucleic acid interactions and low electrostatic binding free energies for all force fields used. RNA was found to be highly dynamic with all force fields, while DNA had relatively more stable conformations with the AMBER force fields compared to that with CHARMM. Furthermore, we performed MD simulations of the complete elongation complex using CHARMM c36m and AMBER ff19sb force fields to compare the dynamics and interactions with the isolated systems. Similarly, strong KOW1 and DNA interactions were observed in the complete elongation complex simulations and DNA was further stabilized by a network of interactions involving SPT5-KOW1, SPT4, and rpb2 of Pol II. Overall, our study showed that the differences between CHARMM and AMBER force fields strongly affect the dynamics of the nucleic acids. CHARMM provides highly flexible DNA, while AMBER largely stabilizes the DNA structure. Although the presence of the entire interaction network stabilized the DNA and decreased the differences in the results from the two force fields, the discrepancies of the force fields for smaller systems may reflect their problems in generating accurate dynamics of nucleic acids.
Collapse
Affiliation(s)
- Adan Gallardo
- Department of Molecular and Cell Biology, University of California, Merced, 5200 North Lake Road, Merced, California 95343, United States
| | - Brandon M Bogart
- Department of Molecular and Cell Biology, University of California, Merced, 5200 North Lake Road, Merced, California 95343, United States
| | - Bercem Dutagaci
- Department of Molecular and Cell Biology, University of California, Merced, 5200 North Lake Road, Merced, California 95343, United States
| |
Collapse
|
5
|
Rieger M, Zacharias M. Nearest-Neighbor dsDNA Stability Analysis Using Alchemical Free-Energy Simulations. J Phys Chem B 2022; 126:3640-3647. [PMID: 35549273 DOI: 10.1021/acs.jpcb.2c01138] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The thermodynamic stability of double-stranded (ds)DNA depends on its sequence. It is influenced by the base pairing and stacking with neighboring bases along DNA molecules. Semiempirical schemes are available that allow us to predict the thermodynamic stability of DNA sequences based on empirically derived nearest-neighbor contributions of base pairs formed in the context of all possible nearest-neighbor base pairs. Current molecular dynamics (MD) simulations allow one to simulate the dynamics of DNA molecules in good agreement with experimentally obtained structures and available data on conformational flexibility. However, the suitability of current force field methods to reproduce dsDNA stability and its sequence dependence has been much less well tested. We have employed alchemical free-energy simulations of whole base pair transversions in dsDNA and in unbound single-stranded partner molecules. Such transversions change the sequence context but not the nucleotide content or base pairing in dsDNA and allow a direct comparison with the empirical nearest-neighbor dsDNA stability model. For the alchemical free-energy changes in the unbound single-stranded (ss)DNA partner molecules, we tested different setups assuming either complete unstacking or unrestrained simulations with partial stacking in the unbound ssDNA. The free-energy simulations predicted nearest-neighbor effects of similar magnitude, as observed experimentally but showed overall limited correlation with experimental data. An inaccurate description of stacking interactions and other possible reasons such as the neglect of electronic polarization effects are discussed. The results indicate the need to improve the realistic description of stacking interactions in current molecular mechanic force fields.
Collapse
Affiliation(s)
- Manuel Rieger
- Physics Department and Center of Protein Assemblies, Technical University of Munich, 85748 Garching, Germany
| | - Martin Zacharias
- Physics Department and Center of Protein Assemblies, Technical University of Munich, 85748 Garching, Germany
| |
Collapse
|
6
|
Cofsky JC, Knott GJ, Gee CL, Doudna JA. Crystal structure of an RNA/DNA strand exchange junction. PLoS One 2022; 17:e0263547. [PMID: 35436289 PMCID: PMC9015157 DOI: 10.1371/journal.pone.0263547] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2022] [Accepted: 03/04/2022] [Indexed: 11/29/2022] Open
Abstract
Short segments of RNA displace one strand of a DNA duplex during diverse processes including transcription and CRISPR-mediated immunity and genome editing. These strand exchange events involve the intersection of two geometrically distinct helix types-an RNA:DNA hybrid (A-form) and a DNA:DNA homoduplex (B-form). Although previous evidence suggests that these two helices can stack on each other, it is unknown what local geometric adjustments could enable A-on-B stacking. Here we report the X-ray crystal structure of an RNA-5'/DNA-3' strand exchange junction at an anisotropic resolution of 1.6 to 2.2 Å. The structure reveals that the A-to-B helical transition involves a combination of helical axis misalignment, helical axis tilting and compression of the DNA strand within the RNA:DNA helix, where nucleotides exhibit a mixture of A- and B-form geometry. These structural principles explain previous observations of conformational stability in RNA/DNA exchange junctions, enabling a nucleic acid architecture that is repeatedly populated during biological strand exchange events.
Collapse
Affiliation(s)
- Joshua C. Cofsky
- Department of Molecular and Cell Biology, University of California, Berkeley, Berkeley, California, United States of America
| | - Gavin J. Knott
- Department of Biochemistry & Molecular Biology, Monash Biomedicine Discovery Institute, Monash University, Clayton, VIC, Australia
| | - Christine L. Gee
- Department of Molecular and Cell Biology, University of California, Berkeley, Berkeley, California, United States of America
- California Institute for Quantitative Biosciences (QB3), University of California, Berkeley, Berkeley, California, United States of America
- Howard Hughes Medical Institute, University of California, Berkeley, Berkeley, California, United States of America
| | - Jennifer A. Doudna
- Department of Molecular and Cell Biology, University of California, Berkeley, Berkeley, California, United States of America
- California Institute for Quantitative Biosciences (QB3), University of California, Berkeley, Berkeley, California, United States of America
- Howard Hughes Medical Institute, University of California, Berkeley, Berkeley, California, United States of America
- Department of Chemistry, University of California, Berkeley, Berkeley, California, United States of America
- Molecular Biophysics and Integrated Bioimaging Division, Lawrence Berkeley National Laboratory, Berkeley, Berkeley, California, United States of America
- Innovative Genomics Institute, University of California, Berkeley, Berkeley, California, United States of America
- Gladstone Institutes, University of California, San Francisco, San Francisco, California, United States of America
| |
Collapse
|
7
|
Grotz KK, Schwierz N. Magnesium Force Fields for OPC Water with Accurate Solvation, Ion-Binding, and Water-Exchange Properties: Successful Transfer from SPC/E. J Chem Phys 2022; 156:114501. [DOI: 10.1063/5.0087292] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Magnesium plays a vital role in a large variety of biological processes. To model such processes by molecular dynamics simulations, researchers rely on accurate force field parameters for Mg2+ and water. OPC is one of the most promising water models yielding an improved description of biomolecules in water. The aim of this work is to provide force field parameters for Mg2+ that lead to accurate simulation results in combination with OPC water. Using twelve different Mg2+ parameter sets, that were previously optimized with different water models, we systematically assess the transferability to OPC based on a large variety of experimental properties. The results show that the Mg2+ parameters for SPC/E are transferable to OPC and closely reproduce the experimental solvation free energy, radius of the first hydration shell, coordination number, activity derivative, and binding affinity toward the phosphate oxygens on RNA. Two optimal parameter sets are presented: MicroMg yields water exchange in OPC on the microsecond timescale in agreement with experiments. NanoMg yields accelerated exchange on the nanosecond timescale and facilitates the direct observation of ion binding events for enhanced sampling purposes.
Collapse
Affiliation(s)
- Kara K. Grotz
- Theoretical Biophysics, Max Planck Institute of Biophysics, Germany
| | - Nadine Schwierz
- Theoretical Biophysics, Max Planck Institute of Biophysics, Germany
| |
Collapse
|
8
|
Grupa U, Liebl K, Zacharias M. Orientation Dependence of DNA Blunt-End Stacking Studied by Free-Energy Simulations. J Phys Chem B 2021; 125:13850-13857. [PMID: 34928161 DOI: 10.1021/acs.jpcb.1c07829] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
DNA blunt ends can associate mediated by stacking interactions between the terminal base pairs that form blunt ends. The blunt end association plays a role in DNA repair and recombination processes and can also be of importance for the design of DNA-based nano-materials. Its function depends on the sequence and on the geometric arrangement that leads to stable interaction. For a stacked state, the relative orientation (twisting) of the base pairs is important. Molecular dynamics and advanced sampling simulations were used to calculate free energy change associated with twist changes of the stacked blunt-end base pairs. The calculations reproduce blunt stacking arrangements found in crystal structures of DNA oligonucleotides as free energy minima. To elucidate the physical origin of the stabilization of certain angular arrangements, the interactions between backbone atoms in the blunt-end stack were switched off in additional free energy calculations. It allows us to decipher the contributions to stacking stabilization due to the nucleobases and the backbone and to analyze the sequence dependence of the angular stacking preferences. Good qualitative agreement was also found for the comparison with quantum mechanical calculations. The results may help in the design of novel DNA-based materials.
Collapse
Affiliation(s)
- Ulrich Grupa
- Center of Functional Protein Assemblies and Physics Department, Technical University of Munich, Ernst-Otto-Fischer-Str. 8, 85748 Garching, Germany
| | - Korbinian Liebl
- Center of Functional Protein Assemblies and Physics Department, Technical University of Munich, Ernst-Otto-Fischer-Str. 8, 85748 Garching, Germany
| | - Martin Zacharias
- Center of Functional Protein Assemblies and Physics Department, Technical University of Munich, Ernst-Otto-Fischer-Str. 8, 85748 Garching, Germany
| |
Collapse
|
9
|
Morzy D, Joshi H, Sandler SE, Aksimentiev A, Keyser UF. Membrane Activity of a DNA-Based Ion Channel Depends on the Stability of Its Double-Stranded Structure. NANO LETTERS 2021; 21:9789-9796. [PMID: 34767378 DOI: 10.1021/acs.nanolett.1c03791] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
DNA nanotechnology has emerged as a promising method for designing spontaneously inserting and fully controllable synthetic ion channels. However, both insertion efficiency and stability of existing DNA-based membrane channels leave much room for improvement. Here, we demonstrate an approach to overcoming the unfavorable DNA-lipid interactions that hinder the formation of a stable transmembrane pore. Our all-atom MD simulations and experiments show that the insertion-driving cholesterol modifications can cause fraying of terminal base pairs of nicked DNA constructs, distorting them when embedded in a lipid bilayer. Importantly, we show that DNA nanostructures with no backbone discontinuities form more stable conductive pores and insert into membranes with a higher efficiency than the equivalent nicked constructs. Moreover, lack of nicks allows design and maintenance of membrane-spanning helices in a tilted orientation within the lipid bilayer. Thus, reducing the conformational degrees of freedom of the DNA nanostructures enables better control over their function as synthetic ion channels.
Collapse
Affiliation(s)
- Diana Morzy
- Cavendish Laboratory, University of Cambridge, JJ Thomson Avenue, Cambridge, CB3 0HE, United Kingdom
| | - Himanshu Joshi
- Department of Physics, University of Illinois at Urbana-Champaign, 1110 West Green Street, Urbana, Illinois 61801, United States
| | - Sarah E Sandler
- Cavendish Laboratory, University of Cambridge, JJ Thomson Avenue, Cambridge, CB3 0HE, United Kingdom
| | - Aleksei Aksimentiev
- Department of Physics, University of Illinois at Urbana-Champaign, 1110 West Green Street, Urbana, Illinois 61801, United States
- Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, 405 North Mathews Avenue, Urbana, Illinois 61801, United States
| | - Ulrich F Keyser
- Cavendish Laboratory, University of Cambridge, JJ Thomson Avenue, Cambridge, CB3 0HE, United Kingdom
| |
Collapse
|
10
|
Svidlov A, Drobotenko M, Basov A, Gerasimenko E, Elkina A, Baryshev M, Nechipurenko Y, Dzhimak S. Influence of Environmental Parameters on the Stability of the DNA Molecule. ENTROPY (BASEL, SWITZERLAND) 2021; 23:1446. [PMID: 34828144 PMCID: PMC8622188 DOI: 10.3390/e23111446] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/25/2021] [Revised: 10/27/2021] [Accepted: 10/29/2021] [Indexed: 12/12/2022]
Abstract
Fluctuations in viscosity within the cell nucleus have wide limits. When a DNA molecule passes from the region of high viscosity values to the region of low values, open states, denaturation bubbles, and unweaving of DNA strands can occur. Stabilization of the molecule is provided by energy dissipation-dissipation due to interaction with the environment. Separate sections of a DNA molecule in a twisted state can experience supercoiling stress, which, among other things, is due to complex entropic effects caused by interaction with a solvent. In this work, based on the numerical solution of a mechanical mathematical model for the interferon alpha 17 gene and a fragment of the Drosophila gene, an analysis of the external environment viscosity influence on the dynamics of the DNA molecule and its stability was carried out. It has been shown that an increase in viscosity leads to a rapid stabilization of the angular vibrations of nitrogenous bases, while a decrease in viscosity changes the dynamics of DNA: the rate of change in the angular deviations of nitrogenous bases increases and the angular deformations of the DNA strands increase at each moment of time. These processes lead to DNA instability, which increases with time. Thus, the paper considers the influence of the external environment viscosity on the dissipation of the DNA nitrogenous bases' vibrational motion energy. Additionally, the study on the basis of the described model of the molecular dynamics of physiological processes at different indicators of the rheological behavior of nucleoplasm will allow a deeper understanding of the processes of nonequilibrium physics of an active substance in a living cell to be obtained.
Collapse
Affiliation(s)
- Alexander Svidlov
- Department of Radiophysics and Nanothechnology, Kuban State University, 350040 Krasnodar, Russia; (A.S.); (M.D.); (A.B.); (A.E.); (M.B.)
- Federal Research Center the Southern Scientific Center of the Russian Academy of Sciences, 344006 Rostov-on-Don, Russia
| | - Mikhail Drobotenko
- Department of Radiophysics and Nanothechnology, Kuban State University, 350040 Krasnodar, Russia; (A.S.); (M.D.); (A.B.); (A.E.); (M.B.)
| | - Alexander Basov
- Department of Radiophysics and Nanothechnology, Kuban State University, 350040 Krasnodar, Russia; (A.S.); (M.D.); (A.B.); (A.E.); (M.B.)
- Department of Fundamental and Clinical Biochemistry, Kuban State Medical University, 350063 Krasnodar, Russia
| | - Eugeny Gerasimenko
- Department of Technology of Fats, Cosmetics, Commodity Science, Processes and Devices Kuban State Technological University, 350072 Krasnodar, Russia;
| | - Anna Elkina
- Department of Radiophysics and Nanothechnology, Kuban State University, 350040 Krasnodar, Russia; (A.S.); (M.D.); (A.B.); (A.E.); (M.B.)
- Federal Research Center the Southern Scientific Center of the Russian Academy of Sciences, 344006 Rostov-on-Don, Russia
- Department of Technology of Fats, Cosmetics, Commodity Science, Processes and Devices Kuban State Technological University, 350072 Krasnodar, Russia;
| | - Mikhail Baryshev
- Department of Radiophysics and Nanothechnology, Kuban State University, 350040 Krasnodar, Russia; (A.S.); (M.D.); (A.B.); (A.E.); (M.B.)
- Federal Research Center the Southern Scientific Center of the Russian Academy of Sciences, 344006 Rostov-on-Don, Russia
- Department of Technology of Fats, Cosmetics, Commodity Science, Processes and Devices Kuban State Technological University, 350072 Krasnodar, Russia;
| | - Yury Nechipurenko
- Engelhardt Institute of Molecular Biology, Russian Academy of Sciences, 119991 Moscow, Russia;
| | - Stepan Dzhimak
- Department of Radiophysics and Nanothechnology, Kuban State University, 350040 Krasnodar, Russia; (A.S.); (M.D.); (A.B.); (A.E.); (M.B.)
- Federal Research Center the Southern Scientific Center of the Russian Academy of Sciences, 344006 Rostov-on-Don, Russia
- Department of Technology of Fats, Cosmetics, Commodity Science, Processes and Devices Kuban State Technological University, 350072 Krasnodar, Russia;
| |
Collapse
|
11
|
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] [MESH Headings] [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
|
12
|
Levintov L, Paul S, Vashisth H. Reaction Coordinate and Thermodynamics of Base Flipping in RNA. J Chem Theory Comput 2021; 17:1914-1921. [PMID: 33594886 DOI: 10.1021/acs.jctc.0c01199] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Base flipping is a key biophysical event involved in recognition of various ligands by ribonucleic acid (RNA) molecules. However, the mechanism of base flipping in RNA remains poorly understood, in part due to the lack of atomistic details on complex rearrangements in neighboring bases. In this work, we applied transition path sampling (TPS) methods to study base flipping in a double-stranded RNA (dsRNA) molecule that is known to interact with RNA-editing enzymes through this mechanism. We obtained an ensemble of 1000 transition trajectories to describe the base-flipping process. We used the likelihood maximization method to determine the refined reaction coordinate (RC) consisting of two collective variables (CVs), a distance and a dihedral angle between nucleotides that form stacking interactions with the flipping base. The free energy profile projected along the refined RC revealed three minima, two corresponding to the initial and final states and one for a metastable state. We suggest that the metastable state likely represents a wobbled conformation of nucleobases observed in NMR studies that is often characterized as the flipped state. The analyses of reactive trajectories further revealed that the base flipping is coupled to a global conformational change in a stem-loop of dsRNA.
Collapse
Affiliation(s)
- Lev Levintov
- Department of Chemical Engineering, University of New Hampshire, Durham 03824, New Hampshire, United States
| | - Sanjib Paul
- Department of Chemistry, New York University, New York 10003, New York, United States
| | - Harish Vashisth
- Department of Chemical Engineering, University of New Hampshire, Durham 03824, New Hampshire, United States
| |
Collapse
|
13
|
Zacharias M. Base-Pairing and Base-Stacking Contributions to Double-Stranded DNA Formation. J Phys Chem B 2020; 124:10345-10352. [PMID: 33156627 DOI: 10.1021/acs.jpcb.0c07670] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
Double-stranded (ds)DNA formation and dissociation are of fundamental biological importance. The negative DNA charge influences the dsDNA stability. However, the base pairing and the stacking between neighboring bases are responsible for the sequence-dependent stability of dsDNA. The stability of a dsDNA molecule can be estimated from empirical nearest-neighbor models based on contributions assigned to base-pair steps along the DNA and additional parameters because of DNA termini. In efforts to separate contributions, it has been concluded that base stacking dominates dsDNA stability, whereas base pairing contributes negligibly. Using a different model for dsDNA formation, we reanalyze dsDNA stability contributions and conclude that base stacking contributes already at the level of separate ssDNAs but that pairing contributions drive the dsDNA formation. The theoretical model also predicts that stability contributions of base-pair steps that contain only guanine/cytosine, mixed steps, and steps with only adenine/thymine follow the order 6:5:4, respectively, as expected based on the formed hydrogen bonds. The model is fully consistent with the available stacking data and the nearest-neighbor dsDNA parameters. It allows assigning a narrowly distributed value for the effective free energy contribution per formed hydrogen bond during dsDNA formation of -0.72 kcal·mol-1 based entirely on the experimental data.
Collapse
Affiliation(s)
- Martin Zacharias
- Physics Department T38, Technical University of Munich, 85748 Garching, Germany
| |
Collapse
|
14
|
Cofsky JC, Karandur D, Huang CJ, Witte IP, Kuriyan J, Doudna JA. CRISPR-Cas12a exploits R-loop asymmetry to form double-strand breaks. eLife 2020; 9:e55143. [PMID: 32519675 PMCID: PMC7286691 DOI: 10.7554/elife.55143] [Citation(s) in RCA: 59] [Impact Index Per Article: 14.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2020] [Accepted: 04/08/2020] [Indexed: 12/17/2022] Open
Abstract
Type V CRISPR-Cas interference proteins use a single RuvC active site to make RNA-guided breaks in double-stranded DNA substrates, an activity essential for both bacterial immunity and genome editing. The best-studied of these enzymes, Cas12a, initiates DNA cutting by forming a 20-nucleotide R-loop in which the guide RNA displaces one strand of a double-helical DNA substrate, positioning the DNase active site for first-strand cleavage. However, crystal structures and biochemical data have not explained how the second strand is cut to complete the double-strand break. Here, we detect intrinsic instability in DNA flanking the RNA-3' side of R-loops, which Cas12a can exploit to expose second-strand DNA for cutting. Interestingly, DNA flanking the RNA-5' side of R-loops is not intrinsically unstable. This asymmetry in R-loop structure may explain the uniformity of guide RNA architecture and the single-active-site cleavage mechanism that are fundamental features of all type V CRISPR-Cas systems.
Collapse
Affiliation(s)
- Joshua C Cofsky
- Department of Molecular and Cell Biology, University of California, BerkeleyBerkeleyUnited States
| | - Deepti Karandur
- Department of Molecular and Cell Biology, University of California, BerkeleyBerkeleyUnited States
- California Institute for Quantitative Biosciences (QB3), University of California, BerkeleyBerkeleyUnited States
- Howard Hughes Medical Institute, University of California, BerkeleyBerkeleyUnited States
| | - Carolyn J Huang
- Department of Molecular and Cell Biology, University of California, BerkeleyBerkeleyUnited States
| | - Isaac P Witte
- Department of Molecular and Cell Biology, University of California, BerkeleyBerkeleyUnited States
| | - John Kuriyan
- Department of Molecular and Cell Biology, University of California, BerkeleyBerkeleyUnited States
- California Institute for Quantitative Biosciences (QB3), University of California, BerkeleyBerkeleyUnited States
- Howard Hughes Medical Institute, University of California, BerkeleyBerkeleyUnited States
- Department of Chemistry, University of California, BerkeleyBerkeleyUnited States
- MBIB Division, Lawrence Berkeley National LaboratoryBerkeleyUnited States
| | - Jennifer A Doudna
- Department of Molecular and Cell Biology, University of California, BerkeleyBerkeleyUnited States
- California Institute for Quantitative Biosciences (QB3), University of California, BerkeleyBerkeleyUnited States
- Howard Hughes Medical Institute, University of California, BerkeleyBerkeleyUnited States
- Department of Chemistry, University of California, BerkeleyBerkeleyUnited States
- MBIB Division, Lawrence Berkeley National LaboratoryBerkeleyUnited States
- Innovative Genomics Institute, University of California, BerkeleyBerkeleyUnited States
- Gladstone Institutes, University of California, San FranciscoSan FranciscoUnited States
| |
Collapse
|
15
|
Mlýnský V, Kührová P, Kühr T, Otyepka M, Bussi G, Banáš P, Šponer J. Fine-Tuning of the AMBER RNA Force Field with a New Term Adjusting Interactions of Terminal Nucleotides. J Chem Theory Comput 2020; 16:3936-3946. [PMID: 32384244 DOI: 10.1021/acs.jctc.0c00228] [Citation(s) in RCA: 29] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Determination of RNA structural-dynamic properties is challenging for experimental methods. Thus, atomistic molecular dynamics (MD) simulations represent a helpful technique complementary to experiments. However, contemporary MD methods still suffer from limitations of force fields (ffs), including imbalances in the nonbonded ff terms. We have recently demonstrated that some improvement of state-of-the-art AMBER RNA ff can be achieved by adding a new term for H-bonding called gHBfix, which increases tuning flexibility and reduces risk of side-effects. Still, the first gHBfix version did not fully correct simulations of short RNA tetranucleotides (TNs). TNs are key benchmark systems due to availability of unique NMR data, although giving too much weight on improving TN simulations can easily lead to overfitting to A-form RNA. Here we combine the gHBfix version with another term called tHBfix, which separately treats H-bond interactions formed by terminal nucleotides. This allows to refine simulations of RNA TNs without affecting simulations of other RNAs. The approach is in line with adopted strategy of current RNA ffs, where the terminal nucleotides possess different parameters for terminal atoms than the internal nucleotides. Combination of gHBfix with tHBfix significantly improves the behavior of RNA TNs during well-converged enhanced-sampling simulations using replica exchange with solute tempering. TNs mostly populate canonical A-form like states while spurious intercalated structures are largely suppressed. Still, simulations of r(AAAA) and r(UUUU) TNs show some residual discrepancies with primary NMR data which suggests that future tuning of some other ff terms might be useful. Nevertheless, the tHBfix has a clear potential to improve modeling of key biochemical processes, where interactions of RNA single stranded ends are involved.
Collapse
Affiliation(s)
- Vojtěch Mlýnský
- Institute of Biophysics of the Czech Academy of Sciences, Kralovopolská 135, 612 65 Brno, Czech Republic
| | - Petra Kührová
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Tomáš Kühr
- Department of Computer Science, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Michal Otyepka
- Institute of Biophysics of the Czech Academy of Sciences, Kralovopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati, SISSA, via Bonomea 265, 34136 Trieste, Italy
| | - Pavel Banáš
- Institute of Biophysics of the Czech Academy of Sciences, Kralovopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences, Kralovopolská 135, 612 65 Brno, Czech Republic
| |
Collapse
|
16
|
Liebl K, Zacharias M. How global DNA unwinding causes non-uniform stress distribution and melting of DNA. PLoS One 2020; 15:e0232976. [PMID: 32413048 PMCID: PMC7228070 DOI: 10.1371/journal.pone.0232976] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2020] [Accepted: 04/24/2020] [Indexed: 12/31/2022] Open
Abstract
DNA unwinding is an important process that controls binding of proteins, gene expression and melting of double-stranded DNA. In a series of all-atom MD simulations on two DNA molecules containing a transcription start TATA-box sequence we demonstrate that application of a global restraint on the DNA twisting dramatically changes the coupling between helical parameters and the distribution of deformation energy along the sequence. Whereas only short range nearest-neighbor coupling is observed in the relaxed case, long-range coupling is induced in the globally restrained case. With increased overall unwinding the elastic deformation energy is strongly non-uniformly distributed resulting ultimately in a local melting transition of only the TATA box segment during the simulations. The deformation energy tends to be stored more in cytidine/guanine rich regions associated with a change in conformational substate distribution. Upon TATA box melting the deformation energy is largely absorbed by the melting bubble with the rest of the sequences relaxing back to near B-form. The simulations allow us to characterize the structural changes and the propagation of the elastic energy but also to calculate the associated free energy change upon DNA unwinding up to DNA melting. Finally, we design an Ising model for predicting the local melting transition based on empirical parameters. The direct comparison with the atomistic MD simulations indicates a remarkably good agreement for the predicted necessary torsional stress to induce a melting transition, for the position and length of the melted region and for the calculated associated free energy change between both approaches.
Collapse
Affiliation(s)
- Korbinian Liebl
- Physics Department T38, Technical University of Munich, Garching, Germany
| | - Martin Zacharias
- Physics Department T38, Technical University of Munich, Garching, Germany
- * E-mail:
| |
Collapse
|
17
|
Islam B, Stadlbauer P, Vorlíčková M, Mergny JL, Otyepka M, Šponer J. Stability of Two-Quartet G-Quadruplexes and Their Dimers in Atomistic Simulations. J Chem Theory Comput 2020; 16:3447-3463. [PMID: 32163706 DOI: 10.1021/acs.jctc.9b01068] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
G-quadruplexes (GQs) are four-stranded noncanonical DNA and RNA architectures that can be formed by guanine-rich sequences. The stability of GQs increases with the number of G-quartets, and three G-quartets generally form stable GQs. However, the stability of two-quartet GQs is an open issue. To understand the intrinsic stability of two-quartet GQ stems, we have carried out a series of unbiased molecular dynamics (MD) simulations (505 μs in total) of two- and four-quartet DNA and RNA GQs, with attention paid mainly to parallel-stranded arrangements. We used AMBER DNA parmOL15 and RNA parmOL3 force fields and tested different ion and water models. Two-quartet parallel-stranded DNA GQs unfolded in all the simulations, while the equivalent RNA GQ was stable in most of the simulations. GQs composed of two stacked units of two-quartet GQs were stable for both DNA and RNA. The simulations suggest that a minimum of three quartets are needed to form an intrinsically stable all-anti parallel-stranded DNA GQ. Parallel two-quartet DNA GQ may exist if substantially stabilized by another molecule or structural element, including multimerization. On the other hand, we predict that isolated RNA two-quartet parallel GQs may form, albeit being weakly stable. We also show that ionic parameters and water models should be chosen with caution because some parameter combinations can cause spurious instability of GQ stems. Some in-so-far unnoticed limitations of force-field description of multiple ions inside the GQs are discussed, which compromise the capability of simulations to fully capture the effect of increase in the number of quartets on the GQ stability.
Collapse
Affiliation(s)
- Barira Islam
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Petr Stadlbauer
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Michaela Vorlíčková
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Jean-Louis Mergny
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Michal Otyepka
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic.,Department of Physical Chemistry, Faculty of Science, Palacký University, tř. 17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| |
Collapse
|
18
|
Lee JY, Kim YJ, Lee C, Lee JG, Yagyu H, Tabata O, Kim DN. Investigating the sequence-dependent mechanical properties of DNA nicks for applications in twisted DNA nanostructure design. Nucleic Acids Res 2019; 47:93-102. [PMID: 30476210 PMCID: PMC6326809 DOI: 10.1093/nar/gky1189] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2018] [Accepted: 11/09/2018] [Indexed: 01/09/2023] Open
Abstract
DNA nick can be used as a design motif in programming the shape and reconfigurable deformation of synthetic DNA nanostructures, but its mechanical properties have rarely been systematically characterized at the level of base sequences. Here, we investigated sequence-dependent mechanical properties of DNA nicks through molecular dynamics simulation for a comprehensive set of distinct DNA oligomers constructed using all possible base-pair steps with and without a nick. We found that torsional rigidity was reduced by 28–82% at the nick depending on its sequence and location although bending and stretching rigidities remained similar to those of regular base-pair steps. No significant effect of a nick on mechanically coupled deformation such as the twist-stretch coupling was observed. These results suggest that the primary structural role of nick is the relaxation of torsional constraint by backbones known to be responsible for relatively high torsional rigidity of DNA. Moreover, we experimentally demonstrated the usefulness of quantified nick properties in self-assembling DNA nanostructure design by constructing twisted DNA origami structures to show that sequence design of nicks successfully controls the twist angle of structures. Our study illustrates the importance as well as the opportunities of considering sequence-dependent properties in structural DNA nanotechnology.
Collapse
Affiliation(s)
- Jae Young Lee
- Department of Mechanical and Aerospace Engineering, Seoul National University, Gwanak-ro 1, Gwanak-gu, Seoul 08826, Republic of Korea
| | - Young-Joo Kim
- Department of Mechanical and Aerospace Engineering, Seoul National University, Gwanak-ro 1, Gwanak-gu, Seoul 08826, Republic of Korea
| | - Chanseok Lee
- Department of Mechanical and Aerospace Engineering, Seoul National University, Gwanak-ro 1, Gwanak-gu, Seoul 08826, Republic of Korea
| | - Jae Gyung Lee
- Department of Mechanical and Aerospace Engineering, Seoul National University, Gwanak-ro 1, Gwanak-gu, Seoul 08826, Republic of Korea
| | - Hiromasa Yagyu
- Department of Mechanical Engineering, Kanto Gakuin University, Yokohama 236-8501, Japan
| | - Osamu Tabata
- Department of Micro Engineering, Kyoto University, Kyoto Daigaku-Katsura, Kyoto 615-8540, Japan
| | - Do-Nyun Kim
- Department of Mechanical and Aerospace Engineering, Seoul National University, Gwanak-ro 1, Gwanak-gu, Seoul 08826, Republic of Korea.,Institute of Advanced Machines and Design, Seoul National University, Gwanak-ro 1, Gwanak-gu, Seoul 08826, Republic of Korea
| |
Collapse
|
19
|
Erlenbach N, Grünewald C, Krstic B, Heckel A, Prisner TF. "End-to-end" stacking of small dsRNA. RNA (NEW YORK, N.Y.) 2019; 25:239-246. [PMID: 30404925 PMCID: PMC6348986 DOI: 10.1261/rna.068130.118] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/17/2018] [Accepted: 11/06/2018] [Indexed: 05/08/2023]
Abstract
PELDOR (pulsed electron-electron double resonance) is an established method to study intramolecular distances and can give evidence for conformational changes and flexibilities. However, it can also be used to study intermolecular interactions as for example oligerimization. Here, we used PELDOR to study the "end-to-end" stacking of small double-stranded (ds) RNAs. For this study, the dsRNA molecules were only singly labeled with the spin label TPA to avoid multispin effects and to measure only the intermolecular stacking interactions. It can be shown that small dsRNAs tend to assemble to rod-like structures due to π-π interactions between the base pairs at the end of the strands. On the one hand, these interactions can influence or complicate measurements aimed at the determining of the structure and dynamics of the dsRNA molecule itself. On the other hand, it can be interesting to study such intermolecular stacking interactions in more detail, as for example their dependence on ion concentration. We quantitatively determined the stacking probability as a function of the monovalent NaCl salt and the dsRNA concentration. From these data, the dissociation constant Kd was deduced and found to depend on the ratio between the NaCl salt and dsRNA concentrations. Additionally, the distances and distance distributions obtained predict a model for the stacking geometry of dsRNAs. Introducing a nucleotide overhangs at one end of the dsRNA molecule restricts the stacking to the other end, leading only to dimer formations. Introducing such an overhang at both ends of the dsRNA molecule fully suppresses stacking, as we demonstrate by PELDOR experiments quantitatively.
Collapse
Affiliation(s)
- Nicole Erlenbach
- Institute of Physical and Theoretical Chemistry, Center of Biomolecular Magnetic Resonance, Goethe University, D-60438 Frankfurt am Main, Germany
| | - Christian Grünewald
- Institute of Organic Chemistry and Chemical Biology, Goethe University, D-60438 Frankfurt am Main, Germany
| | - Bisera Krstic
- Institute of Physical and Theoretical Chemistry, Center of Biomolecular Magnetic Resonance, Goethe University, D-60438 Frankfurt am Main, Germany
| | - Alexander Heckel
- Institute of Organic Chemistry and Chemical Biology, Goethe University, D-60438 Frankfurt am Main, Germany
| | - Thomas F Prisner
- Institute of Physical and Theoretical Chemistry, Center of Biomolecular Magnetic Resonance, Goethe University, D-60438 Frankfurt am Main, Germany
| |
Collapse
|
20
|
Li G, Quan Y, Wang X, Liu R, Bie L, Gao J, Zhang HY. Trinucleotide Base Pair Stacking Free Energy for Understanding TF-DNA Recognition and the Functions of SNPs. Front Chem 2019; 6:666. [PMID: 30713839 PMCID: PMC6345724 DOI: 10.3389/fchem.2018.00666] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2018] [Accepted: 12/21/2018] [Indexed: 01/03/2023] Open
Abstract
Single nucleotide polymorphisms (SNPs) affect base pair stacking, which is the primary factor for maintaining the stability of DNA. However, the mechanism of how SNPs lead to phenotype variations is still unclear. In this work, we connected SNPs and base pair stacking by a 3-mer base pair stacking free energy matrix. The SNPs with large base pair stacking free energy differences led to phenotype variations. A molecular dynamics (MD) simulation was then applied. Our results showed that base pair stacking played an important role in the transcription factor (TF)-DNA interaction. Changes in DNA structure mainly originate from TF-DNA interactions, and with the increased base pair stacking free energy, the structure of DNA approaches its free type, although its binding affinity was increased by the SNP. In addition, quantitative models using base pair stacking features revealed that base pair stacking can be used to predict TF binding specificity. As such, our work combined knowledge from bioinformatics and structural biology and provided a new understanding of the relationship between SNPs and phenotype variations. The 3-mer base pair stacking free energy matrix is useful in high-throughput screening of SNPs and predicting TF-DNA binding affinity.
Collapse
Affiliation(s)
- Gen Li
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan, China
| | - Yuan Quan
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan, China
| | - Xiaocong Wang
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan, China
| | - Rong Liu
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan, China
| | - Lihua Bie
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan, China
| | - Jun Gao
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan, China
| | - Hong-Yu Zhang
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan, China
| |
Collapse
|
21
|
Smith LG, Tan Z, Spasic A, Dutta D, Salas-Estrada LA, Grossfield A, Mathews DH. Chemically Accurate Relative Folding Stability of RNA Hairpins from Molecular Simulations. J Chem Theory Comput 2018; 14:6598-6612. [PMID: 30375860 DOI: 10.1021/acs.jctc.8b00633] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
To benchmark RNA force fields, we compared the folding stabilities of three 12-nucleotide hairpin stem loops estimated by simulation to stabilities determined by experiment. We used umbrella sampling and a reaction coordinate of end-to-end (5' to 3' hydroxyl oxygen) distance to estimate the free energy change of the transition from the native conformation to a fully extended conformation with no hydrogen bonds between non-neighboring bases. Each simulation was performed four times using the AMBER FF99+bsc0+χOL3 force field, and each window, spaced at 1 Å intervals, was sampled for 1 μs, for a total of 552 μs of simulation. We compared differences in the simulated free energy changes to analogous differences in free energies from optical melting experiments using thermodynamic cycles where the free energy change between stretched and random coil sequences is assumed to be sequence-independent. The differences between experimental and simulated ΔΔ G° are, on average, 0.98 ± 0.66 kcal/mol, which is chemically accurate and suggests that analogous simulations could be used predictively. We also report a novel method to identify where replica free energies diverge along a reaction coordinate, thus indicating where additional sampling would most improve convergence. We conclude by discussing methods to more economically perform these simulations.
Collapse
Affiliation(s)
- Louis G Smith
- Department of Biochemistry & Biophysics , University of Rochester , Rochester , New York 14642 , United States.,Center for RNA Biology , University of Rochester , Rochester , New York 14642 , United States
| | - Zhen Tan
- Department of Biochemistry & Biophysics , University of Rochester , Rochester , New York 14642 , United States.,Center for RNA Biology , University of Rochester , Rochester , New York 14642 , United States
| | - Aleksandar Spasic
- Department of Biochemistry & Biophysics , University of Rochester , Rochester , New York 14642 , United States.,Center for RNA Biology , University of Rochester , Rochester , New York 14642 , United States
| | - Debapratim Dutta
- Department of Biochemistry & Biophysics , University of Rochester , Rochester , New York 14642 , United States.,Center for RNA Biology , University of Rochester , Rochester , New York 14642 , United States
| | - Leslie A Salas-Estrada
- Department of Biochemistry & Biophysics , University of Rochester , Rochester , New York 14642 , United States
| | - Alan Grossfield
- Department of Biochemistry & Biophysics , University of Rochester , Rochester , New York 14642 , United States
| | - David H Mathews
- Department of Biochemistry & Biophysics , University of Rochester , Rochester , New York 14642 , United States.,Department of Biostatistics and Computational Biology , University of Rochester , Rochester , New York 14642 , United States.,Center for RNA Biology , University of Rochester , Rochester , New York 14642 , United States
| |
Collapse
|
22
|
Chakraborty D, Hori N, Thirumalai D. Sequence-Dependent Three Interaction Site Model for Single- and Double-Stranded DNA. J Chem Theory Comput 2018; 14:3763-3779. [PMID: 29870236 PMCID: PMC6423546 DOI: 10.1021/acs.jctc.8b00091] [Citation(s) in RCA: 41] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
We develop a robust coarse-grained model for single- and double-stranded DNA by representing each nucleotide by three interaction sites (TIS) located at the centers of mass of sugar, phosphate, and base. The resulting TIS model includes base-stacking, hydrogen bond, and electrostatic interactions as well as bond-stretching and bond angle potentials that account for the polymeric nature of DNA. The choices of force constants for stretching and the bending potentials were guided by a Boltzmann inversion procedure using a large representative set of DNA structures extracted from the Protein Data Bank. Some of the parameters in the stacking interactions were calculated using a learning procedure, which ensured that the experimentally measured melting temperatures of dimers are faithfully reproduced. Without any further adjustments, the calculations based on the TIS model reproduce the experimentally measured salt and sequence-dependence of the size of single-stranded DNA (ssDNA), as well as the persistence lengths of poly(dA) and poly(dT) chains. Interestingly, upon application of mechanical force, the extension of poly(dA) exhibits a plateau, which we trace to the formation of stacked helical domains. In contrast, the force-extension curve (FEC) of poly(dT) is entropic in origin and could be described by a standard polymer model. We also show that the persistence length of double-stranded DNA, formed from two complementary ssDNAs, is consistent with the prediction based on the worm-like chain. The persistence length, which decreases with increasing salt concentration, is in accord with the Odijk-Skolnick-Fixman theory intended for stiff polyelectrolyte chains near the rod limit. Our model predicts the melting temperatures of DNA hairpins with excellent accuracy, and we are able to recover the experimentally known sequence-specific trends. The range of applications, which did not require adjusting any parameter after the initial construction based solely on PDB structures and melting profiles of dimers, attests to the transferability and robustness of the TIS model for ssDNA and dsDNA.
Collapse
Affiliation(s)
- Debayan Chakraborty
- Department of Chemistry, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Naoto Hori
- Department of Chemistry, The University of Texas at Austin, Austin, Texas 78712, United States
| | - D. Thirumalai
- Department of Chemistry, The University of Texas at Austin, Austin, Texas 78712, United States
| |
Collapse
|
23
|
Aziz HA, Abdel-Salam ASG, Al-Obaide MAI, Alobydi HW, Al-Humaish S. Kynurenine 3-Monooxygenase Gene Associated With Nicotine Initiation and Addiction: Analysis of Novel Regulatory Features at 5' and 3'-Regions. Front Genet 2018; 9:198. [PMID: 29951083 PMCID: PMC6008986 DOI: 10.3389/fgene.2018.00198] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/29/2018] [Accepted: 05/17/2018] [Indexed: 11/13/2022] Open
Abstract
Tobacco smoking is widespread behavior in Qatar and worldwide and is considered one of the major preventable causes of ill health and death. Nicotine is part of tobacco smoke that causes numerous health risks and is incredibly addictive; it binds to the α7 nicotinic acetylcholine receptor (α7nAChR) in the brain. Recent studies showed α7nAChR involvement in the initiation and addiction of smoking. Kynurenic acid (KA), a significant tryptophan metabolite, is an antagonist of α7nAChR. Inhibition of kynurenine 3-monooxygenase enzyme encoded by KMO enhances the KA levels. Modulating KMO gene expression could be a useful tactic for the treatment of tobacco initiation and dependence. Since KMO regulation is still poorly understood, we aimed to investigate the 5' and 3'-regulatory factors of KMO gene to advance our knowledge to modulate KMO gene expression. In this study, bioinformatics methods were used to identify the regulatory sequences associated with expression of KMO. The displayed differential expression of KMO mRNA in the same tissue and different tissues suggested the specific usage of the KMO multiple alternative promoters. Eleven KMO alternative promoters identified at 5'-regulatory region contain TATA-Box, lack CpG Island (CGI) and showed dinucleotide base-stacking energy values specific to transcription factor binding sites (TFBSs). The structural features of regulatory sequences can influence the transcription process and cell type-specific expression. The uncharacterized LOC105373233 locus coding for non-coding RNA (ncRNA) located on the reverse strand in a convergent manner at the 3'-side of KMO locus. The two genes likely expressed by a promoter that lacks TATA-Box harbor CGI and two TFBSs linked to the bidirectional transcription, the NRF1, and ZNF14 motifs. We identified two types of microRNA (miR) in the uncharacterized LOC105373233 ncRNA, which are like hsa-miR-5096 and hsa-miR-1285-3p and can target the miR recognition element (MRE) in the KMO mRNA. Pairwise sequence alignment identified 52 nucleotides sequence hosting MRE in the KMO 3' UTR untranslated region complementary to the ncRNA LOC105373233 sequence. We speculate that the identified miRs can modulate the KMO expression and together with alternative promoters at the 5'-regulatory region of KMO might contribute to the development of novel diagnostic and therapeutic algorithm for tobacco smoking.
Collapse
Affiliation(s)
- Hassan A Aziz
- College of Arts and Sciences, Qatar University, Doha, Qatar
| | | | - Mohammed A I Al-Obaide
- School of Medicine, Texas Tech University Health Sciences Center, Amarillo, TX, United States
| | | | | |
Collapse
|
24
|
Yoo J, Aksimentiev A. New tricks for old dogs: improving the accuracy of biomolecular force fields by pair-specific corrections to non-bonded interactions. Phys Chem Chem Phys 2018; 20:8432-8449. [PMID: 29547221 DOI: 10.1039/c7cp08185e] [Citation(s) in RCA: 157] [Impact Index Per Article: 26.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
In contrast to ordinary polymers, the vast majority of biological macromolecules adopt highly ordered three-dimensional structures that define their functions. The key to folding of a biopolymer into a unique 3D structure or to an assembly of several biopolymers into a functional unit is a delicate balance between the attractive and repulsive forces that also makes such self-assembly reversible under physiological conditions. The all-atom molecular dynamics (MD) method has emerged as a powerful tool for studies of individual biomolecules and their functional assemblies, encompassing systems of ever increasing complexity. However, advances in parallel computing technology have outpaced the development of the underlying theoretical models-the molecular force fields, pushing the MD method into an untested territory. Recent tests of the MD method have found the most commonly used molecular force fields to be out of balance, overestimating attractive interactions between charged and hydrophobic groups, which can promote artificial aggregation in MD simulations of multi-component protein, nucleic acid, and lipid systems. One route towards improving the force fields is through the NBFIX corrections method, in which the intermolecular forces are calibrated against experimentally measured quantities such as osmotic pressure by making atom pair-specific adjustments to the non-bonded interactions. In this article, we review development of the NBFIX (Non-Bonded FIX) corrections to the AMBER and CHARMM force fields and discuss their implications for MD simulations of electrolyte solutions, dense DNA systems, Holliday junctions, protein folding, and lipid bilayer membranes.
Collapse
Affiliation(s)
- Jejoong Yoo
- Center for the Physics of Living Cells, Department of Physics and Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, 1110 West Green Street, Urbana, IL 61801, USA. and Center for Self-assembly and Complexity, Institute for Basic Science, Pohang, 37363, Republic of Korea
| | - Aleksei Aksimentiev
- Center for the Physics of Living Cells, Department of Physics and Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, 1110 West Green Street, Urbana, IL 61801, USA.
| |
Collapse
|
25
|
Havrila M, Stadlbauer P, Islam B, Otyepka M, Šponer J. Effect of Monovalent Ion Parameters on Molecular Dynamics Simulations of G-Quadruplexes. J Chem Theory Comput 2017; 13:3911-3926. [PMID: 28657760 DOI: 10.1021/acs.jctc.7b00257] [Citation(s) in RCA: 42] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
G-quadruplexes (GQs) are key noncanonical DNA and RNA architectures stabilized by desolvated monovalent cations present in their central channels. We analyze extended atomistic molecular dynamics simulations (∼580 μs in total) of GQs with 11 monovalent cation parametrizations, assessing GQ overall structural stability, dynamics of internal cations, and distortions of the G-tetrad geometries. Majority of simulations were executed with the SPC/E water model; however, test simulations with TIP3P and OPC water models are also reported. The identity and parametrization of ions strongly affect behavior of a tetramolecular d[GGG]4 GQ, which is unstable with several ion parametrizations. The remaining studied RNA and DNA GQs are structurally stable, though the G-tetrad geometries are always deformed by bifurcated H-bonding in a parametrization-specific manner. Thus, basic 10-μs-scale simulations of fully folded GQs can be safely done with a number of cation parametrizations. However, there are parametrization-specific differences and basic force-field errors affecting the quantitative description of ion-tetrad interactions, which may significantly affect studies of the ion-binding processes and description of the GQ folding landscape. Our d[GGG]4 simulations indirectly suggest that such studies will also be sensitive to the water models. During exchanges with bulk water, the Na+ ions move inside the GQs in a concerted manner, while larger relocations of the K+ ions are typically separated. We suggest that the Joung-Cheatham SPC/E K+ parameters represent a safe choice in simulation studies of GQs, though variation of ion parameters can be used for specific simulation goals.
Collapse
Affiliation(s)
- Marek Havrila
- Institute of Biophysics, Academy of Sciences of the Czech Republic , Královopolská 135, 612 65 Brno, Czech Republic.,CEITEC - Central European Institute of Technology, Masaryk University , Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Petr Stadlbauer
- Institute of Biophysics, Academy of Sciences of the Czech Republic , Královopolská 135, 612 65 Brno, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , 17. listopadu 12, 771 46 Olomouc, Czech Republic
| | - Barira Islam
- Institute of Biophysics, Academy of Sciences of the Czech Republic , Královopolská 135, 612 65 Brno, Czech Republic
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , 17. listopadu 12, 771 46 Olomouc, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic , Královopolská 135, 612 65 Brno, Czech Republic.,CEITEC - Central European Institute of Technology, Masaryk University , Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| |
Collapse
|
26
|
Sedova A, Banavali NK. Geometric Patterns for Neighboring Bases Near the Stacked State in Nucleic Acid Strands. Biochemistry 2017; 56:1426-1443. [PMID: 28187685 DOI: 10.1021/acs.biochem.6b01101] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Structural variation in base stacking has been analyzed frequently in isolated double helical contexts for nucleic acids, but not as often in nonhelical geometries or in complex biomolecular environments. In this study, conformations of two neighboring bases near their stacked state in any environment are comprehensively characterized for single-strand dinucleotide (SSD) nucleic acid crystal structure conformations. An ensemble clustering method is used to identify a reduced set of representative stacking geometries based on pairwise distances between select atoms in consecutive bases, with multiple separable conformational clusters obtained for categories divided by nucleic acid type (DNA/RNA), SSD sequence, stacking face orientation, and the presence or absence of a protein environment. For both DNA and RNA, SSD conformations are observed that are either close to the A-form, or close to the B-form, or intermediate between the two forms, or further away from either form, illustrating the local structural heterogeneity near the stacked state. Among this large variety of distinct conformations, several common stacking patterns are observed between DNA and RNA, and between nucleic acids in isolation or in complex with proteins, suggesting that these might be stable stacking orientations. Noncanonical face/face orientations of the two bases are also observed for neighboring bases in the same strand, but their frequency is much lower, with multiple SSD sequences across categories showing no occurrences of such unusual stacked conformations. The resulting reduced set of stacking geometries is directly useful for stacking-energy comparisons between empirical force fields, prediction of plausible localized variations in single-strand structures near their canonical states, and identification of analogous stacking patterns in newly solved nucleic acid containing structures.
Collapse
Affiliation(s)
- Ada Sedova
- Department of Biomedical Sciences, School of Public Health, State University of New York at Albany , Albany, New York 12222, United States
| | - Nilesh K Banavali
- Laboratory of Computational and Structural Biology, Division of Genetics, Wadsworth Center, New York State Department of Health, CMS 2008, Biggs Laboratory, Empire State Plaza, P.O. Box 509, Albany, New York 12201-0509, United States.,Department of Biomedical Sciences, School of Public Health, State University of New York at Albany , Albany, New York 12222, United States
| |
Collapse
|
27
|
Pinamonti G, Zhao J, Condon DE, Paul F, Noè F, Turner DH, Bussi G. Predicting the Kinetics of RNA Oligonucleotides Using Markov State Models. J Chem Theory Comput 2017; 13:926-934. [PMID: 28001394 DOI: 10.1021/acs.jctc.6b00982] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Abstract
Nowadays different experimental techniques, such as single molecule or relaxation experiments, can provide dynamic properties of biomolecular systems, but the amount of detail obtainable with these methods is often limited in terms of time or spatial resolution. Here we use state-of-the-art computational techniques, namely, atomistic molecular dynamics and Markov state models, to provide insight into the rapid dynamics of short RNA oligonucleotides, to elucidate the kinetics of stacking interactions. Analysis of multiple microsecond-long simulations indicates that the main relaxation modes of such molecules can consist of transitions between alternative folded states, rather than between random coils and native structures. After properly removing structures that are artificially stabilized by known inaccuracies of the current RNA AMBER force field, the kinetic properties predicted are consistent with the time scales of previously reported relaxation experiments.
Collapse
Affiliation(s)
- Giovanni Pinamonti
- Scuola Internazionale Superiore di Studi Avanzati, International School for Advanced Studies , 265 Via Bonomea, I-34136 Trieste, Italy
| | - Jianbo Zhao
- Department of Chemistry, University of Rochester , Rochester, New York 14627, United States
| | - David E Condon
- Department of Chemistry, University of Rochester , Rochester, New York 14627, United States
| | - Fabian Paul
- Department for Mathematics and Computer Science, Freie Universität , Arnimallee 6, Berlin 14195, Germany
| | - Frank Noè
- Department for Mathematics and Computer Science, Freie Universität , Arnimallee 6, Berlin 14195, Germany
| | - Douglas H Turner
- Department of Chemistry, University of Rochester , Rochester, New York 14627, United States
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati, International School for Advanced Studies , 265 Via Bonomea, I-34136 Trieste, Italy
| |
Collapse
|
28
|
Bottaro S, Banáš P, Šponer J, Bussi G. Free Energy Landscape of GAGA and UUCG RNA Tetraloops. J Phys Chem Lett 2016; 7:4032-4038. [PMID: 27661094 DOI: 10.1021/acs.jpclett.6b01905] [Citation(s) in RCA: 56] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/08/2023]
Abstract
We report the folding thermodynamics of ccUUCGgg and ccGAGAgg RNA tetraloops using atomistic molecular dynamics simulations. We obtain a previously unreported estimation of the folding free energy using parallel tempering in combination with well-tempered metadynamics. A key ingredient is the use of a recently developed metric distance, eRMSD, as a biased collective variable. We find that the native fold of both tetraloops is not the global free energy minimum using the AmberχOL3 force field. The estimated folding free energies are 30.2 ± 0.5 kJ/mol for UUCG and 7.5 ± 0.6 kJ/mol for GAGA, in striking disagreement with experimental data. We evaluate the viability of all possible one-dimensional backbone force field corrections. We find that disfavoring the gauche+ region of α and ζ angles consistently improves the existing force field. The level of accuracy achieved with these corrections, however, cannot be considered sufficient by judging on the basis of available thermodynamic data and solution experiments.
Collapse
Affiliation(s)
- Sandro Bottaro
- SISSA, International School for Advanced Studies 265 , Via Bonomea, I-34136 Trieste, Italy
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc , 17. Listopadu 12, 771 46 Olomouc, Czech Republic
- Institute of Biophysics, Academy of Sciences of the Czech Republic , Kralovopolska 135, 612 65 Brno, Czech Republic
| | - Jiří Šponer
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University Olomouc , 17. Listopadu 12, 771 46 Olomouc, Czech Republic
- Institute of Biophysics, Academy of Sciences of the Czech Republic , Kralovopolska 135, 612 65 Brno, Czech Republic
| | - Giovanni Bussi
- SISSA, International School for Advanced Studies 265 , Via Bonomea, I-34136 Trieste, Italy
| |
Collapse
|
29
|
Izadi S, Onufriev AV. Accuracy limit of rigid 3-point water models. J Chem Phys 2016; 145:074501. [PMID: 27544113 PMCID: PMC4991989 DOI: 10.1063/1.4960175] [Citation(s) in RCA: 167] [Impact Index Per Article: 20.9] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2016] [Accepted: 07/19/2016] [Indexed: 11/14/2022] Open
Abstract
Classical 3-point rigid water models are most widely used due to their computational efficiency. Recently, we introduced a new approach to constructing classical rigid water models [S. Izadi et al., J. Phys. Chem. Lett. 5, 3863 (2014)], which permits a virtually exhaustive search for globally optimal model parameters in the sub-space that is most relevant to the electrostatic properties of the water molecule in liquid phase. Here we apply the approach to develop a 3-point Optimal Point Charge (OPC3) water model. OPC3 is significantly more accurate than the commonly used water models of same class (TIP3P and SPCE) in reproducing a comprehensive set of liquid bulk properties, over a wide range of temperatures. Beyond bulk properties, we show that OPC3 predicts the intrinsic charge hydration asymmetry (CHA) of water - a characteristic dependence of hydration free energy on the sign of the solute charge - in very close agreement with experiment. Two other recent 3-point rigid water models, TIP3PFB and H2ODC, each developed by its own, completely different optimization method, approach the global accuracy optimum represented by OPC3 in both the parameter space and accuracy of bulk properties. Thus, we argue that an accuracy limit of practical 3-point rigid non-polarizable models has effectively been reached; remaining accuracy issues are discussed.
Collapse
Affiliation(s)
- Saeed Izadi
- Department of Biomedical Engineering and Mechanics, Virginia Tech, Blacksburg, Virginia 24060, USA
| | - Alexey V Onufriev
- Departments of Computer Science and Physics, Center for Soft Matter and Biological Physics, Virginia Tech, Blacksburg, Virginia 24060, USA
| |
Collapse
|