Rules for modeling signal-transduction systems

Development and Analysis of Multiscale Models for Tuberculosis: From Molecules to Populations

Although infectious disease dynamics are often analyzed at the macro-scale, increasing numbers of drug-resistant infections highlight the importance of within-host modeling that simultaneously solves across multiple scales to effectively respond to epidemics. We review multiscale modeling approaches for complex, interconnected biological systems and discuss critical steps involved in building, analyzing, and applying such models within the discipline of model credibility. We also present our two tools: CaliPro, for calibrating multiscale models (MSMs) to datasets, and tunable resolution, for fine- and coarse-graining sub-models while retaining insights. We include as an example our work simulating infection with Mycobacterium tuberculosis to demonstrate modeling choices and how predictions are made to generate new insights and test interventions. We discuss some of the current challenges of incorporating novel datasets, rigorously training computational biologists, and increasing the reach of MSMs. We also offer several promising future research directions of incorporating within-host dynamics into applications ranging from combinatorial treatment to epidemic response.

Combining Multikinase Tyrosine Kinase Inhibitors Targeting the Vascular Endothelial Growth Factor and Cluster of Differentiation 47 Signaling Pathways Is Predicted to Increase the Efficacy of Antiangiogenic Combination Therapies

Angiogenesis is a critical step in tumor growth, development, and invasion. Nascent tumor cells secrete vascular endothelial growth factor (VEGF) that significantly remodels the tumor microenvironment through interaction with multiple receptors on vascular endothelial cells, including type 2 VEGF receptor (VEGFR2). The complex pathways initiated by VEGF binding to VEGFR2 lead to enhanced proliferation, survival, and motility of vascular endothelial cells and formation of a new vascular network, enabling tumor growth. Antiangiogenic therapies that inhibit VEGF signaling pathways were among the first drugs that targeted stroma rather than tumor cells. Despite improvements in progression-free survival and higher response rates relative to chemotherapy in some types of solid tumors, the impact on overall survival (OS) has been limited, with the majority of tumors eventually relapsing due to resistance or activation of alternate angiogenic pathways. Here, we developed a molecularly detailed computational model of endothelial cell signaling and angiogenesis-driven tumor growth to investigate combination therapies targeting different nodes of the endothelial VEGF/VEGFR2 signaling pathway. Simulations predicted a strong threshold-like behavior in extracellular signal-regulated kinases 1/2 (ERK1/2) activation relative to phosphorylated VEGFR2 levels, as continuous inhibition of at least 95% of receptors was necessary to abrogate phosphorylated ERK1/2 (pERK1/2). Combinations with mitogen-activated protein kinase/ERK kinase (MEK) and spingosine-1-phosphate inhibitors were found to be effective in overcoming the ERK1/2 activation threshold and abolishing activation of the pathway. Modeling results also identified a mechanism of resistance whereby tumor cells could reduce pERK1/2 sensitivity to inhibitors of VEGFR2 by upregulation of Raf, MEK, and sphingosine kinase 1 (SphK1), thus highlighting the need for deeper investigation of the dynamics of the crosstalk between VEGFR2 and SphK1 pathways. Inhibition of VEGFR2 phosphorylation was found to be more effective at blocking protein kinase B, also known as AKT, activation; however, to effectively abolish AKT activation, simulations identified Axl autophosphorylation or the Src kinase domain as potent targets. Simulations also supported activating cluster of differentiation 47 (CD47) on endothelial cells as an effective combination partner with tyrosine kinase inhibitors to inhibit angiogenesis signaling and tumor growth. Virtual patient simulations supported the effectiveness of CD47 agonism in combination with inhibitors of VEGFR2 and SphK1 pathways. Overall, the rule-based system model developed here provides new insights, generates novel hypothesis, and makes predictions regarding combinations that may enhance the OS with currently approved antiangiogenic therapies.

Mechanistic model of MAPK signaling reveals how allostery and rewiring contribute to drug resistance

BRAF is prototypical of oncogenes that can be targeted therapeutically and the treatment of BRAF^V600E melanomas with RAF and MEK inhibitors results in rapid tumor regression. However, drug-induced rewiring generates a drug adapted state thought to be involved in acquired resistance and disease recurrence. In this article, we study mechanisms of adaptive rewiring in BRAF^V600E melanoma cells using an energy-based implementation of ordinary differential equation (ODE) modeling in combination with proteomic, transcriptomic and imaging data. We develop a method for causal tracing of ODE models and identify two parallel MAPK reaction channels that are differentially sensitive to RAF and MEK inhibitors due to differences in protein oligomerization and drug binding. We describe how these channels, and timescale separation between immediate-early signaling and transcriptional feedback, create a state in which the RAS-regulated MAPK channel can be activated by growth factors under conditions in which the BRAF^V600E -driven channel is fully inhibited. Further development of the approaches in this article is expected to yield a unified model of adaptive drug resistance in melanoma.

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.

Computational modeling of RNase, antisense ORF0 RNA, and intracellular compartmentation and their impact on the life cycle of the line retrotransposon

Nearly half of the human genome is occupied by repetitive sequences of ancient virus-like genetic elements. The largest class, comprising 17% of the genome, belong to the type 1 Long INterspersed Elements (LINE-1) and are the only class capable of autonomous propagation in the genome. When epigenetic silencing mechanisms of LINE-1 fail, the proteins encoded by LINE-1 engage in reverse transcription to make new copies of their own or other DNAs that are pasted back into the genome. To elucidate how LINE-1 is dysregulated as a result of carcinogen exposure, we developed a computational model of key elements in the LINE-1 lifecycle, namely, the role of cytosolic ribonuclease (RNase), RNA interference (RNAi) by the antisense ORF0 RNA, and sequestration of LINE-1 products into stress granules and multivesicular structures.

The model showed that when carcinogen exposure is represented as either a sudden increase in LINE-1 mRNA count, or as an increase in mRNA transcription rate, the retrotransposon copy number exhibits a distinct threshold behavior above which LINE-1 enters a positive feedback loop that allows the cDNA copy number to grow exponentially. We also found that most of the LINE-1 RNA was degraded via the RNAase pathway and that neither ORF0 RNAi, nor the sequestration of LINE-1 products into granules and multivesicular structures, played a significant role in regulating the retrotransposon’s life cycle. Several aspects of the prediction agree with experimental results and indicate that the model has significant potential to inform future experiments related to LINE-1 activation.

Model comparison via simplicial complexes and persistent homology

In many scientific and technological contexts, we have only a poor understanding of the structure and details of appropriate mathematical models. We often, therefore, need to compare different models. With available data we can use formal statistical model selection to compare and contrast the ability of different mathematical models to describe such data. There is, however, a lack of rigorous methods to compare different models a priori. Here, we develop and illustrate two such approaches that allow us to compare model structures in a systematic way by representing models as simplicial complexes. Using well-developed concepts from simplicial algebraic topology, we define a distance between models based on their simplicial representations. Employing persistent homology with a flat filtration provides for alternative representations of the models as persistence intervals, which represent model structure, from which the model distances are also obtained. We then expand on this measure of model distance to study the concept of model equivalence to determine the conceptual similarity of models. We apply our methodology for model comparison to demonstrate an equivalence between a positional-information model and a Turing-pattern model from developmental biology, constituting a novel observation for two classes of models that were previously regarded as unrelated.

