Structural basis of pathway-dependent force profiles in stretched DNA

Quantum machine learning corrects classical forcefields: Stretching DNA base pairs in explicit solvent

In order to improve the accuracy of molecular dynamics simulations, classical forcefields are supplemented with a kernel-based machine learning method trained on quantum-mechanical fragment energies. As an example application, a potential-energy surface is generalized for a small DNA duplex, taking into account explicit solvation and long-range electron exchange-correlation effects. A long-standing problem in molecular science is that experimental studies of the structural and thermodynamic behavior of DNA under tension are not well confirmed by simulation; study of the potential energy vs extension taking into account a novel correction shows that leading classical DNA models have excessive stiffness with respect to stretching. This discrepancy is found to be common across multiple forcefields. The quantum correction is in qualitative agreement with the experimental thermodynamics for larger DNA double helices, providing a candidate explanation for the general and long-standing discrepancy between single molecule stretching experiments and classical calculations of DNA stretching. The new dataset of quantum calculations should facilitate multiple types of nucleic acid simulation, and the associated Kernel Modified Molecular Dynamics method (KMMD) is applicable to biomolecular simulations in general. KMMD is made available as part of the AMBER22 simulation software.

The emergence of sequence-dependent structural motifs in stretched, torsionally constrained DNA

The double-helical structure of DNA results from canonical base pairing and stacking interactions. However, variations from steady-state conformations resulting from mechanical perturbations in cells have physiological relevance but their dependence on sequence remains unclear. Here, we use molecular dynamics simulations showing sequence differences result in markedly different structural motifs upon physiological twisting and stretching. We simulate overextension on different sequences of DNA ((AA)12, (AT)12, (CC)12 and (CG)12) with supercoiling densities at 200 and 50 mM salt concentrations. We find that DNA denatures in the majority of stretching simulations, surprisingly including those with over-twisted DNA. GC-rich sequences are observed to be more stable than AT-rich ones, with the specific response dependent on the base pair order. Furthermore, we find that (AT)12 forms stable periodic structures with non-canonical hydrogen bonds in some regions and non-canonical stacking in others, whereas (CG)12 forms a stacking motif of four base pairs independent of supercoiling density. Our results demonstrate that 20-30% DNA extension is sufficient for breaking B-DNA around and significantly above cellular supercoiling, and that the DNA sequence is crucial for understanding structural changes under mechanical stress. Our findings have important implications for the activities of protein machinery interacting with DNA in all cells.

DNA partitions into triplets under tension in the presence of organic cations, with sequence evolutionary age predicting the stability of the triplet phase

Using atomistic simulations, we show the formation of stable triplet structure when particular GC-rich DNA duplexes are extended in solution over a timescale of hundreds of nanoseconds, in the presence of organic salt. We present planar-stacked triplet disproportionated DNA (Σ DNA) as a possible solution phase of the double helix under tension, subject to sequence and the presence of stabilising co-factors. Considering the partitioning of the duplexes into triplets of base pairs as the first step of operation of recombinase enzymes like RecA, we emphasise the structure-function relationship in Σ DNA. We supplement atomistic calculations with thermodynamic arguments to show that codons for 'phase 1' amino acids (those appearing early in evolution) are more likely than a lower entropy GC-rich sequence to form triplets under tension. We further observe that the four amino acids supposed (in the 'GADV world' hypothesis) to constitute the minimal set to produce functional globular proteins have the strongest triplet-forming propensity within the phase 1 set, showing a series of decreasing triplet propensity with evolutionary newness. The weak form of our observation provides a physical mechanism to minimise read frame and recombination alignment errors in the early evolution of the genetic code.

DNA-Based Single-Molecule Electronics: From Concept to Function

Beyond being the repository of genetic information, DNA is playing an increasingly important role as a building block for molecular electronics. Its inherent structural and molecular recognition properties render it a leading candidate for molecular electronics applications. The structural stability, diversity and programmability of DNA provide overwhelming freedom for the design and fabrication of molecular-scale devices. In the past two decades DNA has therefore attracted inordinate amounts of attention in molecular electronics. This review gives a brief survey of recent experimental progress in DNA-based single-molecule electronics with special focus on single-molecule conductance and I–V characteristics of individual DNA molecules. Existing challenges and exciting future opportunities are also discussed.

Effect of mechanical stretching on DNA conductance

Studying the structural and charge transport properties in DNA is important for unraveling molecular scale processes and developing device applications of DNA molecules. Here we study the effect of mechanical stretching-induced structural changes on charge transport in single DNA molecules. The charge transport follows the hopping mechanism for DNA molecules with lengths varying from 6 to 26 base pairs, but the conductance is highly sensitive to mechanical stretching, showing an abrupt decrease at surprisingly short stretching distances and weak dependence on DNA length. We attribute this force-induced conductance decrease to the breaking of hydrogen bonds in the base pairs at the end of the sequence and describe the data with a mechanical model.

The transition mechanism of DNA overstretching: a microscopic view using molecular dynamics

The overstretching transition in torsionally unconstrained DNA is studied by means of atomistic molecular dynamics simulations. The free-energy profile as a function of the length of the molecule is determined through the umbrella sampling technique providing both a thermodynamic and a structural characterization of the transition pathway. The zero-force free-energy profile is monotonic but, in accordance with recent experimental evidence, becomes two-state at high forces. A number of experimental results are satisfactorily predicted: (i) the entropic and enthalpic contributions to the free-energy difference between the basic (B) state and the extended (S) state; (ii) the longitudinal extension of the transition state and (iii) the enthalpic contribution to the transition barrier. A structural explanation of the experimental finding that overstretching is a cooperative reaction characterized by elementary units of approximately 22 base pairs is found in the average distance between adenine/thymine-rich regions along the molecule. The overstretched DNA adopts a highly dynamical and structurally disordered double-stranded conformation which is characterized by residual base pairing, formation of non-native intra-strand hydrogen bonds and effective hydrophobic screening of apolar regions.