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

Molecular Characterization and Genome Mechanical Features of Two Newly Isolated Polyvalent Bacteriophages Infecting Pseudomonas syringae pv. garcae

Coffee plants have been targeted by a devastating bacterial disease, a condition known as bacterial blight, caused by the phytopathogen Pseudomonas syringae pv. garcae (Psg). Conventional treatments of coffee plantations affected by the disease involve frequent spraying with copper- and kasugamycin-derived compounds, but they are both highly toxic to the environment and stimulate the appearance of bacterial resistance. Herein, we report the molecular characterization and mechanical features of the genome of two newly isolated (putative polyvalent) lytic phages for Psg. The isolated phages belong to class Caudoviricetes and present a myovirus-like morphotype belonging to the genuses Tequatrovirus (PsgM02F) and Phapecoctavirus (PsgM04F) of the subfamilies Straboviridae (PsgM02F) and Stephanstirmvirinae (PsgM04F), according to recent bacterial viruses' taxonomy, based on their complete genome sequences. The 165,282 bp (PsgM02F) and 151,205 bp (PsgM04F) genomes do not feature any lysogenic-related (integrase) genes and, hence, can safely be assumed to follow a lytic lifestyle. While phage PsgM02F produced a morphogenesis yield of 124 virions per host cell, phage PsgM04F produced only 12 virions per host cell, indicating that they replicate well in Psg with a 50 min latency period. Genome mechanical analyses established a relationship between genome bendability and virion morphogenesis yield within infected host cells.

A Primer on the oxDNA Model of DNA: When to Use it, How to Simulate it and How to Interpret the Results

The oxDNA model of Deoxyribonucleic acid has been applied widely to systems in biology, biophysics and nanotechnology. It is currently available via two independent open source packages. Here we present a set of clearly documented exemplar simulations that simultaneously provide both an introduction to simulating the model, and a review of the model's fundamental properties. We outline how simulation results can be interpreted in terms of-and feed into our understanding of-less detailed models that operate at larger length scales, and provide guidance on whether simulating a system with oxDNA is worthwhile.

Understanding the paradoxical mechanical response of in-phase A-tracts at different force regimes

A-tracts are A:T rich DNA sequences that exhibit unique structural and mechanical properties associated with several functions in vivo. The crystallographic structure of A-tracts has been well characterized. However, the mechanical properties of these sequences is controversial and their response to force remains unexplored. Here, we rationalize the mechanical properties of in-phase A-tracts present in the Caenorhabditis elegans genome over a wide range of external forces, using single-molecule experiments and theoretical polymer models. Atomic Force Microscopy imaging shows that A-tracts induce long-range (∼200 nm) bending, which originates from an intrinsically bent structure rather than from larger bending flexibility. These data are well described with a theoretical model based on the worm-like chain model that includes intrinsic bending. Magnetic tweezers experiments show that the mechanical response of A-tracts and arbitrary DNA sequences have a similar dependence with monovalent salt supporting that the observed A-tract bend is intrinsic to the sequence. Optical tweezers experiments reveal a high stretch modulus of the A-tract sequences in the enthalpic regime. Our work rationalizes the complex multiscale flexibility of A-tracts, providing a physical basis for the versatile character of these sequences inside the cell.

Bent DNA Bows as Sensing Amplifiers for Detecting DNA-Interacting Salts and Molecules

Due to the central role of DNA, its interactions with inorganic salts and small organic molecules are important. For example, such interactions play important roles in various fundamental cellular processes in living systems and are involved in many DNA-damage related diseases. Strategies to improve the sensitivity of existing techniques for studying DNA interactions with other molecules would be appreciated in situations where the interactions are too weak. Here we report our development and demonstration of bent DNA bows for amplifying, sensing, and detecting the interactions of 14 inorganic salts and small organic molecules with DNA. With the bent DNA bows, these interactions were easily visualized and quantified in gel electrophoresis, which were difficult to measure without bending. In addition, the strength of the interactions of DNA with the various salts/molecules were quantified using the modified Hill equation. This work highlights the amplification effects of the bending elastic energy stored in the DNA bows and the potential use of the DNA bows for quantitatively measuring DNA interactions with small molecules as simple economic methods; it may also pave the way for exploiting the bent DNA bows for other applications such as screening DNA-interacting molecules and drugs.

