CA/TG sequence at the 5' end of oligo(A)-tracts strongly modulates DNA curvature

Molecular stripping underpins derepression of a toxin-antitoxin system

Transcription factors control gene expression; among these, transcriptional repressors must liberate the promoter for derepression to occur. Toxin-antitoxin (TA) modules are bacterial elements that autoregulate their transcription by binding the promoter in a T:A ratio-dependent manner, known as conditional cooperativity. The molecular basis of how excess toxin triggers derepression has remained elusive, largely because monitoring the rearrangement of promoter-repressor complexes, which underpin derepression, is challenging. Here, we dissect the autoregulation of the Salmonella enterica tacAT3 module. Using a combination of assays targeting DNA binding and promoter activity, as well as structural characterization, we determine the essential TA and DNA elements required to control transcription, and we reconstitute a repression-to-derepression path. We demonstrate that excess toxin triggers molecular stripping of the repressor complex off the DNA through multiple allosteric changes causing DNA distortion and ultimately leading to derepression. Thus, our work provides important insight into the mechanisms underlying conditional cooperativity.

Structural basis of gene regulation by the Grainyhead/CP2 transcription factor family

Grainyhead (Grh)/CP2 transcription factors are highly conserved in multicellular organisms as key regulators of epithelial differentiation, organ development and skin barrier formation. In addition, they have been implicated as being tumor suppressors in a variety of human cancers. Despite their physiological importance, little is known about their structure and DNA binding mode. Here, we report the first structural study of mammalian Grh/CP2 factors. Crystal structures of the DNA-binding domains of grainyhead-like (Grhl) 1 and Grhl2 reveal a closely similar conformation with immunoglobulin-like core. Both share a common fold with the tumor suppressor p53, but differ in important structural features. The Grhl1 DNA-binding domain binds duplex DNA containing the consensus recognition element in a dimeric arrangement, supporting parsimonious target-sequence selection through two conserved arginine residues. We elucidate the molecular basis of a cancer-related mutation in Grhl1 involving one of these arginines, which completely abrogates DNA binding in biochemical assays and transcriptional activation of a reporter gene in a human cell line. Thus, our studies establish the structural basis of DNA target-site recognition by Grh transcription factors and reveal how tumor-associated mutations inactivate Grhl proteins. They may serve as points of departure for the structure-based development of Grh/CP2 inhibitors for therapeutic applications.

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.

Probing sequence-specific DNA flexibility in a-tracts and pyrimidine-purine steps by nuclear magnetic resonance (13)C relaxation and molecular dynamics simulations

Sequence-specific DNA flexibility plays a key role in a variety of cellular interactions that are critical for gene packaging, expression, and regulation, yet few studies have experimentally explored the sequence dependence of DNA dynamics that occur on biologically relevant time scales. Here, we use nuclear magnetic resonance (NMR) carbon spin relaxation combined with molecular dynamics (MD) simulations to examine the picosecond to nanosecond dynamics in a variety of dinucleotide steps as well as in varying length homopolymeric A(n)·T(n) repeats (A(n)-tracts, where n = 2, 4, or 6) that exhibit unusual structural and mechanical properties. We extend the NMR spin relaxation time scale sensitivity deeper into the nanosecond regime by using glycerol and a longer DNA duplex to slow overall tumbling. Our studies reveal a structurally unique A-tract core (for n > 3) that is uniformly rigid, flanked by junction steps that show increasing sugar flexibility with A-tract length. High sugar mobility is observed at pyrimidine residues at the A-tract junctions, which is encoded at the dinucleotide level (CA, TG, and CG steps) and increases with A-tract length. The MD simulations reproduce many of these trends, particularly the overall rigidity of A-tract base and sugar sites, and suggest that the sugar-backbone dynamics could involve transitions in sugar pucker and phosphate backbone BI ↔ BII equilibria. Our results reinforce an emerging view that sequence-specific DNA flexibility can be imprinted in dynamics occurring deep within the nanosecond time regime that is difficult to characterize experimentally at the atomic level. Such large-amplitude sequence-dependent backbone fluctuations might flag the genome for specific DNA recognition.

Benzimidazoles: a minor groove-binding ligand-induced stabilization of triple helix

Nonintercalating minor groove-binding ligands netropsin, Hoechst 33258, and DAPI are reported to destabilize the triplex. Ligands with different substitutions on the phenyl ring of bis- and terbenzimidazoles were evaluated for their effect on the stability of C+.GC triplex and Hoogsteen duplex. We found that newly synthesized benzimidazoles stabilize the triplex as shown by fluorescence and melting studies. Modeling studies showed that these molecules bind in the Watson-Crick minor groove of triplex, which can exert a profound impact on the properties of the host triplex. Circular dichroism-binding studies indicate 5.77 base triplets/ligand as an apparent binding site for bis- and 8.66 for terbenzimidazoles. The stabilization of triplex can be attributed to the protonation of nitrogens and amines of benzimidazoles at pH 5.2 that compensate the negative charge of phosphate backbone to reduce the repulsion between the strands resulting in the stabilization.

Distant neighbor base sequence context effects in human nucleotide excision repair of a benzo[a]pyrene-derived DNA lesion

The effects of non-nearest base sequences, beyond the nucleotides flanking a DNA lesion on either side, on nucleotide excision repair (NER) in extracts from human cells were investigated. We constructed two duplexes containing the same minor groove-aligned 10S (+)-trans-anti-B[a]P-N(2)-dG (G*) DNA adduct, derived from the environmental carcinogen benzo[a]pyrene (B[a]P): 5'-C-C-A-T-C-G*-C-T-A-C-C-3' (CG*C-I), and 5'-C-A-C3-A4-C5-G*-C-A-C-A-C-3' (CG*C-II). We used polyacrylamide gel electrophoresis to compare the extent of DNA bending, and molecular dynamics simulations to analyze the structural characteristics of these two DNA duplexes. The NER efficiencies are 1.6(+/-0.2)-fold greater in the case of the CG*C-II than the CG*C-I sequence context in 135-mer duplexes. Gel electrophoresis and self-ligation circularization experiments revealed that the CG*C-II duplex is more bent than the CG*C-I duplex, while molecular dynamics simulations showed that the unique -C3-A4-C5- segment in the CG*C-II duplex plays a key role. The presence of a minor groove-positioned guanine amino group, the Watson-Crick partner to C3, acts as a wedge; facilitated by a highly deformable local -C3-A4- base step, this amino group allows the B[a]P ring system to produce a more enlarged minor groove in CG*C-II than in CG*C-I, as well as a local untwisting and enlarged and flexible Roll only in the CG*C-II sequence. These structural properties fit well with our earlier findings that in the case of the family of minor groove 10S (+)-trans-anti-B[a]P-N(2)-dG lesions, flexible bends and enlarged minor groove widths constitute NER recognition signals, and extend our understanding of sequence context effects on NER to the neighbors that are distant to the lesion.

Estrogen receptor acting in cis enhances WT and mutant p53 transactivation at canonical and noncanonical p53 target sequences

