Fluidification of Entanglements by a DNA Bending Protein.

In spite of the nanoscale and single-molecule insights into nucleoid associated proteins (NAPs), their role in modulating the mesoscale viscoelasticity of entangled DNA has been overlooked so far. By combining microrheology and molecular dynamics simulation, we find that the abundant NAP "integration host factor" (IHF) lowers the viscosity of entangled λDNA 20-fold at physiological concentrations and stoichiometries. Our results suggest that IHF may play a previously unappreciated role in resolving DNA entanglements and in turn may be acting as a "genomic fluidizer" for bacterial genomes.

Integration host factor bends and bridges DNA in a multiplicity of binding modes with varying specificity.

Nucleoid-associated proteins (NAPs) are crucial in organizing prokaryotic DNA and regulating genes. Vital to these activities are complex nucleoprotein structures, however, how these form remains unclear. Integration host factor (IHF) is an Escherichia coli NAP that creates very sharp bends in DNA at sequences relevant to several functions including transcription and recombination, and is also responsible for general DNA compaction when bound non-specifically. We show that IHF-DNA structural multimodality is more elaborate than previously thought, and provide insights into how this drives mechanical switching towards strongly bent DNA. Using single-molecule atomic force microscopy and atomic molecular dynamics simulations we find three binding modes in roughly equal proportions: 'associated' (73° of DNA bend), 'half-wrapped' (107°) and 'fully-wrapped' (147°), only the latter occurring with sequence specificity. We show IHF bridges two DNA double helices through non-specific recognition that gives IHF a stoichiometry greater than one and enables DNA mesh assembly. We observe that IHF-DNA structural multiplicity is driven through non-specific electrostatic interactions that we anticipate to be a general NAP feature for physical organization of chromosomes.

The emergence of sequence-dependent structural motifs in stretched, torsionally constrained DNA.

The double-helical structure of DNA results from canonical base pairing and stacking interactions. However, variations from steady-state conformations resulting from mechanical perturbations in cells have physiological relevance but their dependence on sequence remains unclear. Here, we use molecular dynamics simulations showing sequence differences result in markedly different structural motifs upon physiological twisting and stretching. We simulate overextension on different sequences of DNA ((AA)12, (AT)12, (CC)12 and (CG)12) with supercoiling densities at 200 and 50 mM salt concentrations. We find that DNA denatures in the majority of stretching simulations, surprisingly including those with over-twisted DNA. GC-rich sequences are observed to be more stable than AT-rich ones, with the specific response dependent on the base pair order. Furthermore, we find that (AT)12 forms stable periodic structures with non-canonical hydrogen bonds in some regions and non-canonical stacking in others, whereas (CG)12 forms a stacking motif of four base pairs independent of supercoiling density. Our results demonstrate that 20-30% DNA extension is sufficient for breaking B-DNA around and significantly above cellular supercoiling, and that the DNA sequence is crucial for understanding structural changes under mechanical stress. Our findings have important implications for the activities of protein machinery interacting with DNA in all cells.

SerraNA: a program to determine nucleic acids elasticity from simulation data.

The resistance of DNA to stretch, twist and bend is broadly well estimated by experiments and is important for gene regulation and chromosome packing. However, their sequence-dependence and how bulk elastic constants emerge from local fluctuations is less understood. Here, we present SerraNA, which is an open software that calculates elastic parameters of double-stranded nucleic acids from dinucleotide length up to the whole molecule using ensembles from numerical simulations. The program reveals that global bendability emerge from local periodic bending angles in phase with the DNA helicoidal shape. We apply SerraNA to the whole set of 136 tetra-bp combinations and we observe a high degree of sequence-dependence with differences over 200% for all elastic parameters. Tetramers with TA and CA base-pair steps are especially flexible, while the ones containing AA and AT tend to be the most rigid. Thus, AT-rich motifs can generate extreme mechanical properties, which are critical for creating strong global bends when phased properly. Our results also indicate base mismatches would make DNA more flexible, while protein binding would make it more rigid. SerraNA is a tool to be applied in the next generation of interdisciplinary investigations to further understand what determines the elasticity of DNA.

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/.

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

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.

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.

Correlating fluorescence microscopy, optical and magnetic tweezers to study single chiral biopolymers such as DNA.

Biopolymer topology is critical for determining interactions inside cell environments, exemplified by DNA where its response to mechanical perturbation is as important as biochemical properties to its cellular roles. The dynamic structures of chiral biopolymers exhibit complex dependence with extension and torsion, however the physical mechanisms underpinning the emergence of structural motifs upon physiological twisting and stretching are poorly understood due to technological limitations in correlating force, torque and spatial localization information. We present COMBI-Tweez (Combined Optical and Magnetic BIomolecule TWEEZers), a transformative tool that overcomes these challenges by integrating optical trapping, time-resolved electromagnetic tweezers, and fluorescence microscopy, demonstrated on single DNA molecules, that can controllably form and visualise higher order structural motifs including plectonemes. This technology combined with cutting-edge MD simulations provides quantitative insight into complex dynamic structures relevant to DNA cellular processes and can be adapted to study a range of filamentous biopolymers.

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.

