Computer simulation of DNA supercoiling

DNA supercoiling in bacteria: state of play and challenges from a viewpoint of physics based modeling

DNA supercoiling is central to many fundamental processes of living organisms. Its average level along the chromosome and over time reflects the dynamic equilibrium of opposite activities of topoisomerases, which are required to relax mechanical stresses that are inevitably produced during DNA replication and gene transcription. Supercoiling affects all scales of the spatio-temporal organization of bacterial DNA, from the base pair to the large scale chromosome conformation. Highlighted in vitro and in vivo in the 1960s and 1970s, respectively, the first physical models were proposed concomitantly in order to predict the deformation properties of the double helix. About fifteen years later, polymer physics models demonstrated on larger scales the plectonemic nature and the tree-like organization of supercoiled DNA. Since then, many works have tried to establish a better understanding of the multiple structuring and physiological properties of bacterial DNA in thermodynamic equilibrium and far from equilibrium. The purpose of this essay is to address upcoming challenges by thoroughly exploring the relevance, predictive capacity, and limitations of current physical models, with a specific focus on structural properties beyond the scale of the double helix. We discuss more particularly the problem of DNA conformations, the interplay between DNA supercoiling with gene transcription and DNA replication, its role on nucleoid formation and, finally, the problem of scaling up models. Our primary objective is to foster increased collaboration between physicists and biologists. To achieve this, we have reduced the respective jargon to a minimum and we provide some explanatory background material for the two communities.

DNA fluctuations reveal the size and dynamics of topological domains

DNA supercoiling is a key regulatory mechanism that orchestrates DNA readout, recombination, and genome maintenance. DNA-binding proteins often mediate these processes by bringing two distant DNA sites together, thereby inducing (transient) topological domains. In order to understand the dynamics and molecular architecture of protein-induced topological domains in DNA, quantitative and time-resolved approaches are required. Here, we present a methodology to determine the size and dynamics of topological domains in supercoiled DNA in real time and at the single-molecule level. Our approach is based on quantifying the extension fluctuations-in addition to the mean extension-of supercoiled DNA in magnetic tweezers (MT). Using a combination of high-speed MT experiments, Monte Carlo simulations, and analytical theory, we map out the dependence of DNA extension fluctuations as a function of supercoiling density and external force. We find that in the plectonemic regime, the extension variance increases linearly with increasing supercoiling density and show how this enables us to determine the formation and size of topological domains. In addition, we demonstrate how the transient (partial) dissociation of DNA-bridging proteins results in the dynamic sampling of different topological states, which allows us to deduce the torsional stiffness of the plectonemic state and the kinetics of protein-plectoneme interactions. We expect our results to further the understanding and optimization of magnetic tweezer measurements and to enable quantification of the dynamics and reaction pathways of DNA processing enzymes in the context of physiologically relevant forces and supercoiling densities.

Equilibrium fluctuations of DNA plectonemes

Plectonemes are intertwined helically looped domains which form when a DNA molecule is supercoiled, i.e., over- or underwound. They are ubiquitous in cellular DNA, and their physical properties have attracted significant interest both from the experimental side and from the modeling side. In this paper, we investigate fluctuations of the end-point distance z of supercoiled linear DNA molecules subject to external stretching forces. Our analysis is based on a two-phase model, which describes the supercoiled DNA as composed of a stretched phase and a plectonemic phase. A variety of mechanisms are found to contribute to extension fluctuations, characterized by the variance 〈Δz^{2}〉. We find the dominant contribution to 〈Δz^{2}〉 to originate from phase-exchange fluctuations, the transient shrinking and expansion of plectonemes, which is accompanied by an exchange of molecular length between the two phases. We perform Monte Carlo simulations of the twistable wormlike chain and analyze the fluctuation of various quantities, the results of which are found to agree with the two-phase model predictions. Furthermore, we show that the extension and its variance at high forces are very well captured by the two-phase model, provided that one goes beyond quadratic approximations.

Two-phase dynamics of DNA supercoiling based on DNA polymer physics

DNA supercoils are generated in genome regulation processes such as transcription and replication and provide mechanical feedback to such processes. Under tension, a DNA supercoil can present a coexistence state of plectonemic and stretched phases. Experiments have revealed the dynamic behaviors of plectonemes, e.g., diffusion, nucleation, and hopping. To represent these dynamics with conformational changes, we demonstrated first the fast dynamics on the DNA to reach torque equilibrium within the plectonemic and stretched phases, and then identified the two-phase boundaries as collective slow variables to describe the essential dynamics. According to the timescale separation demonstrated here, we developed a two-phase model on the dynamics of DNA supercoiling, which can capture physiologically relevant events across timescales of several orders of magnitudes. In this model, we systematically characterized the slow dynamics between the two phases and compared the numerical results with those from the DNA polymer physics-based worm-like chain model. The supercoiling dynamics, including the nucleation, diffusion, and hopping of plectonemes, have been well represented and reproduced, using the two-phase dynamic model, at trivial computational costs. Our current developments, therefore, can be implemented to explore multiscale physical mechanisms of the DNA supercoiling-dependent physiological processes.

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.

