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

RNA packaging motor: From structure to quantum mechanical modelling and sequential-stochastic mechanism

The bacteriophages of the Cystoviridae family package their single stranded RNA genomic precursors into empty capsid (procapsids) using a hexameric packaging ATPase motor (P4). This molecular motor shares sequence and structural similarity with RecA-like hexameric helicases. A concerted structural, mutational and kinetic analysis helped to define the mechanical reaction coordinate, i.e. the conformational changes associated with RNA translocation. The results also allowed us to propose a possible scheme of coupling between ATP hydrolysis and translocation which requires the cooperative action of three consecutive subunits. Here, we first test this model by preparing hexamers with defined proportions of wild type and mutant subunits and measuring their activity. Then, we develop a stochastic kinetic model which accounts for the catalytic cooperativity of the P4 hexamer. Finally, we use the available structural information to construct a quantum-chemical model of the chemical reaction coordinate and obtain a detailed description of the electron density changes during ATP hydrolysis. The model explains the results of the mutational analyses and yields new insights into the role of several conserved residues within the ATP binding pocket. These hypotheses will guide future experimental work

COEL: A Web-based Chemistry Simulation Framework

The chemical reaction network (CRN) is a widely used formalism to describe macroscopic behavior of chemical systems. Available tools for CRN modelling and simulation require local access, installation, and often involve local file storage, which is susceptible to loss, lacks searchable structure, and does not support concurrency. Furthermore, simulations are often single-threaded, and user interfaces are non-trivial to use. Therefore there are significant hurdles to conducting efficient and collaborative chemical research. In this paper, we introduce a new enterprise chemistry simulation framework, COEL, which addresses these issues. COEL is the first web-based framework of its kind. A visually pleasing and intuitive user interface, simulations that run on a large computational grid, reliable database storage, and transactional services make COEL ideal for collaborative research and education. COEL's most prominent features include ODE-based simulations of chemical reaction networks and multicompartment reaction networks, with rich options for user interactions with those networks. COEL provides DNA-strand displacement transformations and visualization (and is to our knowledge the first CRN framework to do so), GA optimization of rate constants, expression validation, an application-wide plotting engine, and SBML/Octave/Matlab export. We also present an overview of the underlying software and technologies employed and describe the main architectural decisions driving our development. COEL is available at http://coel-sim.org for selected research teams only. We plan to provide a part of COEL's functionality to the general public in the near future.Comment: 23 pages, 12 figures, 1 tabl

Kinetic Intersections as a Source of Information on Rate Coefficients

C

CFD modelling of two-phase stirred bioreaction systems by segregated solution of the Euler–Euler model

An advanced study of a bioreactor system involving a Navier–Stokes based model has been accomplished. The model allows a more realistic impeller induced flow image to be combined with the Monod bioreaction kinetics reported previously. The time-course of gluconic acid production by Aspergillus niger strain is simulated at kinetic conditions proposed in the literature. The simulation is based on (1) a stepwise solution strategy resolving first the fluid flow field, further imposing oxygen mass transfer and bioreaction with subsequent analysis of flow interactions, and (2) a segregated solution of the model replacing the multiple iterations per grid cell with single iterations. The numerical results are compared with experimental data for the bioreaction dynamics and show satisfactory agreement. The model is used for assessment of the viscosity effect upon the bioreactor performance. A 10-fold viscosity rise results in 2-fold decrease of KLa and 25% decrease of the specific gluconic acid production rate. The model allows better understanding of the mechanism of the important bioprocess

A structured approach for the engineering of biochemical network models, illustrated for signalling pathways

http://dx.doi.org/10.1093/bib/bbn026Quantitative models of biochemical networks (signal transduction cascades, metabolic pathways, gene regulatory circuits) are a central component of modern systems biology. Building and managing these complex models is a major challenge that can benefit from the application of formal methods adopted from theoretical computing science. Here we provide a general introduction to the field of formal modelling, which emphasizes the intuitive biochemical basis of the modelling process, but is also accessible for an audience with a background in computing science and/or model engineering. We show how signal transduction cascades can be modelled in a modular fashion, using both a qualitative approach { Qualitative Petri nets, and quantitative approaches { Continuous Petri Nets and Ordinary Differential Equations. We review the major elementary building blocks of a cellular signalling model, discuss which critical design decisions have to be made during model building, and present ..

The 1999 Center for Simulation of Dynamic Response in Materials Annual Technical Report

Introduction: This annual report describes research accomplishments for FY 99 of the Center for Simulation of Dynamic Response of Materials. The Center is constructing a virtual shock physics facility in which the full three dimensional response of a variety of target materials can be computed for a wide range of compressive, ten- sional, and shear loadings, including those produced by detonation of energetic materials. The goals are to facilitate computation of a variety of experiments in which strong shock and detonation waves are made to impinge on targets consisting of various combinations of materials, compute the subsequent dy- namic response of the target materials, and validate these computations against experimental data

Exploration of Reaction Pathways and Chemical Transformation Networks

For the investigation of chemical reaction networks, the identification of all relevant intermediates and elementary reactions is mandatory. Many algorithmic approaches exist that perform explorations efficiently and automatedly. These approaches differ in their application range, the level of completeness of the exploration, as well as the amount of heuristics and human intervention required. Here, we describe and compare the different approaches based on these criteria. Future directions leveraging the strengths of chemical heuristics, human interaction, and physical rigor are discussed.Comment: 48 pages, 4 figure