In the absence of writhe, DNA relieves torsional stress with localized, sequence-dependent structural failure to preserve B-form

Structural basis of DNA crossover capture by Escherichia coli DNA gyrase

DNA supercoiling must be precisely regulated by topoisomerases to prevent DNA entanglement. The interaction of type IIA DNA topoisomerases with two DNA molecules, enabling the transport of one duplex through the transient double-stranded break of the other, remains elusive owing to structures derived solely from single linear duplex DNAs lacking topological constraints. Using cryo-electron microscopy, we solved the structure of Escherichia coli DNA gyrase bound to a negatively supercoiled minicircle DNA. We show how DNA gyrase captures a DNA crossover, revealing both conserved molecular grooves that accommodate the DNA helices. Together with molecular tweezer experiments, the structure shows that the DNA crossover is of positive chirality, reconciling the binding step of gyrase-mediated DNA relaxation and supercoiling in a single structure.

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.

DNA opening during transcription initiation by RNA polymerase II in atomic detail

RNA polymerase II (RNAP II) synthesizes RNA by reading the DNA code. During transcription initiation, RNAP II opens the double-stranded DNA to expose the DNA template to the active site. The molecular interactions driving and controlling DNA opening are not well understood. We used all-atom steered molecular dynamics simulations to derive a continuous pathway of DNA opening in human RNAP II, involving a 55 Å DNA strand displacement and a nearly 360° DNA helix rotation. To drive such large-scale transitions, we used a combination of RMSD-based collective variables, a newly designed rotational coordinate, and a path collective variable. The simulations reveal extensive interactions of the DNA with three conserved protein loops near the active site, namely with the rudder, fork loop 1, and fork loop 2. According to the simulations, DNA-protein interactions support DNA opening by a twofold mechanism; they catalyze DNA opening by attacking Watson-Crick hydrogen bonds, and they stabilize the open DNA bubble by the formation of a wide set of DNA-protein salt bridges.

The End Restraint Method for Mechanically Perturbing Nucleic Acids In Silico

Far from being a passive information store, the genome is a mechanically dynamic and diverse system in which torsion and tension fluctuate and combine to determine structure and help regulate gene expression. Much of this mechanical perturbation is due to molecular machines such as topoisomerases which must stretch and twist DNA as part of various functions including DNA repair and replication. While the broad-scale mechanical response of nucleic acids to tension and torsion is well characterized, detail at the single base pair level is beyond the limits of even super-resolution imaging. Here, we present a straightforward, flexible, and extensible umbrella-sampling protocol to twist and stretch nucleic acids in silico using the popular biomolecular simulation package Amber-though the principles we describe are applicable also to other packages such as GROMACS. We discuss how to set up the simulation system, decide force fields and solvation models, and equilibrate. We then introduce the torsionally constrained stretching protocol, and finally we present some analysis techniques we have used to characterize structural motif formation. Rather than defining forces or fictional pseudoatoms, we instead define a fixed translation of specified atoms between each umbrella-sampling step, which allows comparison with experiment without needing to estimate applied forces by simply using the fractional end-to-end displacement as a comparison metric. We hope that this easy-to-implement solution will be valuable for interrogating optical and magnetic tweezers data on nucleic acids at base pair resolution.

Sequence-dependent twist-bend coupling in DNA minicircles

Looping of double-stranded DNA molecules with 100-200 base pairs into minicircles, catenanes, and rotaxanes has been suggested as a potential tool for DNA nanotechnologies. However, sharp DNA bending into a minicircle with a diameter of several to ten nanometers occurs with alterations in the DNA helical structure and may lead to defective kink formation that hampers the use of DNA minicircles, catenanes, and rotaxanes in nanoscale DNA applications. Here, we investigated local variations of a helical twist in sharply bent DNA using microsecond-long all-atom molecular dynamics simulations of six different DNA minicircles, focusing on the sequence dependence of the coupling between DNA bending and its helical twist. Twist angles between consecutive base pairs were analyzed at different locations relative to the direction of DNA bending and, among 10 unique dinucleotide steps, we identified four dinucleotide steps with strong twist-bend coupling, the pyrimidine-purine dinucleotide steps of TA/TA, CG/CG, and CA/TG and the purine-purine dinucleotide step of GA/TC. This work suggests the sequence-dependent structural responses of DNA to strong mechanical deformation, providing new molecular-level insights into the structure and stability of sharply bent DNA minicircles for nanoscale applications.

Supercoiling and looping promote DNA base accessibility and coordination among distant sites

DNA in cells is supercoiled and constrained into loops and this supercoiling and looping influence every aspect of DNA activity. We show here that negative supercoiling transmits mechanical stress along the DNA backbone to disrupt base pairing at specific distant sites. Cooperativity among distant sites localizes certain sequences to superhelical apices. Base pair disruption allows sharp bending at superhelical apices, which facilitates DNA writhing to relieve torsional strain. The coupling of these processes may help prevent extensive denaturation associated with genomic instability. Our results provide a model for how DNA can form short loops, which are required for many essential processes, and how cells may use DNA loops to position nicks to facilitate repair. Furthermore, our results reveal a complex interplay between site-specific disruptions to base pairing and the 3-D conformation of DNA, which influences how genomes are stored, replicated, transcribed, repaired, and many other aspects of DNA activity.

Effects of Model Shape, Volume, and Softness of the Capsid for DNA Packaging of phi29

