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.     

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.     

Anisotropic elastic bending models of DNA

Simplified elastic rod models of DNA were developed in which the rigidity of DNA is sequence dependent and asymmetrical, i.e. the bending is facilitated towards the major groove. By subjecting the models to bending load in various directions perpendicular to the longitudinal axis of DNA, the bending deformation and the average conformation of the models can be estimated using finite element methods. Intrinsically curved sequence motifs [(aaaattttgc)_n, (tctctaaaaaatatataaaaa)_n] are found to be curved by this modelling procedure whereas the average conformation of homopolymers and straight motifs [(a)_n, (atctaatctaacacaacaca)_n] show negligible or no curvature. This suggests that sequence dependent asymmetric rigidity of DNA can provide an explanation in itself for intrinsic DNA curvature. The average rigidity of various DNA sequences was calculated and a good correlation was found with such quantities as the free energy change upon the binding of the Cro repressor, the base stacking energy and the thermal fluctuations at room temperature.     

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.     

DNA bending by charged peptides: electrophoretic and spectroscopic analyses

We are testing the idea that placement of fixed charges near one face of the DNA double helix can induce DNA bending by a purely electrostatic mechanism. If stretching forces between DNA phosphates are significant, fixed charges should induce DNA bending by asymmetrically modulating these forces. We have previously tested this hypothesis by adding charged residues to small bZIP DNA binding peptides and monitoring DNA bending using electrophoretic phasing assays. Our results were consistent with an electrostatic model of DNA bending in predicted directions. We now confirm these observations with fluorescence resonance energy transfer (FRET). Using a "U"-shaped DNA probe, we report that DNA bending by charged bZIP peptides is readily detected by FRET. We further show that charged bZIP peptides cause DNA bending rather than DNA twisting.     

Bending and curvature calculations in B-DNA.

A simple program, BEND, has been written to calculate the magnitude of local bending and macroscopic curvature at each point along an arbitrary B-DNA sequence, using any desired bending model that specifies values of twist, roll and tilt as a function of sequence. The program has been used to evaluate six different DNA bending models in three categories. Two are bent non-A-tract models: (a) A new model based on the nucleosome positioning data of Satchwell et al 1986 (J. Mol. Biol. 191, 659-675), (b) The model of Calladine et al 1988 (J. Mol. Biol. 201, 127-137). Three are bent A-tract models: (c) The wedge model of Bolshoy et al 1991 (Proc. Natl. Acad. Sci. USA 88, 2312-2316), (d) The model of Cacchione et al 1989 (Biochem. 28, 8706-8713), (e) A reversed version of model (b). The last is a junction model: (f) The model of Koo & Crothers 1988 (Proc. Natl. Acad. Sci. USA 85, 1763-1767). Although they have widely different assumptions and values for twist, roll and tilt, all six models correctly predict experimental A-tract curvature as measured by gel retardation and cyclization kinetics, but only the new nucleosome positioning model is successful in predicting curvature in regions containing phased GGGCCC sequences. This model--showing local bending at mixed sequence DNA, strong bends at the sequence GGC, and straight, rigid A-tracts--is the only model consistent with both solution data from gel retardation and cyclization kinetics and structural data from x-ray crystallography.     

Direct measurements of the nucleosome-forming preferences of periodic DNA motifs challenge established models

Several periodic motifs have been implicated in facilitating the bending of DNA around the histone core of the nucleosome. For example, di-nucleotides AA/TT/TA and GC at ∼10-bp periods, but offset by 5 bp, are found with higher-than-expected occurrences in aligned nucleosomal DNAs in vitro and in vivo. Additionally, regularly oscillating period-10 trinucleotide motifs non-T, A/T, G and their complements have been implicated in the formation of regular nucleosome arrays. The effects of these periodic motifs on nucleosome formation have not been systematically tested directly by competitive reconstitution assays. We show that, in general, none of these period-10 motifs, except TA, in certain sequence contexts, facilitates nucleosome formation. The influence of periodic TAs on nucleosome formation is appreciable; with some of the 200-bp DNAs out-competing bulk nucleosomal DNA by more than 400-fold. Only the nucleotides immediately flanking TA influence its nucleosome-forming ability. Period-10 TA, when flanked by a pair of permissive nucleotides, facilitates DNA bending through compression of the minor groove. The free energy change for nucleosome formation decreases linearly with the number of consecutive TAs, up to eight. We suggest how these data can be reconciled with previous findings.     

