DNA dynamics can only be understood by taking into account its complex mechanical behavior at different length scales. At the micrometer level, the mechanical properties of single DNA molecules have been well-characterized by polymer models and are commonly quantified by a persistence length of 50 nm (~150 bp). However, at the base pair level (~3.4 Å), the dynamics of DNA involves complex molecular mechanisms that are still being deciphered. Here, we review recent single-molecule experiments and molecular dynamics simulations that are providing novel insights into DNA mechanics from such a molecular perspective. We first discuss recent findings on sequence-dependent DNA mechanical properties, including sequences that resist mechanical stress and sequences that can accommodate strong deformations. We then comment on the intricate effects of cytosine methylation and DNA mismatches on DNA mechanics. Finally, we review recently reported differences in the mechanical properties of DNA and double-stranded RNA, the other double-helical carrier of genetic information. A thorough examination of the recent single-molecule literature permits establishing a set of general ‘rules’ that reasonably explain the mechanics of nucleic acids at the base pair level. These simple rules offer an improved description of certain biological systems and might serve as valuable guidelines for future design of DNA and RNA nanostructures.

Significance: The mechanical properties of nucleic acids regulate multiple biological processes ranging from complex chromosome packing to replication of a plasmid. Single-molecule experiments have reported puzzling differences between the mechanical properties of double-stranded DNA (dsDNA) and double-stranded RNA (dsRNA) subjected to force and torque. This study investigates these differences using constant-force, all-atom, microsecond-long molecular dynamics. We provide a physical mechanism that explains the nonintuitive opposite twist-stretch response of these molecules based on the change of the interstrand distance with the stretching force. Changes in interstrand distance are ultimately related to differences in the chemical structure of dsDNA and dsRNA molecules. The methodology and results shown here open the field to explore larger forces to test experimental measurements, and to challenge the predictions given by our simulations.

Sequence-dependent DNA conformation and flexibility play a fundamental role in the specificity of DNA-protein interactions. Here we quantify the DNA crookedness: a sequence-dependent deformation of DNA that consists of periodic bends of the base pair centers chain. Using extensive $100\mu \mathrm{s}$-long, all-atom molecular dynamics simulations, we found that DNA crookedness and its associated flexibility are bijective, which unveils a one-to-one relation between DNA structure and dynamics. This allowed us to build a predictive model to compute the stretch moduli of different DNA sequences from solely their structure. Sequences with very little crookedness show extremely high stretching stiffness and have been previously shown to form unstable nucleosomes and promote gene expression. Interestingly, the crookedness can be tailored by epigenetic modifications, known to affect gene expression. Our results rationalize the idea that the DNA sequence is not only a chemical code, but also a physical one that allows finely regulating its mechanical properties and, possibly, its 3D arrangement inside the cell.

Large-scale molecular dynamics (MD) simulations and atomic force microscopy (AFM) in liquid are combined to characterize the adsorption of Immunoglobulin G (IgG) antibodies over a hydrophobic surface modeled with a three-layer graphene slab. We consider explicitly the water solvent, simulating systems with massive sizes (up to 770 000 atoms), for four different adsorption orientations. Protocols based on steered MD to speed up the protein diffusion stage and to enhance the dehydration process are combined with long simulation times (>150 ns) in order to make sure that the final adsorption states correspond to actual stable configurations. Our MD results and the AFM images demonstrate that the IgG antibodies are strongly adsorbed, do not unfold, and retain their secondary and tertiary structures upon deposition. Statistical analysis of the AFM images shows that many of the antibodies adopt vertical orientations, even at very small coverages, which expose at least one Fab binding site for recognition events. Single molecule force spectroscopy experiments demonstrate the immunological response of the deposited antibodies by recognizing its specific antigens. The above properties together with the strong anchoring and preservation of the secondary structure, make graphene an excellent candidate for the development of immunosensors.