FUBP1 is a general splicing factor facilitating 3' splice site recognition and splicing of long introns.

Splicing of pre-mRNAs critically contributes to gene regulation and proteome expansion in eukaryotes, but our understanding of the recognition and pairing of splice sites during spliceosome assembly lacks detail. Here, we identify the multidomain RNA-binding protein FUBP1 as a key splicing factor that binds to a hitherto unknown cis-regulatory motif. By collecting NMR, structural, and in vivo interaction data, we demonstrate that FUBP1 stabilizes U2AF2 and SF1, key components at the 3' splice site, through multivalent binding interfaces located within its disordered regions. Transcriptional profiling and kinetic modeling reveal that FUBP1 is required for efficient splicing of long introns, which is impaired in cancer patients harboring FUBP1 mutations. Notably, FUBP1 interacts with numerous U1 snRNP-associated proteins, suggesting a unique role for FUBP1 in splice site bridging for long introns. We propose a compelling model for 3' splice site recognition of long introns, which represent 80% of all human introns.

State- and stimulus-specific dynamics of SMAD signaling determine fate decisions in individual cells.

SMAD-mediated signaling regulates apoptosis, cell cycle arrest, and epithelial-to-mesenchymal transition to safeguard tissue homeostasis. However, it remains elusive how the relatively simple pathway can determine such a broad range of cell fate decisions and how it differentiates between varying ligands. Here, we systematically investigate how SMAD-mediated responses are modulated by various ligands of the transforming growth factor ß (TGFß) family and compare these ligand responses in quiescent and proliferating MCF10A cells. We find that the nature of the phenotypic response is mainly determined by the proliferation status, with migration and cell cycle arrest being dominant in proliferating cells for all tested TGFß family ligands, whereas cell death is the major outcome in quiescent cells. In both quiescent and proliferating cells, the identity of the ligand modulates the strength of the phenotypic response proportional to the dynamics of induced SMAD nuclear-to-cytoplasmic translocation and, as a consequence, the corresponding gene expression changes. Interestingly, the proliferation state of a cell has little impact on the set of genes induced by SMAD signaling; instead, it modulates the relative cellular sensitivity to TGFß superfamily members. Taken together, diversity of SMAD-mediated responses is mediated by differing cellular states, which determine ligand sensitivity and phenotypic effects, while the pathway itself merely serves as a quantitative relay from the cell membrane to the nucleus.

Network switches and their role in circadian clocks.

Circadian rhythms are generated by complex interactions among genes and proteins. Self-sustained ∼24 h oscillations require negative feedback loops and sufficiently strong nonlinearities that are the product of molecular and network switches. Here, we review common mechanisms to obtain switch-like behavior, including cooperativity, antagonistic enzymes, multisite phosphorylation, positive feedback, and sequestration. We discuss how network switches play a crucial role as essential components in cellular circadian clocks, serving as integral parts of transcription-translation feedback loops that form the basis of circadian rhythm generation. The design principles of network switches and circadian clocks are illustrated by representative mathematical models that include bistable systems and negative feedback loops combined with Hill functions. This work underscores the importance of negative feedback loops and network switches as essential design principles for biological oscillations, emphasizing how an understanding of theoretical concepts can provide insights into the mechanisms generating biological rhythms.

Data-based stochastic modeling reveals sources of activity bursts in single-cell TGF-ß signaling.

Cells sense their surrounding by employing intracellular signaling pathways that transmit hormonal signals from the cell membrane to the nucleus. TGF-ß/SMAD signaling encodes various cell fates, controls tissue homeostasis and is deregulated in diseases such as cancer. The pathway shows strong heterogeneity at the single-cell level, but quantitative insights into mechanisms underlying fluctuations at various time scales are still missing, partly due to inefficiency in the calibration of stochastic models that mechanistically describe signaling processes. In this work we analyze single-cell TGF-ß/SMAD signaling and show that it exhibits temporal stochastic bursts which are dose-dependent and whose number and magnitude correlate with cell migration. We propose a stochastic modeling approach to mechanistically describe these pathway fluctuations with high computational efficiency. Employing high-order numerical integration and fitting to burst statistics we enable efficient quantitative parameter estimation and discriminate models that assume noise in different reactions at the receptor level. This modeling approach suggests that stochasticity in the internalization of TGF-ß receptors into endosomes plays a key role in the observed temporal bursting. Further, the model predicts the single-cell dynamics of TGF-ß/SMAD signaling in untested conditions, e.g., successfully reflects memory effects of signaling noise and cellular sensitivity towards repeated stimulation. Taken together, our computational framework based on burst analysis, noise modeling and path computation scheme is a suitable tool for the data-based modeling of complex signaling pathways, capable of identifying the source of temporal noise.

