Coarse-Grain Mechanics of DNA: Part II From Electrons to Oligomers
CECAM-HQ-EPFL, Lausanne, Switzerland
DNA carries genetic information so that understanding processes such as transcription and replication are central to biology. In the simplest classic view, the DNA sequence carries information merely via codons of three bases which encode for an amino acid when transcribed. Recently, this view has broadened to encompass the idea that the detailed chemistry of the four possible bases also provide a sequence-dependent modulation of the physical properties of the DNA, such as nonzero intrinsic curvature, varying intrinsic twist, and variations in stiffness, and that these sequence-dependent variations of the physical properties of the DNA are actually key to the control and regulation of the biological processes involving the DNA. For example, recently Segal et al  proposed that Nature exploits the redundancy of the genetic code (64 codons but only 20 amino acids) to create a "second genetic code" via the sequence-dependent mechanical properties of DNA which control the preferred locations of nucleosomes along the DNA. And the sequence-dependent mechanical properties of DNA is believed to be significant in other biological systems too.
The workshop will be a venue to discuss subjects surrounding the sequence-dependent physical properties of DNA. It is the second in a series of such discussions including the CECAM workshop “Coarse-Grain Mechanics of DNA: Bases to Chromosomes” which was held in Lyon in 2010. The third, to be organized by the Netherlands CECAM node, is anticipated for 2012. This second workshop will be tightly focused on discussing the shorter, biologically pertinent length scales of DNA, i.e. a few hundreds of base-pairs down to single bases and below. We are particularly interested in encouraging discussions between experimentalists and practitioners of modelling of various types. Specifically the participants will include experts in simulations at divers levels of coarse graining, ranging from a) quantum simulations of nucleic acids, b) atomistic molecular dynamics simulations, c) rigid base and base-pair simulations, and d) sequence-dependent continuum mechanics and polymer simulations up to lengths of one or two hundred base pairs. This is the scale of a DNA double-helix persistence length, which is believed to be rather important in biology. This range of length scales will be matched with the participation of experimentalists, from crystallographers at the atomistic level, to practitioners of looping, cyclization and mini-circle experiments, which offer a direct probe of sequence-dependent DNA mechanical properties. On the simulation side one of the primary issues is how one can quantitatively pass information between models at multiple scales. And of course interaction with experimentalists is key to the verification of model parameters at the various scales.
We believe that the undoubted importance of DNA in biology is already a powerful justification for a detailed study of its sequence-dependent physical properties. The mechanical properties of DNA do vary significantly with sequence, and to determine definitively whether and how these modulations are exploited in biology, requires an improved and quantitative understanding of the variations. In addition it should be realised that DNA is increasingly being adopted as a construction material in non-biological nanoscience, due to the inherent self-assembly properties of the double helix. However in the design of artificial nanostructures made from DNA, the sequence-dependence of the duplex DNA mechanical properties has not as yet been exploited, in part due to a lack of quantification of these properties. Accordingly a secondary, minor theme of the workshop will be to discuss such non-biological uses of DNA.
STATE OF THE ART
Atomistic Molecular Dynamics Simulations of Short Oligomers
The year 2011 will be the tenth anniversary of the forming of the ABC (or Ascona B-DNA Consortium) collaboration. The ABC was born at an interdisciplinary workshop organized by John Maddocks. The idea was that by pooling computational resources between groups, a database of entirely consistent simulations of DNA oligomers covering "all possible" sequence-dependence could be constructed. The most recent ABC publication , which is based on a data set of nearly 5microseconds of simulation time of 39 different 18 base-pair oligomers, is only now starting to approach that goal, and of course much was learnt along the way. For example, previous ABC data sets revealed serious limitations of the underlying Amber atomistic force field, which lead to a re-parametrization from quantum simulations . The current data set has no known inconsistencies with experimental data, and has various specific predictions, for example bimodal distribution of twist at certain specific base-pair steps.
Atomistic Molecular Dynamics Simulations of DNA Minicircles and Experiment:
The study of J-factors, or the probability of formation of DNA minicircles, is a classic technique for studying DNA mechanics. And minicircle formation has even been proposed as a way to determine experimentally the sequence-dependent mechanical properties of DNA . In particular some specific sequences were observed to cyclize remarkably efficiently compared to predictions made by standard, sequence-independent theories , although there is still some controversy over the experimental results [6,7]. Motivated by these results there have been two published molecular dynamics simulations of full minicircles [8,9], the first of which observed localized kink deformations of the DNA duplex, remarkably similar to those discussed theoretically thirty years previously . The 3D conformations of slightly longer DNA minicircles can be observed using stereo-cryo-EM [11, 12] although the method is very close to the limits of resolution. Nevertheless after the introduction of new statistical tests for curves in 3D, a statistically significant difference in the fluctuations of two minicircles formed from two different sequences could be observed .
Rigid base and base-pair models:
A standard model of the elasticity of double-helical DNA on the few nm length scale is the rigid base-pair model (rbpm), whose conformation variables are the relative positions and orientations of adjacent base pairs . Corresponding sequence-dependent local elastic potentials have been obtained from (combinations of) all-atom MD simulation and from high-resolution crystal structure data [15,16,17,18] and have been used with some success to model DNA-protein binding  and nucleosome positioning [18, 19]. Nevertheless the analysis of  has shown that at length scales of a few tens of base-pairs, a rbpm with a localized energy is not consistent with molecular dynamics data sets of the type obtained in , while a rigid base model with localized interactions is consistent with molecular dynamics data.
Sequence-dependent continuum mechanics models:
Beyond the 100nm scale, DNA is successfully described by a worm like chain model with homogeneous elastic properties , which can be determined via a systematic coarse-graining procedure from the rbpm . Modifications to model kinking have also been proposed . On the wormlike chain level a systematic treatment of small fluctuations (semi classical approximation) of the statistical mechanics of chains is now available [24,25].Sequence-dependent continuum models ignoring fluctuations have been applied to model cyclization data .
Ralf Everaers (École Normale Supérieure de Lyon) - Organiser
Helmut Schiessel (Leiden University) - Organiser
John H. Maddocks (Swiss Federal Institute of Technology Lausanne (EPFL)) - Organiser