Statistical analysis of DNA duplex structural features

The role of structural parameters in DNA cyclization

BACKGROUND

The intrinsic bendability of DNA plays an important role with relevance for myriad of essential cellular mechanisms. The flexibility of a DNA fragment can be experimentally and computationally examined by its propensity for cyclization, quantified by the Jacobson-Stockmayer J factor. In this study, we use a well-established coarse-grained three-dimensional model of DNA and seven distinct sets of experimentally and computationally derived conformational parameters of the double helix to evaluate the role of structural parameters in calculating DNA cyclization.

RESULTS

We calculate the cyclization rates of 86 DNA sequences with previously measured J factors and lengths between 57 and 325 bp as well as of 20,000 randomly generated DNA sequences with lengths between 350 and 4000 bp. Our comparison with experimental data is complemented with analysis of simulated data.

CONCLUSIONS

Our data demonstrate that all sets of parameters yield very similar results for longer DNA fragments, regardless of the nucleotide sequence, which are in agreement with experimental measurements. However, for DNA fragments shorter than 100 bp, all sets of parameters performed poorly yielding results with several orders of magnitude difference from the experimental measurements. Our data show that DNA cyclization rates calculated using conformational parameters based on nucleosome packaging data are most similar to the experimental measurements. Overall, our study provides a comprehensive large-scale assessment of the role of structural parameters in calculating DNA cyclization rates.

Evaluation of DNA bending models in their capacity to predict electrophoretic migration anomalies of satellite DNA sequences

DNA containing a sequence that generates a local curvature exhibits a pronounced retardation in electrophoretic mobility. Various theoretical models have been proposed to explain relationship between DNA structural features and migration anomaly. Here, we studied the capacity of 15 static wedge-bending models to predict electrophoretic behavior of 69 satellite monomers derived from four divergent families. All monomers exhibited retarded mobility in PAGE corresponding to retardation factors ranging 1.02-1.54. The curvature varied both within and across the groups and correlated with the number, position, and lengths of A-tracts. Two dinucleotide models provided strong correlation between gel mobility and curvature prediction; two trinucleotide models were satisfactory while remaining dinucleotide models provided intermediate results with reliable prediction for subsets of sequences only. In some cases, similarly shaped molecules exhibited relatively large differences in mobility and vice versa. Generally less accurate predictions were obtained in groups containing less homogeneous sequences possessing distinct structural features. In conclusion, relatively universal theoretical models were identified suitable for the analysis of natural sequences known to harbor relatively moderate curvature. These models could be potentially applied to genome wide studies. However, in silico predictions should be viewed in context of experimental measurement of intrinsic DNA curvature.

Molecular dynamics simulations of G-DNA and perspectives on the simulation of nucleic acid structures

The article reviews the application of biomolecular simulation methods to understand the structure, dynamics and interactions of nucleic acids with a focus on explicit solvent molecular dynamics simulations of guanine quadruplex (G-DNA and G-RNA) molecules. While primarily dealing with these exciting and highly relevant four-stranded systems, where recent and past simulations have provided several interesting results and novel insight into G-DNA structure, the review provides some general perspectives on the applicability of the simulation techniques to nucleic acids.

RNA structural motifs that entail hydrogen bonds involving sugar-phosphate backbone atoms of RNA