Computational Phosphorylation Network Reconstruction: An Update on Methods and Resources

Most proteins undergo some form of modification after translation, and phosphorylation is one of the most relevant and ubiquitous post-translational modifications. The succession of protein phosphorylation and dephosphorylation catalyzed by protein kinase and phosphatase, respectively, constitutes a key mechanism of molecular information flow in cellular systems. The protein interactions of kinases, phosphatases, and their regulatory subunits and substrates are the main part of phosphorylation networks. To elucidate the landscape of phosphorylation events has been a central goal pursued by both experimental and computational approaches. Substrate specificity (e.g., sequence, structure) or the phosphoproteome has been utilized in an array of different statistical learning methods to infer phosphorylation networks. In this chapter, different computational phosphorylation network inference-related methods and resources are summarized and discussed.

Atlas: automatic modeling of regulation of bacterial gene expression and metabolism using rule-based languages

MOTIVATION

Cells are complex systems composed of hundreds of genes whose products interact to produce elaborated behaviors. To control such behaviors, cells rely on transcription factors to regulate gene expression, and gene regulatory networks (GRNs) are employed to describe and understand such behavior. However, GRNs are static models, and dynamic models are difficult to obtain due to their size, complexity, stochastic dynamics and interactions with other cell processes.

RESULTS

We developed Atlas, a Python software that converts genome graphs and gene regulatory, interaction and metabolic networks into dynamic models. The software employs these biological networks to write rule-based models for the PySB framework. The underlying method is a divide-and-conquer strategy to obtain sub-models and combine them later into an ensemble model. To exemplify the utility of Atlas, we used networks of varying size and complexity of Escherichia coli and evaluated in silico modifications, such as gene knockouts and the insertion of promoters and terminators. Moreover, the methodology could be applied to the dynamic modeling of natural and synthetic networks of any bacteria.

AVAILABILITY AND IMPLEMENTATION

