Complex Promoter Structure Research Articles

Exact burst-size distributions for gene-expression models with complex promoter structure

In prokaryotic and eukaryotic cells, most genes are transcribed in a bursty fashion on one hand and complex gene regulations may lead to complex promoter structure on the other hand. This raises an unsolved issue: how does promoter structure shape transcriptional bursting kinetics characterized by burst size and frequency? Here we analyze stochastic models of gene transcription, which consider complex regulatory mechanisms. Notably, we develop an efficient method to derive exact burst-size distributions. The analytical results show that if the promoter of a gene contains only one active state, the burst size indeed follows a geometric distribution, in agreement with the previous result derived under certain limiting conditions. However, if it contains a multitude of active states, the burst size in general obeys a non-geometric distribution, which is a linearly weighted sum of geometric distributions. This superposition principle reveals the essential feature of bursting kinetics in complex cases of transcriptional regulation although it seems that there has been no direct experimental confirmation. The derived burst-size distributions not only highlight the importance of promoter structure in regulating bursting kinetics, but can be also used in the exact inference of this kinetics based on experimental data.

Modeling and analytical solution of free energy of complex promoter structure

Based on the biological model of complex promoter structure, this paper establishes a mathematical model, and maps the complex topological structures of promoter into a stochastic transcriptional multi-state model. And in this paper, from the new perspective of non-equilibrium thermodynamics. we combine discretization with continuity and calculate the free energy consumption of the complex promoter. The analytical solution of the free energy consumption of the promoter is obtained, it provides convenience for analyzing the free energy cost of the promoter.

Calculation of Free Energy Consumption in Gene Transcription with Complex Promoter Structure

From the viewpoint of thermodynamics, gene transcription necessarily consumes free energy due to nonequilibrium processes. On the other hand, regulatory molecules present on the core promoter of a gene interact often in a dynamic, highly combinatorial, and possibly energy-dependent manner, leading to a complex promoter structure. This raises the question of how gene transcription with general promoter topology consumes free energy. We propose a biophysically intuitive approach to calculate energy consumption (quantified by the production rate of entropy) of a gene transcription process. Then, we show that the numbers of the ON and OFF states of a promoter can reduce energy consumption of the gene system and the Fano factor of mRNA, and in contrast to other regulatory ways, the cooperative binding of transcription factors to DNA sites always reduces energy consumption but amplifies the mRNA noise. While our proposed approach is general, our obtained qualitative results can in turn be used to the inference of complex promoter structure.

Influence of Complex Promoter Structure on Gene Expression

A gene is often regulated by a variety of transcription factors, leading to complex promoter structure. However, how this structure affects gene expression remains elusive. Here, this paper studies a stochastic gene model with the promoter containing arbitrarily many active and inactive states. First, the authors use the binomial moment method to derive analytical steady-state distributions of the mRNA and protein numbers. Then, the authors analytically investigate how the promoter structure impacts the mean expression levels and the expression noise. Third, numerical simulation finds interesting phenomena, e.g., the common on-off model overestimates the expression noise in contrast to multiple-state models; the multi-on mechanism can reduce the expression noise more than the multi-off mechanism if the mean expression level is kept the same; and multiple exits of transcription can result in multimodal distributions.

Energy Cost and Mean Dwell Times for the Activity of Promoter with Complex Structure

Regulatory molecules present on the core promoter of a gene interact often in a dynamic, highly combinatorial and possibly energy-dependent manner, leading to complex promoter structure and even complex global dynamics. The authors analyze dynamics of an arbitrarily complex promoter from the view of thermodynamics combined with statistic physics. First, the authors formulize transcription factors-mediated promoter kinetics in terms of energy. Then, the authors analyze energetic cost in several representative cases of promoter structure, deriving useful analytical results. Third, the authors derive analytical expressions for mean dwell times of the promoter activity states, experimentally measurable quantities related to the energy cost of promoter dynamics. The overall framework lays a theoretical foundation for analysis of complex promoter kinetics and gene expression dynamics.

Toward predictive R-loop computational biology: genome-scale prediction of R-loops reveals their association with complex promoter structures, G-quadruplexes and transcriptionally active enhancers.