The growing number of high-resolution crystal structures of large RNA molecules provides much information for understanding the principles of structural organization of these complex molecules. Several in-depth analyses of nucleobase-centered RNA structural motifs and backbone conformations have been published based on this information, including a systematic classification of base pairs by Leontis and Westhof. However, hydrogen bonds involving sugar-phosphate backbone atoms of RNA have not been analyzed systematically until recently, although such hydrogen bonds appear to be common both in local and tertiary interactions. Here we review some backbone structural motifs discussed in the literature and analyze a set of eight high-resolution multi-domain RNA structures. The analyzed RNAs are highly structured: among 5372 nucleotides in this set, 89% are involved in at least one "long-range" RNA-RNA hydrogen bond, i.e., hydrogen bonds between atoms in the same residue or sequential residues are ignored. These long-range hydrogen bonds frequently use backbone atoms as hydrogen bond acceptors, i.e., OP1, OP2, O2', O3', O4', or O5', or as a donor (2'OH). A surprisingly large number of such hydrogen bonds are found, considering that neither single-stranded nor double-stranded regions will contain such hydrogen bonds unless additional interactions with other residues exist. Among 8327 long-range hydrogen bonds found in this set of structures, 2811, or about one-third, are hydrogen bonds entailing RNA backbone atoms; they involve 39% of all nucleotides in the structures. The majority of them (2111) are hydrogen bonds entailing ribose hydroxyl groups, which can be used either as a donor or an acceptor; they constitute 25% of all hydrogen bonds and involve 31% of all nucleotides. The phosphate oxygens OP1 or OP2 are used as hydrogen bond acceptors in 12% of all nucleotides, and the ribose ring oxygen O4' and phosphodiester oxygens O3' and O5' are used in 4%, 4%, and 1% of all nucleotides, respectively. Distributions of geometric parameters and some examples of such hydrogen bonds are presented in this report. A novel motif involving backbone hydrogen bonds, the ribose-phosphate zipper, is also identified.

An iterative strategy combining biophysical criteria and duration hidden Markov models for structural predictions of Chlamydia trachomatis sigma66 promoters

BACKGROUND

Promoter identification is a first step in the quest to explain gene regulation in bacteria. It has been demonstrated that the initiation of bacterial transcription depends upon the stability and topology of DNA in the promoter region as well as the binding affinity between the RNA polymerase sigma-factor and promoter. However, promoter prediction algorithms to date have not explicitly used an ensemble of these factors as predictors. In addition, most promoter models have been trained on data from Escherichia coli. Although it has been shown that transcriptional mechanisms are similar among various bacteria, it is quite possible that the differences between Escherichia coli and Chlamydia trachomatis are large enough to recommend an organism-specific modeling effort.

RESULTS

Here we present an iterative stochastic model building procedure that combines such biophysical metrics as DNA stability, curvature, twist and stress-induced DNA duplex destabilization along with duration hidden Markov model parameters to model Chlamydia trachomatis sigma66 promoters from 29 experimentally verified sequences. Initially, iterative duration hidden Markov modeling of the training set sequences provides a scoring algorithm for Chlamydia trachomatis RNA polymerase sigma66/DNA binding. Subsequently, an iterative application of Stepwise Binary Logistic Regression selects multiple promoter predictors and deletes/replaces training set sequences to determine an optimal training set. The resulting model predicts the final training set with a high degree of accuracy and provides insights into the structure of the promoter region. Model based genome-wide predictions are provided so that optimal promoter candidates can be experimentally evaluated, and refined models developed. Co-predictions with three other algorithms are also supplied to enhance reliability.

CONCLUSION

This strategy and resulting model support the conjecture that DNA biophysical properties, along with RNA polymerase sigma-factor/DNA binding collaboratively, contribute to a sequence's ability to promote transcription. This work provides a baseline model that can evolve as new Chlamydia trachomatis sigma66 promoters are identified with assistance from the provided genome-wide predictions. The proposed methodology is ideal for organisms with few identified promoters and relatively small genomes.

Core promoters are predicted by their distinct physicochemical properties in the genome of Plasmodium falciparum

A method is presented to computationally identify core promoters in the Plasmodium falciparum genome using only DNA physicochemical properties.

Little is known about the structure and distinguishing features of core promoters in Plasmodium falciparum. In this work, we describe the first method to computationally identify core promoters in this AT-rich genome. This prediction algorithm uses solely DNA physicochemical properties as descriptors. Our results add to a growing body of evidence that a physicochemical code for eukaryotic genomes plays a crucial role in core promoter recognition.

Single-molecule FRET measures bends and kinks in DNA

We present advances in the use of single-molecule FRET measurements with flexibly linked dyes to derive full 3D structures of DNA constructs based on absolute distances. The resolution obtained by this single-molecule approach harbours the potential to study in detail also protein- or damage-induced DNA bending. If one is to generate a geometric structural model, distances between fixed positions are needed. These are usually not experimentally accessible because of unknown fluorophore-linker mobility effects that lead to a distribution of FRET efficiencies and distances. To solve this problem, we performed studies on DNA double-helices by systematically varying donor acceptor distances from 2 to 10 nm. Analysis of dye-dye quenching and fluorescence anisotropy measurements reveal slow positional and fast orientational fluorophore dynamics, that results in an isotropic average of the FRET efficiency. We use a nonlinear conversion function based on MD simulations that allows us to include this effect in the calculation of absolute FRET distances. To obtain unique structures, we performed a quantitative statistical analysis for the conformational search in full space based on triangulation, which uses the known helical nucleic acid features. Our higher accuracy allowed the detection of sequence-dependent DNA bending by 16 degrees . For DNA with bulged adenosines, we also quantified the kink angles introduced by the insertion of 1, 3 and 5 bases to be 32 degrees +/- 6 degrees , 56 degrees +/- 4 degrees and 73 +/- 2 degrees , respectively. Moreover, the rotation angles and shifts of the helices were calculated to describe the relative orientation of the two arms in detail.

Remeasuring the double helix

DNA is thought to behave as a stiff elastic rod with respect to the ubiquitous mechanical deformations inherent to its biology. To test this model at short DNA lengths, we measured the mean and variance of end-to-end length for a series of DNA double helices in solution, using small-angle x-ray scattering interference between gold nanocrystal labels. In the absence of applied tension, DNA is at least one order of magnitude softer than measured by single-molecule stretching experiments. Further, the data rule out the conventional elastic rod model. The variance in end-to-end length follows a quadratic dependence on the number of base pairs rather than the expected linear dependence, indicating that DNA stretching is cooperative over more than two turns of the DNA double helix. Our observations support the idea of long-range allosteric communication through DNA structure.

Intrinsic curvature associated with the coordinately regulated anthrax toxin gene promoters

The current model for virulence gene regulation in Bacillus anthracis involves several trans-acting factors, the most important of which appears to be the anthrax toxin activator encoded by the atxA gene. AtxA is a positive regulator of the toxin genes pagA, cya and lef, and of a number of other plasmid- and chromosome-encoded genes. The AtxA protein (56 kDa) possesses a predicted winged-helix DNA-binding domain and phosphotransferase system-regulated domains, but the mechanism for positive regulation of AtxA target genes is not known. Sequence similarities in the promoter regions of AtxA-regulated genes are not apparent, and recombinant AtxA binds DNA with a high affinity in a non-specific manner. We hypothesized that the toxin genes possess common structural features or cis-acting elements that are required for positive regulation. We employed deletion analyses to determine the minimal sequences required for atxA-mediated toxin gene expression. In silico modelling and in vitro experiments using double-stranded DNA corresponding to the toxin gene promoter regions indicated significant curvature associated with these regions. These findings suggest that the structural topology of the DNA plays an important role in the control of anthrax toxin gene expression.

Helical coherence of DNA in crystals and solution

The twist, rise, slide, shift, tilt and roll between adjoining base pairs in DNA depend on the identity of the bases. The resulting dependence of the double helix conformation on the nucleotide sequence is important for DNA recognition by proteins, packaging and maintenance of genetic material, and other interactions involving DNA. This dependence, however, is obscured by poorly understood variations in the stacking geometry of the same adjoining base pairs within different sequence contexts. In this article, we approach the problem of sequence-dependent DNA conformation by statistical analysis of X-ray and NMR structures of DNA oligomers. We evaluate the corresponding helical coherence length--a cumulative parameter quantifying sequence-dependent deviations from the ideal double helix geometry. We find, e.g. that the solution structure of synthetic oligomers is characterized by 100-200 A coherence length, which is similar to approximately 150 A coherence length of natural, salmon-sperm DNA. Packing of oligomers in crystals dramatically alters their helical coherence. The coherence length increases to 800-1200 A, consistent with its theoretically predicted role in interactions between DNA at close separations.

Molecular modeling of nucleic acid structure

This unit is the first in a series of four units covering the analysis of nucleic acid structure by molecular modeling. This unit provides an overview of computer simulation of nucleic acids. Topics include the static structure model, computational graphics and energy models, generation of an initial model, and characterization of the overall three-dimensional structure.

Intrinsic curvature of DNA influences LacR-mediated looping

Protein-mediated DNA looping is a common mechanism for regulating gene expression. Loops occur when a protein binds to two operators on the same DNA molecule. The probability of looping is controlled, in part, by the basepair sequence of inter-operator DNA, which influences its structural properties. One structural property is the intrinsic or stress-free curvature. In this article, we explore the influence of sequence-dependent intrinsic curvature by exercising a computational rod model for the inter-operator DNA as applied to looping of the LacR-DNA complex. Starting with known sequences for the inter-operator DNA, we first compute the intrinsic curvature of the helical axis as input to the rod model. The crystal structure of the LacR (with bound operators) then defines the requisite boundary conditions needed for the dynamic rod model that predicts the energetics and topology of the intervening DNA loop. A major contribution of this model is its ability to predict a broad range of published experimental data for highly bent (designed) sequences. The model successfully predicts the loop topologies known from fluorescence resonance energy transfer measurements, the linking number distribution known from cyclization assays with the LacR-DNA complex, the relative loop stability known from competition assays, and the relative loop size known from gel mobility assays. In addition, the computations reveal that highly curved sequences tend to lower the energetic cost of loop formation, widen the energy distribution among stable and meta-stable looped states, and substantially alter loop topology. The inclusion of sequence-dependent intrinsic curvature also leads to nonuniform twist and necessitates consideration of eight distinct binding topologies from the known crystal structure of the LacR-DNA complex.

A triple helix within a pseudoknot is a conserved and essential element of telomerase RNA

Telomerase copies a short template within its integral telomerase RNA onto eukaryotic chromosome ends, compensating for incomplete replication and degradation. Telomerase action extends the proliferative potential of cells, and thus it is implicated in cancer and aging. Nontemplate regions of telomerase RNA are also crucial for telomerase function. However, they are highly divergent in sequence among species, and their roles are largely unclear. Using in silico three-dimensional modeling, constrained by mutational analysis, we propose a three-dimensional model for a pseudoknot in telomerase RNA of the budding yeast Kluyveromyces lactis. Interestingly, this structure includes a U-A.U major-groove triple helix. We confirmed the triple-helix formation in vitro using oligoribonucleotides and showed that it is essential for telomerase function in vivo. While triplex-disrupting mutations abolished telomerase function, triple compensatory mutations that formed pH-dependent G-C.C(+) triples restored the pseudoknot structure in a pH-dependent manner and partly restored telomerase function in vivo. In addition, we identified a novel type of triple helix that is formed by G-C.U triples, which also partly restored the pseudoknot structure and function. We propose that this unusual structure, so far found only in telomerase RNA, provides an essential and conserved telomerase-specific function.

NMR structure of the full-length linear dimer of stem-loop-1 RNA in the HIV-1 dimer initiation site

The packaging signal of HIV-1 RNA contains a stem-loop structure, SL1, which serves as the dimerization initiation site for two identical copies of the genome and is important for packaging of the RNA genome into the budding virion and for overall infectivity. SL1 spontaneously dimerizes via a palindromic hexanucleotide sequence in its apical loop, forming a metastable kissing dimer form. Incubation with nucleocapsid protein causes this form to refold to a thermodynamically stable mature linear dimer. Here, we present an NMR structure of the latter form of the full-length SL1 sequence of the Lai HIV-1 isolate. The structure was refined using nuclear Overhauser effect and residual dipolar coupling data. The structure presents a symmetric homodimer of two RNA strands of 35 nucleotides each; it includes five stems separated by four internal loops. The central palindromic stem is surrounded by two symmetric adenine-rich 1-2 internal loops, A-bulges. All three adenines in each A-bulge are stacked inside the helix, consistent with the solution structures of shorter SL1 constructs determined previously. The outer 4-base pair stems and, proximal to them, purine-rich 1-3 internal loops, or G-bulges, are the least stable parts of the molecule. The G-bulges display high conformational variability in the refined ensemble of structures, despite the availability of many structural restraints for this region. Nevertheless, most conformations share a similar structural motif: a guanine and an adenine from opposite strands form a GA mismatch stacked on the top of the neighboring stem. The two remaining guanines are exposed, one in the minor groove and another in the major groove side of the helix, consistent with secondary structure probing data for SL1. These guanines may be recognized by the nucleocapsid protein, which binds tightly to the G-bulge in vitro.

Formation and positioning of nucleosomes: effect of sequence-dependent long-range correlated structural disorder

The understanding of the long-range correlations (LRC) observed in DNA sequences is still an open and very challenging problem. In this paper, we start reviewing recent results obtained when exploring the scaling properties of eucaryotic, eubacterial and archaeal genomic sequences using the space-scale decomposition provided by the wavelet transform (WT). These results suggest that the existence of LRC up to distances approximately 20-30 kbp is the signature of the nucleosomal structure and dynamics of the chromatin fiber. Actually the LRC are mainly observed in the DNA bending profiles obtained when using some structural coding of the DNA sequences that accounts for the fluctuations of the local double-helix curvature within the nucleosome complex. Because of the approximate planarity of nucleosomal DNA loops, we then study the influence of the LRC structural disorder on the thermodynamical properties of 2D elastic chains submitted locally to mechanical/topological constraint as loops. The equilibrium properties of the one-loop system are derived numerically and analytically in the quite realistic weak-disorder limit. The LRC are shown to favor the spontaneous formation of small loops, the larger the LRC, the smaller the size of the loop. We further investigate the dynamical behavior of such a loop using the mean first passage time (MFPT) formalism. We show that the typical short-time loop dynamics is superdiffusive in the presence of LRC. For displacements larger than the loop size, we use large-deviation theory to derive a LRC-dependent anomalous-diffusion rule that accounts for the lack of disorder self-averaging. Potential biological implications on DNA loops involved in nucleosome positioning and dynamics in eucaryotic chromatin are discussed.

X-ray diffraction "fingerprinting" of DNA structure in solution for quantitative evaluation of molecular dynamics simulation

Solution state x-ray diffraction fingerprinting is demonstrated as a method for experimentally assessing the accuracy of molecular dynamics (MD) simulations. Fourier transforms of coordinate data from MD simulations are used to produce reciprocal space "fingerprints" of atomic pair distance correlations that are characteristic of the ensemble and are the direct numerical analogues of experimental solution x-ray diffraction (SXD). SXD experiments and MD simulations were carried out to test the ability of experiment and simulation to resolve sequence-dependent modifications in helix conformation for B-form DNA. SXD experiments demonstrated that solution-state poly(AT) and poly(A)-poly(T) duplex DNA sequences exist in ensembles close to canonical B-form and B'-form structures, respectively. In contrast, MD simulations analyzed in terms of SXD fingerprints are shown to deviate from experiment, most significantly for poly(A)-poly(T) duplex DNA. Compared with experiment, MD simulation shortcomings were found to include both mismatches in simulated conformer structures and number population within the ensembles. This work demonstrates an experimental approach for quantitatively evaluating MD simulations and other coordinate models to simulate biopolymer structure in solution and suggests opportunities to use solution diffraction data as experimental benchmarks for developing supramolecular force fields optimized for a range of in situ applications.

Role of stacking interactions in the binding sequence preferences of DNA bis-intercalators: insight from thermodynamic integration free energy simulations

The major structural determinant of the preference to bind to CpG binding sites on DNA exhibited by the natural quinoxaline bis-intercalators echinomycin and triostin A, or the quinoline echinomycin derivative, 2QN, is the 2-amino group of guanine (G). However, relocation of this group by means of introduction into the DNA molecule of the 2-aminoadenine (=2,6-diaminopurine, D) base in place of adenine (A) has been shown to lead to a drastic redistribution of binding sites, together with ultratight binding of 2QN to the sequence DTDT. Also, the demethylated triostin analogs, TANDEM and CysMeTANDEM, which bind with high affinity to TpA steps in natural DNA, bind much less tightly to CpI steps, despite the fact that both adenosine and the hypoxanthine-containing nucleoside, inosine (I), provide the same hydrogen bonding possibilities in the minor groove. To study both the increased binding affinity of 2QN for DTDT relative to GCGC sites and the remarkable loss of binding energy between CysMeTANDEM and ICIC compared with ATAT, a series of thermodynamic integration free energy simulations involving conversions between DNA base pairs have been performed. Our results demonstrate that the electrostatic component of the stacking interactions between the heteroaromatic rings of these compounds and the bases that make up the intercalation sites plays a very important role in the modulation of their binding affinities.

Fine mapping of inherent flexibility variation along DNA molecules: validation by atomic force microscopy (AFM) in buffer

Curvature and flexibility are structural properties of central importance to genome function. However, due to the difficulties in finding suitable experimental conditions, methods for studying one without the interference of the other have proven to be difficult. We propose a new approach that provides a measure of inherent flexibility of DNA by taking advantage of two powerful techniques, X-ray crystallography and nuclear magnetic resonance. Both techniques are able to detect local curvature on DNA fragments but, while the first analyzes DNA in the solid state, the second works on DNA in solution. Comparison of the two data sets allowed us to calculate the relative contribution to flexibility of the three rotations and three translations, which relate successive base pair planes for the ten different dinucleotide steps. These values were then used to compute the variation of flexibility along a given nucleotide sequence. This allowed us to validate the method experimentally through comparisons with maps of local fluctuations in DNA molecule trajectory constructed from atomic force microscopy imaging in solution. We conclude that the six dinucleotide-step parameters defined here provide a powerful tool for the exploration of DNA structure and, consequently will make an important contribution to our understanding of DNA-sequence-dependent biological processes.

Solution structure of a DNA duplex containing an alpha-anomeric adenosine: insights into substrate recognition by endonuclease IV

The cytotoxic alpha anomer of adenosine, generated in situ by radicals, must be recognized and repaired to maintain genomic stability. Endonuclease IV (Endo IV), a member of the base excision repair (BER) enzyme family, in addition to acting on abasic sites, has the auxiliary function of removing this mutagenic nucleotide in Escherichia coli. We have employed enzymatic, thermodynamic, and structural studies on DNA duplexes containing a central alpha-anomeric adenosine residue to characterize the role of DNA structure on recognition and catalysis by Endo IV. The enzyme recognizes and cleaves our alphaA-containing DNA duplexes at the site of the modification. The NMR solution structure of the DNA decamer duplex establishes that the single alpha-anomeric adenosine residue is intrahelical and stacks in a reverse Watson-Crick fashion consistent with the slight decrease in thermostability. However, the presence of this lesion confers significant changes to the global duplex conformation, resulting from a kink of the helical axis into the major groove and an opening of the minor groove emanating from the alpha-anomeric site. Interestingly, the conformation of the flanking base-paired segments is not greatly altered from a B-type conformation. The global structural changes caused by this lesion place the DNA along the conformational path leading to the DNA structure observed in the complex. Thus, it appears that the alpha-anomeric lesion facilitates recognition by Endo IV.

NMR structure of an alpha-L-LNA:RNA hybrid: structural implications for RNase H recognition

Alpha-L-LNA (alpha-L-ribo configured locked nucleic acid) is a nucleotide analogue that raises the thermostability of nucleic acid duplexes by up to approximately 4 degrees C per inclusion. We have determined the NMR structure of a nonamer alpha-L-LNA:RNA hybrid with three alpha-L-LNA modifications. The geometry of this hybrid is intermediate between A- and B-type, all nucleobases partake in Watson-Crick base pairing and base stacking, and the global structure is very similar to that of the corresponding unmodified hybrid. The sugar-phosphate backbone is rearranged in the vicinity of the modified nucleotides. As a consequence, the phosphate groups following the modified nucleotides are rotated into the minor groove. It is interesting that the alpha-L-LNA:RNA hybrid, which has an elevation in melting temperature of 17 degrees C relative to the corresponding DNA:RNA hybrid, retains the global structure of this hybrid. To our knowledge, this is the first example of such a substantial increase in melting temperature of a nucleic acid analogue that does not act as an N-type (RNA) mimic. alpha-L-LNA:RNA hybrids are recognised by RNase H with subsequent cleavage of the RNA strand, albeit with slow rates. We attempt to rationalise this impaired enzyme activity from the rearrangement of the sugar-phosphate backbone of the alpha-L-LNA:RNA hybrid.

Non-random features of loop-size chromatin fragmentation

Upon isolation of DNA from normal eukaryotic cells by standard methods involving extensive proteolytic treatment, a rather homogeneous population of loop-size, double-stranded DNA fragments is regularly obtained. These DNA molecules can be efficiently end-labeled by the DNA polymerase I Klenow fragment, as well as by a 3'- to -5'-exonuclease-free Klenow enzyme, but not by terminal transferase (TdT) unless the ends have been filled up by Klenow, suggesting that dominantly 5' protruding termini are generated upon fragmentation. The filled-up termini were used for cloning the distal parts of the approximately 50 kb fragments. BLAST analysis of the sequence of several clones allowed us to determine the sequence of the non-cloned side of the breakpoints. Comparison of 25, 600 bp-long breakpoint sequences demonstrated prevalence of repetitive elements. Consensus motives characteristic of the breakpoint sequences have been identified. Several sequences exhibit peculiar computed conformational characteristics, with sharp transition or center of symmetry located exactly at the breakpoint. Our data collectively suggest that chromatin fragmentation involves nucleolytic cleavages at fragile/hypersensitive sites delimiting loop-size fragments in a non-random manner. Interestingly, the sequence characteristics of the breakpoints are reminiscent of certain breakpoint cluster regions frequently subject to gene rearrangements.

DNA analysis servers: plot.it, bend.it, model.it and IS

The WWW servers at http://www.icgeb.trieste.it/dna/ are dedicated to the analysis of user-submitted DNA sequences; plot.it creates parametric plots of 45 physicochemical, as well as statistical, parameters; bend.it calculates DNA curvature according to various methods. Both programs provide 1D as well as 2D plots that allow localisation of peculiar segments within the query. The server model.it creates 3D models of canonical or bent DNA starting from sequence data and presents the results in the form of a standard PDB file, directly viewable on the user's PC using any molecule manipulation program. The recently established introns server allows statistical evaluation of introns in various taxonomic groups and the comparison of taxonomic groups in terms of length, base composition, intron type etc. The options include the analysis of splice sites and a probability test for exon-shuffling.

NMR solution structure of dsDNA containing a bicyclic D-arabino-configured nucleotide fixed in an O4'-endo sugar conformation

[3.2.0]bcANA is a D-arabino-configured bicyclic nucleotide with a 2'-O,3'-C-methylene bridge. We here present the high-resolution NMR structure of a [3.2.0]bcANA modified dsDNA nonamer with one modified nucleotide incorporated. NOE restraints were obtained by analysis of NOESY cross peak intensities using a full relaxation matrix approach, and subsequently these restraints were incorporated into a simulated annealing scheme for the structure determination. In addition, the furanose ring puckers of the deoxyribose moieties were determined by analysis of COSY cross peaks. The modified duplex adopts a B-like geometry with Watson-Crick base pairing in all base pairs and all glycosidic angles in the anti range. The stacking arrangement of the nucleobases appears to be unperturbed relative to the normal B-like arrangement. The 2'-O,3'-C-methylene bridge of the modified nucleotide is located at the brim of the major groove where it fits well into the B-type duplex framework. The sugar pucker of the [3.2.0]bcANA nucleotide is O4'-endo and this sugar conformation causes a change in the delta backbone angle relative to the C2'-endo deoxyribose sugar pucker. This change is absorbed locally by slight changes in the epsilon and zeta angles of the modified nucleotide. Overall, the [3.2.0]bcANA modifications fits very well into a B-like duplex framework and only small and local perturbations are observed relative to the unmodified dsDNA of identical base sequence.

An assessment of three dinucleotide parameters to predict DNA curvature by quantitative comparison with experimental data

Curved DNA fragments are often found near functionally important sites such as promoters and origins of replication, and hence sequence-dependent DNA curvature prediction is of great utility in genomics and bioinformatics. In light of this, an assessment of three different dinucleotide step parameters (based on gel retardation as well as crystal structure data) is carried out. These parameters (BMHT, LB and CS) are evaluated quantitatively for their ability to predict correctly the experimental results of a large set of nucleic acid sequences containing A-tracts as well as GC-rich motifs. This set contained around 40 synthetic as well as natural sequences whose solution properties have been well characterized experimentally. All three models could account reasonably well for curvature in the various DNA sequences. The CS model, where dinucleotide parameters are calculated from crystal structure data, consistently shows slightly better correlation with experimental data. Our simple analysis also indicates that presently available trinucleotide parameters fail to predict curvature in some of the well-characterized sequences. The study shows that the dinucleotide parameters with some further refinement can be used to predict sequence-dependent curvature correctly in genomic sequences.

Influence of the sequence on elastic properties of long DNA chains

We revisit the results of single-molecule DNA stretching experiments using a rodlike chain (RLC) model that explicitly includes some intrinsic structural disorder induced by the sequence. The investigation of artificial and real genomic sequences shows that the wormlike chain model reproduces quite well the data but with an effective bend stiffness A(eff), which underestimates the true elastic bend stiffness A, independently of the elastic twist stiffness C. Mainly dominated by the amplitude of the structural disorder, this correction seems rather insensitive to the presence of long-range correlations. This RLC model is shown to remarkably fit the experimental data for lambda-DNA when considering A approximately 70+/-10 nm (>A(eff) approximately 50 nm), in good agreement with previous experimental estimates of the "dynamic" persistent length. From the analysis of large human contigs, we speculate about the possible dependence of A(eff) and/or A upon the (G+C) content of the considered sequence.

Accuracy of AFM measurements of the contour length of DNA fragments adsorbed on mica in air and in aqueous buffer

The measurement by atomic force microscope of the contour length of DNA fragments adsorbed on mica has been made as accurate as possible by revisiting the different steps of image acquisition and processing. In air, the DNA helical rise was estimated at 2.97 +/- 0.15 A per base pair (bp) (mean +/- standard deviation) by imaging a 648-bp DNA fragment and 2.95 +/- 0.14 A per bp for a 115-bp fragment. This confirms earlier observations suggesting that drying DNA fragments on mica in the presence of nickel induces limited conformational changes. At this point the exact nature of these conformational changes remains unknown. Simple hypotheses are the transconformation of stretches of the DNA molecules to the A-form of the double helix or alteration of the helix structure at the points of contact between DNA and mica. By contrast, in aqueous buffer, the measured helical rise was 3.14 +/- 0.15 A per bp for the 648-bp fragment and 3.17 +/- 0.13 A per bp for the 1115-bp fragment. Thus, measured helical rises do not depend on the fragment length and are significantly shorter than the 3.38 A per bp measured by crystallography, but close to the 3.18 A per bp found in NMR studies. These findings are discussed with respect to discrepancies in earlier results published in the literature.

Low temperature solution structures and base pair stacking of double helical d(CGTACG)(2)

Solution structures and base pair stacking of a self- complementary DNA hexamer d(CGTACG)(2) have been studied at 5, 10 and 15 degrees C, respectively. The stacking interactions among the center base pair steps of the DNA duplex are found to improve when the terminal base pairs became less stable due to end fraying. A new structural quantity, the stacking sum (Sigma(s)), is introduced to indicate small changes in the stacking overlaps between base pairs. The improvements in the stacking overlaps to maintain the double helical conformation are probably the cause for the observed temperature dependent structural changes in double helical DNA molecule. A detailed analysis of the helical parameters, backbone torsion angles, base orientations and sugar conformations of these structures has been performed.

Molecular dynamics simulation of nucleic acids: successes, limitations, and promise

In the last five years we have witnessed a significant increase in the number publications describing accurate and reliable all-atom molecular dynamics simulations of nucleic acids. This increase has been facilitated by the development of fast and efficient methods for treating the long-range electrostatic interactions, the availability of faster parallel computers, and the development of well-validated empirical molecular mechanical force fields. With these technologies, it has been demonstrated that simulation is not only capable of consistently reproducing experimental observations of sequence specific fine structure of DNA, but also can give detailed insight into prevalent problems in nucleic acid structure, ion association and specific hydration of nucleic acids, polyadenine tract bending, and the subtle environmental dependence of the A-DNA-B-DNA duplex equilibrium. Despite the advances, there are still issues with the methods that need to be resolved through rigorous controlled testing. In general, these relate to deficiencies of the underlying molecular mechanical potentials or applied methods (such as the imposition of true periodicity in Ewald simulations and the need for energy conservation), and significant limits in effective conformational sampling. In this perspective, we provide an overview of our experiences, provide some cautionary notes, and provide recommendations for further study in molecular dynamics simulation of nucleic acids.

1H NMR studies of a 17-mer DNA duplex

Transcription factor 1 (TF1), encoded by the Bacillus subtilis bacteriophage SPO1, is a DNA-binding protein of the HU family. In preparation for a determination of the structure of the DNA-TF1 complex, we have studied the conformation of one core 17-mer duplex d(5'-CACTACTCTTTGTAGTG-3')-d(5'-CACTACAAAGAGTAGTG-3'). NOESY, DQF-COSY and TOCSY spectroscopy provide resonance assignments of non-exchangeable and exchangeable protons, internucleotide and interstrand proton-proton distances, and dihedral angle constraints. Restrained molecular dynamics calculations yield a family of NMR solution structures for which the RMSD is 0.7 A (all atoms). The helical twist is 34.9 degrees for the central 15 bp. Bends toward the major groove are located between the second and fourth base pairs from each end. The G12 x C23 base pair, which is bounded on each side by consecutive A x T pairs, causes a local disturbance to the DNA helix that makes the conformations of the two end segments unsymmetrical. The pyrimidine rings at T9, T10 and T11 experience more extensive rotational movement than the rest of the structure.

Solution structure of the dodecamer d-(CATGGGCC-CATG)2 is B-DNA. Experimental and molecular dynamics study

The DNA duplex d-(CATGGGCCCATG)2 has been studied in solution by FTIR, NMR and CD. The experimental approaches have been complemented by series of large-scale unrestrained molecular dynamics simulation with explicit inclusion of solvent and counterions. Typical proton-proton distances extracted from the NMR spectra and the CD spectra are completely in agreement with slightly modified B-DNA. By molecular dynamics simulation, starting from A-type sugar pucker, a spontaneous repuckering to B-type sugar pucker was observed. Both experimental and theoretical approaches suggest for the dodecamer d-(CATGGGCCCATG)2 under solution conditions puckering of all 2'-deoxyribose residues in the south conformation (mostly C2'-endo) and can exclude significant population of sugars in the north conformation (C3'-endo). NMR, FTIR and CD data are in agreement with a B-form of the dodecamer in solution. Furthermore, the duplex shows a cooperative B-A transition in solution induced by addition of trifluorethanol. This contrasts a recently published crystal structure of the same oligonucleotide found as an intermediate between B- and A-DNA where 23 out of 24 sugar residues were reported to adopt the north (N-type) conformation (C3'-endo) like in A-DNA (Ng, H. L., Kopka, M. L. and Dickerson, R. E., Proc. Natl. Acad. Sci. U S A 97, 2035-2039 (2000)). The simulated structures resemble standard B-DNA. They nevertheless show a moderate shift towards A-type stacking similar to that seen in the crystal, despite the striking difference in sugar puckers between the MD and X-ray structures. This is in agreement with preceding MD reports noticing special stacking features of G-tracts exhibiting a tendency towards the A-type stacking supported by the CD spectra also reflecting the G-tract stacking. MD simulations reveal several noticeable local conformational variations, such as redistribution of helical twist and base pair roll between the central GpC steps and the adjacent G-tract segments, as well as a substantial helical twist variability in the CpA(TpG) steps combined with a large positive base pair roll. These local variations are rather different from those seen in the crystal.