Identifying Physical Causes of Apparent Enhanced Cyclization of Short DNA Molecules with a Coarse-Grained Model

DNA cyclization is a powerful technique to gain insight into the nature of DNA bending. While the wormlike chain model provides a good description of small to moderate bending fluctuations, it is expected to break down for large bending. Recent cyclization experiments on strongly bent shorter molecules indeed suggest enhanced flexibility over and above that expected from the wormlike chain. Here, we use a coarse-grained model of DNA to investigate the subtle thermodynamics of DNA cyclization for molecules ranging from 30 to 210 base pairs. As the molecules get shorter, we find increasing deviations between our computed equilibrium j-factor and the classic wormlike chain predictions of Shimada and Yamakawa for a torsionally aligned looped molecule. These deviations are due to sharp kinking, first at nicks, and only subsequently in the body of the duplex. At the shortest lengths, substantial fraying at the ends of duplex domains is the dominant method of relaxation. We also estimate the dynamic j-factor measured in recent FRET experiments. We find that the dynamic j-factor is systematically larger than its equilibrium counterpart—with the deviation larger for shorter molecules—because not all the stress present in the fully cyclized state is present in the transition state. These observations are important for the interpretation of recent cyclization experiments, suggesting that measured anomalously high j-factors may not necessarily indicate non-WLC behavior in the body of duplexes.

Structural waters in the minor and major grooves of DNA--a major factor governing structural adjustments of the A-T mini-helix

The role of microhydration in structural adjustments of the AT-tract in B-DNA was studied at the B97-D/def2-SV(P) level. The (dA:dT)5 complexes with 10 water molecules in minor and 15 water molecules in major grooves were studied. The obtained network of hydrogen bonds revealed the dependence between the groove width and the types of water patterns. In the minor groove, the following patterns were observed: interstrand one-water bridges similar to that of the Dickerson "water spine" and interstrand two-water bridges. The network of structural waters in the major groove is more diverse than that in the minor groove, which agrees with crystallographic data. As the major groove is wider, it is enriched by water molecules forming two- and three-water bridges. Results suggest the nucleobase-water interactions in both grooves prevent AT-tract twisting and its "collapse" along the minor groove. Whereby, a helix structure with narrow minor and wide major grooves is formed. The structural waters affect the polynucleotide conformation so that it becomes similar to poly(dA)·poly(dT) in fibers and acquires features of the A-tracts in DNA in solution. We suggest that formation of specific water patterns in both grooves is the factor responsible for stabilization of A-tracts with a narrowed minor groove, leading in turn to their strong intrinsic bending in DNA.

Mechanical properties of symmetric and asymmetric DNA A-tracts: implications for looping and nucleosome positioning

A-tracts are functionally important DNA sequences which induce helix bending and have peculiar structural properties. While A-tract structure has been qualitatively well characterized, their mechanical properties remain controversial. A-tracts appear structurally rigid and resist nucleosome formation, but seem flexible in DNA looping. In this work, we investigate mechanical properties of symmetric AnTn and asymmetric A2n tracts for n = 3, 4, 5 using two types of coarse-grained models. The first model represents DNA as an ensemble of interacting rigid bases with non-local quadratic deformation energy, the second one treats DNA as an anisotropically bendable and twistable elastic rod. Parameters for both models are inferred from microsecond long, atomic-resolution molecular dynamics simulations. We find that asymmetric A-tracts are more rigid than the control G/C-rich sequence in localized distortions relevant for nucleosome formation, but are more flexible in global bending and twisting relevant for looping. The symmetric tracts, in contrast, are more rigid than asymmetric tracts and the control, both locally and globally. Our results can reconcile the contradictory stiffness data on A-tracts and suggest symmetric A-tracts to be more efficient in nucleosome exclusion than the asymmetric ones. This would open a new possibility of gene expression manipulation using A-tracts.