Influence of Nucleoid-Associated Proteins on DNA Supercoiling

DNA supercoiling, where the DNA strand forms a writhe to relieve torsional stress, plays a vital role in packaging the genetic material in cells. Experiment, simulation, and theory have all demonstrated how supercoiling emerges due to the over- or underwinding of the DNA strand. Nucleoid-associated proteins (NAPs) help structure DNA in prokaryotes, yet the role that they play in the supercoiling process has not been as thoroughly investigated. We develop a coarse-grained simulation to model DNA supercoiling in the presence of proteins, providing a rigorous physical understanding of how NAPs affect supercoiling behavior. Specifically, we demonstrate how the force and torque necessary to form supercoils are affected by the presence of NAPs. NAPs that bend DNA stabilize the supercoil, thus shifting the transition between extended and supercoiled DNAs. We develop a theory to explain how NAP binding affects DNA supercoiling. This provides insight into how NAPs modulate DNA compaction via a combination of supercoiling and local protein-dependent deformations.

Geometrical structures for radiation biology research as implemented in the TOPAS-nBio toolkit

Computational simulations, such as Monte Carlo track structure simulations, offer a powerful tool for quantitatively investigating radiation interactions within cells. The modelling of the spatial distribution of energy deposition events as well as diffusion of chemical free radical species, within realistic biological geometries, can help provide a comprehensive understanding of the effects of radiation on cells. Track structure simulations, however, generally require advanced computing skills to implement. The TOPAS-nBio toolkit, an extension to TOPAS (TOol for PArticle Simulation), aims to provide users with a comprehensive framework for radiobiology simulations, without the need for advanced computing skills. This includes providing users with an extensive library of advanced, realistic, biological geometries ranging from the micrometer scale (e.g. cells and organelles) down to the nanometer scale (e.g. DNA molecules and proteins). Here we present the geometries available in TOPAS-nBio.

Statistical mechanical model for a closed loop plectoneme with weak helix specific forces

We develop a statistical mechanical framework, based on a variational approximation, to describe closed loop plectonemes. This framework incorporates weak helix structure dependent forces into the determination of the free energy and average structure of a plectoneme. Notably, due to their chiral nature, helix structure dependent forces break the symmetry between left and right handed supercoiling. The theoretical approach, presented here, also provides a systematic way of enforcing the topological constraint of closed loop supercoiling in the variational approximation. At large plectoneme lengths, by considering correlation functions in an expansion in terms of the spatial mean twist density about its thermally averaged value, it can be argued that topological constraint may be approximated by replacing twist and writhe by their thermal averages. A Lagrange multiplier, containing the sum of average twist and writhe, can be added to the free energy to conveniently inforce this result. The average writhe can be calculated through the thermal average of the Gauss' integral in the variational approximation. Furthermore, this approach allows for a possible way to calculate finite size corrections due to the topological constraint. Using interaction energy terms from the mean-field Kornyshev-Leikin theory, for parameter values that correspond to weak helix dependent forces, we calculate the free energy, fluctuation magnitudes and mean geometric parameters for the plectoneme. We see a slight asymmetry, where interestingly, left handed supercoils have a looser structure than right handed ones, although with a lower free energy, unlike what the previous ground state calculations would suggest.

Effect of sequence-dependent rigidity on plectoneme localization in dsDNA

We use Monte-Carlo simulations to study the effect of variable rigidity on plectoneme formation and localization in supercoiled double-stranded DNA. We show that the presence of soft sequences increases the number of plectoneme branches and that the edges of the branches tend to be localized at these sequences. We propose an experimental approach to test our results in vitro, and discuss the possible role played by plectoneme localization in the search process of transcription factors for their targets (promoter regions) on the bacterial genome.

Thermodynamics of long supercoiled molecules: insights from highly efficient Monte Carlo simulations

Supercoiled DNA polymer models for which the torsional energy depends on the total twist of molecules (Tw) are a priori well suited for thermodynamic analysis of long molecules. So far, nevertheless, the exact determination of Tw in these models has been based on a computation of the writhe of the molecules (Wr) by exploiting the conservation of the linking number, Lk=Tw+Wr, which reflects topological constraints coming from the helical nature of DNA. Because Wr is equal to the number of times the main axis of a DNA molecule winds around itself, current Monte Carlo algorithms have a quadratic time complexity, O(L(2)), with respect to the contour length (L) of the molecules. Here, we present an efficient method to compute Tw exactly, leading in principle to algorithms with a linear complexity, which in practice is O(L(1.2)). Specifically, we use a discrete wormlike chain that includes the explicit double-helix structure of DNA and where the linking number is conserved by continuously preventing the generation of twist between any two consecutive cylinders of the discretized chain. As an application, we show that long (up to 21 kbp) linear molecules stretched by mechanical forces akin to magnetic tweezers contain, in the buckling regime, multiple and branched plectonemes that often coexist with curls and helices, and whose length and number are in good agreement with experiments. By attaching the ends of the molecules to a reservoir of twists with which these can exchange helix turns, we also show how to compute the torques in these models. As an example, we report values that are in good agreement with experiments and that concern the longest molecules that have been studied so far (16 kbp).

