Histone modifications within the human X centromere region.

Human centromeres are multi-megabase regions of highly ordered arrays of alpha satellite DNA that are separated from chromosome arms by unordered alpha satellite monomers and other repetitive elements. Complexities in assembling such large repetitive regions have limited detailed studies of centromeric chromatin organization. However, a genomic map of the human X centromere has provided new opportunities to explore genomic architecture of a complex locus. We used ChIP to examine the distribution of modified histones within centromere regions of multiple X chromosomes. Methylation of H3 at lysine 4 coincided with DXZ1 higher order alpha satellite, the site of CENP-A localization. Heterochromatic histone modifications were distributed across the 400-500 kb pericentromeric regions. The large arrays of alpha satellite and gamma satellite DNA were enriched for both euchromatic and heterochromatic modifications, implying that some pericentromeric repeats have multiple chromatin characteristics. Partial truncation of the X centromere resulted in reduction in the size of the CENP-A/Cenp-A domain and increased heterochromatic modifications in the flanking pericentromere. Although the deletion removed approximately 1/3 of centromeric DNA, the ratio of CENP-A to alpha satellite array size was maintained in the same proportion, suggesting that a limited, but defined linear region of the centromeric DNA is necessary for kinetochore assembly. Our results indicate that the human X centromere contains multiple types of chromatin, is organized similarly to smaller eukaryotic centromeres, and responds to structural changes by expanding or contracting domains.

Defining the Role of the Histone Methyltransferase, PR-Set7, in Maintaining the Genome Integrity of Drosophila Melanogaster

The complete and faithful duplication of the genome is essential to ensure normal cell division and organismal development. Eukaryotic DNA replication is initiated at multiple sites termed origins of replication that are activated at different time through S phase. The replication timing program is regulated by the S-phase checkpoint, which signals and repairs replicative stress. Eukaryotic DNA is packaged with histones into chromatin, thus DNA-templated processes including replication are modulated by the local chromatin environment such as post-translational modifications (PTMs) of histones.

One such epigenetic mark, methylation of lysine 20 on histone H4 (H4K20), has been linked to chromatin compaction, transcription, DNA repair and DNA replication. H4K20 can be mono-, di- and tri-methylated. Monomethylation of H4K20 (H4K20me1) is mediated by the cell cycle-regulated histone methyltransferase PR-Set7 and subsequent di-/tri- methylation is catalyzed by Suv4-20. Prior studies have shown that PR-Set7 depletion in mammalian cells results in defective S phase progression and the accumulation of DNA damage, which may be partially attributed to defects in origin selection and activation. Meanwhile, overexpression of mammalian PR-Set7 recruits components of pre-Replication Complex (pre-RC) onto chromatin and licenses replication origins for re-replication. However, these studies were limited to only a handful of mammalian origins, and it remains unclear how PR-Set7 impacts the replication program on a genomic scale. Finally, the methylation substrates of PR-Set7 include both histone (H4K20) and non-histone targets, therefore it is necessary to directly test the role of H4K20 methylation in PR-Set7 regulated phenotypes.

I employed genetic, cytological, and genomic approaches to better understand the role of H4K20 methylation in regulating DNA replication and genome stability in Drosophila melanogaster cells. Depletion of Drosophila PR-Set7 by RNAi in cultured Kc167 cells led to an ATR-dependent cell cycle arrest with near 4N DNA content and the accumulation of DNA damage, indicating a defect in completing S phase. The cells were arrested at the second S phase following PR-Set7 downregulation, suggesting that it was an epigenetic effect that coupled to the dilution of histone modification over multiple cell cycles. To directly test the role of H4K20 methylation in regulating genome integrity, I collaborated with the Duronio Lab and observed spontaneous DNA damage on the imaginal wing discs of third instar mutant larvae that had an alanine substitution on H4K20 (H4K20A) thus unable to be methylated, confirming that H4K20 is a bona fide target of PR-Set7 in maintaining genome integrity.