Molecular dynamics simulation of nucleic acids

We review molecular dynamics simulations of nucleic acids, including those completed from 1995 to 2000, with a focus on the applications and results rather than the methods. After the introduction, which discusses recent advances in the simulation of nucleic acids in solution, we describe force fields for nucleic acids and then provide a detailed summary of the published literature. We emphasize simulations of small nucleic acids ( approximately 6 to 24 mer) in explicit solvent with counterions, using reliable force fields and modern simulation protocols that properly represent the long-range electrostatic interactions. We also provide some limited discussion of simulation in the absence of explicit solvent. Absent from this discussion are results from simulations of protein-nucleic acid complexes and modified DNA analogs. Highlights from the molecular dynamics simulation are the spontaneous observation of A B transitions in duplex DNA in response to the environment, specific ion binding and hydration, and reliable representation of protein-nucleic acid interactions. We close by examining major issues and the future promise for these methods.

Induced fit DNA recognition by a minor groove binding analogue of Hoechst 33258: fluctuations in DNA A tract structure investigated by NMR and molecular dynamics simulations

NMR analysis and molecular dynamics simulations of d(GGTAATTACC)(2) and its complex with a tetrahydropyrimidinium analogue of Hoechst 33258 suggest that DNA minor groove recognition in solution involves a combination of conformational selection and induced fit, rather than binding to a preorganised site. Analysis of structural fluctuations in the bound and unbound states suggests that the degree of induced fit observed is primarily a consequence of optimising van der Waals contacts with the walls of the minor groove resulting in groove narrowing through: (i) changes in base step parameters, including increased helical twist and propeller twist; (ii) changes to the sugar-phosphate backbone conformation to engulf the bound ligand; (iii) suppression of bending modes at the TpA steps. In contrast, the geometrical arrangement of hydrogen bond acceptors on the groove floor appears to be relatively insensitive to DNA conformation (helical twist and propeller twist). We suggest that effective recognition of DNA sequences (in this case an A tract structure) appears to depend to a significant extent on the sequence being flexible enough to be able to adopt the geometrically optimal conformation compatible with the various binding interactions, rather than involving 'lock and key' recognition.

