Special issue: Multiscale simulations of DNA from electrons to nucleosomes.

This editorial for Volume 16, Issue 3 of highlights the three-dimensional structural and dynamic information encoded in DNA sequences and introduces the topics covered in this special issue of the journal on Multiscale Simulations of DNA from Electrons to Nucleosomes. is the official journal of the International Union for Pure and Applied Biophysics (IUPAB 2024). The international scope of the articles in the issue exemplifies the goals of IUPAB to organize worldwide advancements, co-operation, communication, and education in biophysics.

cgNA+web : A Visual Interface to the cgNA+ Sequence-dependent Statistical Mechanics Model of Double-stranded Nucleic Acids.

The sequence-dependent statistical mechanics of double-stranded nucleic acid, or dsNA, is believed to be essential in its biological functions. In turn, the equilibrium statistical mechanics behaviour of dsNA depends strongly both on sequence-dependent perturbations in its ground state shape away from an idealised, uniform, double-helical configuration, and on its fluctuations as governed by its sequence-dependent stiffness. We here describe the cgNA+web browser-based interactive tool for visualising the sequence-dependent ground states of dsNA fragments of arbitrary sequences, as predicted by the underlying cgNA+ coarse-grain model.

Sequence-dependent structural properties of B-DNA: what have we learned in 40 years?

The structure of B-DNA, the physiological form of the DNA molecule, has been a central topic in biology, chemistry and physics. Far from uniform and rigid, the double helix was revealed as a flexible and structurally polymorphic molecule. Conformational changes that lead to local and global changes in the helix geometry are mediated by a complex choreography of base and backbone rearrangements affecting the ability of the B-DNA to recognize ligands and consequently on its functionality.

Mechanics of two filaments in tight orthogonal contact.

Networks of flexible filaments often involve regions of tight contact. Predictively understanding the equilibrium configurations of these systems is challenging due to intricate couplings between topology, geometry, large nonlinear deformations, and friction. Here, we perform an in-depth study of a simple, yet canonical, problem that captures the essence of contact between filaments.

The static and dynamic structural heterogeneities of B-DNA: extending Calladine-Dickerson rules.

We present a multi-laboratory effort to describe the structural and dynamical properties of duplex B-DNA under physiological conditions. By processing a large amount of atomistic molecular dynamics simulations, we determine the sequence-dependent structural properties of DNA as expressed in the equilibrium distribution of its stochastic dynamics. Our analysis includes a study of first and second moments of the equilibrium distribution, which can be accurately captured by a harmonic model, but with nonlocal sequence-dependence.

cgDNAweb: a web interface to the cgDNA sequence-dependent coarse-grain model of double-stranded DNA.

The sequence-dependent statistical mechanical properties of fragments of double-stranded DNA is believed to be pertinent to its biological function at length scales from a few base pairs (or bp) to a few hundreds of bp, e.g. indirect read-out protein binding sites, nucleosome positioning sequences, phased A-tracts, etc.

Analyzing DNA curvature and its impact on the ionic environment: application to molecular dynamics simulations of minicircles.

We propose a method for analyzing the magnitude and direction of curvature within nucleic acids, based on the curvilinear helical axis calculated by Curves+. The method is applied to analyzing curvature within minicircles constructed with varying degrees of over- or under-twisting. Using the molecular dynamics trajectories of three different minicircles, we are able to quantify how curvature varies locally both in space and in time.

Sequence-Dependent Persistence Lengths of DNA.

A Monte Carlo code applied to the cgDNA coarse-grain rigid-base model of B-form double-stranded DNA is used to predict a sequence-averaged persistence length of l = 53.5 nm in the sense of Flory, and of l = 160 bp or 53.5 nm in the sense of apparent tangent-tangent correlation decay.

Analyzing ion distributions around DNA: sequence-dependence of potassium ion distributions from microsecond molecular dynamics.

Microsecond molecular dynamics simulations of B-DNA oligomers carried out in an aqueous environment with a physiological salt concentration enable us to perform a detailed analysis of how potassium ions interact with the double helix. The oligomers studied contain all 136 distinct tetranucleotides and we are thus able to make a comprehensive analysis of base sequence effects. Using a recently developed curvilinear helicoidal coordinate method we are able to analyze the details of ion populations and densities within the major and minor grooves and in the space surrounding DNA.

μABC: a systematic microsecond molecular dynamics study of tetranucleotide sequence effects in B-DNA.

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.

Twist and stretch of helices explained via the Kirchhoff-Love rod model of elastic filaments.

In various single-molecule experiments, a chiral polymer, such as DNA, is simultaneously pulled and twisted. We address an elementary but fundamental question raised by various authors: does the molecule overwind or unwind under tension? We show that within the context of the classic Kirchhoff-Love rod model of elastic filaments, both behaviors are possible, depending on the precise constitutive relations of the polymer. More generally, our analysis provides an effective linear response theory for helical structures that relates axial force and axial torque to axial translation and rotation.

Analyzing ion distributions around DNA.

We present a new method for analyzing ion, or molecule, distributions around helical nucleic acids and illustrate the approach by analyzing data derived from molecular dynamics simulations. The analysis is based on the use of curvilinear helicoidal coordinates and leads to highly localized ion densities compared to those obtained by simply superposing molecular dynamics snapshots in Cartesian space. The results identify highly populated and sequence-dependent regions where ions strongly interact with the nucleic and are coupled to its conformational fluctuations.

Looping probabilities of elastic chains: a path integral approach.