One possible source of DNA damage due to loss of PR-Set7 is reduced origin activity. I used BrdU-seq to profile the genome-wide origin activation pattern. However, I found that deregulation of H4K20 methylation states by manipulating the H4K20 methyltransferases PR-Set7 and Suv4-20 had no impact on origin activation throughout the genome. I then mapped the genomic distribution of DNA damage upon PR-Set7 depletion. Surprisingly, ChIP-seq of the DNA damage marker γ-H2A.v located the DNA damage to late replicating euchromatic regions of the Drosophila genome, and the strength of γ-H2A.v signal was uniformly distributed and spanned the entire late replication domain, implying stochastic replication fork collapse within late replicating regions. Together these data suggest that PR-Set7-mediated monomethylation of H4K20 is critical for maintaining the genomic integrity of late replicating domains, presumably via stabilization of late replicating forks.

In addition to investigating the function of H4K20me, I also used immunofluorescence to characterize the cell cycle regulated chromatin loading of Mcm2-7 complex, the DNA helicase that licenses replication origins, using H4K20me1 level as a proxy for cell cycle stages. In parallel with chromatin spindown data by Powell et al. (Powell et al. 2015), we showed a continuous loading of Mcm2-7 during G1 and a progressive removal from chromatin through S phase.

Quality control autophagy: a joint effort of ubiquitin, protein deacetylase and actin cytoskeleton.

Autophagy has been predominantly studied as a nonselective self-digestion process that recycles macromolecules and produces energy in response to starvation. However, autophagy independent of nutrient status has long been known to exist. Recent evidence suggests that this form of autophagy enforces intracellular quality control by selectively disposing of aberrant protein aggregates and damaged organelles--common denominators in various forms of neurodegenerative diseases. By definition, this form of autophagy, termed quality-control (QC) autophagy, must be different from nutrient-regulated autophagy in substrate selectivity, regulation and function. We have recently identified the ubiquitin-binding deacetylase, HDAC6, as a key component that establishes QC. HDAC6 is not required for autophagy activation per se; rather, it is recruited to ubiquitinated autophagic substrates where it stimulates autophagosome-lysosome fusion by promoting F-actin remodeling in a cortactin-dependent manner. Remarkably, HDAC6 and cortactin are dispensable for starvation-induced autophagy. These findings reveal that autophagosomes associated with QC are molecularly and biochemically distinct from those associated with starvation autophagy, thereby providing a new molecular framework to understand the emerging complexity of autophagy and therapeutic potential of this unique machinery.

HDAC6 and Ubp-M BUZ domains recognize specific C-terminal sequences of proteins.

The BUZ/Znf-UBP domain is a protein module found in the cytoplasmic deacetylase HDAC6, E3 ubiquitin ligase BRAP2/IMP, and a subfamily of ubiquitin-specific proteases. Although several BUZ domains have been shown to bind ubiquitin with high affinity by recognizing its C-terminal sequence (RLRGG-COOH), it is currently unknown whether the interaction is sequence-specific or whether the BUZ domains are capable of binding to proteins other than ubiquitin. In this work, the BUZ domains of HDAC6 and Ubp-M were subjected to screening against a one-bead-one-compound (OBOC) peptide library that exhibited random peptide sequences with free C-termini. Sequence analysis of the selected binding peptides as well as alanine scanning studies revealed that the BUZ domains require a C-terminal Gly-Gly motif for binding. At the more N-terminal positions, the two BUZ domains have distinct sequence specificities, allowing them to bind to different peptides and/or proteins. A database search of the human proteome on the basis of the BUZ domain specificities identified 11 and 24 potential partner proteins for Ubp-M and HDAC6 BUZ domains, respectively. Peptides corresponding to the C-terminal sequences of four of the predicted binding partners (FBXO11, histone H4, PTOV1, and FAT10) were synthesized and tested for binding to the BUZ domains by fluorescence polarization. All four peptides bound to the HDAC6 BUZ domain with low micromolar K(D) values and less tightly to the Ubp-M BUZ domain. Finally, in vitro pull-down assays showed that the Ubp-M BUZ domain was capable of binding to the histone H3-histone H4 tetramer protein complex. Our results suggest that BUZ domains are sequence-specific protein-binding modules, with each BUZ domain potentially binding to a different subset of proteins.

Thermodynamic analysis of a molecular chaperone binding to unfolded protein substrates.

