Agnes Noy

PARMBSC1: A refined force-field for DNA simulations

We present parmbsc1, a force field for DNA atomistic simulation, which has been parameterized from high-level quantum mechanical data and tested for nearly 100 systems (representing a total simulation time of ∼140 μs) covering most of DNA structural space. Parmbsc1 provides high-quality results in diverse systems. Parameters and trajectories are available at http://mmb.irbbarcelona.org/ParmBSC1/.

Theoretical methods for the simulation of nucleic acids

Different theoretical methods for the description of nucleic acid structures are reviewed. Firstly, we introduce the concept of classical force-field in the context of nucleic acid structures, discussing their accuracy. We then examine theoretical approaches to the description of nucleic acids based on: i) a rigid or quasi-rigid description of the molecule, ii) molecular mechanics optimization, and iii) molecular dynamics. Special emphasis is made ion current state of the art molecular dynamics simulations of nucleic acids structures.

Recent advances in the study of nucleic acid flexibility by molecular dynamics

The recent use of molecular dynamics (MD) simulations to study flexibility of nucleic acids has been reviewed from an analysis of the publications appearing in the past two years (from 2005 till date). Despite the existence of some unsolved problems in the methodologies, these years have been witness to major advances in the field. Based on a critical review of the most recent contributions, excitement exists on the expected evolution of the field in the next years.

Theoretical study of large conformational transitions in DNA: the B↔A conformational change in water and ethanol/water

We explore here the possibility of determining theoretically the free energy change associated with large conformational transitions in DNA, like the solvent-induced B⇔A conformational change. We find that a combination of targeted molecular dynamics (tMD) and the weighted histogram analysis method (WHAM) can be used to trace this transition in both water and ethanol/water mixture. The pathway of the transition in the A→B direction mirrors the B→A pathway, and is dominated by two processes that occur somewhat independently: local changes in sugar puckering and global rearrangements (particularly twist and roll) in the structure. The B→A transition is found to be a quasi-harmonic process, which follows closely the first spontaneous deformation mode of B-DNA, showing that a physiologically-relevant deformation is in coded in the flexibility pattern of DNA.

Small DNA circles as probes of DNA topology.

Small DNA circles can occur in Nature, for example as protein-constrained loops, and can be synthesized by a number of methods. Such small circles provide tractable systems for the study of the structure, thermodynamics and molecular dynamics of closed-circular DNA. In the present article, we review the occurrence and synthesis of small DNA circles, and examine their utility in studying the properties of DNA and DNA-protein interactions. In particular, we highlight the analysis of small circles using atomistic simulations.

Data Mining of Molecular Dynamics Trajectories of Nucleic Acids

Abstract Analysis, storage, and transfer of molecular dynamic trajectories are becoming the bottleneck of computer simulations. In this paper we discuss different approaches for data mining and data processing of huge trajectory files generated from molecular dynamic simulations of nucleic acids.

Long-range correlations in the mechanics of small DNA circles under topological stress revealed by multi-scale simulation

Abstract It is well established that gene regulation can be achieved through activator and repressor proteins that bind to DNA and switch particular genes on or off, and that complex metabolic networks determine the levels of transcription of a given gene at a given time. Using three complementary computational techniques to study the sequence-dependence of DNA denaturation within DNA minicircles, we have observed that whenever the ends of the DNA are constrained, information can be transferred over long distances directly by the transmission of mechanical stress through the DNA itself, without any requirement for external signalling factors. Our models combine atomistic molecular dynamics (MD) with coarse-grained simulations and statistical mechanical calculations to span three distinct spatial resolutions and timescale regimes. While they give a consensus view of the non-locality of sequence-dependent denaturation in highly bent and supercoiled DNA loops, each also reveals a unique aspect of long-range informational transfer that occurs as a result of restraining the DNA within the closed loop of the minicircles.

Protein/DNA interactions in complex DNA topologies: expect the unexpected

DNA supercoiling results in compacted DNA structures that can bring distal sites into close proximity. It also changes the local structure of the DNA, which can in turn influence the way it is recognised by drugs, other nucleic acids and proteins. Here, we discuss how DNA supercoiling and the formation of complex DNA topologies can affect the thermodynamics of DNA recognition. We then speculate on the implications for transcriptional control and the three-dimensional organisation of the genetic material, using examples from our own simulations and from the literature. We introduce and discuss the concept of coupling between the multiple length-scales associated with hierarchical nuclear structural organisation through DNA supercoiling and topology.DNA supercoiling results in compacted DNA structures that can bring distal sites into close proximity. It also changes the local structure of the DNA, which can in turn influence the way it is recognised by drugs, other nucleic acids and proteins. Here, we discuss how DNA supercoiling and the formation of complex DNA topologies can affect the thermodynamics of DNA recognition. We then speculate on the implications for transcriptional control and the three-dimensional organisation of the genetic material, using examples from our own simulations and from the literature. We introduce and discuss the concept of coupling between the multiple length-scales associated with hierarchical nuclear structural organisation through DNA supercoiling and topology.

Comparison of molecular contours for measuring writhe in atomistic supercoiled DNA.

DNA molecular center-lines designed from atomistic-resolution structures are compared for the evaluation of the writhe in supercoiled DNA using molecular dynamics simulations of two sets of minicircles with 260 and 336 base pairs. We present a new method called WrLINE that systematically filters out local (i.e., subhelical turn) irregularities using a sliding-window averaged over a single DNA turn and that provides a measure of DNA writhe that is suitable for comparing atomistic resolution data with those obtained from measurements of the global molecular shape. In contrast, the contour traced by the base-pair origins defined by the 3DNA program largely overestimates writhe due to the helical periodicity of DNA. Nonetheless, this local modulation of the molecular axis emerges as an internal mechanism for the DNA to confront superhelical stress, where the adjustment between low and high twist is coupled to a high and low local periodicity, respectively, mimicking the different base-stacking conformational space of A and B canonical DNA forms.

Interference between Triplex and Protein Binding to Distal Sites on Supercoiled DNA

We have explored the interdependence of the binding of a DNA triplex and a repressor protein to distal recognition sites on supercoiled DNA minicircles using MD simulations. We observe that the interaction between the two ligands through their influence on their DNA template is determined by a subtle interplay of DNA mechanics and electrostatics, that the changes in flexibility induced by ligand binding play an important role and that supercoiling can instigate additional ligand-DNA contacts that would not be possible in simple linear DNA sequences.