A refined prediction method for gel retardation of DNA oligonucleotides from dinucleotide step parameters: reconciliation of DNA bending models with crystal structure data

The development and assessment of a prediction method for gel retardation and sequence dependent curvature of DNA based on dinulcleotide step parameters are described. The method is formulated using the Babcock-Olson equations for base pair step geometry (1) and employs Monte Carlo simulated annealing for parameter optimization against experimental data. The refined base pair step parameters define a stuctural construct which, when the width of observed parameter distributions is taken into account, is consistent with the results of DNA oligonucleotide crystal structures. The predictive power of the method is demonstrated and tested via comparisons with DNA bending data on sets of sequences not included in the training set, including A-tracts with and without periodic helix phasing, phased A4T4 and T4A4 motifs, a sequence with a phased GGGCCC motif, some "unconventional" helix phasing sequences, and three short fragments of kinetoplast DNA from Crithidia fasiculata that exhibit significantly different behavior on non-denaturing polyacrylamide gels. The nature of the structural construct produced by the methodology is discussed with respect to static and dynamic models of structure and representations of bending and bendability. An independent theoretical account of sequence dependent chemical footprinting results is provided. Detailed analysis of sequences with A-tract induced axis bending forms the basis for a critical discussion of the applicability of wedge models,junction models and non A-tract, general sequence models for understanding the origin of DNA curvature at the molecular level.