Molecular chaperones are a highly diverse group of proteins that recognize and bind unfolded proteins to facilitate protein folding and prevent nonspecific protein aggregation. The mechanisms by which chaperones bind their protein substrates have been studied for decades. However, there are few reports about the affinity of molecular chaperones for their unfolded protein substrates. Thus, little is known about the relative binding affinities of different chaperones and about the relative binding affinities of chaperones for different unfolded protein substrates. Here we describe the application of SUPREX (stability of unpurified proteins from rates of H-D exchange), an H-D exchange and MALDI-based technique, in studying the binding interaction between the molecular chaperone Hsp33 and four different unfolded protein substrates, including citrate synthase, lactate dehydrogenase, malate dehydrogenase, and aldolase. The results of our studies suggest that the cooperativity of the Hsp33 folding-unfolding reaction increases upon binding with denatured protein substrates. This is consistent with the burial of significant hydrophobic surface area in Hsp33 when it interacts with its substrate proteins. The SUPREX-derived K(d) values for Hsp33 complexes with four different substrates were all found to be within the range of 3-300 nM.

Neuroprotection resulting from insufficiency of RANBP2 is associated with the modulation of protein and lipid homeostasis of functionally diverse but linked pathways in response to oxidative stress.

Oxidative stress is a deleterious stressor associated with a plethora of disease and aging manifestations, including neurodegenerative disorders, yet very few factors and mechanisms promoting the neuroprotection of photoreceptor and other neurons against oxidative stress are known. Insufficiency of RAN-binding protein-2 (RANBP2), a large, mosaic protein with pleiotropic functions, suppresses apoptosis of photoreceptor neurons upon aging and light-elicited oxidative stress, and promotes age-dependent tumorigenesis by mechanisms that are not well understood. Here we show that, by downregulating selective partners of RANBP2, such as RAN GTPase, UBC9 and ErbB-2 (HER2; Neu), and blunting the upregulation of a set of orphan nuclear receptors and the light-dependent accumulation of ubiquitylated substrates, light-elicited oxidative stress and Ranbp2 haploinsufficiency have a selective effect on protein homeostasis in the retina. Among the nuclear orphan receptors affected by insufficiency of RANBP2, we identified an isoform of COUP-TFI (Nr2f1) as the only receptor stably co-associating in vivo with RANBP2 and distinct isoforms of UBC9. Strikingly, most changes in proteostasis caused by insufficiency of RANBP2 in the retina are not observed in the supporting tissue, the retinal pigment epithelium (RPE). Instead, insufficiency of RANBP2 in the RPE prominently suppresses the light-dependent accumulation of lipophilic deposits, and it has divergent effects on the accumulation of free cholesterol and free fatty acids despite the genotype-independent increase of light-elicited oxidative stress in this tissue. Thus, the data indicate that insufficiency of RANBP2 results in the cell-type-dependent downregulation of protein and lipid homeostasis, acting on functionally interconnected pathways in response to oxidative stress. These results provide a rationale for the neuroprotection from light damage of photosensory neurons by RANBP2 insufficiency and for the identification of novel therapeutic targets and approaches promoting neuroprotection.

Modulation of heat shock transcription factor 1 as a therapeutic target for small molecule intervention in neurodegenerative disease.

Neurodegenerative diseases such as Huntington disease are devastating disorders with no therapeutic approaches to ameliorate the underlying protein misfolding defect inherent to poly-glutamine (polyQ) proteins. Given the mounting evidence that elevated levels of protein chaperones suppress polyQ protein misfolding, the master regulator of protein chaperone gene transcription, HSF1, is an attractive target for small molecule intervention. We describe a humanized yeast-based high-throughput screen to identify small molecule activators of human HSF1. This screen is insensitive to previously characterized activators of the heat shock response that have undesirable proteotoxic activity or that inhibit Hsp90, the central chaperone for cellular signaling and proliferation. A molecule identified in this screen, HSF1A, is structurally distinct from other characterized small molecule human HSF1 activators, activates HSF1 in mammalian and fly cells, elevates protein chaperone expression, ameliorates protein misfolding and cell death in polyQ-expressing neuronal precursor cells and protects against cytotoxicity in a fly model of polyQ-mediated neurodegeneration. In addition, we show that HSF1A interacts with components of the TRiC/CCT complex, suggesting a potentially novel regulatory role for this complex in modulating HSF1 activity. These studies describe a novel approach for the identification of new classes of pharmacological interventions for protein misfolding that underlies devastating neurodegenerative disease.

Yeast screens identify the RNA polymerase II CTD and SPT5 as relevant targets of BRCA1 interaction.