p53 is a master regulatory, sequence-specific transcription factor that directly controls expression of over 100 genes in response to various stress signals. Transactivation is generally considered to occur through p53 binding to a consensus response element (RE) composed of two 5'-RRRCWWGYYY-3' decamers. Recently, studying the human angiogenesis-related gene FLT1 we discovered that p53 can mediate limited transactivation at a noncanonical 1/2 site and could synergize with the estrogen receptor (ER) acting in cis at a nearby ER 1/2 site. To address the generality of concerted transactivation by p53 and ER, the 1/2 site in the FLT1 promoter was replaced with a variety of 1/2 sites, as well as canonical weak and strong p53 REs of human target genes. The p53 transactivation of all tested sequences was greatly enhanced by ligand-activated ER acting in cis. Furthermore, enhanced transactivation extends to several cancer-associated p53 mutants with altered function, suggesting ER-dependent mutant p53 activity for at least some REs. The enhanced transactivation was also found with p63 and p73. We propose a general synergistic relationship between p53 family and ER master regulators in transactivation of p53 target canonical and noncanonical REs, which might be poorly responsive to p53 on their own. This relationship greatly expands the transcriptional master network regulated by p53 in terms of genes affected and levels of expression and has implications for the appearance and possible treatments of cancer.

The expanding universe of p53 targets

The p53 tumour suppressor is modified through mutation or changes in expression in most cancers, leading to the altered regulation of hundreds of genes that are directly influenced by this sequence-specific transcription factor. Central to the p53 master regulatory network are the target response element (RE) sequences. The extent of p53 transactivation and transcriptional repression is influenced by many factors, including p53 levels, cofactors and the specific RE sequences, all of which contribute to the role that p53 has in the aetiology of cancer. This Review describes the identification and functionality of REs and highlights the inclusion of non-canonical REs that expand the universe of genes and regulation mechanisms in the p53 tumour suppressor network.

Small local variations in B-form DNA lead to a large variety of global geometries which can accommodate most DNA-binding protein motifs

BACKGROUND

An important question of biological relevance is the polymorphism of the double-helical DNA structure in its free form, and the changes that it undergoes upon protein-binding. We have analysed a database of free DNA crystal structures to assess the inherent variability of the free DNA structure and have compared it with a database of protein-bound DNA crystal structures to ascertain the protein-induced variations.

RESULTS

Most of the dinucleotide steps in free DNA display high flexibility, assuming different conformations in a sequence-dependent fashion. With the exception of the AA/TT and GA/TC steps, which are 'A-phobic', and the GG/CC step, which is 'A-philic', the dinucleotide steps show no preference for A or B forms of DNA. Protein-bound DNA adopts the B-conformation most often. However, in certain cases, protein-binding causes the DNA backbone to take up energetically unfavourable conformations. At the gross structural level, several protein-bound DNA duplexes are observed to assume a curved conformation in the absence of any large distortions, indicating that a series of normal structural parameters at the dinucleotide and trinucleotide level, similar to the ones in free B-DNA, can give rise to curvature at the overall level.

CONCLUSION

The results illustrate that the free DNA molecule, even in the crystalline state, samples a large amount of conformational space, encompassing both the A and the B-forms, in the absence of any large ligands. A-form as well as some non-A, non-B, distorted geometries are observed for a small number of dinucleotide steps in DNA structures bound to the proteins belonging to a few specific families. However, for most of the bound DNA structures, across a wide variety of protein families, the average step parameters for various dinucleotide sequences as well as backbone torsion angles are observed to be quite close to the free 'B-like' DNA oligomer values, highlighting the flexibility and biological significance of this structural form.

Influence of magnesium ion on the binding of p53 DNA-binding domain to DNA-response elements

Site-specific recognition and DNA-binding activity of p53 are crucial for its tumour suppressor function. Previous reports have shown that metal ions can affect the specific recognition and DNA-binding activity of p53DBD. Here we firstly report that magnesium ion can bind to the protein and influence its DNA-binding activity. To elucidate the nature and the effect of metal ions in the reaction chemistry, we utilized endogenous tryptophan fluorescence to quantitate the interaction between p53DBD and metal ions. The K(a) value for the binding of Mg(2+) to the protein is 1.88 x 10(3) M(-1). Analysis of the CD data clearly suggested that the binding of magnesium ion induced a subtle conformational change rather than a radical modification of the overall protein architecture. Based on the results of electrophoretic mobility shift assays and fluorescence experiments, we concluded that the binding of Mg(2+) significantly stimulated the binding of the protein to DNA in a sequence-independent manner, which differed from that of zinc ions in a sequence-specific manner. Based on these results and the fact that Mg(2+) exists at relatively high concentration in the cell, we propose that Mg(2+) is one of potential factors to affect or regulate the transactivation of p53.

Crystal structure of a p53 core tetramer bound to DNA

The tumor suppressor p53 regulates downstream genes in response to many cellular stresses and is frequently mutated in human cancers. Here, we report the use of a crosslinking strategy to trap a tetrameric p53 DNA binding domain (p53DBD) bound to DNA and the X-ray crystal structure of the protein/DNA complex. The structure reveals that two p53DBD dimers bind to B form DNA with no relative twist and that a p53 tetramer can bind to DNA without introducing significant DNA bending. The numerous dimer-dimer interactions involve several strictly conserved residues thus suggesting a molecular basis for p53DBD-DNA binding cooperativity. Surface residue conservation of the p53DBD tetramer bound to DNA highlights possible regions of other p53 domain or p53 cofactor interactions.

Monovalent cation binding by curved DNA molecules containing variable numbers of a-tracts

Monovalent cation binding by DNA A-tracts, runs of four or more contiguous adenine or thymine residues, has been determined for two curved approximately 200 basepair (bp) restriction fragments, one taken from the M13 origin of replication and the other from the VP1 gene of SV40. These two fragments have previously been shown to contain stable, centrally located bends of 44 degrees and 46 degrees , respectively, located within approximately 60 bp "curvature modules" containing four or five irregularly spaced A-tracts. Transient electric birefringence measurements of these two fragments, sequence variants containing reduced numbers of A-tracts in the SV40 curvature module or changes in the residues flanking the A-tracts in the M13 curvature module, have been combined with the free solution electrophoretic mobilities of the same fragments using known equations to estimate the effective charge of each fragment. The effective charge is reduced, on average, by one-third charge for each A-tract in the curvature module, suggesting that each A-tract binds a monovalent cation approximately one-third of the time. Monovalent cation binding to two or more A-tracts is required to observe significant curvature of the DNA helix axis.

Identification of a nisI promoter within the nisABCTIP operon that may enable establishment of nisin immunity prior to induction of the operon via signal transduction

Certain strains of Lactococcus lactis produce the broad-spectrum bacteriocin nisin, which belongs to the lantibiotic class of antimicrobial peptides. The genes encoding nisin are organized in three contiguous operons: nisABTCIP, encoding production and immunity (nisI); nisRK, encoding regulation; and nisFEG, also involved in immunity. Transcription of nisABTCIP and nisFEG requires autoinduction by external nisin via signal transducing by NisRK. This organization poses the intriguing question of how sufficient immunity (NisI) can be expressed when the nisin cluster enters a new cell, before it encounters external nisin. In this study, Northern analysis in both Lactococcus and Enterococcus backgrounds revealed that nisI mRNA was present under conditions when no nisA transcription was occurring, suggesting an internal promoter within the operon. The nisA transcript was significantly more stable than nisI, further substantiating this. Reverse transcriptase PCR analysis revealed that the transcription initiated just upstream from nisI. Fusing this region to a lacZ gene in a promoter probe vector demonstrated that a promoter was present. The transcription start site (TSS) of the nisI promoter was mapped at bp 123 upstream of the nisI translation start codon. Ordered 5' deletions revealed that transcription activation depended on sequences located up to bp -234 from the TSS. The presence of poly(A) tracts and computerized predictions for this region suggested that a high degree of curvature may be required for transcription initiation. The existence of this nisI promoter is likely an evolutionary adaptation of the nisin gene cluster to enable its successful establishment in other cells following horizontal transfer.

