A Refined Prediction Method for Gel Retardation of DNA Oligonucleotides from Dinucleotide Step Parameters: Reconciliation of DNA Bending Models with Crystal Structure Data

Abstract

The development and assessment of a prediction method for gel retardation and sequence dependent curvature of DNA based on dinucleotide 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 structural 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 A₄T₄ and T₄A₄ 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.     

Sequence‐dependent modelling of local DNA bending phenomena: curvature prediction and vibrational analysis

Bending is a local conformational micropolymorphism of DNA in which the original B-DNA structure is only distorted but not extensively modified. Bending can be predicted by simple static geometry models as well as by a recently developed elastic model that incorporate sequence dependent anisotropic bendability (SDAB). The SDAB model qualitatively explains phenomena including affinity of protein binding, kinking, as well as sequence-dependent vibrational properties of DNA. The vibrational properties of DNA segments can be studied by finite element analysis of a model subjected to an initial bending moment. The frequency spectrum is obtained by applying Fourier analysis to the displacement values in the time domain. This analysis shows that the spectrum of the bending vibrations quite sensitively depends on the sequence, for example the spectrum of a curved sequence is characteristically different from the spectrum of straight sequence motifs of identical basepair composition. Curvature distributions are genome-specific, and pronounced differences are found between protein-coding and regulatory regions, respectively, that is, sites of extreme curvature and/or bendability are less frequent in protein-coding regions. A WWW server is set up for the prediction of curvature and generation of 3D models from DNA sequences (http://www.icgeb.trieste.it/dna).This revised version was published online in October 2005 with corrections to the Cover Date.     

Structural Mechanics of DNA Wrapping in the Nucleosome

Experimental X-ray crystal structures and a database of calculated structural parameters of DNA octamers were used in combination to analyse the mechanics of DNA bending in the nucleosome core complex. The 1kx5 X-ray crystal structure of the nucleosome core complex was used to determine the relationship between local structure at the base-step level and the global superhelical conformation observed for nucleosome-bound DNA. The superhelix is characterised by a large curvature (597°) in one plane and very little curvature (10°) in the orthogonal plane. Analysis of the curvature at the level of 10-step segments shows that there is a uniform curvature of 30° per helical turn throughout most of the structure but that there are two sharper kinks of 50° at ± 2 helical turns from the central dyad base pair. The curvature is due almost entirely to the base-step parameter roll. There are large periodic variations in roll, which are in phase with the helical twist and account for 500° of the total curvature. Although variations in the other base-step parameters perturb the local path of the DNA, they make minimal contributions to the total curvature. This implies that DNA bending in the nucleosome is achieved using the roll-slide-twist degree of freedom previously identified as the major degree of freedom in naked DNA oligomers. The energetics of bending into a nucleosome-bound conformation were therefore analysed using a database of structural parameters that we have previously developed for naked DNA oligomers. The minimum energy roll, the roll flexibility force constant and the maximum and minimum accessible roll values were obtained for each base step in the relevant octanucleotide context to account for the effects of conformational coupling that vary with sequence context. The distribution of base-step roll values and corresponding strain energy required to bend DNA into the nucleosome-bound conformation defined by the 1kx5 structure were obtained by applying a constant bending moment. When a single bending moment was applied to the entire sequence, the local details of the calculated structure did not match the experiment. However, when local 10-step bending moments were applied separately, the calculated structure showed excellent agreement with experiment. This implies that the protein applies variable bending forces along the DNA to maintain the superhelical path required for nucleosome wrapping. In particular, the 50° kinks are constraints imposed by the protein rather than a feature of the 1kx5 DNA sequence. The kinks coincide with a relatively flexible region of the sequence, and this is probably a prerequisite for high-affinity nucleosome binding, but the bending strain energy is significantly higher at these points than for the rest of the sequence. In the most rigid regions of the sequence, a higher strain energy is also required to achieve the standard 30° curvature per helical turn. We conclude that matching of the DNA sequence to the local roll periodicity required to achieve bending, together with the increased flexibility required at the kinks, determines the sequence selectivity of DNA wrapping in the nucleosome.     

Structural Basis of Stable Bending in DNA Containing An Tracts. Different Types of Bending

Abstract

Structural determinants of DNA bending of different types have been studied by theoretical conformational analysis of duplexes. Their terminal parts were fixed either in an ordinary low-energy B-like conformation or in “anomalous” conformations with a narrowed minor groove typical of A_n tracts. The anomalous conformations had different negative tilt angles (up to about zero), different propeller twists and minor groove widths. Calculations have been performed for DNA fragments A_nT_m, T_nA_m, A_nGCT_m, A_nCGT_m, T_mGCA_n, T_mCGA_n which are the models of the junction of two anomalous structures on A_n and T_m tracts. On the AT step of the A_nT_m fragment the minor groove can be easily narrowed so that a whole unbent fragment of anomalous structure is formed on A_n T_m. According to our energy estimates, there should not be any reliable bending on A_nT_m. In contrast, in all other cases there was a pronounced roll-like bending into the major groove in the chemical symmetry region. Calculations of the junction between the anomalous and ordinary B-like structure for G_nT_m and C_nA_m have shown that there is an equilibrium bending with a tilt component towards the chain having the anomalous structure at the 5′-end. From our calculations it is impossible to determine precisely the direction of bending, though it can be suggested that the roll component of bending might be directed towards the major groove. The anomalous structure is the main reason of bending; alternations of pyrimidines and purines can modulate the value and the direction of equilibrium bending (only the value in the case of self-complemantary fragments).

The results are consistent with the experimental data and promote a better understanding of the problem of DNA bending.     

Molecular mechanics calculations of dA12·dT12 and of the curved molecule d(GCTCGAAAAA)4·d(TTTTTCGAGC)4

Using the AMBER software package (Weiner and Kollman 1981) substantially modified for electrostatic contributions, the structural energies of the double-stranded oligonucleotides dA₁₂·dT₁₂ and d(GCTCGAAAAA)₄·d(TTTTTCGAGC)₄ were minimized. Using various starting structures for the molecule dA₁₂·dT₁₂, one final structure is obtained which possesses the experimentally determined properties of poly(dA)·poly(dT). This structure is an A-form-B-form-hybrid structure similar to that of Arnott et al. (1983). The dA-strand is similar to an A-form while the dT-strand is similar to normal B-form. This structure and separately optimized B-form sequence stretches were used to construct the double-stranded fragment d(GCTCGAAAAA)₄ which again was optimized. This sequence, when imbedded in a DNA fragment as contiguous repeats, shows a gel migration anomaly which has been interpreted as stable curvature of the DNA (Diekmann 1986). The calculated structure of this sequence indeed has a curved helix axis and is discussed as a model for curved DNA. A theoretical formalism is presented which allows one to calculate the structural parameters of any nucleic acid double helix in two different geometrical representations. This formalism is used to determine the parameters of the base-pair orientations of the curved structure in terms of wedge as well as cylindrical parameters. In the structural model presented here, the curvature of the helix axis results from an alternation of two different DNA structures in which the base-pairs possess different angles with the helix axis (cylinder tilt). Resulting from geometric restraints, a negative cylinder tilt angle correlates strongly with the closing of the minor groove (wedge roll). The blocks with different structure are not exactly coincident with the dA₅-blocks and the B-DNA stretches. Within the dA₅ block, base-pair tilt and wedge roll adopt large values which proceed into the 3 flanking B-DNA sequence by about one base-pair. These properties of the structure calculated here are discussed in terms of different models explaining DNA curvature.     

Triplet repeat DNA structures and human genetic disease: dynamic mutations from dynamic DNA

Fourteen genetic neurodegenerative diseases and three fragile sites have been associated with the expansion of (CTG)_n•(CAG)_n, (CGG)_n•(CCG)_n, or (GAA)_n•(TTC)_n repeat tracts. Different models have been proposed for the expansion of triplet repeats, most of which presume the formation of alternative DNA structures in repeat tracts. One of the most likely structures, slipped strand DNA, may stably and reproducibly form within triplet repeat sequences. The propensity to form slipped strand DNA is proportional to the length and homogeneity of the repeat tract. The remarkable stability of slipped strand DNA may, in part, be due to loop-loop interactions facilitated by the sequence complementarity of the loops and the dynamic structure of three-way junctions formed at the loop-outs.     

Survey of plant short tandem DNA repeats

Length variations in simple sequence tandem repeats are being given increased attention in plant genetics. Some short tandem repeats (STRs) from a few plant species, mainly those at the dinucleotide level, have been demonstrated to show polymorphisms and Mendelian inheritance. In the study reported here a search for all of the possible STRs ranging from mononucleotide up to tetranucleotide repeats was carried out on EMBL and GenBank DNA sequence databases of 3026 kb nuclear DNA and 1268 kb organelle DNA in 54 and 28 plant species (plus algae), respectively. An extreme rareness of STRs (4 STRs in 1268 kb DNA) was detected in organelle compared with nuclear DNA sequences. In nuclear DNA sequences, (AT)_n sequences were the most abundant followed by (A)_n · (T)_n, (AG)_n · (CT)_n, (AAT)_n · (ATT)_n, (AAC)_n · (GTT), (AGC)_n · (GCT)_n, (AAG)_n · (CTT)_n, (AATT)_n · (TTAA)_n, (AAAT)_n · (ATTT)_n and (AC)_n · (GT)_n sequences. A total of 130 STRs were found, including 49 (AT)_n sequences in 31 species, giving an average of 1 STR every 23.3 kb and 1 (AT)_n STR every 62 kb. An abundance comparable to that for the dinucleotide repeat was observed for the tri- and tetranucleotide repeats together. On average, there was 1 STR every 64.6 kb DNA in monocotyledons versus 1 every 21.2 kb DNA in dicotyledons. The fraction of STRs that contained G-C basepairs increased as the G+C contents went up from dicotyledons, monocotyledons to algae. While STRs of mono-, di- and tetranucleotide repeats were all located in non coding regions, 57% of the trinucleotide STRs containing G-C basepairs resided in coding regions.     

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.     

GATC motifs may alter the conformation of DNA depending on sequence context and N6-adenine methylation status: possible implications for DNA-protein recognition

As part of our analysis of the role of a uniquely clustered set of dam methylation sites (the motif GATC) within the origin of DNA replication in Escherichia coli, we have studied the effect of GATCs in various methylation states on the intrinsic curvature of DNA. We have designed a set of DNA linkers and used commercially available linkers containing GATC motifs. The linkers were ligated and the electrophoretic mobility of the resulting multimers in different states of methylation was tested relative to reference fragments. We report that properly phased GATCs in certain sequence environments modulate DNA curvature and that these effects may be enhanced by N⁶-adenine methylation of the GATCs. These structural alterations may in turn affect DNA-protein interactions, especially those involving proteins that rely on both primary sequence and structure for recognition. We present an example, where introduction of a GATC within an integration host factor (IHF) binding site, which does not alter the consensus sequence, reduces the binding affinity of the protein for the modified site. Received: 16 December 1997 / Accepted: 24 February 1998     

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.     

Stathmin and Interfacial Microtubule Inhibitors Recognize a Naturally Curved Conformation of Tubulin Dimers

Tubulin is able to switch between a straight microtubule-like structure and a curved structure in complex with the stathmin-like domain of the RB3 protein (T₂RB3). GTP hydrolysis following microtubule assembly induces protofilament curvature and disassembly. The conformation of the labile tubulin heterodimers is unknown. One important question is whether free GDP-tubulin dimers are straightened by GTP binding or if GTP-tubulin is also curved and switches into a straight conformation upon assembly. We have obtained insight into the bending flexibility of tubulin by analyzing the interplay of tubulin-stathmin association with the binding of several small molecule inhibitors to the colchicine domain at the tubulin intradimer interface, combining structural and biochemical approaches. The crystal structures of T₂RB3 complexes with the chiral R and S isomers of ethyl-5-amino-2-methyl-1,2-dihydro-3-phenylpyrido[3,4-b]pyrazin-7-yl-carbamate, show that their binding site overlaps with colchicine ring A and that both complexes have the same curvature as unliganded T₂RB3. The binding of these ligands is incompatible with a straight tubulin structure in microtubules. Analytical ultracentrifugation and binding measurements show that tubulin-stathmin associations (T₂RB3, T₂Stath) and binding of ligands (R, S, TN-16, or the colchicine analogue MTC) are thermodynamically independent from one another, irrespective of tubulin being bound to GTP or GDP. The fact that the interfacial ligands bind equally well to tubulin dimers or stathmin complexes supports a bent conformation of the free tubulin dimers. It is tempting to speculate that stathmin evolved to recognize curved structures in unassembled and disassembling tubulin, thus regulating microtubule assembly.     

Electroporation of curved lipid membranes in ionic strength gradients

A thermodynamic theory for the membrane electroporation of curved membranes such as those of lipid vesicles and cylindrical membrane tubes has been developed. The theory covers in particular the observation that electric pore formation and shape deformation of vesicles and cells are dependent on the salt concentration of the suspending solvent. It is shown that transmembrane salt gradients can appreciably modify the electrostatic part of Helfrich's spontaneous curvature, elastic bending rigidity and Gaussian curvature modulus of charged membranes. The Gibbs reaction energy of membrane electroporation can be explicitely expressed in terms of salt gradient-dependent contributions of bending, the ionic double layers and electric surface potentials and dielectric polarisation of aqueous pores. In order to cover the various physical contribution to the chemical process of electroporation-resealing, we have introduced a generalised chemophysical potential covering all generalised forces and generalised displacements in terms of a transformed Gibbs energy formalism. Comparison with, and analysis of, the data of electrooptical relaxation kinetic studies show that the Gibbs reaction energy terms can be directly determined from turbidity dichroism (Planck's conservative dichroism). The approach also quantifies the electroporative cross-membrane material exchange such as electrolyte release, electrohaemolysis of red blood cells or uptake of drugs and dyes and finally gene DNA by membrane electroporation.     

