Modeling DNA Dynamics under Steady Deforming Forces and Torques

Recombination-independent recognition of DNA homology for meiotic silencing in Neurospora crassa

The pairing of homologous chromosomes represents a critical step of meiosis in nearly all sexually reproducing species. In many organisms, pairing involves chromosomes that remain apparently intact. The mechanistic nature of homology recognition at the basis of such pairing is unknown. Using "meiotic silencing by unpaired DNA" (MSUD) as a model process, we demonstrate the existence of a cardinally different approach to DNA homology recognition in meiosis. The main advantage of MSUD over other experimental systems lies in its ability to identify any relatively short DNA fragment lacking a homologous allelic partner. Here, we show that MSUD does not rely on the canonical mechanism of meiotic recombination, yet it is promoted by REC8, a conserved component of the meiotic cohesion complex. We also show that certain patterns of interspersed homology are recognized as pairable during MSUD. Such patterns need to be colinear and must contain short tracts of sequence identity spaced apart at 21 or 22 base pairs. By using these periodicity values as a guiding parameter in all-atom molecular modeling, we discover that homologous DNA molecules can pair by forming quadruplex-based contacts with an interval of 2.5 helical turns. This process requires right-handed plectonemic coiling and additional conformational changes in the intervening double-helical segments. Our results 1) reconcile genetic and biophysical evidence for the existence of direct homologous double-stranded DNA (dsDNA)-dsDNA pairing, 2) identify a role for this process in initiating RNA interference, and 3) suggest that chromosomes can be cross-matched by a precise mechanism that operates on intact dsDNA molecules.

Discovery and design of soft polymeric bio-inspired materials with multiscale simulations and artificial intelligence

Establishing the “Materials 4.0” paradigm requires intimate knowledge of the virtual space in materials design.

Weak nanoscale chaos and anomalous relaxation in DNA

Anomalous nonexponential relaxation in hydrated biomolecules is commonly attributed to the complexity of the free-energy landscapes, similarly to polymers and glasses. It was found recently that the hydrogen-bond breathing of terminal DNA base pairs exhibits a slow power-law relaxation attributable to weak Hamiltonian chaos, with parameters similar to experimental data. Here, the relationship is studied between this motion and spectroscopic signals measured in DNA with a small molecular photoprobe inserted into the base-pair stack. To this end, the earlier computational approach in combination with an analytical theory is applied to the experimental DNA fragment. It is found that the intensity of breathing dynamics is strongly increased in the internal base pairs that flank the photoprobe, with anomalous relaxation quantitatively close to that in terminal base pairs. A physical mechanism is proposed to explain the coupling between the relaxation of base-pair breathing and the experimental response signal. It is concluded that the algebraic relaxation observed experimentally is very likely a manifestation of weakly chaotic dynamics of hydrogen-bond breathing in the base pairs stacked to the photoprobe and that the weak nanoscale chaos can represent an ubiquitous hidden source of nonexponential relaxation in ultrafast spectroscopy.

Understanding the mechanical response of double-stranded DNA and RNA under constant stretching forces using all-atom molecular dynamics

Multiple biological processes involve the stretching of nucleic acids (NAs). Stretching forces induce local changes in the molecule structure, inhibiting or promoting the binding of proteins, which ultimately affects their functionality. Understanding how a force induces changes in the structure of NAs at the atomic level is a challenge. Here, we use all-atom, microsecond-long molecular dynamics to simulate the structure of dsDNA and dsRNA subjected to stretching forces up to 20 pN. We determine all of the elastic constants of dsDNA and dsRNA and provide an explanation for three striking differences in the mechanical response of these two molecules: the threefold softer stretching constant obtained for dsRNA, the opposite twist-stretch coupling, and its nontrivial force dependence. The lower dsRNA stretching resistance is linked to its more open structure, whereas the opposite twist-stretch coupling of both molecules is due to the very different evolution of molecules' interstrand distance with the stretching force. A reduction of this distance leads to overwinding in dsDNA. In contrast, dsRNA is not able to reduce its interstrand distance and can only elongate by unwinding. Interstrand distance is directly correlated with the slide base-pair parameter and its different behavior in dsDNA and dsRNA traced down to changes in the sugar pucker angle of these NAs.

Algebraic Statistics of Poincaré Recurrences in a DNA Molecule

