Intrinsically bent DNA in a eukaryotic transcription factor recognition sequence potentiates transcription activation

Deciphering the dynamic code: DNA recognition by transcription factors in the ever-changing genome

Transcription factors (TFs) intricately navigate the vast genomic landscape to locate and bind specific DNA sequences for the regulation of gene expression programs. These interactions occur within a dynamic cellular environment, where both DNA and TF proteins experience continual chemical and structural perturbations, including epigenetic modifications, DNA damage, mechanical stress, and post-translational modifications (PTMs). While many of these factors impact TF-DNA binding interactions, understanding their effects remains challenging and incomplete. This review explores the existing literature on these dynamic changes and their potential impact on TF-DNA interactions.

Small extrachromosomal circular DNAs, microDNA, produce short regulatory RNAs that suppress gene expression independent of canonical promoters

Interest in extrachromosomal circular DNA (eccDNA) molecules has increased recently because of their widespread presence in normal cells across every species ranging from yeast to humans, their increased levels in cancer cells and their overlap with oncogenic and drug-resistant genes. However, the majority of eccDNA (microDNA) in mammalian tissues and cell lines are too small to carry protein coding genes. We have tested functional capabilities of microDNA by creating artificial microDNA molecules mimicking known microDNA sequences and have discovered that they express functional small regulatory RNA including microRNA and novel si-like RNA. MicroDNA are transcribed in vitro and in vivo independent of a canonical promoter sequence. MicroDNA that carry miRNA genes form transcripts that are processed by the endogenous RNA-interference pathway into mature miRNA molecules, which repress a luciferase reporter gene as well as endogenous mRNA targets of the miRNA. Further, microDNA that contain sequences of exons repress the endogenous gene from which the microDNA were derived through the formation of novel si-like RNA. We also show that endogenous microDNA associate with RNA polymerases subunits, POLR2H and POLR3F. Together, these results suggest that microDNA may modulate gene expression through the production of both known and novel regulatory small RNA.

Photochemical analysis of structural transitions in DNA liquid crystals reveals differences in spatial structure of DNA molecules organized in liquid crystalline form

The anisotropic shape of DNA molecules allows them to form lyotropic liquid crystals (LCs) at high concentrations. This liquid crystalline arrangement is also found in vivo (e.g., in bacteriophage capsids, bacteria or human sperm nuclei). However, the role of DNA liquid crystalline organization in living organisms still remains an open question. Here we show that in vitro, the DNA spatial structure is significantly changed in mesophases compared to non-organized DNA molecules. DNA LCs were prepared from pBluescript SK (pBSK) plasmid DNA and investigated by photochemical analysis of structural transitions (PhAST). We reveal significant differences in the probability of UV-induced pyrimidine dimer photoproduct formation at multiple loci on the DNA indicative of changes in major groove architecture.

Molecular organization of the 25S-18S rDNA IGS of Fagus sylvatica and Quercus suber: a comparative analysis

The 35S ribosomal DNA (rDNA) units, repeated in tandem at one or more chromosomal loci, are separated by an intergenic spacer (IGS) containing functional elements involved in the regulation of transcription of downstream rRNA genes. In the present work, we have compared the IGS molecular organizations in two divergent species of Fagaceae, Fagus sylvatica and Quercus suber, aiming to comprehend the evolution of the IGS sequences within the family. Self- and cross-hybridization FISH was done on representative species of the Fagaceae. The IGS length variability and the methylation level of 18 and 25S rRNA genes were assessed in representatives of three genera of this family: Fagus, Quercus and Castanea. The intergenic spacers in Beech and Cork Oak showed similar overall organizations comprising putative functional elements needed for rRNA gene activity and containing a non-transcribed spacer (NTS), a promoter region, and a 5′-external transcribed spacer. In the NTS: the sub-repeats structure in Beech is more organized than in Cork Oak, sharing some short motifs which results in the lowest sequence similarity of the entire IGS; the AT-rich region differed in both spacers by a GC-rich block inserted in Cork Oak. The 5′-ETS is the region with the higher similarity, having nonetheless different lengths. FISH with the NTS-5′-ETS revealed fainter signals in cross-hybridization in agreement with the divergence between genera. The diversity of IGS lengths revealed variants from ∼2 kb in Fagus, and Quercus up to 5.3 kb in Castanea, and a lack of correlation between the number of variants and the number of rDNA loci in several species. Methylation of 25S Bam HI site was confirmed in all species and detected for the first time in the 18S of Q. suber and Q. faginea. These results provide important clues for the evolutionary trends of the rDNA 25S-18S IGS in the Fagaceae family.

Construction of a genome-scale structural map at single-nucleotide resolution

Few methods are available for mapping the local structure of DNA throughout a genome. The hydroxyl radical cleavage pattern is a measure of the local variation in solvent-accessible surface area of duplex DNA, and thus provides information on the local shape and structure of DNA. We report the construction of a relational database, ORChID (OH Radical Cleavage Intensity Database), that contains extensive hydroxyl radical cleavage data produced from two DNA libraries. We have used the ORChID database to develop a set of algorithms that are capable of predicting the hydroxyl radical cleavage pattern of a DNA sequence of essentially any length, to high accuracy. We have used the prediction algorithm to produce a structural map of the 30 Mb of the ENCODE regions of the human genome.

Interaction of estrogen receptor alpha with proliferating cell nuclear antigen

