BioNetGen 2.2: Advances in Rule-Based Modeling

BioNetGen is an open-source software package for rule-based modeling of complex biochemical systems. Version 2.2 of the software introduces numerous new features for both model specification and simulation. Here, we report on these additions, discussing how they facilitate the construction, simulation, and analysis of larger and more complex models than previously possible.Comment: 3 pages, 1 figure, 1 supplementary text file. Supplementary text includes a brief discussion of the RK-PLA along with a performance analysis, two tables listing all new actions/arguments added in BioNetGen 2.2, and the "BioNetGen Quick Reference Guide". Accepted for publication in Bioinformatic

An Introduction to Rule-based Modeling of Immune Receptor Signaling

Cells process external and internal signals through chemical interactions. Cells that constitute the immune system (e.g., antigen presenting cell, T-cell, B-cell, mast cell) can have different functions (e.g., adaptive memory, inflammatory response) depending on the type and number of receptor molecules on the cell surface and the specific intracellular signaling pathways activated by those receptors. Explicitly modeling and simulating kinetic interactions between molecules allows us to pose questions about the dynamics of a signaling network under various conditions. However, the application of chemical kinetics to biochemical signaling systems has been limited by the complexity of the systems under consideration. Rule-based modeling (BioNetGen, Kappa, Simmune, PySB) is an approach to address this complexity. In this chapter, by application to the Fc$\varepsilon$RI receptor system, we will explore the origins of complexity in macromolecular interactions, show how rule-based modeling can be used to address complexity, and demonstrate how to build a model in the BioNetGen framework. Open source BioNetGen software and documentation are available at http://bionetgen.org.Comment: 5 figure

Compact Representation of Photosynthesis Dynamics by Rule-based Models (Full Version)

Traditional mathematical models of photosynthesis are based on mass action kinetics of light reactions. This approach requires the modeller to enumerate all the possible state combinations of the modelled chemical species. This leads to combinatorial explosion in the number of reactions although the structure of the model could be expressed more compactly. We explore the use of rule-based modelling, in particular, a simplified variant of Kappa, to compactly capture and automatically reduce existing mathematical models of photosynthesis. Finally, the reduction procedure is implemented in BioNetGen language and demonstrated on several ODE models of photosynthesis processes. This is an extended version of the paper published in proceedings of 5th International Workshop on Static Analysis and Systems Biology (SASB) 2014.Comment: SASB 2014 full pape

Advances in Rule-based Modeling: Compartments, Energy, and Hybrid Simulation, with Application to Sepsis and Cell Signaling

Biological systems are commonly modeled as reaction networks, which describe the system at the resolution of biochemical species. Cellular systems, however, are governed by events at a finer scale: local interactions among macromolecular domains. The multi-domain structure of macromolecules, combined with the local nature of interactions, can lead to a combinatorial explosion that pushes reaction network methods to their limits. As an alternative, rule-based models (RBMs) describe the domain-based structure and local interactions found in biological systems. Molecular complexes are represented by graphs: functional domains as vertices, macromolecules as groupings of vertices, and molecular bonding as edges. Reaction rules, which describe classes of reactions, govern local modifications to molecular graphs, such as binding, post-translational modification, and degradation. RBMs can be transformed to equivalent reaction networks and simulated by differential or stochastic methods, or simulated directly with a network-free approach that avoids the problem of combinatorial complexity. Although RBMs and network-free methods resolve many problems in systems modeling, challenges remain. I address three challenges here: (i) managing model complexity due to cooperative interactions, (ii) representing biochemical systems in the compartmental setting of cells and organisms, and (iii) reducing the memory burden of large-scale network-free simulations. First, I present a general theory of energy-based modeling within the BioNetGen framework. Free energy is computed under a pattern-based formalism, and contextual variations within reaction classes are enumerated automatically. Next, I extend the BioNetGen language to permit description of compartmentalized biochemical systems, with treatment of volumes, surfaces and transport. Finally, a hybrid particle/population method is developed to reduce memory requirements of network-free simulations. All methods are implemented and available as part of BioNetGen. The remainder of this work presents an application to sepsis and inflammation. A multi-organ model of peritoneal infection and systemic inflammation is constructed and calibrated to experiment. Extra-corporeal blood purification, a potential treatment for sepsis, is explored in silico. Model simulations demonstrate that removal of blood cytokines and chemokines is a sufficient mechanism for improved survival in sepsis. However, differences between model predictions and the latest experimental data suggest directions for further exploration

Domain-Oriented Reduction of Rule-Based Network Models

The coupling of membrane-bound receptors to transcriptional regulators and other effector functions is mediated by multi-domain proteins that form complex assemblies. The modularity of protein interactions lends itself to a rule-based description, in which species and reactions are generated by rules that encode the necessary context for an interaction to occur, but also can produce a combinatorial explosion in the number of chemical species that make up the signaling network. We have shown previously that exact network reduction can be achieved using hierarchical control relationships between sites/domains on proteins to dissect multi-domain proteins into sets of non-interacting sites, allowing the replacement of each “full” (progenitor) protein with a set of derived auxiliary (offspring) proteins. The description of a network in terms of auxiliary proteins that have fewer sites than progenitor proteins often greatly reduces network size. We describe here a method for automating domain-oriented model reduction and its implementation as a module in the BioNetGen modeling package. It takes as input a standard BioNetGen model and automatically performs the following steps: 1) detecting the hierarchical control relationships between sites; 2) building up the auxiliary proteins; 3) generating a raw reduced model; and 4) cleaning up the raw model to provide the correct mass balance for each chemical species in the reduced network. We tested the performance of this module on models representing portions of growth factor receptor and immunoreceptor-mediated signaling networks, and confirmed its ability to reduce the model size and simulation cost by at least one or two orders of magnitude. Limitations of the current algorithm include the inability to reduce models based on implicit site dependencies or heterodimerization, and loss of accuracy when dynamics are computed stochastically