Mathematical model of alternative mechanism of telomere length maintenance

Biopolymer length regulation is a complex process that involves a large number of biological, chemical, and physical subprocesses acting simultaneously across multiple spatial and temporal scales. An illustrative example important for genomic stability is the length regulation of telomeres-nucleoprotein structures at the ends of linear chromosomes consisting of tandemly repeated DNA sequences and a specialized set of proteins. Maintenance of telomeres is often facilitated by the enzyme telomerase but, particularly in telomerase-free systems, the maintenance of chromosomal termini depends on alternative lengthening of telomeres (ALT) mechanisms mediated by recombination. Various linear and circular DNA structures were identified to participate in ALT, however, dynamics of the whole process is still poorly understood. We propose a chemical kinetics model of ALT with kinetic rates systematically derived from the biophysics of DNA diffusion and looping. The reaction system is reduced to a coagulation-fragmentation system by quasi-steady-state approximation. The detailed treatment of kinetic rates yields explicit formulas for expected size distributions of telomeres that demonstrate the key role played by the J factor, a quantitative measure of bending of polymers. The results are in agreement with experimental data and point out interesting phenomena: an appearance of very long telomeric circles if the total telomere density exceeds a critical value (excess mass) and a nonlinear response of the telomere size distributions to the amount of telomeric DNA in the system. The results can be of general importance for understanding dynamics of telomeres in telomerase-independent systems as this mode of telomere maintenance is similar to the situation in tumor cells lacking telomerase activity. Furthermore, due to its universality, the model may also serve as a prototype of an interaction between linear and circular DNA structures in various settings.

Bending of short DNA helices

In their recent Science paper, Vafabakhsh and Ha claim that DNA duplexes at the range of 100 bp experience anomalous flexibility, much greater than the flexibility of large DNA molecules. ( 1) However, careful reevaluation of their data leads to the conclusion that the presented data do not warrant the authors' claim.

Impact of methylation on the physical properties of DNA

There is increasing evidence for the presence of an alternative code imprinted in the genome that might contribute to gene expression regulation through an indirect reading mechanism. In mammals, components of this coarse-grained regulatory mechanism include chromatin structure and epigenetic signatures, where d(CpG) nucleotide steps are key players. We report a comprehensive experimental and theoretical study of d(CpG) steps that provides a detailed description of their physical characteristics and the impact of cytosine methylation on these properties. We observed that methylation changes the physical properties of d(CpG) steps, having a dramatic effect on enriched CpG segments, such as CpG islands. We demonstrate that methylation reduces the affinity of DNA to assemble into nucleosomes, and can affect nucleosome positioning around transcription start sites. Overall, our results suggest a mechanism by which the basic physical properties of the DNA fiber can explain parts of the cellular epigenetic regulatory mechanisms.

Bullied no more: when and how DNA shoves proteins around

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

Quantitative methods for measuring DNA flexibility in vitro and in vivo

The double-helical DNA biopolymer is particularly resistant to bending and twisting deformations. This property has important implications for DNA folding in vitro and for the packaging and function of DNA in living cells. Among the outstanding questions in the field of DNA biophysics are the underlying origin of DNA stiffness and the mechanisms by which DNA stiffness is overcome within cells. Exploring these questions requires experimental methods to quantitatively measure DNA bending and twisting stiffness both in vitro and in vivo. Here, we discuss two classical approaches: T4 DNA ligase-mediated DNA cyclization kinetics and lac repressor-mediated DNA looping in Escherichia coli. We review the theoretical basis for these techniques and how each can be applied to quantitate biophysical parameters that describe the DNA polymer. We then show how we have modified these methods and applied them to quantitate how apparent DNA physical properties are altered in vitro and in vivo by sequence-nonspecific architectural DNA-binding proteins such as the E. coli HU protein and eukaryotic HMGB proteins.

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 curvature and flexibility in vitro and in vivo

