Twist-bend coupling, twist waves, and the shape of DNA loops

Coarse-grained modeling of DNA-RNA hybrids

We introduce oxNA, a new model for the simulation of DNA-RNA hybrids that is based on two previously developed coarse-grained models-oxDNA and oxRNA. The model naturally reproduces the physical properties of hybrid duplexes, including their structure, persistence length, and force-extension characteristics. By parameterizing the DNA-RNA hydrogen bonding interaction, we fit the model's thermodynamic properties to experimental data using both average-sequence and sequence-dependent parameters. To demonstrate the model's applicability, we provide three examples of its use-calculating the free energy profiles of hybrid strand displacement reactions, studying the resolution of a short R-loop, and simulating RNA-scaffolded wireframe origami.

Influence of temperature on bend, twist and twist-bend coupling of dsDNA

The temperature-dependent bend and twist elasticities of dsDNA, as well as their couplings, were explored through all-atom molecular dynamics simulations. Three rotational parameters, tilt, roll, and twist, were employed to assess the bend and twist elasticities through their stiffness matrix. Our analysis indicates that the bend and twist stiffnesses decrease as the temperature rises, primarily owing to entropic influences stemming from thermodynamic fluctuations. Furthermore, the couplings between these rotational parameters also exhibit a decline with increasing temperature, although the roll-twist coupling displays greater strength than the tilt-roll and tilt-twist couplings, attributed to its more robust correction component. We elucidated the influence of temperature on bend and twist elasticities based on the comparisons between various models and existing data.

Bending Unwinds DNA

In biology, DNA is often tightly bent to small radii. Solely based on the groove asymmetry, a 30-year-old theoretical paper predicted that such bending should unwind DNA, but this effect has not been directly experimentally quantified so far. We developed a ligation-based assay with nicked DNA circles of variable length, thereby decoupling the twist-dependent ligation efficiency from the large bending strain which dominates conventional circularization assays. We demonstrate that tightly bent DNA indeed unwinds to over 11 base pairs/turn, exactly as predicted. Our discovery requires reassessing the molecular mechanisms and energetics of all processes where DNA is tightly bent or relaxed again, including DNA packaging, gene regulation and expression.

Two Different Ways of Stress Release in Supercoiled DNA Minicircles under DNA Nick

It is generally believed that DNA nick is an effective way to release stress in supercoiled DNA, resulting from the twisting motion that individual strands rotate around the axis of the DNA helix. Here, we use MD simulations based on the oxDNA model to investigate the relaxation of 336 bp supercoiled minicircular DNA under DNA nick. Our simulations show that stress release, characterized by the abrupt decrease in linking number, may be induced by two types of DNA motion depending on the nick position. Except for the twisting motion, there is a writhing motion, that is, double strands collectively rotating with one plectoneme removal, which may occur in the process of DNA relaxation with the nick position in the loop region. Moreover, the writhing motion is more likely to occur in the DNA with relatively high hardness, such as C-G pairs. Our simulation results uncover the relationship between structural transformation, stress release, and DNA motion during the dynamic process under DNA nick, indicating the influence of nick position on the relaxation of the supercoiled DNA.

Temperature dependence of DNA elasticity: An all-atom molecular dynamics simulation study

We used all-atom molecular dynamics simulation to investigate the elastic properties of double-stranded DNA (dsDNA). We focused on the influences of temperature on the stretch, bend, and twist elasticities, as well as the twist-stretch coupling, of the dsDNA over a wide range of temperature. The results showed that the bending and twist persistence lengths, together with the stretch and twist moduli, decrease linearly with temperature. However, the twist-stretch coupling behaves in a positive correction and enhances as the temperature increases. The potential mechanisms of how temperature affects dsDNA elasticity and coupling were investigated by using the trajectories from atomistic simulation, in which thermal fluctuations in structural parameters were analyzed in detail. We analyzed the simulation results by comparing them with previous simulation and experimental data, which are in good agreement. The prediction about the temperature dependence of dsDNA elastic properties provides a deeper understanding of DNA elasticities in biological environments and potentially helps in the further development of DNA nanotechnology.