Double-stranded DNA is under extreme confinement when packed in phage phi29 with osmotic pressures approaching 60 atm and densities near liquid crystalline. The shape of the capsid determined from experiment is elongated. We consider the effects of the capsid shape and volume on the DNA distribution. We propose simple models for the capsid of phage phi29 to capture volume, shape, and wall flexibility, leading to an accurate DNA density profile. The effect of the packaging motor twisting the DNA on the resulting density distribution has been explored. We find packing motor induced twisting leads to a greater numbers of defects formed. The emergence of defects such as bubbles or large roll angles along the DNA shows a sequence dependence, and the resulting flexibility leads to an inhomogeneous distribution of defects occurring more often at TpA steps and AT-rich regions. In conjunction with capsid elongation, this has effects on the global DNA packing structures.

Specifically bound BZIP transcription factors modulate DNA supercoiling transitions

Torsional stress on DNA, introduced by molecular motors, constitutes an important regulatory mechanism of transcriptional control. Torsional stress can modulate specific binding of transcription factors to DNA and introduce local conformational changes that facilitate the opening of promoters and nucleosome remodelling. Using all-atom microsecond scale molecular dynamics simulations together with a torsional restraint that controls the total twist of a DNA fragment, we address the impact of torsional stress on DNA complexation with a human BZIP transcription factor, MafB. We gradually over- and underwind DNA alone and in complex with MafB by 0.5° per dinucleotide step, starting from the relaxed state to a maximum of 5° per dinucleotide step, monitoring the evolution of the protein-DNA contacts at different degrees of torsional strain. Our computations show that MafB changes the DNA sequence-specific response to torsional stress. The dinucleotide steps that are susceptible to absorbing most of the torsional stress become more torsionally rigid, as they are involved in protein-DNA contacts. Also, the protein undergoes substantial conformational changes to follow the stress-induced DNA deformation, but mostly maintains the specific contacts with DNA. This results in a significant asymmetric increase of free energy of DNA twisting transitions, relative to free DNA, where overtwisting is more energetically unfavourable. Our data suggest that specifically bound BZIP factors could act as torsional stress insulators, modulating the propagation of torsional stress along the chromatin fibre, which might promote cooperative binding of collaborative DNA-binding factors.

How global DNA unwinding causes non-uniform stress distribution and melting of DNA

DNA unwinding is an important process that controls binding of proteins, gene expression and melting of double-stranded DNA. In a series of all-atom MD simulations on two DNA molecules containing a transcription start TATA-box sequence we demonstrate that application of a global restraint on the DNA twisting dramatically changes the coupling between helical parameters and the distribution of deformation energy along the sequence. Whereas only short range nearest-neighbor coupling is observed in the relaxed case, long-range coupling is induced in the globally restrained case. With increased overall unwinding the elastic deformation energy is strongly non-uniformly distributed resulting ultimately in a local melting transition of only the TATA box segment during the simulations. The deformation energy tends to be stored more in cytidine/guanine rich regions associated with a change in conformational substate distribution. Upon TATA box melting the deformation energy is largely absorbed by the melting bubble with the rest of the sequences relaxing back to near B-form. The simulations allow us to characterize the structural changes and the propagation of the elastic energy but also to calculate the associated free energy change upon DNA unwinding up to DNA melting. Finally, we design an Ising model for predicting the local melting transition based on empirical parameters. The direct comparison with the atomistic MD simulations indicates a remarkably good agreement for the predicted necessary torsional stress to induce a melting transition, for the position and length of the melted region and for the calculated associated free energy change between both approaches.

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.

BIOPHYSICS MEETS GENE THERAPY: HOW EXPLORING SUPERCOILING-DEPENDENT STRUCTURAL CHANGES IN DNA LED TO THE DEVELOPMENT OF MINIVECTOR DNA

Supercoiling affects every aspect of DNA function (replication, transcription, repair, recombination, etc.), yet the vast majority of studies on DNA and crystal structures of the molecule utilize short linear duplex DNA, which cannot be supercoiled. To study how supercoiling drives DNA biology, we developed and patented methods to make milligram quantities of tiny supercoiled circles of DNA called minicircles. We used a collaborative and multidisciplinary approach, including computational simulations (both atomistic and coarse-grained), biochemical experimentation, and biophysical methods to study these minicircles. By determining the three-dimensional conformations of individual supercoiled DNA minicircles, we revealed the structural diversity of supercoiled DNA and its highly dynamic nature. We uncovered profound structural changes, including sequence-specific base-flipping (where the DNA base flips out into the solvent), bending, and denaturing in negatively supercoiled minicircles. Counterintuitively, exposed DNA bases emerged in the positively supercoiled minicircles, which may result from inside-out DNA (Pauling-like, or "P-DNA"). These structural changes strongly influence how enzymes interact with or act on DNA. We hypothesized that, because of their small size and lack of bacterial sequences, these small supercoiled DNA circles may be efficient at delivering DNA into cells for gene therapy applications. "Minivectors," as we named them for this application, have proven to have therapeutic potential. We discovered that minivectors efficiently transfect a wide range of cell types, including many clinically important cell lines that are refractory to transfection with conventional plasmid vectors. Minivectors can be aerosolized for delivery to lungs and transfect human cells in culture to express RNA or genes. Importantly, minivectors demonstrate no obvious vector-associated toxicity. Minivectors can be repeatedly delivered and are long-lasting without integrating into the genome. Requests from colleagues around the world for minicircle and minivector DNA revealed a demand for our invention. We successfully obtained start-up funding for Twister Biotech, Inc. to help fulfill this demand, providing DNA for those who needed it, with a long-term goal of developing human therapeutics. In summary, what started as a tool for studying DNA structure has taken us in new and unanticipated directions.