Comparative binding of p53 to its promoter and DNA recognition elements

Tumor suppressor p53 is a transcription factor that transactivates a wide range of genes, including those in DNA repair, cell cycle arrest, apoptosis and its own degradation. To estimate the role of selectivity in binding to its promoters, we measured the binding affinities of a tetrameric p53 construct (p53CT) in vitro with 20 of its recognition elements from a variety of representative genes. The binding of full length p53 to four representative sequences exactly paralleled the affinities to p53CT. The binding of p53 to different recognition elements was co-operative and the affinities varied by up to 50-fold. p53 bound with high affinity to the recognition elements of all the genes involved in cell cycle arrest and some of the genes in apoptosis. All of the lower affinity-binding sites were in genes involved in apoptosis. Our quantitative-binding data were in agreement with published cell-based assays. The regulation of p53 activity is in part determined through the specificity of its DNA-binding interactions.

Structural and transcriptional features of Bombus terrestris satellite DNA and their potential involvement in the differentiation process

A unique satellite DNA family was characterized in the genome of the bumble bee, Bombus terrestris. Sequence analysis revealed that it contains two wide palindromes of about 160 and 190 bp, respectively, that span 75% of the repeated unit. One feature of this satellite DNA is that it accounts for different amounts of genomic DNA in males and females. The DNA curvature and bendability were determined by migration on PAGE and by computer analysis. It has been correlated with the presence of dA/dT stretches repeated in phase with the helix turn and with the presence of the deformable dinucleotide CA-TG embedded in some of these A-T-rich regions. Transcription of the satellite DNA was also analyzed by Northern blot hybridization and RT-PCR. Multimeric transcripts spanning several satellite DNA units were found in RNA samples from males, workers, and queens. These transcripts resulted from a specific transcription occurring on one DNA strand in the embryos or on both DNA strands in imagoes. The involvement of DNA curvature in the organization of the satellite DNA and the function of the satellite transcripts is discussed.

Permanganate oxidation reactions of DNA: perspective in biological studies

KMnO4 has been well known as a powerful chemical probe for numerous applications in biological fields, particularly for those used in conformational studies of DNA. The KMnO4 assay provides essential information for understanding biochemical processes and detecting aberrant DNA, which is associated with many genetic diseases. Elegant examples are sequencing techniques, foot-printing assays for transcriptional studies, an interference method for hormone receptor binding assays as well as DNA conformational studies of Z-DNA, Z-Z junctions, hairpins, curvatures, short nucleotide base repeats, binding of intercalators and groove binders, etc. Recently, KMnO4 has been successfully applied to detect single base changes and mutations in DNA (chemical cleavage of mismatch method, CCM) as well as other types of base damage (8-oxoguanine and thymine dimers). This paper aims to review the usefulness and limitations of the permanganate oxidation reaction used in various biological studies of DNA.

Molecular dynamics simulations of DNA curvature and flexibility: helix phasing and premelting

Recent studies of DNA axis curvature and flexibility based on molecular dynamics (MD) simulations on DNA are reviewed. The MD simulations are on DNA sequences up to 25 base pairs in length, including explicit consideration of counterions and waters in the computational model. MD studies are described for ApA steps, A-tracts, for sequences of A-tracts with helix phasing. In MD modeling, ApA steps and A-tracts in aqueous solution are essentially straight, relatively rigid, and exhibit the characteristic features associated with the B'-form of DNA. The results of MD modeling of A-tract oligonucleotides are validated by close accord with corresponding crystal structure results and nuclear magnetic resonance (NMR) nuclear Overhauser effect (NOE) and residual dipolar coupling (RDC) structures of d(CGCGAATTCGCG) and d(GGCAAAAAACGG). MD simulation successfully accounts for enhanced axis curvature in a set of three sequences with phased A-tracts studied to date. The primary origin of the axis curvature in the MD model is found at those pyrimidine/purine YpR "flexible hinge points" in a high roll, open hinge conformational substate. In the MD model of axis curvature in a DNA sequence with both phased A-tracts and YpR steps, the A-tracts appear to act as positioning elements that make the helix phasing more precise, and key YpR steps in the open hinge state serve as curvature elements. Our simulations on a phased A-tract sequence as a function of temperature show that the MD simulations exhibit a premelting transition in close accord with experiment, and predict that the mechanism involves a B'-to-B transition within A-tracts coupled with the prediction of a transition in key YpR steps from the high roll, open hinge, to a low roll, closed hinge substate. Diverse experimental observations on DNA curvature phenomena are examined in light of the MD model with no serious discrepancies. The collected MD results provide independent support for the "non-A-tract model" of DNA curvature. The "junction model" is indicated to be a special case of the non-A-tract model when there is a Y base at the 5' end of an A-tract. In accord with crystallography, the "ApA wedge model" is not supported by MD.

Spectroscopic study of permanganate oxidation reactions of oligonucleotides containing single base mismatches

The electrophoretic gel-based chemical cleavage of the mismatch method gives an incomplete view of the DNA conformational changes induced by a single base mismatch. This spectroscopic study investigates the permanganate oxidation reactions with matched and mismatched DNA under constant and variable temperature conditions. The results, which include the oxidation levels, reaction patterns with isosbestic points, color changes, thermal spectra, spectroscopy derivative, and gel separation and melting temperatures, provide a fundamental background for identification of oligonucleotides containing single base mismatches by chemical means.

Minor Groove Binding DNA Ligands with Expanded A/T Sequence Length Recognition, Selective Binding to Bent DNA Regions and Enhanced Fluorescent Properties