Robust Ordering of Anaphase Events by Adaptive Thresholds and Competing Degradation Pathways.

The splitting of chromosomes in anaphase and their delivery into the daughter cells needs to be accurately executed to maintain genome stability. Chromosome splitting requires the degradation of securin, whereas the distribution of the chromosomes into the daughter cells requires the degradation of cyclin B. We show that cells encounter and tolerate variations in the abundance of securin or cyclin B. This makes the concurrent onset of securin and cyclin B degradation insufficient to guarantee that early anaphase events occur in the correct order. We uncover that the timing of chromosome splitting is not determined by reaching a fixed securin level, but that this level adapts to the securin degradation kinetics. In conjunction with securin and cyclin B competing for degradation during anaphase, this provides robustness to the temporal order of anaphase events. Our work reveals how parallel cell-cycle pathways can be temporally coordinated despite variability in protein concentrations.

In vitro iCLIP-based modeling uncovers how the splicing factor U2AF2 relies on regulation by cofactors.

Alternative splicing generates distinct mRNA isoforms and is crucial for proteome diversity in eukaryotes. The RNA-binding protein (RBP) U2AF2 is central to splicing decisions, as it recognizes 3' splice sites and recruits the spliceosome. We establish "in vitro iCLIP" experiments, in which recombinant RBPs are incubated with long transcripts, to study how U2AF2 recognizes RNA sequences and how this is modulated by trans-acting RBPs. We measure U2AF2 affinities at hundreds of binding sites and compare in vitro and in vivo binding landscapes by mathematical modeling. We find that trans-acting RBPs extensively regulate U2AF2 binding in vivo, including enhanced recruitment to 3' splice sites and clearance of introns. Using machine learning, we identify and experimentally validate novel trans-acting RBPs (including FUBP1, CELF6, and PCBP1) that modulate U2AF2 binding and affect splicing outcomes. Our study offers a blueprint for the high-throughput characterization of in vitro mRNP assembly and in vivo splicing regulation.

Exon Definition Facilitates Reliable Control of Alternative Splicing in the RON Proto-Oncogene.

Alternative splicing is a key step in eukaryotic gene expression that allows for the production of multiple transcript and protein isoforms from the same gene. Even though splicing is perturbed in many diseases, we currently lack insights into regulatory mechanisms promoting its precision and efficiency. We analyze high-throughput mutagenesis data obtained for an alternatively spliced exon in the proto-oncogene RON and determine the functional units that control this splicing event. Using mathematical modeling of distinct splicing mechanisms, we show that alternative splicing is based in RON on a so-called "exon definition" mechanism. Here, the recognition of the adjacent exons by the spliceosome is required for removal of an intron. We use our model to analyze the differences between the exon and intron definition scenarios and find that exon definition prevents the accumulation of deleterious, partially spliced retention products during alternative splicing regulation. Furthermore, it modularizes splicing control, as multiple regulatory inputs are integrated into a common net input, irrespective of the location and nature of the corresponding cis-regulatory elements in the pre-messenger RNA. Our analysis suggests that exon definition promotes robust and reliable splicing outcomes in RON splicing.

Estrogen-dependent control and cell-to-cell variability of transcriptional bursting.

Cellular decision-making and environmental adaptation are dependent upon a heterogeneous response of gene expression to external cues. Heterogeneity arises in transcription from random switching between transcriptionally active and inactive states, resulting in bursts of RNA synthesis. Furthermore, the cellular state influences the competency of transcription, thereby globally affecting gene expression in a cell-specific manner. We determined how external stimuli interplay with cellular state to modulate the kinetics of bursting. To this end, single-cell dynamics of nascent transcripts were monitored at the endogenous estrogen-responsive GREB1 locus. Stochastic modeling of gene expression implicated a two-state promoter model in which the estrogen stimulus modulates the frequency of transcriptional bursting. The cellular state affects transcriptional dynamics by altering initiation and elongation kinetics and acts globally, as GREB1 alleles in the same cell correlate in their transcriptional output. Our results suggest that cellular state strongly affects the first step of the central dogma of gene expression, to promote heterogeneity in the transcriptional output of isogenic cells.

Cell-specific responses to the cytokine TGFß are determined by variability in protein levels.