The ability of estrogen receptor alpha (ERalpha) to modulate gene expression is influenced by the recruitment of a host of co-regulatory proteins to target genes. To further understand how estrogen-responsive genes are regulated, we have isolated and identified proteins associated with ERalpha when it is bound to DNA containing the consensus estrogen response element (ERE). One of the proteins identified in this complex, proliferating cell nuclear antigen (PCNA), is required for DNA replication and repair. We show that PCNA interacts with ERalpha in the absence and in the presence of DNA, enhances the interaction of ERalpha with ERE-containing DNA, and associates with endogenous estrogen-responsive genes. Interestingly, rather than altering hormone responsiveness of endogenous, estrogen-responsive genes, PCNA increases the basal expression of these genes. Our studies suggest that in addition to serving as a platform for the recruitment of DNA replication and repair proteins, PCNA may serve as a platform for transcription factors involved in regulating gene expression.

The deoxyribonucleic acid repair protein flap endonuclease-1 modulates estrogen-responsive gene expression

The ligand-occupied estrogen receptor alpha (ERalpha) initiates changes in gene expression through its interaction with target DNA. The capacity of ERalpha to modulate gene expression is influenced by the association of the receptor with a variety of coregulatory proteins. To further understand the role of these coregulatory proteins in ERalpha-mediated transcription, we have isolated and identified proteins associated with ERalpha when it is bound to the consensus estrogen response element. One of the proteins identified in this complex, flap endonuclease-1 (FEN-1), is required for DNA replication and repair. We show that FEN-1 interacts directly with ERalpha and enhances the interaction of ERalpha with estrogen response element-containing DNA. More importantly, chromatin immunoprecipitation and RNA interference assays demonstrate that endogenously expressed FEN-1 associates with the native pS2 gene in MCF-7 cells and influences estrogen-responsive gene expression. Interestingly, estrogen differentially regulates expression of FEN-1 in mouse uterine epithelial, stromal, and myometrial cells. Together, our studies help to elucidate the functional consequence of the ERalpha-FEN-1 interaction and increase our understanding of the elaborate regulatory mechanisms that drive estrogen-responsive gene expression and DNA repair.

Eukaryotic promoter prediction based on relative entropy and positional information

The eukaryotic promoter prediction is one of the most important problems in DNA sequence analysis, but also a very difficult one. Although a number of algorithms have been proposed, their performances are still limited by low sensitivities and high false positives. We present a method for improving the performance of promoter regions prediction. We focus on the selection of most effective features for different functional regions in DNA sequences. Our feature selection algorithm is based on relative entropy or Kullback-Leibler divergence, and a system combined with position-specific information for promoter regions prediction is developed. The results of testing on large genomic sequences and comparisons with the PromoterInspector and Dragon Promoter Finder show that our algorithm is efficient with higher sensitivity and specificity in predicting promoter regions.

Regulation of transcription by synthetic DNA-bending agents

Gene expression is regulated by a complex interplay between binding and the three-dimensional arrangement of transcription factors with RNA polymerase and DNA. Previous studies have supported a direct role for DNA bending and conformation in gene expression, which suggests that agents that induce bends in DNA might be able to control gene expression. To test this hypothesis, we examined the effect of triple-helix-forming oligonucleotide (TFO) bending agents on the transcription of luciferase in an in vitro transcriptional/translational system. We find that transcription is regulated only by a TFO that induces a bend in the DNA. Related TFOs that do not induce bends in DNA have no effect on transcription. Reporter expression can be increased by as much as 80 % or decreased by as much as 50 % depending on the phasing of the upstream bend relative to the promoter. We interpret the results as follows: when the bend is positioned such that the upstream DNA is curved toward the RNA polymerase on the same DNA face, transcription is enhanced. When the upstream DNA is curved away, transcription is attenuated. These results support the hypothesis that DNA-bending agents might have the capability to regulate gene expression, thereby opening up a previously undervalued avenue in research on the artificial control of gene expression.

The importance and identification of regulatory polymorphisms and their mechanisms of action

The search for the genetic variations underlying all human phenotypes is in its infancy but must be one of the long term goals of the scientific community. There is evidence that most, if not all human phenotypes, including illnesses are influenced by the genetic makeup of the individual. There are an estimated 11 million human genetic polymorphisms with a minor allele frequency >1% and possibly many times that number of rare sequence variants. The proportion of these sequence variants which have any functional effect is unknown but it is likely that the majority of those which influence illness lie outside of the amino acid coding regions of genes, and affect the regulation of gene expression--these are called rSNPs. Recent research suggests that about 50% of genes have one or more common rSNPs associated with them and probably most if not all genes have an rSNP within the human population. In the long term, determining which polymorphisms are potentially functional must be done bio-informatically using algorithms based upon experimental data. However, at the current time, the limited data that has been obtained does not allow the creation of such an algorithm. In vitro studies suggest that a large proportion of rSNPs lie within the core and proximal promoter regions of genes but it is not clear how the majority of these influence transcription, as they do not appear to be within any known transcription factor binding sites. However, promoter regions possess a number of sequence-dependent characteristics which make them distinct from the rest of the genome, namely stability, curvature and flexibility. Subtle changes to these features may underlie the mechanisms by which many polymorphisms exert their function.

A hybrid neural network system for prediction and recognition of promoter regions in human genome

