DNA-Topology Simplification by Topoisomerases.

The topological properties of DNA molecules, supercoiling, knotting, and catenation, are intimately connected with essential biological processes, such as gene expression, replication, recombination, and chromosome segregation. Non-trivial DNA topologies present challenges to the molecular machines that process and maintain genomic information, for example, by creating unwanted DNA entanglements. At the same time, topological distortion can facilitate DNA-sequence recognition through localized duplex unwinding and longer-range loop-mediated interactions between the DNA sequences. Topoisomerases are a special class of essential enzymes that homeostatically manage DNA topology through the passage of DNA strands. The activities of these enzymes are generally investigated using circular DNA as a model system, in which case it is possible to directly assay the formation and relaxation of DNA supercoils and the formation/resolution of knots and catenanes. Some topoisomerases use ATP as an energy cofactor, whereas others act in an ATP-independent manner. The free energy of ATP hydrolysis can be used to drive negative and positive supercoiling or to specifically relax DNA topologies to levels below those that are expected at thermodynamic equilibrium. The latter activity, which is known as topology simplification, is thus far exclusively associated with type-II topoisomerases and it can be understood through insight into the detailed non-equilibrium behavior of type-II enzymes. We use a non-equilibrium topological-network approach, which stands in contrast to the equilibrium models that are conventionally used in the DNA-topology field, to gain insights into the rates that govern individual transitions between topological states. We anticipate that our quantitative approach will stimulate experimental work and the theoretical/computational modeling of topoisomerases and similar enzyme systems.

Loop-closure kinetics reveal a stable, right-handed DNA intermediate in Cre recombination.

In Cre site-specific recombination, the synaptic intermediate is a recombinase homotetramer containing a pair of loxP DNA target sites. The enzyme system's strand-exchange mechanism proceeds via a Holliday-junction (HJ) intermediate; however, the geometry of DNA segments in the synapse has remained highly controversial. In particular, all crystallographic structures are consistent with an achiral, planar Holliday-junction (HJ) structure, whereas topological assays based on Cre-mediated knotting of plasmid DNAs are consistent with a right-handed chiral junction. We use the kinetics of loop closure involving closely spaced (131-151 bp) loxP sites to investigate the in-aqueo ensemble of conformations for the longest-lived looped DNA intermediate. Fitting the experimental site-spacing dependence of the loop-closure probability, J, to a statistical-mechanical theory of DNA looping provides evidence for substantial out-of-plane HJ distortion, which unequivocally stands in contrast to the square-planar intermediate geometry from Cre-loxP crystal structures and those of other int-superfamily recombinases. J measurements for an HJ-isomerization-deficient Cre mutant suggest that the apparent geometry of the wild-type complex is consistent with temporal averaging of right-handed and achiral structures. Our approach connects the static pictures provided by crystal structures and the natural dynamics of macromolecules in solution, thus advancing a more comprehensive dynamic analysis of large nucleoprotein structures and their mechanisms.

Kinetic pathways of topology simplification by Type-II topoisomerases in knotted supercoiled DNA.

The topological state of covalently closed, double-stranded DNA is defined by the knot type $K$ and the linking-number difference $\Delta Lk$ relative to unknotted relaxed DNA. DNA topoisomerases are essential enzymes that control the topology of DNA in all cells. In particular, type-II topoisomerases change both $K$ and $\Delta Lk$ by a duplex-strand-passage mechanism and have been shown to simplify the topology of DNA to levels below thermal equilibrium at the expense of ATP hydrolysis. It remains a key question how small enzymes are able to preferentially select strand passages that result in topology simplification in much larger DNA molecules. Using numerical simulations, we consider the non-equilibrium dynamics of transitions between topological states $(K,\Delta Lk)$ in DNA induced by type-II topoisomerases. For a biological process that delivers DNA molecules in a given topological state $(K,\Delta Lk)$ at a constant rate we fully characterize the pathways of topology simplification by type-II topoisomerases in terms of stationary probability distributions and probability currents on the network of topological states $(K,\Delta Lk)$. In particular, we observe that type-II topoisomerase activity is significantly enhanced in DNA molecules that maintain a supercoiled state with constant torsional tension. This is relevant for bacterial cells in which torsional tension is maintained by enzyme-dependent homeostatic mechanisms such as DNA-gyrase activity.

DNA cyclization and looping in the wormlike limit: Normal modes and the validity of the harmonic approximation.