Plectoneme tip bubbles: coupled denaturation and writhing in supercoiled DNA

We predict a novel conformational regime for DNA, where denaturation bubbles form at the tips of plectonemes, and study its properties using coarse-grained simulations. For negative supercoiling, this regime lies between bubble-dominated and plectoneme-dominated phases, and explains the broad transition between the two observed in experiment. Tip bubbles cause localisation of plectonemes within thermodynamically weaker AT-rich sequences, and can greatly suppress plectoneme diffusion by a pinning mechanism. They occur for supercoiling densities and forces that are typically encountered for DNA in vivo, and may be exploited for biological control of genomic processes.

Efficient deformation algorithm for plasmid DNA simulations

BACKGROUND

Plasmid DNA molecules are closed circular molecules that are widely used in life sciences, particularly in gene therapy research. Monte Carlo methods have been used for several years to simulate the conformational behavior of DNA molecules. In each iteration these simulation methods randomly generate a new trial conformation, which is either accepted or rejected according to a criterion based on energy calculations and stochastic rules. These simulation trials are generated using a method based on crankshaft motion that, apart from some slight improvements, has remained the same for many years.

RESULTS

In this paper, we present a new algorithm for the deformation of plasmid DNA molecules for Monte Carlo simulations. The move underlying our algorithm preserves the size and connectivity of straight-line segments of the plasmid DNA skeleton. We also present the results of three experiments comparing our deformation move with the standard and biased crankshaft moves in terms of acceptance ratio of the trials, energy and temperature evolution, and average displacement of the molecule. Our algorithm can also be used as a generic geometric algorithm for the deformation of regular polygons or polylines that preserves the connections and lengths of their segments.

CONCLUSION

Compared with both crankshaft moves, our move generates simulation trials with higher acceptance ratios and smoother deformations, making it suitable for real-time visualization of plasmid DNA coiling. For that purpose, we have adopted a DNA assembly algorithm that uses nucleotides as building blocks.

Lac repressor mediated DNA looping: Monte Carlo simulation of constrained DNA molecules complemented with current experimental results

Tethered particle motion (TPM) experiments can be used to detect time-resolved loop formation in a single DNA molecule by measuring changes in the length of a DNA tether. Interpretation of such experiments is greatly aided by computer simulations of DNA looping which allow one to analyze the structure of the looped DNA and estimate DNA-protein binding constants specific for the loop formation process. We here present a new Monte Carlo scheme for accurate simulation of DNA configurations subject to geometric constraints and apply this method to Lac repressor mediated DNA looping, comparing the simulation results with new experimental data obtained by the TPM technique. Our simulations, taking into account the details of attachment of DNA ends and fluctuations of the looped subsegment of the DNA, reveal the origin of the double-peaked distribution of RMS values observed by TPM experiments by showing that the average RMS value for anti-parallel loop types is smaller than that of parallel loop types. The simulations also reveal that the looping probabilities for the anti-parallel loop types are significantly higher than those of the parallel loop types, even for loops of length 600 and 900 base pairs, and that the correct proportion between the heights of the peaks in the distribution can only be attained when loops with flexible Lac repressor conformation are taken into account. Comparison of the in silico and in vitro results yields estimates for the dissociation constants characterizing the binding affinity between O1 and Oid DNA operators and the dimeric arms of the Lac repressor.

Determining DNA supercoiling enthalpy by isothermal titration calorimetry

DNA supercoiling plays a critical role in certain essential DNA transactions, such as DNA replication, recombination, and transcription. For this reason, exploring energetics of DNA supercoiling is fundamentally important for understanding its biological functions. In this paper, using a unique property of DNA intercalators, such as ethidium bromide and daunorubicin, which bind to supercoiled, nicked, and relaxed DNA templates with different DNA-binding enthalpies, we determined DNA supercoiling enthalpy of plasmid pXXZ6, a 4.5 kb plasmid to be about 11.5 kcal/mol per linking number change. This determination allowed us to partition the DNA supercoiling free energy into enthalpic and entropic contributions where the unfavorable DNA supercoiling free energy exclusively originated from the large positive supercoiling enthalpy and was compensated by a large, favorable entropy term (TΔS).

Chiral electrostatics breaks the mirror symmetry of DNA supercoiling

DNA supercoiling plays a fundamental role in regulating cellular activity and in the packaging of genetic material. In this communication, we analyse the effect of attractive chiral forces on the conformation of a closed circular DNA molecule, arising due to the helical patterns of charges on the DNA. We propose a model for closed loop DNA which uses the results of the recent theory of electrostatic interactions of a braid of two free-ended DNA molecules. Our model reproduces the known features of DNA supercoiling in an environment of low ionic strength. In high salt conditions, and in the presence of counterions that have high affinity to the DNA grooves, helix-specific forces significantly affect the conformation of the molecule by favouring a state characterized by a central left-handed braided section where there is close contact between distant portions of the loop. In such an environment we predict a previously unexplored possibility that nicked or topologically relaxed DNA molecules adopt a writhed state. This prediction suggests an alternative explanation for experiments in which it was assumed that the most stable topoisomer is always an open circle. Our results also give the first plausible explanation for the occurrence of tightly interwound molecules observed in cryo-electron microscopy and atomic force microscopy in a high ionic strength environment. We suggest several new experiments to test the predictions of this theory.