The statistics of Poincaré recurrences is studied for the base-pair breathing dynamics of an all-atom DNA molecule in a realistic aqueous environment with thousands of degrees of freedom. It is found that at least over five decades in time the decay of recurrences is described by an algebraic law with the Poincaré exponent close to β=1.2. This value is directly related to the correlation decay exponent ν=β-1, which is close to ν≈0.15 observed in the time resolved Stokes shift experiments. By applying the virial theorem we analyze the chaotic dynamics in polynomial potentials and demonstrate analytically that an exponent β=1.2 is obtained assuming the dominance of dipole-dipole interactions in the relevant DNA dynamics. Molecular dynamics simulations also reveal the presence of strong low frequency noise with the exponent η=1.6. We trace parallels with the chaotic dynamics of symplectic maps with a few degrees of freedom characterized by the Poincaré exponent β~1.5.

On the origin of thermal untwisting of DNA

In aqueous solutions, the helical twist of DNA decreases with temperature. This phenomenon was noticed and studied experimentally several decades ago, but its physical origin remains elusive. The present paper shows that the thermal untwisting can be predicted from the specific properties of the torsional elasticity of the double helix revealed in recent computational studies. The temperature coefficient of untwisting estimated using coarse-grained models fitted to all-atom MD data accounts for the experimental results nearly quantitatively. The agreement is further improved with the computed torsional rigidity scaled to remove the discrepancy from experiment. The results confirm that the torsional rigidity of DNA is strongly anharmonic. They indicate that for random DNA, its value grows with small twisting and decreases with untwisting.

Torque transfer coefficient in DNA under torsional stress

In recent years, significant progress in understanding the properties of supercoiled DNA has been obtained due to nanotechniques that made stretching and twisting of single molecules possible. Quantitative interpretation of such experiments requires accurate knowledge of torques inside manipulated DNA. This paper argues that it is not possible to transfer the entire magnitudes of external torques to the twisting stress of the double helix, and that a reducing torque transfer coefficient (TTC < 1) should always be assumed. This assertion agrees with simple physical intuition and is supported by the results of all-atom molecular dynamics (MD) simulations. According to MD, the TTCs around 0.8 are observed in nearly optimal conditions. Reaching higher values requires special efforts and it should be difficult in practice. The TTC can be partially responsible for the persistent discrepancies between the twisting rigidity of DNA measured by different methods.

Local elasticity of strained DNA studied by all-atom simulations

Genomic DNA is constantly subjected to various mechanical stresses arising from its biological functions and cell packaging. If the local mechanical properties of DNA change under torsional and tensional stress, the activity of DNA-modifying proteins and transcription factors can be affected and regulated allosterically. To check this possibility, appropriate steady forces and torques were applied in the course of all-atom molecular dynamics simulations of DNA with AT- and GC-alternating sequences. It is found that the stretching rigidity grows with tension as well as twisting. The torsional rigidity is not affected by stretching, but it varies with twisting very strongly, and differently for the two sequences. Surprisingly, for AT-alternating DNA it passes through a minimum with the average twist close to the experimental value in solution. For this fragment, but not for the GC-alternating sequence, the bending rigidity noticeably changes with both twisting and stretching. The results have important biological implications and shed light on earlier experimental observations.

Keep It Flexible: Driving Macromolecular Rotary Motions in Atomistic Simulations with GROMACS

We describe a versatile method to enforce the rotation of subsets of atoms, e.g., a protein subunit, in molecular dynamics (MD) simulations. In particular, we introduce a "flexible axis" technique that allows realistic flexible adaptions of both the rotary subunit as well as the local rotation axis during the simulation. A variety of useful rotation potentials were implemented for the GROMACS 4.5 MD package. Application to the molecular motor F(1)-ATP synthase demonstrates the advantages of the flexible axis approach over the established fixed axis rotation technique.

Anharmonic torsional stiffness of DNA revealed under small external torques

DNA supercoiling plays an important role in a variety of cellular processes. The torsional stress related to supercoiling may also be involved in gene regulation through the local structure and dynamics of the double helix. To check this possibility, steady torsional stress was applied in the course of all-atom molecular dynamics simulations of two DNA fragments with different base pair sequences. For one fragment, the torsional stiffness significantly varied with small twisting. The effect is traced to sequence-specific asymmetry of local torsional fluctuations, and it should be small in long random DNA due to compensation. In contrast, the stiffness of special short sequences can change significantly, which gives a simple possibility of gene regulation via probabilities of strong fluctuations. These results have important implications for the role of DNA twisting in complexes with transcription factors.