Length scale dependence of DNA mechanical properties.

Although mechanical properties of DNA are well characterized at the kilobase-pair range, a number of recent experiments have suggested that DNA is more flexible at shorter length scales, which correspond to the regime that is crucial for cellular processes such as DNA packaging and gene regulation. Here, we perform a systematic study of the effective elastic properties of DNA at different length scales by probing the conformation and fluctuations of DNA from the single base-pair level up to four helical turns, using trajectories from atomistic simulation. We find evidence that supports cooperative softening of the stretch modulus and identify the essential modes that give rise to this effect. The bend correlation exhibits modulations that reflect the helical periodicity, while it yields a reasonable value for the effective persistence length, and the twist modulus undergoes a smooth crossover--from a relatively smaller value at the single base-pair level to the bulk value--over half a DNA turn.

Atomistic Molecular Dynamics Simulations of DNA in Complex 3D Arrangements for Comparison with Lower Resolution Structural Experiments.

Atomic-level computer simulations are a very useful tool for describing the structure and dynamics of complex biomolecules such as DNA and for providing detail at a resolution where experimental techniques cannot arrive. Molecular dynamics (MD) simulations of mechanically distorted DNA caused by agents like supercoiling and protein binding are computationally challenging due to the large size of the associated systems and timescales. However, nowadays they are achievable thanks to the efficient usage of GPU and to the improvements of continuum solvation models. This together with the concurrent improvements in the resolution of single-molecule experiments, such as atomic force microscopy (AFM), makes possible the convergence between the two. Here we present detailed protocols for doing so: for performing molecular dynamics (MD) simulations of DNA adopting complex three-dimensional arrangements and for comparing the outcome of the calculations with single-molecule experimental data with a lower resolution than atomic.

Rolling circle RNA synthesis catalyzed by RNA.

RNA-catalyzed RNA replication is widely considered a key step in the emergence of life's first genetic system. However, RNA replication can be impeded by the extraordinary stability of duplex RNA products, which must be dissociated for re-initiation of the next replication cycle. Here, we have explored rolling circle synthesis (RCS) as a potential solution to this strand separation problem. We observe sustained RCS by a triplet polymerase ribozyme beyond full-length circle synthesis with strand displacement yielding concatemeric RNA products. Furthermore, we show RCS of a circular Hammerhead ribozyme capable of self-cleavage and re-circularization. Thus, all steps of a viroid-like RNA replication pathway can be catalyzed by RNA alone. Finally, we explore potential RCS mechanisms by molecular dynamics simulations, which indicate a progressive build-up of conformational strain upon RCS with destabilization of nascent strand 5'- and 3'-ends. Our results have implications for the emergence of RNA replication and for understanding the potential of RNA to support complex genetic processes.

Many organisms today rely on a trio of molecules for their survival: DNA, to store their genetic information; proteins, to conduct the biological processes required for growth or replication; and RNA, to mainly act as an intermediary between DNA and proteins. Yet, how these inanimate molecules first came together to form a living system remains unclear. Circumstantial evidence suggests that the first lifeforms relied to a much greater exrtent on RNA to conduct all necessary biological processes. There is no trace of this 'RNA world' today, but molecular 'fossils' may exist in current biology. Viroids, for example, are agents which can infect and replicate inside plant cells. They are formed of nothing but a circular strand of RNA that serves not only as genetic storage but also as ribozymes (RNA-based enzymes). Viroids need proteins from the host plant to replicate, but scientists have been able to engineer ribozymes that can copy complex RNA strands. This suggests that viroid-like replication could be achieved using only RNA. Kristoffersen et al. put this idea to the test and showed that it is possible to use RNA enzymatic activity alone to carry out all the steps of a viroid-like copying mechanism. This process included copying a viroid-like RNA circle with RNA, followed by trimming the copy to the right size and reforming the circle. These two latter steps could be carried out by a ribozyme that could itself be encoded on the RNA circle. A computer simulation indicated that RNA synthesis on the circle caused increasing tension that could ease some of the barriers to replication. These results increase our understanding of how RNA copying by RNA could be possible. This may lead to developing molecular models of a primordial RNA-based replication, which could be used to investigate early genetic systems and may have potential applications in synthetic biology.

Structural interplay between DNA-shape protein recognition and supercoiling: The case of IHF.

The integration host factor (IHF) is a prominent example of indirect readout as it imposes one of the strongest bends on relaxed linear DNA. However, the relation between IHF and torsionally constrained DNA, as occurs physiologically, remains unclear. By using atomistic molecular dynamics simulations on DNA minicircles, we reveal, for the first time, the reciprocal influence between a DNA-bending protein and supercoiling. On one hand, the increased curvature of supercoiled DNA enhances wrapping around IHF making the final complex topologically dependent. On the other hand, IHF acts as a 'supercoiling relief' factor by compacting relaxed DNA loops and, when supercoiled, it pins the position of plectonemes in a unique and specific manner. In addition, IHF restrains under- or overtwisted DNA depending on whether the complex is formed in negatively or positively supercoiled DNA, becoming effectively a 'supercoiling buffer'. We finally provide evidence of DNA bridging by IHF and reveal that these bridges divide DNA into independent topological domains. We anticipate that the crosstalk detected here between the 'active' DNA and the multifaceted IHF could be common to other DNA-protein complexes relying on the deformation of DNA.

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.