R-loops are three-stranded RNA:DNA hybrid structures essential for many normal and pathobiological processes. Previously, we generated a quantitative R-loop forming sequence (RLFS) model, quantitative model of R-loop-forming sequences (QmRLFS) and predicted ∼660 000 RLFSs; most of them located in genes and gene-flanking regions, G-rich regions and disease-associated genomic loci in the human genome. Here, we conducted a comprehensive comparative analysis of these RLFSs using experimental data and demonstrated the high performance of QmRLFS predictions on the nucleotide and genome scales. The preferential co-localization of RLFS with promoters, U1 splice sites, gene ends, enhancers and non-B DNA structures, such as G-quadruplexes, provides evidence for the mechanical linkage between DNA tertiary structures, transcription initiation and R-loops in critical regulatory genome regions. We introduced and characterized an abundant class of reverse-forward RLFS clusters highly enriched in non-B DNA structures, which localized to promoters, gene ends and enhancers. The RLFS co-localization with promoters and transcriptionally active enhancers suggested new models for in cis and in trans regulation by RNA:DNA hybrids of transcription initiation and formation of 3D-chromatin loops. Overall, this study provides a rationale for the discovery and characterization of the non-B DNA regulatory structures involved in the formation of the RNA:DNA interactome as the basis for an emerging quantitative R-loop biology and pathobiology.

Toward predictive R-loop computational biology: genome-scale prediction of R-loops reveals their association with complex promoter structures, G-quadruplexes and transcriptionally active enhancers

R-loops are three-stranded RNA:DNA hybrid structures essential for many normal and pathobiological processes. Previously, we generated a quantitative R-loop forming sequence (RLFS) model, quantitative model of R-loop-forming sequences (QmRLFS) and predicted ∼660 000 RLFSs; most of them located in genes and gene-flanking regions, G-rich regions and disease-associated genomic loci in the human genome. Here, we conducted a comprehensive comparative analysis of these RLFSs using experimental data and demonstrated the high performance of QmRLFS predictions on the nucleotide and genome scales. The preferential co-localization of RLFS with promoters, U1 splice sites, gene ends, enhancers and non-B DNA structures, such as G-quadruplexes, provides evidence for the mechanical linkage between DNA tertiary structures, transcription initiation and R-loops in critical regulatory genome regions. We introduced and characterized an abundant class of reverse-forward RLFS clusters highly enriched in non-B DNA structures, which localized to promoters, gene ends and enhancers. The RLFS co-localization with promoters and transcriptionally active enhancers suggested new models for in cis and in trans regulation by RNA:DNA hybrids of transcription initiation and formation of 3D-chromatin loops. Overall, this study provides a rationale for the discovery and characterization of the non-B DNA regulatory structures involved in the formation of the RNA:DNA interactome as the basis for an emerging quantitative R-loop biology and pathobiology.

TSSPlant: a new tool for prediction of plant Pol II promoters.

Our current knowledge of eukaryotic promoters indicates their complex architecture that is often composed of numerous functional motifs. Most of known promoters include multiple and in some cases mutually exclusive transcription start sites (TSSs). Moreover, TSS selection depends on cell/tissue, development stage and environmental conditions. Such complex promoter structures make their computational identification notoriously difficult. Here, we present TSSPlant, a novel tool that predicts both TATA and TATA-less promoters in sequences of a wide spectrum of plant genomes. The tool was developed by using large promoter collections from ppdb and PlantProm DB. It utilizes eighteen significant compositional and signal features of plant promoter sequences selected in this study, that feed the artificial neural network-based model trained by the backpropagation algorithm. TSSPlant achieves significantly higher accuracy compared to the next best promoter prediction program for both TATA promoters (MCC≃0.84 and F1-score≃0.91 versus MCC≃0.51 and F1-score≃0.71) and TATA-less promoters (MCC≃0.80, F1-score≃0.89 versus MCC≃0.29 and F1-score≃0.50). TSSPlant is available to download as a standalone program at http://www.cbrc.kaust.edu.sa/download/.

Third promoter and differential regulation of mouse and human fatty acid translocase/CD36 genes