For much of the last three decades, Monte Carlo-simulation methods have been the standard approach for accurately calculating the cyclization probability, J, or J factor, for DNA models having sequence-dependent bends or inhomogeneous bending flexibility. Within the last 10 years approaches based on harmonic analysis of semi-flexible polymer models have been introduced, which offer much greater computational efficiency than Monte Carlo techniques. These methods consider the ensemble of molecular conformations in terms of harmonic fluctuations about a well-defined elastic-energy minimum. However, the harmonic approximation is only applicable for small systems, because the accessible conformation space of larger systems is increasingly dominated by anharmonic contributions. In the case of computed values of the J factor, deviations of the harmonic approximation from the exact value of J as a function of DNA length have not been characterized. Using a recent, numerically exact method that accounts for both anharmonic and harmonic contributions to J for wormlike chains of arbitrary size, we report here the apparent error that results from neglecting anharmonic behavior. For wormlike chains having contour lengths less than four times the persistence length, the error in J arising from the harmonic approximation is generally small, amounting to free energies less than the thermal energy, kB T. For larger systems, however, the deviations between harmonic and exact J values increase approximately linearly with size.

Free-energy calculations for semi-flexible macromolecules: applications to DNA knotting and looping.

We present a method to obtain numerically accurate values of configurational free energies of semiflexible macromolecular systems, based on the technique of thermodynamic integration combined with normal-mode analysis of a reference system subject to harmonic constraints. Compared with previous free-energy calculations that depend on a reference state, our approach introduces two innovations, namely, the use of internal coordinates to constrain the reference states and the ability to freely select these reference states. As a consequence, it is possible to explore systems that undergo substantially larger fluctuations than those considered in previous calculations, including semiflexible biopolymers having arbitrary ratios of contour length L to persistence length P. To validate the method, high accuracy is demonstrated for free energies of prime DNA knots with L/P = 20 and L/P = 40, corresponding to DNA lengths of 3000 and 6000 base pairs, respectively. We then apply the method to study the free-energy landscape for a model of a synaptic nucleoprotein complex containing a pair of looped domains, revealing a bifurcation in the location of optimal synapse (crossover) sites. This transition is relevant to target-site selection by DNA-binding proteins that occupy multiple DNA sites separated by large linear distances along the genome, a problem that arises naturally in gene regulation, DNA recombination, and the action of type-II topoisomerases.

The thermodynamics of DNA loop formation, from J to Z.

The formation of DNA loops is a ubiquitous theme in biological processes, including DNA replication, recombination and repair, and gene regulation. These loops are mediated by proteins bound at specific sites along the contour of a single DNA molecule, in some cases many thousands of base pairs apart. Loop formation incurs a thermodynamic cost that is a sensitive function of the length of looped DNA as well as the geometry and elastic properties of the DNA-bound protein. The free energy of DNA looping is logarithmically related to a generalization of the Jacobson-Stockmayer factor for DNA cyclization, termed the J factor. In the present article, we review the thermodynamic origins of this quantity, discuss how it is measured experimentally and connect the macroscopic interpretation of the J factor with a statistical-mechanical description of DNA looping and cyclization.

Denaturation transition of stretched DNA.

In the last two decades, single-molecule force measurements using optical and magnetic tweezers and atomic force spectroscopy have dramatically expanded our knowledge of nucleic acids and proteins. These techniques characterize the force on a biomolecule required to produce a given molecular extension. When stretching long DNA molecules, the observed force-extension relationship exhibits a characteristic plateau at approximately 65 pN where the DNA may be extended to almost twice its B-DNA length with almost no increase in force. In the present review, I describe this transition in terms of the Poland-Scheraga model and summarize recent related studies.

Non-equilibrium Casimir force between vibrating plates.

We study the fluctuation-induced, time-dependent force between two plates confining a correlated fluid which is driven out of equilibrium mechanically by harmonic vibrations of one of the plates. For a purely relaxational dynamics of the fluid we calculate the fluctuation-induced force generated by the vibrating plate on the plate at rest. The time-dependence of this force is characterized by a positive lag time with respect to the driving. We obtain two distinctive contributions to the force, one generated by diffusion of stress in the fluid and another related to resonant dissipation in the cavity. The relation to the dynamic Casimir effect of the electromagnetic field and possible experiments to measure the time-dependent Casimir force are discussed.

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.

Single DNA denaturation and bubble dynamics.

While the Watson-Crick double-strand is the thermodynamically stable state of DNA in a wide range of temperature and salt conditions, even at physiological conditions local denaturation bubbles may open up spontaneously due to thermal activation. By raising the ambient temperature, titration, or by external forces in single molecule setups bubbles proliferate until full denaturation of the DNA occurs. Based on the Poland-Scheraga model we investigate both the equilibrium transition of DNA denaturation and the dynamics of the denaturation bubbles with respect to recent single DNA chain experiments for situations below, at, and above the denaturation transition. We also propose a new single molecule setup based on DNA constructs with two bubble zones to measure the bubble coalescence and extract the physical parameters relevant to DNA breathing. Finally we consider the interplay between denaturation bubbles and selectively single-stranded DNA binding proteins.