Solution measurement of DNA curvature in papillomavirus E2 binding sites

‘Indirect readout’ refers to the proposal that proteins can recognize the intrinsic three-dimensional shape or flexibility of a DNA binding sequence apart from direct protein contact with DNA base pairs. The differing affinities of human papillomavirus (HPV) E2 proteins for different E2 binding sites have been proposed to reflect indirect readout. DNA bending has been observed in X-ray structures of E2 protein–DNA complexes. X-ray structures of three different E2 DNA binding sites revealed differences in intrinsic curvature. DNA sites with intrinsic curvature in the direction of protein-induced bending were bound more tightly by E2 proteins, supporting the indirect readout model. We now report solution measurements of intrinsic DNA curvature for three E2 binding sites using a sensitive electrophoretic phasing assay. Measured E2 site curvature agrees well the predictions of a dinucleotide model and supports an indirect readout hypothesis for DNA recognition by HPV E2.     

E1 recognition sequences in the bovine papillomavirus type 1 origin of DNA replication: interaction between half sites of the inverted repeats.

The E1 protein encoded by bovine papillomavirus type 1 (BPV-1) is required for viral DNA replication, and it binds site specifically to an A/T-rich palindromic sequence within the viral origin of replication. The protein is targeted to this site through cooperative interactions and binding with the virus-encoded E2 protein. To explore the nature of the E1 binding site, we inserted a series of homologous DNA linkers at the center of dyad symmetry within the E1 recognition palindrome. The effects of these modifications indicated that the E1 recognition palindrome can be separated into functional half sites. The series of insertions manifest a phasing relationship with respect to the wild-type BPV-1 genome in that greater biological activity was measured when full integral turns of the DNA helix separated the palindrome than when the separations were half-turns. This phasing pattern of activity was observed to occur in a variety of biological phenotypes, including transformation efficiency, stable plasmid copy number in cell lines established from pooled foci, and transient replication of full-length viral genomes. For replication reporter constructs where E1 and E2 are supplied in trans by the respective expression vectors, distance between the half sites seems to play a major role, yet the phasing relationships are measurable. DNase I protection studies showed that E1 bound very poorly to the construct containing a 5-bp linker, and binding was close to the wild-type level for the 10-bp insertion, consistent with a requirement for a phasing function between half sites with a modulus of 10 bp. Binding to the 15- and 20-bp insertion mutants was weak, but only for the 20-bp insertions was protection over both halves of the palindrome measurable. As it had been previously reported that the 18-bp palindrome contains sufficient nucleotide sequence information for E1 binding, we speculate that a minimal E1 recognition motif is presented in each half site. A comparison between this sequence and that of an upstream region that also binds E1 (the E2RE1 region) revealed a common pentanucleotide motif of APyAAPy. Mutants with substitutions of the ATAAT elements within E2RE1 failed to bind E1 protein. We present models for how repeats of the pentanucleotide sequence may coordinate E1 binding at the dyad symmetry axis of the origin and compare the DNA sequence organization of BPV-1 with those of the simian virus 40 and polyomaviruses at their origins of DNA replication.     

64 Hoogsteen or not Hoogsteen? Iodine-125 radioprobing of the p53-induced DNA deformations

Radioprobing is suitable for tracing the DNA and RNA trajectories in nucleoprotein complexes in solution. The method is based on the analysis of the single-strand breaks produced by decay of iodine-125 incorporated in the C5 position of cytosine (Karamychev et al., 1999, 2012). Here, we used radioprobing to study the conformation of DNA in complex with the DNA binding domain (DBD) of the tumor-suppressor protein p53. Two recently crystallized DNA-p53 DBD complexes have different conformations of the CATG motifs: one with the Hoogsteen A:T pairs (Kitayner et al., 2010) and the other with the Watson–Crick pairs (Chen et al., 2010). The two complexes differ in the sequence of the central YYY|RRR junction: the first one has the C|G step and the second has the T|A step. Thus, it is interesting to apply the radioprobing method to the two DNA sequences used in crystallography to see if the local changes (T|A to C|G) in the center of the p53 response element would produce significant distortions in the CATG motifs. To this aim, the iodine-containing cytosine was incorporated in the duplexes containing p53-binding sites, in one of the two CATG motifs and the frequencies of DNA breaks were analyzed. Frequencies of breaks are negatively correlated with the iodine–sugar distances, thus, one can evaluate the changes in these distances upon DNA binding to a protein. The radioprobing distances obtained for both DNA sequences proved to be consistent with the Watson–Crick structure observed by Chen et al. (2010). We did not find any evidence of the Hoogsteen A:T base pair formation in the DNA-p53 DBD complexes in solution using our radioprobing method. The most significant changes in the break frequency distributions were detected in the central segment of the p53 binding site, YYY|RRR, which are consistent with an increase in DNA twisting in this region and local DNA bending and sliding (Nagaich et al., 1999). We interpret these p53-induced DNA deformations in the context of p53 binding to nucleosomal DNA (Sahu et al., 2010).     