DNA minor groove ligands provide a paradigm for double-stranded DNA recognition, where common structural motifs provide a crescent shape that matches the helix turn. Since minor groove ligands are useful in medicine, new ligands with improved binding properties based on the structural information about DNA-ligand complexes could be useful in developing new drugs. Here, two new synthetic analogues of AT specific Hoechst 33258 5-(4-methylpiperazin-1-yl)-2-[2'-(3,4-dimethoxyphenyl)-5'-benzimidazolyl] benzimidazole (DMA) and 5-(4-methylpiperazin-1-yl)-2-[2'[2''-(4-hydroxy-3-methoxyphenyl)-5' '-benzimidazolyl]-5'-benzimidazolyl] benzimidazole (TBZ) were evaluated for their DNA binding properties. Both analogues are bisubstituted on the phenyl ring. DMA contains two ortho positioned methoxy groups, and TBZ contains a phenolic group at C-4 and a methoxy group at C-3. Fluorescence yield upon DNA binding increased 100-fold for TBZ and 16-fold for DMA. Like the parent compound, the new ligands showed low affinity to GC-rich (K approximately 4 x 10(7) M(-1)) relative to AT-rich sequences (K approximately 5 x 10(8) M(-1)), and fluorescence lifetime and anisotropy studies suggest two distinct DNA-ligand complexes. Binding studies indicate expanded sequence recognition for TBZ (8-10 AT base pairs) and tighter binding (DeltaT(m) of 23 degrees C for d (GA(5)T(5)C). Finally, EMSA and equilibrium binding titration studies indicate that TBZ preferentially binds highly hydrated duplex domains with altered A-tract conformations d (GA(4)T(4)C)(2) (K= 3.55 x 10(9) M(-1)) and alters its structure over d (GT(4)A(4)C)(2) (K = 3.3 x 10(8) M(-1)) sequences. Altered DNA structure and higher fluorescence output for the bound fluorophore are consistent with adaptive binding and a constrained final complex. Therefore, the new ligands provide increased sequence and structure selective recognition and enhanced fluorescence upon minor groove binding, features that can be useful for further development as probes for chromatin structure stability.

Analysis of DNA bending by transient electric birefringence

Transient electric birefringence has been used to analyze DNA bending in six restriction fragments containing 171, 174, 207, 263, 289, and 471 bp in three different low ionic strength buffers. The target fragments contain sequences corresponding to the apparent bend centers in pUC19 and Litmus 28, previously identified by the circular permutation assay (Strutz, K.; Stellwagen, N. C. Electrophoresis 1996, 17, 989-995). The target fragments migrate anomalously slowly in polyacrylamide gels and exhibit birefringence relaxation times that are shorter than those of restriction fragments of the same size, taken from nonbent regions of the same plasmids. Apparent bend angles ranging from 30 degrees to 41 degrees were calculated for the target fragments by tau-ratio method. The bend angles of four of the target fragments were independent of temperature from 4 degrees C to 20 degrees C, but decreased when the temperature was increased to 37 degrees C. The bend angles of the other two target fragments were independent of temperature over the entire range examined, 4 degrees -37 degrees C. Hence, the thermal stability of sequence-dependent bends in random-sequence DNA is variable. The bend angles of five of the six target fragments were independent of the presence or absence of Mg2+ ions in the solution, indicating most of the target fragments were stably bent or curved, rather than anisometrically flexible. Restriction fragments containing 219 and 224 bp, with sequences somewhat offset from the sequence of the 207 bp fragment, were also studied. Comparison of the tau-ratios of these overlapping fragments allowed both the bend angle and bend position to be independently determined. These methods should be useful for analyzing sequence-dependent bending in other random-sequence DNAs.

A-tract polarity dominate the curvature in flanking sequences

It is well-known, but little understood, that the nucleotide sequences between phased A(4-6)-tracts (at 10-11 bp intervals) have only a slight effect on overall curvature. To explore this phenomenon, we have examined the gel-migration properties of sequences containing both A-tracts as well as G-tracts (i.e., sequences of the form G(n)C(m) or C(n)G(m), n + m > 4) in various relative positioning. We show that the composite bend of these sequences depends on their relative arrangement. When G-tracts are placed between two A-tracts, such that both tracts are repeated in phase to themselves (e.g., G(5)A(6)G(5)A(5)), or adjacent to the 3'-side of A-tracts (e.g., A(6)G(5)N(10)), they have minimal influence on the extent of bending of the composite sequence. When G-tracts are placed one helical repeat away from A-tracts (e.g., G(5)N(5)A(6)N(6)), or are adjacent only to the 5'-side of A-tracts (e.g., G(5)A(6)N(10)) their influence on the composite bend is larger. The differential behavior of AG- versus GA-tracts means that A-tracts influence their flanking sequences in a polar manner. Whereas they suppress, or make constant, the intrinsic bending characteristics of any sequence placed immediately 3' to them (and hence by definition any sequence placed between two phased A-tracts), sequences adjoining them on their 5'-side are free to modulate the overall curvature. We interpret these results as evidence for the dominant nature of the unique and nonuniform structure adopted by tracts of four adenines or more. The effects of A-tracts extend at least five base pairs into the adjoining 3' region. This is further evidence for the complexity of DNA structure and the inadequacy of simple nearest-neighbor models to explain all its manifestations.

Differential transactivation by the p53 transcription factor is highly dependent on p53 level and promoter target sequence

Little is known about the mechanisms that regulate differential transactivation by p53. We developed a system in the yeast Saccharomyces cerevisiae that addresses p53 transactivation capacity from 26 different p53 response elements (REs) under conditions where all other factors, such as chromatin, are kept constant. The system relies on a tightly regulated promoter (rheostatable) that can provide for a broad range of p53 expression. The p53 transactivation capacity toward each 20- to 22-bp-long RE could be ranked by using a simple phenotypic assay. Surprisingly, there was as much as a 1,000-fold difference in transactivation. There was no correlation between the functional rank and statistical predictions of binding energy of the REs. Instead we found that the central sequence element in an RE greatly affects p53 transactivation capacity, possibly because of DNA structural properties. Our results suggest that intrinsic DNA binding affinity and p53 protein levels are important contributors to p53-induced differential transactivation. These results are also relevant to understanding the regulation by other families of transcription factors that recognize several sequence-related response elements and/or have tightly regulated expression. We found that p53 had weak activity towards half the apoptotic REs. In addition, p53 alleles associated with familial breast cancer, previously classified as wild type, showed subtle differences in transactivation capacity towards several REs.

Evidence for a common molecular basis for sequence recognition of N3-guanine and N3-adenine DNA adducts involving the covalent bonding reaction of (+)-CC-1065

The antitumor antibiotic (+)-CC-1065 can alkylate N3 of guanine in certain sequences. A previous high-field 1H NMR study on the (+)-CC-1065d[GCGCAATTG*CGC]2 adduct (* indicates the drug alkylation site) showed that drug modification on N3 of guanine results in protonation of the cross-strand cytosine [Park, H.-J.; Hurley, L. H. J. Am. Chem. Soc. 1997, 119, 629.]. In this contribution we describe a further analysis of the NMR data sets together with restrained molecular dynamics. This study provides not only a solution structure of the (+)-CC-1065(N3-guanine) DNA duplex adduct but also new insight into the molecular basis for the sequence-specific interaction between (+)-CC-1065 and N3-guanine in the DNA duplex. On the basis of NOESY data, we propose that the narrow minor groove at the 7T8T step and conformational kinks at the junctions of 16C17A and 18A19T are both related to DNA bending in the drugDNA adduct. Analysis of the one-dimensional 1H NMR (in H2O) data and rMD trajectories strongly suggests that hydrogen bonding linkages between the 8-OH group of the (+)-CC-1065 A-subunit and the 9G10C phosphate via a water molecule are present. All the phenomena observed here in the (+)-CC-1065(N3-guanine) adduct at 5'-AATTG* are reminiscent of those obtained from the studies on the (+)-CC-1065(N3-adenine) adduct at 5'-AGTTA*, suggesting that (+)-CC-1065 takes advantage of the conformational flexibility of the 5'-TPu step to entrap the bent structure required for the covalent bonding reaction. This study reveals a common molecular basis for (+)-CC-1065 alkylation at both 5'-TTG* and 5'-TTA*, which involves a trapping out of sequence-dependent DNA conformational flexibility as well as sequence-dependent general acid and general base catalysis by duplex DNA.

Base sequence effects in bending induced by bulky carcinogen-DNA adducts: experimental and computational analysis

The covalent binding of bulky mutagenic or carcinogenic compounds to DNA can lead to bending, which could significantly alter the interactions of DNA with critical replication and transcription proteins. The impact of adducts derived from the highly reactive bay region enantiomeric (+)- and (-)-anti-7,8-diol-9,10-epoxide derivatives of benzo[a]pyrene (BPDE) are of interest because the (+)-7R,8S,9S,10R-anti-BPDE enantiomer is highly tumorigenic in rodents, while the (-)-7S,8R,9R,10S-anti-BPDE enantiomer is not. Both (+)- and (-)-anti-BPDE bind covalently with DNA predominantly by trans addition at the exocyclic amino group of guanine to yield 10S (+)- and 10R (-)-trans-anti-[BP]-N(2)-dG adducts. We have synthesized a number of different oligonucleotides with single (+)- and (-)-trans-anti-[BP]-N(2)-dG adducts (G) in the base sequence context XG*Y, where X and Y are different DNA bases. The G* residues were positioned at or close to the center of 11 base pair ( approximately 1 helical turn) or 16 base pair ( approximately 1.5 turns) duplexes. All bases, except for X and Y and their partners, were identical. These sequences were self-ligated with T4 ligase to form multimers that yield a ladder of bands upon electrophoresis in native polyacrylamide gels. The extent of bending in each oligonucleotide was assessed by monitoring the decrease in gel mobilities of these linear, self-ligated oligomers, relative to unmodified oligonucleotides of the same base sequence. The extent of global bending was then estimated using a sequence-specific three-dimensional model from which the values of the base-pair step parameter roll adjacent to the lesion site could be extracted. We find that (+)-trans-anti-[BP]-N(2)-dG adducts are considerably more bent than the (-) isomers regardless of sequence and that A-T base pairs flanking the [BP]-N(2)-dG lesion site allow for local flexibility consistent with adduct conformational heterogeneity. Interestingly, the fit of computed versus observed gel mobilities using classical reptation treatments requires enhancement of unmodified DNA flexibility in gels, compared to aqueous salt solution. The differences in bending between the two stereoisomeric adduct duplexes and the observed base sequence context effects may play a significant role in the differential processing of these lesions by cellular replication, transcription, and repair enzymes.

Slow motion in the CAA*TTG sequence of a DNA decamer duplex studied by NMR

In DNA duplexes, pyrimidine-purine steps are believed to be flexible or conformationally unstable. Indeed, several DNA crystal structures exhibit a multitude of conformations for CpA*TpG steps. The question arises of whether this structural flexibility is accompanied by dynamical flexibility, i.e., a question pertaining to the energy barrier between conformations. Except for TpA steps, slow motions on the microsecond-to-millisecond time scale have not been detected in duplexes until now. In the present study, such slow motion was investigated by 1H, 13C, and 15N NMR relaxation measurements on a DNA decamer d(CATTTGCATC)*d(GATGCAAATG). The DNA decamer was enriched with 15% 13C and 98% 15N isotopes for each adenosine and guanosine residue. Three lines of evidence support the notion of slow motion in the CAA*TTG moiety. Analysis of (15)N relaxation showed that the order parameter, S2, of guanosine imino NH groups was about 0.8, similar to that of CH groups for this oligomer. The strong temperature dependence of guanosine NH S2 in the CAA*TTG sequence indicated the presence of a large-amplitude motion. Signals of adenosine H8 protons in the CAA*TTG sequence were broadened in 2D 1H NOESY spectra, which also suggested the existence of slow motion. As well as being smaller than for other adenine residues, the 1H T2 values exhibited a magnetic field strength dependence for all adenosine H8 signals in the ATTTG*CAAAT region, suggesting slow motions more pronounced at the first adenosine in the CAA*TTG sequence but extending over the CAAAT*ATTTG region. This phenomenon was further examined by the pulse field strength dependence of the 1H, 13C, and 15N T1rho values. 1H and 13C T1rho values showed a pulse field strength dependence, but 15N T1rho did not. Assuming a two-site exchange process, an exchange time constant of 20-300 micros was estimated for the first adenosine in the CAA sequence. The exact nature of this motion remains unknown.

Conformation and rigidity of DNA microcircles containing waf1 response element for p53 regulatory protein

The tumor-suppressor activity of p53 is closely related to its DNA-binding properties. It binds a number of DNA response-elements and it is likely that these share a common structural feature. Here, we present a new, general method to determine the absolute twist of flexible DNA promoter sequences based on direct imaging of the topology of microcircles containing the sequences. We have used magnetically driven dynamic force microscopy ("MacMode" AFM) to observe, in solution, the conformation of 168 base-pair DNA microcircles, each containing four equally spaced copies of the waf1/cip1/p21 p53 response-element. Analysis of the images showed that the microcircles are markedly puckered with a small excess of negatively writhed molecules. The average measured values of writhe are 0.109+/-0.013 (for 60 positively writhed molecules) and -0.098+/-0.011 (for 65 negatively writhed molecules). These values lead directly to a difference in linking number for the positively and negatively writhed molecules prior to ligation, from which we derive a twist mismatch of 178 degrees (overtwist). This is 44.5 degrees for each 42-mer precursor containing a single waf1/cip1/p21 p53 response-element, in good agreement with the range of values deduced by indirect biochemical techniques. The two values of writhe may also be used to determine the ratio of the bending (B) to twisting (C) rigidity, yielding B/C=0.23. This is about one-third of the value for long, random-sequence DNA, suggesting that the waf1/cip1/p21 p53 response-element is extremely flexible, a result that is also consistent with indirect biochemical experiments. These results support the idea, proposed by us earlier, that torsional stress may play a role in the regulation of p53 binding through modulation of twist at the binding site.

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.

Effect of neighboring bases on base-pair stacking orientation: a molecular dynamics study

It is generally believed that base-pair stacking interaction in DNA double helix is one of the strongest interactions that governs sequence directed structural variability. However, X-ray crystal structures of some base-paired doublet sequences have been seen to adopt different structures when flanked by different base-pairs. DNA crystal database, however, is still too small to make good statistical inference about effect of such flanking residues. Influence of neighboring residue on the local helical geometry of a base-paired doublet in B-DNA has been investigated here using molecular dynamics simulation. We have generated ensembles of structures for d(CA).d(TG) and d(AA).d(TT) base-paired doublets located at the centers of d(CGCGCAAAGCG).d(CGCTTTGCGCG) and d(CGCGAAAACGCG).d(CGCGTTTTCGCG) sequences along with their analogs by varying the bases either at 5'- or 3'- position to the central doublet. Comparison of base paired doublet parameters for the ensembles of structures show that stacking geometry of d(CA).d(TG) doublet depends on some of the flanking base-pairs. On the other hand d(AA).d(TT) doublet remains nearly unperturbed when the flanking residues are altered.

A frequent TG deletion near the polyadenylation signal of the human HEXB gene: occurrence of an irregular DNA structure and conserved nucleotide sequence motif in the 3' untranslated region

While screening for new mutations in the HEXB gene, which encodes the beta-subunit of beta-hexosaminidase, a TG deletion (deltaTG) was found in the 3' untranslated region (3'UTR) of the gene, 7 bp upstream from the polyadenylation signal. Examination of DNA samples of 145 unrelated Argentinean individuals from different racial backgrounds showed that the deltaTG allele was present with a frequency of approximately 0.1, compared with the wild-type (WT) allele. The deletion was not associated with infantile or variant forms of Sandhoff disease when present in combination with a deleterious allele. Total Hex and Hex B enzymatic activities measured in individuals heterozygous for deltaTG and a null allele, IVS-2 + 1G-->A (G-->A), were approximately 30% lower than the activities of G-->A/WT individuals. Analysis of the HEXB mRNA from leukocytes of deltaTG/WT individuals by RT-PCR of the 3'UTR showed that the deltaTG allele is present at lower level than the WT allele. By polyacrylamide gel electrophoresis, it was determined that a PCR fragment containing the +TG version of the 3'UTR of the HEXB gene had an irregular structure. On inspection of genes containing a TG dinucleotide upstream from the polyadenylation signal we found that this dinucleotide was part of a conserved sequence (TGTTTT) immersed in a A/T-rich region. This sequence arrangement was present in more than 40% analyzed eukaryotic mRNAs, including in the human, mouse and cat HEXB genes. The significance of the TG deletion in reference to Sandhoff disease as well as the possible functional role of the consensus sequence and the DNA structure of the 3'UTR are considered.