Intrinsically bent DNA in the promoter regions of the yeast GAAL1-10 and GAL80 genes

Circular permutation analysis has detected fairly strong sites of intrinsic DNA bending on the promoter regions of the yeast GAL1-10 and GAL80 genes. These bends lie in functionally suggestive locations. On the promoter of the GAL1-10 structural genes, strong bends bracket nucleosome B, which lies between the UAS(G) and the GAL1 TATA. These intrinsic bends could help position nucleosome B. Nucleosome B plus two other promoter nucleosomes protect the TATA and start site elements in the inactive state of expression but are completely disrupted (removed) when GAL1-10 expression is induced. The strongest intrinsic bend ( approximately 70 degrees ) lies at the downstream edge of nucleosome B; this places it approximately 30 base pairs upstream of the GAL1 TATA, a position that could allow it to be involved in GAL1 activation in several ways, including the recruitment of a yeast HMG protein that is required for the normally robust level of GAL1 expression in the induced state (Paull, T., Carey, M., and Johnson, R. (1996) Genes Dev. 10, 2769-2781). On the regulatory gene GAL80, the single bend lies in the non-nucleosomal hypersensitive region, between a GAL80-specific far upstream promoter element and the more gene-proximal promoter elements. GAL80 promoter region nucleosomes contain no intrinsically bent DNA.