Optimization of Monte Carlo trial moves for protein simulations

Closed rigid-body rotations of residue segments under bond-angle restraints are simple and effective Monte Carlo moves for searching the conformational space of proteins. The efficiency of these moves is examined here as a function of the number of moving residues and the magnitude of their displacement. It is found that the efficiency of folding and equilibrium simulations can be significantly improved by tailoring the distribution of the number of moving residues to the simulation temperature. In general, simulations exploring compact conformations are more efficient when the average number of moving residues is smaller. It is also demonstrated that the moves do not require additional restrictions on the magnitude of the rotation displacements and perform much better than other rotation moves that do not restrict the bond angles a priori. As an example, these results are applied to the replica exchange method. By assigning distributions that generate a smaller number of moving residues to lower temperature replicas, the simulation times are decreased as long as the higher temperature replicas are effective.

Energetics at the DNA supercoiling transition

Twisting a DNA molecule held under constant tension is accompanied by a transition from a linear to a plectonemic DNA configuration, in which part of the applied twist is absorbed in a superhelical structure. Recent experiments revealed the occurrence of an abrupt extension change at the onset of this transition. To elucidate its origin we study this abrupt DNA shortening using magnetic tweezers. We find that it strongly depends on the length of the DNA molecule and the ionic strength of the solution. This behavior can be well understood in the framework of a model in which the energy per writhe for the initial plectonemic loop is larger than for subsequent turns of the superhelix. By quantitative data analysis, relevant plectoneme energies and other parameters were extracted, providing good agreement with a simple theory. As a direct confirmation of the initial-loop model, we find that for a kinked DNA molecule the abrupt extension change occurs at significantly lower twist than the subsequent superhelix formation. This should allow pinning of the plectoneme position within supercoiled DNA if a kinked substrate is used, and enable the detection of enzymes and proteins which, themselves, bend or kink DNA.

Gel mobilities of linking-number topoisomers and their dependence on DNA helical repeat and elasticity

Agarose-gel electrophoresis has been used for more than thirty years to characterize the linking-number (Lk) distribution of closed-circular DNA molecules. Although the physical basis of this technique remains poorly understood, the gel-electrophoretic behavior of covalently closed DNAs has been used to determine the local unwinding of DNA by proteins and small-molecule ligands, characterize supercoiling-dependent conformational transitions in duplex DNA, and to measure helical-repeat changes due to shifts in temperature and ionic strength. Those results have been analyzed by assuming that the absolute mobility of a particular topoisomer is mainly a function of the integral number of superhelical turns, and thus a slowly varying function of plasmid molecular weight. In examining the mobilities of Lk topoisomers for a series of plasmids that differ incrementally in size over more than one helical turn, we found that the size-dependent agarose-gel mobility of individual topoisomers with identical values of Lk (but different values of the excess linking number, DeltaLk) vary dramatically over a duplex turn. Our results suggest that a simple semi-empirical relationship holds between the electrophoretic mobility of linking-number topoisomers and their average writhe in solution.

Coupling of twist and writhe in short DNA loops

While bending and twist can be treated as independent degrees of freedom for linear DNA molecules, the loop closure constraint introduces a coupling between these variables in circular DNA. We performed Monte Carlo simulations of wormlike rods with both bending and twist rigidity in order to study the coupling between the writhe and twist distributions for various DNA lengths. We find that for sufficiently short DNA, the writhe distribution differs from that of a model with bending energy only. We show that the factorization approximation introduced by previous researchers coincides, within numerical accuracy, with our simulation results, and conclude that the closure constraint is fully accounted for by the White-Fuller relation. Experimental tests of our results for short DNA plasmids are proposed.

Elasticity and electrostatics of plectonemic DNA

We present a self-contained theory for the mechanical response of DNA in single molecule experiments. Our model is based on a one-dimensional continuum description of the DNA molecule and accounts both for its elasticity and for DNA-DNA electrostatic interactions. We consider the classical loading geometry used in experiments where one end of the molecule is attached to a substrate and the other one is pulled by a tensile force and twisted by a given number of turns. We focus on configurations relevant to the limit of a large number of turns, which are made up of two phases, one with linear DNA and the other one with superhelical DNA. The model takes into account thermal fluctuations in the linear phase and electrostatic interactions in the superhelical phase. The values of the torsional stress, of the supercoiling radius and angle, and key features of the experimental extension-rotation curves, namely the slope of the linear region and thermal buckling threshold, are predicted. They are found in good agreement with experimental data.

CRANKITE: A fast polypeptide backbone conformation sampler

BACKGROUND

CRANKITE is a suite of programs for simulating backbone conformations of polypeptides and proteins. The core of the suite is an efficient Metropolis Monte Carlo sampler of backbone conformations in continuous three-dimensional space in atomic details.

METHODS

In contrast to other programs relying on local Metropolis moves in the space of dihedral angles, our sampler utilizes local crankshaft rotations of rigid peptide bonds in Cartesian space.