BRCA1 has been implicated in numerous DNA repair pathways that maintain genome integrity, however the function responsible for its tumor suppressor activity in breast cancer remains obscure. To identify the most highly conserved of the many BRCA1 functions, we screened the evolutionarily distant eukaryote Saccharomyces cerevisiae for mutants that suppressed the G1 checkpoint arrest and lethality induced following heterologous BRCA1 expression. A genome-wide screen in the diploid deletion collection combined with a screen of ionizing radiation sensitive gene deletions identified mutants that permit growth in the presence of BRCA1. These genes delineate a metabolic mRNA pathway that temporally links transcription elongation (SPT4, SPT5, CTK1, DEF1) to nucleopore-mediated mRNA export (ASM4, MLP1, MLP2, NUP2, NUP53, NUP120, NUP133, NUP170, NUP188, POM34) and cytoplasmic mRNA decay at P-bodies (CCR4, DHH1). Strikingly, BRCA1 interacted with the phosphorylated RNA polymerase II (RNAPII) carboxy terminal domain (P-CTD), phosphorylated in the pattern specified by the CTDK-I kinase, to induce DEF1-dependent cleavage and accumulation of a RNAPII fragment containing the P-CTD. Significantly, breast cancer associated BRCT domain defects in BRCA1 that suppressed P-CTD cleavage and lethality in yeast also suppressed the physical interaction of BRCA1 with human SPT5 in breast epithelial cells, thus confirming SPT5 as a relevant target of BRCA1 interaction. Furthermore, enhanced P-CTD cleavage was observed in both yeast and human breast cells following UV-irradiation indicating a conserved eukaryotic damage response. Moreover, P-CTD cleavage in breast epithelial cells was BRCA1-dependent since damage-induced P-CTD cleavage was only observed in the mutant BRCA1 cell line HCC1937 following ectopic expression of wild type BRCA1. Finally, BRCA1, SPT5 and hyperphosphorylated RPB1 form a complex that was rapidly degraded following MMS treatment in wild type but not BRCA1 mutant breast cells. These results extend the mechanistic links between BRCA1 and transcriptional consequences in response to DNA damage and suggest an important role for RNAPII P-CTD cleavage in BRCA1-mediated cancer suppression.

Risk of ovarian cancer and inherited variants in relapse-associated genes.

BACKGROUND: We previously identified a panel of genes associated with outcome of ovarian cancer. The purpose of the current study was to assess whether variants in these genes correlated with ovarian cancer risk. METHODS AND FINDINGS: Women with and without invasive ovarian cancer (749 cases, 1,041 controls) were genotyped at 136 single nucleotide polymorphisms (SNPs) within 13 candidate genes. Risk was estimated for each SNP and for overall variation within each gene. At the gene-level, variation within MSL1 (male-specific lethal-1 homolog) was associated with risk of serous cancer (p = 0.03); haplotypes within PRPF31 (PRP31 pre-mRNA processing factor 31 homolog) were associated with risk of invasive disease (p = 0.03). MSL1 rs7211770 was associated with decreased risk of serous disease (OR 0.81, 95% CI 0.66-0.98; p = 0.03). SNPs in MFSD7, BTN3A3, ZNF200, PTPRS, and CCND1A were inversely associated with risk (p<0.05), and there was increased risk at HEXIM1 rs1053578 (p = 0.04, OR 1.40, 95% CI 1.02-1.91). CONCLUSIONS: Tumor studies can reveal novel genes worthy of follow-up for cancer susceptibility. Here, we found that inherited markers in the gene encoding MSL1, part of a complex that modifies the histone H4, may decrease risk of invasive serous ovarian cancer.

A novel, non-apoptotic role for Scythe/BAT3: a functional switch between the pro- and anti-proliferative roles of p21 during the cell cycle.