Superstructure-Dependent Loading of DNA Origami Nanostructures with a Groove-Binding Drug

DNA origami nanostructures are regarded as powerful and versatile vehicles for targeted drug delivery. So far, DNA origami-based drug delivery strategies mostly use intercalation of the therapeutic molecules between the base pairs of the DNA origami's double helices for drug loading. The binding of nonintercalating drugs to DNA origami nanostructures, however, is less studied. Therefore, in this work, we investigate the interaction of the drug methylene blue (MB) with different DNA origami nanostructures under conditions that result in minor groove binding. We observe a noticeable effect of DNA origami superstructure on the binding affinity of MB. In particular, non-B topologies as for instance found in designs using the square lattice with 10.67 bp/turn may result in reduced binding affinity because groove binding efficiency depends on groove dimensions. Also, mechanically flexible DNA origami shapes that are prone to structural fluctuations may exhibit reduced groove binding, even though they are based on the honeycomb lattice with 10.5 bp/turn. This can be attributed to the induction of transient over- and underwound DNA topologies by thermal fluctuations. These issues should thus be considered when designing DNA origami nanostructures for drug delivery applications that employ groove-binding drugs.

Probing hyper-negatively supercoiled mini-circles with nucleases and DNA binding proteins

It is well accepted that the introduction of negative supercoils locally unwinds the DNA double helix, influencing thus the activity of proteins. Despite the use of recent methods of molecular dynamics simulations to model the DNA supercoiling-induced DNA deformation, the precise extent and location of unpaired bases induced by the negative supercoiling have never been investigated at the nucleotide level. Our goals in this study were to use radiolabeled double-stranded DNA mini-circles (dsMCs) to locate the unpaired bases on dsMCs whose topology ranged from relaxed to hyper-negatively supercoiled states, and to characterize the binding of proteins involved in the DNA metabolism. Our results show that the Nuclease SI is nearly ten times more active on hyper-negatively supercoiled than relaxed DNA. The structural changes responsible for this stimulation of activity were mapped for the first time with a base pair resolution and shown to be subtle and distributed along the entire sequence. As divalent cations modify the DNA topology, our binding studies were conducted with or without magnesium. Without magnesium, the dsMCs topoisomers mostly differ by their twist. Under these conditions, the Escherichia coli topoisomerase I weakly binds relaxed dsMCs and exhibits a stronger binding on negatively and hyper-negatively supercoiled dsMCs than relaxed dsMCs, with no significant difference in the binding activity among the supercoiled topoisomers. For the human replication protein A (hRPA), the more negatively supercoiled is the DNA, the better the binding, illustrating the twist-dependent binding activity for this protein. The presence of magnesium permits the dsMCs to writhe upon introduction of negative supercoiling and greatly modifies the binding properties of the hRPA and Escherichia coli SSB on dsMCs, indicating a magnesium-dependent DNA binding behavior. Finally, our experiments that probe the topology of the DNA in the hRPA-dsMC complexes show that naked and hRPA-bound dsMCs have the same topology.

Sequence-dependent response of DNA to torsional stress: a potential biological regulation mechanism

Torsional restraints on DNA change in time and space during the life of the cell and are an integral part of processes such as gene expression, DNA repair and packaging. The mechanical behavior of DNA under torsional stress has been studied on a mesoscopic scale, but little is known concerning its response at the level of individual base pairs and the effects of base pair composition. To answer this question, we have developed a geometrical restraint that can accurately control the total twist of a DNA segment during all-atom molecular dynamics simulations. By applying this restraint to four different DNA oligomers, we are able to show that DNA responds to both under- and overtwisting in a very heterogeneous manner. Certain base pair steps, in specific sequence environments, are able to absorb most of the torsional stress, leaving other steps close to their relaxed conformation. This heterogeneity also affects the local torsional modulus of DNA. These findings suggest that modifying torsional stress on DNA could act as a modulator for protein binding via the heterogeneous changes in local DNA structure.

Unwinding Induced Melting of Double-Stranded DNA Studied by Free Energy Simulations

DNA unwinding plays a major role in many biological processes, such as replication, transcription, and repair. It can lead to local melting and strand separation and can serve as a key mechanism to promote access to the separate strands of a double-stranded DNA. While DNA unwinding has been investigated extensively by DNA cyclization and single-molecule studies on a length-scale of kilo base pairs, it is neither fully understood at the base pair level nor at the level of molecular interactions. By employing a torque acting on the termini of DNA oligonucleotides during molecular dynamics free energy simulations, we locally unwind the central part of a DNA beyond an elastic (harmonic) regime. The simulations reproduce experimental results on the twist elasticity in the harmonic regime (characterized by a mostly quadratic free energy change with respect to changes in twist) and a deformation up to 7° was found as a limit of the harmonic response. Beyond this limit the free energy increase per twist change dropped dramatically coupled to local base pair disruptions and significant deformation of the nucleic acid backbone structure. Restriction of the DNA bending flexibility resulted in a stiffer harmonic response and an earlier onset of the anharmonic response. Whereas local melting with a complete disruption of base pairing and flipping of nucleotides was observed in case of an AT rich central segment strong backbone changes and changes in the stacking arrangements were observed in case of a GC rich segment. Unrestrained MD simulations starting from locally melted DNA reformed regular B-DNA after 50-300 ns simulation time. The simulations may have important implications for understanding DNA recognition processes coupled with significant structural alterations.