RESULTS

The sampler allows fast simulation and analysis of secondary structure formation and conformational changes for proteins of average length.

Shapes of semiflexible polymer rings

The shape of semiflexible polymer rings is studied over their whole range of flexibility. Investigating the joint distribution of asphericity and the nature of asphericity as well as their respective averages, we find two distinct shape regimes depending on the flexibility of the polymer. For a small perimeter to persistence length the fluctuating rings exhibit only planar, elliptical configurations. At higher flexibilities three-dimensional, crumpled structures arise. Analytic calculations confirm the qualitative behavior of the averaged shape parameters and the elliptical shape in the stiff regime.

Fluctuating semiflexible polymer ribbon constrained to a ring

Twist stiffness and an asymmetric bending stiffness of a polymer or a polymer bundle is captured by the elastic ribbon model. We investigate the effects a ring geometry induces to a thermally fluctuating ribbon, finding bend-bend coupling in addition to twist-bend coupling. Furthermore, due to the geometric constraint the polymer's effective bending stiffness increases. A new parameter for experimental investigations of polymer bundles is proposed: the mean square diameter of a ribbonlike ring, which is determined analytically in the semiflexible limit. Monte Carlo simulations are performed which affirm the model's prediction up to high flexibility.

Torque and dynamics of linking number relaxation in stretched supercoiled DNA

In micromechanical studies of DNA, plectonemically supercoiled domains are often used as sources of constant torque. These torques are not easily measured and are instead usually estimated. Here, coexisting extended and supercoiled DNA domains are analyzed, and closed-form expressions for the dependence of extension and torque on force and linking number are presented. When there are coexisting domains of plectonemic and extended DNA, the torque depends only on force, with no dependence on linking number. However, torque depends on force in a manner more complex than a simple power law, involving the free energy of the extended and plectonemic DNA. A simple strategy is described for measurement of the free energies of both extended and plectonemic DNA without reference to specific microscopic polymer models. Applications of the theory to analysis of relaxation of supercoiling by enzymes which permit friction-controlled rotational relaxation of linking number is also presented. Such enzymes must display a breaking of symmetry between relaxations driven by equal magnitude but opposite direction torques.

Exhaustive Metropolis Monte Carlo sampling and analysis of polyalanine conformations adopted under the influence of hydrogen bonds

We propose a novel Metropolis Monte Carlo procedure for protein modeling and analyze the influence of hydrogen bonding on the distribution of polyalanine conformations. We use an atomistic model of the polyalanine chain with rigid and planar polypeptide bonds, and elastic alpha carbon valence geometry. We adopt a simplified energy function in which only hard-sphere repulsion and hydrogen bonding interactions between the atoms are considered. Our Metropolis Monte Carlo procedure utilizes local crankshaft moves and is combined with parallel tempering to exhaustively sample the conformations of 16-mer polyalanine. We confirm that Flory's isolated-pair hypothesis (the steric independence between the dihedral angles of individual amino acids) does not hold true in long polypeptide chains. In addition to 3(10)- and alpha-helices, we identify a kink stabilized by 2 hydrogen bonds with a shared acceptor as a common structural motif. Varying the strength of hydrogen bonds, we induce the helix-coil transition in the model polypeptide chain. We compare the propensities for various hydrogen bonding patterns and determine the degree of cooperativity of hydrogen bond formation in terms of the Hill coefficient. The observed helix-coil transition is also quantified according to Zimm-Bragg theory.

DNA Mechanics

We review the history of DNA mechanics and its analysis. We evaluate several methods to analyze the structures of superhelical DNA molecules, each predicated on the assumption that DNA can be modeled with reasonable accuracy as an extended, linearly elastic polymer. Three main approaches are considered: mechanical equilibrium methods, which seek to compute minimum energy conformations of topologically constrained molecules; statistical mechanical methods, which seek to compute the Boltzmann distribution of equilibrium conformations that arise in a finite temperature environment; and dynamic methods, which seek to compute deterministic trajectories of the helix axis by solving equations of motion. When these methods include forces of self-contact, which prevent strand passage and preserve the topological constraint, each predicts plectonemically interwound structures. On the other hand, the extent to which these mechanical methods reliably predict energetic and thermodynamic properties of superhelical molecules is limited, in part because of their inability to account explicitly for interactions involving solvent. Monte Carlo methods predict the entropy associated with supercoiling to be negative, in conflict with a body of experimental evidence that finds it is large and positive, as would be the case if superhelical deformations significantly disrupt the ordering of ambient solvent molecules. This suggests that the large-scale conformational properties predicted by elastomechanical models are not the only ones determining the energetics and thermodynamics of supercoiling. Moreover, because all such models that preserve the topological constraint correctly predict plectonemic interwinding, despite these and other limitations, this constraint evidently dominates energetic and thermodynamic factors in determining supercoil geometry. Therefore, agreement between predicted structures and structures obtained experimentally, for example, by electron microscopy, does not in itself provide evidence for the correctness or completeness of any given model of DNA mechanics.

Braiding DNA: experiments, simulations, and models