We consider an elastic chain at thermodynamic equilibrium with a heat bath, and derive an approximation to the probability density function, or pdf, governing the relative location and orientation of the two ends of the chain. Our motivation is to exploit continuum mechanics models for the computation of DNA looping probabilities, but here we focus on explaining the novel analytical aspects in the derivation of our approximation formula. Accordingly, and for simplicity, the current presentation is limited to the illustrative case of planar configurations.

A systematic molecular dynamics study of nearest-neighbor effects on base pair and base pair step conformations and fluctuations in B-DNA.

It is well recognized that base sequence exerts a significant influence on the properties of DNA and plays a significant role in protein-DNA interactions vital for cellular processes. Understanding and predicting base sequence effects requires an extensive structural and dynamic dataset which is currently unavailable from experiment. A consortium of laboratories was consequently formed to obtain this information using molecular simulations.

Bending modes of DNA directly addressed by cryo-electron microscopy of DNA minicircles.

We use cryo-electron microscopy (cryo-EM) to study the 3D shapes of 94-bp-long DNA minicircles and address the question of whether cyclization of such short DNA molecules necessitates the formation of sharp, localized kinks in DNA or whether the necessary bending can be redistributed and accomplished within the limits of the elastic, standard model of DNA flexibility. By comparing the shapes of covalently closed, nicked and gapped DNA minicircles, we conclude that 94-bp-long covalently closed and nicked DNA minicircles do not show sharp kinks while gapped DNA molecules, containing very flexible single-stranded regions, do show sharp kinks. We corroborate the results of cryo-EM studies by using Bal31 nuclease to probe for the existence of kinks in 94-bp-long minicircles.

Structural motifs of biomolecules.

Biomolecular structures are assemblies of emergent anisotropic building modules such as uniaxial helices or biaxial strands. We provide an approach to understanding a marginally compact phase of matter that is occupied by proteins and DNA. This phase, which is in some respects analogous to the liquid crystal phase for chain molecules, stabilizes a range of shapes that can be obtained by sequence-independent interactions occurring intra- and intermolecularly between polymeric molecules.

Kinking occurs during molecular dynamics simulations of small DNA minicircles.

Recent experiments on minicircle formation suggest that a conformational mechanism other than smooth deformation may be playing a role in enhancing DNA flexibility. Both local base unpairing and kink formation have been suggested as possible explanations. Although kinks within isolated DNA were proposed 30 years ago, they have, until now, only been observed within DNA complexed with proteins.

3D reconstruction and comparison of shapes of DNA minicircles observed by cryo-electron microscopy.

We use cryo-electron microscopy to compare 3D shapes of 158 bp long DNA minicircles that differ only in the sequence within an 18 bp block containing either a TATA box or a catabolite activator protein binding site. We present a sorting algorithm that correlates the reconstructed shapes and groups them into distinct categories. We conclude that the presence of the TATA box sequence, which is believed to be easily bent, does not significantly affect the observed shapes.

Helices.

Helices are among the simplest shapes that are observed in the filamentary and molecular structures of nature. The local mechanical properties of such structures are often modeled by a uniform elastic potential energy dependent on bending and twist, which is what we term a rod model. Our first result is to complete the semi-inverse classification, initiated by Kirchhoff, of all infinite, helical equilibria of inextensible, unshearable uniform rods with elastic energies that are a general quadratic function of the flexures and twist.

3-D shape estimation of DNA molecules from stereo cryo-electron micro-graphs using a projection-steerable snake.

We introduce a three-dimensional (3-D) parametric active contour algorithm for the shape estimation of DNA molecules from stereo cryo-electron micrographs. We estimate the shape by matching the projections of a 3-D global shape model with the micrographs; we choose the global model as a 3-D filament with a B-spline skeleton and a specified radial profile. The active contour algorithm iteratively updates the B-spline coefficients, which requires us to evaluate the projections and match them with the micrographs at every iteration.

Measuring limits of telomere movement on nuclear envelope.

The dynamic behavior of the decondensed chromatin can be monitored by real-time fluorescence confocal microscopy. It can be observed that different chromosomal sites enjoy different degrees of freedom during a certain period, exploring larger or smaller portions of nuclear volume. Here we measure the accessible surface for two chromosomal sites (yeast telomeres Tel3R and Tel6R) that both exhibit strong preferential association with the nuclear membrane in galactose-containing media, but differ significantly in gene activity.

Molecular dynamics simulations of the 136 unique tetranucleotide sequences of DNA oligonucleotides. II: sequence context effects on the dynamical structures of the 10 unique dinucleotide steps.

Molecular dynamics (MD) simulations including water and counterions on B-DNA oligomers containing all 136 unique tetranucleotide basepair steps are reported. The objective is to obtain the calculated dynamical structure for at least two copies of each case, use the results to examine issues with regard to convergence and dynamical stability of MD on DNA, and determine the significance of sequence context effects on all unique dinucleotide steps. This information is essential to understand sequence effects on DNA structure and has implications on diverse problems in the structural biology of DNA.

Molecular dynamics simulations of the 136 unique tetranucleotide sequences of DNA oligonucleotides. I. Research design and results on d(CpG) steps.

We describe herein a computationally intensive project aimed at carrying out molecular dynamics (MD) simulations including water and counterions on B-DNA oligomers containing all 136 unique tetranucleotide base sequences. This initiative was undertaken by an international collaborative effort involving nine research groups, the "Ascona B-DNA Consortium" (ABC). Calculations were carried out on the 136 cases imbedded in 39 DNA oligomers with repeating tetranucleotide sequences, capped on both ends by GC pairs and each having a total length of 15 nucleotide pairs.