Elucidating the Role of Topological Constraint on the Structure of Overstretched DNA Using Fluorescence Polarization Microscopy.

The combination of DNA force spectroscopy and polarization microscopy of fluorescent DNA intercalator dyes can provide valuable insights into the structure of DNA under tension. These techniques have previously been used to characterize S-DNA-an elongated DNA conformation that forms when DNA overstretches at forces ≥ 65 pN. In this way, it was deduced that the base pairs of S-DNA are highly inclined, relative to those in relaxed (B-form) DNA. However, it is unclear whether and how topological constraints on the DNA may influence the base-pair inclinations under tension. Here, we apply polarization microscopy to investigate the impact of DNA pulling geometry, torsional constraint, and negative supercoiling on the orientations of intercalated dyes during overstretching. In contrast to earlier predictions, the pulling geometry (namely, whether the DNA molecule is stretched via opposite strands or the same strand) is found to have little influence. However, torsional constraint leads to a substantial reduction in intercalator tilting in overstretched DNA, particularly in AT-rich sequences. Surprisingly, the extent of intercalator tilting is similarly reduced when the DNA molecule is negatively supercoiled up to a critical supercoiling density (corresponding to ∼70% reduction in the linking number). We attribute these observations to the presence of P-DNA (an overwound DNA conformation). Our results suggest that intercalated DNA preferentially flanks regions of P-DNA rather than those of S-DNA and also substantiate previous suggestions that P-DNA forms predominantly in AT-rich sequences.

Base-pair resolution analysis of the effect of supercoiling on DNA flexibility and major groove recognition by triplex-forming oligonucleotides.

In the cell, DNA is arranged into highly-organised and topologically-constrained (supercoiled) structures. It remains unclear how this supercoiling affects the detailed double-helical structure of DNA, largely because of limitations in spatial resolution of the available biophysical tools. Here, we overcome these limitations, by a combination of atomic force microscopy (AFM) and atomistic molecular dynamics (MD) simulations, to resolve structures of negatively-supercoiled DNA minicircles at base-pair resolution. We observe that negative superhelical stress induces local variation in the canonical B-form DNA structure by introducing kinks and defects that affect global minicircle structure and flexibility. We probe how these local and global conformational changes affect DNA interactions through the binding of triplex-forming oligonucleotides to DNA minicircles. We show that the energetics of triplex formation is governed by a delicate balance between electrostatics and bonding interactions. Our results provide mechanistic insight into how DNA supercoiling can affect molecular recognition, that may have broader implications for DNA interactions with other molecular species.

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. In this sense, the knowledge obtained from the sequence-dependent structural properties of B-DNA has always been thought crucial to rationalize how ligands and, most notably, proteins recognize B-DNA and modulate its activity, i.e. the structural basis of gene regulation. Honouring the anniversary of the first high-resolution X-ray structure of a B-DNA molecule, in this contribution, we present the most important discoveries of the last 40 years on the sequence-dependent structural and dynamical properties of B-DNA, from the early beginnings to the current frontiers in the field.

Diversification of DNA-Binding Specificity by Permissive and Specificity-Switching Mutations in the ParB/Noc Protein Family.

Specific interactions between proteins and DNA are essential to many biological processes. Yet, it remains unclear how the diversification in DNA-binding specificity was brought about, and the mutational paths that led to changes in specificity are unknown. Using a pair of evolutionarily related DNA-binding proteins, each with a different DNA preference (ParB [Partitioning Protein B] and Noc [Nucleoid Occlusion Factor], which both play roles in bacterial chromosome maintenance), we show that specificity is encoded by a set of four residues at the protein-DNA interface. Combining X-ray crystallography and deep mutational scanning of the interface, we suggest that permissive mutations must be introduced before specificity-switching mutations to reprogram specificity and that mutational paths to new specificity do not necessarily involve dual-specificity intermediates. Overall, our results provide insight into the possible evolutionary history of ParB and Noc and, in a broader context, might be useful for understanding the evolution of other classes of DNA-binding proteins.

The impact of monovalent ion force field model in nucleic acids simulations.

Different classical models for monovalent ions (typically used to neutralize proteins or nucleic acids) are available in the literature and are widely used in molecular dynamics simulations without a great knowledge of their quality, consistency with the macromolecular force field and impact on the global simulation results. In this paper the ability of several of the most popular ion models to reproduce both quantum mechanics and experimental results is examined. Artefacts due to the use of incorrect ion models in molecular dynamics simulations of concentrated solutions of NaCl and KCl in water and of a short DNA duplex in 500 mM aqueous solutions of NaCl and KCl have been analyzed. Our results allow us to discuss the robustness and reliability of different ion models and to highlight the source of potential errors arising from non-optimal models. However, it is also found that the structural and dynamic characteristics of DNA (as an example of a heavily charged macromolecule) in near-physiological conditions are quite independent of the ion model used, providing support to most already-published simulations of macromolecules.

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.