A quantitative model of transcription factor-activated gene expression - PubMed (original) (raw)
. 2008 Nov;15(11):1192-8.
doi: 10.1038/nsmb.1500. Epub 2008 Oct 12.
Affiliations
- PMID: 18849996
- PMCID: PMC2696132
- DOI: 10.1038/nsmb.1500
A quantitative model of transcription factor-activated gene expression
Harold D Kim et al. Nat Struct Mol Biol. 2008 Nov.
Abstract
A challenge facing biology is to develop quantitative, predictive models of gene regulation. Eukaryotic promoters contain transcription factor binding sites of differing affinity and accessibility, but we understand little about how these variables combine to generate a fine-tuned, quantitative transcriptional response. Here we used the PHO5 promoter in budding yeast to quantify the relationship between transcription factor input and gene expression output, termed the gene-regulation function (GRF). A model that captures variable interactions between transcription factors, nucleosomes and the promoter faithfully reproduced the observed quantitative changes in the GRF that occur upon altering the affinity of transcription factor binding sites, and implicates nucleosome-modulated accessibility of transcription factor binding sites in increasing the diversity of gene expression profiles. This work establishes a quantitative framework that can be applied to predict GRFs of other eukaryotic genes.
Figures
Figure 1
Measuring the GRF from the tetracycline-regulated gene expression system. A variant of a tetracycline-controlled reverse transactivator is constitutively expressed under the control of the MYO2 promoter (PMYO2). When bound to doxycycline, it can activate the TETO7 promoter (PTETO7, containing seven repeats of the tetO operator sequence) and drive the expression of PHO4-YFP, which in turn activates the expression of CFP from a PHO5 promoter variant (PPHO5*). Pho4SA1-4PA6 used in this study is a constitutively nuclear form of Pho4 that has serine-to-alanine substitutions at phosphorylation sites 1, 2, 3 and 4, and a proline-to-alanine substitution at site 6 (ref. 16). Nhp6a, a chromatin protein in the nucleus, was also tagged with RFP to mark the nucleus and to serve as an internal standard for intensity normalization. The chromatin architectures of the PHO5 promoter variants are shown, with nucleosomes (Nuc) −3, −2 and −1 drawn as ellipses. Nucleosome −1 is the closest to the PHO5 open reading frame and occludes the TATA box, which is depicted as a white square with the letter T. Low-affinity (●) and high-affinity (▲) Pho4 binding sites differ by 1 bp out of a 6-bp core sequence (CACGTT versus CACGTG). Variant 4 shown inside the box is the wild-type PHO5 promoter, with a low-affinity site (UASp1) in the exposed region between nucleosome −3 and −2, and a high-affinity site (UASp2) under nucleosome −2. Variants 1, 2, 4, 8, 9 and 12 are labeled as LX, LL, LH, HX, HL and HH, respectively, where H, L and X are the high-affinity sites, the low-affinity site and no site, respectively.
Figure 2
Raw data of fluorescence intensities from single cells and the fits to GRFs. YFP and CFP intensities were obtained from single cells, normalized to RFP intensity and plotted as dashed lines in arbitrary normalized units (NU). As an example, the GRFs of variant 4 and variant 9 are shown in gold and blue, respectively. Variant 4 (LH) is the wild-type strain, and variant 9 (HL) is obtained by swapping UASp1 and UASp2. The Hill curves fit to the data are shown as solid lines in corresponding colors.
Figure 3
The diversity of GRFs and the relationship between the maximum expression level and the nucleosome occupancy over the TATA box region. (a) The GRF of each strain is represented by its maximum expression level on the x axis and its induction threshold on the y axis (NU, normalized unit). The labeling of the promoter variants is the same as in Figure 1. Error bars are one s.d. from a minimum of five different measurements. The variants are grouped by their similarity to one of the four representative variants—LX, LH, HX and HH, colored purple, red, blue and brown, respectively. The distinct features in the promoter architecture of each group are also highlighted in the corresponding group color. The segregation of variants by the horizontal and vertical dashed lines suggests that the induction threshold is determined by the affinity of the exposed binding site (LX and LH versus HX and HH), whereas the maximum expression level depends strongly on the affinity of the binding site under nucleosome −2 (LX and HX versus LH and HH). (b) The normalized histone H3 occupancy was measured by ChIP and plotted against the maximum expression level for each strain. The error bars are one s.d. of six different points, coming from three independent measurements normalized in two different ways (Methods). The inverse relationship supports the idea that the nucleosome occupancy over the TATA box region limits gene expression.
Figure 4
Quantitative models of Pho4-dependent chromatin remodeling. (a) This minimal model describes a promoter composed of one exposed Pho4 binding site and one adjacent nucleosome containing the TATA box. The presence or absence of Pho4 and the nucleosome (Nuc) are marked by O or X in the table, the combination of which defines the four states in this model. The orange pentagon represents Pho4, and all other symbols are as defined in Figure 1. The two states on the right are considered to be transcriptionally active as the general transcription machinery can access the TATA box region. The transitional frequencies of association, dissociation, remodeling and reassembly are described by kassoc∗, _k_dissoc, _k_remod and _k_reass, respectively. (b) This model expands on the minimal model in a to describe the variants LX, LL, LH, HX, HL and HH. The six states on the right comprise the transcriptionally active fraction. The rate constants in this model are for the association of Pho4 to DNA (kassoc∗), dissociation from the exposed (kdissocexp) or nucleosomal region (kdissocnuc), remodeling of nucleosome −2 or −1 (_k_remod) and the reassembly of nucleosome −2 or −1 (_k_reass). The arrow scheme for these rate constants is shown on the left. The entire model with three layers of squares is the model with nucleosome −2, and the model with the inner two layers highlighted in green corresponds to the model without nucleosome −2.
Figure 5
Comparison between the data and the model prediction. (a) We globally fitted six threshold values and six maximum expression levels from variants LX, LL, LH, HX, HL and HH using two different models—those with and without nucleosome −2. The experimental data (blue squares with error bars) are compared to the results from the model with nucleosome −2 (purple circles) and the model without nucleosome −2 (red circles). The promoter variants used are listed on the left axis. (b) Gene expression landscapes from the models with and without nucleosome −2. The maximum expression level and the induction threshold are plotted as a function of k¯dissocnuc (x axis) and k¯dissocexp (y axis) as heat maps. The other parameters were set to the values obtained from the global fit (Supplementary Discussion). The color scale is shown next to each heat map (AU, arbitrary unit). k¯dissocnuc and k¯dissocexp values of the six variants determined by the global fitting are plotted on each heat map: LX (●), LL (▼), LH (◄), HX (■), HL (▲) and HH (►). For these variants, k¯dissocexp can be either _k̄_H or _k̄_L, and k¯dissocnuc can be _k̄_H, _k̄_Lor infinity (∞). _k̄_H and _k̄_L are 0.17 and 1.27, respectively, with nucleosome −2, and 0.34 and 2.76, respectively, without nucleosome −2. For convenience, k¯dissocnuc = infinity is plotted as the maximum k¯dissocnuc value shown in the heat map. Nucleosome −2, by occluding one of the transcription factor binding sites, breaks the anti-correlated change of the maximum expression level and the induction threshold, thus increasing the variety of gene expression profiles.
Similar articles
- Nucleosomes Are Essential for Proper Regulation of a Multigated Promoter in Saccharomyces cerevisiae.
Yarrington RM, Goodrum JM, Stillman DJ. Yarrington RM, et al. Genetics. 2016 Feb;202(2):551-63. doi: 10.1534/genetics.115.183715. Epub 2015 Dec 1. Genetics. 2016. PMID: 26627840 Free PMC article. - Genome-wide analysis of transcriptional dependence and probable target sites for Abf1 and Rap1 in Saccharomyces cerevisiae.
Yarragudi A, Parfrey LW, Morse RH. Yarragudi A, et al. Nucleic Acids Res. 2007;35(1):193-202. doi: 10.1093/nar/gkl1059. Epub 2006 Dec 7. Nucleic Acids Res. 2007. PMID: 17158163 Free PMC article. - Nucleosome repositioning underlies dynamic gene expression.
Nocetti N, Whitehouse I. Nocetti N, et al. Genes Dev. 2016 Mar 15;30(6):660-72. doi: 10.1101/gad.274910.115. Epub 2016 Mar 10. Genes Dev. 2016. PMID: 26966245 Free PMC article. - Transcription factors vs nucleosomes: regulation of the PHO5 promoter in yeast.
Svaren J, Hörz W. Svaren J, et al. Trends Biochem Sci. 1997 Mar;22(3):93-7. doi: 10.1016/s0968-0004(97)01001-3. Trends Biochem Sci. 1997. PMID: 9066259 Review. - The yeast PHO5 promoter: from single locus to systems biology of a paradigm for gene regulation through chromatin.
Korber P, Barbaric S. Korber P, et al. Nucleic Acids Res. 2014;42(17):10888-902. doi: 10.1093/nar/gku784. Epub 2014 Sep 4. Nucleic Acids Res. 2014. PMID: 25190457 Free PMC article. Review.
Cited by
- MONITTR allows real-time imaging of transcription and endogenous proteins in C. elegans.
Liu X, Chang Z, Sun P, Cao B, Wang Y, Fang J, Pei Y, Chen B, Zou W. Liu X, et al. J Cell Biol. 2025 Jan 6;224(1):e202403198. doi: 10.1083/jcb.202403198. Epub 2024 Oct 14. J Cell Biol. 2025. PMID: 39400293 - Live-cell analysis of IMPDH protein levels during yeast colony growth provides insights into the regulation of GTP synthesis.
Shand EL, Sweeney K, Sundling KE, McClean MN, Brow DA. Shand EL, et al. mBio. 2024 Aug 14;15(8):e0102124. doi: 10.1128/mbio.01021-24. Epub 2024 Jun 28. mBio. 2024. PMID: 38940616 Free PMC article. - Single-molecule chromatin configurations link transcription factor binding to expression in human cells.
Doughty BR, Hinks MM, Schaepe JM, Marinov GK, Thurm AR, Rios-Martinez C, Parks BE, Tan Y, Marklund E, Dubocanin D, Bintu L, Greenleaf WJ. Doughty BR, et al. bioRxiv [Preprint]. 2024 Feb 4:2024.02.02.578660. doi: 10.1101/2024.02.02.578660. bioRxiv. 2024. PMID: 38352517 Free PMC article. Preprint. - PGE2 alters chromatin through H2A.Z-variant enhancer nucleosome modification to promote hematopoietic stem cell fate.
Sporrij A, Choudhuri A, Prasad M, Muhire B, Fast EM, Manning ME, Weiss JD, Koh M, Yang S, Kingston RE, Tolstorukov MY, Clevers H, Zon LI. Sporrij A, et al. Proc Natl Acad Sci U S A. 2023 May 9;120(19):e2220613120. doi: 10.1073/pnas.2220613120. Epub 2023 May 1. Proc Natl Acad Sci U S A. 2023. PMID: 37126722 Free PMC article. - Minimal synthetic enhancers reveal control of the probability of transcriptional engagement and its timing by a morphogen gradient.
Alamos S, Reimer A, Westrum C, Turner MA, Talledo P, Zhao J, Luu E, Garcia HG. Alamos S, et al. Cell Syst. 2023 Mar 15;14(3):220-236.e3. doi: 10.1016/j.cels.2022.12.008. Epub 2023 Jan 24. Cell Syst. 2023. PMID: 36696901 Free PMC article.
References
- Cyert MS. Regulation of nuclear localization during signaling. J Biol Chem. 2001;276:20805–20808. - PubMed
Publication types
MeSH terms
Substances
Grants and funding
- R01 GM051377-15/GM/NIGMS NIH HHS/United States
- Howard Hughes Medical Institute/United States
- GM51377/GM/NIGMS NIH HHS/United States
- R01 GM051377-16/GM/NIGMS NIH HHS/United States
- R01 GM051377/GM/NIGMS NIH HHS/United States
LinkOut - more resources
Full Text Sources
Molecular Biology Databases