This paper proposes a high specificity and sensitivity algorithm called PromPredictor for recognizing promoter regions in the human genome. PromPredictor extracts compositional features and CpG islands information from genomic sequence, feeding these features as input for a hybrid neural network system (HNN) and then applies the HNN for prediction. It combines a novel promoter recognition model, coding theory, feature selection and dimensionality reduction with machine learning algorithm. Evaluation on Human chromosome 22 was approximately 66% in sensitivity and approximately 48% in specificity. Comparison with two other systems revealed that our method had superior sensitivity and specificity in predicting promoter regions. PromPredictor is written in MATLAB and requires Matlab to run. PromPredictor is freely available at http://www.whtelecom.com/Prompredictor.htm.

Structural properties of promoters: similarities and differences between prokaryotes and eukaryotes

During the process of transcription, RNA polymerase can exactly locate a promoter sequence in the complex maze of a genome. Several experimental studies and computational analyses have shown that the promoter sequences apparently possess some special properties, such as unusual DNA structures and low stability, which make them distinct from the rest of the genome. But most of these studies have been carried out on a particular set of promoter sequences or on promoter sequences from similar organisms. To examine whether the promoters from a wide variety of organisms share these special properties, we have carried out an analysis of sets of promoters from bacteria, vertebrates and plants. These promoters were analyzed with respect to the prediction of three different properties, such as DNA curvature, bendability and stability, which are relevant to transcription. All the promoter sequences are predicted to share certain features, such as stability and bendability profiles, but there are significant differences in DNA curvature profiles and nucleotide composition between the different organisms. These similarities and differences are correlated with some of the known facts about transcription process in the promoters from the three groups of organisms.

Left-handedly curved DNA regulates accessibility to cis-DNA elements in chromatin

There is little information on chromatin structure that allows access of trans-acting transcription factors. Logically, the target DNA elements become accessible by either exposing themselves towards the environment on the surface of the nucleosome, or making the regulatory region free of the nucleosome. Here, we demonstrate that curved DNA that mimics a negative supercoil can play both roles in the promoter region. By constructing 35 reporter plasmids and using in vivo assay systems, we scrutinized the relationships between upstream DNA geometry, nucleosome positioning and promoter activity. When the left-handedly curved DNA was linked to the herpes simplex virus thymidine kinase (HSV tk) promoter at a specific rotational phase and distance, the curved DNA attracted the nucleosome and the TATA box was thereby left in the linker DNA with its minor groove facing outwards, which led to the activation of transcription. Neither planar curving, nor right-handedly curved DNA nor straight DNA had this effect. Our results seem to provide a clue for solving the problem of why curved DNA is often located near transcriptional control regions.

Differential regulation of the human progesterone receptor gene through an estrogen response element half site and Sp1 sites

The progesterone receptor (PR) gene is regulated by estrogen in normal reproductive tissues and in MCF-7 human breast cancer cells. Although it is generally thought that estrogen responsiveness is mediated by interaction of the ligand-occupied estrogen receptor (ER) with estrogen response elements (EREs) in target genes, the human progesterone receptor (PR) gene lacks a palindromic ERE. Promoter A of the PR gene does, however, contain an ERE half site upstream of two adjacent Sp1 sites from +571 to +595, the +571 ERE/Sp1 site. We have examined the individual contributions of the ERE half site and the two Sp1 sites in regulating estrogen responsiveness. Transient transfection assays demonstrated that both Sp1 sites were critical for estrogen-mediated activation of the PR gene. Interestingly, rather than decreasing transcription, mutations in the ERE half site increased transcription substantially suggesting that this site plays a role in limiting transcription. Chromatin immunoprecipitation assays demonstrated that Sp1 was associated with the +571 ERE/Sp1 site in the endogenous PR gene in the absence and in the presence of estrogen, but that ERalpha was only associated with this region of the PR gene after MCF-7 cells had been treated with estrogen. Our studies provide evidence that effective regulation of transcription through the +571 ERE/Sp1 site requires the binding of ERalpha and Sp1 to their respective cis elements and the appropriate interaction of ERalpha and Sp1 with other coregulatory proteins and transcription factors.

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.

HMG proteins and DNA flexibility in transcription activation

The relative stiffness of naked DNA is evident from measured values of longitudinal persistence length (approximately 150 bp) and torsional persistence length (approximately 180 bp). These parameters predict that certain arrangements of eukaryotic transcription activator proteins in gene promoters should be much more effective than others in fostering protein-protein interactions with the basal RNA polymerase II transcription apparatus. Thus, if such interactions require some kind of DNA looping, DNA loop energies should depend sensitively on helical phasing of protein binding sites, loop size, and intrinsic DNA curvature within the loop. Using families of artificial transcription templates where these parameters were varied, we were surprised to find that the degree of transcription activation by arrays of Gal4-VP1 transcription activators in HeLa cell nuclear extract was sensitive only to the linear distance separating a basal promoter from an array of bound activators on DNA templates. We now examine the hypothesis that this unexpected result is due to factors in the extract that act to enhance apparent DNA flexibility. We demonstrate that HeLa cell nuclear extract is rich in a heat-resistant activity that dramatically enhances apparent DNA longitudinal and torsional flexibility. Recombinant mammalian high-mobility group 2 (HMG-2) protein can substitute for this activity. We propose that the abundance of HMG proteins in eukaryotic nuclei provides an environment in which DNA is made sufficiently flexible to remove many constraints on protein binding site arrangements that would otherwise limit efficient transcription activation to certain promoter geometries.