DNA encounters topological problems in vivo because of its extended double-helical structure. As a consequence, the semiconservative mechanism of DNA replication leads to the formation of DNA braids or catenanes, which have to be removed for the completion of cell division. To get a better understanding of these structures, we have studied the elastic behavior of two braided nicked DNA molecules using a magnetic trap apparatus. The experimental data let us identify and characterize three regimes of braiding: a slightly twisted regime before the formation of the first crossing, followed by genuine braids which, at large braiding number, buckle to form plectonemes. Two different approaches support and quantify this characterization of the data. First, Monte Carlo (MC) simulations of braided DNAs yield a full description of the molecules' behavior and their buckling transition. Second, modeling the braids as a twisted swing provides a good approximation of the elastic response of the molecules as they are intertwined. Comparisons of the experiments and the MC simulations with this analytical model allow for a measurement of the diameter of the braids and its dependence upon entropic and electrostatic repulsive interactions. The MC simulations allow for an estimate of the effective torsional constant of the braids (at a stretching force F = 2 pN): C(b) approximately 48 nm (as compared with C approximately 100 nm for a single unnicked DNA). Finally, at low salt concentrations and for sufficiently large number of braids, the diameter of the braided molecules is observed to collapse to that of double-stranded DNA. We suggest that this collapse is due to the partial melting and fraying of the two nicked molecules and the subsequent right- or left-handed intertwining of the stretched single strands.

Computational analysis of DNA gyrase action

DNA gyrase introduces negative supercoiling into circular DNA by catalyzing the passage of one DNA segment through another. The efficiency of the reaction is many times higher than that of other topological transformations. We analyze, by a computer simulation, the reaction selectivity for a model of DNA gyrase action that assumes existence of a free loop between the G- and T- DNA segments participating in the reaction. A popular model of this type assumed that the selectivity can be provided by the conformation of the DNA segment wrapped around the enzyme into the right-handed helix turn (G-segment). We simulated the distribution of the reaction products for this model. Equilibrium sets of DNA conformations with one segment of the double helix wrapped around the enzyme were constructed. From these sets we selected conformations that had a second segment properly juxtaposed with the first one. Assuming that the juxtapositions result in the strand-passing reaction, we calculated the reaction products for all these conformations. The results show that different products have to be formed if the enzyme acts according to the model. This conclusion can be extended for any model with a free loop between the G- and T-segments. An alternative model that is consistent with the major experimental observations and the computational analysis, is suggested.

Effects of small neutral osmolytes on the supercoiling free energy and intrinsic twist of p30? DNA

Both theory and experiments are employed to investigate the effects of small neutral osmolytes on the average intrinsic twist (l0), the torsion and bending elastic constants, and the twist energy parameter (ET) that governs the supercoiling free energy. The experimental data for ethylene glycol and acetamide at 37 degrees C suggest, and are interpreted in terms of, a model wherein the DNA exhibits an equilibrium between two distinct conformational states that possess different numbers of bound water molecules and exhibit different intrinsic twists and torsion and bending elastic constants. Expressions are derived to relate the effective ET and l0 to the equilibrium constant, water activity (aw), and number (n) of bound water molecules released per cooperative domain undergoing the two-state transition. The variations of l0 and ET with -ln(aw) are similar for acetamide and ethylene glycol at 37 degrees C. Fitting the theory to those data yields the range n = 103-125 for ethylene glycol and n = 71-113 for acetamide, depending on the assumed value of ET for the dehydrated state. The cooperative domain size of the two-state transition is estimated to exceed about 25-30 base pairs (bp). Between 0 and 19.4 w/v % ethylene glycol, the torsion elastic constant, measured by time-resolved fluorescence polarization anisotropy (FPA), increases by 1.37-fold, whereas the measured ET decreases by 1.15-fold over that same range. The implied decrease in bending rigidity over that range is by a factor of about 0.7. The variations of l0 and ET with increasing -ln(aw) due to added ethylene glycol at 37 degrees C are far smaller than the corresponding variations observed previously at 14 and 15 degrees C. However, at 21 degrees C, upon adding either ethylene glycol or acetamide, l0 and ET initially decline steeply with increasing -ln(aw), with slopes possibly comparable to those seen at 14 and 15 degrees C, but then flatten out and follow curves similar to those at 37 degrees C. Possible origins of such mixed behavior are discussed. The effects of betaine at both 37 and 21 degrees C differ qualitatively and quantitatively in various respects from those of ethylene glycol and acetamide. Upon adding sucrose, l0 initially jumps to higher plateaus at both 37 and 21 degrees C, but its effects on ET cannot be reliably assessed, due to the limited range of -ln(aw).

Investigation of viral DNA packaging using molecular mechanics models

