Strong bending of the DNA double helix.

During the past decade, the issue of strong bending of the double helix has attracted a lot of attention. Here, we overview the major experimental and theoretical developments in the field sorting out reliably established facts from speculations and unsubstantiated claims. Theoretical analysis shows that sharp bends or kinks have to facilitate strong bending of the double helix. It remains to be determined what is the critical curvature of DNA that prompts the appearance of the kinks. Different experimental and computational approaches to the problem are analyzed. We conclude that there is no reliable evidence that any anomalous behavior of the double helix happens when DNA fragments in the range of 100 bp are circularized without torsional stress. The anomaly starts at the fragment length of about 70 bp when sharp bends or kinks emerge in essentially every molecule. Experimental data and theoretical analysis suggest that kinks may represent openings of isolated base pairs, which had been experimentally detected in linear DNA molecules. The calculation suggests that although the probability of these openings in unstressed DNA is close to 10(-5), it increases sharply in small DNA circles reaching 1 open bp per circle of 70 bp.

Temperature dependence of DNA persistence length.

We have determined the temperature dependence of DNA persistence length, a, using two different methods. The first approach was based on measuring the j-factors of short DNA fragments at various temperatures. Fitting the measured j-factors by the theoretical equation allowed us to obtain the values of a for temperatures between 5°C and 42°C. The second approach was based on measuring the equilibrium distribution of the linking number between the strands of circular DNA at different temperatures. The major contribution into the distribution variance comes from the fluctuations of DNA writhe in the nicked circular molecules which are specified by the value of a. The computation-based analysis of the measured variances was used to obtain the values of a for temperatures up to 60°C. We found a good agreement between the results obtained by these two methods. Our data show that DNA persistence length strongly depends on temperature and accounting for this dependence is important in quantitative comparison between experimental results obtained at different temperatures.

Sequence dependence of DNA bending rigidity.

For many aspects of DNA-protein interaction, it is vital to know how DNA bending rigidity (or persistence length, a) depends on its sequence. We addressed this problem using the method based on cyclization of short DNA fragments, which allows very accurate determination of a. Our approach was based on assigning specific values of a to each of 10 distinct dinucleotide steps. We prepared DNA fragments, each about 200 bp in length, with various quasi-periodic sequences, measured their cyclization efficiencies (j factors), and fitted the data by the theoretical equation to obtain the values of a for each fragment. From these data, we obtained a set of a for the dinucleotide steps. To test this set, we used it to design DNA sequences that should correspond to very low and very high values of a, prepared the corresponding fragments, and determined their values of a experimentally. The measured and calculated values of a were very close to one another, confirming that we have found the correct solution to this long-standing problem. The same experimental data also allowed us to determine the sequence dependence of DNA helical repeat.

DNA supercoiling helps to unlink sister duplexes after replication.

DNA supercoiling is one of the mechanisms that can help unlinking of newly replicated DNA molecules. Although DNA topoisomerases, which catalyze the strand passing of DNA segments through one another, make the unlinking problem solvable in principle, it remains difficult to complete the process that enables the separation of the sister duplexes. A few different mechanisms were developed by nature to solve the problem. Some of the mechanisms are very intuitive while the others, like topology simplification by type II DNA topoisomerases and DNA supercoiling, are not so evident. A computer simulation and analysis of linked sister plasmids formed in Escherichia coli cells with suppressed topoisomerase IV suggests an insight into the latter mechanism.

Unlinking of supercoiled DNA catenanes by type IIA topoisomerases.

It was found recently that DNA catenanes, formed during replication of circular plasmids, become positively (+) supercoiled, and the unlinking of such catenanes by type IIA topoisomerases proceeds much more efficiently than the unlinking of negatively (-) supercoiled catenanes. In an attempt to explain this striking finding we studied, by computer simulation, conformational properties of supercoiled DNA catenanes. Although the simulation showed that conformational properties of (+) and (-) supercoiled replication catenanes are very different, these properties per se do not give any advantage to (+) supercoiled over (-) supercoiled DNA catenanes for unlinking. An advantage became evident, however, when we took into account the established features of the enzymatic reaction catalyzed by the topoisomerases. The enzymes create a sharp DNA bend in the first bound DNA segment and allow for the transport of the second segment only from inside the bend to its outside. We showed that in (-) supercoiled DNA catenanes this protein-bound bent segment becomes nearly inaccessible for segments of the other linked DNA molecule, inhibiting the unlinking.