Fatty acid translocase (FAT)/CD36 is a multi-functional membrane protein participating in uptake of long-chain fatty acids and oxidized low density lipoprotein, playing an important role in lipid homeostasis. We reported the dual promoter structures of the mouse and human FAT/CD36 genes in relation to the tissue-specific expression and regulation by peroxisome proliferator-activated receptors (PPARs) (Sato et.al. J. Biol. Chem. 277, 15703-15711, 2002). However, these structures alone could not explain a large induction of the mRNA by a PPAR ligand in the mouse intestine. To learn more about the tissue-specific transcriptional regulation of the FAT/CD36 gene for lipid homeostasis at the whole body level, we cloned the cDNA for the induced mRNA in the mouse intestine and identified third promoter close to the first common exon of the gene. The human gene also has essentially the same gene organization. RT-PCR and Northern blot analyses showed that the three promoters are differentially used in various mouse tissues and the cultured human cells. The physiological significance of these complex promoter structures was suggested by the differential responsiveness to a peroxisome proliferator-activated receptor ligand and the differences in translational efficiencies or stabilities of the mRNA among the isoforms.

T-bet Binding to Newly Identified Target Gene Promoters Is Cell Type-independent but Results in Variable Context-dependent Functional Effects

Recently developed target gene identification strategies based upon the chromatin immunoprecipitation assay provide a powerful method to determine the localization of transcription factor binding within mammalian genomes. However, in many cases, it is unclear if the binding capacity of a transcription factor correlates with an obligate role in gene regulation in diverse contexts. It is therefore important to carefully examine the relationship between transcription factor binding and its ability to functionally regulate gene expression. T-bet is a T-box transcription factor expressed in several hematopoietic cell types. By utilizing a chromatin immunoprecipitation assay coupled to genomic microarray technology approach, we identified numerous promoters, including CXCR3, IL2Rbeta, and CCL3, that are bound by T-bet in B cells. Most surprisingly, the ability of T-bet to associate with the target promoters is not dependent upon the cell type background. Several of the promoters appear to be functionally regulated by T-bet. However, we could not detect a functional consequence for T-bet association with many of the identified promoters in overexpression studies or an examination of wild type and T-bet-/- primary B, CD4+, and CD8+ T cells. Thus, there is a high variability in the functional consequences, if any, that result from the association of T-bet with individual target promoters.

Global RNA half-life analysis in Escherichia coli reveals positional patterns of transcript degradation.

Subgenic-resolution oligonucleotide microarrays were used to study global RNA degradation in wild-type Escherichia coli MG1655. RNA chemical half-lives were measured for 1036 open reading frames (ORFs) and for 329 known and predicted operons. The half-life of total mRNA was 6.8 min under the conditions tested. We also observed significant relationships between gene functional assignments and transcript stability. Unexpectedly, transcription of a single operon (tdcABCDEFG) was relatively rifampicin-insensitive and showed significant increases 2.5 min after rifampicin addition. This supports a novel mechanism of transcription for the tdc operon, whose promoter lacks any recognizable sigma binding sites. Probe by probe analysis of all known and predicted operons showed that the 5' ends of operons degrade, on average, more quickly than the rest of the transcript, with stability increasing in a 3' direction, supporting and further generalizing the current model of a net 5' to 3' directionality of degradation. Hierarchical clustering analysis of operon degradation patterns revealed that this pattern predominates but is not exclusive. We found a weak but highly significant correlation between the degradation of adjacent operon regions, suggesting that stability is determined by a combination of local and operon-wide stability determinants. The 16 ORF dcw gene cluster, which has a complex promoter structure and a partially characterized degradation pattern, was studied at high resolution, allowing a detailed and integrated description of its abundance and degradation. We discuss the application of subgenic resolution DNA microarray analysis to study global mechanisms of RNA transcription and processing.

Both phyA and phyB mediate light-imposed repression of PHYA gene expression in Arabidopsis.

The negatively photoregulated PHYA gene has a complex promoter structure in Arabidopsis, with three active transcription start sites. To identify the photoreceptors responsible for regulation of this gene, and to assess the relative roles of the three transcription start sites, we analyzed the changes in PHYA transcript levels in wild-type and photoreceptor mutant seedlings under various irradiation conditions. Continuous far-red or red light exposures each induced a significant decline in transcript levels in wild-type etiolated seedlings. Analysis of mutants specifically lacking either phyA or phyB protein demonstrated that these phytochromes are required for the negative regulation induced by far-red and red light, respectively. Ribonuclease protection experiments showed further that this negative regulation is confined almost exclusively to the shortest, most abundant PHYA transcript, and occurs predominantly in shoots. By contrast, both of the other minor transcripts in shoots, and all three transcripts in roots, exhibit near constitutive expression. This complex expression pattern indicates that the PHYA gene is subject to regulation by multiple signals, including environmental, developmental, and organ-specific signals.

