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.

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.

Decoding a cancer-relevant splicing decision in the RON proto-oncogene using high-throughput mutagenesis.

Mutations causing aberrant splicing are frequently implicated in human diseases including cancer. Here, we establish a high-throughput screen of randomly mutated minigenes to decode the cis-regulatory landscape that determines alternative splicing of exon 11 in the proto-oncogene MST1R (RON). Mathematical modelling of splicing kinetics enables us to identify more than 1000 mutations affecting RON exon 11 skipping, which corresponds to the pathological isoform RON∆165. Importantly, the effects correlate with RON alternative splicing in cancer patients bearing the same mutations. Moreover, we highlight heterogeneous nuclear ribonucleoprotein H (HNRNPH) as a key regulator of RON splicing in healthy tissues and cancer. Using iCLIP and synergy analysis, we pinpoint the functionally most relevant HNRNPH binding sites and demonstrate how cooperative HNRNPH binding facilitates a splicing switch of RON exon 11. Our results thereby offer insights into splicing regulation and the impact of mutations on alternative splicing in cancer.

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.

Actin filament elasticity and retrograde flow shape the force-velocity relation of motile cells.

Cells migrate through a crowded environment during processes such as metastasis or wound healing, and must generate and withstand substantial forces. The cellular motility responses to environmental forces are represented by their force-velocity relation, which has been measured for fish keratocytes but remains unexplained. Even pN opposing forces slow down lamellipodium motion by three orders of magnitude. At larger opposing forces, the retrograde flow of the actin network accelerates until it compensates for polymerization, and cell motion stalls. Subsequently, the lamellipodium adapts to the stalled state. We present a mechanism quantitatively explaining the cell's force-velocity relation and its changes upon application of drugs that hinder actin polymerization or actomyosin-based contractility. Elastic properties of filaments, close to the lamellipodium leading edge, and retrograde flow shape the force-velocity relation. To our knowledge, our results shed new light on how these migratory responses are regulated, and on the mechanics and structure of the lamellipodium.

Modeling morphodynamic phenotypes and dynamic regimes of cell motion.

Many cellular processes and signaling pathways converge onto cell morphology and cell motion, which share important components. The mechanisms used for propulsion could also be responsible for shape changes, if they are capable of generating the rich observed variety of dynamic regimes. Additionally, the analysis of cell shape changes in space and time promises insight into the state of the cytoskeleton and signaling pathways controlling it. While this has been obvious for some time by now, little effort has been made to systematically and quantitatively explore this source of information. First pioneering experimental work revealed morphodynamic phenotypes which can be associated with dynamic regimes like oscillations and excitability. Here, we review the current state of modeling of morphodynamic phenotypes, the experimental results and discuss the ideas on the mechanisms driving shape changes which are suggested by modeling.

Active colloidal suspensions exhibit polar order under gravity.

Recently, the steady sedimentation profile of a dilute suspension of chemically powered colloids was studied experimentally [J. Palacci et al., Phys. Rev. Lett. 105, 088304 (2010)]. It was found that the sedimentation length increases quadratically with the swimming speed of the active Brownian particles. Here we investigate theoretically the sedimentation of self-propelled particles undergoing translational and rotational diffusion. We find that the measured increase of the sedimentation length is coupled to a partial alignment of the suspension with the mean swimming direction oriented against the gravitational field. We suggest realistic parameter values to observe this polar order. Furthermore, we find that the dynamics of the active suspension can be derived from a generalized free energy functional.

Modeling of protrusion phenotypes driven by the actin-membrane interaction.

We propose a mathematical model for simulating the leading-edge dynamics of a migrating cell from the interplay among elastic properties, architecture of the actin cytoskeleton, and the mechanics of the membrane. Our approach is based on the description of the length and attachment dynamics of actin filaments in the lamellipodium network. It is used to determine the total force exerted on the membrane at each position along the leading edge and at each time step. The model reproduces the marked state switches in protrusion morphodynamics found experimentally between epithelial cells in control conditions and cells expressing constitutively active Rac, a signaling molecule involved in the regulation of lamellipodium network assembly. The model also suggests a mechanistic explanation of experimental distortions in protrusion morphodynamics induced by deregulation of Arp2/3 and cofilin activity.

Leading-edge-gel coupling in lamellipodium motion.

We present a model for actin-based motility that combines the dynamics of the semiflexible region at the leading edge of the lamellipodium with actomyosin gel properties in the bulk described by the theory of active polar gels. We calculate the velocity of the lamellipodium determined by the interaction of the gel and adhesion with forces in the semiflexible region. The stationary concave force-velocity relation of the model reproduces experimental results. We suggest that it is determined by retrograde flow at small forces and gel formation and retrograde flow at large ones. The variety of dynamic regimes of the semiflexible region reproducing experimentally observed morphodynamics is conserved when we couple the leading edge to the gel.

Dynamic regimes and bifurcations in a model of actin-based motility.

Propulsion by actin polymerization is widely used in cell motility. Here, we investigate a model of the brush range of an actin gel close to a propelled object, describing the force generation and the dynamics of the propagation velocity. We find transitions between stable steady states and relaxation oscillations when the attachment rate of actin filaments to the obstacle is varied. The oscillations set in at small values of the attachment rate via a homoclinic bifurcation. A second transition from a stable steady state to relaxation oscillations, found for higher values of the attachment rate, occurs via a supercritical Hopf bifurcation. The behavior of the model near the second transition is similar that of a system undergoing a canard explosion. Consequently, we observe excitable dynamics also. The model further exhibits bistability between stationary states or stationary states and limit cycles. Therefore, the brush of actin filament ends appears to have a much richer dynamics than was assumed until now.

Activity dynamics in nonlocal interacting neural fields.

We study the activity of a synaptically coupled neuronal network consisting of an excitatory and an inhibitory layer with isotropic connections and nonlinear interactions. Using the mathematical model of Wilson and Cowan in two spatial dimensions, we first discuss a spatial hysteresis phenomenon. Then we analyze special traveling wave solutions with stationary shape. We establish existence conditions, derive analytic expressions of the particular solutions and their velocity, and finally present numerical simulations.