Theoretical models of DNA topology simplification by type IIA DNA topoisomerases.

It was discovered 12 years ago that type IIA topoisomerases can simplify DNA topology--the steady-state fractions of knots and links created by the enzymes are many times lower than the corresponding equilibrium fractions. Though this property of the enzymes made clear biological sense, it was not clear how small enzymes could selectively change the topology of very large DNA molecules, since topology is a global property and cannot be determined by a local DNA-protein interaction. A few models, suggested to explain the phenomenon, are analyzed in this review. We also consider experimental data that both support and contravene these models.

Kinking the double helix by bending deformation.

DNA bending and torsional deformations that often occur during its functioning inside the cell can cause local disruptions of the regular helical structure. The disruptions created by negative torsional stress have been studied in detail, but those caused by bending stress have only been analyzed theoretically. By probing the structure of very small DNA circles, we determined that bending stress disrupts the regular helical structure when the radius of DNA curvature is smaller than 3.5 nm. First, we developed an efficient method to obtain covalently closed DNA minicircles. To detect structural disruptions in the minicircles we treated them by single-strand-specific endonucleases. The data showed that the regular DNA structure is disrupted by bending deformation in the 64-65-bp minicircles, but not in the 85-86-bp minicircles. Our results suggest that strong DNA bending initiates kink formation while preserving base pairing.

Determining protein-induced DNA bending in force-extension experiments: theoretical analysis.

Computer simulations were used to investigate the possibility of determining protein-induced DNA bend angles by measuring the extension of a single DNA molecule. Analysis of the equilibrium sets of DNA conformations showed that shortening of DNA extension by a single protein-induced DNA bend can be as large as 35 nm. The shortening has a maximum value at the extending force of approximately 0.1 pN. At this force, the DNA extension experiences very large fluctuations that dramatically complicate the measurement. Using Brownian dynamics simulation of a DNA molecule extended by force, we were able to estimate the observation time needed to obtain the desired accuracy of the extension measurement. Also, the simulation revealed large fluctuations of the force, acting on the attached magnetic bead from the stretched DNA molecule.

Theoretical analysis of disruptions in DNA minicircles.

Under sufficient bending stress, which appears in DNA minicircles and small DNA loops, the double helix experiences local disruptions of its regular structure. We developed a statistical-mechanical treatment of the disruptions in DNA minicircles, studied experimentally by Du et al. The model of disruptions used in our Monte Carlo simulation of minicircle conformations specifies these conformations by three parameters: DNA bend angle at the disruption, theta(d); local DNA unwinding caused by the disruption; and the free energy associated with the disruption in the unstressed double helix, G(d). The model is applicable to any structural type of disruption, kinks or opening of single basepairs. The simulation shows that accounting for both torsional and bending deformation associated with the disruptions is very important for proper analysis. We obtained a relationship between values of G(d) and theta(d) under which the simulation results are compatible with the experimental data. The relationship suggests that the free energy of basepair opening, which includes flipping out both bases, is significantly higher than the generally accepted value. The model is also applied to the analysis of j-factors of very short DNA fragments.

Simulation of DNA catenanes.

DNA catenanes are important objects in biology, foremost as they appear during replication of circular DNA molecules. In this review we analyze how conformational properties of DNA catenanes can be studied by computer simulation. We consider classification of catenanes, their topological invariants and the methods of calculation of these invariants. We briefly analyze the DNA model and the simulation procedure used to sample the equilibrium conformational ensemble of catenanes with a particular topology. We consider how to avoid direct simulation of many DNA molecules when we need to account for the linking-unlinking process. The simulation methods and their comparisons with experiments are illustrated by some examples. We also describe an approach that allows simulating the steady state fraction of DNA catenanes created by type II topoisomerases.

Protein-induced local DNA bends regulate global topology of recombination products.