The cytokine TGFß provides important information during embryonic development, adult tissue homeostasis, and regeneration. Alterations in the cellular response to TGFß are involved in severe human diseases. To understand how cells encode the extracellular input and transmit its information to elicit appropriate responses, we acquired quantitative time-resolved measurements of pathway activation at the single-cell level. We established dynamic time warping to quantitatively compare signaling dynamics of thousands of individual cells and described heterogeneous single-cell responses by mathematical modeling. Our combined experimental and theoretical study revealed that the response to a given dose of TGFß is determined cell specifically by the levels of defined signaling proteins. This heterogeneity in signaling protein expression leads to decomposition of cells into classes with qualitatively distinct signaling dynamics and phenotypic outcome. Negative feedback regulators promote heterogeneous signaling, as a SMAD7 knock-out specifically affected the signal duration in a subpopulation of cells. Taken together, we propose a quantitative framework that allows predicting and testing sources of cellular signaling heterogeneity.

Sharing of Phosphatases Promotes Response Plasticity in Phosphorylation Cascades.

Sharing of positive or negative regulators between multiple targets is frequently observed in cellular signaling cascades. For instance, phosphatase sharing between multiple kinases is ubiquitous within the MAPK pathway. Here we investigate how such phosphatase sharing could shape robustness and evolvability of the phosphorylation cascade. Through modeling and evolutionary simulations, we demonstrate that 1) phosphatase sharing dramatically increases robustness of a bistable MAPK response, and 2) phosphatase-sharing cascades evolve faster than nonsharing cascades. This faster evolution is particularly pronounced when evolving from a monostable toward a bistable phenotype, whereas the transition speed of a population from a bistable to monostable response is not affected by phosphatase sharing. This property may enable the phosphatase-sharing design to adapt better in a changing environment. Analysis of the respective mutational landscapes reveal that phosphatase sharing reduces the number of limiting mutations required for transition from monostable to bistable responses, hence facilitating a faster transition to such response types. Taken together, using MAPK cascade as an example, our study offers a general theoretical framework to explore robustness and evolutionary plasticity of signal transduction cascades.

Modelling Systemic Iron Regulation during Dietary Iron Overload and Acute Inflammation: Role of Hepcidin-Independent Mechanisms.

Systemic iron levels must be maintained in physiological concentrations to prevent diseases associated with iron deficiency or iron overload. A key role in this process plays ferroportin, the only known mammalian transmembrane iron exporter, which releases iron from duodenal enterocytes, hepatocytes, or iron-recycling macrophages into the blood stream. Ferroportin expression is tightly controlled by transcriptional and post-transcriptional mechanisms in response to hypoxia, iron deficiency, heme iron and inflammatory cues by cell-autonomous and systemic mechanisms. At the systemic level, the iron-regulatory hormone hepcidin is released from the liver in response to these cues, binds to ferroportin and triggers its degradation. The relative importance of individual ferroportin control mechanisms and their interplay at the systemic level is incompletely understood. Here, we built a mathematical model of systemic iron regulation. It incorporates the dynamics of organ iron pools as well as regulation by the hepcidin/ferroportin system. We calibrated and validated the model with time-resolved measurements of iron responses in mice challenged with dietary iron overload and/or inflammation. The model demonstrates that inflammation mainly reduces the amount of iron in the blood stream by reducing intracellular ferroportin transcription, and not by hepcidin-dependent ferroportin protein destabilization. In contrast, ferroportin regulation by hepcidin is the predominant mechanism of iron homeostasis in response to changing iron diets for a big range of dietary iron contents. The model further reveals that additional homeostasis mechanisms must be taken into account at very high dietary iron levels, including the saturation of intestinal uptake of nutritional iron and the uptake of circulating, non-transferrin-bound iron, into liver. Taken together, our model quantitatively describes systemic iron metabolism and generated experimentally testable predictions for additional ferroportin-independent homeostasis mechanisms.

Correlated receptor transport processes buffer single-cell heterogeneity.

Cells typically vary in their response to extracellular ligands. Receptor transport processes modulate ligand-receptor induced signal transduction and impact the variability in cellular responses. Here, we quantitatively characterized cellular variability in erythropoietin receptor (EpoR) trafficking at the single-cell level based on live-cell imaging and mathematical modeling. Using ensembles of single-cell mathematical models reduced parameter uncertainties and showed that rapid EpoR turnover, transport of internalized EpoR back to the plasma membrane, and degradation of Epo-EpoR complexes were essential for receptor trafficking. EpoR trafficking dynamics in adherent H838 lung cancer cells closely resembled the dynamics previously characterized by mathematical modeling in suspension cells, indicating that dynamic properties of the EpoR system are widely conserved. Receptor transport processes differed by one order of magnitude between individual cells. However, the concentration of activated Epo-EpoR complexes was less variable due to the correlated kinetics of opposing transport processes acting as a buffering system.

