An In-Depth Look at DNA Crystals through the Prism of Molecular Dynamics Simulations

X-ray crystallography is the primary tool for biomolecular structural determination. However, contacts formed through the crystal lattice are known to affect structures, especially for small and flexible molecules such as DNA oligomers, by introducing significant structural changes in comparison to solution. Furthermore, why molecules crystallize in certain symmetry groups, which role crystallization additives play, and whether they are just innocuous and unspecific crystallization catalysts remain unclear. By using one of the currently best-performing DNA force fields and applying significant computational effort, we described the nature of intermolecular forces that stabilize B-DNA crystals in various symmetry groups and solvent environments with an unprecedented level of detail. We showed a tight coupling between the lattice stability and the type of crystallization additives and that certain symmetry groups are stable only in the presence of a specific additive. Additives and crystal contacts induce small but non-negligible changes in the physical properties of DNA

Allosterism and signal transfer in DNA

We analysed the basic mechanisms of signal transmission in DNA and the origins of the allostery exhibited by systems such as the ternary complex BAMHI-DNA-GRDBD. We found that perturbation information generated by a primary protein binding event travels as a wave to distant regions of DNA following a hopping mechanism. However, such a structural perturbation is transient and does not lead to permanent changes in the DNA geometry and interaction properties at the secondary binding site. The BAMHI-DNA-GRDBD allosteric mechanism does not occur through any traditional models: direct (protein-protein), indirect (reorganization of the secondary site) readout or solvent-release. On the contrary, it is generated by a subtle and less common entropy-mediated mechanism, which might have an important role to explain other DNA-mediated cooperative effects

The structural impact of DNA mismatches

© 2015 © The Author(s) 2015. Published by Oxford University Press on behalf of Nucleic Acids Research. The structure and dynamics of all the transversion and transition mismatches in three different DNA environments have been characterized by molecular dynamics simulations and NMR spectroscopy. We found that the presence of mismatches produced significant local structural alterations, especially in the case of purine transversions. Mismatched pairs often show promiscuous hydrogen bonding patterns, which interchange among each other in the nanosecond time scale. This therefore defines flexible base pairs, where breathing is frequent, and where distortions in helical parameters are strong, resulting in significant alterations in groove dimension. Even if the DNA structure is plastic enough to absorb the structural impact of the mismatch, local structural changes can be propagated far from the mismatch site, following the expected through-backbone and a previously unknown through-space mechanism. The structural changes related to the presence of mismatches help to understand the different susceptibility of mismatches to the action of repairing proteins.Peer Reviewe

Multiscale simulation of DNA

DNA is not only among the most important molecules in life, but a meeting point for biology, physics and chemistry, being studied by numerous techniques. Theoretical methods can help in gaining a detailed understanding of DNA structure and function, but their practical use is hampered by the multiscale nature of this molecule. In this regard, the study of DNA covers a broad range of different topics, from sub-Angstrom details of the electronic distributions of nucleobases, to the mechanical properties of millimeter-long chromatin fibers. Some of the biological processes involving DNA occur in femtoseconds, while others require years. In this review, we describe the most recent theoretical methods that have been considered to study DNA, from the electron to the chromosome, enriching our knowledge on this fascinating molecule

Molecular basis or arginine and lysine DNA sequence-dependent thermo-stability modulation

We have used a variety of theoretical and experimental techniques to study the role of four basic amino acids-Arginine, Lysine, Ornithine and L-2,4-Diaminobutyric acid-on the structure, flexibility and sequence-dependent stability of DNA. We found that the presence of organic ions stabilizes the duplexes and significantly reduces the difference in stability between AT- and GC-rich duplexes with respect to the control conditions. This suggests that these amino acids, ingredients of the primordial soup during abiogenesis, could have helped to equalize the stability of AT- and GC-rich DNA oligomers, facilitating a general noncatalysed self-replication of DNA. Experiments and simulations demonstrate that organic ions have an effect that goes beyond the general electrostatic screening, involving specific interactions along the grooves of the double helix. We conclude that organic ions, largely ignored in the DNA world, should be reconsidered as crucial structural elements far from mimics of small inorganic cations

How B-DNA dynamics decipher sequence-selective protein recognition

The rules governing sequence-specific DNA-protein recognition are under a long-standing debate regarding the prevalence of base versus shape readout mechanisms to explain sequence specificity and of the conformational selection versus induced fit binding paradigms to explain binding-related conformational changes in DNA. Using a combination of atomistic simulations on a subset of representative sequences and mesoscopic simulations at the protein-DNA interactome level, we demonstrate the prevalence of the shape readout model in determining sequence-specificity and of the conformational selection paradigm in defining the general mechanism for binding-related conformational changes in DNA. Our results suggest that the DNA uses a double mechanism to adapt its structure to the protein: it moves along the easiest deformation modes to approach the bioactive conformation, while final adjustments require localized rearrangements at the base-pair step and backbone level. Our study highlights the large impact of B-DNA dynamics in modulating DNA-protein binding