Intrinsic DNA bends: an organizer of local chromatin structure for transcription

DNA with a curved trajectory of its helix axis is called bent DNA, or curved DNA. Interestingly, biologically important DNA regions often contain this structure, irrespective of the origin of DNA. In the last decade, considerable progress has been made in clarifying one role of bent DNA in prokaryotic transcription and its mechanism of action. However, the role of bent DNA in eukaryotic transcription remains unclear. Our recent study raises the possibility that bent DNA is implicated in the "functional packaging" of transcriptional regulatory regions into chromatin. In this article, I review recent progress in bent DNA research in eukaryotic transcription, and summarize the history of bent DNA research and several subjects relevant to this theme. Finally, I propose a hypothesis that bent DNA structures that mimic a negative supercoil, or have a right-handed superhelical writhe, organize local chromatin infrastructure to help the very first interaction between cis-acting DNA elements and activators that trigger transcription.

Intrinsically bent DNA in the promoter regions of the yeast GAAL1-10 and GAL80 genes

Circular permutation analysis has detected fairly strong sites of intrinsic DNA bending on the promoter regions of the yeast GAL1-10 and GAL80 genes. These bends lie in functionally suggestive locations. On the promoter of the GAL1-10 structural genes, strong bends bracket nucleosome B, which lies between the UAS(G) and the GAL1 TATA. These intrinsic bends could help position nucleosome B. Nucleosome B plus two other promoter nucleosomes protect the TATA and start site elements in the inactive state of expression but are completely disrupted (removed) when GAL1-10 expression is induced. The strongest intrinsic bend ( approximately 70 degrees ) lies at the downstream edge of nucleosome B; this places it approximately 30 base pairs upstream of the GAL1 TATA, a position that could allow it to be involved in GAL1 activation in several ways, including the recruitment of a yeast HMG protein that is required for the normally robust level of GAL1 expression in the induced state (Paull, T., Carey, M., and Johnson, R. (1996) Genes Dev. 10, 2769-2781). On the regulatory gene GAL80, the single bend lies in the non-nucleosomal hypersensitive region, between a GAL80-specific far upstream promoter element and the more gene-proximal promoter elements. GAL80 promoter region nucleosomes contain no intrinsically bent DNA.

A-Tract bending: insights into experimental structures by computational models

While solution structures of adenine tract (A-tract) oligomers have indicated a unique bend direction equivalent to negative global roll (commonly termed "minor-groove bending"), crystallographic data have not unambiguously characterized the bend direction; nevertheless, many features are shared by all A-tract crystal and solution structures (e.g. propeller twisting, narrow minor grooves, and localized water spines). To examine the origin of bending and to relate findings to the crystallographic and solution data, we analyze molecular dynamics trajectories of two solvated A-tract dodecamers: 1D89, d(CGCGA(6)CG), and 1D98, d(CGCA(6)GCG), using a new general global bending framework for analyzing bent DNA and DNA/protein complexes. It is significant that the crystallographically-based initial structures are converted from dissimilar to similar bend directions equivalent to negative global roll, with the average helical-axis bend ranging from 10.5 degrees to 14.1 degrees. The largest bend occurs as positive roll of 12 degrees on the 5' side of the A-tracts (supporting a junction model) and is reinforced by gradual curvature at each A-tract base-pair (bp) step (supporting a wedge model). The precise magnitude of the bend is subtly sequence dependent (consistent with a curved general sequence model). The conversion to negative global roll only requires small local changes at each bp, accumulated over flexible moieties both outside and inside the A-tract. In contrast, the control sequence 1BNA, d(CGCGA(2)TTCGCG), bends marginally (only 6.9 degrees ) with no preferred direction. The molecular features that stabilize the bend direction in the A-tract dodecamers include propeller twisting of AT base-pairs, puckering differences between A and T deoxyriboses, a narrow minor groove, and a stable water spine (that extends slightly beyond the A-tract, with lifetimes approaching 0.2 ns). The sugar conformations, in particular, are proposed as important factors that support bent DNA. It is significant that all these curvature-stabilizing features are also observed in the crystallographic structures, but yield overall different bending paths, largely due to the effects of sequences outside the A-tract. These results merge structural details reported for A-tract structures by experiment and theory and lead to structural and dynamic insights into sequence-dependent DNA flexibility, as highlighted by the effect of an A-tract variant of a TATA-box element on bending and flexibility required for TBP binding.

Characterization of the selective staining of DNA on gels using two fluorochromes

Ethidium bromide (EB) and 4'-6-diamidino-2-phenylindole (DAPI) are both well-known fluorochromes for detecting DNA fragments. EB binds to DNA by intercalation and DAPI binds in the DNA minor groove. We previously developed a staining method using both EB and DAPI that is selective for AT-rich DNA fragments. Using this double-staining method, AT-rich DNA fragments are visualized as bluish-white fluorescent bands. To further characterize this method, a series of synthetic DNA fragments were designed with systematic variations in the length, AT content, and DNA sequence pattern. The staining properties of these fragments were determined in the presence of DAPI and EB, and the following results were obtained. (i) In a series of fragments with three AT base pairs followed by one GC base pair, the stained DNA fragments exhibited different fluorescent colors and varied from bluish (more DAPI staining) to pinkish (less DAPI staining) in the order 5'-AAA-3', 5'-AAT-3', 5'-ATA-3', 5'-TTA-3'. (ii) In fragments with constant AT content, the blue fluorescent color increased with increasing number of A (or T) nucleotides, due to increased DAPI binding. The blue color was saturated when the number of A (or T) nucleotides was 12 or greater. (iii) The fluorescent color of the stained DNA fragments changed in the order of red-orange, pink, pinkish-white, white, bluish-white, blue as the AT content increased from 0 to 100%. Thus, the fluorescent color of DNA fragments stained with DAPI and EB depends on base composition and nucleotide sequence, suggesting that individual stained DNA fragments may have characteristic and specific fluorescent colors. The fluorescent color emitted by specific stained DNA fragments in the presence of EB and DAPI can be analyzed with a high degree of sensitivity and resolution using the XYZ colorimetric system.