Influence of DNA sequence on the structure of minicircles under torsional stress

The sequence dependence of the conformational distribution of DNA under various levels of torsional stress is an important unsolved problem. Combining theory and coarse-grained simulations shows that the DNA sequence and a structural correlation due to topology constraints of a circle are the main factors that dictate the 3D structure of a 336 bp DNA minicircle under torsional stress. We found that DNA minicircle topoisomers can have multiple bend locations under high torsional stress and that the positions of these sharp bends are determined by the sequence, and by a positive mechanical correlation along the sequence. We showed that simulations and theory are able to provide sequence-specific information about individual DNA minicircles observed by cryo-electron tomography (cryo-ET). We provided a sequence-specific cryo-ET tomogram fitting of DNA minicircles, registering the sequence within the geometric features. Our results indicate that the conformational distribution of minicircles under torsional stress can be designed, which has important implications for using minicircle DNA for gene therapy.

Phage-like packing structures with mean field sequence dependence

Packing of double-stranded DNA in phages must overcome both electrostatic repulsions and the problem of persistence length. We consider coarse-grained models with the ability to kink and with randomly generated disorder. We show that the introduction of kinking into configurations of the DNA polymer packaged within spherical confinement results in significant reductions of the overall energies and pressures. We use a kink model which has the ability to deform every 24 bp, close to the average length predicted from phage sequence. The introduction of such persistence length defects even with highly random packing models increases the local nematic ordering of the packed DNA polymer segments. Such local ordering allowed by kinking not only reduces the total bending energy of confined DNA due to nonlinear elasticity but also reduces the electrostatic component of the energy and pressure. We show that a broad ensemble of polymer configurations is consistent with the structural data. © 2016 Wiley Periodicals, Inc.

DNA Shape versus Sequence Variations in the Protein Binding Process

The binding process of a protein with a DNA involves three stages: approach, encounter, and association. It has been known that the complexation of protein and DNA involves mutual conformational changes, especially for a specific sequence association. However, it is still unclear how the conformation and the information in the DNA sequences affects the binding process. What is the extent to which the DNA structure adopted in the complex is induced by protein binding, or is instead intrinsic to the DNA sequence? In this study, we used the multiscale simulation method to explore the binding process of a protein with DNA in terms of DNA sequence, conformation, and interactions. We found that in the approach stage the protein can bind both the major and minor groove of the DNA, but uses different features to locate the binding site. The intrinsic conformational properties of the DNA play a significant role in this binding stage. By comparing the specific DNA with the nonspecific in unbound, intermediate, and associated states, we found that for a specific DNA sequence, ∼40% of the bending in the association forms is intrinsic and that ∼60% is induced by the protein. The protein does not induce appreciable bending of nonspecific DNA. In addition, we proposed that the DNA shape variations induced by protein binding are required in the early stage of the binding process, so that the protein is able to approach, encounter, and form an intermediate at the correct site on DNA.

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.

Sequence Affects the Cyclization of DNA Minicircles

Understanding how the sequence of a DNA molecule affects its dynamic properties is a central problem affecting biochemistry and biotechnology. The process of cyclizing short DNA, as a critical step in molecular cloning, lacks a comprehensive picture of the kinetic process containing sequence information. We have elucidated this process by using coarse-grained simulations, enhanced sampling methods, and recent theoretical advances. We are able to identify the types and positions of structural defects during the looping process at a base-pair level. Correlations along a DNA molecule dictate critical sequence positions that can affect the looping rate. Structural defects change the bending elasticity of the DNA molecule from a harmonic to subharmonic potential with respect to bending angles. We explore the subelastic chain as a possible model in loop formation kinetics. A sequence-dependent model is developed to qualitatively predict the relative loop formation time as a function of DNA sequence.

Structural diversity of supercoiled DNA

By regulating access to the genetic code, DNA supercoiling strongly affects DNA metabolism. Despite its importance, however, much about supercoiled DNA (positively supercoiled DNA, in particular) remains unknown. Here we use electron cryo-tomography together with biochemical analyses to investigate structures of individual purified DNA minicircle topoisomers with defined degrees of supercoiling. Our results reveal that each topoisomer, negative or positive, adopts a unique and surprisingly wide distribution of three-dimensional conformations. Moreover, we uncover striking differences in how the topoisomers handle torsional stress. As negative supercoiling increases, bases are increasingly exposed. Beyond a sharp supercoiling threshold, we also detect exposed bases in positively supercoiled DNA. Molecular dynamics simulations independently confirm the conformational heterogeneity and provide atomistic insight into the flexibility of supercoiled DNA. Our integrated approach reveals the three-dimensional structures of DNA that are essential for its function.

DNA supercoiling strongly affects its metabolism. By electron cryo-tomography, biochemical assays and molecular dynamics simulations, here the authors show that supercoiled DNA minicircles adopt unique and wide distributions of three-dimensional conformations, many with disrupted base pairs.

Twist-induced defects of the P-SSP7 genome revealed by modeling the cryo-EM density

We consider the consequences of assuming that DNA inside of phages can be approximated as a strongly nonlinear persistence length polymer. Recent cryo-EM experiments find a hole in the density map of P-SSP7 phage, located in the DNA segment filling the portal channel of the phage. We use experimentally derived structural constraints with coarse-grained simulation techniques to consider contrasting model interpretations of reconstructed density in the portal channel. The coarse-grained DNA models used are designed to capture the effects of torsional strain and electrostatic environment. Our simulation results are consistent with the interpretation that the vacancy or hole in the experimental density map is due to DNA strain leading to strand separation. We further demonstrate that a moderate negative twisting strain is able to account for the strand separation. This effect of nonlinear persistence length may be important in other aspects of phage DNA packing.