Sequence directed flexibility of DNA and the role of cross-strand hydrogen bonds

Persistence length and torsional rigidity for different B-DNA sequences have been calculated by analysing crystal structure database. The values of these parameters for mixed sequence DNA are in good agreement with those estimated by others. Persistence lengths for the homopolymeric sequences, namely poly(dA).poly(dT) and poly(dG).poly(dC), are significantly large compared to those of others as expected from the inability of these sequences to form nucleosome under normal conditions. The heteropolymeric sequences poly(dA-dC).poly(dG-dT) and poly(dG-dC).poly(dG-dC), on the other hand, have smaller persistence lengths. This implies larger flexibility of the d(AC).d(GT), d(CA).d(TG), d(GC).d(GC) and d(CG).d(CG) doublets, some of which constitute the genetic disease forming triplet repeats d(CTG).d(CAG) and d(CGG).d(CCG). Thus it is expected that these triplet repeat sequences are also flexible and wrap around the histone octamer efficiently. Persistence length calculations also indicate larger flexibility for these triplet repeat sequences. Furthermore, our computations reveal that the rigidity of a given DNA sequence is controlled by its ability to form cross-strand bifurcated hydrogen bonds between the successive base pairs. Molecular orbital calculations suggest that these hydrogen bonds are generally extended with bond lengths around 3A.

