Investigation of early events in Fc epsilon RI-mediated signaling using a detailed mathematical model

bnglViz: online visualization of rule-based models

MOTIVATION

Rule-based modeling is a powerful method to describe and simulate interactions among multi-site molecules and multi-molecular species, accounting for the internal connectivity of molecules in chemical species. This modeling technique is implemented in BioNetGen software that is used by various tools and software frameworks, such as BioNetGen stand-alone software, NFSim simulation engine, Virtual Cell simulation and modeling framework, SmolDyn and PySB software tools. These tools exchange models using BioNetGen scripting language (BNGL). Until now, there was no online visualization of such rule-based models. Modelers and researchers reading the manuscripts describing rule-based models had to learn BNGL scripting or master one of these tools to understand the models.

RESULTS

Here, we introduce bnglViz, an online platform for visualizing BNGL files as graphical cartoons, empowering researchers to grasp the nuances of rule-based models swiftly and efficiently, and making the exploration of complex biological systems more accessible than ever before. The produced visualizations can be used as supplemental figures in publications or as a way to annotate BNGL models on web repositories.

AVAILABILITY AND IMPLEMENTATION

Available at https://bnglviz.github.io/.

Catalyst: Fast and flexible modeling of reaction networks

We introduce Catalyst.jl, a flexible and feature-filled Julia library for modeling and high-performance simulation of chemical reaction networks (CRNs). Catalyst supports simulating stochastic chemical kinetics (jump process), chemical Langevin equation (stochastic differential equation), and reaction rate equation (ordinary differential equation) representations for CRNs. Through comprehensive benchmarks, we demonstrate that Catalyst simulation runtimes are often one to two orders of magnitude faster than other popular tools. More broadly, Catalyst acts as both a domain-specific language and an intermediate representation for symbolically encoding CRN models as Julia-native objects. This enables a pipeline of symbolically specifying, analyzing, and modifying CRNs; converting Catalyst models to symbolic representations of concrete mathematical models; and generating compiled code for numerical solvers. Leveraging ModelingToolkit.jl and Symbolics.jl, Catalyst models can be analyzed, simplified, and compiled into optimized representations for use in numerical solvers. Finally, we demonstrate Catalyst's broad extensibility and composability by highlighting how it can compose with a variety of Julia libraries, and how existing open-source biological modeling projects have extended its intermediate representation.

Spatial Stochastic Model of the Pre-B Cell Receptor

Survival and proliferation of immature B lymphocytes requires expression and tonic signaling of the pre-B cell receptor (pre-BCR). This low level, ligand-independent signaling is likely achieved through frequent, but short-lived, homo interactions. Tonic signaling is also central in the pathology of precursor B acute lymphoblastic leukemia (B-ALL). In order to understand how repeated, transient events can lead to sustained signaling and to assess the impact of receptor accumulation induced by the membrane landscape, we developed a spatial stochastic model of receptor aggregation and downstream signaling events. Our rule- and agent-based model builds on previous mature BCR signaling models and incorporates novel parameters derived from single particle tracking of pre-BCR on surfaces of two different B-ALL cell lines, 697 and Nalm6. Live cell tracking of receptors on the two cell lines revealed characteristic differences in their dimer dissociation rates and diffusion coefficients. We report here that these differences affect pre-BCR aggregation and consequent signal initiation events. Receptors on Nalm6 cells, which have a lower off-rate and lower diffusion coefficient, more frequently form higher order oligomers than pre-BCR on 697 cells, resulting in higher levels of downstream phosphorylation in the Nalm6 cell line.

Comparison of rule- and ordinary differential equation-based dynamic model of DARPP-32 signalling network

Dynamic modelling has considerably improved our understanding of complex molecular mechanisms. Ordinary differential equations (ODEs) are the most detailed and popular approach to modelling the dynamics of molecular systems. However, their application in signalling networks, characterised by multi-state molecular complexes, can be prohibitive. Contemporary modelling methods, such as rule- based (RB) modelling, have addressed these issues. The advantages of RB modelling over ODEs have been presented and discussed in numerous reviews. In this study, we conduct a direct comparison of the time courses of a molecular system founded on the same reaction network but encoded in the two frameworks. To make such a comparison, a set of reactions that underlie an ODE model was manually encoded in the Kappa language, one of the RB implementations. A comparison of the models was performed at the level of model specification and dynamics, acquired through model simulations. In line with previous reports, we confirm that the Kappa model recapitulates the general dynamics of its ODE counterpart with minor differences. These occur when molecules have multiple sites binding the same interactor. Furthermore, activation of these molecules in the RB model is slower than in the ODE one. As reported for other molecular systems, we find that, also for the DARPP-32 reaction network, the RB representation offers a more expressive and flexible syntax that facilitates access to fine details of the model, easing model reuse. In parallel with these analyses, we report a refactored model of the DARPP-32 interaction network that can serve as a canvas for the development of more complex dynamic models to study this important molecular system.

Docking of Syk to FcεRI is enhanced by Lyn but limited in duration by SHIP1

The high-affinity immunoglobulin E (IgE) receptor, FcεRI, is the primary immune receptor found on mast cells and basophils. Signal initiation is classically attributed to phosphorylation of FcεRI β− and γ-subunits by the Src family kinase (SFK) Lyn, followed by the recruitment and activation of the tyrosine kinase Syk. FcεRI signaling is tuned by the balance between Syk-driven positive signaling and the engagement of inhibitory molecules, including SHIP1. Here, we investigate the mechanistic contributions of Lyn, Syk, and SHIP1 to the formation of the FcεRI signalosome. Using Lyn-deficient RBL-2H3 mast cells, we found that another SFK can weakly monophosphorylate the γ-subunit, yet Syk still binds the incompletely phosphorylated immunoreceptor tyrosine-based activation motifs (ITAMs). Once recruited, Syk further enhances γ-phosphorylation to propagate signaling. In contrast, the loss of SHIP1 recruitment indicates that Lyn is required for phosphorylation of the β-subunit. We demonstrate two noncanonical Syk binding modes, trans γ-bridging and direct β-binding, that can support signaling when SHIP1 is absent. Using single particle tracking, we reveal a novel role of SHIP1 in regulating Syk activity, where the presence of SHIP1 in the signaling complex acts to increase the Syk:receptor off-rate. These data suggest that the composition and dynamics of the signalosome modulate immunoreceptor signaling activities.

Modularisation of published and novel models toward a complex KIR2DL4 pathway in pbNK cell