Experimental phase diagram of negatively supercoiled DNA measured by magnetic tweezers and fluorescence

The most common form of DNA is the well-known B-structure of double-helix DNA. Many processes in the cell, however, exert force and torque, inducing structural changes to the DNA that are vital to biological function. Virtually all DNA in cells is in a state of negative supercoiling, with a DNA structure that is complex. Using magnetic tweezers combined with fluorescence imaging, we here study DNA structure as a function of negative supercoiling at the single-molecule level. We classify DNA phases based on DNA length as a function of supercoiling, down to a very high negative supercoiling density σ of -2.5, and forces up to 4.5 pN. We characterize plectonemes using fluorescence imaging. DNA bubbles are visualized by the binding of fluorescently labelled RPA, a eukaryotic single-strand-binding protein. The presence of Z-DNA, a left-handed form of DNA, is probed by the binding of Zα77, the minimal binding domain of a Z-DNA-binding protein. Without supercoiling, DNA is in the relaxed B-form. Upon going toward negative supercoiling, plectonemic B-DNA is being formed below 0.6 pN. At higher forces and supercoiling densities down to about -1.9, a mixed state occurs with plectonemes, multiple bubbles and left-handed L-DNA. Around σ = -1.9, a buckling transition occurs after which the DNA end-to-end length linearly decreases when applying more negative turns, into a state that we interpret as plectonemic L-DNA. By measuring DNA length, Zα77 binding, plectoneme and ssDNA visualisation, we thus have mapped the co-existence of many DNA structures and experimentally determined the DNA phase diagram at (extreme) negative supercoiling.

Modeling DNA thermodynamics under torsional stress

Negatively twisted DNA is essential to many biological functions. Due to torsional stress, duplex DNA can have local, sequence-dependent structural defects. In this work, a thermodynamic model of DNA was built to qualitatively predict the local sequence-dependent mechanical instabilities under torsional stress. The results were compared to both simulation of a coarse-grained model and experiment results. By using the Kirkwood superposition approximation, we built an analytical model to represent the free energy difference ΔW of a hydrogen-bonded basepair between the B-form helical state and the basepair opened (or locally melted) state, within a given sequence under torsional stress. We showed that ΔW can be well approximated by two-body interactions with its nearest-sequence-neighbor basepairs plus a free energy correction due to long-range correlations. This model is capable of rapidly predicting the position and thermodynamics of local defects in a given sequence. The result qualitatively matches with an in vitro experiment for a long DNA sequence (>4000 basepairs). The 12 parameters used in this model can be further quantitatively refined when more experimental data are available.

Close encounters with DNA

Over the past ten years, the all-atom molecular dynamics method has grown in the scale of both systems and processes amenable to it and in its ability to make quantitative predictions about the behavior of experimental systems. The field of computational DNA research is no exception, witnessing a dramatic increase in the size of systems simulated with atomic resolution, the duration of individual simulations and the realism of the simulation outcomes. In this topical review, we describe the hallmark physical properties of DNA from the perspective of all-atom simulations. We demonstrate the amazing ability of such simulations to reveal the microscopic physical origins of experimentally observed phenomena. We also discuss the frustrating limitations associated with imperfections of present atomic force fields and inadequate sampling. The review is focused on the following four physical properties of DNA: effective electric charge, response to an external mechanical force, interaction with other DNA molecules and behavior in an external electric field.

Triple helical DNA in a duplex context and base pair opening

It is fundamental to explore in atomic detail the behavior of DNA triple helices as a means to understand the role they might play in vivo and to better engineer their use in genetic technologies, such as antigene therapy. To this aim we have performed atomistic simulations of a purine-rich antiparallel triple helix stretch of 10 base triplets flanked by canonical Watson–Crick double helices. At the same time we have explored the thermodynamic behavior of a flipping Watson–Crick base pair in the context of the triple and double helix. The third strand can be accommodated in a B-like duplex conformation. Upon binding, the double helix changes shape, and becomes more rigid. The triple-helical region increases its major groove width mainly by oversliding in the negative direction. The resulting conformations are somewhere between the A and B conformations with base pairs remaining almost perpendicular to the helical axis. The neighboring duplex regions maintain a B DNA conformation. Base pair opening in the duplex regions is more probable than in the triplex and binding of the Hoogsteen strand does not influence base pair breathing in the neighboring duplex region.

Communication: Origin of the contributions to DNA structure in phages

Cryo electron microscopy (cryo-EM) data of the interior of phages show ordering of the interior DNA that has been interpreted as a nearly perfectly ordered polymer. We show surface-induced correlations, excluded volume, and electrostatic forces are sufficient to predict most of the major features of the current structural data for DNA packaged within viral capsids without additional ordering due to elastic bending forces for the polymer. Current models assume highly-ordered, even spooled, hexagonally packed conformations based on interpretation of cryo-EM density maps. We show herein that the surface induced packing of short (6mer), unconnected DNA polymer segments is the only necessary ingredient in creating ringed densities consistent with experimental density maps. This implies the ensemble of possible conformations of polymeric DNA within the capsid that are consistent with cryo-EM data may be much larger than implied by traditional interpretations where such rings can only result from highly-ordered spool-like conformations. This opens the possibility of a more disordered, entropically-driven view of phage packaging thermodynamics. We also show the electrostatics of the DNA contributes a large portion of the internal hydrostatic and osmotic pressures of a phage virion, suggesting that nonlinear elastic anomalies might reduce the overall elastic bending enthalpy of more disordered conformations to have allowable free energies.

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.