Effect of Flanking Sequences on the Right- to Left- Handed Transition of a (dA-dT)n Tract in Supercoiled DNA

Abstract

Alternating (dA-dT)_n sequences in supercoiled DNA may undergo a transition to a left-handed conformation in the presence of Ni²⁺ ions and high NaCl concentration (Nejedlý, K., Klysik, J. and Pale?ek, E, FEBS Lett. 243, 313–317 (1989)). In this work we have found that ionic conditions necessary for the B-to-Z transition are strongly dependent on the sequences flanking the (dA-dT)_n tract. In particular, the presence of 5′- homopyrimidine (C3) and 3′-homopurine (G4) blocks adjacent to the tract were found to facilitate the transition to the left- handed form. Within a constant sequence context it was found that the ionic strength required to promote the transition was inversely proportional to the length of the (dA-dT) _n sequence.     

Base Sequence Effects in Double Helical DNA. I. Potential Energy Estimates of Local Base Morphology

Abstract

A series of potential energy calculations have been carried out to estimate base sequence dependent structural differences in B-DNA. Attention has been focused on the simplest dimeric fragments that can be used to build long chains, computing the energy as a function of the orientation and displacement of the 16 possible base pair combinations within the double helix. Calculations have been performed, for simplicity, on free base pairs rather than complete nucleotide units. Conformational preferences and relative flexibilities are reported for various combinations of the roll, tilt, twist, lateral displacement, and propeller twist of individual residues. The predictions are compared with relevant experimental measures of conformation and flexibility, where available. The energy surfaces are found to fit into two distinct categories, some dimer duplexes preferring to bend in a symmetric fashion and others in a skewed manner. The effects of common chemical substitutions (uracil for thymine, 5-methyl cytosine for cytosine, and hypoxanthine for guanine) on the preferred arrangements of neighboring residues are also examined, and the interactions of the sugar-phosphate backbone are included in selected cases. As a first approximation, long range interactions between more distant neighbors, which may affect the local chain configuration, are ignored. A rotational isomeric state scheme is developed to describe the average configurations of individual dimers and is used to develop a static picture of overall double helical structure. The ability of the energetic scheme to account for documented examples of intrinsic B-DNA curvature is presented, and some new predictions of sequence directed chain bending are offered.     