KIR2DL4 is an interesting receptor expressed on the peripheral blood natural killer (pbNK) cell as it can be either activating or inhibitory depending on the amino acid residues in the domain. This model uses mathematical modelling to investigate the downstream effects of natural killer cells’ activation (KIR2DL4) receptor after stimulation by key ligand (HLA-G) on pbNK cells. Development of this large pathway is based on a comprehensive qualitative description of pbNKs’ intracellular signalling pathways leading to chemokine and cytotoxin secretion, obtained from the KEGG database (https://www.genome.jp/pathway/hsa04650). From this qualitative description we built a quantitative model for the pathway, reusing existing curated models where possible and implementing new models as needed. This model employs a composite approach for generating modular models. The approach allows for the construction of large-scale complex model by combining component of sub-models that can be modified individually. This large pathway consists of two published sub-models; the Ca²⁺ model and the NFAT model, and a newly built FCεRIγ sub-model. The full pathway was fitted to published dataset and fitted well to one of two secreted cytokines. The model can be used to predict the production of IFNγ and TNFα cytokines.•

Development of pathway and mathematical model

•

Reusing existing curated models and implementing new models

•

Model optimization and analysis

A multiscale model of the regulation of aquaporin 2 recycling

The response of cells to their environment is driven by a variety of proteins and messenger molecules. In eukaryotes, their distribution and location in the cell are regulated by the vesicular transport system. The transport of aquaporin 2 between membrane and storage region is a crucial part of the water reabsorption in renal principal cells, and its malfunction can lead to Diabetes insipidus. To understand the regulation of this system, we aggregated pathways and mechanisms from literature and derived three models in a hypothesis-driven approach. Furthermore, we combined the models to a single system to gain insight into key regulatory mechanisms of Aquaporin 2 recycling. To achieve this, we developed a multiscale computational framework for the modeling and simulation of cellular systems. The analysis of the system rationalizes that the compartmentalization of cAMP in renal principal cells is a result of the protein kinase A signalosome and can only occur if specific cellular components are observed in conjunction. Endocytotic and exocytotic processes are inherently connected and can be regulated by the same protein kinase A signal.

BiPSim: a flexible and generic stochastic simulator for polymerization processes

Detailed whole-cell modeling requires an integration of heterogeneous cell processes having different modeling formalisms, for which whole-cell simulation could remain tractable. Here, we introduce BiPSim, an open-source stochastic simulator of template-based polymerization processes, such as replication, transcription and translation. BiPSim combines an efficient abstract representation of reactions and a constant-time implementation of the Gillespie’s Stochastic Simulation Algorithm (SSA) with respect to reactions, which makes it highly efficient to simulate large-scale polymerization processes stochastically. Moreover, multi-level descriptions of polymerization processes can be handled simultaneously, allowing the user to tune a trade-off between simulation speed and model granularity. We evaluated the performance of BiPSim by simulating genome-wide gene expression in bacteria for multiple levels of granularity. Finally, since no cell-type specific information is hard-coded in the simulator, models can easily be adapted to other organismal species. We expect that BiPSim should open new perspectives for the genome-wide simulation of stochastic phenomena in biology.

CLUE: Exact maximal reduction of kinetic models by constrained lumping of differential equations

MOTIVATION

Detailed mechanistic models of biological processes can pose significant challenges for analysis and parameter estimations due to the large number of equations used to track the dynamics of all distinct configurations in which each involved biochemical species can be found. Model reduction can help tame such complexity by providing a lower-dimensional model in which each macro-variable can be directly related to the original variables.

RESULTS

We present CLUE, an algorithm for exact model reduction of systems of polynomial differential equations by constrained linear lumping. It computes the smallest dimensional reduction as a linear mapping of the state space such that the reduced model preserves the dynamics of user-specified linear combinations of the original variables. Even though CLUE works with nonlinear differential equations, it is based on linear algebra tools, which makes it applicable to high-dimensional models. Using case studies from the literature, we show how CLUE can substantially lower model dimensionality and help extract biologically intelligible insights from the reduction.

AVAILABILITY

An implementation of the algorithm and relevant resources to replicate the experiments herein reported are freely available for download at https://github.com/pogudingleb/CLUE.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Applications of personalised signalling network models in precision oncology

As our ability to provide in-depth, patient-specific characterisation of the molecular alterations within tumours rapidly improves, it is becoming apparent that new approaches will be required to leverage the power of this data and derive the full benefit for each individual patient. Systems biology approaches are beginning to emerge within this field as a potential method of incorporating large volumes of network level data and distilling a coherent, clinically-relevant prediction of drug response. However, the initial promise of this developing field is yet to be realised. Here we argue that in order to develop these precise models of individual drug response and tailor treatment accordingly, we will need to develop mathematical models capable of capturing both the dynamic nature of drug-response signalling networks and key patient-specific information such as mutation status or expression profiles. We also review the modelling approaches commonly utilised within this field, and outline recent examples of their use in furthering the application of systems biology for a precision medicine approach to cancer treatment.

Parallel Tempering with Lasso for model reduction in systems biology

Systems Biology models reveal relationships between signaling inputs and observable molecular or cellular behaviors. The complexity of these models, however, often obscures key elements that regulate emergent properties. We use a Bayesian model reduction approach that combines Parallel Tempering with Lasso regularization to identify minimal subsets of reactions in a signaling network that are sufficient to reproduce experimentally observed data. The Bayesian approach finds distinct reduced models that fit data equivalently. A variant of this approach that uses Lasso to perform selection at the level of reaction modules is applied to the NF-κB signaling network to test the necessity of feedback loops for responses to pulsatile and continuous pathway stimulation. Taken together, our results demonstrate that Bayesian parameter estimation combined with regularization can isolate and reveal core motifs sufficient to explain data from complex signaling systems.

Cells respond to diverse environmental cues using complex networks of interacting proteins and other biomolecules. Mathematical and computational models have become invaluable tools to understand these networks and make informed predictions to rationally perturb cell behavior. However, the complexity of detailed models that try to capture all known biochemical elements of signaling networks often makes it difficult to determine the key regulatory elements that are responsible for specific cell behaviors. Here, we present a Bayesian computational approach, PTLasso, to automatically extract minimal subsets of detailed models that are sufficient to explain experimental data. The method simultaneously calibrates and reduces models, and the Bayesian approach samples globally, allowing us to find alternate mechanistic explanations for the data if present. We demonstrate the method on both synthetic and real biological data and show that PTLasso is an effective method to isolate distinct parts of a larger signaling model that are sufficient for specific data.

Pleione: A tool for statistical and multi-objective calibration of Rule-based models

Mathematical models based on Ordinary Differential Equations (ODEs) are frequently used to describe and simulate biological systems. Nevertheless, such models are often difficult to understand. Unlike ODE models, Rule-Based Models (RBMs) utilise formal language to describe reactions as a cumulative number of statements that are easier to understand and correct. They are also gaining popularity because of their conciseness and simulation flexibility. However, RBMs generally lack tools to perform further analysis that requires simulation. This situation arises because exact and approximate simulations are computationally intensive. Translating RBMs into ODEs is commonly used to reduce simulation time, but this technique may be prohibitive due to combinatorial explosion. Here, we present the software called Pleione to calibrate RBMs. Parameter calibration is essential given the incomplete experimental determination of reaction rates and the goal of using models to reproduce experimental data. The software distributes stochastic simulations and calculations and incorporates equivalence tests to determine the fitness of RBMs compared with data. The primary features of Pleione were thoroughly tested on a model of gene regulation in Escherichia coli. Pleione yielded satisfactory results regarding calculation time and error reduction for multiple simulators, models, parameter search strategies, and computing infrastructures.

PyBioNetFit and the Biological Property Specification Language

In systems biology modeling, important steps include model parameterization, uncertainty quantification, and evaluation of agreement with experimental observations. To help modelers perform these steps, we developed the software PyBioNetFit, which in addition supports checking models against known system properties and solving design problems. PyBioNetFit introduces Biological Property Specification Language (BPSL) for the formal declaration of system properties. BPSL allows qualitative data to be used alone or in combination with quantitative data. PyBioNetFit performs parameterization with parallelized metaheuristic optimization algorithms that work directly with existing model definition standards: BioNetGen Language (BNGL) and Systems Biology Markup Language (SBML). We demonstrate PyBioNetFit's capabilities by solving various example problems, including the challenging problem of parameterizing a 153-parameter model of cell cycle control in yeast based on both quantitative and qualitative data. We demonstrate the model checking and design applications of PyBioNetFit and BPSL by analyzing a model of targeted drug interventions in autophagy signaling.

Combinatorial diversity of Syk recruitment driven by its multivalent engagement with FcεRIγ

Syk/Zap70 family kinases are essential for signaling via multichain immune-recognition receptors such as tetrameric (αβγ2) FcεRI. Syk activation is generally attributed to cis binding of its tandem SH2 domains to dual phosphotyrosines within FcεRIγ-ITAMs (immunoreceptor tyrosine-based activation motifs). However, the mechanistic details of Syk docking on γ homodimers are unresolved. Here, we estimate that multivalent interactions for WT Syk improve cis-oriented binding by three orders of magnitude. We applied molecular dynamics (MD), hybrid MD/worm-like chain polymer modeling, and live cell imaging to evaluate relative binding and signaling output for all possible cis and trans Syk-FcεRIγ configurations. Syk binding is likely modulated during signaling by autophosphorylation on Y130 in interdomain A, since a Y130E phosphomimetic form of Syk is predicted to lead to reduced helicity of interdomain A and alter Syk's bias for cis binding. Experiments in reconstituted γ-KO cells, whose γ subunits are linked by disulfide bonds, as well as in cells expressing monomeric ITAM or hemITAM γ-chimeras, support model predictions that short distances between γ ITAM pairs are required for trans docking. We propose that the full range of docking configurations improves signaling efficiency by expanding the combinatorial possibilities for Syk recruitment, particularly under conditions of incomplete ITAM phosphorylation.

Efficient anticorrelated variance reduction for stochastic simulation of biochemical reactions

We investigate the computational challenge of improving the accuracy of the stochastic simulation estimation by inducing negative correlation through the anticorrelated variance reduction technique. A direct application of the technique to the stochastic simulation algorithm (SSA), employing the inverse transformation, is not efficient for simulating large networks because its computational cost is similar to the sum of independent simulation runs. We propose in this study a new algorithm that employs the propensity bounds of reactions, introduced recently in their rejection-based SSA, to correlate and synchronise the trajectories during the simulation. The selection of reaction firings by our approach is exact due to the rejection-based mechanism. In addition, by applying the anticorrelated variance technique to select reaction firings, our approach can induce substantial correlation between realisations, hence reducing the variance of the estimator. The computational advantage of our rejection-based approach in comparison with the traditional inverse transformation is that it only needs to maintain a single data structure storing propensity bounds of reactions, which is updated infrequently, hence achieving better performance.

MCell-R: A Particle-Resolution Network-Free Spatial Modeling Framework

Spatial heterogeneity can have dramatic effects on the biochemical networks that drive cell regulation and decision-making. For this reason, a number of methods have been developed to model spatial heterogeneity and incorporated into widely used modeling platforms. Unfortunately, the standard approaches for specifying and simulating chemical reaction networks become untenable when dealing with multistate, multicomponent systems that are characterized by combinatorial complexity. To address this issue, we developed MCell-R, a framework that extends the particle-based spatial Monte Carlo simulator, MCell, with the rule-based model specification and simulation capabilities provided by BioNetGen and NFsim. The BioNetGen syntax enables the specification of biomolecules as structured objects whose components can have different internal states that represent such features as covalent modification and conformation and which can bind components of other molecules to form molecular complexes. The network-free simulation algorithm used by NFsim enables efficient simulation of rule-based models even when the size of the network implied by the biochemical rules is too large to enumerate explicitly, which frequently occurs in detailed models of biochemical signaling. The result is a framework that can efficiently simulate systems characterized by combinatorial complexity at the level of spatially resolved individual molecules over biologically relevant time and length scales.

Multicellular Models Bridging Intracellular Signaling and Gene Transcription to Population Dynamics

Cell signaling and gene transcription occur at faster time scales compared to cellular death, division, and evolution. Bridging these multiscale events in a model is computationally challenging. We introduce a framework for the systematic development of multiscale cell population models. Using message passing interface (MPI) parallelism, the framework creates a population model from a single-cell biochemical network model. It launches parallel simulations on a single-cell model and treats each stand-alone parallel process as a cell object. MPI mediates cell-to-cell and cell-to-environment communications in a server-client fashion. In the framework, model-specific higher level rules link the intracellular molecular events to cellular functions, such as death, division, or phenotype change. Cell death is implemented by terminating a parallel process, while cell division is carried out by creating a new process (daughter cell) from an existing one (mother cell). We first demonstrate these capabilities by creating two simple example models. In one model, we consider a relatively simple scenario where cells can evolve independently. In the other model, we consider interdependency among the cells, where cellular communication determines their collective behavior and evolution under a temporally evolving growth condition. We then demonstrate the framework’s capability by simulating a full-scale model of bacterial quorum sensing, where the dynamics of a population of bacterial cells is dictated by the intercellular communications in a time-evolving growth environment.

A spatio-temporal model reveals self-limiting FcɛRI cross-linking by multivalent antigens

Aggregation of cell surface receptor proteins by multivalent antigens is an essential early step for immune cell signalling. A number of experimental and modelling studies in the past have investigated multivalent ligand-mediated aggregation of IgE receptors (FcɛRI) in the plasma membrane of mast cells. However, understanding of the mechanisms of FcɛRI aggregation remains incomplete. Experimental reports indicate that FcɛRI forms relatively small and finite-sized clusters when stimulated by a multivalent ligand. By contrast, modelling studies have shown that receptor cross-linking by a trivalent ligand may lead to the formation of large receptor superaggregates that may potentially give rise to hyperactive cellular responses. In this work, we have developed a Brownian dynamics-based spatio-temporal model to analyse FcɛRI aggregation by a trivalent antigen. Unlike the existing models, which implemented non-spatial simulation approaches, our model explicitly accounts for the coarse-grained site-specific features of the multivalent species (molecules and complexes). The model incorporates membrane diffusion, steric collisions and sub-nanometre-scale site-specific interaction of the time-evolving species of arbitrary structures. Using the model, we investigated temporal evolution of the species and their diffusivities. Consistent with a recent experimental report, our model predicted sharp decay in species mobility in the plasma membrane in response receptor cross-linking by a multivalent antigen. We show that, due to such decay in the species mobility, post-stimulation receptor aggregation may become self-limiting. Our analysis reveals a potential regulatory mechanism suppressing hyperactivation of immune cells in response to multivalent antigens.

A Critical Comparison of Rejection-Based Algorithms for Simulation of Large Biochemical Reaction Networks

The rejection-based simulation technique has been applying to improve the computational efficiency of the stochastic simulation algorithm (SSA) in simulating large reaction networks, which are required for a thorough understanding of biological systems. We compare two recently proposed simulation methods, namely the composition-rejection algorithm (SSA-CR) and the rejection-based SSA (RSSA), aiming for this purpose. We discuss the right interpretation of the rejection-based technique used in these algorithms in order to make an informed choice when dealing with different aspects of biochemical networks. We provide the theoretical analysis as well as the detailed runtime comparison of these algorithms on concrete biological models. We highlight important factors that are omitted in previous analysis of these algorithms. The numerical comparison shows that for reaction networks where the search cost is expensive then SSA-CR is more efficient, and for reaction networks where the update cost is dominant, often the case in practice, then RSSA should be the choice.

An Overview of Network-Based and -Free Approaches for Stochastic Simulation of Biochemical Systems

Stochastic simulation has been widely used to model the dynamics of biochemical reaction networks. Several algorithms have been proposed that are exact solutions of the chemical master equation, following the work of Gillespie. These stochastic simulation approaches can be broadly classified into two categories: network-based and -free simulation. The network-based approach requires that the full network of reactions be established at the start, while the network-free approach is based on reaction rules that encode classes of reactions, and by applying rule transformations, it generates reaction events as they are needed without ever having to derive the entire network. In this study, we compare the efficiency and limitations of several available implementations of these two approaches. The results allow for an informed selection of the implementation and methodology for specific biochemical modeling applications.

Timescale Separation of Positive and Negative Signaling Creates History-Dependent Responses to IgE Receptor Stimulation

The high-affinity receptor for IgE expressed on the surface of mast cells and basophils interacts with antigens, via bound IgE antibody, and triggers secretion of inflammatory mediators that contribute to allergic reactions. To understand how past inputs (memory) influence future inflammatory responses in mast cells, a microfluidic device was used to precisely control exposure of cells to alternating stimulatory and non-stimulatory inputs. We determined that the response to subsequent stimulation depends on the interval of signaling quiescence. For shorter intervals of signaling quiescence, the second response is blunted relative to the first response, whereas longer intervals of quiescence induce an enhanced second response. Through an iterative process of computational modeling and experimental tests, we found that these memory-like phenomena arise from a confluence of rapid, short-lived positive signals driven by the protein tyrosine kinase Syk; slow, long-lived negative signals driven by the lipid phosphatase Ship1; and slower degradation of Ship1 co-factors. This work advances our understanding of mast cell signaling and represents a generalizable approach for investigating the dynamics of signaling systems.

Automated visualization of rule-based models

Frameworks such as BioNetGen, Kappa and Simmune use "reaction rules" to specify biochemical interactions compactly, where each rule specifies a mechanism such as binding or phosphorylation and its structural requirements. Current rule-based models of signaling pathways have tens to hundreds of rules, and these numbers are expected to increase as more molecule types and pathways are added. Visual representations are critical for conveying rule-based models, but current approaches to show rules and interactions between rules scale poorly with model size. Also, inferring design motifs that emerge from biochemical interactions is an open problem, so current approaches to visualize model architecture rely on manual interpretation of the model. Here, we present three new visualization tools that constitute an automated visualization framework for rule-based models: (i) a compact rule visualization that efficiently displays each rule, (ii) the atom-rule graph that conveys regulatory interactions in the model as a bipartite network, and (iii) a tunable compression pipeline that incorporates expert knowledge and produces compact diagrams of model architecture when applied to the atom-rule graph. The compressed graphs convey network motifs and architectural features useful for understanding both small and large rule-based models, as we show by application to specific examples. Our tools also produce more readable diagrams than current approaches, as we show by comparing visualizations of 27 published models using standard graph metrics. We provide an implementation in the open source and freely available BioNetGen framework, but the underlying methods are general and can be applied to rule-based models from the Kappa and Simmune frameworks also. We expect that these tools will promote communication and analysis of rule-based models and their eventual integration into comprehensive whole-cell models.

Interaction of DJ-1 with Lyn is essential for IgE-mediated stimulation of human mast cells

BACKGROUND

DJ-1 is a redox-sensitive protein with multiple roles in cell homeostasis, levels of which are altered in patients with mast cell (MC)-related disorders. However, whether DJ-1 can regulate human MC function is unknown.

OBJECTIVE

We sought to investigate the potential role of DJ-1 in the responses of human MCs to antigen stimulation.

METHODS

DJ-1 was silenced in human CD34⁺-derived MCs and in the LAD2 MC line by using lentiviral short hairpin RNA constructs. Release of β-hexosaminidase, prostaglandin D₂, and GM-CSF and changes in reactive oxygen species levels were measured after FcεRI engagement. Enzymatic assays, sucrose density gradient centrifugation, immunoprecipitation, dot and Western blotting, and confocal imaging were performed for signaling, cellular localization, and coassociation studies.

RESULTS

DJ-1 knockdown substantially reduced mediator release, as well as Lyn kinase and spleen tyrosine kinase activation and signaling through mechanisms that appeared largely unrelated to DJ-1 antioxidant activity. Following FcεRI activation, nonoxidized rather than oxidized DJ-1 translocated to lipid rafts, where it associated with Lyn, an interaction that appeared critical for maximal Lyn activation and initiation of signaling. Using purified recombinant proteins, we demonstrated that DJ-1 directly bound to Lyn but not to other Src kinases, and this interaction was specific for human but not mouse proteins. In addition, DJ-1 reduced Src homology 2 domain-containing phosphatase 2 phosphatase activity by scavenging reactive oxygen species, thus preventing spleen tyrosine kinase dephosphorylation and perpetuating MC signaling.

CONCLUSION

We demonstrate a novel role for DJ-1 in the early activation of Lyn by FcεRI, which is essential for human MC responses and provides the basis for an alternative target in allergic disease therapy.

Maximal aggregation of polynomial dynamical systems

Ordinary differential equations (ODEs) with polynomial derivatives are a fundamental tool for understanding the dynamics of systems across many branches of science, but our ability to gain mechanistic insight and effectively conduct numerical evaluations is critically hindered when dealing with large models. Here we propose an aggregation technique that rests on two notions of equivalence relating ODE variables whenever they have the same solution (backward criterion) or if a self-consistent system can be written for describing the evolution of sums of variables in the same equivalence class (forward criterion). A key feature of our proposal is to encode a polynomial ODE system into a finitary structure akin to a formal chemical reaction network. This enables the development of a discrete algorithm to efficiently compute the largest equivalence, building on approaches rooted in computer science to minimize basic models of computation through iterative partition refinements. The physical interpretability of the aggregation is shown on polynomial ODE systems for biochemical reaction networks, gene regulatory networks, and evolutionary game theory.

Stochastic simulation of biochemical reactions with partial-propensity and rejection-based approaches

We present in this paper a new exact algorithm for improving performance of exact stochastic simulation algorithm. The algorithm is developed on concepts of the partial-propensity and the rejection-based approaches. It factorizes the propensity bounds of reactions and groups factors by common reactant species for selecting next reaction firings. Our algorithm provides favorable computational advantages for simulating of biochemical reaction networks by reducing the cost for selecting the next reaction firing to scale with the number of chemical species and avoiding expensive propensity updates during the simulation. We present the details of our new algorithm and benchmark it on concrete biological models to demonstrate its applicability and efficiency.

Simple biochemical networks allow accurate sensing of multiple ligands with a single receptor

Cells use surface receptors to estimate concentrations of external ligands. Limits on the accuracy of such estimations have been well studied for pairs of ligand and receptor species. However, the environment typically contains many ligands, which can bind to the same receptors with different affinities, resulting in cross-talk. In traditional rate models, such cross-talk prevents accurate inference of concentrations of individual ligands. In contrast, here we show that knowing the precise timing sequence of stochastic binding and unbinding events allows one receptor to provide information about multiple ligands simultaneously and with a high accuracy. We show that such high-accuracy estimation of multiple concentrations can be realized with simple structural modifications of the familiar kinetic proofreading biochemical network diagram. We give two specific examples of such modifications. We argue that structural and functional features of real cellular biochemical sensory networks in immune cells, such as feedforward and feedback loops or ligand antagonism, sometimes can be understood as solutions to the accurate multi-ligand estimation problem.

Cells live in chemically complex environments with many different chemical ligands around them. Can cells estimate concentrations of more ligands than they have receptor types? In this paper, we show that, surprisingly, the answer is “yes”, and the estimation can be implemented with simple biochemical components already present in many cells. Therefore, cells may “know” a lot more about their environment and thus may be able to implement more complex and accurate response strategies than was previously thought.

Design of an embedded inverse-feedforward biomolecular tracking controller for enzymatic reaction processes

Feedback control is widely used in chemical engineering to improve the performance and robustness of chemical processes. Feedback controllers require a 'subtractor' that is able to compute the error between the process output and the reference signal. In the case of embedded biomolecular control circuits, subtractors designed using standard chemical reaction network theory can only realise one-sided subtraction, rendering standard controller design approaches inadequate. Here, we show how a biomolecular controller that allows tracking of required changes in the outputs of enzymatic reaction processes can be designed and implemented within the framework of chemical reaction network theory. The controller architecture employs an inversion-based feedforward controller that compensates for the limitations of the one-sided subtractor that generates the error signals for a feedback controller. The proposed approach requires significantly fewer chemical reactions to implement than alternative designs, and should have wide applicability throughout the fields of synthetic biology and biological engineering.

Demystifying the cytokine network: Mathematical models point the way

Cytokines provide the means by which immune cells communicate with each other and with parenchymal cells. There are over one hundred cytokines and many exist in families that share receptor components and signal transduction pathways, creating complex networks. Reductionist approaches to understanding the role of specific cytokines, through the use of gene-targeted mice, have revealed further complexity in the form of redundancy and pleiotropy in cytokine function. Creating an understanding of the complex interactions between cytokines and their target cells is challenging experimentally. Mathematical and computational modeling provides a robust set of tools by which complex interactions between cytokines can be studied and analyzed, in the process creating novel insights that can be further tested experimentally. This review will discuss and provide examples of the different modeling approaches that have been used to increase our understanding of cytokine networks. This includes discussion of knowledge-based and data-driven modeling approaches and the recent advance in single-cell analysis. The use of modeling to optimize cytokine-based therapies will also be discussed.

On the rejection-based algorithm for simulation and analysis of large-scale reaction networks

Stochastic simulation for in silico studies of large biochemical networks requires a great amount of computational time. We recently proposed a new exact simulation algorithm, called the rejection-based stochastic simulation algorithm (RSSA) [Thanh et al., J. Chem. Phys. 141(13), 134116 (2014)], to improve simulation performance by postponing and collapsing as much as possible the propensity updates. In this paper, we analyze the performance of this algorithm in detail, and improve it for simulating large-scale biochemical reaction networks. We also present a new algorithm, called simultaneous RSSA (SRSSA), which generates many independent trajectories simultaneously for the analysis of the biochemical behavior. SRSSA improves simulation performance by utilizing a single data structure across simulations to select reaction firings and forming trajectories. The memory requirement for building and storing the data structure is thus independent of the number of trajectories. The updating of the data structure when needed is performed collectively in a single operation across the simulations. The trajectories generated by SRSSA are exact and independent of each other by exploiting the rejection-based mechanism. We test our new improvement on real biological systems with a wide range of reaction networks to demonstrate its applicability and efficiency.

Symbolic computation of differential equivalences

Ordinary differential equations (ODEs) are widespread in many natural sciences including chemistry, ecology, and systems biology, and in disciplines such as control theory and electrical engineering. Building on the celebrated molecules-as-processes paradigm, they have become increasingly popular in computer science, with high-level languages and formal methods such as Petri nets, process algebra, and rule-based systems that are interpreted as ODEs. We consider the problem of comparing and minimizing ODEs automatically. Influenced by traditional approaches in the theory of programming, we propose differential equivalence relations. We study them for a basic intermediate language, for which we have decidability results, that can be targeted by a class of high-level specifications. An ODE implicitly represents an uncountable state space, hence reasoning techniques cannot be borrowed from established domains such as probabilistic programs with finite-state Markov chain semantics. We provide novel symbolic procedures to check an equivalence and compute the largest one via partition refinement algorithms that use satisfiability modulo theories. We illustrate the generality of our framework by showing that differential equivalences include (i) well-known notions for the minimization of continuous-time Markov chains (lumpability), (ii)~bisimulations for chemical reaction networks recently proposed by Cardelli et al., and (iii) behavioral relations for process algebra with ODE semantics. With a prototype implementation we are able to detect equivalences in biochemical models from the literature that cannot be reduced using competing automatic techniques.

Nanodomains in early and later phases of FcɛRI signalling

Our long-term efforts to elucidate receptor-mediated signalling in immune cells, particularly transmembrane signalling initiated by FcɛRI, the receptor for IgE in mast cells, led us unavoidably to contemplate the role of the heterogeneous plasma membrane. Our early investigations with fluorescence microscopy revealed co-redistribution of certain lipids and signalling components with antigen-cross-linked IgE-FcɛRI and pointed to participation of ordered membrane domains in the signalling process. With a focus on this function, we have worked along with others to develop diverse and increasingly sophisticated tools to analyse the complexity of membrane structure that facilitates regulation and targeting of signalling events. The present chapter describes how initial membrane interactions of clustered IgE-FcɛRI lead to downstream cellular responses and how biochemical information integrated with nanoscale resolution spectroscopy and imaging is providing mechanistic insights at the level of molecular complexes.

Modeling the T cell immune response: a fascinating challenge

The immune system is designed to protect the organism from infection and to repair damaged tissue. An effective response requires recognition of the threat, the appropriate effector mechanism to clear the pathogen and a return to homeostasis with minimal damage to self-tissues. T cells play a central role in orchestrating the immune response at all stages of the response and have been the subject of intense study by both experimental immunologists and modelers. This review examines some of the more critical questions in T cell biology and describes the latest attempts to address those questions using approaches that combine mathematical modeling and experiments.

Phosphorylation site dynamics of early T-cell receptor signaling

In adaptive immune responses, T-cell receptor (TCR) signaling impacts multiple cellular processes and results in T-cell differentiation, proliferation, and cytokine production. Although individual protein-protein interactions and phosphorylation events have been studied extensively, we lack a systems-level understanding of how these components cooperate to control signaling dynamics, especially during the crucial first seconds of stimulation. Here, we used quantitative proteomics to characterize reshaping of the T-cell phosphoproteome in response to TCR/CD28 co-stimulation, and found that diverse dynamic patterns emerge within seconds. We detected phosphorylation dynamics as early as 5 s and observed widespread regulation of key TCR signaling proteins by 30 s. Development of a computational model pointed to the presence of novel regulatory mechanisms controlling phosphorylation of sites with central roles in TCR signaling. The model was used to generate predictions suggesting unexpected roles for the phosphatase PTPN6 (SHP-1) and shortcut recruitment of the actin regulator WAS. Predictions were validated experimentally. This integration of proteomics and modeling illustrates a novel, generalizable framework for solidifying quantitative understanding of a signaling network and for elucidating missing links.

An Interaction Library for the FcεRI Signaling Network

Antigen receptors play a central role in adaptive immune responses. Although the molecular networks associated with these receptors have been extensively studied, we currently lack a systems-level understanding of how combinations of non-covalent interactions and post-translational modifications are regulated during signaling to impact cellular decision-making. To fill this knowledge gap, it will be necessary to formalize and piece together information about individual molecular mechanisms to form large-scale computational models of signaling networks. To this end, we have developed an interaction library for signaling by the high-affinity IgE receptor, FcεRI. The library consists of executable rules for protein–protein and protein–lipid interactions. This library extends earlier models for FcεRI signaling and introduces new interactions that have not previously been considered in a model. Thus, this interaction library is a toolkit with which existing models can be expanded and from which new models can be built. As an example, we present models of branching pathways from the adaptor protein Lat, which influence production of the phospholipid PIP₃ at the plasma membrane and the soluble second messenger IP₃. We find that inclusion of a positive feedback loop gives rise to a bistable switch, which may ensure robust responses to stimulation above a threshold level. In addition, the library is visualized to facilitate understanding of network circuitry and identification of network motifs.

Exact hybrid particle/population simulation of rule-based models of biochemical systems

Detailed modeling and simulation of biochemical systems is complicated by the problem of combinatorial complexity, an explosion in the number of species and reactions due to myriad protein-protein interactions and post-translational modifications. Rule-based modeling overcomes this problem by representing molecules as structured objects and encoding their interactions as pattern-based rules. This greatly simplifies the process of model specification, avoiding the tedious and error prone task of manually enumerating all species and reactions that can potentially exist in a system. From a simulation perspective, rule-based models can be expanded algorithmically into fully-enumerated reaction networks and simulated using a variety of network-based simulation methods, such as ordinary differential equations or Gillespie's algorithm, provided that the network is not exceedingly large. Alternatively, rule-based models can be simulated directly using particle-based kinetic Monte Carlo methods. This “network-free” approach produces exact stochastic trajectories with a computational cost that is independent of network size. However, memory and run time costs increase with the number of particles, limiting the size of system that can be feasibly simulated. Here, we present a hybrid particle/population simulation method that combines the best attributes of both the network-based and network-free approaches. The method takes as input a rule-based model and a user-specified subset of species to treat as population variables rather than as particles. The model is then transformed by a process of “partial network expansion” into a dynamically equivalent form that can be simulated using a population-adapted network-free simulator. The transformation method has been implemented within the open-source rule-based modeling platform BioNetGen, and resulting hybrid models can be simulated using the particle-based simulator NFsim. Performance tests show that significant memory savings can be achieved using the new approach and a monetary cost analysis provides a practical measure of its utility.

Rule-based modeling is a modeling paradigm that addresses the problem of combinatorial complexity in biochemical systems. The key idea is to specify only those components of a biological macromolecule that are directly involved in a biochemical transformation. Until recently, this “pattern-based” approach greatly simplified the process of model building but did nothing to improve the performance of model simulation. This changed with the introduction of “network-free” simulation methods, which operate directly on the compressed rule set of a rule-based model rather than on a fully-enumerated set of reactions and species. However, these methods represent every molecule in a system as a particle, limiting their use to systems containing less than a few million molecules. Here, we describe an extension to the network-free approach that treats rare, complex species as particles and plentiful, simple species as population variables, while retaining the exact dynamics of the model system. By making more efficient use of computational resources for species that do not require the level of detail of a particle representation, this hybrid particle/population approach can simulate systems much larger than is possible using network-free methods and is an important step towards realizing the practical simulation of detailed, mechanistic models of whole cells.

The duration of T cell stimulation is a critical determinant of cell fate and plasticity

Variations in T cell receptor (TCR) signal strength, as indicated by differential activation of downstream signaling pathways, determine the fate of naïve T cells after encounter with antigen. Low-strength signals favor differentiation into regulatory T (T(reg)) cells containing the transcription factor Foxp3, whereas high-strength signals favor generation of interleukin-2-producing T helper (T(H)) cells. We constructed a logic circuit model of TCR signaling pathways, a major feature of which is an incoherent feed-forward loop involving both TCR-dependent activation of Foxp3 and its inhibition by mammalian target of rapamycin (mTOR), leading to the transient appearance of Foxp3(+) cells under T(H) cell-generating conditions. Experiments confirmed this behavior and the prediction that the immunosuppressive cytokine TGF-β (transforming growth factor-β) could generate T(reg) cells even during continued Akt-mTOR signaling. We predicted that sustained mTOR activity could suppress FOXP3 expression upon TGF-β removal, suggesting a possible mechanism for the experimentally observed instability of Foxp3(+) cells. Our model predicted, and experiments confirmed, that transient stimulation of cells with high-dose antigen generated T(H), T(reg), and nonactivated cells in proportions depending on the duration of TCR stimulation. Experimental analysis of cells after antigen removal identified three populations that correlated with these T cell fates. Further analysis of simulations implicated a negative feedback loop involving Foxp3, the phosphatase PTEN, and Akt-mTOR in determining fate. These results suggest that there is a critical time after TCR stimulation during which heterogeneity in the differentiating population of cells leads to increased plasticity of cell fate.

Efficient stochastic simulation of chemical kinetics networks using a weighted ensemble of trajectories

We apply the "weighted ensemble" (WE) simulation strategy, previously employed in the context of molecular dynamics simulations, to a series of systems-biology models that range in complexity from a one-dimensional system to a system with 354 species and 3680 reactions. WE is relatively easy to implement, does not require extensive hand-tuning of parameters, does not depend on the details of the simulation algorithm, and can facilitate the simulation of extremely rare events. For the coupled stochastic reaction systems we study, WE is able to produce accurate and efficient approximations of the joint probability distribution for all chemical species for all time t. WE is also able to efficiently extract mean first passage times for the systems, via the construction of a steady-state condition with feedback. In all cases studied here, WE results agree with independent "brute-force" calculations, but significantly enhance the precision with which rare or slow processes can be characterized. Speedups over "brute-force" in sampling rare events via the Gillespie direct Stochastic Simulation Algorithm range from ~10(12) to ~10(18) for characterizing rare states in a distribution, and ~10(2) to ~10(4) for finding mean first passage times.

Stochastic transitions in a bistable reaction system on the membrane

Transitions between steady states of a multi-stable stochastic system in the perfectly mixed chemical reactor are possible only because of stochastic switching. In realistic cellular conditions, where diffusion is limited, transitions between steady states can also follow from the propagation of travelling waves. Here, we study the interplay between the two modes of transition for a prototype bistable system of kinase-phosphatase interactions on the plasma membrane. Within microscopic kinetic Monte Carlo simulations on the hexagonal lattice, we observed that for finite diffusion the behaviour of the spatially extended system differs qualitatively from the behaviour of the same system in the well-mixed regime. Even when a small isolated subcompartment remains mostly inactive, the chemical travelling wave may propagate, leading to the activation of a larger compartment. The activating wave can be induced after a small subdomain is activated as a result of a stochastic fluctuation. Such a spontaneous onset of activity is radically more probable in subdomains characterized by slower diffusion. Our results show that a local immobilization of substrates can lead to the global activation of membrane proteins by the mechanism that involves stochastic fluctuations followed by the propagation of a semi-deterministic travelling wave.

Systems model of T cell receptor proximal signaling reveals emergent ultrasensitivity

Receptor phosphorylation is thought to be tightly regulated because phosphorylated receptors initiate signaling cascades leading to cellular activation. The T cell antigen receptor (TCR) on the surface of T cells is phosphorylated by the kinase Lck and dephosphorylated by the phosphatase CD45 on multiple immunoreceptor tyrosine-based activation motifs (ITAMs). Intriguingly, Lck sequentially phosphorylates ITAMs and ZAP-70, a cytosolic kinase, binds to phosphorylated ITAMs with differential affinities. The purpose of multiple ITAMs, their sequential phosphorylation, and the differential ZAP-70 affinities are unknown. Here, we use a systems model to show that this signaling architecture produces emergent ultrasensitivity resulting in switch-like responses at the scale of individual TCRs. Importantly, this switch-like response is an emergent property, so that removal of multiple ITAMs, sequential phosphorylation, or differential affinities abolishes the switch. We propose that highly regulated TCR phosphorylation is achieved by an emergent switch-like response and use the systems model to design novel chimeric antigen receptors for therapy.

Recognition of antigen by the T cell antigen receptor (TCR) is a central event in the initiation of adaptive immune responses and for this reason the TCR has been extensively studied. Multiple studies performed over the past 20 years have revealed intriguing findings that include the observation that the TCR has multiple phosphorylation sites that are sequentially phosphorylated by the kinase Lck and that ZAP-70, a cytosolic kinase, binds to these sites with different affinities. The purpose of multiple sites, their sequential phosphorylation by Lck, and the differential binding affinities of ZAP-70 are unknown. Using a novel mechanistic model that incorporates a high level of molecular detail, we find, unexpectedly, that all factors are critical for producing ultrasensitivity (switch-like response) and therefore this signaling architecture exhibits systems-level emergent ultrasensitivity. We use the model to study existing therapeutic chimeric antigen receptors and in the design of novel ones. The work also has direct implications to the study of many other immune receptors.

Single-cell measurements of IgE-mediated FcεRI signaling using an integrated microfluidic platform

Heterogeneity in responses of cells to a stimulus, such as a pathogen or allergen, can potentially play an important role in deciding the fate of the responding cell population and the overall systemic response. Measuring heterogeneous responses requires tools capable of interrogating individual cells. Cell signaling studies commonly do not have single-cell resolution because of the limitations of techniques used such as Westerns, ELISAs, mass spectrometry, and DNA microarrays. Microfluidics devices are increasingly being used to overcome these limitations. Here, we report on a microfluidic platform for cell signaling analysis that combines two orthogonal single-cell measurement technologies: on-chip flow cytometry and optical imaging. The device seamlessly integrates cell culture, stimulation, and preparation with downstream measurements permitting hands-free, automated analysis to minimize experimental variability. The platform was used to interrogate IgE receptor (FcεRI) signaling, which is responsible for triggering allergic reactions, in RBL-2H3 cells. Following on-chip crosslinking of IgE-FcεRI complexes by multivalent antigen, we monitored signaling events including protein phosphorylation, calcium mobilization and the release of inflammatory mediators. The results demonstrate the ability of our platform to produce quantitative measurements on a cell-by-cell basis from just a few hundred cells. Model-based analysis of the Syk phosphorylation data suggests that heterogeneity in Syk phosphorylation can be attributed to protein copy number variations, with the level of Syk phosphorylation being particularly sensitive to the copy number of Lyn.

Leveraging modeling approaches: reaction networks and rules

We have witnessed an explosive growth in research involving mathematical models and computer simulations of intracellular molecular interactions, ranging from metabolic pathways to signaling and gene regulatory networks. Many software tools have been developed to aid in the study of such biological systems, some of which have a wealth of features for model building and visualization, and powerful capabilities for simulation and data analysis. Novel high-resolution and/or high-throughput experimental techniques have led to an abundance of qualitative and quantitative data related to the spatiotemporal distribution of molecules and complexes, their interactions kinetics, and functional modifications. Based on this information, computational biology researchers are attempting to build larger and more detailed models. However, this has proved to be a major challenge. Traditionally, modeling tools require the explicit specification of all molecular species and interactions in a model, which can quickly become a major limitation in the case of complex networks - the number of ways biomolecules can combine to form multimolecular complexes can be combinatorially large. Recently, a new breed of software tools has been created to address the problems faced when building models marked by combinatorial complexity. These have a different approach for model specification, using reaction rules and species patterns. Here we compare the traditional modeling approach with the new rule-based methods. We make a case for combining the capabilities of conventional simulation software with the unique features and flexibility of a rule-based approach in a single software platform for building models of molecular interaction networks.

Monovalent and multivalent ligation of the B cell receptor exhibit differential dependence upon Syk and Src family kinases

The Src and Syk families of kinases are two distinct sets of kinases that play critical roles in initiating membrane-proximal B cell receptor (BCR) signaling. However, unlike in other lymphocytes, such as T cells, the "division of labor" between Src family kinases (SFKs) and Syk in B cells is not well separated because both Syk and SFKs can phosphorylate immunoreceptor tyrosine-based activation motifs (ITAMs) present in proteins comprising the BCR. To understand why B cells require both SFKs and Syk for activation, we investigated the roles of both families of kinases in BCR signaling with computational modeling and in vitro experiments. Our computational model suggested that positive feedback enabled Syk to substantially compensate for the absence of SFKs when spatial clustering of BCRs was induced by multimeric ligands. We confirmed this prediction experimentally. In contrast, when B cells were stimulated by monomeric ligands that failed to produce BCR clustering, both Syk and SFKs were required for complete and rapid BCR activation. Our data suggest that SFKs could play a pivotal role in increasing BCR sensitivity to monomeric antigens of pathogens and in mediating a rapid response to soluble multimeric antigens of pathogens that can induce spatial BCR clustering.

A mechanistic model of early FcεRI signaling: lipid rafts and the question of protection from dephosphorylation

We present a model of the early events in mast cell signaling mediated by FcεRI where the plasma membrane is composed of many small ordered lipid domains (rafts), surrounded by a non-order region of lipids consisting of the remaining plasma membrane. The model treats the rafts as transient structures that constantly form and breakup, but that maintain a fixed average number per cell. The rafts have a high propensity for harboring Lyn kinase, aggregated, but not unaggregated receptors, and the linker for the activation of T cells (LAT). Phosphatase activity in the rafts is substantially reduced compared to the nonraft region. We use the model to analyze published experiments on the rat basophilic leukemia (RBL)-2H3 cell line that seem to contradict the notion that rafts offer protection. In these experiments IgE was cross-linked with a multivalent antigen and then excess monovalent hapten was added to break-up cross-links. The dephosphorylation of the unaggregated receptor (nonraft associated) and of LAT (raft associated) were then monitored in time and found to decay at similar rates, leading to the conclusion that rafts offer no protection from dephosphorylation. In the model, because the rafts are transient, a protein that is protected while in a raft will be subject to dephosphorylation when the raft breaks up and the protein finds itself in the nonraft region of the membrane. We show that the model is consistent with the receptor and LAT dephosphorylation experiments while still allowing rafts to enhance signaling by providing substantial protection from phosphatases.