A simple molecular mechanics model has been used to investigate optimal spool-like packing conformations of double-stranded DNA molecules in viral capsids with icosahedral symmetry. The model represents an elastic segmented chain by using one pseudoatom for each ten basepairs (roughly one turn of the DNA double helix). Force constants for the various terms in the energy function were chosen to approximate known physical properties, and a radial restraint was used to confine the DNA into a sphere with a volume corresponding to that of a typical bacteriophage capsid. When the DNA fills 90% of the spherical volume, optimal packaging is obtained for coaxially spooled models, but this result does not hold when the void volume is larger. When only 60% of the spherical volume is filled with DNA, the lowest energy structure has two layers, with a coiled core packed at an angle to an outer coaxially spooled shell. This relieves bending strain associated with tight curvature near the poles in a model with 100% coaxial spooling. Interestingly, the supercoiling density of these models is very similar to typical values observed in plasmids in bacterial cells. Potential applications of the methodology are also discussed.

Computational analysis of the chiral action of type II DNA topoisomerases

It was found recently that bacterial type II DNA topoisomerase, topo IV, is much more efficient in relaxing (+) DNA supercoiling than (-) supercoiling. This means that the DNA-enzyme complex is chiral. This chirality can appear upon binding the first segment that participates in the strand passing reaction (G segment) or only after the second segment (T segment) joins the complex. The former possibility is analyzed here. We assume that upon binding the enzyme, the G segment forms a part of left-handed helical turn. This model is an extension of the hairpin model introduced earlier to explain simplification of DNA topology by these enzymes. Using statistical-mechanical simulation of DNA properties, we estimated different consequences of the model: (1) relative rates of relaxation of (+) and (-) supercoiling by the enzyme; (2) the distribution of positions of the G segment in supercoiled molecules; (3) steady-state distribution of knots in circular molecules created by the topoisomerase; (4) the variance of topoisomer distribution created by the enzyme; (5) the effect of (+) and (-) supercoiling on the binding topo II with G segment. The simulation results are capable of explaining nearly all available experimental data, at least semiquantitatively. A few predictions obtained in the model analysis can be tested experimentally.

Surface-induced DNA superhelicity

Certain biopolymers, such as DNA, have a double-stranded twisted structure and frequently exhibit a supercoiled conformation. Over the past decade, extensive conformational analyses of different biopolymers have been performed using atomic force microscopy. In this technique, a necessary step in sample preparation is the adsorption of molecules on the surface, which could affect the chain conformation. Using a Monte Carlo simulation, we studied the adsorption process of circular semiflexible twisted double-stranded polymer chains on a solid surface with an emphasis on the conformational properties. We found that the conformation of an adsorbed chain strongly depends on the number of double-helical turns of the chain. Chains with weak twisting adsorbed on the surface assumes higher superhelical twisting, whereas the coiled state is seen in the bulk solution. After the adsorption, double-helical turns are accumulated in the adsorbed ("train") sections giving more conformational freedom to the nonadsorbed ("loop") sections. Chains with strong twisting show small conformational changes with adsorption. In both cases, superhelicity shows the opposite sign of writhing.

Time-lapse imaging of conformational changes in supercoiled DNA by scanning force microscopy

Most of the scanning force microscopy (SFM) images of supercoiled DNA on untreated mica thus far reported have not shown tight plectonemic structure seen by electron microscopy, but instead less coiled molecules and sometimes a partly "condensed" state with intimate chain-chain interactions. By observing time-lapse images of conformational changes of DNA induced by decreasing ionic strength of imaging buffer in solution SFM, we could show that the process of water rinsing, an indispensable step for preparation of dried samples, may be responsible for some of the conformational anomalies in the images previously reported. We have studied several protocols to observe supercoiled DNA molecules by SFM and discuss the merits and the demerits. Images obtained following uranyl acetate treatment may be ideal for the detection of DNA damage, as the supercoiled and nicked forms are easily distinguishable.

Kinetic proofreading can explain the supression of supercoiling of circular DNA molecules by type-II topoisomerases

The enzymes that pass DNA through DNA so as to remove entanglements, adenosine-triphosphate-hydrolyzing type-II topoisomerases, are able to suppress the probability of self-entanglements (knots) and mutual entanglements (links) between approximately 10 kb plasmids, well below the levels expected, given the assumption that the topoisomerases pass DNA segments at random by thermal motion. This implies that a 10-nm type-II topoisomerase can somehow sense the topology of a large DNA. We previously introduced a "kinetic proofreading" model which supposes the enzyme to require two successive collisions in order to allow exchange of DNA segments, and we showed how it could quantitatively explain the reduction in knotting and linking complexity. Here we show how the same model quantitatively explains the reduced variance of the double-helix linking number (supercoiling) distribution observed experimentally.

Multimerization-cyclization of DNA fragments as a method of conformational analysis

Ligation of short DNA fragments results in the formation of linear and circular multimers of various lengths. The distribution of products in such a reaction is often used to evaluate fragment bending caused by specific chemical modification, by bound ligands or by the presence of irregular structural elements. We have developed a more rigorous quantitative approach to the analysis of such experimental data based on determination of j-factors for different multimers from the distribution of the reaction products. j-Factors define the effective concentration of one end of a linear chain in the vicinity of the other end. To extract j-factors we assumed that kinetics of the reaction is described by a system of differential equations where j-factors appear as coefficients. The assumption was confirmed by comparison with experimental data obtained here for DNA fragments containing A-tracts. At the second step of the analysis j-factors are used to determine conformational parameters of DNA fragments: the equilibrium bend angle, the bending rigidity of the fragment axis, and the total twist of the fragments. This procedure is based on empirical equations that connect the conformational parameters with the set of j-factors. To obtain the equations, we computed j-factors for a large array of conformational parameters that describe model fragments. The approach was tested on both simulated and actual experimental data for DNA fragments containing A-tracts. A-tract DNA bend angle determined here is in good agreement with previously published data. We have established a set of experimental conditions necessary for the data analysis to be successful.