BACKGROUND: Scythe/BAT3 is a member of the BAG protein family whose role in apoptosis has been extensively studied. However, since the developmental defects observed in Bat3-null mouse embryos cannot be explained solely by defects in apoptosis, we investigated whether BAT3 is also involved in cell-cycle progression. METHODS/PRINCIPAL FINDINGS: Using a stable-inducible Bat3-knockdown cellular system, we demonstrated that reduced BAT3 protein level causes a delay in both G1/S transition and G2/M progression. Concurrent with these changes in cell-cycle progression, we observed a reduction in the turnover and phosphorylation of the CDK inhibitor p21, which is best known as an inhibitor of DNA replication; however, phosphorylated p21 has also been shown to promote G2/M progression. Our findings indicate that in Bat3-knockdown cells, p21 continues to be synthesized during cell-cycle phases that do not normally require p21, resulting in p21 protein accumulation and a subsequent delay in cell-cycle progression. Finally, we showed that BAT3 co-localizes with p21 during the cell cycle and is required for the translocation of p21 from the cytoplasm to the nucleus during the G1/S transition and G2/M progression. CONCLUSION: Our study reveals a novel, non-apoptotic role for BAT3 in cell-cycle regulation. By maintaining a low p21 protein level during the G1/S transition, BAT3 counteracts the inhibitory effect of p21 on DNA replication and thus enables the cells to progress from G1 to S phase. Conversely, during G2/M progression, BAT3 facilitates p21 phosphorylation by cyclin A/Cdk2, an event required for G2/M progression. BAT3 modulates these pro- and anti-proliferative roles of p21 at least in part by regulating cyclin A abundance, as well as p21 translocation between the cytoplasm and the nucleus to ensure that it functions in the appropriate intracellular compartment during each phase of the cell cycle.

Transcription factors MYOCD, SRF, Mesp1 and SMARCD3 enhance the cardio-inducing effect of GATA4, TBX5, and MEF2C during direct cellular reprogramming.

Transient overexpression of defined combinations of master regulator genes can effectively induce cellular reprogramming: the acquisition of an alternative predicted phenotype from a differentiated cell lineage. This can be of particular importance in cardiac regenerative medicine wherein the heart lacks the capacity to heal itself, but simultaneously contains a large pool of fibroblasts. In this study we determined the cardio-inducing capacity of ten transcription factors to actuate cellular reprogramming of mouse embryonic fibroblasts into cardiomyocyte-like cells. Overexpression of transcription factors MYOCD and SRF alone or in conjunction with Mesp1 and SMARCD3 enhanced the basal but necessary cardio-inducing effect of the previously reported GATA4, TBX5, and MEF2C. In particular, combinations of five or seven transcription factors enhanced the activation of cardiac reporter vectors, and induced an upregulation of cardiac-specific genes. Global gene expression analysis also demonstrated a significantly greater cardio-inducing effect when the transcription factors MYOCD and SRF were used. Detection of cross-striated cells was highly dependent on the cell culture conditions and was enhanced by the addition of valproic acid and JAK inhibitor. Although we detected Ca(2+) transient oscillations in the reprogrammed cells, we did not detect significant changes in resting membrane potential or spontaneously contracting cells. This study further elucidates the cardio-inducing effect of the transcriptional networks involved in cardiac cellular reprogramming, contributing to the ongoing rational design of a robust protocol required for cardiac regenerative therapies.

Copper signaling axis as a target for prostate cancer therapeutics.

Previously published reports indicate that serum copper levels are elevated in patients with prostate cancer and that increased copper uptake can be used as a means to image prostate tumors. It is unclear, however, to what extent copper is required for prostate cancer cell function as we observed only modest effects of chelation strategies on the growth of these cells in vitro. With the goal of exploiting prostate cancer cell proclivity for copper uptake, we developed a "conditional lethal" screen to identify compounds whose cytotoxic actions were manifested in a copper-dependent manner. Emerging from this screen was a series of dithiocarbamates, which, when complexed with copper, induced reactive oxygen species-dependent apoptosis of malignant, but not normal, prostate cells. One of the dithiocarbamates identified, disulfiram (DSF), is an FDA-approved drug that has previously yielded disappointing results in clinical trials in patients with recurrent prostate cancer. Similarly, in our studies, DSF alone had a minimal effect on the growth of prostate cancer tumors when propagated as xenografts. However, when DSF was coadministered with copper, a very dramatic inhibition of tumor growth in models of hormone-sensitive and of castrate-resistant disease was observed. Furthermore, we determined that prostate cancer cells express high levels of CTR1, the primary copper transporter, and additional chaperones that are required to maintain intracellular copper homeostasis. The expression levels of most of these proteins are increased further upon treatment of androgen receptor (AR)-positive prostate cancer cell lines with androgens. Not surprisingly, robust CTR1-dependent uptake of copper into prostate cancer cells was observed, an activity that was accentuated by activation of AR. Given these data linking AR to intracellular copper uptake, we believe that dithiocarbamate/copper complexes are likely to be effective for the treatment of patients with prostate cancer whose disease is resistant to classical androgen ablation therapies.

