Parmbsc1, a new force field for DNA simulations, was broadly tested on nearly 100 DNA systems and overcame simulation artifacts that affected previous force fields. We present parmbsc1, a force field for DNA atomistic simulation, which has been parameterized from high-level quantum mechanical data and tested for nearly 100 systems (representing a total simulation time of ∼140 μs) covering most of DNA structural space. Parmbsc1 provides high-quality results in diverse systems. Parameters and trajectories are available at http://mmb.irbbarcelona.org/ParmBSC1/ .

We present the results of microsecond molecular dynamics simulations carried out by the ABC group of laboratories on a set of B-DNA oligomers containing the 136 distinct tetranucleotide base sequences. We demonstrate that the resulting trajectories have extensively sampled the conformational space accessible to B-DNA at room temperature. We confirm that base sequence effects depend strongly not only on the specific base pair step, but also on the specific base pairs that flank each step. Beyond sequence effects on average helical parameters and conformational fluctuations, we also identify tetranucleotide sequences that oscillate between several distinct conformational substates. By analyzing the conformation of the phosphodiester backbones, it is possible to understand for which sequences these substates will arise, and what impact they will have on specific helical parameters.

The electric polarizability of DNA, represented by the dielectric constant, is a key intrinsic property that modulates DNA interaction with effector proteins. Surprisingly, it has so far remained unknown owing to the lack of experimental tools able to access it. Here, we experimentally resolved it by detecting the ultraweak polarization forces of DNA inside single T7 bacteriophages particles using electrostatic force microscopy. In contrast to the common assumption of low-polarizable behavior like proteins (ε_r ∼ 2–4), we found that the DNA dielectric constant is ∼8, considerably higher than the value of ∼3 found for capsid proteins. State-of-the-art molecular dynamic simulations confirm the experimental findings, which result in sensibly decreased DNA interaction free energy than normally predicted by Poisson–Boltzmann methods. Our findings reveal a property at the basis of DNA structure and functions that is needed for realistic theoretical descriptions, and illustrate the synergetic power of scanning probe microscopy and theoretical computation techniques. Significance The strength of DNA–DNA and DNA–ligand electrostatic interactions crucially depends on the electric polarizability of DNA, represented by its dielectric constant. This has remained unknown owing to the lack of experimental techniques able to measure it. Here, we experimentally determined the dielectric constant of double-stranded DNA in a native condensed state inside a single bacteriophage as well as the dielectric constants of the protein shell and tail that compose the viral capsid using scanning force microscopy. We supported the experimental data by theoretically determining the DNA dielectric constant using atomistic simulations. Both approaches yield a dielectric constant of DNA around 8, sensibly higher than commonly assumed, thus revealing a DNA intrinsic property essential for realistic computational description 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.

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.