cgDNA: a software package for the prediction of sequence-dependent coarse-grain free energies of B-form DNA

John Maddocks, Marco Pasi, Daiva Petkeviciute
Oxford Univ Press, 2014
Journal paper

Abstract

cgDNA is a package for the prediction of sequence-dependent configuration-space free energies for B-form DNA at the coarse-grain level of rigid bases. For a fragment of any given length and sequence, cgDNA calculates the configuration of the associated free energy minimizer, i.e. the relative positions and orientations of each base, along with a stiffness matrix, which together govern differences in free energies. The model predicts non-local (i.e. beyond base-pair step) sequence dependence of the free energy minimizer. Configurations can be input or output in either the Curves+definition of the usual helical DNA structural variables, or as a PDB file of coordinates of base atoms. We illustrate the cgDNA package by comparing predictions of free energy minimizers from (a) the cgDNA model, (b) time-averaged atomistic molecular dynamics (or MD) simulations, and (c) NMR or Xray experimental observation, for (i) the Dickerson-Drew dodecamer and (ii) three oligomers containing A-tracts. The cgDNA predictions are rather close to those of the MD simulations, but many orders of magnitude faster to compute. Both the cgDNA and MD predictions are in reasonable agreement with the available experimental data. Our conclusion is that cgDNA can serve as a highly efficient tool for studying structural variations in B-form DNA over a wide range of sequences.

Official source

https://infoscience.epfl.ch/record/205664?ln=en

About this result

This page is automatically generated and may contain information that is not correct, complete, up-to-date, or relevant to your search query. The same applies to every other page on this website. Please make sure to verify the information with EPFL's official sources.

Related concepts

Related publications

A DNA Coarse-Grain Rigid Base Model and Parameter Estimation from Molecular Dynamics Simulations

Daiva Petkeviciute

Sequence dependent mechanics of DNA is believed to play a central role in the functioning of the cell through the expression of genetic information. Nucleosome positioning, gene regulation, DNA looping and packaging within the cell are only some of the processes that are believed to be at least partially governed by mechanical laws. Therefore it is important to understand how the sequence of DNA affects its mechanical properties. For exploring the mechanical properties of DNA, various discrete and continuum models have been, and continue to be, developed. A large family of these models, including the model considered in this work, assume that bases or base pairs of DNA are rigid bodies. The most standard are rigid base pair models, with parameters either obtained directly from experimental data or from Molecular Dynamics (MD) simulations. The drawback of current experimental data, such as crystal structures, is that only small ensembles of configurations are available for a small number of sequences. In contrast, MD simulations allow a much more detailed view of a larger number of DNA sequences. However, the drawback is that the results of these simulations depend on the choice of the simulation protocol and force field parameters. MD simulations also have sequence length limitations and are currently too intensive for (linear) molecules longer than a few tens of base pairs. The only way to simulate longer sequences is to construct a coarse-grain model. The goal of this work is to construct a small parameter set that can model a sequence- dependent equilibrium probability distribution for rigid base configurations of a DNA oligomer with any given sequence of any length. The model parameter sets previously available were for rigid base pair models ignoring all the couplings beyond nearest neighbour interactions. However it was shown in previous work, that this standard model of rigid base pair nearest neighbour interactions is inconsistent with a (then) large scale MD simulation of a single oligomer [36]. In contrast we here show that a rigid base nearest neighbour, dimer sequence dependent model is a quite good fit to many MD simulations of different duration and se- quence. In fact a hierarchy of rigid base models with different interaction range and length of sequence-dependence is discussed, and it is concluded that the nearest neighbour, dimer based model is a good compromise between accuracy and complexity of the model. A full parameter set for this model is estimated. An interesting feature is that despite the dimer dependence of the parameter set, due to the phenomenon of frustration, our model predicts non local changes in the oligomer shape as a function of local changes in the sequence, down to the level of a point mutation.

EPFL2012

A sequence-dependent coarse-grain model of B-DNA with explicit description of bases and phosphate groups parametrised from large scale Molecular Dynamics simulations

Alessandro Samuele Patelli