DNA rings with multiple energy minima

Within the context of DNA rings, we analyze the relationship between intrinsic shape and the existence of multiple stable equilibria, either nicked or cyclized with the same link. A simple test, based on a perturbation expansion of symmetry breaking within a continuum elastic rod model, provides good predictions of the occurrence of such multiple equilibria. The reliability of these predictions is verified by direct computation of nicked and cyclized equilibria for several thousand DNA minicircles with lengths of 200 and 900 bp. Furthermore, our computations of equilibria for nicked rings predict properties of the equilibrium distribution of link, as calculated by much more computationally intensive Monte Carlo simulations.

Equilibrium distributions of topological states in circular DNA: interplay of supercoiling and knotting

Two variables define the topological state of closed double-stranded DNA: the knot type, K, and DeltaLk, the linking number difference from relaxed DNA. The equilibrium distribution of probabilities of these states, P(DeltaLk, K), is related to two conditional distributions: P(DeltaLk|K), the distribution of DeltaLk for a particular K, and P(K|DeltaLk) and also to two simple distributions: P(DeltaLk), the distribution of DeltaLk irrespective of K, and P(K). We explored the relationships between these distributions. P(DeltaLk, K), P(DeltaLk), and P(K|DeltaLk) were calculated from the simulated distributions of P(DeltaLk|K) and of P(K). The calculated distributions agreed with previous experimental and theoretical results and greatly advanced on them. Our major focus was on P(K|DeltaLk), the distribution of knot types for a particular value of DeltaLk, which had not been evaluated previously. We found that unknotted circular DNA is not the most probable state beyond small values of DeltaLk. Highly chiral knotted DNA has a lower free energy because it has less torsional deformation. Surprisingly, even at |DeltaLk| > 12, only one or two knot types dominate the P(K|DeltaLk) distribution despite the huge number of knots of comparable complexity. A large fraction of the knots found belong to the small family of torus knots. The relationship between supercoiling and knotting in vivo is discussed.

Effect of magnesium on cruciform extrusion in supercoiled DNA

Recently, it was reported that Mg2+greatly facilitates cruciform extrusion in the short palindromes of supercoiled DNA, thereby allowing the formation of cruciform structures in vivo. Because of the potential biological importance of this phenomenon, we undertook a broader study of the effect of Mg2+on a cruciform extrusion in supercoiled DNA. The method of two-dimensional gel electrophoresis was used to detect the cruciform extrusion both in the absence and in the presence of these ions. Our results show that Mg2+shifts the cruciform extrusion in the d(CCC(AT)16GGG) palindrome to a higher, rather than to a lower level of supercoiling. In order to study possible sequence-specific properties of the short palindromes for which the unusual cruciform extrusion in the presence Mg2+was reported, we constructed a plasmid with a longer palindromic region. This region bears the same sequences in the hairpin loops and four-arm junction as the short palindrome, except that the short stems of the hairpins are extended. The extension allowed us to overcome the limitation of our experimental approach which cannot be used for very short palindromes. Our results show that Mg2+also shifts the cruciform extrusion in this palindrome to a higher level of supercoiling. These data suggest that cruciform extrusion in the short palindromes at low supercoiling is highly improbable. We performed a thermodynamic analysis of the effect of Mg2+on cruciform extrusion. The treatment accounted for the effect of Mg2+on both free energy of supercoiling and the free energy of cruciform structure per se. Our analysis showed that although the level of supercoiling required for the cruciform extrusion is not reduced by Mg2+, the ions reduce the free energy of the cruciform structure.

Electrostatic-undulatory theory of plectonemically supercoiled DNA

We present an analytical calculation of the electrostatic interaction in a plectonemic supercoil within the Poisson-Boltzmann approximation. Undulations of the supercoil strands arising from thermal motion couple nonlinearly with the electrostatic interaction, giving rise to a strong enhancement of the bare interaction. In the limit of fairly tight winding, the free energy of a plectonemic supercoil may be split into an elastic contribution containing the bending and torsional energies and an electrostatic-undulatory free energy. The total free energy of the supercoil is minimized according to an iterative scheme, which utilizes the special symmetry inherent in the usual elastic free energy of the plectoneme. The superhelical radius, opening angle, and undulation amplitudes in the radius and pitch are obtained as a function of the specific linking difference and the concentration of monovalent salt. Our results compare favorably with the experimental values for these parameters of Boles et al. (1990. J. Mol. Biol. 213:931-951). In particular, we confirm the experimental observation that the writhe is a virtually constant fraction of the excess linking number over a wide range of superhelical densities. Another important prediction is the ionic strength dependence of the plectonemic parameters, which is in reasonable agreement with the results from computer simulations.