DNA Topological Context Affects Access to Eukaryotic DNA Topoisomerase I

Abstract

We have analyzed the reactivity of a 217 base pair segment of the intrinsically curved Crithidia fasciculata kinetoplast DNA towards eukaryotic DNA topoisomerase I. The substrates were open [linear fragment and nicked circle] and closed minidomains [closed relaxed circle and circles with linking differences of ?1 and ?2], We interpreted the results with the aid of a model that was used to predict the structures of the topoisomers. The modelling shows that the ΔLk(?l) form is unusually compact because of the curvature in the DNA. To determine the role of sequence-directed curvature in both the experimental and modeling studies, controls were examined in which the curved Crithidia sequence was replaced by an uncurved sequence obtained from the plasmid pBR322.

Reactivity of the Crithidia DNA [as analyzed both by the cleavage and the topoisomerization reactions] markedly varied among the DNA forms: (i) the hierarchy of overall reactivity observed is: linear fragment > nicked circular, closed circular [ΔLk(O)], interwound [ΔLk(?2)] > bent interwound [ΔLk(?l)]; (ii) the intensity of several cleavage positions differs among DNA forms.

The results show that eukaryotic DNA topoisomerase I is very sensitive to the conformation of the substrates and that its reactivity is modulated by the variation of the compactness of the DNA molecule. The C. fasciculata sequence contains a highly curved segment that determines the conformation of the closed circle in a complex way.     

