1
|
Fukal J, Zgarbová M, Jurečka P, Šebera J, Sychrovský V. Probabilistic Interpretation of NMR J-Couplings Determines BI-BII State Equilibria in DNA. J Chem Theory Comput 2022; 18:6989-6999. [PMID: 36206364 DOI: 10.1021/acs.jctc.2c00733] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Interpretation of 3JP,H3' NMR scalar spin-spin coupling constants in DNA becomes more reliable by including distinct structural states such as BI and BII, using the weighted-static or, better still, the recently implemented adiabatic-MD (Ad-MD) method. The calculation method employs an adiabatic ("Ad") dependence of 3JP,H3' coupling on NMR-assigned torsion angle, ε, weighted by P(ε) probability distribution calculated by molecular dynamics (MD). Ad-MD calculations enable cross-validation of the bsc1, OL15, and OL21 force fields and various parametrizations of the Karplus equation describing the dependence of 3JP,H3' coupling on ε torsion (KE). The mean absolute deviation of Ad-MD 3JP,H3' couplings from the experimental values in Dickerson-Drew DNA is comparable to the scatter of 3JP,H3' couplings among four separate NMR experiments. A commonly accepted assumption of homogeneity of one kind of structure-dynamic state within DNA (BI or BII) is questionable because the principal characteristics of relevant P(ε) probabilities (shapes and positioning) vary with DNA sequence. The theory outlined in the present work sets limits to future reparameterization of MD force fields, as relevant to NMR data.
Collapse
Affiliation(s)
- Jiří Fukal
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic.,Department of Physical Chemistry, Faculty of Science, Palacký University Olomouc, 17. listopadu 12, 77146 Olomouc, Czech Republic
| | - Marie Zgarbová
- Department of Physical Chemistry, Faculty of Science, Palacký University Olomouc, 17. listopadu 12, 77146 Olomouc, Czech Republic
| | - Petr Jurečka
- Department of Physical Chemistry, Faculty of Science, Palacký University Olomouc, 17. listopadu 12, 77146 Olomouc, Czech Republic
| | - Jakub Šebera
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic
| | - Vladimír Sychrovský
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic.,Department of Electrotechnology, Electrical Engineering, Czech Technical University, Technická 2, 166 27 Praha 6, Czech Republic
| |
Collapse
|
2
|
Savelyev A, MacKerell AD. Differential Deformability of the DNA Minor Groove and Altered BI/BII Backbone Conformational Equilibrium by the Monovalent Ions Li(+), Na(+), K(+), and Rb(+) via Water-Mediated Hydrogen Bonding. J Chem Theory Comput 2015; 11:4473-85. [PMID: 26575937 DOI: 10.1021/acs.jctc.5b00508] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Recently, we reported the differential impact of the monovalent cations Li(+), Na(+), K(+), and Rb(+) on DNA conformational properties. These were identified from variations in the calculated solution-state X-ray DNA spectra as a function of the ion type in solvation buffer in MD simulations using our recently developed polarizable force field based on the classical Drude oscillator. Changes in the DNA structure were found to mainly involve variations in the minor groove width. Because minor groove dimensions vary significantly in protein-DNA complexes and have been shown to play a critical role in both specific and nonspecific DNA readout, understanding the origins of the observed differential DNA modulation by the first-group monovalent ions is of great biological importance. In the present study, we show that the primary microscopic mechanism for the phenomenon is the formation of water-mediated hydrogen bonds between solvated cations located inside the minor groove and simultaneously to two DNA strands, a process whose intensity and impact on DNA structure depends on both the type of ion and the DNA sequence. Additionally, it is shown that the formation of such ion-DNA hydrogen bond complexes appreciably modulates the conformation of the backbone by increasing the population of the BII substate. Notably, the differential impact of the ions on DNA conformational behavior is only predicted by the Drude polarizable model for DNA with virtually no effect observed from MD simulations utilizing the additive CHARMM36 model. Analysis of dipole moments of the water shows the Drude SWM4 model to possess high sensitivity to changes in the local environment, which indicates the important role of electronic polarization in the salt-dependent conformational properties. This also suggests that inclusion of polarization effects is required to model even relatively simple biological systems, such as DNA, in various ionic solutions.
Collapse
Affiliation(s)
- Alexey Savelyev
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland , Baltimore, Maryland 21201, United States
| | - Alexander D MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland , Baltimore, Maryland 21201, United States
| |
Collapse
|
3
|
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
|
4
|
van der Vaart A. Coupled binding-bending-folding: The complex conformational dynamics of protein-DNA binding studied by atomistic molecular dynamics simulations. Biochim Biophys Acta Gen Subj 2014; 1850:1091-1098. [PMID: 25161164 DOI: 10.1016/j.bbagen.2014.08.009] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2014] [Revised: 08/14/2014] [Accepted: 08/18/2014] [Indexed: 12/21/2022]
Abstract
BACKGROUND Protein-DNA binding often involves dramatic conformational changes such as protein folding and DNA bending. While thermodynamic aspects of this behavior are understood, and its biological function is often known, the mechanism by which the conformational changes occur is generally unclear. By providing detailed structural and energetic data, molecular dynamics simulations have been helpful in elucidating and rationalizing protein-DNA binding. SCOPE OF REVIEW This review will summarize recent atomistic molecular dynamics simulations of the conformational dynamics of DNA and protein-DNA binding. A brief overview of recent developments in DNA force fields is given as well. MAJOR CONCLUSIONS Simulations have been crucial in rationalizing the intrinsic flexibility of DNA, and have been instrumental in identifying the sequence of binding events, the triggers for the conformational motion, and the mechanism of binding for a number of important DNA-binding proteins. GENERAL SIGNIFICANCE Molecular dynamics simulations are an important tool for understanding the complex binding behavior of DNA-binding proteins. With recent advances in force fields and rapid increases in simulation time scales, simulations will become even more important for future studies. This article is part of a Special Issue entitled Recent developments of molecular dynamics.
Collapse
Affiliation(s)
- Arjan van der Vaart
- Department of Chemistry, University of South Florida, 4202 East Fowler Avenue CHE 205, Tampa, FL 33620, USA.
| |
Collapse
|
5
|
Jalili S, Karami L, Schofield J. Study of base pair mutations in proline-rich homeodomain (PRH)-DNA complexes using molecular dynamics. EUROPEAN BIOPHYSICS JOURNAL: EBJ 2013; 42:427-40. [PMID: 23385423 DOI: 10.1007/s00249-013-0892-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/04/2012] [Revised: 11/11/2012] [Accepted: 01/21/2013] [Indexed: 11/26/2022]
Abstract
Proline-rich homeodomain (PRH) is a regulatory protein controlling transcription and gene expression processes by binding to the specific sequence of DNA, especially to the sequence 5'-TAATNN-3'. The impact of base pair mutations on the binding between the PRH protein and DNA is investigated using molecular dynamics and free energy simulations to identify DNA sequences that form stable complexes with PRH. Three 20-ns molecular dynamics simulations (PRH-TAATTG, PRH-TAATTA and PRH-TAATGG complexes) in explicit solvent water were performed to investigate three complexes structurally. Structural analysis shows that the native TAATTG sequence forms a complex that is more stable than complexes with base pair mutations. It is also observed that upon mutation, the number and occupancy of the direct and water-mediated hydrogen bonds decrease. Free energy calculations performed with the thermodynamic integration method predict relative binding free energies of 0.64 and 2 kcal/mol for GC to AT and TA to GC mutations, respectively, suggesting that among the three DNA sequences, the PRH-TAATTG complex is more stable than the two mutated complexes. In addition, it is demonstrated that the stability of the PRH-TAATTA complex is greater than that of the PRH-TAATGG complex.
Collapse
Affiliation(s)
- Seifollah Jalili
- Department of Chemistry, K. N. Toosi University of Technology, P.O. Box 15875-4416, Tehran, Iran.
| | | | | |
Collapse
|
6
|
Foloppe N, Guéroult M, Hartmann B. Simulating DNA by molecular dynamics: aims, methods, and validation. Methods Mol Biol 2013; 924:445-468. [PMID: 23034759 DOI: 10.1007/978-1-62703-017-5_17] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/01/2023]
Abstract
The structure and dynamics of the B-DNA double helix involves subtle sequence-dependent effects which are decisive for its function, but difficult to characterize. These structural and dynamic effects can be addressed by simulations of DNA sequences in explicit solvent. Here, we present and discuss the state-of-art of B-DNA molecular dynamics simulations with the major force fields in use today. We explain why a critical analysis of the MD trajectories is required to assess their reliability, and estimate the value and limitations of these models. Overall, simulations of DNA bear great promise towards deciphering the structural and physical subtleties of this biopolymer, where much remains to be understood.
Collapse
|
7
|
Hart K, Foloppe N, Baker CM, Denning EJ, Nilsson L, MacKerell AD. Optimization of the CHARMM additive force field for DNA: Improved treatment of the BI/BII conformational equilibrium. J Chem Theory Comput 2012; 8:348-362. [PMID: 22368531 PMCID: PMC3285246 DOI: 10.1021/ct200723y] [Citation(s) in RCA: 390] [Impact Index Per Article: 32.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
The B-form of DNA can populate two different backbone conformations: BI and BII, defined by the difference between the torsion angles ε and ζ (BI = ε-ζ < 0 and BII = ε-ζ > 0). BI is the most populated state, but the population of the BII state, which is sequence dependent, is significant and accumulating evidence shows that BII affects the overall structure of DNA, and thus influences protein-DNA recognition. This work presents a reparametrization of the CHARMM27 additive nucleic acid force field to increase the sampling of the BII form in MD simulations of DNA. In addition, minor modifications of sugar puckering were introduced to facilitate sampling of the A form of DNA under the appropriate environmental conditions. Parameter optimization was guided by quantum mechanical data on model compounds, followed by calculations on several DNA duplexes in the condensed phase. The selected optimized parameters were then validated against a number of DNA duplexes, with the most extensive tests performed on the EcoRI dodecamer, including comparative calculations using the Amber Parm99bsc0 force field. The new CHARMM model better reproduces experimentally observed sampling of the BII conformation, including sampling as a function of sequence. In addition, the model reproduces the A form of the 1ZF1 duplex in 75 % ethanol, and yields a stable Z-DNA conformation of duplex (GTACGTAC) in its crystal environment. The resulting model, in combination with a recent reoptimization of the CHARMM27 force field for RNA, will be referred to as CHARMM36.
Collapse
Affiliation(s)
- Katarina Hart
- Department of Biosciences and Nutrition, Center for Biosciences, Karolinska Institutet, SE-141 83 HUDDINGE, Sweden
| | | | - Christopher M. Baker
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, 20 Penn Street, Baltimore, Maryland 21201, USA
| | - Elizabeth J. Denning
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, 20 Penn Street, Baltimore, Maryland 21201, USA
| | - Lennart Nilsson
- Department of Biosciences and Nutrition, Center for Biosciences, Karolinska Institutet, SE-141 83 HUDDINGE, Sweden
| | - Alexander D. MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, 20 Penn Street, Baltimore, Maryland 21201, USA
| |
Collapse
|
8
|
Beierlein FR, Kneale GG, Clark T. Predicting the effects of basepair mutations in DNA-protein complexes by thermodynamic integration. Biophys J 2011; 101:1130-8. [PMID: 21889450 DOI: 10.1016/j.bpj.2011.07.003] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2011] [Revised: 06/28/2011] [Accepted: 07/05/2011] [Indexed: 10/17/2022] Open
Abstract
Thermodynamically rigorous free energy methods in principle allow the exact computation of binding free energies in biological systems. Here, we use thermodynamic integration together with molecular dynamics simulations of a DNA-protein complex to compute relative binding free energies of a series of mutants of a protein-binding DNA operator sequence. A guanine-cytosine basepair that interacts strongly with the DNA-binding protein is mutated into adenine-thymine, cytosine-guanine, and thymine-adenine. It is shown that basepair mutations can be performed using a conservative protocol that gives error estimates of ∼10% of the change in free energy of binding. Despite the high CPU-time requirements, this work opens the exciting opportunity of being able to perform basepair scans to investigate protein-DNA binding specificity in great detail computationally.
Collapse
Affiliation(s)
- Frank R Beierlein
- Computer-Chemie-Centrum, Friedrich-Alexander-Universität Erlangen-Nürnberg, Erlangen, Germany
| | | | | |
Collapse
|
9
|
Allnér O, Nilsson L. Nucleotide modifications and tRNA anticodon-mRNA codon interactions on the ribosome. RNA (NEW YORK, N.Y.) 2011; 17:2177-2188. [PMID: 22028366 PMCID: PMC3222130 DOI: 10.1261/rna.029231.111] [Citation(s) in RCA: 37] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/11/2011] [Accepted: 09/14/2011] [Indexed: 05/31/2023]
Abstract
We have carried out molecular dynamics simulations of the tRNA anticodon and mRNA codon, inside the ribosome, to study the effect of the common tRNA modifications cmo(5)U34 and m(6)A37. In tRNA(Val), these modifications allow all four nucleotides to be successfully read at the wobble position in a codon. Previous data suggest that entropic effects are mainly responsible for the extended reading capabilities, but detailed mechanisms have remained unknown. We have performed a wide range of simulations to elucidate the details of these mechanisms at the atomic level and quantify their effects: extensive free energy perturbation coupled with umbrella sampling, entropy calculations of tRNA (free and bound to the ribosome), and thorough structural analysis of the ribosomal decoding center. No prestructuring effect on the tRNA anticodon stem-loop from the two modifications could be observed, but we identified two mechanisms that may contribute to the expanded decoding capability by the modifications: The further reach of the cmo(5)U34 allows an alternative outer conformation to be formed for the noncognate base pairs, and the modification results in increased contacts between tRNA, mRNA, and the ribosome.
Collapse
Affiliation(s)
- Olof Allnér
- Department of Biosciences and Nutrition, Center for Biosciences, Karolinska Institutet, SE-141 83 Huddinge, Sweden
| | - Lennart Nilsson
- Department of Biosciences and Nutrition, Center for Biosciences, Karolinska Institutet, SE-141 83 Huddinge, Sweden
| |
Collapse
|