Structural properties of hybrid triplex of polycation deoxyribonucleic S-methylthiourea (DNmt) strands with a complementary DNA strand, probed by nanosecond molecular dynamics

The replacement of phosphodiester linkages of the polyanion DNA with S-methylthiourea linkers provides the polycation deoxyribonucleic S-methylthiourea (DNmt). Molecular dynamics studies to 1,220 ps of the hybrid triplex formed from octameric DNmt strands d(Tmt)8 with a complementary DNA oligomer strand d(Ap)8 have been carried out with explicit water solvent and Na+Cl- counterions under periodic boundary conditions using the CHARMM force field and the Ewald summation method. The Watson-Crick and Hoogsteen hydrogen-bonding patterns of the A/T tracts remained intact without any structural restraints for triplex structures throughout the simulation. The duplex portion of the triplex structure equilibrated at a B-DNA conformation in terms of the helical rise and other helical parameters. The dynamic structures of the DNmt x DNA x DNmt triplex were determined by examining histograms from the last 800 ps of the dynamics run. These included the hydrogen-bonding pattern (sequence recognition), three-centered bifurcating occurrences, minor groove width variations, and bending of tracts for the hybrid triplex structures. Together with the Watson-Crick hydrogen-bondings, the strong Hoogsteen hydrogen-bondings, the partially maintained three-centered bifurcatings in the Watson-Crick pair, and the medium-strength three-centered bifurcatings in the Hoogsteen pair suggest that the hybrid triplex is energetically favorable as compared to a duplex with similar base stacking, van der Waals interactions, and helical parameters. This is in agreement with our previously reported thermodynamic study, in which only triplex structures were observed in solution. The bending angle measured between the local axis vectors of the first and last helical axis segments is about 20 degrees for the Watson-Crick portion of the averaged structure. Propeller twist (associated with three-centered hydrogen-bonding) up to -30 degrees, native to DNA AT base pairing, was also observed for the triplex structure. The sugar pseudorotation phase angles and the ring rotation angles for the DNA strand are within the C3'-endo domain and C2'-endo domain for the DNmt strand. Water spines are observed in both minor and major grooves throughout the dynamics run. The molecular dynamics simulations of the structural properties of DNmt x DNA x DNmt hybrid triplex is compared to the DNG x DNA x DNG hybrid triplex (In DNG the -O-(PO2-)-O- linkers in DNA is replaced by -NH-C(=N+H2)-NH-).

Promoter sequences and algorithmical methods for identifying them

This paper presents a survey of currently available mathematical models and algorithmical methods for trying to identify promoter sequences. The methods concern both searching in a genome for a previously defined consensus and extracting a consensus from a set of sequences. Such methods were often tailored for either eukaryotes or prokaryotes although this does not preclude use of the same method for both types of organisms. The survey therefore covers all methods; however, emphasis is placed on prokaryotic promoter sequence identification. Illustrative applications of the main extracting algorithms are given for three bacteria.

DNA damage by nitrogen mustard in a gene containing multiple Sp1-binding sites

The human cytochrome c(1) gene TATA-less promoter contains 10 Sp1-binding elements that regulate the activation of transcription of this gene. Quantitative PCR was used to show that nitrogen mustard induces DNA lesions within this Sp1-binding region following exposure of HeLa cells to clinical levels of the drug. Alkylation of the cytochrome c(1) gene in HeLa cells increased with reaction time up to 4 h following exposure to nitrogen mustard, with 50% of the lesions (approximately 0.8/kb) forming within 1 h. An Sp1 competition assay showed that nitrogen mustard inhibited the binding of Sp1 to the promoter region of the cytochrome c(1) gene in HeLa cells. These results show that nitrogen mustard-induced damage to Sp1-binding sites may contribute to the toxicity of this compound by interfering with the activation of specific genes.

Promoter upstream bent DNA activates the transcription of the Clostridium perfringens phospholipase C gene in a low temperature-dependent manner

The phospholipase C gene (plc) of Clostridium perfringens possesses three phased A-tracts forming bent DNA upstream of the promoter. An in vitro transcription assay involving C.perfringens RNA polymerase (RNAP) showed that the phased A-tracts have a stimulatory effect on the plc promoter, and that the effect is proportional to the number of A-tracts, and more prominent at lower temperature. A gel retardation assay and hydroxyl radical footprinting revealed that the phased A-tracts facilitate the formation of the RNAP-plc promoter complex through extension of the contact region. The upstream (UP) element of the Escherichia coli rrnB P1 promoter stimulated the downstream promoter activity temperature independently, differing from the phased A-tracts. When the UP element was placed upstream of the plc promoter, low temperature-dependent stimulation was observed, although this effect was less prominent than that of the phased A-tracts. These results suggest that both the phased A-tracts and UP element cause low temperature-dependent activation of the plc promoter through a similar mechanism, and that the more efficient low temperature-dependent activation by the phased A-tracts may be due to an increase in the bending angle at a lower temperature.