Code, models and tutorials are available online (https://github.com/networkbiolab/atlas).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Inferring phenomenological models of first passage processes

Biochemical processes in cells are governed by complex networks of many chemical species interacting stochastically in diverse ways and on different time scales. Constructing microscopically accurate models of such networks is often infeasible. Instead, here we propose a systematic framework for building phenomenological models of such networks from experimental data, focusing on accurately approximating the time it takes to complete the process, the First Passage (FP) time. Our phenomenological models are mixtures of Gamma distributions, which have a natural biophysical interpretation. The complexity of the models is adapted automatically to account for the amount of available data and its temporal resolution. The framework can be used for predicting behavior of FP systems under varying external conditions. To demonstrate the utility of the approach, we build models for the distribution of inter-spike intervals of a morphologically complex neuron, a Purkinje cell, from experimental and simulated data. We demonstrate that the developed models can not only fit the data, but also make nontrivial predictions. We demonstrate that our coarse-grained models provide constraints on more mechanistically accurate models of the involved phenomena.

Exploring the mechanistic and temporal regulation of LRP6 endocytosis in canonical WNT signaling

Endocytosis plays a pivotal regulatory role in canonical WNT signaling. Internalization of the low-density lipoprotein receptor-related protein 6 (LRP6) receptor complex can either promote or attenuate canonical WNT signaling, depending on the employed internalization pathway. Detailed analysis of the mechanism of LRP6 internalization and its temporal regulation is crucial for understanding the different cellular responses to WNT stimulation under varying conditions and in various cell types. Here, we elucidate the mechanisms involved in the internalization of LRP6 and re-evaluate existing, partly contradicting, theories on the regulation of LRP6 receptor internalization. We utilize a computational approach that aims at finding a set of mechanisms that accounts for the temporal dynamics of LRP6 receptor internalization upon WNT stimulation. Starting with a simple simulation model, we successively extend and probe the model's behavior based on quantitative measurements. The final model confirms that LRP6 internalization is clathrin independent in vertebrates, is not restricted to microdomains, and that signalosome formation delays LRP6 internalization within the microdomains. These findings partly revise the current understanding of LRP6 internalization in vertebrates.

LogicNet: probabilistic continuous logics in reconstructing gene regulatory networks

Background

Gene Regulatory Networks (GRNs) have been previously studied by using Boolean/multi-state logics. While the gene expression values are usually scaled into the range [0, 1], these GRN inference methods apply a threshold to discretize the data, resulting in missing information. Most of studies apply fuzzy logics to infer the logical gene-gene interactions from continuous data. However, all these approaches require an a priori known network structure.

Results

Here, by introducing a new probabilistic logic for continuous data, we propose a novel logic-based approach (called the LogicNet) for the simultaneous reconstruction of the GRN structure and identification of the logics among the regulatory genes, from the continuous gene expression data. In contrast to the previous approaches, the LogicNet does not require an a priori known network structure to infer the logics. The proposed probabilistic logic is superior to the existing fuzzy logics and is more relevant to the biological contexts than the fuzzy logics. The performance of the LogicNet is superior to that of several Mutual Information-based and regression-based tools for reconstructing GRNs.

Conclusions

The LogicNet reconstructs GRNs and logic functions without requiring prior knowledge of the network structure. Moreover, in another application, the LogicNet can be applied for logic function detection from the known regulatory genes-target interactions. We also conclude that computational modeling of the logical interactions among the regulatory genes significantly improves the GRN reconstruction accuracy.

NERDSS: A Nonequilibrium Simulator for Multibody Self-Assembly at the Cellular Scale

Currently, a significant barrier to building predictive models of cellular self-assembly processes is that molecular models cannot capture minutes-long dynamics that couple distinct components with active processes, whereas reaction-diffusion models cannot capture structures of molecular assembly. Here, we introduce the nonequilibrium reaction-diffusion self-assembly simulator (NERDSS), which addresses this spatiotemporal resolution gap. NERDSS integrates efficient reaction-diffusion algorithms into generalized software that operates on user-defined molecules through diffusion, binding and orientation, unbinding, chemical transformations, and spatial localization. By connecting the fast processes of binding with the slow timescales of large-scale assembly, NERDSS integrates molecular resolution with reversible formation of ordered, multisubunit complexes. NERDSS encodes models using rule-based formatting languages to facilitate model portability, usability, and reproducibility. Applying NERDSS to steps in clathrin-mediated endocytosis, we design multicomponent systems that can form lattices in solution or on the membrane, and we predict how stochastic but localized dephosphorylation of membrane lipids can drive lattice disassembly. The NERDSS simulations reveal the spatial constraints on lattice growth and the role of membrane localization and cooperativity in nucleating assembly. By modeling viral lattice assembly and recapitulating oscillations in protein expression levels for a circadian clock model, we illustrate the adaptability of NERDSS. NERDSS simulates user-defined assembly models that were previously inaccessible to existing software tools, with broad applications to predicting self-assembly in vivo and designing high-yield assemblies in vitro.

ModelBricks-modules for reproducible modeling improving model annotation and provenance

Most computational models in biology are built and intended for "single-use"; the lack of appropriate annotation creates models where the assumptions are unknown, and model elements are not uniquely identified. Simply recreating a simulation result from a publication can be daunting; expanding models to new and more complex situations is a herculean task. As a result, new models are almost always created anew, repeating literature searches for kinetic parameters, initial conditions and modeling specifics. It is akin to building a brick house starting with a pile of clay. Here we discuss a concept for building annotated, reusable models, by starting with small well-annotated modules we call ModelBricks. Curated ModelBricks, accessible through an open database, could be used to construct new models that will inherit ModelBricks annotations and thus be easier to understand and reuse. Key features of ModelBricks include reliance on a commonly used standard language (SBML), rule-based specification describing species as a collection of uniquely identifiable molecules, association with model specific numerical parameters, and more common annotations. Physical bricks can vary substantively; likewise, to be useful the structure of ModelBricks must be highly flexible-it should encapsulate mechanisms from single reactions to multiple reactions in a complex process. Ultimately, a modeler would be able to construct large models by using multiple ModelBricks, preserving annotations and provenance of model elements, resulting in a highly annotated model. We envision the library of ModelBricks to rapidly grow from community contributions. Persistent citable references will incentivize model creators to contribute new ModelBricks.

Rule-Based Modeling Using Wildcards in the Smoldyn Simulator

Many biological molecules exist in multiple variants, such as proteins with different posttranslational modifications, DNAs with different sequences, and phospholipids with different chain lengths. Representing these variants as distinct species, as most biochemical simulators do, leads to the problem that the number of species, and chemical reactions that interconvert them, typically increase combinatorially with the number of ways that the molecules can vary. This can be alleviated by "rule-based modeling methods," in which software generates the chemical reaction network from relatively simple "rules." This chapter presents a new approach to rule-based modeling. It is based on wildcards that match to species names, much as wildcards can match to file names in computer operating systems. It is much simpler to use than the formal rule-based modeling approaches developed previously but can lead to unintended consequences if not used carefully. This chapter demonstrates rule-based modeling with wildcards through examples for signaling systems, protein complexation, polymerization, nucleic acid sequence copying and mutation, the "SMILES" chemical notation, and others. The method is implemented in Smoldyn, a spatial and stochastic biochemical simulator, for both generate-first and on-the-fly expansion, meaning whether the reaction network is generated before or during the simulation.

Using Python for Spatially Resolved Modeling with Simmune

Mechanistic models are an important tool to gain insights about the quantitative behavior of cell-biological signal transduction networks. Here we show how Simmune can be used in conjunction with IPython to create repeatable, self-contained analyses of signal transduction processes in spatially inhomogeneous environments.

Mapping Tyrosine Kinase Receptor Dimerization to Receptor Expression and Ligand Affinities

Tyrosine kinase receptor (RTK) ligation and dimerization is a key mechanism for translating external cell stimuli into internal signaling events. This process is critical to several key cell and physiological processes, such as in angiogenesis and embryogenesis, among others. While modulating RTK activation is a promising therapeutic target, RTK signaling axes have been shown to involve complicated interactions between ligands and receptors both within and across different protein families. In angiogenesis, for example, several signaling protein families, including vascular endothelial growth factors and platelet-derived growth factors, exhibit significant cross-family interactions that can influence pathway activation. Computational approaches can provide key insight to detangle these signaling pathways but have been limited by the sparse knowledge of these cross-family interactions. Here, we present a framework for studying known and potential non-canonical interactions. We constructed generalized models of RTK ligation and dimerization for systems of two, three and four receptor types and different degrees of cross-family ligation. Across each model, we developed parameter-space maps that fully determine relative pathway activation for any set of ligand-receptor binding constants, ligand concentrations and receptor concentrations. Therefore, our generalized models serve as a powerful reference tool for predicting not only known ligand: Receptor axes but also how unknown interactions could alter signaling dimerization patterns. Accordingly, it will drive the exploration of cross-family interactions and help guide therapeutic developments across processes like cancer and cardiovascular diseases, which depend on RTK-mediated signaling.

Efficiently Encoding Complex Biochemical Models with the Multistate Model Builder (MSMB)

Biologists seek to create increasingly complex molecular regulatory network models. Writing such a model is a creative effort that requires flexible analysis tools and better modeling languages than offered by many of today's biochemical model editors. Our Multistate Model Builder (MSMB) supports multistate models created using different modeling styles that suit the modeler rather than the software. MSMB defines a simple but powerful syntax to describe multistate species. Our syntax reduces the number of reactions needed to encode the model, thereby reducing the cognitive load involved with model creation. MSMB gives extensive feedback during all stages of model creation. Users can activate error notifications, and use these notifications as a guide toward a consistent, syntactically correct model. Any consistent model can be exported to SBML or COPASI formats. We show the effectiveness of MSMB's multistate syntax through realistic models of cell cycle regulation and mRNA transcription. MSMB is an open-source project implemented in Java and it uses the COPASI API. Complete information and the installation package can be found at http://copasi.org/Projects/ .

Robustness of Nutrient Signaling Is Maintained by Interconnectivity Between Signal Transduction Pathways

Systems biology approaches provide means to study the interplay between biological processes leading to the mechanistic understanding of the properties of complex biological systems. Here, we developed a vector format rule-based Boolean logic model of the yeast S. cerevisiae cAMP-PKA, Snf1, and the Snf3-Rgt2 pathway to better understand the role of crosstalk on network robustness and function. We identified that phosphatases are the common unknown components of the network and that crosstalk from the cAMP-PKA pathway to other pathways plays a critical role in nutrient sensing events. The model was simulated with known crosstalk combinations and subsequent analysis led to the identification of characteristics and impact of pathway interconnections. Our results revealed that the interconnections between the Snf1 and Snf3-Rgt2 pathway led to increased robustness in these signaling pathways. Overall, our approach contributes to the understanding of the function and importance of crosstalk in nutrient signaling.

A Step-by-Step Guide to Using BioNetFit

BioNetFit is a software tool designed for solving parameter identification problems that arise in the development of rule-based models. It solves these problems through curve fitting (i.e., nonlinear regression). BioNetFit is compatible with deterministic and stochastic simulators that accept BioNetGen language (BNGL)-formatted files as inputs, such as those available within the BioNetGen framework. BioNetFit can be used on a laptop or stand-alone multicore workstation as well as on many Linux clusters, such as those that use the Slurm Workload Manager to schedule jobs. BioNetFit implements a metaheuristic population-based global optimization procedure, an evolutionary algorithm (EA), to minimize a user-defined objective function, such as a residual sum of squares (RSS) function. BioNetFit also implements a bootstrapping procedure for determining confidence intervals for parameter estimates. Here, we provide step-by-step instructions for using BioNetFit to estimate the values of parameters of a BNGL-encoded model and to define bootstrap confidence intervals. The process entails the use of several plain-text files, which are processed by BioNetFit and BioNetGen. In general, these files include (1) one or more EXP files, which each contains (experimental) data to be used in parameter identification/bootstrapping; (2) a BNGL file containing a model section, which defines a (rule-based) model, and an actions section, which defines simulation protocols that generate GDAT and/or SCAN files with model predictions corresponding to the data in the EXP file(s); and (3) a CONF file that configures the fitting/bootstrapping job and that defines algorithmic parameter settings.

In Silico Evolution of Signaling Networks Using Rule-Based Models: Bistable Response Dynamics

One of the ultimate goals in biology is to understand the design principles of biological systems. Such principles, if they exist, can help us better understand complex, natural biological systems and guide the engineering of de novo ones. Toward deciphering design principles, in silico evolution of biological systems with proper abstraction is a promising approach. Here, we demonstrate the application of in silico evolution combined with rule-based modeling for exploring design principles of cellular signaling networks. This application is based on a computational platform, called BioJazz, which allows in silico evolution of signaling networks with unbounded complexity. We provide a detailed introduction to BioJazz architecture and implementation and describe how it can be used to evolve and/or design signaling networks with defined dynamics. For the latter, we evolve signaling networks with switch-like response dynamics and demonstrate how BioJazz can result in new biological insights into network structures that can endow bistable response dynamics. This example also demonstrated both the power of BioJazz in evolving and designing signaling networks and its limitations at the current stage of development.

Soft Set Theory for Decision Making in Computational Biology under Incomplete Information

The study of biological systems is complex and of great importance. There exist numerous approaches to signal transduction processes, including symbolic modeling of cellular adaptation. The use of formal methods for computational systems biology eases the analysis of cellular models and the establishment of the causes and consequences of certain cellular situations associated to diseases. In this paper, we define an application of logic modeling with rewriting logic and soft set theory. Our approach to decision making with soft sets offers a novel strategy that complements standard strategies. We implement a metalevel strategy to control and guide the rewriting process of the Maude rewriting engine. In particular, we adapt mathematical methods to capture imprecision, vagueness, and uncertainty in the available data. Using this new strategy, we propose an extension in the biological symbolic models of Pathway Logic. Our ultimate aim is to automatically determine the rules that are most appropriate and adjusted to reality in dynamic systems using decision making with incomplete soft sets.

Spatial Rule-Based Simulations: The SRSim Software

SRSim combines rule-based reaction network models with spatial particle simulations allowing to simulate the dynamics of large molecular complexes changing according to a set of chemical reaction rules. As the rule can contain patterns of molecular complexes and specific states of certain binding sites, a combinatorially complex or even infinitely sized reaction network can be defined. Particles move in a three-dimensional space according to molecular dynamics implemented by LAMMPS, while the BioNetGen language is used to formulate reaction rules. Geometric information is added in a specific XML format. The simulation protocol is exemplified by two different variants of polymerization as well as a toy model of DNA helix formation. SRSim is open source and available for download.

Intrinsic limits of information transmission in biochemical signalling motifs

All living things have evolved to sense changes in their environment in order to respond in adaptive ways. At the cellular level, these sensing systems generally involve receptor molecules at the cell surface, which detect changes outside the cell and relay those changes to the appropriate response elements downstream. With the advent of experimental technologies that can track signalling at the single-cell level, it has become clear that many signalling systems exhibit significant levels of 'noise,' manifesting as differential responses of otherwise identical cells to the same environment. This noise has a large impact on the capacity of cell signalling networks to transmit information from the environment. Application of information theory to experimental data has found that all systems studied to date encode less than 2.5 bits of information, with the majority transmitting significantly less than 1 bit. Given the growing interest in applying information theory to biological data, it is crucial to understand whether the low values observed to date represent some sort of intrinsic limit on information flow given the inherently stochastic nature of biochemical signalling events. In this work, we used a series of computational models to explore how much information a variety of common 'signalling motifs' can encode. We found that the majority of these motifs, which serve as the basic building blocks of cell signalling networks, can encode far more information (4-6 bits) than has ever been observed experimentally. In addition to providing a consistent framework for estimating information-theoretic quantities from experimental data, our findings suggest that the low levels of information flow observed so far in living system are not necessarily due to intrinsic limitations. Further experimental work will be needed to understand whether certain cell signalling systems actually can approach the intrinsic limits described here, and to understand the sources and purpose of the variation that reduces information flow in living cells.

Modeling the Self-Assembly of Protein Complexes through a Rigid-Body Rotational Reaction-Diffusion Algorithm

The reaction-diffusion equations provide a powerful framework for modeling nonequilibrium, cell-scale dynamics over the long time scales that are inaccessible by traditional molecular modeling approaches. Single-particle reaction-diffusion offers the highest resolution technique for tracking such dynamics, but it has not been applied to the study of protein self-assembly due to its treatment of reactive species as single-point particles. Here, we develop a relatively simple but accurate approach for building rigid structure and rotation into single-particle reaction-diffusion methods, providing a rate-based method for studying protein self-assembly. Our simplifying assumption is that reactive collisions can be evaluated purely on the basis of the separations between the sites, and not their orientations. The challenge of evaluating reaction probabilities can then be performed using well-known equations based on translational diffusion in both 3D and 2D, by employing an effective diffusion constant we derive here. We show how our approach reproduces both the kinetics of association, which is altered by rotational diffusion, and the equilibrium of reversible association, which is not. Importantly, the macroscopic kinetics of association can be predicted on the basis of the microscopic parameters of our structurally resolved model, allowing for critical comparisons with theory and other rate-based simulations. We demonstrate this method for efficient, rate-based simulations of self-assembly of clathrin trimers, highlighting how formation of regular lattices impacts the kinetics of association.

Computational modeling of synergistic interaction between αVβ3 integrin and VEGFR2 in endothelial cells: Implications for the mechanism of action of angiogenesis-modulating integrin-binding peptides

Cooperation between VEGFR2 and integrin αVβ3 is critical for neovascularization in wound healing, cardiovascular ischemic diseases, ocular diseases, and tumor angiogenesis. In the present study, we developed a rule-based computational model to investigate the potential mechanism by which the Src-induced integrin association with VEGFR2 enhances VEGFR2 activation. Simulations demonstrated that the main function of integrin is to reduce the degradation of VEGFR2 and hence stabilize the activation signal. In addition, receptor synthesis rate and recruitment from internal compartment were found to be sensitive determinants of the activation state of VEGFR2. The model was then applied to simulate the effect of integrin-binding peptides such as tumstatin and cilengitide on VEGFR2 signaling. Further, computational modeling proposed potential molecular mechanisms for the angiogenesis-modulating activity of other integrin-binding peptides. The model highlights the complexity of the crosstalk between αVβ3 integrin and VEGFR2 and the necessity of utilizing models to elucidate potential mechanisms in angiogenesis-modulating peptide therapy.

Computer Simulation of TSP1 Inhibition of VEGF-Akt-eNOS: An Angiogenesis Triple Threat

The matricellular protein thrombospondin-1 (TSP1) is a potent inhibitor of angiogenesis. Specifically, TSP1 has been experimentally shown to inhibit signaling downstream of vascular endothelial growth factor (VEGF). The molecular mechanism of this inhibition is not entirely clear. We developed a detailed computational model of VEGF signaling to Akt-endothelial nitric oxide synthase (eNOS) to investigate the quantitative molecular mechanism of TSP1 inhibition. The model demonstrated that TSP1 acceleration of VEGFR2 degradation is sufficient to explain the inhibition of VEGFR2 and eNOS phosphorylation. However, Akt inhibition requires TSP1-induced phosphatase recruitment to VEGFR2. The model was then utilized to test various strategies for the rescue of VEGF signaling to Akt and eNOS. Inhibiting TSP1 was predicted to be not as effective as CD47 depletion in rescuing signaling to Akt. The model further predicts that combination strategy involving depletion of CD47 and inhibition of TSP1 binding to CD47 is necessary for effective recovery of signaling to eNOS. In all, computational modeling offers insight to molecular mechanisms involving TSP1 interaction with VEGF signaling and provides strategies for rescuing angiogenesis by targeting TSP1-CD47 axis.

Biophysical attributes that affect CaMKII activation deduced with a novel spatial stochastic simulation approach

Calcium/calmodulin-dependent protein kinase II (CaMKII) holoenzymes play a critical role in decoding Ca2+ signals in neurons. Understanding how this occurs has been the focus of numerous studies including many that use models. However, CaMKII is notoriously difficult to simulate in detail because of its multi-subunit nature, which causes a combinatorial explosion in the number of species that must be modeled. To study the Ca2+-calmodulin-CaMKII reaction network with detailed kinetics while including the effect of diffusion, we have customized an existing stochastic particle-based simulator, Smoldyn, to manage the problem of combinatorial explosion. With this new method, spatial and temporal aspects of the signaling network can be studied without compromising biochemical details. We used this new method to examine how calmodulin molecules, both partially loaded and fully loaded with Ca2+, choose pathways to interact with and activate CaMKII under various Ca2+ input conditions. We found that the dependence of CaMKII phosphorylation on Ca2+ signal frequency is intrinsic to the network kinetics and the activation pattern can be modulated by the relative amount of Ca2+ to calmodulin and by the rate of Ca2+ diffusion. Depending on whether Ca2+ influx is saturating or not, calmodulin molecules could choose different routes within the network to activate CaMKII subunits, resulting in different frequency dependence patterns. In addition, the size of the holoenzyme produces a subtle effect on CaMKII activation. The more extended the subunits are organized, the easier for calmodulin molecules to access and activate the subunits. The findings suggest that particular intracellular environmental factors such as crowding and calmodulin availability can play an important role in decoding Ca2+ signals and can give rise to distinct CaMKII activation patterns in dendritic spines, Ca2+ channel nanodomains and cytoplasm.

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.

Integrative biological simulation praxis: Considerations from physics, philosophy, and data/model curation practices

Integrative biological simulations have a varied and controversial history in the biological sciences. From computational models of organelles, cells, and simple organisms, to physiological models of tissues, organ systems, and ecosystems, a diverse array of biological systems have been the target of large-scale computational modeling efforts. Nonetheless, these research agendas have yet to prove decisively their value among the broader community of theoretical and experimental biologists. In this commentary, we examine a range of philosophical and practical issues relevant to understanding the potential of integrative simulations. We discuss the role of theory and modeling in different areas of physics and suggest that certain sub-disciplines of physics provide useful cultural analogies for imagining the future role of simulations in biological research. We examine philosophical issues related to modeling which consistently arise in discussions about integrative simulations and suggest a pragmatic viewpoint that balances a belief in philosophy with the recognition of the relative infancy of our state of philosophical understanding. Finally, we discuss community workflow and publication practices to allow research to be readily discoverable and amenable to incorporation into simulations. We argue that there are aligned incentives in widespread adoption of practices which will both advance the needs of integrative simulation efforts as well as other contemporary trends in the biological sciences, ranging from open science and data sharing to improving reproducibility.

Information processing in the NF-κB pathway

The NF-κB pathway is known to transmit merely 1 bit of information about stimulus level. We combined experimentation with mathematical modeling to elucidate how information about TNF concentration is turned into a binary decision. Using Kolmogorov-Smirnov distance, we quantified the cell’s ability to discern 8 TNF concentrations at each step of the NF-κB pathway, to find that input discernibility decreases as signal propagates along the pathway. Discernibility of low TNF concentrations is restricted by noise at the TNF receptor level, whereas discernibility of high TNF concentrations it is restricted by saturation/depletion of downstream signaling components. Consequently, signal discernibility is highest between 0.03 and 1 ng/ml TNF. Simultaneous exposure to TNF or LPS and a translation inhibitor, cycloheximide, leads to prolonged NF-κB activation and a marked increase of transcript levels of NF-κB inhibitors, IκBα and A20. The impact of cycloheximide becomes apparent after the first peak of nuclear NF-κB translocation, meaning that the NF-κB network not only relays 1 bit of information to coordinate the all-or-nothing expression of early genes, but also over a longer time course integrates information about other stimuli. The NF-κB system should be thus perceived as a feedback-controlled decision-making module rather than a simple information transmission channel.

Mathematical models of microbial growth and metabolism: a whole-organism perspective

We review the principles underpinning the development of mathematical models of the metabolic activities of micro-organisms. Such models are important to understand and chart the substantial contributions made by micro-organisms to geochemical cycles, and also to optimise the performance of bioreactors that exploit the biochemical capabilities of these organisms. We advocate an approach based on the principle of dynamic allocation. We survey the biological background that motivates this approach, including nutrient assimilation, the regulation of gene expression, and the principles of microbial growth. In addition, we discuss the classic models of microbial growth as well as contemporary approaches. The dynamic allocation theory generalises these classic models in a natural manner and is readily amenable to the additional information provided by transcriptomics and proteomics approaches. Finally, we touch upon these organising principles in the context of the transition from the free-living unicellular mode of life to multicellularity.

Derivative-Free Optimization of Rate Parameters of Capsid Assembly Models from Bulk in Vitro Data

The assembly of virus capsids proceeds by a complicated cascade of association and dissociation steps, the great majority of which cannot be directly experimentally observed. This has made capsid assembly a rich field for computational models, but there are substantial obstacles to model inference for such systems. Here, we describe progress on fitting kinetic rate constants defining capsid assembly models to experimental data, a difficult data-fitting problem because of the high computational cost of simulating assembly trajectories, the stochastic noise inherent to the models, and the limited and noisy data available for fitting. We evaluate the merits of data-fitting methods based on derivative-free optimization (DFO) relative to gradient-based methods used in prior work. We further explore the advantages of alternative data sources through simulation of a model of time-resolved mass spectrometry data, a technology for monitoring bulk capsid assembly that can be expected to provide much richer data than previously used static light scattering approaches. The results show that advances in both the data and the algorithms can improve model inference. More informative data sources lead to high-quality fits for all methods, but DFO methods show substantial advantages on less informative data sources that better represent current experimental practice.

Multisite Phosphorylation Modulates the T Cell Receptor ζ-Chain Potency but not the Switchlike Response

Multisite phosphorylation is ubiquitous in cellular signaling and is thought to provide signaling proteins with additional regulatory mechanisms. Indeed, mathematical models have revealed a large number of mechanisms by which multisite phosphorylation can produce switchlike responses. The T cell antigen receptor (TCR) is a multisubunit receptor on the surface of T cells that is a prototypical multisite substrate as it contains 20 sites that are distributed on 10 conserved immunoreceptor tyrosine-based activation motifs (ITAMs). The TCR ζ-chain is a homodimer subunit that contains six ITAMs (12 sites) and exhibits a number of properties that are predicted to be sufficient for a switchlike response. We have used cellular reconstitution to systematically study multisite phosphorylation of the TCR ζ-chain. We find that multisite phosphorylation proceeds by a nonsequential random mechanism, and find no evidence that multiple ITAMs modulate a switchlike response but do find that they alter receptor potency and maximum phosphorylation. Modulation of receptor potency can be explained by a reduction in molecular entropy of the disordered ζ-chain upon phosphorylation. We further find that the tyrosine kinase ZAP-70 increases receptor potency but does not modulate the switchlike response. In contrast to other multisite proteins, where phosphorylations act in strong concert to modulate protein function, we suggest that the multiple ITAMs on the TCR function mainly to amplify subsequent signaling.

Multivalent Molecules as Modulators of RNA Granule Size and Composition

RNA granules are ensembles of specific RNA and protein molecules that mediate localized translation in eukaryotic cells. The mechanisms for formation and selectivity of RNA granules are unknown. Here we present a model for assembly of one type of RNA granule based on experimentally measured binding interactions among three core multivalent molecular components necessary for such assembly: specific RNA molecules that contain a cis-acting sequence called the A2 response element (A2RE), hnRNP A2 proteins that bind specifically (with high affinity) to A2RE sequences or nonspecifically (with lower affinity) to other RNA sequences, and heptavalent protein cytoskeleton-associated protein 5 (CKAP5, an alternative name for TOG protein) that binds both hnRNP A2 molecules and RNA. Non-A2RE RNA molecules (RNA without the A2RE sequence) that may be recruited to the granules through nonspecific interactions are also considered in the model. Modeling multivalent molecular interactions in granules is challenging because of combinatorial complexity in the number of potential molecular complexes among these core components and dynamic changes in granule composition and structure in response to changes in local intracellular environment. We use a hybrid modeling approach (deterministic-stochastic-statistical) that is appropriate when the overall compositions of multimolecular ensembles are of greater importance than the specific interactions among individual molecular components. Modeling studies titrating the concentrations of various granule components and varying effective site pair affinities and RNA valency demonstrate that interactions between multivalent components (TOG and RNA) are modulated by a bivalent adaptor molecule (hnRNP A2). Formation and disruption of granules, as well as RNA selectivity in granule composition are regulated by distinct concentration regimes of A2. Our results suggest that granule assembly is tightly controlled by multivalent molecular interactions among RNA molecules, adaptor proteins, and scaffold proteins.

Computational investigation of sphingosine kinase 1 (SphK1) and calcium dependent ERK1/2 activation downstream of VEGFR2 in endothelial cells

Vascular endothelial growth factor (VEGF) is a powerful regulator of neovascularization. VEGF binding to its cognate receptor, VEGFR2, activates a number of signaling pathways including ERK1/2. Activation of ERK1/2 is experimentally shown to involve sphingosine kinase 1 (SphK1) activation and its calcium-dependent translocation downstream of ERK1/2. Here we construct a rule-based computational model of signaling downstream of VEGFR2, by including SphK1 and calcium positive feedback mechanisms, and investigate their consequences on ERK1/2 activation. The model predicts the existence of VEGF threshold in ERK1/2 activation that can be continuously tuned by cellular concentrations of SphK1 and sphingosine 1 phosphate (S1P). The computer model also predicts powerful effects of perturbations in plasma and ER calcium pump rates and the current through the CRAC channels on ERK1/2 activation dynamics, highlighting the critical role of intracellular calcium in shaping the pERK1/2 signal. The model is then utilized to simulate anti-angiogenic therapeutic interventions targeting VEGFR2-ERK1/2 axis. Simulations indicate that monotherapies that exclusively target VEGFR2 phosphorylation, VEGF, or VEGFR2 are ineffective in shutting down signaling to ERK1/2. By simulating therapeutic strategies that target multiple nodes of the pathway such as Raf and SphK1, we conclude that combination therapy should be much more effective in blocking VEGF signaling to EKR1/2. The model has important implications for interventions that target signaling pathways in angiogenesis relevant to cancer, vascular diseases, and wound healing.

Vascular endothelial growth factor (VEGF) signaling is a potent regulator of angiogenesis, the growth and development of new vessels out of a preexisting vascular network. Angiogenesis requires enhanced survival, proliferation, and motility of the vascular endothelial cells. Crucial signaling endpoints in VEGF-mediated angiogenic response include elevation in intracellular calcium and the activation of the proteins ERK1 and 2 (ERK1/2). In this study, we have developed a novel computer model for the activation of ERK1/2 and calcium downstream of VEGF receptor type 2 (VEGFR2). Our model is the first of its kind to incorporate and investigate the consequences of calcium elevation and the role of a cellular lipid modifier known as sphingosine kinase 1 (SphK1). We also utilize the model to simulate therapeutic strategies targeting VEGF signaling to ERK1/2 indicating inefficiency of single therapies known as tyrosine kinase inhibitors (TKI) that target receptor phosphorylation. Computer simulations indicate that combination therapy is essential for effective blockade of this important pathway. Our results have important implications for human diseases such as cancer where plethora of anti-VEGF therapies are currently employed. Overall, our computer model sheds new light on a complex feedback involving SphK1 and calcium that radically alters the response of cells to VEGF.

Modeling of Receptor Tyrosine Kinase Signaling: Computational and Experimental Protocols

The advent of systems biology has convincingly demonstrated that the integration of experiments and dynamic modelling is a powerful approach to understand the cellular network biology. Here we present experimental and computational protocols that are necessary for applying this integrative approach to the quantitative studies of receptor tyrosine kinase (RTK) signaling networks. Signaling by RTKs controls multiple cellular processes, including the regulation of cell survival, motility, proliferation, differentiation, glucose metabolism, and apoptosis. We describe methods of model building and training on experimentally obtained quantitative datasets, as well as experimental methods of obtaining quantitative dose-response and temporal dependencies of protein phosphorylation and activities. The presented methods make possible (1) both the fine-grained modeling of complex signaling dynamics and identification of salient, course-grained network structures (such as feedback loops) that bring about intricate dynamics, and (2) experimental validation of dynamic models.

Importins promote high-frequency NF-κB oscillations increasing information channel capacity

Background

Importins and exportins influence gene expression by enabling nucleocytoplasmic shuttling of transcription factors. A key transcription factor of innate immunity, NF-κB, is sequestered in the cytoplasm by its inhibitor, IκBα, which masks nuclear localization sequence of NF-κB. In response to TNFα or LPS, IκBα is degraded, which allows importins to bind NF-κB and shepherd it across nuclear pores. NF-κB nuclear activity is terminated when newly synthesized IκBα enters the nucleus, binds NF-κB and exportin which directs the complex to the cytoplasm. Although importins/exportins are known to regulate spatiotemporal kinetics of NF-κB and other transcription factors governing innate immunity, the mechanistic details of these interactions have not been elucidated and mathematically modelled.

Results

Based on our quantitative experimental data, we pursue NF-κB system modelling by explicitly including NF-κB–importin and IκBα–exportin binding to show that the competition between importins and IκBα enables NF-κB nuclear translocation despite high levels of IκBα. These interactions reduce the effective relaxation time and allow the NF-κB regulatory pathway to respond to recurrent TNFα pulses of 45-min period, which is about twice shorter than the characteristic period of NF-κB oscillations. By stochastic simulations of model dynamics we demonstrate that randomly appearing, short TNFα pulses can be converted to essentially digital pulses of NF-κB activity, provided that intervals between input pulses are not shorter than 1 h.

Conclusions

By including interactions involving importin-α and exportin we bring the modelling of spatiotemporal kinetics of transcription factors to a more mechanistic level. Basing on the analysis of the pursued model we estimated the information transmission rate of the NF-κB pathway as 1 bit per hour.

Reviewers

This article was reviewed by Marek Kimmel, James Faeder and William Hlavacek.

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-016-0164-z) contains supplementary material.