Targeting of high mobility group-14/-17 proteins in chromatin is independent of DNA sequence

Chromosomal proteins high mobility group (HMG)-14 and HMG-17 are nucleosomal-binding proteins that unfold the chromatin fiber and enhance transcription from chromatin templates. Their intracellular organization is dynamic and related to both cell cycle and transcription. Here we examine possible mechanisms for targeting HMG-14/-17 to specific regions in chromatin. Chromatin immunoprecipitation assays indicate that HMG-17 protein is not preferentially associated with chromatin regions containing transcriptionally active genes, or any type of specific DNA. We used a modification of the random amplified polymorphic DNA method to analyze DNA in various HMG-14/-17.nucleosome complexes. We found that although HMG-14 or HMG-17 proteins preferentially associate with core particles in which the DNA has a low frequency of CG dinucleotides, the genome does not contain consensus sequences that serve as specific targeting sites for the binding of either HMG-14 or HMG-17 proteins to nucleosomes. We used size exclusion and ion exchange chromatography to demonstrate that nuclei contain a large portion of HMG-17 associated with other proteins in a multiprotein complex. We suggest that these complexes regulate the dynamic organization of HMG-14/-17 in the nucleus and serve to target the proteins to specific sites in chromatin.

A statistical mechanical model for predicting B-DNA curvature and flexibility