Molecular dynamics simulations of an oligonucleotide duplex with adenine tracts phased by a full helix turn

A theoretical model of a DNA oligonucleotide duplex featuring A-tracts phased by a full helix turn is developed based on molecular dynamics computer simulation. The extent to which this model agrees with relevant experimental data on axis bending and the relationship of A-tracts to bending and other aspects of helix morphology is investigated. Specifically, a series of nanosecond-level molecular dynamics (MD) simulations have been carried out for the 25 bp duplex d(ATAGGCAAAAAATAGGCAAAAATGG) at various concentrations of saline solution. A 30 base-pair sequence composed of three 10 bp repeats of the BamHI recognition sequence ligated together, d(CGGGATCCCG. CGGGATCCCG.CGGGATCCCG), was simulated as a control. The MD was carried out using the AMBER 4.1 suite of programs, and utilized the Cornell et al. force-field with the electrostatic boundary conditions treated by the particle-mesh Ewald summation protocol. The MD results show that at a concentration of 60 mM KCl, 10 mM MgCl2 added salt plus minimal neutralizing cations, the MD model exhibits concerted axis bending to the extent of 15.5 degrees per A-tract. This compares favorably with the bending per turn of 17 to 21 degrees inferred from cyclization experiments. The MD model also exhibits a progressive 5' to 3' narrowing of the minor-groove region of A-tracts, a feature inferred from DNA footprinting experiments. Analysis of the dynamic structure of the MD models shows that the origin of the bending follows a junction-type bending model with an admixture of mixed sequence effects, with A-tracts relatively straight, as in oligonucleotide crystal structures of sequences containing A-tracts. The results are shown to be sensitive to environmental conditions: MD on d(ATAGGCAAAAAATAGGCAAAAATGG) in neutralizing Na+ buffer results in markedly reduced curvature, and the removal of Mg2+ measurably affects bending. Carrying out the simulations at experimental salt conditions appears to be essential to obtain an accurate account of the experimentally observed bending.

The ligand insertion hypothesis in the genomic action of steroid hormones

Gene regulation by steroids is tightly coupled to hormone concentration and stereochemistry. A key step is binding of hormones to receptors which interact with consensus DNA sequences known as hormone response elements (HREs). The specificity and strength of hormone binding do not correlate well with hormonal activity suggesting an additional step involving recognition of ligand by the gene. Stereospecific fit of hormones between base pairs and correlation of fit with hormonal activity led to the proposal that such recognition involves insertion of hormone into DNA. Here, the feasibility of insertion was investigated using computer models of the glucocorticoid receptor DNA binding domain bound to its HRE. The site reported to accommodate glucocorticoids was found in the HRE and was exposed to permit unwinding at this locus. The resulting cavity in the unwound DNA/receptor interface fit cortisol remarkably well; cortisol formed hydrogen bonds to both the receptor and DNA. Current experimental evidence is generally consistent with ligand binding domains of receptors undergoing a conformational change which facilitates transfer of the ligand into the unwound DNA/receptor interface. We propose this step is rate limiting and alterations in receptor, DNA or hormone which attenuate insertion impair hormonal regulation of gene function.

Prebending the estrogen response element destabilizes binding of the estrogen receptor DNA binding domain

Binding of many eukaryotic transcription regulatory proteins to their DNA recognition sequences results in conformational changes in DNA. To test the effect of altering DNA topology by prebending a transcription factor binding site, we examined the interaction of the estrogen receptor (ER) DNA binding domain (DBD) with prebent estrogen response elements (EREs). When the ERE in minicircle DNA was prebent toward the major groove, which is in the same direction as the ER-induced DNA bend, there was no significant effect on ER DBD binding relative to the linear counterparts. However, when the ERE was bent toward the minor groove, in a direction that opposes the ER-induced DNA bend, there was a four- to eightfold reduction in ER DBD binding. Since reduced binding was also observed with the ERE in nicked circles, the reduction in binding was not due to torsional force induced by binding of ER DBD to the prebent ERE in covalently closed minicircles. To determine the mechanism responsible for reduced binding to the prebent ERE, we examined the effect of prebending the ERE on the association and dissociation of the ER DBD. Binding of the ER DBD to ERE-containing minicircles was rapid when the EREs were prebent toward either the major or minor groove of the DNA (k(on) of 9.9 x 10(6) to 1.7 x 10(7) M(-1) s(-1)). Prebending the ERE toward the minor groove resulted in an increase in k(off) of four- to fivefold. Increased dissociation of the ER DBD from the ERE is, therefore, the major factor responsible for reduced binding of the ER DBD to an ERE prebent toward the minor groove. These data provide the first direct demonstration that the interaction of a eukaryotic transcription factor with its recognition sequence can be strongly influenced by altering DNA topology through prebending the DNA.

Supercoiling-dependent flexibility of adenosine-tract-containing DNA detected by a topological method