A multi-scale model of hepcidin promoter regulation reveals factors controlling systemic iron homeostasis.

Systemic iron homeostasis involves a negative feedback circuit in which the expression level of the peptide hormone hepcidin depends on and controls the iron blood levels. Hepcidin expression is regulated by the BMP6/SMAD and IL6/STAT signaling cascades. Deregulation of either pathway causes iron-related diseases such as hemochromatosis or anemia of inflammation. We quantitatively analyzed how BMP6 and IL6 control hepcidin expression. Transcription factor (TF) phosphorylation and reporter gene expression were measured under co-stimulation conditions, and the promoter was perturbed by mutagenesis. Using mathematical modeling, we systematically analyzed potential mechanisms of cooperative and competitive promoter regulation by the transcription factors, and experimentally validated the model predictions. Our results reveal that hepcidin cross-regulation primarily occurs by combinatorial transcription factor binding to the promoter, whereas signaling crosstalk is insignificant. We find that the presence of two BMP-responsive elements enhances the steepness of the promoter response towards the iron-sensing BMP signaling axis, which promotes iron homeostasis in vivo. IL6 co-stimulation reduces the promoter sensitivity towards the BMP signal, because the SMAD and STAT transcription factors compete for recruiting RNA polymerase to the transcription start site. This may explain why inflammatory signals disturb iron homeostasis in anemia of inflammation. Taken together, our results reveal why the iron homeostasis circuit is sensitive to perturbations implicated in disease.

Determinants of cell-to-cell variability in protein kinase signaling.

Cells reliably sense environmental changes despite internal and external fluctuations, but the mechanisms underlying robustness remain unclear. We analyzed how fluctuations in signaling protein concentrations give rise to cell-to-cell variability in protein kinase signaling using analytical theory and numerical simulations. We characterized the dose-response behavior of signaling cascades by calculating the stimulus level at which a pathway responds ('pathway sensitivity') and the maximal activation level upon strong stimulation. Minimal kinase cascades with gradual dose-response behavior show strong variability, because the pathway sensitivity and the maximal activation level cannot be simultaneously invariant. Negative feedback regulation resolves this trade-off and coordinately reduces fluctuations in the pathway sensitivity and maximal activation. Feedbacks acting at different levels in the cascade control different aspects of the dose-response curve, thereby synergistically reducing the variability. We also investigated more complex, ultrasensitive signaling cascades capable of switch-like decision making, and found that these can be inherently robust to protein concentration fluctuations. We describe how the cell-to-cell variability of ultrasensitive signaling systems can be actively regulated, e.g., by altering the expression of phosphatase(s) or by feedback/feedforward loops. Our calculations reveal that slow transcriptional negative feedback loops allow for variability suppression while maintaining switch-like decision making. Taken together, we describe design principles of signaling cascades that promote robustness. Our results may explain why certain signaling cascades like the yeast pheromone pathway show switch-like decision making with little cell-to-cell variability.

Robustness of signal transduction pathways.

Signal transduction pathways transduce information about the outside of the cell to the nucleus, regulating gene expression and cell fate. To reliably inform the cell about its surroundings, information transfer has to be robust against typical perturbation that a cell experiences. Robustness of several mammalian signaling pathways has been studied recently by quantitative experimentation and using mathematical modeling. Here, we review these studies, and describe the emerging concepts of robustness and the underlying mechanisms.

Negative feedback in the bone morphogenetic protein 4 (BMP4) synexpression group governs its dynamic signaling range and canalizes development.

What makes embryogenesis a robust and canalized process is an important question in developmental biology. A bone morphogenetic protein (BMP) morphogen gradient plays a key role in embryonic development, and we are beginning to understand how the self-regulating properties of its signaling circuitry ensure robust embryonic patterning. An unexplored question is why the BMP signaling circuit is organized as a modular synexpression group, with a prevalence of feedback inhibitors. Here, we provide evidence from direct experimentation and mathematical modeling that the synexpressed feedback inhibitors BAMBI, SMAD6, and SMAD7 (i) expand the dynamic BMP signaling range essential for proper embryonic patterning and (ii) reduce interindividual phenotypic and molecular variability in Xenopus embryos. Thereby, negative feedback linearizes signaling responses and confers robust patterning, thus promoting canalized development. The presence of negative feedback inhibitors in other growth factor synexpression groups suggests that these properties may constitute a general principle.

