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.

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.

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.

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.

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.

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.

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.

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

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.

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.

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.

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.

Atomistic Molecular Dynamics Simulations of DNA Minicircle Topoisomers: A Practical Guide to Setup, Performance, and Analysis.

While DNA supercoiling is ubiquitous in vivo, the structure of supercoiled DNA is more challenging to study experimentally than simple linear sequences because the DNA must have a closed topology in order to sustain superhelical stress. DNA minicircles, which are closed circular double-stranded DNA sequences typically containing between 60 and 500 base pairs, have proven to be useful biochemical tools for the study of supercoiled DNA mechanics. We present detailed protocols for constructing models of DNA minicircles in silico, for performing atomistic molecular dynamics (MD) simulations of supercoiled minicircle DNA, and for analyzing the results of the calculations. These simulations are computationally challenging due to the large system sizes. However, improvements in parallel computing software and hardware promise access to improve conformational sampling and simulation timescales. Given the concurrent improvements in the resolution of experimental techniques such as atomic force microscopy (AFM) and cryo-electron microscopy, the study of DNA minicircles will provide a more complete understanding of both the structure and the mechanics of supercoiled 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/.

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.