Salicylic acid biosynthesis is enhanced and contributes to increased biotrophic pathogen resistance in Arabidopsis hybrids.

Heterosis, the phenotypic superiority of a hybrid over its parents, has been demonstrated for many traits in Arabidopsis thaliana, but its effect on defence remains largely unexplored. Here, we show that hybrids between some A. thaliana accessions show increased resistance to the biotrophic bacterial pathogen Pseudomonas syringae pv. tomato (Pst) DC3000. Comparisons of transcriptomes between these hybrids and their parents after inoculation reveal that several key salicylic acid (SA) biosynthesis genes are significantly upregulated in hybrids. Moreover, SA levels are higher in hybrids than in either parent. Increased resistance to Pst DC3000 is significantly compromised in hybrids of pad4 mutants in which the SA biosynthesis pathway is blocked. Finally, increased histone H3 acetylation of key SA biosynthesis genes correlates with their upregulation in infected hybrids. Our data demonstrate that enhanced activation of SA biosynthesis in A. thaliana hybrids may contribute to their increased resistance to a biotrophic bacterial pathogen.

Predictive Modeling of Loss-of-Heterozygosity Events in Yeast

During mitotic cell cycles, DNA experiences many types of endogenous and exogenous damaging agents that could potentially cause double strand breaks (DSB). In S. cerevisiae, DSBs are primarily repaired by mitotic recombination and as a result, could lead to loss-of-heterozygosity (LOH). Genetic recombination can happen in both meiosis and mitosis. While genome-wide distribution of meiotic recombination events has been intensively studied, mitotic recombination events have not been mapped unbiasedly throughout the genome until recently. Methods for selecting mitotic crossovers and mapping the positions of crossovers have recently been developed in our lab. Our current approach uses a diploid yeast strain that is heterozygous for about 55,000 SNPs, and employs SNP-Microarrays to map LOH events throughout the genome. These methods allow us to examine selected crossovers and unselected mitotic recombination events (crossover, noncrossover and BIR) at about 1 kb resolution across the genome. Using this method, we generated maps of spontaneous and UV-induced LOH events. In this study, we explore machine learning and variable selection techniques to build a predictive model for where the LOH events occur in the genome.

Randomly from the yeast genome, we simulated control tracts resembling the LOH tracts in terms of tract lengths and locations with respect to single-nucleotide-polymorphism positions. We then extracted roughly 1,100 features such as base compositions, histone modifications, presence of tandem repeats etc. and train classifiers to distinguish control tracts and LOH tracts. We found interesting features of good predictive values. We also found that with the current repertoire of features, the prediction is generally better for spontaneous LOH events than UV-induced LOH events.

Different Mechanisms Confer Gradual Control and Memory at Nutrient- and Stress-Regulated Genes in Yeast.

Cells respond to environmental stimuli by fine-tuned regulation of gene expression. Here we investigated the dose-dependent modulation of gene expression at high temporal resolution in response to nutrient and stress signals in yeast. The GAL1 activity in cell populations is modulated in a well-defined range of galactose concentrations, correlating with a dynamic change of histone remodeling and RNA polymerase II (RNAPII) association. This behavior is the result of a heterogeneous induction delay caused by decreasing inducer concentrations across the population. Chromatin remodeling appears to be the basis for the dynamic GAL1 expression, because mutants with impaired histone dynamics show severely truncated dose-response profiles. In contrast, the GRE2 promoter operates like a rapid off/on switch in response to increasing osmotic stress, with almost constant expression rates and exclusively temporal regulation of histone remodeling and RNAPII occupancy. The Gal3 inducer and the Hog1 mitogen-activated protein (MAP) kinase seem to determine the different dose-response strategies at the two promoters. Accordingly, GAL1 becomes highly sensitive and dose independent if previously stimulated because of residual Gal3 levels, whereas GRE2 expression diminishes upon repeated stimulation due to acquired stress resistance. Our analysis reveals important differences in the way dynamic signals create dose-sensitive gene expression outputs.