How accurate are accurate force-fields for B-DNA?

Last generation of force-fields are raising expectations on the quality of molecular dynamics (MD) simulations of DNA, as well as to the belief that theoretical models can substitute experimental ones in several cases. However these claims are based on limited benchmarks, where MD simulations have shown the ability to reproduce already existing 'experimental models', which in turn, have an unclear accuracy to represent DNA conformation in solution. In this work we explore the ability of different force-fields to predict the structure of two new B-DNA dodecamers, determined herein by means of 1H nuclear magnetic resonance (NMR). The study allowed us to check directly for experimental NMR observables on duplexes previously not solved, and also to assess the reliability of 'experimental structures'. We observed that technical details in the annealing procedures can induce non-negligible local changes in the final structures. We also found that while not all theoretical simulations are equally reliable, those obtained using last generation of AMBER force-fields (BSC1 and BSC0OL15) show predictive power in the multi-microsecond timescale and can be safely used to reproduce global structure of DNA duplexes and fine sequence-dependent details

The Impact of the HydroxyMethylCytosine epigenetic signature on DNA structure and function.

Funder: Programa de Desarrollo de las Ciencias BasicasFunder: Institució Catalana de Recerca i Estudis AvancatsFunder: Sistema Nacional de Investigadores, Agencia Nacional de Investigación e Innovación, UruguayFunder: Government of SpainWe present a comprehensive, experimental and theoretical study of the impact of 5-hydroxymethylation of DNA cytosine. Using molecular dynamics, biophysical experiments and NMR spectroscopy, we found that Ten-Eleven translocation (TET) dioxygenases generate an epigenetic variant with structural and physical properties similar to those of 5-methylcytosine. Experiments and simulations demonstrate that 5-methylcytosine (mC) and 5-hydroxymethylcytosine (hmC) generally lead to stiffer DNA than normal cytosine, with poorer circularization efficiencies and lower ability to form nucleosomes. In particular, we can rule out the hypothesis that hydroxymethylation reverts to unmodified cytosine physical properties, as hmC is even more rigid than mC. Thus, we do not expect dramatic changes in the chromatin structure induced by differences in physical properties between d(mCpG) and d(hmCpG). Conversely, our simulations suggest that methylated-DNA binding domains (MBDs), associated with repression activities, are sensitive to the substitution d(mCpG) ➔ d(hmCpG), while MBD3 which has a dual activation/repression activity is not sensitive to the d(mCpG) d(hmCpG) change. Overall, while gene activity changes due to cytosine methylation are the result of the combination of stiffness-related chromatin reorganization and MBD binding, those associated to 5-hydroxylation of methylcytosine could be explained by a change in the balance of repression/activation pathways related to differential MBD binding

The structural impact of DNA mismatches

ABSTRACT The structure and dynamics of all the transversion and transition mismatches in three different DNA environments have been characterized by molecular dynamics simulations and NMR spectroscopy. We found that the presence of mismatches produced significant local structural alterations, especially in the case of purine transversions. Mismatched pairs often show promiscuous hydrogen bonding patterns, which interchange among each other in the nanosecond time scale. This therefore defines flexible base pairs, where breathing is frequent, and where distortions in helical parameters are strong, resulting in significant alterations in groove dimension. Even if the DNA structure is plastic enough to absorb the structural impact of the mismatch, local structural changes can be propagated far from the mismatch site, following the expected through-backbone and a previously unknown through-space mechanism. The structural changes related to the presence of mismatches help to understand the different susceptibility of mismatches to the action of repairing proteins

The Role of Unconventional Hydrogen Bonds in Determining BII Propensities in B-DNA

An accurate understanding of DNA backbone transitions is likely to be the key for elucidating the puzzle of the intricate sequence-dependent mechanical properties that govern most of the biologically relevant functions of the double helix. One factor believed to be important in indirect recognition within protein–DNA complexes is the combined effect of two DNA backbone torsions (ε and ζ) which give rise to the well-known BI/BII conformational equilibrium. In this work we explain the sequence-dependent BII propensity observed in RpY steps (R = purine; Y = pyrimidine) at the tetranucleotide level with the help of a previously undetected C–H···O contact between atoms belonging to adjacent bases. Our results are supported by extensive multimicrosecond molecular dynamics simulations from the Ascona B-DNA Consortium, high-level quantum mechanical calculations, and data mining of the experimental structures deposited in the Protein Data Bank