The tyrosine family of recombinases produces two smaller DNA circles when acting on circular DNA harboring two recombination sites in head-to-tail orientation. If the substrate is supercoiled, these circles can be unlinked or form multiply linked catenanes. The topological complexity of the products varies strongly even for similar recombination systems. This dependence has been solved here. Our computer simulation of the synapsis showed that the bend angles, phi, created in isolated recombination sites by protein binding before assembly of the full complex, determine the product topology. To verify the validity of this theoretical finding we measured the values of phi for Cre/loxP and Flp/FRT systems. The measurement was based on cyclization of the protein-bound short DNA fragments in solution. Despite the striking similarity of the synapses for these recombinases, action of Cre on head-to-tail target sites produces mainly unlinked circles, while that of Flp yields multiply linked catenanes. In full agreement with theoretical expectations we found that the values of phi for these systems are very different, close to 35 degrees and 80 degrees, respectively. Our findings have general implications in how small protein machines acting locally on large DNA molecules exploit statistical properties of their substrates to bring about directed global changes in topology.

DNA melting and energetics of the double helix.

Studying melting and energetics of the DNA double helix has been one of the major topics of molecular biophysics over the past six decades. The main objective of this article is to overview the current state of the field and to emphasize that there are still serious gaps in our understanding of the issue. We start with a concise description of the commonly accepted theoretical model of the DNA melting. We then concentrate on studies devoted to the comparison with experiment of theoretically predicted melting profiles of DNAs with known sequences. For long DNA molecules, such comparison is significant from the basic-science viewpoint while an accurate theoretical description of melting of short duplexes is necessary for various very important applications in biotechnology. Several sets of DNA melting parameters, proposed within the framework of the nearest neighbor model, are compared and analyzed. The analysis leads to a conclusion that in case of long DNA molecules the consensus set of nearest neighbor parameters describes well the experimental melting profiles. Unexpectedly, for short DNA duplexes the same set of parameters hardly yields any improvement as compared to the simplest model, which completely ignores the effect of heterogeneous stacking. Possible causes of this striking observation are discussed. We then overview the issue of separation of base-pairing and base-stacking contributions into the double helix stability. The recent experimental attempts to solve the problem are extensively analyzed. It is concluded that the double helix is essentially stabilized by stacking interaction between adjacent base pairs. Base pairing between complementary pairs does not appreciably contribute into the duplex stability. In the final section of the article, kinetic aspects of the DNA melting phenomenon are discussed. The main emphasis is made on the hysteresis effects often observed in melting of long DNA molecules. It is argued that the phenomenon can be well described via an accurate theoretical treatment of the random-walk model of melting kinetics of an isolated helical segment in DNA.

Disentangling DNA molecules.

The widespread circular form of DNA molecules inside cells creates very serious topological problems during replication. Due to the helical structure of the double helix the parental strands of circular DNA form a link of very high order, and yet they have to be unlinked before the cell division. DNA topoisomerases, the enzymes that catalyze passing of one DNA segment through another, solve this problem in principle. However, it is very difficult to remove all entanglements between the replicated DNA molecules due to huge length of DNA comparing to the cell size. One strategy that nature uses to overcome this problem is to create the topoisomerases that can dramatically reduce the fraction of linked circular DNA molecules relative to the corresponding fraction at thermodynamic equilibrium. This striking property of the enzymes means that the enzymes that interact with DNA only locally can access their topology, a global property of circular DNA molecules. This review considers the experimental studies of the phenomenon and analyzes the theoretical models that have been suggested in attempts to explain it. We describe here how various models of enzyme action can be investigated computationally. There is no doubt at the moment that we understand basic principles governing enzyme action. Still, there are essential quantitative discrepancies between the experimental data and the theoretical predictions. We consider how these discrepancies can be overcome.

Contribution of the intrinsic curvature to measured DNA persistence length.

