CA runs increase DNA flexibility in the complex of lambda Cro protein with the OR3 site

Molecular Characterization and Genome Mechanical Features of Two Newly Isolated Polyvalent Bacteriophages Infecting Pseudomonas syringae pv. garcae

Coffee plants have been targeted by a devastating bacterial disease, a condition known as bacterial blight, caused by the phytopathogen Pseudomonas syringae pv. garcae (Psg). Conventional treatments of coffee plantations affected by the disease involve frequent spraying with copper- and kasugamycin-derived compounds, but they are both highly toxic to the environment and stimulate the appearance of bacterial resistance. Herein, we report the molecular characterization and mechanical features of the genome of two newly isolated (putative polyvalent) lytic phages for Psg. The isolated phages belong to class Caudoviricetes and present a myovirus-like morphotype belonging to the genuses Tequatrovirus (PsgM02F) and Phapecoctavirus (PsgM04F) of the subfamilies Straboviridae (PsgM02F) and Stephanstirmvirinae (PsgM04F), according to recent bacterial viruses' taxonomy, based on their complete genome sequences. The 165,282 bp (PsgM02F) and 151,205 bp (PsgM04F) genomes do not feature any lysogenic-related (integrase) genes and, hence, can safely be assumed to follow a lytic lifestyle. While phage PsgM02F produced a morphogenesis yield of 124 virions per host cell, phage PsgM04F produced only 12 virions per host cell, indicating that they replicate well in Psg with a 50 min latency period. Genome mechanical analyses established a relationship between genome bendability and virion morphogenesis yield within infected host cells.

Ultrafast differential flexibility of Cro-protein binding domains of two operator DNAs with different sequences

The crucial ultrafast domain fluctuation of the operator DNA O_R3 over O_R2 upon complexation with the repressor Cro-protein dimer has been investigated.

Unusual sequence effects on nucleotide excision repair of arylamine lesions: DNA bending/distortion as a primary recognition factor

The environmental arylamine mutagens are implicated in the etiology of various sporadic human cancers. Arylamine-modified dG lesions were studied in two fully paired 11-mer duplexes with a -G*CN- sequence context, in which G* is a C8-substituted dG adduct derived from fluorinated analogs of 4-aminobiphenyl (FABP), 2-aminofluorene (FAF) or 2-acetylaminofluorene (FAAF), and N is either dA or dT. The FABP and FAF lesions exist in a simple mixture of 'stacked' (S) and 'B-type' (B) conformers, whereas the N-acetylated FAAF also samples a 'wedge' (W) conformer. FAAF is repaired three to four times more efficiently than FABP and FAF. A simple A- to -T polarity swap in the G*CA/G*CT transition produced a dramatic increase in syn-conformation and resulted in 2- to 3-fold lower nucleotide excision repair (NER) efficiencies in Escherichia coli. These results indicate that lesion-induced DNA bending/thermodynamic destabilization is an important DNA damage recognition factor, more so than the local S/B-conformational heterogeneity that was observed previously for FAF and FAAF in certain sequence contexts. This work represents a novel 3'-next flanking sequence effect as a unique NER factor for bulky arylamine lesions in E. coli.

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.

Sequence-specific physical properties of African green monkey alpha-satellite DNA contribute to centromeric heterochromatin formation

Satellite DNA, a major component of eukaryotic centromeric heterochromatin, is potentially associated with the processes ensuring the faithful segregation of the genetic material during cell division. Structural properties of alpha-satellite DNA (AS) from African green monkey (AGM) were studied. Atomic force microscopy imaging showed smaller end-to-end distances of AS fragments than would be expected for the persistence length of random sequence DNA. The apparent persistence length of the AS was determined as 35nm. Gel-electrophoresis indicated only a weak contribution of intrinsic curvature to the DNA conformations suggesting an additional contribution of an elevated bending flexibility to the reduced end-to-end distances. Next, the force-extension behavior of the naked AS and in complex with nucleosomes was studied using optical tweezers. The naked AS showed a reduced overstretching transition force (-18% the value determined for random DNA) and higher forces required to straighten the DNA. Finally, reconstituted AS nucleosomes disrupted at significantly higher forces as compared with random DNA nucleosomes which is probably due to structural properties of the AS which stabilize the nucleosomes. The data support that the AS plays a role in the formation of centromeric heterochromatin due to specific structural properties and suggest that a relatively higher mechanical stability of nucleosomes is important in AGM-AS chromatin.

Specific Interactions of Divalent Metal Ions with a DNA Duplex Containing the d(CA)n/(GT)nTandem Repeat

Divalent metal ions are essential for maintaining functional states of the DNA molecule. Their participation in DNA structure is modulated by the base sequence and varies depending on the nature of the ion. The present investigation addresses the interaction of Ca2+ ions with a tandem repeat of two CA dinucleotides, (CA)2/(TG)2. The binding of Ca2+ to the repeat is monitored by nuclear magnetic resonance (NMR) spectroscopy using chemical shift mapping. Parallel experiments monitor binding of Mg2+ ions to the repeat as well as binding of each ion to a DNA duplex in which the (CA)2/(TG)2 repeat is eliminated. The results reveal that the direction and the magnitude of chemical shift changes induced by Ca2+ ions in the NMR spectra of the repeat are different from those induced by Mg2+ ions. The differences between the two cations are significantly diminished by the elimination of the (CA)2/(TG)2 repeat. These findings suggest a specific interaction of Ca2+ ions with the (CA)2/(TG)2 motif. The specificity of the interaction resides in the two A-T base pairs of the repeat, and it involves the major groove of the first A-T base pair and both grooves of the second A-T base pair.

Indirect readout of DNA sequence at the primary-kink site in the CAP-DNA complex: recognition of pyrimidine-purine and purine-purine steps