Rigid Base Biasing in Molecular Dynamics Enables Enhanced Sampling of DNA Conformations

All-atom simulations have become increasingly popular to study conformational and dynamical properties of nucleic acids as they are accurate and provide high spatial and time resolutions. This high resolution, however, comes at a heavy computational cost, and, within the time scales of simulations, nucleic acids weakly fluctuate around their ideal structure exploring a limited set of conformations. We introduce the RBB-NA algorithm (available as a package in the Open Source Library PLUMED), which is capable of controlling rigid base parameters in all-atom simulations of nucleic acids. With suitable biasing potentials, this algorithm can "force" a DNA or RNA molecule to assume specific values of the six rotational (tilt, roll, twist, buckle, propeller, opening) and/or the six translational parameters (shift, slide, rise, shear, stretch, stagger). The algorithm enables the use of advanced sampling techniques to probe the structure and dynamics of locally strongly deformed nucleic acids. We illustrate its performance showing some examples in which DNA is strongly twisted, bent, or locally buckled. In these examples, RBB-NA reproduces well the unconstrained simulations data and other known features of DNA mechanics, but it also allows one to explore the anharmonic behavior characterizing the mechanics of nucleic acids in the high deformation regime.

Three-dimensional loop extrusion

Loop extrusion convincingly describes how certain structural maintenance of chromosome (SMC) proteins mediate the formation of large DNA loops. Yet most of the existing computational models cannot reconcile recent in vitro observations showing that condensins can traverse each other, bypass large roadblocks, and perform steps longer than their own size. To fill this gap, we propose a three-dimensional (3D) "trans-grabbing" model for loop extrusion, which not only reproduces the experimental features of loop extrusion by one SMC complex but also predicts the formation of so-called Z-loops via the interaction of two or more SMCs extruding along the same DNA substrate. By performing molecular dynamics simulations of this model, we discover that the experimentally observed asymmetry in the different types of Z-loops is a natural consequence of the DNA tethering in vitro. Intriguingly, our model predicts this bias to disappear in the absence of tethering and a third type of Z-loop, which has not yet been identified in experiments, to appear. Our model naturally explains roadblock bypassing and the appearance of steps larger than the SMC size as a consequence of non-contiguous DNA grabbing. Finally, this study is the first, to our knowledge, to address how Z-loops and bypassing might occur in a way that is broadly consistent with existing cis-only 1D loop extrusion models.

Sequence-dependent twist-bend coupling in DNA minicircles

Looping of double-stranded DNA molecules with 100-200 base pairs into minicircles, catenanes, and rotaxanes has been suggested as a potential tool for DNA nanotechnologies. However, sharp DNA bending into a minicircle with a diameter of several to ten nanometers occurs with alterations in the DNA helical structure and may lead to defective kink formation that hampers the use of DNA minicircles, catenanes, and rotaxanes in nanoscale DNA applications. Here, we investigated local variations of a helical twist in sharply bent DNA using microsecond-long all-atom molecular dynamics simulations of six different DNA minicircles, focusing on the sequence dependence of the coupling between DNA bending and its helical twist. Twist angles between consecutive base pairs were analyzed at different locations relative to the direction of DNA bending and, among 10 unique dinucleotide steps, we identified four dinucleotide steps with strong twist-bend coupling, the pyrimidine-purine dinucleotide steps of TA/TA, CG/CG, and CA/TG and the purine-purine dinucleotide step of GA/TC. This work suggests the sequence-dependent structural responses of DNA to strong mechanical deformation, providing new molecular-level insights into the structure and stability of sharply bent DNA minicircles for nanoscale applications.

Design Approaches and Computational Tools for DNA Nanostructures

Designing a structure in nanoscale with desired shape and properties has been enabled by structural DNA nanotechnology. Design strategies in this research field have evolved to interpret various aspects of increasingly more complex nanoscale assembly and to realize molecular-level functionality by exploring static to dynamic characteristics of the target structure. Computational tools have naturally been of significant interest as they are essential to achieve a fine control over both shape and physicochemical properties of the structure. Here, we review the basic design principles of structural DNA nanotechnology together with its computational analysis and design tools.