Kinetics and thermodynamics of salt-dependent T7 gene 2.5 protein binding to single- and double-stranded DNA.

Bacteriophage T7 gene 2.5 protein (gp2.5) is a single-stranded DNA (ssDNA)-binding protein that has essential roles in DNA replication, recombination and repair. However, it differs from other ssDNA-binding proteins by its weaker binding to ssDNA and lack of cooperative ssDNA binding. By studying the rate-dependent DNA melting force in the presence of gp2.5 and its deletion mutant lacking 26 C-terminal residues, we probe the kinetics and thermodynamics of gp2.5 binding to ssDNA and double-stranded DNA (dsDNA). These force measurements allow us to determine the binding rate of both proteins to ssDNA, as well as their equilibrium association constants to dsDNA. The salt dependence of dsDNA binding parallels that of ssDNA binding. We attribute the four orders of magnitude salt-independent differences between ssDNA and dsDNA binding to nonelectrostatic interactions involved only in ssDNA binding, in contrast to T4 gene 32 protein, which achieves preferential ssDNA binding primarily through cooperative interactions. The results support a model in which dimerization interactions must be broken for DNA binding, and gp2.5 monomers search dsDNA by 1D diffusion to bind ssDNA. We also quantitatively compare the salt-dependent ssDNA- and dsDNA-binding properties of the T4 and T7 ssDNA-binding proteins for the first time.

Denaturation transition of stretched DNA.

We generalize the Poland-Scheraga model to consider DNA denaturation in the presence of an external stretching force. We demonstrate the existence of a force-induced DNA denaturation transition and obtain the temperature-force phase diagram. The transition is determined by the loop exponent c, for which we find the new value c = 4 nu-1/2 such that the transition is second order with c = 1.85 < 2 in d = 3. We show that a finite stretching force F destabilizes DNA, corresponding to a lower melting temperature T(F), in agreement with single-molecule DNA stretching experiments.

Probing molecular free energy landscapes by periodic loading.

Single molecule pulling experiments provide information about interactions in biomolecules that cannot be obtained by any other method. However, the reconstruction of the molecule's free energy profile from the experimental data is still a challenge, in particular, for the unstable barrier regions. We propose a new method for obtaining the full profile by introducing a periodic ramp and using Jarzynski's relation for obtaining equilibrium quantities from nonequilibrium data. Our simulated experiments show that this method delivers significant more accurate data than previous methods, under the constraint of equal experimental effort.

Entropy loss in long-distance DNA looping.

The entropy loss due to the formation of one or multiple loops in circular and linear DNA chains is calculated from a scaling approach in the limit of long chain segments. The analytical results allow us to obtain a fast estimate for the entropy loss for a given configuration. Numerical values obtained for some examples suggest that the entropy loss encountered in loop closure in typical genetic switches may become a relevant factor in comparison to both k(B)T and typical bond energies in biopolymers, which has to be overcome by the released bond energy between the looping contact sites.

Tightness of slip-linked polymer chains.

We study the interplay between entropy and topological constraints for a polymer chain in which sliding rings (slip links) enforce pair contacts between monomers. These slip links divide a closed ring polymer into a number of subloops which can exchange length among each other. In the ideal chain limit, we find the joint probability density function for the sizes of segments within such a slip-linked polymer chain (paraknot). A particular segment is tight (small in size) or loose (of the order of the overall size of the paraknot) depending on both the number of slip links it incorporates and its competition with other segments. When self-avoiding interactions are included, scaling arguments can be used to predict the statistics of segment sizes for certain paraknot configurations.

Equilibrium shapes of flat knots.

We study the equilibrium shapes of prime and composite knots confined to two dimensions. Using scaling arguments we show that, due to self-avoiding effects, the topological details of prime knots are localized on a small portion of the larger ring polymer. Within this region, the original knot configuration can assume a hierarchy of contracted shapes, the dominating one given by just one small loop. This hierarchy is investigated in detail for the flat trefoil knot, and corroborated by Monte Carlo simulations.

Correlation functions near modulated and rough surfaces.

In a system with long-ranged correlations, the behavior of correlation functions is sensitive to the presence of a boundary. We show that surface deformations strongly modify this behavior as compared to a flat surface. The modified near surface correlations can be measured by scattering probes. To determine these correlations, we develop a perturbative calculation in the deformations in height from a flat surface. Detailed results are given for a regularly patterned surface, as well as for a self-affinely rough surface with roughness exponent zeta. By combining this perturbative calculation in height deformations with the field-theoretic renormalization-group approach, we also estimate the values of critical exponents governing the behavior of the decay of correlation functions near a self-affinely rough surface. We find that for the interacting theory, a large enough zeta can lead to a different surface critical behavior. We also provide scaling relations between roughness induced critical exponents for thermodynamic surface quantities.