It has been more than 50 years since the elucidation of the structure of double-helical DNA. Despite active research and progress in DNA biology and biochemistry, much remains to be learned in the field of DNA biophysics. Predicting the sequence-dependent curvature and flexibility of DNA is difficult. Applicability of the conventional worm-like chain polymer model of DNA has been challenged. The fundamental forces responsible for the remarkable resistance of DNA to bending and twisting remain controversial. The apparent 'softening' of DNA measured in vivo in the presence of kinking proteins and superhelical strain is incompletely understood. New methods and insights are being applied to these problems. This review places current work on DNA biophysics in historical context and illustrates the ongoing interplay between theory and experiment in this exciting field.

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.

First-principles calculation of DNA looping in tethered particle experiments

We calculate the probability of DNA loop formation mediated by regulatory proteins such as Lac repressor (LacI), using a mathematical model of DNA elasticity. Our model is adapted to calculating quantities directly observable in tethered particle motion (TPM) experiments, and it accounts for all the entropic forces present in such experiments. Our model has no free parameters; it characterizes DNA elasticity using information obtained in other kinds of experiments. It assumes a harmonic elastic energy function (or wormlike chain type elasticity), but our Monte Carlo calculation scheme is flexible enough to accommodate arbitrary elastic energy functions. We show how to compute both the 'looping J factor' (or equivalently, the looping free energy) for various DNA construct geometries and LacI concentrations, as well as the detailed probability density function of bead excursions. We also show how to extract the same quantities from recent experimental data on TPM, and then compare to our model's predictions. In particular, we present a new method to correct observed data for finite camera shutter time and other experimental effects. Although the currently available experimental data give large uncertainties, our first-principles predictions for the looping free energy change are confirmed to within about 1 k(B)T, for loops of length around 300 basepairs. More significantly, our model successfully reproduces the detailed distributions of bead excursion, including their surprising three-peak structure, without any fit parameters and without invoking any alternative conformation of the LacI tetramer. Indeed, the model qualitatively reproduces the observed dependence of these distributions on tether length (e.g., phasing) and on LacI concentration (titration). However, for short DNA loops (around 95 basepairs) the experiments show more looping than is predicted by the harmonic-elasticity model, echoing other recent experimental results. Because the experiments we study are done in vitro, this anomalously high looping cannot be rationalized as resulting from the presence of DNA-bending proteins or other cellular machinery. We also show that it is unlikely to be the result of a hypothetical 'open' conformation of the LacI tetramer.

Poly(dA:dT) tracts: major determinants of nucleosome organization