牛乳腺核基质附着区多重顺式作用序列表明其功能的多样性

高等真核细胞的染色体DNA通过基质结合区(MAR)不时地与核基质特异性结合而组织成一种空间环状结构。为了研究以DNA套环形式附着于核基质上的DNA序列的特性,从处于泌乳期的乳腺组织中克隆了多个MAR DNA序列。体外结合实验表明,这些序列能够同核基质蛋白共结合成不溶性的复合物,这些复合物可较容易的通过离心去除。其中,两个MAR序列中包含有TL、CA—和GA—阻断以及ATTA基序。这两个序列中含有多个复制／转录因子的结合位点、增强子基序、多个完全的和非完全的反向重复序列以及潜在的DNA弯曲核心序列样结构。同一DNA序列中存在不同元件的组合可能说明在控制一系列细胞的发育过程中,它们可能发挥有正的或负的调控元件的功能。     

DNA unwinding produced by site-specific intrastrand cross-links of the antitumor drug cis-diamminedichloroplatinum(II)

The DNA unwinding produced by specific adducts of the antitumor drug cis-diamminedichloroplatinum(II) has been quantitatively determined. Synthetic DNA duplex oligonucleotides of varying lengths with two base pair cohesive ends were synthesized and characterized that contained site-specific intrastrand N7-purine/N7-purine cross-links. Included are cis-[Pt(NH3)2[d(GpG)]], cis-[Pt(NH3)2(d(ApG)]], and cis-[Pt(NH3)2[d(GpTpG)]] adducts, respectively referred to as cis-GG, cis-AG, and cis-GTG. Local DNA distortions at the site of platination were amplified by polymerization of these monomers and quantitatively evaluated by using polyacrylamide gel electrophoresis. The extent of DNA unwinding was determined by systematically varying the interplatinum distance, or phasing, in polymers containing the adducts. The multimer that migrates most slowly gives the optimal phasing for cooperative bending, from which the degree of unwinding can be obtained. We find that the cis-GG and cis-AG adducts both unwind DNA by 13 degrees, while the cis-GTG adduct unwinds DNA by 23 degrees. In addition, experiments are presented that support previous studies revealing that a hinge joint forms at the sites of platination in DNA molecules containing trans-GTG adducts. On the basis of an analysis of the present and other published studies of site-specifically modified DNA, we propose that local duplex unwinding is a major determinant in the recognition of DNA damage by the Escherichia coli (A)BC excinuclease. In addition, local duplex unwinding of 13 degrees and bending by 35 degrees are shown to correlate well with the recognition of platinated DNA by a previously identified damage recognition protein (DRP) in human cells.     

Long-range organization and sequence-directed curvature of Xenopus laevis satellite 1 DNA.

We have investigated the long-range organization and the intrinsic curvature of satellite 1 DNA, an unusual tandemly-repeated DNA family of Xenopus laevis presenting sequence homologies to SINEs. PFGE was used in combination with frequent-cutter restriction enzymes not likely to cut within satellite 1 DNA and revealed that almost all the repeating units are tandemly organized to form large arrays (200 kb to 2 Mb) that are marked by restriction length polymorphism and contain intra-array domains of sequence variation. Besides that, we have analysed the secondary structure of satellite 1 DNA by computer modelling. Theoretical maps of curvature obtained from three independent models of DNA bending (the dinucleotide wedge model of Trifonov, the junction model of Crothers and the model of de Santis) showed that satellite 1 DNA is intrinsically curved and these results were confirmed experimentally by polyacrylamide gel electrophoresis. Moreover, we observed that this bending element is highly conserved among all the members of the satellite 1 DNA family that are accessible to analysis. A potential genetic role for satellite 1 DNA based on this unusual structural feature is discussed.