The catabolite activator protein (CAP) bends DNA in the CAP-DNA complex, typically introducing a sharp DNA kink, with a roll angle of approximately 40 degrees and a twist angle of approximately 20 degrees, between positions 6 and 7 of the DNA half-site, 5'-A1A2A3T4G5T6G7A8T9C10T11 -3' ("primary kink"). In previous work, we showed that CAP recognizes the nucleotide immediately 5' to the primary-kink site, T6, through an "indirect-readout" mechanism involving sequence effects on energetics of primary-kink formation. Here, to understand further this example of indirect readout, we have determined crystal structures of CAP-DNA complexes containing each possible nucleotide at position 6. The structures show that CAP can introduce a DNA kink at the primary-kink site with any nucleotide at position 6. The DNA kink is sharp with the consensus pyrimidine-purine step T6G7 and the non-consensus pyrimidine-purine step C6G7 (roll angles of approximately 42 degrees, twist angles of approximately 16 degrees ), but is much less sharp with the non-consensus purine-purine steps A6G7 and G6G7 (roll angles of approximately 20 degrees, twist angles of approximately 17 degrees). We infer that CAP discriminates between consensus and non-consensus pyrimidine-purine steps at positions 6-7 solely based on differences in the energetics of DNA deformation, but that CAP discriminates between the consensus pyrimidine-purine step and non-consensus purine-purine steps at positions 6-7 both based on differences in the energetics of DNA deformation and based on qualitative differences in DNA deformation. The structures further show that CAP can achieve a similar, approximately 46 degrees per DNA half-site, overall DNA bend through a sharp DNA kink, a less sharp DNA kink, or a smooth DNA bend. Analysis of these and other crystal structures of CAP-DNA complexes indicates that there is a large, approximately 28 degrees per DNA half-site, out-of-plane component of CAP-induced DNA bending in structures not constrained by end-to-end DNA lattice interactions and that lattice contacts involving CAP tend to involve residues in or near biologically functional surfaces.

A highly distinctive mechanical property found in the majority of human promoters and its transcriptional relevance

A recent study revealed that TATA boxes and initiator sequences have a common anomalous mechanical property, i.e. they comprise distinctive flexible and rigid sequences when compared with the other parts of the promoter region. In the present study, using the flexibility parameters from two different models, we calculated the average flexibility profiles of 1004 human promoters that do not contain canonical promoter elements, such as a TATA box, initiator (Inr) sequence, downstream promoter element or a GC box, and those of 382 human promoters that contain the GC box only. Here, we show that they have a common characteristic mechanical property that is strikingly similar to those of the TATA box-containing or Inr-containing promoters. Their most interesting feature is that the TATA- or Inr-corresponding region lies in the several nucleotides around the transcription start site. We have also found that a dinucleotide step from −1 to +1 (transcription start site) has a slight tendency to adopt CA that is known to be flexible. We also demonstrate that certain synthetic DNA fragments designed to mimic the average mechanical property of these 1386 promoters can drive transcription. This distinctive mechanical property may be the hallmark of a promoter.

Structure of DNA sequence d-TGATCA by two-dimensional nuclear magnetic resonance spectroscopy and restrained molecular dynamics

The 5' d-TpG 3' element is a part of DNA sequences involved in regulation of gene expression and is also a site for intercalation of several anticancer drugs. Solution conformation of DNA duplex d-TGATCA containing this element has been investigated by two-dimensional NMR spectroscopy. Using a total of 12 torsional angles and 121 distance constraints, structural refinement has been carried out by restrained molecular dynamics (rMDs) in vacuum up to 100 ps. The structure is characterized by a large positive roll at TpG/CpA base pair step and large negative propeller twist for AT and TA base pairs. The backbone torsional angle, gamma(O5'-C5'-C4'-C3'), of T1 residue adopts a trans-conformation which is corroborated by short intra nucleotide T1H6-T1H5' (3.7A) distance in nuclear overhauser effect spectroscopy (NOESY) spectra while the backbone torsional angle, beta(P-O5'-C5'-C4'), exists in trans as well as gauche state for T1 and C5 residues. There is evidence of significant flexibility of the sugar-phosphate backbone with rapid inter-conversion between two different conformers at TpG/CpA base pair step. The base sequence dependent variations and local structural heterogeneity have important implications in specific recognition of DNA by ligands.

Intrinsic curvature in the VP1 gene of SV40: comparison of solution and gel results

DNA restriction fragments that are stably curved are usually identified by polyacrylamide gel electrophoresis because curved fragments migrate more slowly than normal fragments containing the same number of basepairs. In free solution, curved DNA molecules can be identified by transient electric birefringence (TEB) because they exhibit rotational relaxation times that are faster than those of normal fragments of the same size. In this article, the results observed in free solution and in polyacrylamide gels are compared for a highly curved 199-basepair (bp) restriction fragment taken from the VP1 gene in Simian Virus 40 (SV40) and various sequence mutants and insertion derivatives. The TEB method of overlapping fragments was used to show that the 199-bp fragment has an apparent bend angle of 46 +/- 2 degrees centered at sequence position 1922 +/- 2 bp. Four unphased A- and T-tracts and a mixed A3T4-tract occur within a span of approximately 60 bp surrounding the apparent bend center; for brevity, this 60-bp sequence element is called a curvature module. Modifying any of the A- or T-tracts in the curvature module by site-directed mutagenesis decreases the curvature of the fragment; replacing all five A- and T-tracts by random-sequence DNA causes the 199-bp mutant to adopt a normal conformation, with normal electrophoretic mobilities and birefringence relaxation times. Hence, stable curvature in this region of the VP1 gene is due to the five unphased A- and T- tracts surrounding the apparent bend center. Discordant solution and gel results are observed when long inverted repeats are inserted within the curvature module. These insertion derivatives migrate anomalously slowly in polyacrylamide gels but have normal, highly flexible conformations in free solution. Discordant solution and gel results are not observed if the insert does not contain a long inverted repeat or if the long inverted repeat is added to the 199-bp fragment outside the curvature module. The results suggest that long inverted repeats can form hairpins or cruciforms when they are located within a region of the helix backbone that is intrinsically curved, leading to large mobility anomalies in polyacrylamide gels. Hairpin/cruciform formation is not observed in free solution, presumably because of rapid conformational exchange. Hence, DNA restriction fragments that migrate anomalously slowly in polyacrylamide gels are not necessarily stably curved in free solution.

Structure of DNA Hexamer Sequence d-CGATCG by Two-dimensional Nuclear Magnetic Resonance Spectroscopy and Restrained Molecular Dynamics

