DNA curving and bending in protein-DNA recognition

From straight to curved: A historical perspective of DNA shape

DNA is the macromolecule responsible for storing the genetic information of a cell and it has intrinsic properties such as deformability, stability and curvature. DNA Curvature plays an important role in gene transcription and, consequently, in the subsequent production of proteins, a fundamental process of cells. With recent advances in bioinformatics and theoretical biology, it became possible to analyze and understand the involvement of DNA Curvature as a discriminatory characteristic of gene-promoting regions. These regions act as sites where RNAp (ribonucleic acid-polymerase) binds to initiate transcription. This review aims to describe the formation of Curvature, as well as highlight its importance in predicting promoters. Furthermore, this article provides the potential of DNA Curvature as a distinguishing feature for promoter prediction tools, as well as outlining the calculation procedures that have been described by other researchers. This work may support further studies directed towards the enhancement of promoter prediction software.

COL1A1 proximal promoter topology regulates its transcriptional response to transforming growth factor β

OBJECTIVE

The COL1A1 proximal promoter contains two GC-rich regions and two inverted CCAAT boxes. The transcription factors Sp1 and CBF bind to the GC sequence at -122 to -115 bp and the inverted CCAAT box at -101 to -96 bp, respectively, and stimulate COL1A1 transcriptional activity.

METHODS