Length-scale-dependent elasticity in DNA from coarse-grained and all-atom models

We investigate the influence of nonlocal couplings on the torsional and bending elasticities of DNA. Such couplings have been observed in the past by several simulation studies. Here, we use a description of DNA conformations based on the variables tilt, roll, and twist. Our analysis of both coarse-grained (oxDNA) and all-atom models indicates that these share strikingly similar features: there are strong off-site couplings for tilt-tilt and twist-twist, while they are much weaker in the roll-roll case. By developing an analytical framework to estimate bending and torsional persistence lengths in nonlocal DNA models, we show how off-site interactions generate a length-scale-dependent elasticity. Based on the simulation-generated elasticity data, the theory predicts a significant length-scale-dependent effect on torsional fluctuations but only a modest effect on bending fluctuations. These results are in agreement with experiments probing DNA mechanics from single base pair to kilobase pair scales.

Bistability induced by a spontaneous twisting rate for a two-dimensional intrinsically curved filament

We find that a moderate intrinsic twisting rate (ITR) can induce a bistable state for a force-free two-dimensional intrinsically curved filament. There are two different configurations of equal energy in a bistable state so that the filament is clearly different from its three-dimensional counterpart. The smaller the ITR or the larger the intrinsic curvature (IC), the clearer the distinction between two isoenergetic configurations and the longer the filament. In bistable states, the relationship between length and ITR is approximately a hyperbola and relationship between IC and critical ITR is approximately linear. Thermal fluctuation can result in a shift between two isoenergetic configurations, but large bending and twisting rigidities can prevent the shift and maintain the filament in one of these two configurations. Moreover, a filament can have a metastable state and at a finite temperature such a filament has the similar property as that of a filament with bistable state.

Coarse-grained modelling of DNA plectoneme pinning in the presence of base-pair mismatches

Damaged or mismatched DNA bases result in the formation of physical defects in double-stranded DNA. In vivo, defects in DNA must be rapidly and efficiently repaired to maintain cellular function and integrity. Defects can also alter the mechanical response of DNA to bending and twisting constraints, both of which are important in defining the mechanics of DNA supercoiling. Here, we use coarse-grained molecular dynamics (MD) simulation and supporting statistical-mechanical theory to study the effect of mismatched base pairs on DNA supercoiling. Our simulations show that plectoneme pinning at the mismatch site is deterministic under conditions of relatively high force (>2 pN) and high salt concentration (>0.5 M NaCl). Under physiologically relevant conditions of lower force (0.3 pN) and lower salt concentration (0.2 M NaCl), we find that plectoneme pinning becomes probabilistic and the pinning probability increases with the mismatch size. These findings are in line with experimental observations. The simulation framework, validated with experimental results and supported by the theoretical predictions, provides a way to study the effect of defects on DNA supercoiling and the dynamics of supercoiling in molecular detail.

A review on nonlinear DNA physics

The study and the investigation of structural and dynamical properties of complex systems have attracted considerable interest among scientists in general and physicists and biologists in particular. The present review paper represents a broad overview of the research performed over the nonlinear dynamics of DNA, devoted to some different aspects of DNA physics and including analytical, quantum and computational tools to understand nonlinear DNA physics. We review in detail the semi-discrete approximation within helicoidal Peyrard-Bishop model and show that localized modulated solitary waves, usually called breathers, can emerge and move along the DNA. Since living processes occur at submolecular level, we then discuss a quantum treatment to address the problem of how charge and energy are transported on DNA and how they may play an important role for the functioning of living cells. While this problem has attracted the attention of researchers for a long time, it is still poorly understood how charge and energy transport can occur at distances comparable to the size of macromolecules. Here, we review a theory based on the mechanism of 'self-trapping' of electrons due to their interaction with mechanical (thermal) oscillation of the DNA structure. We also describe recent computational models that have been developed to capture nonlinear mechanics of DNA in vitro and in vivo, possibly under topological constraints. Finally, we provide some conjectures on potential future directions for this field.

SerraNA: a program to determine nucleic acids elasticity from simulation data

AT-rich motifs can generate extreme mechanical properties, which are critical for creating strong global bends when phased properly.