A Quantitative Measure of DNA Curvature Enabling the Comparison of Predicted Structures

Abstract

A growing body of data indicates that the equilibrium structures of some DNA fragments are curved and that curvature is sequence-directed. We describe a quantitative measure of DNA curvature that can be used for evaluating and comparing current proposed models for the molecular basis of DNA curvature. We demonstrate that this measure, in conjunction with any given prediction model, enables both the comparison of experimental data to predictions and the scanning of nucleotide sequence databases for potential curved regions.     

A theoretical model for the prediction of sequence-dependent nucleosome thermodynamic stability

A theoretical model for predicting nucleosome thermodynamic stability in terms of DNA sequence is advanced. The model is based on a statistical mechanical approach, which allows the calculation of the canonical ensemble free energy involved in the competitive nucleosome reconstitution. It is based on the hypothesis that nucleosome stability mainly depends on the bending and twisting elastic energy to transform the DNA intrinsic superstructure into the nucleosomal structure. The ensemble average free energy is calculated starting from the intrinsic curvature, obtained by integrating the dinucleotide step deviations from the canonical B-DNA and expressed in terms of a Fourier series, in the framework of first-order elasticity. The sequence-dependent DNA flexibility is evaluated from the differential double helix thermodynamic stability. A large number of free-energy experimental data, obtained in different laboratories by competitive nucleosome reconstitution assays, are successfully compared to the theoretical results. They support the hypothesis that the stacking energies are the major factor in DNA rigidity and could be a measure of DNA stiffness. A dual role of DNA intrinsic curvature and flexibility emerges in the determination of nucleosome stability. The difference between the experimental and theoretical (elastic) nucleosome-reconstitution free energy for the whole pool of investigated DNAs suggests a significant role for the curvature-dependent DNA hydration and counterion interactions, which appear to destabilize nucleosomes in highly curved DNAs. This model represents an attempt to clarify the main features of the nucleosome thermodynamic stability in terms of physical-chemical parameters and suggests that in molecular systems with a large degree of complexity, the average molecular properties dominate over the local features, as in a statistical ensemble.     