Homopolymeric stretches of deoxyadenosine nucleotides (A's) on one strand of double-stranded DNA, referred to as poly(dA:dT) tracts or A-tracts, are overabundant in eukaryotic genomes. They have unusual structural, dynamic, and mechanical properties, and may resist sharp bending. Such unusual material properties, together with their overabundance in eukaryotes, raised the possibility that poly(dA:dT) tracts might function in eukaryotes to influence the organization of nucleosomes at many genomic regions. Recent genome-wide studies strongly confirm these ideas and suggest that these tracts play major roles in chromatin organization and genome function. Here we review what is known about poly(dA:dT) tracts and how they work.

Computational analysis of looping of a large family of highly bent DNA by LacI

Sequence-dependent intrinsic curvature of DNA influences looping by regulatory proteins such as LacI and NtrC. Curvature can enhance stability and control shape, as observed in LacI loops formed with three designed sequences with operators bracketing an A-tract bend. We explore geometric, topological, and energetic effects of curvature with an analysis of a family of highly bent sequences, using the elastic rod model from previous work. A unifying straight-helical-straight representation uses two phasing parameters to describe sequences composed of two straight segments that flank a common helically supercoiled segment. We exercise the rod model over this two-dimensional space of phasing parameters to evaluate looping behaviors. This design space is found to comprise two subspaces that prefer parallel versus anti-parallel binding topologies. The energetic cost of looping varies from 4 to 12 kT. Molecules can be designed to yield distinct binding topologies as well as hyperstable or hypostable loops and potentially loops that can switch conformations. Loop switching could be a mechanism for control of gene expression. Model predictions for linking numbers and sizes of LacI-DNA loops can be tested using multiple experimental approaches, which coupled with theory could address whether proteins or DNA provide the observed flexibility of protein-DNA loops.

Effects of the nucleoid protein HU on the structure, flexibility, and ring-closure properties of DNA deduced from Monte Carlo simulations

The histone-like HU (heat unstable) protein plays a key role in the organization and regulation of the Escherichia coli genome. The nonspecific nature of HU binding to DNA complicates analysis of the mechanism by which the protein contributes to the looping of DNA. Conventional models of the looping of HU-bound duplexes attribute the changes in biophysical properties of DNA brought about by the random binding of protein to changes in the effective parameters of an ideal helical wormlike chain. Here, we introduce a novel Monte Carlo approach to study the effects of nonspecific HU binding on the configurational properties of DNA directly. We randomly decorated segments of an ideal double-helical DNA with HU molecules that induce the bends and other structural distortions of the double helix find in currently available X-ray structures. We find that the presence of HU at levels approximating those found in the cell reduces the persistence length by roughly threefold compared with that of naked DNA. The binding of protein has particularly striking effects on the cyclization properties of short duplexes, altering the dependence of ring closure on chain length in a way that cannot be mimicked by a simple wormlike model and accumulating at higher-than-expected levels on successfully closed chains. Moreover, the uptake of protein on small minicircles depends on chain length, taking advantage of the HU-induced deformations of DNA structure to facilitate ligation. Circular duplexes with bound HU show much greater propensity than protein-free DNA to exist as negatively supercoiled topoisomers, suggesting a potential role of HU in organizing the bacterial nucleoid. The local bending and undertwisting of DNA by HU, in combination with the number of bound proteins, provide a structural rationale for the condensation of DNA and the observed expression levels of reporter genes in vivo.

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.

Key experimental approaches in DNA nanotechnology

DNA nanotechnology combines unusual DNA motifs with sticky-ended cohesion to build polyhedral objects, topological targets, nanomechanical devices, and both crystalline and aperiodic arrays. The goal of DNA nanotechnology is control of the structure of macroscopic matter on the finest possible scale. Applications are expected to arise in the areas of X-ray crystallography, nanoelectronics, nanorobotics, and DNA-based computation. DNA and its close molecular relatives appear extremely well suited for these goals. This overview covers the generation of new DNA motifs, construction methods (synthesis, hybridization, phosphorylation, ligation), and a variety of methods for characterization of motifs, devices, and arrays. Finally, the use of DNA nanotechnology as a tool in biochemistry is discussed.

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.

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.

Template-free, polymerase-free DNA polymerization

A combination of a DNA ligase and a restriction endonuclease provides a DNA polymerase activity, which might suggest a novel strategy for polymer synthesis.

Cyclization of short DNA fragments and bending fluctuations of the double helix

Cloutier and Widom [Cloutier, T. E. & Widom, J. (2004) Mol. Cell 14, 355-362] recently reported that the cyclization efficiency of short DNA fragments, about 100 bp in length, exceeds theoretical expectations by three orders of magnitude. In an effort to resolve this discrepancy, we tried modifying the theory. We investigated how the distribution of the angles between adjacent base pairs of the double helix affects the cyclization efficiency. We found that only the incorporation of sharp kinks in the angle distribution provides the desired increase of the cyclization efficiency. We did not find a model, however, that fits all cyclization data for DNA fragments of different lengths. Therefore, we carefully reinvestigated the cyclization of 100-bp DNA fragments experimentally and found their cyclization efficiency to be in remarkable agreement with the traditional model of DNA bending. We also found an explanation for the discrepancy between our results and those of Cloutier and Widom.

Gapped DNA and cyclization of short DNA fragments

We use the cyclization of small DNA molecules, approximately 200 bp in length, to study conformational properties of DNA fragments with single-stranded gaps. The approach is extremely sensitive to DNA conformational properties and, being complemented by computations, allows a very accurate determination of the fragment's conformational parameters. Sequence-specific nicking endonucleases are used to create the 4-nt-long gap. We determined the bending rigidity of the single-stranded region in the gapped DNA. We found that the gap of 4 nt in length makes all torsional orientations of DNA ends equally probable. Our results also show that the gap has isotropic bending rigidity. This makes it very attractive to use gapped DNA in the cyclization experiments to determine DNA conformational properties, since the gap eliminates oscillations of the cyclization efficiency with the DNA length. As a result, the number of measurements is greatly reduced in the approach, and the analysis of the data is greatly simplified. We have verified our approach on DNA fragments containing well-characterized intrinsic bends caused by A-tracts. The obtained experimental results and theoretical analysis demonstrate that gapped-DNA cyclization is an exceedingly sensitive and accurate approach for the determination of DNA bending.

Intrinsic curvature in the VP1 gene of SV40: comparison of solution and gel results

DNA restriction fragments that are stably curved are usually identified by polyacrylamide gel electrophoresis because curved fragments migrate more slowly than normal fragments containing the same number of basepairs. In free solution, curved DNA molecules can be identified by transient electric birefringence (TEB) because they exhibit rotational relaxation times that are faster than those of normal fragments of the same size. In this article, the results observed in free solution and in polyacrylamide gels are compared for a highly curved 199-basepair (bp) restriction fragment taken from the VP1 gene in Simian Virus 40 (SV40) and various sequence mutants and insertion derivatives. The TEB method of overlapping fragments was used to show that the 199-bp fragment has an apparent bend angle of 46 +/- 2 degrees centered at sequence position 1922 +/- 2 bp. Four unphased A- and T-tracts and a mixed A3T4-tract occur within a span of approximately 60 bp surrounding the apparent bend center; for brevity, this 60-bp sequence element is called a curvature module. Modifying any of the A- or T-tracts in the curvature module by site-directed mutagenesis decreases the curvature of the fragment; replacing all five A- and T-tracts by random-sequence DNA causes the 199-bp mutant to adopt a normal conformation, with normal electrophoretic mobilities and birefringence relaxation times. Hence, stable curvature in this region of the VP1 gene is due to the five unphased A- and T- tracts surrounding the apparent bend center. Discordant solution and gel results are observed when long inverted repeats are inserted within the curvature module. These insertion derivatives migrate anomalously slowly in polyacrylamide gels but have normal, highly flexible conformations in free solution. Discordant solution and gel results are not observed if the insert does not contain a long inverted repeat or if the long inverted repeat is added to the 199-bp fragment outside the curvature module. The results suggest that long inverted repeats can form hairpins or cruciforms when they are located within a region of the helix backbone that is intrinsically curved, leading to large mobility anomalies in polyacrylamide gels. Hairpin/cruciform formation is not observed in free solution, presumably because of rapid conformational exchange. Hence, DNA restriction fragments that migrate anomalously slowly in polyacrylamide gels are not necessarily stably curved in free solution.

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.

Inducing and modulating anisotropic DNA bends by pseudocomplementary peptide nucleic acids

DNA bending is significant for various DNA functions in the cell. Here, we demonstrate that pseudocomplementary peptide nucleic acids (pcPNAs) represent a class of versatile, sequence-specific DNA-bending agents. The occurrence of anisotropic DNA bends induced by pcPNAs is shown by gel electrophoretic phasing analysis. The magnitude of DNA bending is determined by circular permutation assay and by electron microscopy, with good agreement of calculated mean values between both methods. Binding of a pair of 10-meric pcPNAs to its target DNA sequence results in moderate DNA bending with a mean value of 40-45 degrees, while binding of one self-pc 8-mer PNA to target DNA yields a somewhat larger average value of the induced DNA bend. Both bends are found to be in phase when the pcPNA target sites are separated by distances of half-integer numbers of helical turns of regular duplex DNA, resulting in an enhanced DNA bend with an average value in the range of 80-90 degrees. The occurrence of such a sharp bend within the DNA double helix is confirmed and exploited through efficient formation of 170-bp-long DNA minicircles by means of dimerization of two bent DNA fragments. The pcPNAs offer two main advantages over previously designed classes of nonnatural DNA-bending agents: they have very mild sequence limitations while targeting duplex DNA and they can easily be designed for a chosen target sequence, because their binding obeys the principle of complementarity. We conclude that pcPNAs are promising tools for inducing bends in DNA at virtually any chosen site.

The flexibility of DNA double crossover molecules

Double crossover molecules are DNA structures containing two Holliday junctions connected by two double helical arms. There are several types of double crossover molecules, differentiated by the relative orientations of their helix axes, parallel or antiparallel, and by the number of double helical half-turns (even or odd) between the two crossovers. They are found as intermediates in meiosis and they have been used extensively in structural DNA nanotechnology for the construction of one-dimensional and two-dimensional arrays and in a DNA nanomechanical device. Whereas the parallel double helical molecules are usually not well behaved, we have focused on the antiparallel molecules; antiparallel molecules with an even number of half-turns between crossovers (termed DAE molecules) produce a reporter strand when ligated, facilitating their characterization in a ligation cyclization assay. Hence, we have estimated the flexibility of antiparallel DNA double crossover molecules by means of ligation-closure experiments. We are able to show that these molecules are approximately twice as rigid as linear duplex DNA.

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.

Experimental evaluation of the Liu-Beveridge dinucleotide step model of DNA structure

Methods for predicting DNA curvature have many possible applications. Dinucleotide step models describe DNA shape by characterization of helical twist, deflection angles and the direction of deflection for nearest neighbor base pairs. Liu and Beveridge have extended previous applications of dinucleotide step models with the development and qualitative validation of a predictive method for sequence-dependent DNA curvature (the LB model). We tested whether the LB model accurately predicts experimentally deduced curvature angles and helical repeat parameters for DNA sequences not in its training set, particularly when challenged with quantitative data and subtle sequence phasings. We examined a series of 17 well-characterized DNA sequences to compare electrophoretic and computational results. The LB model is superior to two other models in the prediction of helical repeat parameters. We observed a strong linear correlation between curvature magnitudes predicted using the LB model and those determined by electrophoretic ligation ladder experiments, although the LB model somewhat underestimated apparent curvature. With longer electrophoretic phasing probes the LB model slightly overestimated gel mobility anomalies, with modest deviations in predicted helical repeat parameters. Overall, our analyses suggest that the LB model provides reasonably accurate predictions for the electrophoretic behavior of DNA.

Mechanism of topology simplification by type II DNA topoisomerases

Type II DNA topoisomerases actively reduce the fractions of knotted and catenated circular DNA below thermodynamic equilibrium values. To explain this surprising finding, we designed a model in which topoisomerases introduce a sharp bend in DNA. Because the enzymes have a specific orientation relative to the bend, they act like Maxwell's demon, providing unidirectional strand passage. Quantitative analysis of the model by computer simulations proved that it can explain much of the experimental data. The required sharp DNA bend was demonstrated by a greatly increased cyclization of short DNA fragments from topoisomerase binding and by direct visualization with electron microscopy.