Intrinsically bent DNA sequences have been implicated in the activation of transcription by promoting juxtaposition of DNA sequences near the terminal loop of a superhelical domain. We have developed a novel topological assay for DNA looping based on lambda integrative recombination to study the effects of intrinsically bent DNA sequences on the tertiary structure of negatively supercoiled DNA. Remarkably, the localization of adenosine-tract (A-tract) sequences in the terminal loop of a supercoiled plasmid is independent of the extent of intrinsic bending. The results suggest that A-tract-containing sequences have other properties that organize the structure of superhelical domains apart from intrinsic bending and may explain the lack of conservation in the degree of A-tract-dependent bending among DNA sequences located upstream of bacterial promoters.

The transcription promoter of the spliced leader gene from Trypanosoma cruzi

A putative promoter element responsible for transcription of the spliced leader (SL) gene of Trypanosoma cruzi was identified by overlapping deletion and linker scanning analyses of the upstream flanking sequences using the bacterial chloramphenicol acetyltransferase (CAT) gene as a reporter in transient transfections of cultured epimastigotes. Deletion or substitution of a proximal sequence element (PSE) between positions -53 and -40 relative to the transcription start point eliminated CAT gene expression. Comparison of SL genes from several strains of T. cruzi revealed two alternative sequence patterns for the putative SL PSE, both composed of a short run of purines followed by a run of pyrimidines. Moreover, an examination of these sequences supports the subdivision of T. cruzi isolates into two divergent groups. Double-stranded oligonucleotides containing the sequence of the PSE exhibited specific gel mobility shifts after incubation with T. cruzi nuclear extracts, suggesting that a transcription factor binds this site. Finally, experiments designed to increase the level of CAT expression from the SL promoter suggest that it is not a strong promoter in cultured T. cruzi epimastigotes.

Distamycin prolongs E-selectin expression by interacting with a specific NF-kappaB-HMG-I(Y) binding site in the promoter

The E-selectin cell adhesion protein plays a critical role in mediating adherence of leukocytes to endothelium at sites of inflammation. Cytokine-induced E-selectin expression on the surface of endothelial cells is transient; mRNA expression peaks at 3-4 h after induction and returns to basal levels within 24 h. The mechanism for this transcriptional down-modulation is not known. Promoter binding factors responsible for induced gene expression include NF-kappaB, which binds at three sites within the E-selectin promoter, and HMG-I(Y), which binds to the A/T-rich core found at the centre of these binding sites. Distamycin is an antibiotic that also binds A/T-rich DNA and inhibits HMG-I(Y) DNA binding. To study the role of HMG-I(Y) in E-selectin expression, we have examined the effect of distamycin on the cytokine-induced E-selectin expression cycle. We found that distamycin prolonged E-selectin expression, both by sustaining mRNA transcription and by extending the transcript's half-life. The distamycin effect on transcription was mediated through one of the three NF-kappaB-HMG-I(Y) binding sites (NF-kappaBII) within the promoter. This suggests that the NF-kappaB-HMG-I(Y) complex interacting at the NF-kappaBII site plays a role not only in cytokine induction of E-selectin expression, but also in its down-modulation.

DNA bending is induced by binding of the glucocorticoid receptor DNA binding domain and progesterone receptors to their response element

Circular permutation analysis was used to determine the degree of DNA bending induced by binding of the glucocorticoid receptor (GR) DNA binding domain (DBD), the human progesterone receptor (PR) DBD, PR-A:A and PR-B:B homodimers, and PR-A:B heterodimers to the glucocorticoid response element/progesterone response element (GRE/PRE). The bending angles induced by the GR DBD and the PR DBD were approximately 28 degrees and 25 degrees, respectively. The PR-B:B and PR-A:A homodimers and the PR-A:B heterodimers all induced similar DNA bending angles of 72-77 degrees. The substantially greater DNA bend induced by full-length PR compared to the PR DBD indicates that sequences outside the classic zinc finger DNA binding domain may play an important role in the interaction of PR with the GRE/PRE. Because PR-A:A and PR-B:B homodimers and the PR-A:B heterodimers induce similar DNA bends, the different abilities of the PR-A and PR-B isoforms to activate transcription are not due to differences in their abilities to distort DNA structure.

Characterization of the cAMP response element of the cystic fibrosis transmembrane conductance regulator gene promoter

A dominant negative inhibitor of the cAMP-dependent protein kinase has been shown to inhibit the basal expression of the cystic fibrosis transmembrane conductance regulator (CFTR) gene in the human colon carcinoma cell line, T84. A functional cAMP response element (CRE) was localized at -48 in the CFTR promoter, and we have analyzed the interactions of this regulatory region with transcription factors. An adjacent inverted CCAAT element (Y box) at position -60 was also investigated. Mutation of the CRE or the Y box decreases the activity of the promoter in transient transfections of T84 or JEG-3 cells. Electrophoretic mobility shift assays demonstrate that CRE-binding protein (CREB) binds to the CFTR CRE with high affinity and independently of the adjacent Y box and that the CFTR CRE binds CREB and activating transcription factor-1 in nuclear extracts of T84 and CaLu-3 cells. In transient transfections of JEG-3 cells, activation of the CFTR promoter is blocked by a dominant negative CREB mutant. The CFTR CRE will also drive cAMP-mediated expression when placed upstream of a heterologous basal promoter. These results demonstrate that CFTR is a bona fide CRE-dependent gene, and we suggest that CFTR expression levels in vivo may be responsive to hormones or drugs that activate the cAMP-dependent protein kinase system.