A- and T-tract-mediated intrinsic curvature in native DNA between the binding site of the upstream activator NtrC and the nifLA promoter of Klebsiella pneumoniae facilitates transcription

The nif promoters of Klebsiella pneumoniae must be activated by proteins bound to upstream sequences which are thought to interact with the sigma54-RNA polymerase holoenzyme by DNA looping. NifA is the activator for most of the promoters, and integration host factor (IHF) mediates the DNA looping. While NtrC is the activator for the nifLA promoter, no IHF appears to be involved. There are two A tracts and one T tract between the upstream enhancer and the nifLA promoter. This DNA segment exhibits anomalous electrophoretic mobility, suggesting intrinsic sequence-induced curvature in the DNA. On the one hand, mutation of the A tracts or T tract individually or together, or deletion of the A tracts and the T tract reduces the anomaly; on the other hand, creation of two additional A tracts enhances the anomaly. Intrinsic curvature in the DNA has been confirmed by circular permutation analysis after cloning the DNA fragment in the vector pBend 2 and also by electron microscopy. Computer simulation with the DNA base sequence is also suggestive of intrinsic curvature. A transcriptional fusion with the Escherichia coli lacZ gene of the DNA fragment containing the nifLA promoter and the wild-type or the mutated upstream sequences was constructed, and in vivo transcription in K. pneumoniae and E. coli was monitored. There was indeed very good correlation between the extent of intrinsic curvature of the DNA and transcription from the promoter, suggesting that DNA curvature due to the A tracts and the T tract was necessary for transcription in vivo from the nifLA promoter of K. pneumoniae.

Cytosine methylation in a CpG sequence leads to enhanced reactivity with Benzo[a]pyrene diol epoxide that correlates with a conformational change

Benzo[a]pyrene (B[a]P) is a widespread environmental carcinogen that must be activated by cellular metabolism to a diol epoxide form (BPDE) before it reacts with DNA. It has recently been shown that BPDE preferentially modifies the guanine in methylated 5'-CpG-3' sequences in the human p53 gene, providing one explanation for why these sites are mutational hot spots. Using purified duplex oligonucleotides containing identical methylated and unmethylated CpG sequences, we show here that BPDE preferentially modified the guanine in hemimethylated or fully methylated CpG sequences, producing between 3- and 8-fold more modification at this site. Analysis of this reaction using shorter duplex oligonucleotides indicated that it was the level of the (+)-trans isomer that was specifically increased. To determine if there were conformational differences between the methylated and unmethylated B[a]P-modified DNA sequences that may be responsible for this enhanced reactivity, a native polyacrylamide gel electrophoresis analysis was carried out using DNA containing isomerically pure B[a]P-DNA adducts. These experiments showed that each adduct resulted in an altered gel mobility in duplex DNA but that only the presence of a (+)-trans isomer and a methylated C 5' to the adduct resulted in a significant gel mobility shift compared with the unmethylated case.

p53-induced DNA bending and twisting: p53 tetramer binds on the outer side of a DNA loop and increases DNA twisting

DNA binding activity of p53 is crucial for its tumor suppressor function. Our recent studies have shown that four molecules of the DNA binding domain of human p53 (p53DBD) bind the response elements with high cooperativity and bend the DNA. By using A-tract phasing experiments, we find significant differences between the bending and twisting of DNA by p53DBD and by full-length human wild-type (wt) p53. Our data show that four subunits of p53DBD bend the DNA by 32-36 degrees, whereas wt p53 bends it by 51-57 degrees. The directionality of bending is consistent with major groove bends at the two pentamer junctions in the consensus DNA response element. More sophisticated phasing analyses also demonstrate that p53DBD and wt p53 overtwist the DNA response element by approximately 35 degrees and approximately 70 degrees, respectively. These results are in accord with molecular modeling studies of the tetrameric complex. Within the constraints imposed by the protein subunits, the DNA can assume a range of conformations resulting from correlated changes in bend and twist angles such that the p53-DNA tetrameric complex is stabilized by DNA overtwisting and bending toward the major groove at the CATG tetramers. This bending is consistent with the inherent sequence-dependent anisotropy of the duplex. Overall, the four p53 moieties are placed laterally in a staggered array on the external side of the DNA loop and have numerous interprotein interactions that increase the stability and cooperativity of binding. The novel architecture of the p53 tetrameric complex has important functional implications including possible p53 interactions with chromatin.

Nuclear magnetic resonance structure of d(GCATATGATAG). d(CTATCATATGC): a consensus sequence for promoters recognized by sigma K RNA polymerase