Ab initio determination of coarse-grained interactions in double-stranded DNA

We derive the coarse-grained interactions between DNA nucleotides from ab initio total-energy calculations based on density functional theory (DFT). The interactions take into account base and sequence specificity, and are decomposed into physically distinct contributions that include hydrogen bonding, stacking interactions, backbone, and backbone-base interactions. The interaction energies of each contribution are calculated from DFT for a wide range of configurations and are fitted by simple analytical expressions for use in the coarse-grained model, which reduces each nucleotide into two sites. This model is not derived from experimental data, yet it successfully reproduces the stable B-DNA structure and gives good predictions for the persistence length. It may be used to realistically probe dynamics of DNA strands in various environments at the μs time scale and the μm length scale.

CAG/CTG repeats alter the affinity for the histone core and the positioning of DNA in the nucleosome

Trinucleotide repeats (TNRs) occur throughout the genome, and their expansion has been linked to several neurodegenerative disorders, including Huntington's disease. TNRs have been studied using both oligonucleotides and plasmids; however, less is know about how repetitive DNA responds to genomic packaging. Here, we investigate the behavior of CAG/CTG repeats incorporated into nucleosome core particles, the most basic unit of chromatin packaging. To assess the general interaction between CAG/CTG repeats and the histone core, we determined the efficiency with which various TNR-containing DNA substrates form nucleosomes, revealing that even short CAG/CTG tracts are robust incorporators. However, the presence of the Huntington gene flanking sequence (htt) decreases the rate of incorporation. Enzymatic and chemical probing revealed repositioning of the DNA in the nucleosome as the number of CAG/CTG repeats increased, regardless of the flanking sequence. Notably, the periodicity of the repeat tract remained unchanged as a function of length and is consistently 10.7 bp per helical turn. In contrast, the periodicity of the nonrepetitive flanking sequence varies and is smaller than the repeat tract at ~10.0-10.5 bp per turn. Furthermore, while the CAG/CTG repeats remain as a canonical duplex in the nucleosome, nucleosome formation causes kinking in a secondary repeat tract in the htt gene, comprised of CCG/CGG repeats. This work highlights the innate ability of CAG/CTG repeats to incorporate and to position in nucleosomes and how that behavior is modulated by the htt flanking sequence. In addition, it illuminates the differences in packaging of healthy and diseased length repeat tracts within the genome.

DNA cruciform arms nucleate through a correlated but asynchronous cooperative mechanism

Inverted repeat (IR) sequences in DNA can form noncanonical cruciform structures to relieve torsional stress. We use Monte Carlo simulations of a recently developed coarse-grained model of DNA to demonstrate that the nucleation of a cruciform can proceed through a cooperative mechanism. First, a twist-induced denaturation bubble must diffuse so that its midpoint is near the center of symmetry of the IR sequence. Second, bubble fluctuations must be large enough to allow one of the arms to form a small number of hairpin bonds. Once the first arm is partially formed, the second arm can rapidly grow to a similar size. Because bubbles can twist back on themselves, they need considerably fewer bases to resolve torsional stress than the final cruciform state does. The initially stabilized cruciform therefore continues to grow, which typically proceeds synchronously, reminiscent of the S-type mechanism of cruciform formation. By using umbrella sampling techniques, we calculate, for different temperatures and superhelical densities, the free energy as a function of the number of bonds in each cruciform arm along the correlated but asynchronous nucleation pathways we observed in direct simulations.

Bullied no more: when and how DNA shoves proteins around

The predominant protein-centric perspective in protein-DNA-binding studies assumes that the protein drives the interaction. Research focuses on protein structural motifs, electrostatic surfaces and contact potentials, while DNA is often ignored as a passive polymer to be manipulated. Recent studies of DNA topology, the supercoiling, knotting, and linking of the helices, have shown that DNA has the capability to be an active participant in its transactions. DNA topology-induced structural and geometric changes can drive, or at least strongly influence, the interactions between protein and DNA. Deformations of the B-form structure arise from both the considerable elastic energy arising from supercoiling and from the electrostatic energy. Here, we discuss how these energies are harnessed for topology-driven, sequence-specific deformations that can allow DNA to direct its own metabolism.

Modelling Nucleic Acid Structure and Flexibility: From Atomic to Mesoscopic Scale

This chapter surveys some of the recent developments in coarse-grained modelling of nucleic acids. We first discuss models based on pseudoatoms, effective spherical particles representing groups of atoms. A major part of the chapter is devoted to models in which bases or base pairs are represented as independent, interacting rigid bodies. Two popular definitions of internal coordinates, as implemented in the programs 3DNA and Curves+, are outlined from a common perspective. Recently developed rigid base and basepair models with nonlocal quadratic interactions are presented. A statistical mechanical description of the models on their full phase space yields exact relations between model parameters and expected values of some state functions. We estimated shape and stiffness parameters for nonlocal rigid base and basepair models of a DNA oligomer containing A-tract. The parameterization is based on atomic-resolution molecular dynamics simulation data. We found that the rigid base model is consistent with a local interaction pattern, while interactions in the rigid basepair model are visibly non-local, in agreement with earlier findings. Differences in shape and stiffness parameters obtained using Curves+ and 3DNA coordinates are found to be small for structures within the B-DNA family. Anharmonic effects, coarser models, and other approaches to describe nucleic acid structure and flexibility are discussed.