Complex transcriptional control of the streptokinase gene of Streptococcus equisimilis H46A.

On the Streptococcus equisimilis H46A chromosome, the divergent coding sequences of the genes for the plasminogen activator streptokinase (skc) and a leucine-rich protein (lrp), the function of which is unknown, are separated by a 328 bp intrinsically bent DNA region rich in AT tracts. To begin to understand the expression control of these two genes, we mapped their transcriptional initiation sites by S1 nuclease analysis and studied the influence of the bent intergenic region on promoter strength, using promoter-reporter gene fusions of skc' and lrp' to 'lacZ from Escherichia coli. The major transcriptional start sites, in both S. equisimilis and E. coli, mapped 22 bases upstream of the ATG start site of lrp (G), and 24 and 32 bases upstream of the translational initiation codon of skc (A and G, respectively), indicating the existence of two overlapping canonical skc promoters arranged in tandem on opposite faces of the helix. The reporter gene fusions were cloned in E. coli on a vector containing a 1.1 kb fragment of the S. equisimilis dexB gene, thus allowing promoter strength to be measured in multiple plasmid-form copies in the heterologous host and in single-copy genomic form following integration into the skc region of the homologous host. In S. equisimilis, skc'-'lacZ was expressed about 200-fold more strongly than the corresponding lrp'-'lacZ fusion. In contrast, in E. coli, the corresponding levels of expression differed by only about 11-fold. Deletion of the 202 bp bent region upstream of the skc and lrp core promoters caused a 13-fold decrease in skc promoter activity in S. equisimilis but did not alter lrp promoter strength in this host. In contrast, when studied in E. coli, this deletion did not alter the strength of the skc-double promoter and even increased by 2.4- to 3-fold the activity of the lrp promoter. This comparative promoter analysis shows that skc has a complex promoter structure, the activity of which in the homologous genomic environment specifically depends on sequences upstream of the two core promoters. Thus, the skc promoter structure resembles that of an array of promoters involved in a transcriptional switch; however, the nature of the potential switch factor(s) remains unknown.

Promoter-proximal pausing of RNA polymerase II defines a general rate-limiting step after transcription initiation.

We have shown previously that the majority of RNA polymerase II complexes initiated at the c-myc gene are paused in the promoter-proximal region, similar to observations in the Drosophila hsp70 gene. Our analyses define the TATA box or initiator sequences in the c-myc gene as necessary components for the establishment of paused RNA polymerase II. Deletion of upstream sequences or even the TATA box does not influence significantly the degree of transcriptional initiation or pausing. Deletion of both the TATA box and sequences at the transcription initiation site, however, abolishes transcriptional pausing of transcription complexes but still allows synthesis of full-length RNA. Further analyses with synthetic promoter constructs reveal that the simple combination of upstream activator with TATA consensus sequences or initiator sequences act synergistically to recruit high levels of RNA polymerase II complexes. Only a minor fraction of these complexes escapes into regions further downstream. Several different trans-activation domains fused to GAL4-DNA-binding domains, including strong activators such as VP16, do not eliminate promoter-proximal pausing of RNA polymerase. Thus, we conclude that pausing of RNA polymerase II is a common phenomenon in eukaryotic transcription and does not require complex promoter structures. Further analyses reveal that enhancers have a modest influence on transcription initiation and on release of transcription complexes out of the pause site but may function primarily to increase the elongation competence of transcription complexes.

Identification of sequences responsible for transcriptional regulation of the strongly expressed methanoi oxidase-encoding gene in Hansenula polymorpha

The methylotrophic yeasts have been the subject of intensive studies, because of their highly regulated methanol metabolism and the biogenesis of peroxisomes. We investigated the 5' regulatory region of the MOX gene from the yeast, Hansenula polymorpha, encoding the peroxisomal methanol oxidase, the key enzyme of methanol metabolism. This tightly regulated yeast promoter of approximately 1.5 kb is unusually large, and also of remarkable strength under inducing conditions, belonging to the strongest yeast promoters yet described. Deletion analyses revealed a complex promoter structure composed of several sequence elements with positive and negative regulatory effects on reporter gene expression and a pronounced cooperation between the elements. Specific binding of several factors was detected in vitro by gel retardation and DNase I footprinting experiments. On the basis of deletion data, two binding sites could be identified as upstream activation sequences ( UAS1 and UAS2) and one binding site as an upstream repressing sequence ( URS1).