The three-dimensional structure of d(GCATATGATAG).d(CTATCATATGC), from the promoter region of a gene regulating sporulation in Bacillus subtilis mother cells, was determined utilizing two-dimensional nuclear Overhauser effect (2D NOE) and double-quantum-filtered COSY (2QF-COSY) spectra. To minimize the effect of methods used to obtain restraints and refine structure, several variables were studied. Interproton distance bounds were calculated very conservatively by running the complete relaxation matrix program MARDIGRAS hundreds of times using 2D NOE spectra for exchangeable and for nonexchangeable protons at different mixing times, assuming different overall correlation times and different starting structures. The 435 distance restraints were used with two different structural refinement methods: restrained molecular dynamics (rMD) and restrained Monte Carlo calculations (rMC). Refinement using different procedures and starting structures resulted in essentially the same structure (<0.8 A rmsd), indicating that the structure is defined by experimental restraints and not the refinement method or variables used. R factors indicate the structures fit the experimental NOE data very well. Some helical parameters, notably large negative X displacement, are characteristic of A-DNA, but others are characteristic of B-DNA. As with TG.CA steps in other duplex DNA sequences studied in our laboratory, the two TG.CA steps have a positive roll, with T6-G7 exhibiting the largest, and consequently a bent helix axis. The converged structure represents a time-averaged structure. However, multiple conformations, especially in deoxyriboses, were evident from vicinal coupling constants obtained from quantitative simulations of 2QF-COSY cross-peaks and from persistent inconsistencies in experimental distances due to nonlinear conformational averaging.

Small circles of helical DNA obtained on the basis of transcriptional pause sites sequence

The 21-base pair synthetic DNA duplexes with basic 'pause-motif site ('CATGC') were ligated head-to-tail to produce linear and circular multimers. This also was done from other closely related sequences. Electrophoretic mobilities of the linear multimers in polyacrylamide gels were determined under the standard and modified conditions. We revealed that small linear multimers (approximately 90 bp) were characterized by comparable value of gel retardation relative to the well known curved DNA, while longer multimers (130 to approximately 170 bp) had only slightly expressed mobility anomaly. Nevertheless these multimers containing nontruncated 'pause-motif were capable of cyclization, in particular, formation of unusually small circles while truncated ones were not. We conclude that basic 'pause-motif site increases the closure ability while the multimers based on truncated 'pause motif fail to curve into the small circles. We tend to explain this situation as a result of intrinsic bending as well as the influence of the thermal fluctuations of DNA, the latter most probably can be associated with 'pause motif'. We have estimated the equilibrial and maximal bend angles per 10.5 bp to be 12 degrees to approximately 16 degrees and 32 degrees accordingly under experimental conditions of our study.

DNA binding discrimination of the murine DNA cytosine-C5 methyltransferase

Mammalian DNA cytosine-C5 methyltransferase modifies the CpG dinucleotide in the context of many different genomic sequences. A rigorous DNA binding assay was developed for the murine enzyme and used to define how sequences flanking the CpG dinucleotide affect the stability of the enzyme:DNA complex. Oligonucleotides containing a single CpG site form reversible 1:1 complexes with the enzyme that are sequence-specific. A guanine/cytosine-rich 30 base-pair sequence, a mimic of the GC-box cis-element, bound threefold more tightly than an adenine/thymine-rich sequence, a mimic of the cyclic AMP responsive element. However, the binding discrimination between hemi- and unmethylated forms of these DNA substrates was small, as we previously observed at the K(m)DNA level (Biochemistry, 35, 7308-7315 (1996)). Single-stranded substrates are bound much more weakly than double-stranded DNA forms. An in vitro screening method was used to select for CpG flanking sequence preferences of the DNA methyltransferase from a large, divergent population of DNA substrates. After five iterative rounds of increasing selective pressure, guanosine/cytosine-rich sequences were abundant and contributed to binding stabilization for at least 12 base-pairs on either side of a central CpG. Our results suggest a read-out of sequence-dependent conformational features, such as helical flexibility, minor groove dimensions and critical phosphate orientation and mobility, rather than interactions with specific bases over the course of two complete helical turns. Thus, both studies reveal a preference for guanosine/cytosine deoxynucleotides flanking the cognate CpG. The enzyme specificity for similar sequences in the genome may contribute to the in vivo functions of this vital enzyme.

Sequence determinants of the intrinsic bend in the cyclic AMP response element

The cyclic AMP response element (CRE site, ATGACGTCAT) is the DNA target for transcription factors whose activities are regulated by cyclic AMP (1). Recently, we discovered that the CRE site is bent by 10-13 degrees toward the major groove (2). Little or no bend is detected in the related AP-1 site (ATGACTCAT), which differs from the CRE site by loss of a single, central, C.G base pair (2, 3). Here we describe experiments designed to identify which base pairs within the CRE site induce the bent structure in an attempt to understand the origins of the dramatically different conformations of the CRE and AP-1 sites. Our data indicate that the intrinsic CRE bend results from distortion within the TGA sequence found in each CRE half site (ATGAC). These two TGA sequences are located in phase with one another in the CRE sequence but are not (completely) in phase in the AP-1 sequence. This difference in phasing leads to the overall difference in bend as detected by gel (2) and cyclization methods (S. C. Hockings, J. D. Kahn, and D. M. Crothers, unpublished results; M. A. Fabian and A. Schepartz, unpublished results). Our results confirm earlier predictions of altered structure within TG steps, provide insight into the structural reorganizations induced in DNA by bZIP proteins, and lead to a revision of the relationship between the structures of the free and bZIP-bound forms of the CRE and AP-1 sites.

Sequence dependence and characteristics of bends induced by site-specific polynuclear aromatic carcinogen-deoxyguanosine lesions in oligonucleotides

The tumorigenic metabolite of benzo[a]pyrene, the (+)-7R,8S,9S,10R enantiomer, and the nontumorigenic mirror-image isomer, (-)-7S,8R,9R, 10S, of r7,t8-dihydroxy-t9,10-epoxy-7,8,9, 10-tetrahydrobenzo[a]pyrene (anti-BPDE) bind covalently to the exocyclic amino group of deoxyguanosine (N2-dG) in native DNA. These adducts can cause structural perturbations such as DNA bends, which in turn may influence the cellular processing of these lesions. The characteristics of bends in site-specifically modified oligodeoxyribonucleotide duplexes induced by single (+)- and (-)-anti-[BP]-N2-dG lesions were examined by self-ligation and gel electrophoresis techniques. The modified residues (dG*) were centrally positioned in the 11-mer oligonucleotide d(CACAXG*XACAC) complexed with the natural complementary strands, with X = T or C, or in oligonucleotides 16 or 22 base pairs long with the same centrally positioned 11-mer. Among the four stereochemically distinct lesions, the 10S(+)-trans-anti-[BP]-N2-dG adducts were significantly more bent than any of the other three stereoisomeric adducts and were selected for detailed studies. In the TG*T sequence context (X = T), the retardation factor RL (apparent length of multimer/sequence length) is approximately independent of the phasing (distance, in base pairs, between the lesions) of the adducts with respect to the helical repeat (10.5 base pairs/helix turn). In contrast, in the CG*C sequence context (X = C), RL is markedly lower in the case of ligated 16-mers than in the case of ligated 11-mer duplexes. The dependence of RL on the phasing of the bends as a function of the helical repeat, indicate that the bends associated with (+)-trans-anti-[BP]-N2-dG lesions are relatively rigid in the d(...CG*C...).d(...GCG...) sequences, and flexible in the d(...TG*T...).d(...ACA...) sequence context. These differences are attributed to the orientations of the pyrenyl residues on the 5'-side of the modified deoxyguanosine residues in the minor groove and to the intrinsic roll and tilt characteristics of DNA dinucleotide steps CG, GC, TG, and GT. The influence of flanking bases on the extent and character of DNA bending suggest that base sequence effects may be important in the cellular processing of (+)-trans-anti-[BP]-N2-dG lesions.