We introduce a sequence-dependent coarse-grain model of double-stranded DNA with an explicit description of both the bases and the phosphate groups as interacting rigid-bodies. The model parameters are trained on extensive, state-of-the-art large scale molecular dynamics (MD) simulations. The model paradigm relies on three main approximations: 1) nucleic acid bases and phosphate groups are rigid, 2) interactions are nearest-neighbour and can be modelled with a quadratic energy, 3) model parameters have dimer sequence dependence. For an arbitrary sequence, the model predicts a sequence-dependent Gaussian equilibrium probability distribution. The parameter set comprises dimer-based elements, which are used to reconstruct mean configurations, called ground-states, which can have strong non-local sequence dependence, and precision matrices, or stiffness matrices, for any sequence of any length. This prediction step is sufficiently efficient that it is straightforward to construct probability density functions for millions of fragments each of length a few hundred base-pairs. The estimation of a parameter set consists in minimising the sum of Kullback-Leibler divergences between Gaussians predicted by the model and analogous Gaussians estimated directly from MD simulations of a training library of sequences. The training library comprises a short list of short palindromic DNA sequences. We designed the palindromic library using an ad hoc algorithm to include multiple instances of all independent tetramer sub-sequences. We exploit palindromic symmetry properties to study the convergence of the statistics extracted from MD simulations of palindromes and to define palindromically symmetrised estimators of first and second centred moments. The computation of the parameter set is delicate and needs the use of sophisticated numerics. We present an efficient and reliable procedure for estimating a complete parameter set which involves a generalisation of the classic Fisher information matrix and its relationship to the relative entropy, or Kullback-Leibler divergence. The model is a computationally efficient tool that allows the study of the mechanical properties of double-stranded DNA of arbitrary length and sequence. We use the model to study the sequence-dependent rigidity of DNA and we compute sequence-dependent apparent and dynamic persistence lengths. The explicit treatment of the phosphate group also allows computation of sequence-dependent grooves widths. Moreover, with fine-grained representation of predicted ground-states, we can also study sequence-dependence of sugar puckering modes and BI-BII backbone conformations.

EPFL2019

The interaction between neutral particles and the plasma plays a key role in determining the dynamics of the tokamak boundary that, in turn, significantly impact the overall performance of the device. Leveraging the work in Wersal and Ricci [Nucl. Fusion 55, 123014 (2015)], the present thesis describes the development and implementation in the GBS code of a mass-conserving multi-component self-consistent model to simulate the interplay between neutrals and plasma in the tokamak boundary. The simulation results are shown to be a useful tool to disentangle the physics at play in the tokamak boundary, and in particular for the interpretation of gas puff imaging (GPI).Developed in the last decade, the GBS code [Ricci et al, Plasma Phys. Control. Fusion 54, 124047 (2012)] allows for self-consistent three-dimensional numerical simulations of the turbulent plasma and neutral dynamics in the tokamak boundary. In GBS, a set of Braginskii equations in the drift limit describes the plasma time evolution, and the neutrals are modelled by solving a kinetic equation using the method of characteristics. While GBS enables simulations in arbitrary magnetic configurations, here we focus on limited plasmas. We first describe the geometrical operators and proper boundary conditions to ensure that mass conservation is satisfied. In comparison to the previous non-mass-conserving model of GBS, the mass-conserving simulations capture more accurately the sharp transition of the plasma and neutral quantities between the edge and scrape-off layer regions. In addition, we show that mass conserving simulations allow for reliable quantitative studies of particle fluxes in the tokamak boundary.A multi-component model of the neutral-plasma interaction is then developed by extending the single-component model to the description of a deuterium plasma that includes electrons, D+ ions, D atoms, D2 molecules and D2+ ions. The molecular dynamics is introduced through a set of drift-reduced Braginskii equations for the D2+ species and considering a kinetic equation for D2 molecules, in addition to the kinetic equation for D atoms, thus resulting in a coupled system of kinetic equations for the atomic and molecular neutral distribution functions. The first multi-component GBS simulations show that, in the sheath-limited regime under consideration, most of the D2 molecules cross the last closed flux surface (LCFS) and are dissociated or ionized in the edge region, thus giving rise to sources of D atoms inside the LCFS. This leads to an inward radial shift of the peak of the plasma source due to ionization of D atoms with respect to the single-component simulations.The multi-component model is applied to the simulation of GPI diagnostics, where the presence of a molecular gas puff is simulated self-consistently with the plasma and neutral dynamics of the tokamak boundary. The injected molecules interact with the boundary plasma, resulting in the emission of light in the D alpha wavelength that can be measured to infer the turbulent properties of the plasma. The simulated mechanisms underlying the light emission, which include the excitation of D atoms and dissociation of both D2 and D2+, provide a reliable tool for the interpretation of GPI experimental measurements. The impact of neutral fluctuations on the D alpha emission rate is investigated, as well as the correlation between the D alpha emission and the plasma and neutral quantities.