Rule-based modeling with Virtual Cell

Rule-based modeling is invaluable when the number of possible species and reactions in a model become too large to allow convenient manual specification. The popular rule-based software tools BioNetGen and NFSim provide powerful modeling and simulation capabilities at the cost of learning a complex scripting language which is used to specify these models. Here, we introduce a modeling tool that combines new graphical rule-based model specification with existing simulation engines in a seamless way within the familiar Virtual Cell (VCell) modeling environment. A mathematical model can be built integrating explicit reaction networks with reaction rules. In addition to offering a large choice of ODE and stochastic solvers, a model can be simulated using a network free approach through the NFSim simulation engine.

AVAILABILITY AND IMPLEMENTATION

Available as VCell (versions 6.0 and later) at the Virtual Cell web site (http://vcell.org/). The application installs and runs on all major platforms and does not require registration for use on the user's computer. Tutorials are available at the Virtual Cell website and Help is provided within the software. Source code is available at Sourceforge.

CONTACT

vcell_support@uchc.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Non-canonical antagonism of PI3K by the kinase Itpkb delays thymocyte β-selection and renders it Notch-dependent

β-selection is the most pivotal event determining αβ T cell fate. Here, surface-expression of a pre-T cell receptor (pre-TCR) induces thymocyte metabolic activation, proliferation, survival and differentiation. Besides the pre-TCR, β-selection also requires co-stimulatory signals from Notch receptors - key cell fate determinants in eukaryotes. Here, we show that this Notch-dependence is established through antagonistic signaling by the pre-TCR/Notch effector, phosphoinositide 3-kinase (PI3K), and by inositol-trisphosphate 3-kinase B (Itpkb). Canonically, PI3K is counteracted by the lipid-phosphatases Pten and Inpp5d/SHIP-1. In contrast, Itpkb dampens pre-TCR induced PI3K/Akt signaling by producing IP₄, a soluble antagonist of the Akt-activating PI3K-product PIP₃. Itpkb^-/- thymocytes are pre-TCR hyperresponsive, hyperactivate Akt, downstream mTOR and metabolism, undergo an accelerated β-selection and can develop to CD4⁺CD8⁺ cells without Notch. This is reversed by inhibition of Akt, mTOR or glucose metabolism. Thus, non-canonical PI3K-antagonism by Itpkb restricts pre-TCR induced metabolic activation to enforce coincidence-detection of pre-TCR expression and Notch-engagement.

DOI:http://dx.doi.org/10.7554/eLife.10786.001

T cells defend our body against cancer and infectious agents such as viruses. However, they can also cause rheumatoid arthritis and other autoimmune diseases by attacking healthy tissue. T cells recognize target cells via receptor proteins on their surface. To maximize the variety of infections and cancers our immune system can recognize, we generate millions of T cells with different T cell receptors every day.

To ensure T cells work correctly, T cell receptors are tested at various checkpoints. The first checkpoint involves a process called beta (β) selection, during which T cells produce their first T cell receptor – the so-called pre-T cell receptor. This receptor causes T cells to divide and mature, and sets their future identity or “fate”. To complete β-selection, T cells must also receive signals from another surface receptor – one that belongs to the Notch family, which determines cell fate in many different tissues.

The Notch receptor and the pre-T cell receptor both activate an enzyme called PI3K – a key mediator of β-selection. But the pre-T cell receptor also activates another enzyme called Itpkb that is required for T cell development. Westernberg, Conche et al. have now investigated how these different proteins and signaling processes work and interact during β-selection, using mice that lack several immune genes, including the gene that produces Itpkb.

The results of the experiments show that during β-selection, Itpkb limits the ability of PI3K to activate some of its key target proteins. This “dampened” PI3K signaling ensures that both the pre-T cell receptor and the Notch receptor must be activated to trigger T cell maturation. Without Itpkb, β-selection can occur in the absence of Notch signaling.

As Notch signaling is important for determining the fate of many different cell types, Westernberg, Conche et al.’s findings raise the possibility that Itpkb might also regulate cell fate determination in other tissues. Moreover, Itpkb may suppress tumor development, because excessive PI3K signaling drives many cancers.

DOI:http://dx.doi.org/10.7554/eLife.10786.002

Computational Modeling, Formal Analysis, and Tools for Systems Biology

As the amount of biological data in the public domain grows, so does the range of modeling and analysis techniques employed in systems biology. In recent years, a number of theoretical computer science developments have enabled modeling methodology to keep pace. The growing interest in systems biology in executable models and their analysis has necessitated the borrowing of terms and methods from computer science, such as formal analysis, model checking, static analysis, and runtime verification. Here, we discuss the most important and exciting computational methods and tools currently available to systems biologists. We believe that a deeper understanding of the concepts and theory highlighted in this review will produce better software practice, improved investigation of complex biological processes, and even new ideas and better feedback into computer science.

Use of mechanistic models to integrate and analyze multiple proteomic datasets

Proteins in cell signaling networks tend to interact promiscuously through low-affinity interactions. Consequently, evaluating the physiological importance of mapped interactions can be difficult. Attempts to do so have tended to focus on single, measurable physicochemical factors, such as affinity or abundance. For example, interaction importance has been assessed on the basis of the relative affinities of binding partners for a protein of interest, such as a receptor. However, multiple factors can be expected to simultaneously influence the recruitment of proteins to a receptor (and the potential of these proteins to contribute to receptor signaling), including affinity, abundance, and competition, which is a network property. Here, we demonstrate that measurements of protein copy numbers and binding affinities can be integrated within the framework of a mechanistic, computational model that accounts for mass action and competition. We use cell line-specific models to rank the relative importance of protein-protein interactions in the epidermal growth factor receptor (EGFR) signaling network for 11 different cell lines. Each model accounts for experimentally characterized interactions of six autophosphorylation sites in EGFR with proteins containing a Src homology 2 and/or phosphotyrosine-binding domain. We measure importance as the predicted maximal extent of recruitment of a protein to EGFR following ligand-stimulated activation of EGFR signaling. We find that interactions ranked highly by this metric include experimentally detected interactions. Proteins with high importance rank in multiple cell lines include proteins with recognized, well-characterized roles in EGFR signaling, such as GRB2 and SHC1, as well as a protein with a less well-defined role, YES1. Our results reveal potential cell line-specific differences in recruitment.

Insights into the Role of Chemokines, Damage-Associated Molecular Patterns, and Lymphocyte-Derived Mediators from Computational Models of Trauma-Induced Inflammation

SIGNIFICANCE

Traumatic injury elicits a complex, dynamic, multidimensional inflammatory response that is intertwined with complications such as multiple organ dysfunction and nosocomial infection. The complex interplay between inflammation and physiology in critical illness remains a challenge for translational research, including the extrapolation to human disease from animal models.

RECENT ADVANCES

Over the past decade, we and others have attempted to decipher the biocomplexity of inflammation in these settings of acute illness, using computational models to improve clinical translation. In silico modeling has been suggested as a computationally based framework for integrating data derived from basic biology experiments as well as preclinical and clinical studies.

CRITICAL ISSUES

Extensive studies in cells, mice, and human blunt trauma patients have led us to suggest (i) that while an adequate level of inflammation is required for healing post-trauma, inflammation can be harmful when it becomes self-sustaining via a damage-associated molecular pattern/Toll-like receptor-driven feed-forward circuit; (ii) that chemokines play a central regulatory role in driving either self-resolving or self-maintaining inflammation that drives the early activation of both classical innate and more recently recognized lymphoid pathways; and (iii) the presence of multiple thresholds and feedback loops, which could significantly affect the propagation of inflammation across multiple body compartments.

FUTURE DIRECTIONS

These insights from data-driven models into the primary drivers and interconnected networks of inflammation have been used to generate mechanistic computational models. Together, these models may be used to gain basic insights as well as serving to help define novel biomarkers and therapeutic targets.

Automated adaptive inference of phenomenological dynamical models

Dynamics of complex systems is often driven by large and intricate networks of microscopic interactions, whose sheer size obfuscates understanding. With limited experimental data, many parameters of such dynamics are unknown, and thus detailed, mechanistic models risk overfitting and making faulty predictions. At the other extreme, simple ad hoc models often miss defining features of the underlying systems. Here we develop an approach that instead constructs phenomenological, coarse-grained models of network dynamics that automatically adapt their complexity to the available data. Such adaptive models produce accurate predictions even when microscopic details are unknown. The approach is computationally tractable, even for a relatively large number of dynamical variables. Using simulated data, it correctly infers the phase space structure for planetary motion, avoids overfitting in a biological signalling system and produces accurate predictions for yeast glycolysis with tens of data points and over half of the interacting species unobserved.

Mechanistic modelling of dynamical phenomena with many degrees of freedom runs the risk of overfitting and making faulty predictions, whereas ad hoc models may miss defining features. Here the authors develop an approach to construct dynamical models that adapt their complexity to the amount of available data.

Modeling for (physical) biologists: an introduction to the rule-based approach

Models that capture the chemical kinetics of cellular regulatory networks can be specified in terms of rules for biomolecular interactions. A rule defines a generalized reaction, meaning a reaction that permits multiple reactants, each capable of participating in a characteristic transformation and each possessing certain, specified properties, which may be local, such as the state of a particular site or domain of a protein. In other words, a rule defines a transformation and the properties that reactants must possess to participate in the transformation. A rule also provides a rate law. A rule-based approach to modeling enables consideration of mechanistic details at the level of functional sites of biomolecules and provides a facile and visual means for constructing computational models, which can be analyzed to study how system-level behaviors emerge from component interactions.

BioJazz: in silico evolution of cellular networks with unbounded complexity using rule-based modeling

Systems biologists aim to decipher the structure and dynamics of signaling and regulatory networks underpinning cellular responses; synthetic biologists can use this insight to alter existing networks or engineer de novo ones. Both tasks will benefit from an understanding of which structural and dynamic features of networks can emerge from evolutionary processes, through which intermediary steps these arise, and whether they embody general design principles. As natural evolution at the level of network dynamics is difficult to study, in silico evolution of network models can provide important insights. However, current tools used for in silico evolution of network dynamics are limited to ad hoc computer simulations and models. Here we introduce BioJazz, an extendable, user-friendly tool for simulating the evolution of dynamic biochemical networks. Unlike previous tools for in silico evolution, BioJazz allows for the evolution of cellular networks with unbounded complexity by combining rule-based modeling with an encoding of networks that is akin to a genome. We show that BioJazz can be used to implement biologically realistic selective pressures and allows exploration of the space of network architectures and dynamics that implement prescribed physiological functions. BioJazz is provided as an open-source tool to facilitate its further development and use. Source code and user manuals are available at: http://oss-lab.github.io/biojazz and http://osslab.lifesci.warwick.ac.uk/BioJazz.aspx.

Model reduction for stochastic CaMKII reaction kinetics in synapses by graph-constrained correlation dynamics

A stochastic reaction network model of Ca(2+) dynamics in synapses (Pepke et al PLoS Comput. Biol. 6 e1000675) is expressed and simulated using rule-based reaction modeling notation in dynamical grammars and in MCell. The model tracks the response of calmodulin and CaMKII to calcium influx in synapses. Data from numerically intensive simulations is used to train a reduced model that, out of sample, correctly predicts the evolution of interaction parameters characterizing the instantaneous probability distribution over molecular states in the much larger fine-scale models. The novel model reduction method, 'graph-constrained correlation dynamics', requires a graph of plausible state variables and interactions as input. It parametrically optimizes a set of constant coefficients appearing in differential equations governing the time-varying interaction parameters that determine all correlations between variables in the reduced model at any time slice.