The intrinsic curvature of a 51 bp K-DNA fragment of Leishmania tarentolae: a molecular model

DNA intrinsic structure and curvature is a subject of debate because of the importance of these attributes in processes such as DNA packaging, transcription, and gene regulation. X-ray crystallography of DNA single crystals has provided a wealth of information about the local, short range conformational features of DNA. On the other hand, gel electrophoresis analysis of DNA has not only uncovered the macroscopic curvature of DNA but it also provides most of the available data on DNA intrinsic curvature. However, gel electrophoresis can not identify features of DNA structure at the nucleotide or atomic level. In order to address the problem of DNA intrinsic curvature in an attempt to bridge the gap between X-ray crystallography and gel electrophoresis, we use the computational method of molecular dynamics (MD). In this study, we report the results of 2.0 ns MD simulations on a 51 bp fragment of the K-DNA of Leishmania tarentolae containing several A-tracts. The K-DNA double helix is very stable and remains in an intermediate state between the canonical A and B forms of the duplex. The magnitude of global curvature (75 degrees) agrees well with the experimental estimate (72 degrees) available. Analysis of local (every base triplet) and sublocal (every helix turn) curvature shows that the 51 bp K-DNA fragment has curvature features also present in the Wedge, Junction and Calladine's models of DNA intrinsic curvature. We further characterize the flexibility of individual nucleotides in the molecule and find the sugar flexibility within the A-tracts to be strongly correlated with the pattern of A-tract cleavage by the hydroxyl radical. Differential curvature and flexibility at the 5' and 3'junctions between A-tracts and general-sequence DNA are found to modulate the global curvature of the K-DNA fragment.

DNA bending is essential for the site-specific recognition of DNA response elements by the DNA binding domain of the tumor suppressor protein p53

We have used circular permutation assays to determine the extent and location of the DNA bend induced by the DNA binding domain of human wild type p53 (p53DBD) upon binding to several naturally occurring DNA response elements. We have found that p53DBD binding induces axial bending in all of the response elements investigated. In particular, response elements having a d(CATG) sequence at the junction of two consensus pentamers in each half-site favor highly bent complexes (bending angle is approximately 50 degrees ), whereas response elements having d(CTTG) bases at this position are less bent (bending angles from approximately 37 to approximately 25 degrees ). Quantitative electrophoretic mobility shift assays of different complexes show a direct correlation between the DNA bending angle and the binding affinity of the p53DBD with the response elements, i.e. the greater the stability of the complex, the more the DNA is bent by p53DBD binding. The study provides evidence that the energetics of DNA bending, as determined by the presence or absence of flexible sites in the response elements, may contribute significantly to the overall binding affinity of the p53DBD for different sequences. The results therefore suggest that both the structure and the stability of the p53-DNA complex may vary with different response elements. This variability may be correlated with variability in p53 function.

Architectural accommodation in the complex of four p53 DNA binding domain peptides with the p21/waf1/cip1 DNA response element

High resolution chemical footprinting and cross-linking experiments have provided a basis for elucidating the overall architecture of the complex between the core DNA binding domain of p53 (p53DBD, amino acids 98-309) and the p21/waf1/cip1 DNA response element implicated in the G1/S phase cell cycle checkpoint. These studies complement both a crystal structure and earlier biophysical studies and provide the first direct experimental evidence that four subunits of p53DBD bind to the response element in a regular staggered array having pseudodyad symmetry. The invariant guanosines in the highly conserved C(A/T)|(T/A)G parts of the consensus half-sites are critical to the p53DBD-DNA binding. Molecular modeling of the complex using the observed peptide-DNA contacts shows that when four subunits of p53DBD bind the response element, the DNA has to bend approximately 50 degrees to relieve steric clashes among different subunits, consistent with recent DNA cyclization studies. The overall lateral arrangement of the four p53 subunits with respect to the DNA loop comprises a novel nucleoprotein assembly that has not been reported previously in other complexes. We suggest that this kind of nucleoprotein superstructure may be important for p53 binding to response elements packed in chromatin and for subsequent transactivation of p53-mediated genes.

DNA mobility anomalies are determined primarily by polyacrylamide gel concentration, not gel pore size

The dependence of DNA mobility anomalies on gel pore size has been studied in polyacrylamide gels with a wide variety of compositions, using molecular weight ladders containing multiple copies of normal (12B) and anomalously slowly migrating (12A) 147-base pair restriction fragments from plasmid pBR322 as the migrating probe molecules. If the gel pore size is increased by decreasing the total acrylamide concentration (%T) at constant cross-linker ratio (%C), the usual method of increasing gel pore size, the mobility anomalies decrease with increasing gel pore radius as though the 12A multimers were retarded by a sieving mechanism. However, the decrease in the mobility anomalies is independent of whether the apparent gel pore radius is larger or smaller than the DNA radius of gyration, suggesting that gel pore size is not the controlling variable. If the acrylamide concentration is held constant and the gel pore size is increased by decreasing %C at constant %T, the mobility anomalies of the largest 12A multimers (6 mers and higher) decrease with increasing gel pore radius, because of sieving effects, until the effective gel pore radius becomes approximately equal to the DNA radius of gyration, after which the mobility anomalies level off and become independent of gel pore size. The mobility anomalies exhibited by 5-mers and smaller multimers of fragment 12A are independent of gel pore radius in all gels with constant %T. Similar results are observed with a molecular weight ladder containing phased A-tracts from the kinetoplast bending locus. Since the anomalous electrophoretic mobilities depend primarily on the total acrylamide concentration in the gel, and not on the apparent gel pore radius, increases in the magnitude of the mobility anomalies with increasing gel concentration (and decreasing gel pore radius) cannot be taken as evidence for DNA curvature.

Cruciform structures and functions

In this paper, a structure-function analysis of B-DNA self-fitting is reviewed in the light of recent oligonucleotide crystal structures. Their crystal packings provided a high-resolution view of B-DNA helices closely and specifically fitted by groove-backbone interaction, a natural and biologically relevant manner to assemble B-DNA helices. In revealing that new properties of the DNA molecule emerge during condensation, these crystallographic studies have pointed to the biological importance of DNA—DNA interactions.

DNA structural forms

In the years that have passed since the publication of Wolfram Saenger's classic book on nucleic acid structure (Saenger, 1984), a considerable amount of new data has been accumulated on the range of conformations which can be adopted by DNA. Many unusual species have joined the DNA zoo, including new varieties of two, three and four stranded helices. Much has been learnt about intrinsic DNA curvature, dynamics and conformational transitions and many types of damaged or deformed DNA have been investigated. In this article, we will try to summarise this progress, pointing out the scope of the various experimental techniques used to study DNA structure, and, where possible, trying to discern the rules which govern the behaviour of this subtle macromolecule. The article is divided into six major sections which begin with a general discussion of DNA structure and then present successively, B-DNA, DNA deformations, A-DNA, Z-DNA and DNARNA hybrids. An extensive set of references is included and should serve the reader who wishes to delve into greater detai.