To further define the regulatory mechanisms controlling COL1A1 expression by Sp1 and CBF, we introduced 2, 4, 6, or 8 thymidine nucleotides (T-tracts) at position -111 bp of the COL1A1 gene promoter to increase the physical distance between these two binding sites and examined in vitro the transcriptional activities of the resulting constructs and their response to TGF-β1.`.

RESULTS

Insertion of 2 or 4 nucleotides decreased COL1A1 promoter activity by up to 70%. Furthermore, the expected increase in COL1A1 transcription in response to TGF-β1 was abolished. Computer modeling of the modified DNA structure indicated that increasing the physical distance between the Sp1 and CBF binding sites introduces a rotational change in the DNA topology that disrupts the alignment of Sp1 and CBF binding sites and likely alters protein-protein interactions among these transcription factors or their associated co-activators.

CONCLUSION

The topology of the COL1A1 proximal promoter is crucial in determining the transcriptional activity of the gene and its response to the stimulatory effects of TGF-β1.

C5' omitted DNA enhances bendability and protein binding

Protein-DNA interactions are of great biological importance. The specificity and strength of these intimate contacts are crucial in the proper functioning of a cell, wherein the role of DNA dynamic bendability has been a matter of discussion. We relate DNA bendability to protein binding by introducing some simple modifications in the DNA structure. We removed C5' carbon in first modified structure and the second has an additional carbon between C3' and 3'-OH, hereby pronounced as C(-) and C(+) nucleic acids respectively. We observed that C(+) nucleic acid retains B-DNA duplex as seen by means of 500 ns long molecular dynamics (MD) simulations, structural and energetic calculations, while C(-) nucleic acid attains a highly bend structure. We transferred these observations to a protein-DNA system in order to monitor as to what extent the bendability enhances the protein binding. The energetics of binding is explored by performing 100 ns long MD simulations on control and modified DNA-protein complexes followed by running MM-PBSA/GBSA calculations on the resultant structures. It is observed that C(+) nucleic acid has protein binding in close correspondence to the control system (∼-14 kcal/mol) due to their relatable structure, while the C(-) nucleic acid displayed high binding to the protein (∼-18 kcal/mol). DelPhi based calculations reveal that the high binding could be the result of enhanced electrostatic interactions caused by exposed bases in the bend structure for protein recognition. Such modified oligonucleotides, due to their improved binding to protein and resistance to nuclease degradation, have a great therapeutic value.

DNA polymorphism in crystals: three stable conformations for the decadeoxynucleotide d(GCATGCATGC)

High-resolution structures of DNA fragments determined using X-ray crystallography or NMR have provided descriptions of a veritable alphabet of conformations. They have also shown that DNA is a flexible molecule, with some sequences capable of adopting two different structures. Here, the first example is presented of a DNA fragment that can assume three different and distinct conformations in crystals. The decanucleotide d(GCATGCATGC) was previously reported to assume a single-stranded double-fold structure. In one of the two crystal structures described here the decamer assumes both the double-fold conformation and, simultaneously, the more conventional B-type double-helical structure. In the other crystal the sequence assumes the A-type double-helical conformation. These results, taken together with CD spectra, which were recorded as the decamer was titrated against four metal ions and spermine, indicate that the molecule may exist as a mixed population of structures in solution. Small differences in the environmental conditions, such as the concentration of metal ion, may decide which of these crystallizes out. The results also support the idea that it may be possible for DNA to change its structure to suit the binding requirements of proteins or drugs.

DNA and RNA folds in transcription complex as evidenced by iodine-125 radioprobing

Abstract Folding of the DNA and RNA strands in an arrested T7 RNA polymerase (RNAP) transcription complex was studied by radioprobing, a novel method based on an analysis of the strand breaks produced by decay of the iodine-125 incorporated at the C5 position of cytosine. (125)I-labeled cytosines were incorporated into transcripts at different positions relative to the site of the arrest. The intensities of the DNA breaks inversely correlate with the distances from the (125)I decay site, so the radioprobing data provide information about the spatial RNA/DNA folding during transcription. We found that the yield of DNA strand breaks is significantly higher in the template than the non-template strand. This is consistent with local opening of the DNA duplex and formation of a hybrid between RNA and the template DNA strand. Our data demonstrate that the RNA-DNA hybrid has a nonuni form A-like structure. When the (125)I is incorporatedseven nucleotides from the active center of RNAP, the yield of strand breaks is substantially lower than if (125)I is positioned at the ends of the hybrid. Consequently, the DNA and RNA strands are located closer to each other at the ends of the hybrid and somehow separated in the middle. Surprisingly, the (125)I-induced breaks were detected in both DNA strands upstream from the transcription "bubble" indicating that DNA and RNA are closely associated outside the RNAP cleft. Thus, radioprobing data imply that the RNA/DNA fold in the complex with T7 RNAP is more complicated than had been anticipated by the existing models. Based on the present data, we suggest a sterically feasible model explaining how formation of the long RNA-DNA hybrid can result in the initiation-to-elongation switch in the T7 transcription complex. According to this model, the topological linkage between the DNA and RNA strands provides the necessary stability for the elongation complex, while permitting movement of the polymerase along the DNA duplex.

Oxygen-aromatic contacts in intra-strand base pairs: analysis of high-resolution DNA crystal structures and quantum chemical calculations

Three-dimensional structures of biomolecules are stabilized by a large number of non-covalent interactions and some of them such as van der Waals, electrostatic and hydrogen bond interactions are well characterized. Delocalized π-electron clouds of aromatic residues are known to be involved in cation-π, CH-π, OH-π and π-π interactions. In proteins, many examples have been found in which the backbone carbonyl oxygen of one residue makes close contact with the aromatic center of aromatic residues. Quantum chemical calculations suggest that such contacts may provide stability to the protein secondary structures. In this study, we have systematically analyzed the experimentally determined high-resolution DNA crystal structures and identified 91 examples in which the aromatic center of one base is in close contact (<3.5Ǻ) with the oxygen atom of preceding (Group-I) or succeeding base (Group-II). Examples from Group-I are overwhelmingly observed and cytosine or thymine is the preferred base contributing oxygen atom in Group-I base pairs. A similar analysis of high-resolution RNA structures surprisingly did not yield many examples of oxygen-aromatic contact of similar type between bases. Ab initio quantum chemical calculations on compounds based on DNA crystal structures and model compounds show that interactions between the bases in base pairs with oxygen-aromatic contacts are energetically favorable. Decomposition of interaction energies indicates that dispersion forces are the major cause for energetically stable interaction in these base pairs. We speculate that oxygen-aromatic contacts in intra-strand base pairs in a DNA structure may have biological significance.

A novel indirect sequence readout component in the E. coli cyclic AMP receptor protein operator

The cyclic AMP receptor protein (CRP) from Escherichia coli has been extensively studied for several decades. In particular, a detailed characterization of CRP interaction with DNA has been obtained. The CRP dimer recognizes a consensus sequence AANTGTGANNNNNNTCACANTT through direct amino acid nucleobase interactions in the major groove of the two operator half-sites. Crystal structure analyses have revealed that the interaction results in two strong kinks at the TG/CA steps closest to the 6-base-pair spacer (N6). This spacer exhibits high sequence variability among the more than 100 natural binding sites in the E. coli genome, but the exact role of the N6 region in CRP interaction has not previously been systematic examined. Here we employ an in vitro selection system based on a randomized N6 spacer region to demonstrate that CRP binding to the lacP1 site may be enhanced up to 14-fold or abolished by varying the N6 spacer sequences. Furthermore, on the basis of sequence analysis and uranyl (UO2(2+)) probing data, we propose that the underlying mechanism relies on N6 deformability.

Recognition of single- and double-stranded oligonucleotides by bovine serum albumin via nonspecific interactions

The complexes formed by bovine serum albumin (BSA) with single-stranded oligonucleotide (ss-oligo) or double-stranded oligonucleotide (ds-oligo) were investigated by laser light scattering, zeta potential analysis, and atomic force microscopy. It was found that BSA was able to recognize ss-oligo and ds-oligo upon forming complexes in HCOOH-HCOONa buffer at pH 3.0. When oligonucleotide was added dropwise to BSA, BSA formed a complex with ss-oligo but not with ds-oligo in the studied charge ratio. When BSA was added to oligonucleotides, BSA formed complexes with both ss-oligo and ds-oligo but via different paths: the BSA/ds-oligo underwent two processes, heavy precipitation followed by reentry, with increasing BSA/oligo charge ratio, whereas BSA/ss-oligo underwent only aggregation process, but with a charge reversal occurred at BSA/oligo charge ratio about 0.1. Moreover, the complex formed by BSA and ds-oligo showed a kinetics much slower than that of BSA and ss-oligo. We attributed the big difference upon complexation to the physical nature of oligonucleotides as well as the conformational change of BSA under severe conditions. The differentiation of ss-oligo from ds-oligo by BSA via nonspecific interactions gained insight in the recognition of DNA or RNA by specific protein (enzyme) under physiological conditions.

Sequential XylS-CTD binding to the Pm promoter induces DNA bending prior to activation

XylS protein, a member of the AraC family of transcriptional regulators, comprises a C-terminal domain (CTD) involved in DNA binding and an N-terminal domain required for effector binding and protein dimerization. In the absence of benzoate effectors, the N-terminal domain behaves as an intramolecular repressor of the DNA binding domain. To date, the poor solubility properties of the full-length protein have restricted XylS analysis to genetic approaches in vivo. To characterize the molecular consequences of XylS binding to its operator, we used a recombinant XylS-CTD variant devoid of the N-terminal domain. The resulting protein was soluble and monomeric in solution and activated transcription from its cognate promoter in an effector-independent manner. XylS binding sites in the Pm promoter present an intrinsic curvature of 35 degrees centered at position -42 within the proximal site. Gel retardation and DNase footprint analysis showed XylS-CTD binding to Pm occurred sequentially: first a XylS-CTD monomer binds to the proximal site overlapping the RNA polymerase binding sequence to form complex I. This first event increased Pm bending to 50 degrees and was followed by the binding of the second monomer, which further increased the observed global curvature to 98 degrees. This generated a concomitant shift in the bending center to a region centered at position -51 when the two sites were occupied (complex II). We propose a model in which DNA structure and binding sequences strongly influence XylS binding events previous to transcription activation.

An FTIR investigation of flanking sequence effects on the structure and flexibility of DNA binding sites

Fourier transform infrared (FTIR) spectroscopy and a library of FTIR marker bands have been used to examine the structure and relative flexibilities conferred by different flanking sequences on the EcoRI binding site. This approach allowed us to examine unique peaks and subtle changes in the spectra of d(AAAGAATTCTTT)(2), d(TTCGAATTCGAA)(2), and d(CGCGAATTCGCG)(2) and thereby identify local changes in base pairing, base stacking, backbone conformation, glycosidic bond rotation, and sugar puckering in the studied sequences. The changes in flanking sequences induce differences in the sugar puckers, glycosidic bond rotation, and backbone conformations. Varying levels of local flexibility are observed within the sequences in agreement with previous biological activity assays. The results also provide supporting evidence for the presence of a splay in the G(4)-C(9) base pair of the EcoRI binding site and a potential pocket of flexibility at the G(4) cleavage site that have been proposed in the literature. In sum, we have demonstrated that FTIR is a powerful methodology for studying the effect of flanking sequences on DNA structure and flexibility, for it can provide information about the local structure of the nucleic acid and the overall relative flexibilities conferred by different flanking sequences.

Upstream curved sequences in E. coli are related to the regulation of transcription initiation

The advancement in Escherichia coli genome research has made the information regarding transcription start sites of many genes available. A study relying on the availability of transcription start locations was performed. The first question addressed was what an average DNA curvature profile upstream of genes would look like when these genes are aligned by transcription start sites in comparison to alignment by translation start sites. Since it was hypothesized that curvature plays a role in transcription regulation, the expectation was that curvature measurements relative to transcription starts, rather than translation, should strengthen the signal. Our study justified this expectation. The second question aimed to clarify the relation between DNA curvature and promoter strength. Through clustering based on DNA curvature profiles along promoter regions, a strong positive correlation between the promoter strength and the curved DNA was found. The third question dealt with dinucleotide periodicity in E. coli to see whether a periodicity pattern specific to promoter regions exists. Such unknown pattern might shed new light on transcription regulation mechanisms in E. coli. A sequence periodicity of about 11 bp is characteristic to the whole E. coli genome, and is especially well-expressed in intergenic regions. Here it was shown that regions of the size of about 100-150 bp centered 70-100 bp upstream to transcription starts carry hidden periodicity with a period of about 10.3 bp.

Measuring spatial preferences at fine-scale resolution identifies known and novel cis-regulatory element candidates and functional motif-pair relationships

Transcriptional regulation is mediated by the collective binding of proteins called transcription factors to cis-regulatory elements. A handful of factors are known to function at particular distances from the transcription start site, although the extent to which this occurs is not well understood. Spatial dependencies can also exist between pairs of binding motifs, facilitating factor-pair interactions. We sought to determine to what extent spatial preferences measured at high-scale resolution could be utilized to predict cis-regulatory elements as well as motif-pairs binding interacting proteins. We introduce the ‘motif positional function’ model which predicts spatial biases using regression analysis, differentiating noise from true position-specific overrepresentation at single-nucleotide resolution. Our method predicts 48 consensus motifs exhibiting positional enrichment within human promoters, including fourteen motifs without known binding partners. We then extend the model to analyze distance preferences between pairs of motifs. We find that motif-pairs binding interacting factors often co-occur preferentially at multiple distances, with intervals between preferred distances often corresponding to the turn of the DNA double-helix. This offers a novel means by which to predict sequence elements with a collective role in gene regulation.

Parameterization of azole-bridged dinuclear platinum anticancer drugs via a QM/MM force matching procedure

Azole-bridged diplatinum compounds are promising new anticancer drugs designed to induce small distortions upon DNA alkylation, able to circumvent resistance problems of existing platinum drugs. Hybrid quantum classical (QM/MM) molecular dynamics (MD) simulations of different azole-bridged platinum drugs have recently revealed the nature of the local deformations at the DNA binding site. However, the description of global slow converging rearrangements cannot be addressed by QM/MM MD due to the short time scale accessible. Extensive classical MD simulations are therefore mandatory to describe accurately the structural distortions in the DNA double helix. This issue is now addressed by developing a new set of accurate force field parameters of the platinated moiety via a recently proposed force matching procedure of the classical forces to ab initio forces obtained from QM/MM trajectories. The accuracy of our force field parameters is validated by comparison of structural properties from classical MD and hybrid QM/MM simulations. The structural characteristics of the Pt-lesion are well reproduced during classical MD compared with QM/MM simulations and available experimental data. The global distortions in the DNA duplex upon binding of dinuclear Pt-compounds are very small and rather opposite to those induced by cisplatin. Thus, the force match approach significantly extends the potentialities of molecular simulations in the study of anticancer drugs and of the interactions with their biological targets.

Energetics of the protein-DNA-water interaction

Background

To understand the energetics of the interaction between protein and DNA we analyzed 39 crystallographically characterized complexes with the HINT (Hydropathic INTeractions) computational model. HINT is an empirical free energy force field based on solvent partitioning of small molecules between water and 1-octanol. Our previous studies on protein-ligand complexes demonstrated that free energy predictions were significantly improved by taking into account the energetic contribution of water molecules that form at least one hydrogen bond with each interacting species.

Results

An initial correlation between the calculated HINT scores and the experimentally determined binding free energies in the protein-DNA system exhibited a relatively poor r² of 0.21 and standard error of ± 1.71 kcal mol^-1. However, the inclusion of 261 waters that bridge protein and DNA improved the HINT score-free energy correlation to an r² of 0.56 and standard error of ± 1.28 kcal mol^-1. Analysis of the water role and energy contributions indicate that 46% of the bridging waters act as linkers between amino acids and nucleotide bases at the protein-DNA interface, while the remaining 54% are largely involved in screening unfavorable electrostatic contacts.

Conclusion

This study quantifies the key energetic role of bridging waters in protein-DNA associations. In addition, the relevant role of hydrophobic interactions and entropy in driving protein-DNA association is indicated by analyses of interaction character showing that, together, the favorable polar and unfavorable polar/hydrophobic-polar interactions (i.e., desolvation) mostly cancel.

A-tract clusters may facilitate DNA packaging in bacterial nucleoid

Molecular mechanisms of bacterial chromosome packaging are still unclear, as bacteria lack nucleosomes or other apparent basic elements of DNA compaction. Among the factors facilitating DNA condensation may be a propensity of the DNA molecule for folding due to its intrinsic curvature. As suggested previously, the sequence correlations in genome reflect such a propensity [Trifonov and Sussman (1980) Proc. Natl Acad. Sci. USA, 77, 3816-3820]. To further elaborate this concept, we analyzed positioning of A-tracts (the sequence motifs introducing the most pronounced DNA curvature) in the Escherichia coli genome. First, we observed that the A-tracts are over-represented and distributed 'quasi-regularly' throughout the genome, including both the coding and intergenic sequences. Second, there is a 10-12 bp periodicity in the A-tract positioning indicating that the A-tracts are phased with respect to the DNA helical repeat. Third, the phased A-tracts are organized in approximately 100 bp long clusters. The latter feature was revealed with the help of a novel approach based on the Fourier series expansion of the A-tract distance autocorrelation function. Since the A-tracts introduce local bends of the DNA duplex and these bends accumulate when properly phased, the observed clusters would facilitate DNA looping. Also, such clusters may serve as binding sites for the nucleoid-associated proteins that have affinities for curved DNA (such as HU, H-NS, Hfq and CbpA). Therefore, we suggest that the approximately 100 bp long clusters of the phased A-tracts constitute the 'structural code' for DNA compaction by providing the long-range intrinsic curvature and increasing stability of the DNA complexes with architectural proteins.

Influence of DNA stiffness in protein–DNA recognition

Protein-DNA recognition plays an essential role in the regulation of gene expression. The protein-DNA binding specificity is based on direct atomic contacts between protein and DNA and/or the conformational properties of DNA. In this work, we have analyzed the influence of DNA stiffness (E) to the specificity of protein-DNA complexes. The average DNA stiffness parameters for several protein-DNA complexes have been computed using the structure based sequence dependent stiffness scale. The relationship between DNA stiffness and experimental protein-DNA binding specificity has been brought out. We have investigated the importance of DNA stiffness with the aid of experimental free energy changes (DeltaDeltaG) due to binding in several protein-DNA complexes, such as, ETS proteins, 434, lambda, Mnt and trp repressors, 434 cro protein, EcoRV endonuclease V and zinc fingers. We found a correlation in the range 0.65-0.97 between DeltaDeltaG and E in these examples. Further, we have qualitatively analyzed the effect of mutations in the target sequence of lambda repressor and we observed that the DNA stiffness could correctly identify 70% of the correct bases among the considered nine positions.

Analysis of scanning force microscopy images of protein-induced DNA bending using simulations

Bending of DNA is a feature essential to the function of many DNA-binding proteins. Bending angles can be estimated with a variety of techniques, but most directly from images obtained using scanning force microscopy (SFM). Direct measurement of the bending angle using a tangent method often produces angles that deviate significantly from values obtained using other techniques. Here, we describe the application of SFM in combination with simulations of DNA as a means to estimate protein-induced bending angles in a reliable and unbiased fashion. In this manner, we were able to obtain accurate estimates for the bending angles induced by nuclear factor I, octamer-binding transcription factor 1, the human XPC-Rad23B complex and integration host factor [correction]

New elements of the termination of transcription in prokaryotes

We report here that phased runs of adenines and thymines are very frequent in the neighborhood of 3' of the coding regions of Escherichia coli and Bacillus subtilis. These findings suggest that the DNA curvature could affect transcription termination either directly, through contacts with RNA polymerase, or indirectly, via contacts with some regulatory proteins.

Thermodynamics of DNA binding and distortion by the hyperthermophile chromatin protein Sac7d

Sac7d is a hyperthermophile chromatin protein which binds non-specifically to the minor groove of duplex DNA and induces a sharp kink of 66 degrees with intercalation of valine and methionine side-chains. We have utilized the thermal stability of Sac7d and the lack of sequence specificity to define the thermodynamics of DNA binding over a wide temperature range. The binding affinity for poly(dGdC) was moderate at 25 degrees C (Ka = 3.5(+/-1.6) x 10(6) M(-1)) and increased by nearly an order of magnitude from 10 degrees C to 80 degrees C. The enthalpy of binding was unfavorable at 25 degrees C, and decreased linearly from 5 degrees C to 60 degrees C. A positive binding heat at 25 degrees C is attributed in part to the energy of distorting DNA, and ensures that the temperature of maximal binding affinity (75.1+/-5.6 degrees C) is near the growth temperature of Sulfolobus acidocaldarius. Truncation of the two intercalating residues to alanine led to a decreased ability to bend and unwind DNA at 25 degrees C with a small decrease in binding affinity. The energy gained from intercalation is slightly greater than the free energy penalty of bending duplex DNA. Surprisingly, reduced distortion from the double alanine substitution did not lead to a significant decrease in the heat of binding at 25 degrees C. In addition, an anomalous positive DeltaCp of binding was observed for the double alanine mutant protein which could not be explained by the change in polar and apolar accessible surface areas. Both the larger than expected binding enthalpy and the positive heat capacity can be explained by a temperature dependent structural transition in the protein-DNA complex with a Tm of 15-20 degrees C and a DeltaH of 15 kcal/mol. Data are discussed which indicate that the endothermic transition in the complex is consistent with DNA distortion.

Sequence periodicity of Escherichia coli is concentrated in intergenic regions

BACKGROUND

Sequence periodicity with a period close to the DNA helical repeat is a very basic genomic property. This genomic feature was demonstrated for many prokaryotic genomes. The Escherichia coli sequences display the period close to 11 base pairs.

RESULTS

Here we demonstrate that practically only ApA/TpT dinucleotides contribute to overall dinucleotide periodicity in Escherichia coli. The noncoding sequences reveal this periodicity much more prominently compared to protein-coding sequences. The sequence periodicity of ApC/GpT, ApT and GpC dinucleotides along the Escherichia coli K-12 is found to be located as well mainly within the intergenic regions.

CONCLUSIONS

The observed concentration of the dinucleotide sequence periodicity in the intergenic regions of E. coli suggests that the periodicity is a typical property of prokaryotic intergenic regions. We suppose that this preferential distribution of dinucleotide periodicity serves many biological functions; first of all, the regulation of transcription.

Atomic force microscopy reveals kinks in the p53 response element DNA

p53 is a 53 kDa nuclear phosphoprotein. Its function as a tumor suppressor critically lies in its ability to recognize its target DNA response elements as a tetramer. Here, we report the structural theme intrinsic to the response element DNA that governs this recognition phenomenon. The intrinsic flexibility or dynamic bending between two distinctly different, but naturally occurring p53 response elements has been compared by ring closure. Results show that DNA binding sites containing helically phased d(CATG.CATG) tetra-nucleotide sequences at the centers of quasi-dyad symmetry in each half-response site are more intrinsically flexible (i.e. preferentially bent under axial stress) than their d(CTTG.CTTG) counterparts. Intriguingly, p53 binding sites containing these more flexible d(CATG.CATG) sequence elements also exhibit a stronger tendency for tetrameric binding of the p53 DNA binding domain peptide. Examination of the shapes of DNA microcircles obtained by circularization of oligomers constructed from such flexible p53 target DNA sequences in tandem using MacMode atomic force microscopy directly revealed sequence-specific kinks in solution. The tetra-nucleotide sequence d(CATG.CATG) is highly conserved in most functional p53 response elements. Consequently, we propose that the sequence-specific kinks originating from d(CATG.CATG) sequences could be a common structural theme in p53 response elements and as evident from the results reported here, could be a determinant of binding site recognition by the p53 protein and the subsequent stability of the p53-DNA complex.

ADAPT: a molecular mechanics approach for studying the structural properties of long DNA sequences

We describe an original approach to determining sequence-structure relationships for DNA. This approach, termed ADAPT, combines all-atom molecular mechanics with a multicopy algorithm to build nucleotides that contain all four standard bases in variable proportions. These nucleotides enable us to search very rapidly for base sequences that energetically favor chosen types of DNA deformation or chosen DNA-protein or DNA-ligand interactions. Sequences satisfying the chosen criteria can be found by energy minimization, combinatorial sequence searching, or genome scanning, in a manner similar to the threading approaches developed for protein structure prediction. In the latter case, we are able to analyze roughly 2000 base pairs per second. Applications of the method to DNA allomorphic transitions, DNA deformation, and specific DNA interactions are presented.

Relationship between whole proteome aminoacid composition and static DNA curvature

To study possible relationships between an organism's genomic DNA curvature and the aminoacid composition of its proteome, every peptidic sequence from fully determined genomes was retrotranslated using the E. coli codon preferences, and the curvature profiles of the resulting DNA sequences were calculated and compared. A clear interdependence between these two variables was observed, as each retrotranslated proteome presented a distinctive, statistically significant DNA curvature profile biased toward its natural DNA curvature profile. In addition, by comparing the profiles arising from real and randomly permuted proteomes, we also found a position-dependent contribution of the peptidic sequence to DNA curvature. The implications of these results support the idea of a possible selection toward a specific global curvature of genomes.

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.

The global intrinsic curvature of archaeal and eubacterial genomes is mostly contained in their dinucleotide composition and is probably not an adaptation

Until now, the genomic DNA of all eubacteria analyzed has been hyper-curved, its global intrinsic curvature being higher than that of a random sequence. In contrast, that rule failed for archaea or eukaryotes, which could be either hypo- or hyper-curved. The existence of the rule suggested that, at least for eubacteria, global intrinsic curvature is adaptive. However, the present results from analyzing 21 eubacterial and six archaeal genomes argue against adaptation. First, there are two eubacterial exceptions to the former rule. More significantly, we found that the dinucleotide composition of the genome alone (which lacks all sequence information) is enough to determine the genome curvature. Additional evidence against adaptation came from showing that the global curvature of bacterial genomes could not have evolved under either of two complementary models of curvature selection: (i) that curvature is selected locally from unbiased variability; (ii) that curvature is established globally through the selection of a curvature-altering mutational bias. We found that the observed relationship between curvature and dinucleotide composition is incompatible with model (i). We also found that, contrary to the predictions of model (ii), the dinucleo-tide compositions of bacterial genomes were not statistically special in their curvature-related properties (when compared to stochastically generated dinucleotide compositions).

Inducible expression of the chromosomal cdiA from Citrobacter diversus NF85, encoding an ambler class A beta-lactamase, is under similar genetic control to the chromosomal ampC, encoding an ambler class C enzyme, from Citrobacter freundii OS60

This study aimed to characterize the molecular basis of beta-lactamase induction in Citrobacter diversus. The chromosomal beta-lactamase encoding region from C. diversus, strain NF85, was cloned and expressed in Escherichia coli. The cloned region was sequenced and open-reading frames encoding a class A beta-lactamase, designated cdiA, and a putative LysR-type transcriptional regulator protein, divergently transcribed from the beta-lactamase gene and designated cdiR, were identified. The nucleotide sequence of the NF85 cdiA was identical to that of the published C. diversus ULA27 ampC sequence. A putative helix-turn-helix DNA-binding motif was located at the N-terminus of CdiR, and homology with enterobacterial AmpR proteins was noted. CdiR was demonstrated to bind to the C. diversus cdiAR intergenic region but not to the C. freundii ampCR intergenic region. A putative CdiR binding motif was identified in the cdiAR intergenic region. The cloned cdiAR region was inducible in E. coli strains SNO3 and HfrH. The inducible phenotype was dependent on the E. coli ampD and ampG gene products. We conclude that the molecular basis of inducible cdiA expression in C. diversus is similar to that of C. freundii ampC.

31P NMR analysis of the DNA conformation induced by protein binding SRY/DNA complexes

Complexes of the HMG box protein SRY with two duplexes of 8 and 14 base pairs have been studied by 31P NMR and complete assignment of all phosphorus signals of the bound DNA duplexes are presented. While for the free DNA, all 31P signals display limited spectral dispersion (< 0.8 p.p.m.) for the bound duplexes, 31P resonances are spread over 2 p.p.m. Based on the previously published 3D structure of hSRY-HMG, with the 8 mer it is demonstrated that the upfield shifted resonances correspond to the site of partial intercalation of an isoleucine side chain into the DNA. Moreover, the observation of significant difference in linewidths between the two duplexes allows to estimate lifetime of the complexes from 31P-31P 2D exchange experiments.

Insertional mutation of orfD of the DCW cluster of Streptococcus pneumoniae attenuates virulence

Mutational analysis of a 5.5 kb fragment of the genome Streptococcus pneumoniae led to the identification of a putative new virulence gene, designated orfD. Insertion mutagenesis of flanking genes on the fragment suggested that the corresponding gene products were required for in vitro growth. In contrast, insertion mutation of orfD did not alter in vitro growth or the transformability pattern of the mutated strain. However, it did reduce bacterial growth in mice and attenuated virulence in an intraperitoneal model of infection. orfD is flanked by orfC (63 codons) and ftsL (105 codons) and all three genes are upstream of pbpx. orfC showed no similarity with other known proteins. ftsL of S. pneumoniae exhibits minimal sequence similarity with ftsL of E. coli, but shares 16% identical residues with the ftsL homologue encoded by ylld of B. subtilis. Also, ftsL of S. pneumoniae has a predicted topology similar to that described for ftsL of E. coli. Putative promoters with an extended -10 box could be identified upstream of both orfC or orfD. The four open reading frames (including pbpx) are orientated in the same direction, and polycistronic transcription could theoretically start at either promoter. Interestingly, this region shows organizational and sequence homologies with genes controlling division and cell wall biosynthesis (DCW) in other bacteria. The attenuation of virulence in the orfD insertion mutant might be due to the loss of function of the orfD gene product or to an altered level of expression of downstream genes.

DNA structure and polymerase fidelity

The accuracy of DNA replication results from both the intrinsic DNA polymerase fidelity and the DNA sequence. Although the recent structural studies on polymerases have brought new insights on polymerase fidelity, the role of DNA sequence and structure is less well understood. Here, the analysis of the crystal structures of hotspots for polymerase slippage including (CA)n and (A)n tracts in different intermolecular contexts reveals that, in the B-form, these sequences share common structural alterations which may explain the high rate of replication errors. In particular, a two-faced "Janus-like" structure with shifted base-pairs in the major groove but an apparent normal geometry in the minor groove constitutes a molecular decoy specifically suitable to mislead the polymerases. A model of the rat polymerase beta bound to this structure suggests that an altered conformation of the nascent template-primer duplex can interfere with correct nucleotide incorporation by affecting the geometry of the active site and breaking the rules of base-pairing, while at the same time escaping enzymatic mechanisms of error discrimination which scan for the correct geometry of the minor groove.In contrast, by showing that the A-form greatly attenuates the sequence-dependent structural alterations in hotspots, this study suggests that the A-conformation of the nascent template-primer duplex at the vicinity of the polymerase active site will contribute to fidelity. The A-form may play the role of a structural buffer which preserves the correct geometry of the active site for all sequences. The detailed comparison of the conformation of the nascent template-primer duplex in the available crystal structures of DNA polymerase-DNA complexes shows that polymerase beta, the least accurate enzyme, is unique in binding to a B-DNA duplex even close to its active site. This model leads to several predictions which are discussed in the light of published experimental data.

Damage increases the flexibility of duplex DNA

It is proposed that much of the recognition of specific types of damaged DNAs is based on accessible structural features, while much of the recognition of damaged DNAs, as a class, is based on flexibility. The more flexible a DNA the faster its diffusion rate. The diffusion rates of each member of a series of damaged duplex DNAs has been found to be significantly faster than that of the corresponding undamaged duplex DNA. The damaged sites studied include apurinic and apyrimidinic a basic sites, thymine glycol and urea. The presence of mismatched sites also increases the diffusion. Thus, damaged DNAs appear to have sufficient flexibility for recognition and the flexibility may allow damaged sites to act as a universal joint or hinge that allows distant sites on the DNA to come together.

The role of DNA-protein salt bridges in molecular recognition: a model study

A theoretical study is presented of the influence of salt bridges between cationic side chains and DNA phosphates on DNA conformation and flexibility. The DNA sequence studied is that of the catabolite activator protein binding oligomer from the crystallized complex. The effect of salt bridges is modeled by neutralization of net phosphate charges for the groups involved in such interactions in the crystallized complex. Energy-optimized conformations are obtained by molecular mechanics using the JUMNA program. Base sequence dependence is studied by moving the phosphate neutralization pattern along the sequence and also by point mutations. Normal mode analysis is used to evaluate DNA flexibility. The results obtained show that the free oligomer is already precurved in the direction favored by the protein, and the effect of phosphate neutralization is principally to increase this curvature. This effect is, however, strongly sequence dependent. In addition, it is shown that oligomer flexibility cannot be explained by a simple superposition of the properties of successive dinucleotide steps, strong long-range coupling effects are observed. In all the cases examined, phosphate neutralization, however, leads to a reduction in oligomer flexibility.

DNA bending induced by the archaebacterial histone-like protein MC1

The conformational changes induced by the binding of the histone-like protein MC1 to DNA duplexes have been analyzed by dark-field electron microscopy and polyacrylamide gel electrophoresis. Visualisation of the DNA molecules by electron microscopy reveals that the binding of MC1 induces sharp kinks. Linear DNA duplexes (176 bp) which contained a preferential site located at the center were used for quantitative analysis. Measurements of the angle at the center of all duplexes, at a fixed DNA concentration, as a function of the MC1 concentration, were very well fitted by a simple model of an isotropic flexible junction and an equilibrium between the two conformations of DNA with bound or unbound MC1. This model amounts to double-folded Gaussian distributions and yields an equilibrium deflection angle of theta0=116 degrees for the DNA with bound MC1. It allowed measurements of the fraction of DNA with bound MC1 to be taken as a function of MC1 concentrations and yields an equilibrium dissociation constant of Kd=100 nM. It shows that the flexibility of DNA is reduced by the binding of MC1 and the formation of a kink. The equilibrium dissociation constant value was corroborated by gel electrophoresis. Control of the model by the computation of the reduced chi2 shows that the measurements are consistent and that electron microscopy can be used to quantify precisely the DNA deformations induced by the binding of a protein to a preferential site.

Characterization of CorR, a transcriptional activator which is required for biosynthesis of the phytotoxin coronatine

Coronatine (COR) is a plasmid-encoded phytotoxin synthesized by several pathovars of phytopathogenic Pseudomonas syringae. The COR biosynthetic gene cluster in P. syringae pv. glycinea PG4180 is encoded by a 32-kb region which contains both the structural and regulatory genes needed for COR synthesis. The regulatory region contains three genes: corP, corS, and corR. corS is thought to function as a histidine protein kinase, whereas corP and corR show relatedness to response regulators of the two-component regulatory paradigm. In the present study, we investigated whether CorR is a positive activator of COR gene expression. We also studied whether CorR specifically binds the DNA region located upstream of cfl, a gene located at the 5' end of the gene cluster encoding coronafacic acid, the polyketide portion of COR. Complementation analysis with a corR mutant, PG4180.P2, and transcriptional fusions to a promoterless glucuronidase gene (uidA) indicated that CorR functions as a positive regulator of COR gene expression. Deletion analysis of the 5' end of the cfl upstream region was used to define the minimal region required for COR gene expression. A 360-bp DNA fragment located over 500 bp upstream from the cfl transcriptional start site was used in DNase I protection assays to define the specific bases bound by CorR. An area extending from -704 to -650 with respect to the cfl transcriptional start site was protected by DNase I footprinting, indicating a rather large area of protection. This area was also conserved in the promoter region for cmaA, which encodes a transcript containing genes for coronamic acid synthesis, another intermediate in the COR biosynthetic pathway. The results obtained in the current study suggest that both the coronafacic acid and the coronamic acid structural genes are controlled by CorR, a positive activator of COR gene expression.

DNA curvature and phosphate neutralization: an important aspect of specific protein binding

A theoretical study is presented of the influence of salt bridges between protein cationic side chains and DNA phosphates on DNA conformation and flexibility. Two DNA sequences are studied containing respectively the HNF3 and CAP binding sites. The effect of salt bridges is modelled by the neutralisation of net phosphate charges for the groups involved in such interactions in the complex. Energy optimised conformations are obtained by molecular mechanics calculations using the JUMNA program. Base sequence dependence is studied by moving the phosphate neutralisation pattern along the sequence, while normal mode analysis is used to evaluate DNA flexibility. The results show that phosphate neutralisation has a strong influence on DNA conformation. For the HNF3 binding sequence, the free oligomer is bent in direction very different from that observed in the protein complex. Phosphate neutralisation changes this direction by 45 degrees to within only 4 degrees of the direction in the complex, without changing the bending angle. For the CAP binding sequence, the free oligomer is already intrinsically curved in the direction favoured by the protein, but phosphate neutralisation increases the bending angle. For both oligomers studied these effects are strongly sequence dependent. It is also shown that oligomer flexibility cannot be explained by a simple superposition of the properties of successive dinucleotide steps. Important long range coupling effects are observed. However, for both sequence studied, phosphate neutralisation however leads to a reduction in oligomer flexibility.

Oxygen-insensitive nitroreductases: analysis of the roles of nfsA and nfsB in development of resistance to 5-nitrofuran derivatives in Escherichia coli

Nitroheterocyclic and nitroaromatic compounds constitute an enormous range of chemicals whose potent biological activity has significant human health and environmental implications. The biological activity of nitro-substituted compounds is derived from reductive metabolism of the nitro moiety, a process catalyzed by a variety of nitroreductase activities. Resistance of bacteria to nitro-substituted compounds is believed to result primarily from mutations in genes encoding oxygen-insensitive nitroreductases. We have characterized the nfsA and nfsB genes of a large number of nitrofuran-resistant mutants of Escherichia coli and have correlated mutation with cell extract nitroreductase activity. Our studies demonstrate that first-step resistance to furazolidone or nitrofurazone results from an nfsA mutation, while the increased resistance associated with second-step mutants is a consequence of an nfsB mutation. Inferences made from mutation about the structure-function relationships of NfsA and NfsB are discussed, especially with regard to the identification of flavin mononucleotide binding sites. We show that expression of plasmid-carried nfsA and nfsB genes in resistant mutants restores sensitivity to nitrofurans. Among the 20 first-step and 53 second-step mutants isolated in this study, 65 and 49%, respectively, contained insertion sequence elements in nfsA and nfsB. IS1 integrated in both genes, while IS30 and IS186 were found only in nfsA and IS2 and IS5 were observed only in nfsB. Insertion hot spots for IS30 and IS186 are indicated in nfsA, and a hot spot for IS5 insertion is evident in nfsB. We discuss potential regional and sequence-specific determinants for insertion sequence element integration in nfsA and nfsB.

A novel class of supercoil-independent nuclease hypersensitive site is comprised of alternative DNA structures that flank eukaryotic genes

The cell makes a fundamental distinction between genes and non-gene sequences, which mechanistically underlies the process of gene regulation. Here, we describe the properties of a novel class of genetic sites that reproducibly flank and delineate the coding regions of the eukaryotic genes tested. Defined in vitro reaction conditions that include altered solvation and elevated temperature rendered the sites hypersensitive to nuclease cleavage. Consequently, the complete coding regions of the Drosophila genes tested were quantitatively excised from genomic DNA or genomic clones by this treatment. Identical reaction products were generated from linear or supercoiled DNA substrates. Chemical modification and fine-structure analysis of several cleavage sites flanking Drosophila genes showed that the cleavage sites were stable nucleic acid structures that contained specific arrangements of paired and unpaired nucleotides. The locations and properties of the cleavage sites did not correspond to previously known nuclease hypersensitive sites nor to known alternative DNA structures. Thus, they appear to represent a new class of genetic site. In a deletion analysis, the minimal sequence information necessary to direct in vitro nuclease cleavage 3' to the Drosophila GART gene co-localized with the signal required for termination of transcription in vivo. The data suggest that a novel class of DNA site with distinct structural properties encodes biological information by marking the boundaries of at least some gene expression units in organisms as diverse as Plasmodium and Drosophila.

A role for H-NS in the regulation of the virF gene of Shigella and enteroinvasive Escherichia coli

We have investigated the role of H-NS, one of the major components of the bacterial nucleoid, in the expression of the virF gene present on the large virulence plasmid of Shigella and enteroinvasive Escherichia coli in response to different environmental conditions. VirF is an AraC-like protein which activates at least two promoters, virB and virG, both repressed by H-NS. Band shift experiments reveal that the affinity of H-NS for the virF and virB promoters is comparable, while the affinity for the virG promoter is higher. Polyacrylamide gel electrophoresis of three DNA fragments containing the virF, the virB and the VirG promoters demonstrates, in agreement with computer predictions, that they have an intrinsically curved structure, confirming the preference of H-NS for bent DNA. In vivo transcriptional analysis of virF mRNA shows that H-NS negatively controls the expression of virF at 30 degrees C. The expression of a virF-lacZ translational fusion in E.coli wild type and in an hns-defective derivative grown at 30 degrees or 37 degrees C and at pH 6.0 or 7.0 indicates that, in the absence of H-NS, virF expression becomes insensitive to temperature and to limited pH changes. Our results strongly suggest that H-NS controls virF expression by binding to the virF promoter and by repressing its expression at low temperature and at low pH.

Clues and consequences of DNA bending in transcription

This review attempts to substantiate the notion that nonlinear DNA structures allow prokaryotic cells to evolve complex signal integration devices that, to some extent, parallel the transduction cascades employed by higher organisms to control cell growth and differentiation. Regulatory cascades allow the possibility of inserting additional checks, either positive or negative, in every step of the process. In this context, the major consequence of DNA bending in transcription is that promoter geometry becomes a key regulatory element. By using DNA bending, bacteria afford multiple metabolic control levels simply through alteration of promoter architecture, so that positive signals favor an optimal constellation of protein-protein and protein-DNA contacts required for activation. Additional effects of regulated DNA bending in prokaryotic promoters include the amplification and translation of small physiological signals into major transcriptional responses and the control of promoter specificity for cognate regulators.

A novel bend of DNA CIT: changeable bending-center sites of an intrinsic curvature under temperature conditions

We found a novel DNA curvature, which has changeable bending-center sites of an intrinsic curvature under temperature conditions (CIT) in the cyanobacterium strain Microcystis aeruginosa K-81. Circular permutation analyses (CPA) for CIT under different temperature conditions (4-50 degrees C) revealed that the changeable bending-center sites are located in the 5'-upstream region (-141 to -184) of the psbA2 gene, encoding the D1 protein homolog for photosynthesis. The nucleotide sequence around the bending center contains several dT (deoxy thymine) tracts, which seem to be a pivotal determinant for CIT.

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.

Action at a distance for negative control of transcription of the glpD gene encoding sn-glycerol 3-phosphate dehydrogenase of Escherichia coli K-12

Aerobic sn-glycerol 3-phosphate dehydrogenase is a cytoplasmic membrane-associated respiratory enzyme encoded by the glpD gene of Escherichia coli. The glpD operon is tightly controlled by cooperative binding of the glp repressor to tandem operators (O(D)1 and O(D)2) that cover the -10 promoter element and 30 bp downstream of the transcription start site. In this work, two additional operators were identified within the glpD structural gene at positions 568 to 587 (0(D)3) and 609 to 628 (0(D)4). The two internal operators bound the glp repressor in the presence or absence of the tandem operators (O(D)1 and O(D)2) in vitro, as shown by DNase I footprinting. To assess a potential regulatory role for the two internal operators in vivo, a glpD-lacZ transcriptional fusion containing all four operators was constructed. The response of this fusion to the glp repressor was compared with those of fusion constructs in which O(D)3 and O(D)4 were inactivated by either deletion or site-directed mutagenesis. It was found that the repression conferred by binding of the glp repressor to O(D)1 and O(D)2 was increased five- to sevenfold upon introduction of the internal operators. A regulatory role for HU was suggested when it was found that repressor-mediated control of glpD transcription was increased fourfold in strains containing HU compared with that of strains deficient in HU. The effect of HU was apparent only in the presence of all four glpD operators. The results suggest that glpD is controlled by formation of a repression loop between the tandem and internal operators. HU may assist repression by bending the DNA to facilitate loop formation.

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.

Identification of Simian virus 40 promoter DNA sequences capable of conferring restriction endonuclease hypersensitivity

The simian virus 40 (SV40) DNA sequences found in the enhancer domain, nucleotides (nt) 103 to 177, and the early domain, nt 5149 to 5232, of the SV40 promoter have been analyzed for their ability to confer restriction endonuclease hypersensitivity in SV40 chromatin by using an SV40-based recombinant reporter system. The reporter system consists of a polylinker of various unique restriction endonuclease recognition sequences introduced into SV40 at nt 2666. We observed that the introduction of the enhancer domain at one end of the reporter and the early domain at the other end of the reporter resulted in a 20% increase in nuclease sensitivity within the reporter. In the enhancer domain, an element capable of conferring hypersensitivity was found between nt 114 and 124 with the sequence 5'CTGACTAATTG3', which has previously been shown to be the SV40 AP-1 binding site. In the early domain, an element capable of conferring hypersensitivity was localized to nt 5164 to 5187 and had the sequence 5'CATTTGCAAAGCTTTTTGCAAAAGC3'.

Molecular modelling of (A4T4NN)n and (T4A4NN)n: sequence elements responsible for curvature

The molecular modelling program JUMNA has been used to investigate the origins of the strikingly different curvature of the two sequences, (A4T4NN)n and (T4A4NN)n. Gel electrophoresis and cyclisation studies have shown that only the former of these two sequences is significantly curved. By developing novel superhelical symmetry constraints we were able to study the energetic and structural aspects of polymeric DNA having a controlled curvature. The results obtained (which do not take into account specific hydration effects) correlate well with the experimental data and offer a molecular level explanation of curvature. Although curvature is found to be initiated by specific dinucleotide junctions, deformations spread to surrounding dinucleotide steps and, moreover, sequence effects beyond the dinucleotide level are observed.

Examining the contribution of a dA+dT element to the conformation of Escherichia coli integration host factor-DNA complexes

DNA binding proteins that induce structural changes in DNA are common in both prokaryotes and eukaryotes. Integration host factor (IHF) is a multi-functional DNA binding and bending protein of Escherichia coli that can mediate protein-protein and protein-DNA interactions by bending DNA. Previously we have shown that the presence of a dA+dT element 5'-proximal to an IHF consensus sequence can affect the binding of IHF to a particular site. In this study the contribution of various sequence elements to the formation of IHF-DNA complexes was examined. We show that IHF bends DNA more when it binds to a site containing a dA+dT element upstream of its core consensus element than to a site lacking a dA+dT element. We demonstrate that IHF can be specifically crosslinked to DNA with binding sites either containing or lacking this dA+dT element. These results indicate the importance of flanking DNA and a dA+dT element in the binding and bending of a site by IHF.

The effects of sequence context on DNA curvature

Recent experiments have exposed significant discrepancies between experimental data and predictive models for DNA structure. These results strongly suggest that DNA structural parameters incorporated in the models are not always sufficient to account for the influence of sequence context and of specific ion effects. In an attempt to evaluate these two effects, we have investigated repetitive DNA sequences with the sequence motif GAGAG.CTCTC located in different helical phasing arrangements with respect to poly(A) tracts and GGGCCC.GGGCCC sequence motifs. Methods used are ligase-mediated cyclization and gel mobility experiments along with DNase I cutting and chemical probe studies. The results provide new evidence for curvature in poly(A) tracts. They also show that the sequence context in which bending and flexible sequence elements are found is an important aspect of sequence-dependent DNA conformation. Although dinucleotide models generally have good predictive power, this work demonstrates that in some instances sequence elements larger than the dinucleotide must be taken into account, and hence it provides a starting point for the appropriate modification and refinement of existing structural models for DNA.