Thermodynamic Framework of the Interaction between Protein and Solvent Drives Protein Folding

Abstract

We use internal coordinate molecular mechanics calculations to study the impact of abasic sites on the conformation and the mechanics of the DNA double helix. Abasic sites, which are common mutagenic lesions, are shown to locally modify both the groove geometry and the curvature of DNA in a sequence dependent manner. By controlled twisting and bending, it is also shown that these lesions modify the deformability of the duplex, generally increasing its flexibility, but again to an extent which depends on the nature of the abasic site and on the surrounding base sequence. Both the conformational and mechanical influence of this type of DNA damage may be significant for recognition and repair mechanisms.     

Computer simulation of the interaction of α3hclix of λ Cro protein with DNA and origin of sequence specific recognition

The conformation of the α3 helix of Cro protein (residues 27–36) of bacteriophageλ is optimised by the damped least square minimization technique, with the steric constraint that Cα atom positions should match the crystallographic data available to date. On the basis of minimization of total interaction and conformation energy, models for complexes of this peptide sequence with heptanucleotide duplexes from native and altered O_R3 operator are obtained in the major groove of B DNA. Analysis of the energetics for 3 sequences of the DNA show that binding strength is derived mainly from the interaction of side chains of the peptide with DNA. Sequence specificity (maximum difference in binding energy for different DNA sequences) is due to hydrogen bonding interaction. A small amount of sequence specificity is derived from non-bonded interaction also. Stereochemical aspects of peptide DNA interaction and their role in DNA recognition are discussed in this paper.