Reverse engineering a hierarchical regulatory network downstream of oncogenic KRAS.

RAS mutations are highly relevant for progression and therapy response of human tumours, but the genetic network that ultimately executes the oncogenic effects is poorly understood. Here, we used a reverse-engineering approach in an ovarian cancer model to reconstruct KRAS oncogene-dependent cytoplasmic and transcriptional networks from perturbation experiments based on gene silencing and pathway inhibitor treatments. We measured mRNA and protein levels in manipulated cells by microarray, RT-PCR and western blot analysis, respectively. The reconstructed model revealed complex interactions among the transcriptional and cytoplasmic components, some of which were confirmed by double pertubation experiments. Interestingly, the transcription factors decomposed into two hierarchically arranged groups. To validate the model predictions, we analysed growth parameters and transcriptional deregulation in the KRAS-transformed epithelial cells. As predicted by the model, we found two functional groups among the selected transcription factors. The experiments thus confirmed the predicted hierarchical transcription factor regulation and showed that the hierarchy manifests itself in downstream gene expression patterns and phenotype.

Modeling Cellular Signaling Variability Based on Single-Cell Data: The TGFß-SMAD Signaling Pathway.

Nongenetic heterogeneity is key to cellular decisions, as even genetically identical cells respond in very different ways to the same external stimulus, e.g., during cell differentiation or therapeutic treatment of disease. Strong heterogeneity is typically already observed at the level of signaling pathways that are the first sensors of external inputs and transmit information to the nucleus where decisions are made. Since heterogeneity arises from random fluctuations of cellular components, mathematical models are required to fully describe the phenomenon and to understand the dynamics of heterogeneous cell populations. Here, we review the experimental and theoretical literature on cellular signaling heterogeneity, with special focus on the TGFß/SMAD signaling pathway.

Position-dependent effects of RNA-binding proteins in the context of co-transcriptional splicing.

Alternative splicing is an important step in eukaryotic mRNA pre-processing which increases the complexity of gene expression programs, but is frequently altered in disease. Previous work on the regulation of alternative splicing has demonstrated that splicing is controlled by RNA-binding proteins (RBPs) and by epigenetic DNA/histone modifications which affect splicing by changing the speed of polymerase-mediated pre-mRNA transcription. The interplay of these different layers of splicing regulation is poorly understood. In this paper, we derived mathematical models describing how splicing decisions in a three-exon gene are made by combinatorial spliceosome binding to splice sites during ongoing transcription. We additionally take into account the effect of a regulatory RBP and find that the RBP binding position within the sequence is a key determinant of how RNA polymerase velocity affects splicing. Based on these results, we explain paradoxical observations in the experimental literature and further derive rules explaining why the same RBP can act as inhibitor or activator of cassette exon inclusion depending on its binding position. Finally, we derive a stochastic description of co-transcriptional splicing regulation at the single-cell level and show that splicing outcomes show little noise and follow a binomial distribution despite complex regulation by a multitude of factors. Taken together, our simulations demonstrate the robustness of splicing outcomes and reveal that quantitative insights into kinetic competition of co-transcriptional events are required to fully understand this important mechanism of gene expression diversity.

High-throughput mutagenesis identifies mutations and RNA-binding proteins controlling CD19 splicing and CART-19 therapy resistance.

Following CART-19 immunotherapy for B-cell acute lymphoblastic leukaemia (B-ALL), many patients relapse due to loss of the cognate CD19 epitope. Since epitope loss can be caused by aberrant CD19 exon 2 processing, we herein investigate the regulatory code that controls CD19 splicing. We combine high-throughput mutagenesis with mathematical modelling to quantitatively disentangle the effects of all mutations in the region comprising CD19 exons 1-3. Thereupon, we identify ~200 single point mutations that alter CD19 splicing and thus could predispose B-ALL patients to developing CART-19 resistance. Furthermore, we report almost 100 previously unknown splice isoforms that emerge from cryptic splice sites and likely encode non-functional CD19 proteins. We further identify cis-regulatory elements and trans-acting RNA-binding proteins that control CD19 splicing (e.g., PTBP1 and SF3B4) and validate that loss of these factors leads to pervasive CD19 mis-splicing. Our dataset represents a comprehensive resource for identifying predictive biomarkers for CART-19 therapy.