Frontiers in molecular dynamics simulations of DNA

It has been known for decades that DNA is extremely flexible and polymorphic, but our knowledge of its accessible conformational space remains limited. Structural data, primarily from X-ray diffraction studies, is sparse in comparison to the manifold configurations possible, and direct experimental examinations of DNA's flexibility still suffer from many limitations. In the face of these shortcomings, molecular dynamics (MD) is now an essential tool in the study of DNA. It affords detailed structural and dynamical insights, which explains its recent transition from a small number of highly specialized laboratories to a large variety of groups dealing with challenging biological problems. MD is now making an irreversible journey to the mainstream of research in biology, with the attendant opportunities and challenges. But given the speed with which MD studies of DNA have spread, the roots remain somewhat shallow: in many cases, there is a lack of deep knowledge about the foundations, strengths, and limits of the technique. In this Account, we discuss how MD has become the most important source of structural and flexibility data on DNA, focusing on advances since 2007 of atomistic MD in the description of DNA under near-physiological conditions and highlighting the possibilities and shortcomings of the technique. The evolution in the field over the past four years is a prelude to the ongoing revolution. The technique has gained in robustness and predictive power, which when coupled with the spectacular improvements in software and hardware has enabled the tackling of systems of increasing complexity. Simulation times of microseconds have now been achieved, with even longer times when specialized hardware is used. As a result, we have seen the first real-time simulation of large conformational transitions, including folding and unfolding of short DNA duplexes. Noteworthy advances have also been made in the study of DNA-ligand interactions, and we predict that a global thermodynamic and kinetic picture of the binding landscape of DNA will become available in a few years. MD will become a crucial tool in areas such as biomolecular engineering and synthetic biology. MD has also been shown to be an excellent source of parameters for mesoscopic models of DNA flexibility. Such models can be refined through atomistic MD simulations on small duplexes and then applied to the study of entire chromosomes. Recent evidence suggests that MD-derived elastic models can successfully predict the position of regulatory regions in DNA and can help advance our understanding of nucleosome positioning and chromatin plasticity. If these results are confirmed, MD simulations can become the ultimate tool to decipher a physical code that can contribute to gene regulation. We are entering the golden age of MD simulations of DNA. Undoubtedly, the expectations are high, but the challenges are also enormous. These include the need for more accurate potential energy functionals and for longer and more complex simulations in more realistic systems. The joint research effort of several groups will be crucial for adapting the technique to the requirements of the coming decade.

Supercoiled Minivector DNA resists shear forces associated with gene therapy delivery

Supercoiled DNAs varying from 281 to 5302 bp were subjected to shear forces generated by aerosolization or sonication. DNA shearing strongly correlated with length. Typical sized plasmids (≥ 3000 bp) degraded rapidly. DNAs 2000-3000 bp persisted ~10 min. Even in the absence of condensing agents, supercoiled DNA <1200 bp survived nebulization, and increased forces of sonication were necessary to shear it. Circular vectors were considerably more resistant to shearing than linear vectors of the same length. DNA supercoiling afforded additional protection. These results show the potential of shear-resistant Minivector DNAs to overcome one of the major challenges associated with gene therapy delivery.

DNA movies and panspermia

There are several ways that our species might try to send a message to another species separated from us by space and/or time. Synthetic biology might be used to write an epitaph to our species, or simply “Kilroy was here”, in the genome of a bacterium via the patterns of either (1) the codons to exploit Life's non-equilibrium character or (2) the bases themselves to exploit Life's quasi-equilibrium character. We suggest here how DNA movies might be designed using such patterns. We also suggest that a search for mechanisms to create and preserve such patterns might lead to a better understanding of modern cells. Finally, we argue that the cutting-edge microbiology and synthetic biology needed for the Kilroy project would put origin-of-life studies in the vanguard of research.

Cooperative kinking at distant sites in mechanically stressed DNA

In cells, DNA is routinely subjected to significant levels of bending and twisting. In some cases, such as under physiological levels of supercoiling, DNA can be so highly strained, that it transitions into non-canonical structural conformations that are capable of relieving mechanical stress within the template. DNA minicircles offer a robust model system to study stress-induced DNA structures. Using DNA minicircles on the order of 100 bp in size, we have been able to control the bending and torsional stresses within a looped DNA construct. Through a combination of cryo-EM image reconstructions, Bal31 sensitivity assays and Brownian dynamics simulations, we have been able to analyze the effects of biologically relevant underwinding-induced kinks in DNA on the overall shape of DNA minicircles. Our results indicate that strongly underwound DNA minicircles, which mimic the physical behavior of small regulatory DNA loops, minimize their free energy by undergoing sequential, cooperative kinking at two sites that are located about 180° apart along the periphery of the minicircle. This novel form of structural cooperativity in DNA demonstrates that bending strain can localize hyperflexible kinks within the DNA template, which in turn reduces the energetic cost to tightly loop DNA.

Underwound DNA under tension: structure, elasticity, and sequence-dependent behaviors

DNA melting under torsion plays an important role in a wide variety of cellular processes. In the present Letter, we have investigated DNA melting at the single-molecule level using an angular optical trap. By directly measuring force, extension, torque, and angle of DNA, we determined the structural and elastic parameters of torsionally melted DNA. Our data reveal that under moderate forces, the melted DNA assumes a left-handed structure as opposed to an open bubble conformation and is highly torsionally compliant. We have also discovered that at low forces melted DNA properties are highly dependent on DNA sequence. These results provide a more comprehensive picture of the global DNA force-torque phase diagram.