Solution conformation of self-complementary DNA duplex d-CGATCG, containing 5' d-CpG 3' site for intercalation of anticancer drug, daunomycin and adriamycin, has been investigated by nuclear magnetic resonance (NMR) spectroscopy. Complete resonance assignments of all the protons (except some H5'/H5" protons) have been obtained following standard procedures based on double quantum filtered correlation spectroscopy (dQF COSY) and two-dimensional nuclear Overhauser effect (NOE) spectra. Analysis of sums of coupling constants in one-dimensional NMR spectra, cross peak patterns in dQF COSY spectra and inter proton distances shows that the DNA sequence assumes a conformation close to the B-DNA family. The deoxyribose sugar conformation is in dynamic equilibrium with predominantly S-type conformer and a minor N-type conformer with N<-->S equilibrium varying with temperature. At 325 K, the mole fraction of the N-conformer increases for some of the residues by approximately 9%. Using a total of 10 spin-spin coupling constants and 112 NOE intensities, structural refinement has been carried out using Restrained Molecular Dynamics (rMD) with different starting structures, potential functions and rMD protocols. It is observed that pseudorotation phase angle of deoxyribose sugar for A3 and T4 residues is approximately 180 degrees and approximately 120 degrees, respectively while all other residues are close to C2'endo-conformation. A large propeller twist (approximately -18 degrees) and smallest twist angle (approximately 31 degrees) at A3pT4 step, in the middle of the sequence, a wider (12 A) and shallower (3.0 A) major groove with glycosidic bond rotation as high anti at both the ends of hexanucleotide are observed. The structure shows base-sequence dependent variations and hence strong local structural heterogeneity, which may have implications in ligand binding.

Long-distance radical cation reactions in DNA three-way junctions: inter-arm interaction and migration through the junction

DNA three-way junctions (TWJ) are branched molecules having three 'arms'. We studied long-distance radical cation migration in these assemblies by incorporating anthraquinone (AQ) groups linked by a covalent tether to one strand of one arm of the TWJ. Excitation of the AQ at 350 nm results in one-electron oxidation of the DNA, which generates a base radical cation. This leads to relatively inefficient (compared with duplex DNA) strand cleavage at guanines following piperidine treatment of the irradiated samples. When the AQ is linked to the 5'-terminus of arm III by a flexible tether, gel electrophoretic analysis shows that strand cleavage occurs at the guanines in all three arms. We also investigated a TWJ in which the anthraquinone is specifically intercalated in arm III. In this case, a different pattern of strand cleavage is detected. We conclude that there are at least two mechanisms for long-distance radical cation migration in TWJs: (i) by inefficient charge hopping through the junction; (ii) by a through-space, cross-arm interaction when the AQ is on a flexible tether.

Transcriptional regulation of the CRLR gene in human microvascular endothelial cells by hypoxia

Adrenomedullin is a 52 amino acid peptide that shows a remarkable range of effects on the vasculature that include inter alia, vasodilatation, regulation of permeability, inhibition of endothelial cell apoptosis, and promotion of angiogenesis. Recently the G-protein coupled receptor (GPCR) calcitonin receptor-like receptor (CRLR), and receptor activity modifying proteins (RAMPs) have become recognized as integral components of the adrenomedullin signaling system. However, mechanisms of regulation of CRLR expression are still largely unknown. This is in part due to lack of information on the gene promoter. In this study we have determined the transcriptional start of human CRLR cDNA by 5'-RACE and cloned the proximal 5'-flanking region of the gene by PCR. The 2318 bp genomic fragment contains the basal promoter of human CRLR, including potential TATA-boxes and several GC boxes. Regulatory elements binding known transcription factors, such as Sp-1, Pit-1, glucocorticoid receptor, and hypoxia-inducible factor-1 alpha (HIF-1alpha) were also identified. When cloned into reporter gene vectors, the genomic fragment showed significant promoter activity, indicating that the 5'-flanking region isolated by PCR contains the gene promoter of human CRLR. Of significance is that the cloned promoter fragments were activated by hypoxia when transfected in primary microvascular endothelial cells. Site-directed mutagenesis of the consensus hypoxia-response element (HRE) in the 5'-flanking region abolished such a response. We also demonstrated by semi-quantitative RT-PCR that transcription of the gene is activated by hypoxia in microvascular endothelial cells. In contrast, expression of RAMPs 1, 2, and 3 was unaffected by low oxygen tension. We conclude that simultaneous transcriptional up-regulation of CRLR and its ligand adrenomedullin in endothelial cells could lead to a potent survival loop and therefore might play a significant role in vascular responses to hypoxia and ischemia.

Structural analysis of hemicatenated DNA loops

BACKGROUND

We have previously isolated a stable alternative DNA structure, which was formed in vitro by reassociation of the strands of DNA fragments containing a 62 bp tract of the CA-microsatellite poly(CA).poly(TG). In the model which was proposed for this structure the double helix is folded into a loop, the base of the loop consists of a DNA junction in which one of the strands of one duplex passes between the two strands of the other duplex, forming a DNA hemicatenane in a hemiknot structure. The hemiknot DNA structures obtained with long CA/TG inserts have been imaged by AFM allowing us to directly visualize the loops.

RESULTS

Here we have analyzed this structure with several different techniques: high-resolution gel electrophoresis, probing by digestion with single stranded DNA-specific nucleases or with DNase I, modification with chemicals specific for unpaired bases, and atomic force microscopy. The data show a change in DNA structure localized to the CA/TG sequence and allow us to better understand the structure of this alternative conformation and the mechanism of its formation.

CONCLUSIONS

The present work is in good agreement with the model of hemicatenated DNA loop proposed previously. In the presence of protein HMGB1, shifted reassociation of the strands of DNA fragments containing a tract of the poly(CA).poly(TG) microsatellite leads to the formation of DNA loops maintained at their base by a hemicatenated junction located within the repetitive sequence. No mobility of the junction along the DNA molecule could be detected under the conditions used. The novel possibility to prepare DNA hemicatenanes should be useful to further study this alternative DNA structure and its involvement in replication or recombination.

Visualization of hemiknot DNA structure with an atomic force microscope

The hemiknot, a novel type of DNA structure in which a loop is stabilized by threading one end of the duplex through another, has been studied in this paper. The hemiknot was obtained by reassociation of a DNA fragment with (CA/TG)n inserts of different lengths. Slow and fast migrating products were purified by gel electrophoresis and imaged by atomic force microscopy (AFM) using the aminopropylsilatrane-mica technique for sample preparation. Slow migrating product was characterized by the formation of small blobs for the short insert (60 bp) and clear loops and other morphologies for the long insert (188 bp). These structural features were found in almost 100% of the molecules of the slow migrating sample and were not present in the control sample. Measurements showed that the location of the blobs coincided with the positions of the inserts. The sample with the 188 bp insert in the 573 bp fragment had large structural irregularities. The majority of the molecules (77%) had asymmetrically located loops. The location of the loop in the molecules correlated well with the position of the insert in the fragment. The measured sizes of the loops were in agreement with the insert size. Altogether, these data support the hypothesis for hemiknot formation suggested earlier. In addition to looped structures, other morphologies of the hemiknot were identified in AFM images. Possible models for hemiknot formation and structure are discussed.

Indirect readout of DNA sequence at the primary-kink site in the CAP-DNA complex: DNA binding specificity based on energetics of DNA kinking

The catabolite activator protein (CAP) makes no direct contact with the consensus base-pair T:A at position 6 of the DNA half-site 5'-A(1)A(2)A(3)T(4)G(5)T(6)G(7)A(8)T(9)C(10)T(11)-3' but, nevertheless, exhibits strong specificity for T:A at position 6. Binding of CAP results in formation of a sharp DNA kink, with a roll angle of approximately 40 degrees and a twist angle of approximately 20 degrees, between positions 6 and 7 of the DNA half-site. The consensus base-pair T:A at position 6 and the consensus base-pair G:C at position 7 form a T:A/G:C step, which is known to be associated with DNA flexibility. It has been proposed that specificity for T:A at position 6 is a consequence of formation of the DNA kink between positions 6 and 7, and of effects of the T:A(6)/G:C(7) step on the geometry of DNA kinking, or the energetics of DNA kinking. In this work, we determine crystallographic structures of CAP-DNA complexes having the consensus base-pair T:A at position 6 or the non-consensus base-pair C:G at position 6. We show that complexes containing T:A or C:G at position 6 exhibit similar overall DNA bend angles and local geometries of DNA kinking. We infer that indirect readout in this system does not involve differences in the geometry of DNA kinking but, rather, solely differences in the energetics of DNA kinking. We further infer that the main determinant of DNA conformation in this system is protein-DNA interaction, and not DNA sequence.

Multimerization-cyclization of DNA fragments as a method of conformational analysis

Ligation of short DNA fragments results in the formation of linear and circular multimers of various lengths. The distribution of products in such a reaction is often used to evaluate fragment bending caused by specific chemical modification, by bound ligands or by the presence of irregular structural elements. We have developed a more rigorous quantitative approach to the analysis of such experimental data based on determination of j-factors for different multimers from the distribution of the reaction products. j-Factors define the effective concentration of one end of a linear chain in the vicinity of the other end. To extract j-factors we assumed that kinetics of the reaction is described by a system of differential equations where j-factors appear as coefficients. The assumption was confirmed by comparison with experimental data obtained here for DNA fragments containing A-tracts. At the second step of the analysis j-factors are used to determine conformational parameters of DNA fragments: the equilibrium bend angle, the bending rigidity of the fragment axis, and the total twist of the fragments. This procedure is based on empirical equations that connect the conformational parameters with the set of j-factors. To obtain the equations, we computed j-factors for a large array of conformational parameters that describe model fragments. The approach was tested on both simulated and actual experimental data for DNA fragments containing A-tracts. A-tract DNA bend angle determined here is in good agreement with previously published data. We have established a set of experimental conditions necessary for the data analysis to be successful.

Structure and dynamics of three-way DNA junctions: atomic force microscopy studies

We have used atomic force microscopy (AFM) to study the conformation of three-way DNA junctions, intermediates of DNA replication and recombination. Immobile three-way junctions with one hairpin arm (50, 27, 18 and 7 bp long) and two relatively long linear arms were obtained by annealing two partially homologous restriction fragments. Fragments containing inverted repeats of specific length formed hairpins after denaturation. Three-way junctions were obtained by annealing one strand of a fragment from a parental plasmid with one strand of an inverted repeat-containing fragment, purified from gels, and examined by AFM. The molecules are clearly seen as three-armed molecules with one short arm and two flexible long arms. The AFM analysis revealed two important features of three-way DNA junctions. First, three-way junctions are very dynamic structures. This conclusion is supported by a high variability of the inter-arm angle detected on dried samples. The mobility of the junctions was observed directly by imaging the samples in liquid (AFM in situ). Second, measurements of the angle between the arms led to the conclusion that three-way junctions are not flat, but rather pyramid-like. Non-flatness of the junction should be taken into account in analysis of the AFM data.

Isolation and characterization of the promoter region of the rat vasopressin V1b receptor gene

Regulation of pituitary vasopressin V1b receptors plays a critical role in regulating pituitary adrenocorticotropic hormone (ACTH) secretion during adaptation to stress. The objective of this study was to isolate the promoter regulatory region of the V1b receptor gene to better understand the molecular mechanisms involved in V1b receptor regulation. Screening of a rat genomic library using probes directed to the coding region and to the 5'UTR of the rat V1b receptor resulted in the isolation of several clones containing the 5'upstream regions of the V1b receptor cDNA. Sequencing of an 11.2 Kb fragment revealed 8.2 Kb upsteam of the reported cDNA sequence, which contains a putative promoter regulatory region. The 3' end of the clone contained 1472 base pairs corresponding to the recognized cDNA sequence, followed by 1506 bp of unknown sequence located at the end of the sixth transmembrane domain, probably corresponding to an intron, characteristic of these family of receptors. An additional 161 bp intron was found in the 5'UTR, similar to that described in the rat oxytocin receptor gene. 5'RACE and RNase protection analysis mapped two major putative transcription start points at -830 and -861 bp from the starting methionine. Analysis of the putative promoter region showed no indication of a proximal TATA box, but the presence of a CACA box, a GAGA box, several AP-1 and AP-2 sites and a cluster of Sp1 sites upstream of the AP-2 sites. A luciferase construct containing a 2.1-kb of putative promoter, and part of the 5'UTR including the first intron, showed promoter activity when transfected into COS-7, CHO and PC12 cell lines but not in AtT-20 cells. A similar construct without the intron and distal 5'UTR sequence has no promoter activity in the same cell lines. In summary, the V1b receptor gene contains at least 3 exons and 2 introns. The 5'flanking sequence contains several potential sites for transcriptional regulation, and induced luciferace activity only in constructs containing intron 1, suggesting that the latter is important for receptor gene activation. The data provide bases for future analysis of the regulatory elements controlling V1b receptor transcription.

Regulation of vasopressin V1b receptors in the anterior pituitary gland of the rat

Vasopressin secreted by parvocellular neurones of the hypothalamic paraventricular nucleus modulates pituitary adrenocorticotrophic hormone (ACTH) secretion by acting upon vasopressin V1b type receptors in the pituitary corticotroph coupled to phospholipase C. Regulation of V1b receptors contributes to the adaptation of the hypothalamic-pituitary-adrenal (HPA) axis to stress, as evidenced by the correlation between vasopressin receptor number and pituitary ACTH responsiveness. V1b receptor upregulation during chronic stress is associated with elevated circulating glucocorticoids and vasopressin expression in parvocellular neurones, suggesting that these factors control V1b receptor expression. Removal of circulating glucocorticoids by adrenalectomy causes sustained vasopressin receptor downregulation, but reduces V1b receptor mRNA only transiently. The latter effect is not mediated by increased corticotrophin-releasing hormone (CRH) and vasopressin release, since it is not prevented by lesions of the hypothalamic paraventricular nucleus. Adrenalectomy causes sustained V1b receptor loss in Brattleboro rats, which lack hypothalamic vasopressin, suggesting that vasopressin mediates V1b receptor mRNA recovery. Exogenous glucocorticoid administration downregulates pituitary vasopressin binding but increases V1b receptor mRNA and facilitates coupling of the receptor to phospholipase C, effects which may contribute to the refractoriness of vasopressin actions to glucocorticoid feedback. The lack of parallelism between changes in pituitary vasopressin binding and V1b receptor mRNA levels during manipulation of the HPA axis indicates that V1b receptor content depends on post-transcriptional mechanisms rather than steady-state V1b receptor mRNA levels. These studies suggest that interaction between glucocorticoids and vasopressin plays an important role in regulating V1b receptor mRNA expression during alterations of the HPA axis. In addition, the recent characterization of a major part of the V1b receptor gene provides a basis for studying the molecular mechanisms regulating the V1b receptor.

DNA distortion as a factor in nucleosome positioning

In a previous report we constructed a synthetic DNA sequence that directed the deposition of histone octamers to a single site, and it was proposed that DNA distortion was involved in the positioning effect. In the present study we utilized the chemical probe potassium permanganate to identify sites of DNA distortion in the synthetic positioning sequence. A permanganate hypersite was identified 15 bp from the nucleosome pseudo-dyad at a site known to display DNA distortion in the mature nucleosome. The sequence of the site contained a TA step flanked by an oligo-pyrimidine tract. A series of substitutions were made in the region of the permanganate hypersite and the resulting constructs tested for affinity for histone octamers and translational positioning in in vitro studies. The results revealed that either a single base substitution at the TA step or in the adjacent homopolymeric tract dramatically affected affinity and positioning activity. The rotational orientation of the permanganate-sensitive sequence was shown to be important for functions, since altering the orientation of the site in a positioning fragment reduced positioning activity and octamer affinity, while altering the rotational orientation of the sequence in a non-positioning fragment had the opposite effects. A reconstituted 5 S rDNA positioning sequence from Lytechinus variegatus was also shown to display a permanganate hypersite 16 bp from its pseudo-dyad.

Intermolecular contacts between the lambda-Cro repressor and the operator DNA characterized by nuclear magnetic resonance spectroscopy

The specific interaction between lambda phage Cro repressor and the DNA fragment bearing the consensus sequence of operators has been studied using nuclear magnetic resonance (NMR). Using both 15N- and 13C/15N- labeled lambda-Cro in complex with unlabeled DNA, chemical shift assignments of the lambda-Cro-DNA complex were obtained using heteronuclear NMR experiments. Inter-molecular contacts between the protein and DNA were identified using heteronuclear filtered NOESY experiments. The inter-molecular contacts were supplemented with intra-protein and intra-DNA NOE constraints to dock lambda-Cro to the bent B-form DNA using restrained molecular dynamics. The structure of one of the subunits of lambda-Cro in the complex is essentially the same as that of the unbound form. In the complex, inter-molecular NOEs were observed between the "helix-turn-helix" region comprising the alpha2 and alpha3 helices of the lambda-Cro protein and the major groove of the DNA. The methyl group of Thr17 forms a hydrophobic contact with the methyl group of the thymine at base pair 1 in the DNA, and Val25 and Ala29 make hydrophobic contacts with the methyl group of the thymine at base pair 5. The presence and the absence of these contacts can explain the difference in the affinity of lambda-Cro to several variants of the operator sequence.

Environmental influences on DNA curvature

DNA curvature plays an important role in many biological processes. To study environmental influences on DNA curvature we compared the anomalous migration on polyacrylamide gels of ligation ladders of 11 specifically-designed oligonucleotides. At low temperatures (25 degrees C and below) most of the sequences exhibited a degree of anomalous migration. Increased temperature had a significant effect on the anomalous migration (curvature) of some sequences but limited effects on others; at 50 degrees C only 1 sequence migrated anomalously. Mg2+ had a strong influence on the migration of certain sequences, whilst spermine enhanced the anomalous migration of a different set of sequences. Sequences with a GGC motif exhibited greater curvature than predicted by the presently-used angles for the nearest-neighbour wedge model and are especially sensitive to Mg2+. The data have implications for models for DNA curvature and for environmentally-sensitive DNA conformations in the regulation of gene expression.

Structure determination and analysis of helix parameters in the DNA decamer d(CATGGCCATG)2 comparison of results from NMR and crystallography

The solution structure of the DNA decamer (CATGGCCATG)2 has been determined by NMR spectroscopy and restrained molecular dynamic and distance geometry calculations. The restrainted data set includes interproton distances and torsion angles for the deoxyribose sugar ring which were obtained by nuclear Overhauser enhancement intensities and quantitative simulation of cross-peaks from double quantum filtered correlation spectroscopy. The backbone torsion angles were constrained using experimental data from NOE cross-peaks, 1H-1H and 1H-31P-coupling constants. The NMR structure and the crystal structure of the DNA decamer deviates from the structure of the canonical form of B-DNA in a number of observable characteristics. Particularly, both structures display a specific pattern of stacking interaction in the central GGC base triplet. Furthermore, a specific local conformation of the TG/CA base-pair step is present in NMR and crystal structure, highlighting the unusually high flexibility of this DNA duplex part. The solution structure of the TG/CA base-pair step obtained by our high resolution NMR study is characterized by a positive roll angle, whereas in crystal this base-pair step tends to adopt remarkably high twist angles.

A promoter-associated polymorphic repeat modulates PAX-6 expression in human brain

The PAX-6 gene plays a critical role in neurodevelopment and brain plasticity. While transcription of human PAX-6 is regulated by alternate usage of two distinct promoters termed A and B, expression in adult human brain is primarily controlled by promoter B. We now report that a novel polymorphic dinucleotide repeat sequence with the structure (AC)m(AG)n is located approximately 1 kb upstream of the transcription initiation site associated with promoter B. PCR-based systematic screening for length variations in a caucasian population showed a skewed distribution of multiple alleles containing between 24 and 36 repeat units. In 217 unrelated individuals, the frequency of alleles in the range between 25 and 29 repeats was 90%, with the 26 repeat allele alone accounting for 50%; the heterozygosity rate was 65%. Variants of this PAX-6 gene-linked polymorphic region (PAX-6LPR) had different transcriptional efficiencies when fused to a luciferase reporter gene and transfected into Cos-7 cells. Promoter activity of variants with >/=29 repeats was 4- to 9-fold higher than that of the 26 repeat allele. The influence of the PAX-6LPR on PAX-6 expression was confirmed in postmortem cerebellum from individuals with different genotypes. mRNA levels were 2-fold higher in genotypes with long alleles compared to those with short alleles. Allelic variation in PAX-6 expression may be a determinant of interindividual differences in brain plasticity and function.

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.

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.

Insertion of telomere repeat sequence decreases plasmid DNA condensation by cobalt (III) hexaammine

Telomere repeat sequence (TRS) DNA is found at the termini of most eukaryotic chromosomes. The sequences are highly repetitive and G-rich (e.g., [C(1-3)A/TG(1-3)]n for the yeast Saccharomyces cerevisiae) and are packaged into nonnucleosomal protein-DNA structures in vivo. We have used total intensity light scattering and electron microscopy to monitor the effects of yeast TRS inserts on in vitro DNA condensation by cobalt (III) hexaammine. Insertion of 72 bp of TRS into a 3.3-kb plasmid depresses condensation as seen by light scattering and results in a 22% decrease in condensate thickness as measured by electron microscopy. Analysis of toroidal condensate dimensions suggests that the growth stages of condensation are inhibited by the presence of a TRS insert. The depression in total light scattering intensity is greater when the plasmid is linearized with the TRS at an end (39-49%) than when linearized with the TRS in the interior (18-22%). Circular dichroism of a 95-bp fragment containing the TRS insert gives a spectrum that is intermediate between the A-form and B-form, and the anomalous condensation behavior of the TRS suggests a noncanonical DNA structure. We speculate that under conditions in which the plasmid DNA condenses, the telomeric insert assumes a helical geometry that is similar to the A-form and is incompatible with packing into the otherwise B-form lattice of the condensate interior.

DNA length and concentration dependencies of anisotropic phase transitions of DNA solutions

Critical concentrations for the isotropic to cholesteric phase transitions of double-stranded DNA fragments in simple buffered saline (0.1 M NaCl) solutions were determined as a function of DNA contour length ranging from approximately 50 nm to 2700 nm, by solid-state 31P NMR spectroscopy and polarized light microscopy. As expected for semirigid chains, the critical concentrations decrease sharply with increasing DNA length near the persistence length in the range from 50 to 110 nm, and approach a plateau when the contour length exceeds 190 nm. The biphasic region is substantially wider than observed for xanthan, another semirigid polyelectrolyte approximately twice as stiff as DNA, primarily because of low critical concentrations for first appearance of the anisotropic phase, C(i)*, in DNA samples > or =110 nm (320 base pairs) long. The limiting C(i)* for DNA > or =490 nm long is exceptionally low (only 13 mg/ml) and is substantially lower than the C(i)* of approximately 40 mg/ml reported for the stiffer xanthan polyelectrolyte. The much higher values of the critical concentrations, C(a)*, for the disappearance of the isotropic DNA phase (> or =67 mg/ml) are modestly higher than those observed for xanthan and are predicted reasonably well by a theory that has been applied to other semirigid polymers, if a DNA persistence length in the consensus range of 50-100 nm is assumed. By contrast, the broad biphasic region and low C(i)* values of DNA fragments > or =190 nm long could only be reconciled with theory by assuming persistence lengths of 200-400 nm. The latter discrepancies are presumed to reflect some combination of deficiencies in current theory as applied to chiral, strong polyelectrolytes such as DNA, and sequence-dependent variations in DNA properties such as flexibility, curvature, or interaction potential. The propensity of DNA to spontaneously self-order at low concentrations well in the physiological range may have biological significance.

Does TATA matter? A structural exploration of the selectivity determinants in its complexes with TATA box-binding protein

The binding of the TATA box-binding protein (TBP) to a TATA sequence in DNA is essential for eukaryotic basal transcription. TBP binds in the minor groove of DNA, causing a large distortion of the DNA helix. Given the apparent stereochemical equivalence of AT and TA basepairs in the minor groove, DNA deformability must play a significant role in binding site selection, because not all AT-rich sequences are bound effectively by TBP. To gain insight into the precise role that the properties of the TATA sequence have in determining the specificity of the DNA substrates of TBP, the solution structure and dynamics of seven DNA dodecamers have been studied by using molecular dynamics simulations. The analysis of the structural properties of basepair steps in these TATA sequences suggests a reason for the preference for alternating pyrimidine-purine (YR) sequences, but indicates that these properties cannot be the sole determinant of the sequence specificity of TBP. Rather, recognition depends on the interplay between the inherent deformability of the DNA and steric complementarity at the molecular interface.

Effects of natural and synthetic polyamines on the conformation of an oligodeoxyribonucleotide with the estrogen response element

We studied the effects of natural and synthetic polyamines on the conformation of an oligodeoxyribonucleotide (ODN1) harboring the estrogen response element (ERE) by circular dichroism (CD) spectroscopy and polyacrylamide gel electrophoresis. Putrescine and spermidine had no marked effect on the CD spectrum of ODN1. In contrast, spermine provoked and stabilized two characteristic changes in the CD spectrum. The first change was indicated by an increase in the intensity of the CD band at 280 nm at 0.5 mM spermine in Tris-HCl buffer containing 50 mM NaCl. This change appears to be related to changes in base tilt and conformational alterations similar to A-DNA. At 1-2 mM spermine, the CD spectrum was characterized by a loss of positive bands at 220 and 270 nm. This change might have contributions from polyamine-induced condensation/aggregation of DNA. Spectral measurements were also conducted in Tris-HCl buffer containing 150 mM NaCl to minimize contributions from condensation and aggregation of ODN1. Under these conditions, CD spectral changes were retained by (ODN1), although the magnitude of the change was diminished. In contrast, a control oligdeoxyribonucleotide (ODN2) having similar base composition did not show any significant change in the CD spectrum in the presence of 150 mM NaCl and 2 mM spermine. The changes in the CD spectrum of ODN1 were highly sensitive to polyamine structure, as evidenced by experiments using spermine analogs with altered number of -CH2- groups separating the amino and imino groups. Electrophoretic mobility shift analysis further showed ODN1 stabilization by spermine and its analogs. These data demonstrate the ability of an ODN containing ERE to undergo conformational transitions in the presence of polyamines and suggest a possible mechanism for polyamine-mediated alterations in the interaction of estrogen receptor with ERE.

Distribution of sequence-dependent curvature in genomic DNA sequences

The distribution of inherent, sequence-dependent curvature was calculated for a number of prokaryotic (M. genitalium, H. influenzae, M. jannaschii), viral (adenovirus 2, equine herpes virus 1), phage (M13, lambda), eukaryotic (S. cerevisiae) and mitochondrial genomes as well as E. coli and human genomic fragments. The genomic averages are in the range of 6-8 degrees/helical turn and only about 20% of DNA is curved less than 3 degrees/helical turn. The prokaryotes and phages appear to have a consistently higher frequency of curved DNA in their genomes than the other genomes tested. Long, highly curved segments, similar to artificially designed curved DNA, are apparently absent from the genomes. Short, curved segments, differing in G+C content may provide environmentally modulated conformational signals for gene regulation. A WWW-server was constructed for the prediction of curved sites from DNA sequences (http://icgeb.trieste.it/dna/curve_it.html/)..

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.

Novel DNA binding specificities of a putative herpesvirus bZIP oncoprotein

Marek's disease virus is a highly oncogenic herpesvirus that can cause T lymphomas and peripheral nerve demyelination in chickens. meq, a candidate oncogene of Marek's disease virus, encodes a basic leucine zipper (bZIP) transcription factor which contains a large proline-rich domain in its C terminus. On the basis of its bZIP structural homology, meq is perhaps the only member of the jun-fos gene family completely viral in origin. We previously showed that Meq's C-terminal domain has potent transactivation activity and that its bZIP domain can dimerize with itself and with c-Jun also. In an effort to identify viral and cellular targets of Meq, we have determined the optimal binding sites for Meq-Jun heterodimers and Meq-Meq homodimers. By a PCR-based approach using cyclic amplification of selected targets, Meq-Jun heterodimers were found to optimally bind tetradecanoylphorbol acetate response element (TRE) and cyclic AMP response element (CRE) consensus sequences. This result was consistent with the results of our previous functional analysis implicating Meq-Jun heterodimers in the transactivation of the Meq promoter through a TRE- or CRE-like sequence. Interestingly, Meq-Meq homodimers were found to bind two distinct motif elements. The first [GAGTGATG AC(G)TCATC] has a consensus which includes a TRE or CRE core flanked by additional nucleotides critical for tight binding. Methylation interference and mutational analyses confirmed the importance of the flanking residues. The sequences of a subset of TRE and CRE sites selected by Meq-Meq are closely related to the binding motif of Maf, another bZIP oncoprotein. The second putative Meq binding site (RACACACAY) bears a completely different consensus not shared by other bZIP proteins. Binding to this consensus sequence also requires secondary structure characteristics associated with DNA bending. CACA motifs are known to promote DNA curvature and function in a number of special biological processes. Our results lend further weight to the increasing importance of DNA bending in transcriptional regulation and provide a baseline for the identification of Meq-responsive targets.

Distribution of bending propensity in DNA sequences

Local bending propensity and curvature of DNA can be characterized using a vector description of DNA bendability, based on a set of parameters derived from deoxyribonuclease I (DNase I) cleavage experiments. Two characteristics-arithmetic and vector averages of bendability-were successfully used to predict experimentally known bendable, rigid and curved segments in DNA. A characteristic distribution of bendability is conserved in evolutionarily related kinetoplast sequences. An analysis of the M. genitalium and H. influenzae genomes as well as fragments of human and yeast genomes shows, on the other hand, that highly curved segments--similar to artificially designed curved oligonucleotides--are extremely rare in natural DNA.

Bending and torsional flexibility of G/C-rich sequences as determined by cyclization assays

The structural polymorphism of DNA is a vital aspect of its biological function. However, it has become increasingly apparent in recent years that DNA polymorphism is a complicated, multidimensional phenomenon that includes not only static sequence-directed structures but dynamic effects as well, including influences of counterions and sequence context. In order to address some of these additional factors that govern DNA conformation, we have used T4 ligase-mediated cyclization to investigate bending in a series of DNA sequences containing the GGGCCC.GGGCCC motif in different sequence contexts including various helical phasings with (A)5-tracts. We present evidence for curvature in GGGCCC.GGGCCC and (A)5-tract motifs in the presence of physiological levels of Mg2+ and show that these motifs curve through similar but oppositely directed bending angles under these ionic strength conditions. Although these two sequence motifs appear to bend similarly, our results suggest significant differences in stiffness and stability of curvature between them. We also show that under the same experimental conditions, the CTAG-CTAG sequence element possesses unusual torsional flexibility and that this appears to be associated with the central TA.TA dinucleotide. The results underscore the need to include sequence context and specific ion effects as well as a dynamic basis in more complete predictive models for functionally related DNA polymorphism.

Sequence determinants of DNA bending in the ilvlH promoter and regulatory region of Escherichia coli

Previous studies have shown that the promoter/regulatory region of the ilvlH operon displays intrinsic curvature, with the bend center located at position -120 relative to the transcription start site. In this report, a 57 bp sequence spanning the bend center was mutagenized in vitro in order to study the relationship between nucleotide sequence and curvature measured by electrophoresis. The strategy used for analyzing the results consisted of determining the strengths of the relationships between electrophoretic anomaly and predicted curvature calculated by computer programs that differ in wedge angle composition. The results revealed that programs which assume that bending occurs only at AA/TT display good predictive value, with correlation coefficients between electrophoretic anomaly and predicted curvature as high as 0.93. In contrast, a program which assumes that bending occurs at all 16 dinucleotide steps exhibited lower predictive value, while there were no significant relationships between the experimental data and curvature calculated by a program that was based on all non-AA/TT wedge values. These results show that the complete wedge model which incorporates values for all dinucleotide steps does not adequately describe the electrophoretic data in this report.

Sequence dependence of base-pair opening in a DNA dodecamer containing the CACA/GTGT sequence motif

Proton nuclear magnetic resonance spectroscopy is used to characterize the kinetics and energetics of base-pair opening in two self-complementary DNA dodecamer duplexes: [d(CGCACATGTGCG)]2 and [d(CGCAGATCTGCG)]2. The first dodecamer contains two symmetrical CACA/GTGT motifs; in the second dodecamer, each motif is interrupted by a change of the central C.G base pair to a G.C base pair. The opening rates and the equilibrium constants for formation of the open state of each base pair are obtained from the dependence of the imino proton exchange rates on the concentration of ammonia catalyst. The results indicate that the opening rates of the central three base pairs in the CACA/GTGT motif are 3-8-fold larger than the corresponding ones in the CAGA/GTCT sequence. The activation enthalpies and entropies, and the standard enthalpy and entropy changes for formation of the open state, are obtained from the temperature dependence of the opening rates and equilibrium constants, respectively. The results reveal that enthalpy/entropy compensation exists, for all base pairs in both dodecamers, in activation as well as in the equilibria between closed and open states. As a result, the opening rates and equilibrium constants for opening are maintained, in both dodecamers, within a relatively narrow range of values. Nevertheless, large sequence-induced variations are observed for the activation enthalpies and the standard enthalpy changes for opening. The A.T base pair located between the C.G base pairs in the CACA/GTGT motif has a negative enthalpy change for formation of the activated state during opening. This is the first case in which a negative activation enthalpy is observed for opening of a Watson-Crick base pair in DNA.(ABSTRACT TRUNCATED AT 250 WORDS)

Structure of three-way DNA junctions. 2. Effects of extra bases and mismatches

The structure of three-way DNA junctions, containing two linear double helices (arms) and a hairpin as a third arm, was studied by means of a cyclization technique. In addition to branched molecules containing perfect base-pairing in helical parts, three-way junctions with mismatches and extra non-complementary nucleotides (bulges) at junction points were studied. Molecules thus designed were ligated at identical conditions and their geometry was compared through the analysis of the efficiency of circle formation. The analysis showed that irregularities in base pairing listed above dramatically change the static and dynamic structural characteristics of the three-way junctions. All mismatches facilitate the kink between linear arms, but quantitatively, the effect depends on the position of the mismatch. The effect is maximal for GG-mismatch placed at the hairpin junction point. The results for bulges are of different kind, and they lead us to conclude that the three-way DNA junction with unpaired nucleotides adopts a T-like geometry with an angle around 90 degrees between arms containing the bulge and two other arms coaxially stacked. Broad distribution of circles indicates that this T-form geometry of bulge-containing junction is more flexible than initial pyramidal structure predominantly due to high mobility of the third arm.

Structure of three-way DNA junctions. 1. Non-planar DNA geometry

Three-way junctions were obtained by annealing two synthetic DNA-oligomers. One of the strands contains a short palindrome sequence, leading to the formation of a hairpin with four base pairs in the stem and four bases in the loop. Another strand is complementary to the linear arms of the first hairpin-containing strand. Both strands were annealed to form a three-way branched structure with sticky ends on the linear arms. The branched molecules were ligated, and the ligation mixture was analysed on a two-dimensional gel in conditions which separated linear and circular molecules. Analysis of 2D-electrophoresis data shows that circular molecules with high mobility are formed. Formation of circular molecules is indicative of bends between linear arms. We estimate the magnitude of the angle between linear arms from the predominant size of the circular molecules formed. When the junction-to-junction distance is 20-21 bp, trimers and tetramers are formed predominately, giving an angle between linear arms as small as 60-90 degrees. Rotation of the hairpin position in the three-way junction allowed us to measure angles between other arms, yielding similar values. These results led us to conclude that the three-way DNA junction possesses a non-planar pyramidal geometry with 60-90 degrees between the arms. Computer modeling of the three-way junction with 60 degrees pyramidal geometry showed a predominantly B-form structure with local distortions at the junction points that diminish towards the ends of the helices. The size distributions of circular molecules are rather broad indicating a dynamic flexibility of three-way DNA junctions.

Studies of DNA bending and flexibility using gel electrophoresis

Gel electrophoretic methods have become established as primary tools in the study and elucidation of sequence-directed curvature both in free DNA and in the operator DNA of several site-specific nucleoprotein complexes. Results using them have been generally consistent with physical methods sensitive to DNA structure and conformation in those instances where direct comparisons can be made, and in a number of cases, gel methods have provided unique information not presently available from other techniques. Two basic strategies have been used: one based upon anomalous gel mobility effects; and a second based upon cyclization properties of curved DNA. Within each of these categories, various approaches have been used, some of which can lead, in favorable cases, to quantitative estimation of bending angles. In this review, the various gel-based methods that have been used to date are critically discussed and the qualitative and quantitative information that can be obtained from them is evaluated. A number of possible structural models for DNA curvature are described and a distinction is drawn between static or fixed bending and bending due to anisotropic flexibility at specific sequence loci. The importance and roles of gel electrophoretic methods in providing experimental approaches to this question are discussed. It is suggested that both static curvature and anisotropic flexibility in operator DNA may provide much of the basis for indirect readout of sequence information by specific site-binding regulatory proteins.