The persistence length of DNA, a, depends both on the intrinsic curvature of the double helix and on the thermal fluctuations of the angles between adjacent base-pairs. We have evaluated two contributions to the value of a by comparing measured values of a for DNA containing a generic sequence and for an "intrinsically straight" DNA. In each 10 bp segment of the intrinsically straight DNA an initial sequence of five bases is repeated in the sequence of the second five bases, so any bends in the first half of the segment are compensated by bends in the opposite direction in the second half. The value of a for the latter DNA depends, to a good approximation, on thermal fluctuations only; there is no intrinsic curvature. The values of a were obtained from measurements of the cyclization efficiency for short DNA fragments, about 200 bp in length. This method determines the persistence length of DNA with exceptional accuracy, due to the very strong dependence of the cyclization efficiency of short fragments on the value of a. We find that the values of a for the two types of DNA fragment are very close and conclude that the contribution of the intrinsic curvature to a is at least 20 times smaller than the contribution of thermal fluctuations. The relationship between this result and the angles between adjacent base-pairs, which specify the intrinsic curvature, is analyzed.

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.

Dynamics of DNA loop capture by the SfiI restriction endonuclease on supercoiled and relaxed DNA.

The SfiI endonuclease is a prototype for DNA looping. It binds two copies of its recognition sequence and, if Mg(2+) is present, cuts both concertedly. Looping was examined here on supercoiled and relaxed forms of a 5.5 kb plasmid with three SfiI sites: sites 1 and 2 were separated by 0.4 kb, and sites 2 and 3 by 2.0 kb. SfiI converted this plasmid directly to the products cut at all three sites, though DNA species cleaved at one or two sites were formed transiently during a burst phase. The burst revealed three sets of doubly cut products, corresponding to the three possible pairings of sites. The equilibrium distribution between the different loops was evaluated from the burst phases of reactions initiated by adding MgCl(2) to SfiI bound to the plasmid. The short loop was favored over the longer loops, particularly on supercoiled DNA. The relative rates for loop capture were assessed after adding SfiI to solutions containing the plasmid and MgCl(2). On both supercoiled and relaxed DNA, the rate of loop capture across 0.4 kb was only marginally faster than over 2.0 kb or 2.4 kb. The relative strengths and rates of looping were compared to computer simulations of conformational fluctuations in DNA. The simulations concurred broadly with the experimental data, though they predicted that increasing site separations should cause a shallower decline in the equilibrium constants than was observed but a slightly steeper decline in the rates for loop capture. Possible reasons for these discrepancies are discussed.

Symmetric DNA sites are functionally asymmetric within Flp and Cre site-specific DNA recombination synapses.

Flp and Cre-mediated recombination on symmetrized FRT and loxP sites, respectively, in circular plasmid substrates yield both DNA inversion and deletion. However, upon sequestering three negative supercoils outside the recombination complex using the resII-resIII synapse formed by Tn3 resolvase and the LER synapse formed by phage Mu transposase in the case of Flp and Cre, respectively, the reactions are channeled towards inversion at the expense of deletion. The inversion product is a trefoil, its unique topology being conferred by the external resolvase or LER synapse. Thus, Flp and Cre assign their symmetrized substrates a strictly antiparallel orientation with respect to strand cleavage and exchange. These conclusions are supported by the product profiles from tethered parallel and antiparallel native FRT sites in dilution and competition assays. Furthermore, the observed recombination bias favoring deletion over inversion in a nicked circular substrate containing two symmetrized FRT sites is consistent with the predictions from Monte Carlo simulations based on antiparallel synapsis of the DNA partners.

Topological Behavior of Plasmid DNA.

The discovery of the B-form structure of DNA by Watson and Crick led to an explosion of research on nucleic acids in the fields of biochemistry, biophysics, and genetics. Powerful techniques were developed to reveal a myriad of different structural conformations that change B-DNA as it is transcribed, replicated, and recombined and as sister chromosomes are moved into new daughter cell compartments during cell division. This article links the original discoveries of superhelical structure and molecular topology to non-B form DNA structure and contemporary biochemical and biophysical techniques. The emphasis is on the power of plasmids for studying DNA structure and function. The conditions that trigger the formation of alternative DNA structures such as left-handed Z-DNA, inter- and intra-molecular triplexes, triple-stranded DNA, and linked catenanes and hemicatenanes are explained. The DNA dynamics and topological issues are detailed for stalled replication forks and for torsional and structural changes on DNA in front of and behind a transcription complex and a replisome. The complex and interconnected roles of topoisomerases and abundant small nucleoid association proteins are explained. And methods are described for comparing in vivo and in vitro reactions to probe and understand the temporal pathways of DNA and chromosome chemistry that occur inside living cells.