Local elasticity of strained DNA studied by all-atom simulations

Genomic DNA is constantly subjected to various mechanical stresses arising from its biological functions and cell packaging. If the local mechanical properties of DNA change under torsional and tensional stress, the activity of DNA-modifying proteins and transcription factors can be affected and regulated allosterically. To check this possibility, appropriate steady forces and torques were applied in the course of all-atom molecular dynamics simulations of DNA with AT- and GC-alternating sequences. It is found that the stretching rigidity grows with tension as well as twisting. The torsional rigidity is not affected by stretching, but it varies with twisting very strongly, and differently for the two sequences. Surprisingly, for AT-alternating DNA it passes through a minimum with the average twist close to the experimental value in solution. For this fragment, but not for the GC-alternating sequence, the bending rigidity noticeably changes with both twisting and stretching. The results have important biological implications and shed light on earlier experimental observations.

Atomistic simulations reveal bubbles, kinks and wrinkles in supercoiled DNA

Although DNA is frequently bent and supercoiled in the cell, much of the available information on DNA structure at the atomistic level is restricted to short linear sequences. We report atomistic molecular dynamics (MD) simulations of a series of DNA minicircles containing between 65 and 110 bp which we compare with a recent biochemical study of structural distortions in these tight DNA loops. We have observed a wealth of non-canonical DNA structures such as kinks, denaturation bubbles and wrinkled conformations that form in response to bending and torsional stress. The simulations show that bending alone is sufficient to induce the formation of kinks in circles containing only 65 bp, but we did not observe any defects in simulations of larger torsionally relaxed circles containing 110 bp over the same MD timescales. We also observed that under-winding in minicircles ranging in size from 65 to 110 bp leads to the formation of single stranded bubbles and wrinkles. These calculations are used to assess the ability of atomistic MD simulations to determine the structure of bent and supercoiled DNA.

Distinguishing the roles of Topoisomerases I and II in relief of transcription-induced torsional stress in yeast rRNA genes

To better understand the role of topoisomerase activity in relieving transcription-induced supercoiling, yeast genes encoding rRNA were visualized in cells deficient for either or both of the two major topoisomerases. In the absence of both topoisomerase I (Top1) and topoisomerase II (Top2) activity, processivity was severely impaired and polymerases were unable to transcribe through the 6.7-kb gene. Loss of Top1 resulted in increased negative superhelical density (two to six times the normal value) in a significant subset of rRNA genes, as manifested by regions of DNA template melting. The observed DNA bubbles were not R-loops and did not block polymerase movement, since genes with DNA template melting showed no evidence of slowed elongation. Inactivation of Top2, however, resulted in characteristic signs of slowed elongation in rRNA genes, suggesting that Top2 alleviates transcription-induced positive supercoiling. Together, the data indicate that torsion in front of and behind transcribing polymerase I has different consequences and different resolution. Positive torsion in front of the polymerase induces supercoiling (writhe) and is largely resolved by Top2. Negative torsion behind the polymerase induces DNA strand separation and is largely resolved by Top1.

Direct observation of multiple pathways of single-stranded DNA stretching

We observed multiple pathways of stretching single-stranded polydeoxynucleotides, poly(dA). Poly(dA) has been shown to undergo unique transitions under mechanical force, and such transitions were attributed to the stacking characteristics of poly(dA). Using single-molecule manipulation studies, we found that poly(dA) has two stretching pathways at high forces. The previously observed pathway has a free energy that is less than what is expected of single-stranded DNA with a random sequence, indicating the existence of a novel conformation of poly(dA) at large extensions. We also observed stepwise transitions between the two pathways by pulling the molecule with constant force, and found that the transitions are cooperative. These results suggest that the unique mechanical property of poly(dA) may play an important role in biological processes such as gene expression.

Action at hooked or twisted-hooked DNA juxtapositions rationalizes unlinking preference of type-2 topoisomerases

The mathematical basis of the hypothesis that type-2 topoisomerases recognize and act at specific DNA juxtapositions has been investigated by coarse-grained lattice polymer models, showing that selective segment passages at hooked juxtapositions can result in dramatic reductions in catenane and knot populations. The lattice modeling approach is here extended to account for the narrowing of variance of linking number (Lk) of DNA circles by type-2 topoisomerases. In general, the steady-state variance of Lk resulting from selective segment passages at a specific juxtaposition geometry j is inversely proportional to the average linking number, Lk(j), of circles with the given juxtaposition. Based on this formulation, we demonstrate that selective segment passages at hooked juxtapositions reduce the variance of Lk. The dependence of this effect on model DNA circle size is remarkably similar to that observed experimentally for type-2 topoisomerases, which appear to be less capable in narrowing Lk variance for small DNA circles than for larger DNA circles. This behavior is rationalized by a substantial cancellation of writhe in small circles with hook-like juxtapositions. During our simulations, we uncovered a twisted variation of the hooked juxtaposition that has an even more dramatic effect on Lk variance narrowing than the hooked juxtaposition. For an extended set of juxtapositions, we detected a significant correlation between the Lk narrowing potential and the logarithmic decatenating and unknotting potentials for a given juxtaposition, a trend reminiscent of scaling relations observed with experimental measurements on type-2 topoisomerases from a variety of organisms. The consistent agreement between theory and experiment argues for type-2 topoisomerase action at hooked or twisted-hooked DNA juxtapositions.