In simple synthetic promoters YY1-induced DNA bending is important in transcription activation and repression

Depending on promoter context, YY1 can activate or repress transcription, or provide a site for transcription initiation. To investigate whether the ability of YY1 to induce DNA bending influenced its ability to activate and repress transcription, simple synthetic promoters were constructed in which the YY1 binding site was inserted between the TATA box and either the NF1 or AP1 recognition sequences. In transient transfections of COS cells, the NF1YY1TATA and NF1RYY1TATA promoters exhibited a dramatic 15-20-fold increase in correctly initiated transcription. These promoters exhibited even larger 60-80-fold increases in transcription in HeLa cells. Neither multiple copies of the YY1 binding site alone, nor placement of a YY1 site upstream of the NF1 site activated transcription. Deletion of 4 bp between the NF1 and YY1 sites, which changes the phase of the DNA bends, abolished the 16-fold activation of transcription by NF1YY1TATA. Insertion of the YY1 site between the AP1 site and the TATA box decreased transcription approximately 3-fold. Replacing the YY1 binding site with an intrinsic DNA bending sequence mimicked this transcription repression. Sequences of similar length which do not bend DNA fail to repress AP1-mediated transcription. Gel mobility shift assays were used to show that binding of YY1 to its recognition sequence did not repress binding of AP1 to its recognition sequences. Our data indicate that YY1-induced DNA bending may activate and repress transcription by changing the spatial relationships between transcription activators and components of the basal transcription apparatus.

Multiple factors regulating the expression of human thromboxane synthase gene

Characterization of the 5.5 kb promoter of human thromboxane synthase (TS) gene revealed a proximal positive regulatory sequence (PPRS, -90 to -25 bp) and several distal repressive elements. The maximal promoter activity was found to reside within the first 285 bp, approximately 75% of which was contributed by the PPRS. The sequence between -365 and -665 bp exerted a strong repressive effect (approximately 55%) on reporter gene expression independent of orientation and position, consistent with properties expected for a silencer. The sequence upstream of -665 bp to -5.5 kb contains mainly repressive elements which further reduce the promoter activity by 30%. The 65 bp PPRS worked in an orientation-independent, but position-dependent, manner and could be further divided into two independent elements, PPRS1 (-90 to -50 bp) and PPRS2 (-50 to -25 bp). While similar nuclear factor(s) from different cell types interact with PPRS2, those interacting with PPRS1 exhibit cell specificity. Internal sequence deletion and oligonucleotide competition established that a binding sequence for NF-E2 in PPRS1 (-60 tgctgattcat -50) was important for enhancing TS promoter activity in HL-60 cells. The presence of NF-E2 mRNA in HL-60 cells was demonstrated by reverse-transcription PCR amplification of the cDNA and Northern blot analysis. A 9-fold transactivation of luciferase (luc) reporter gene expression had been detected when NF-E2 cDNA was co-expressed with a TS promoter/luc construct. Despite the fact that NF-E2 and the cis-elements could alter the efficiency of TS transcription, they were not sufficient for restricting cell-specific TS expression. Analysis of the methylation status at the TS promoter in several human cell lines reveals cell-specific patterns of methylation that might correlate with TS expression. Taken together, these results suggest that the expression of human TS gene is modulated by multiple factors including cis-elements, trans-activator(s), and possibly genomic methylation.

Design of artificial sequence-specific DNA bending ligands

Proteins that bend DNA are important regulators of biological processes. Sequence-specific DNA bending ligands have been designed that bind two noncontiguous sites in the major groove and induce a bend in the DNA. An oligonucleotide containing pyrimidine segments separated by a central variable linker domain simultaneously binds by triple helix formation two 15-bp purine tracts separated by 10 bp. Bend angles of 61 degrees, 50 degrees, and 38 degrees directed towards the minor groove were quantitated by phasing analysis for linkers of four, five, and six T residues, respectively. The design and synthesis of nonnatural architectural factors may provide a new class of reagents for use in biology and human medicine.

The cytoplasmic and N-terminal transmembrane domains of cytochrome P450 contain independent signals for retention in the endoplasmic reticulum

Microsomal cytochrome P450 is inserted into the membrane of the endoplasmic reticulum (ER) by its N-terminal signal/anchor sequence which also functions as an ER retention signal. To analyze further potential retention signals of cytochrome P450, topological domains of cytochrome P450 2C1 or 2C2, epidermal growth factor receptor, a plasma membrane protein, and bacterial alkaline phosphatase, a secreted protein were exchanged. The N-terminal signal/anchor of cytochrome P450 2C1 functioned as an ER retention signal when placed at the N terminus of several reporter proteins but not when fused at the C terminus of the extracellular domain of epidermal growth factor receptor, with or without a heterologous cytoplasmic domain. Chimeric proteins in which the cytoplasmic domain of cytochrome P450 2C2 was substituted for that of epidermal growth factor receptor were retained in the ER indicating that an independent retention signal is present in the cytoplasmic part of cytochrome P450 2C2. These chimeras were enzymatically active which argues against misfolding as the primary cause of retention. The ER retention signal of the cytoplasmic domain could not be localized to a single amino acid segment by deletion analysis. These results show that cytochrome P450 2C2 contains redundant, complex ER retention signals in its cytoplasmic and N-terminal hydrophobic domains and that the function of the N-terminal signal is context-dependent.