A statistical mechanical model taking into account the symmetric twisting, tilting, sliding fluctuations and asymmetric rolling fluctuations has been proposed to predict the macroscopic curvature and flexibility of B-DNA. Based on the statistical data of structural parameters of double helix in nucleic acid database and the related theoretical analysis, the equilibrium angular parameters (Omega, rho and tau) describing the orientation of successive base-pair planes, the translation parameters (D(y)) along the long axis of neighboring base-pair step and the corresponding force constants are arranged for ten dimers appropriately. Under the assumption of independent angular parameters, independent base-pair steps and a simple energy function, we can calculate the macroscopic curvature and the flexibility of DNA sequences through the transformation matrix <R> and the Boltzmann ensemble average. The predictions on curvature and flexibility of DNA have been compared with the corresponding experimental data. The agreement is remarkably good. It is demonstrated that the lowering of the temperature does increase the DNA curvature.

NMR structure of a DNA duplex containing nucleoside analog 1-(2'-deoxy-beta-D-ribofuranosyl)-3-nitropyrrole and the structure of the unmodified control

The three-dimensional structures of two DNA duplexes d(CATGAGTAC). d(GTACXCATG) (1) and d(CATGAGTAC).d(GTACTCATG) (2), where X represents 1-(2'-deoxy-beta-D-ribofuranosyl)-3-nitropyrrole, were solved using high resolution nuclear magnetic resonance spectroscopy and restrained molecular dynamics. Good convergence was observed between final structures derived from A- and B-form starting geometries for both 1 and 2. Structures of 1 and 2 are right-handed duplexes within the B-form conformational regime. Furthermore, the structures of 1 and 2 are highly similar, with differences in the structures localized to the vicinity of residue 14 (X versus T). The pyrrole group of 1 is in the syn conformation and it is displaced towards the major groove. Furthermore, unlike T14 in 2, the base of X14 has reduced pi-pi stacking interactions with C13 and C15 and the nitro group of X14 is pointing out of the major groove. The structures presented here establish the basis of the thermal data of DNA duplexes containing X and should be informative during the design of improved wild card nucleobase analogs.

Yeast chromatin structure and regulation of GAL gene expression

Yeast genomic DNA is covered by nucleosome cores spaced by short, discrete length linkers. The short linkers, reinforced by novel histone properties, create a number of unique and dynamic nucleosome structural features in vivo: permanent unpeeling of DNA from the ends of the core, an inability to bind even full 147 bp core DNA lengths, and facility to undergo a conformational transition that resembles the changes found in active chromatin. These features probably explain how yeast can maintain most of its genome in a transcribable state and avoid large-scale packaging away of inactive genes. The GAL genes provide a closely regulated system in which to study gene-specific chromatin structure. GAL structural genes are inactive without galactose but are highly transcribed in its presence; the expression patterns of the regulatory genes can account for many of the features of GAL structural gene control. In the inactive state, GAL genes demonstrate a characteristic promoter chromosomal organization; the major upstream activation sequence (UASG) elements lie in open, hypersensitive regions, whereas the TATA and transcription start sites are in nucleosomes. This organization helps implement gene regulation in this state and may benefit the organism. Induction of GAL expression triggers Gal4p-dependent upstream nucleosome disruption. Disruption is transient and can readily be reversed by a Gal80p-dependent nucleosome deposition process. Both are sensitive to the metabolic state of the cell. Induction triggers different kinds of nucleosome changes on the coding sequences, perhaps reflecting the differing roles of nucleosomes on coding versus promoter regions. GAL gene activation is a complex process involving multiple Gal4p activities, numerous positive and negative cofactors, and the histone tails. DNA bending and chromosomal architecture of the promoter regions may also play a role in GAL regulation. Regulator-mediated competition between nucleosomes and the TATA binding protein complex for the TATA region is probably a central aspect of GAL regulation and a focal point for the numerous factors and processes that contribute to it.

Bending and adaptability to proteins of the cAMP DNA-responsive element: molecular dynamics contrasted with NMR

DNA bending is assumed to play a crucial role during recognition of the cAMP-responsive element (CRE) by transcription factors. However, diverging results have been obtained for the bending direction of the unbound double helix. The refined NMR structures present a bend directed toward the minor groove, while biochemical methods conclude that there is a bend toward the major groove. The present 10-ns molecular dynamics (MD) simulation of d(GAGATGACGTCATCTC)(2), which contains the octamer CRE in its center, was carried out with AMBER in explicit water and counterions. It shows that CRE is a flexible segment, although it is bent slightly toward the major groove (7 degrees -8 degrees ) on the average. The MD structure agrees with both the biochemical results and unrefined NMR data. The divergence with the NMR refined structures suggests an improper electrostatic parameterization in the refinement software. The malleability of the central CpG is certainly the major contribution to the curving of the whole CRE segment in both the unbound and bound states. Comparison with the crystal structure of CRE bound to GCN4 shows that the deformation induced by the protein is concentrated mainly on the CpG step, rendering the bound structure of CRE closer to the structure of the 12-0 tetradecanoylphorbol-beta-acetate-responsive element.

Aqueous trifluorethanol solutions simulate the environment of DNA in the crystalline state

We took 28 fragments of DNA whose crystal structures were known and used CD spectroscopy to search for conditions stabilising the crystal structures in solution. All 28 fragments switched into their crystal structures in 60-80% aqueous trifluorethanol (TFE) to indicate that the crystals affected the conformation of DNA like the concentrated TFE. The fragments crystallising in the B-form also underwent cooperative TFE-induced changes that took place within the wide family of B-form structures, suggesting that the aqueous and crystal B-forms differed as well. Spermine and magnesium or calcium cations, which were contained in the crystallisation buffers, promoted or suppressed the TFE-induced changes of several fragments to indicate that the crystallisation agents can decide which of the possible structures is adopted by the DNA fragment in the crystal.

A modified version of the Cornell et al. force field with improved sugar pucker phases and helical repeat

We have examined some subtle parameter modifications to the Cornell et al. force field, which has proven quite successful in reproducing nucleic acid properties, but whose C2'-endo sugar pucker phase and helical repeat for B DNA appear to be somewhat underestimated. Encouragingly, the addition of a single V2 term involving the atoms C(sp3)-O-(sp3)-C(sp3)-N(sp2), which can be nicely rationalized because of the anomeric effect (lone pairs on oxygen are preferentially oriented relative to the electron withdrawing N), brings the sugar pucker phase of C2'-endo sugars to near perfect agreement with ab initio calculations (W near 162 degrees). Secondly, the use of high level ab initio calculations on entire nucleosides (in contrast to smaller model systems necessitated in 1994-95 by computer limitations) lets one improve the chi torsional potential for nucleic acids. Finally, the O(sp3)-C(sp3)- C(sp3)-O(sp3) V2 torsional potential has been empirically adjusted to reproduce the ab initio calculated relative energy of C2'-endo and C3'-endo nucleosides. These modifications are tested in molecular dynamics simulations of mononucleosides (to assess sugar pucker percentages) and double helices of DNA and RNA (to assess helical and sequence specific structural properties). In both areas, the modified force field leads to improved agreement with experimental data.

Unconventional helical phasing of repetitive DNA motifs reveals their relative bending contributions

A novel, multiple DNA phasing analysis is described in which three sequence motifs associated with bent DNA are clustered together in oligomers of identical base composition, but with different phasing relationships of these motifs to each other. Synthetic oligonucleotides containing different combinations of AAAAA(A), GGGCCC and GAGAG sequence motifs were ligated and analyzed by gel mobility and cyclization experiments to determine their global curvature. These assays were used to obtain relative bending contributions of the analyzed sequence motifs. The experimental results also provide a rigorous test of predictive models for DNA bending. We report, using molecular modeling, that none of the most widely used dinucleotide (nearest neighbor) models can accurately describe the conformational properties of these DNA sequences and that more complex models, at least at the trinucleotide level, are required.