Leveraging modeling approaches: reaction networks and rules

Optimal control of agent-based models via surrogate modeling

This paper describes and validates an algorithm to solve optimal control problems for agent-based models (ABMs). For a given ABM and a given optimal control problem, the algorithm derives a surrogate model, typically lower-dimensional, in the form of a system of ordinary differential equations (ODEs), solves the control problem for the surrogate model, and then transfers it back to the original ABM. It applies to quite general ABMs and offers several options for the ODE structure, depending on what information about the ABM is to be used. There is a broad range of applications for such an algorithm, since ABMs are used widely in the life sciences, such as ecology, epidemiology, and biomedicine and healthcare, areas where optimal control is an important purpose for modeling, such as for medical digital twin technology.

Biosystem models, generated from a complex rule/reaction/influence network and from two functionality prototypes

In this work we have further developed the Direct Computer Mapping (DCM) based modelling and simulation methodology. A unified, transition-based representation of complex rule, reaction and influence networks has been introduced and two prototypes (one general state- and another general transition-prototype) have been developed for the unified functional modelling of the state and transition nodes. Starting from the network and from the functional prototypes, an automatic generation method of the graphically editable and extensible GraphML description of biosystem models has been elaborated. The new developments have been implemented in the improved kernel of DCM models. The applied knowledge representation makes possible the unified generation and execution of the balance-based quantitative and influence- or rule-based qualitative, as well as optionally time-driven, multiscale biosystem models. Application of the developed methodology has been illustrated by the improved implementation of the formerly studied and upgraded example biosystem model for combining the detailed, quantitative p53/miR34a signalling system with the pathological model through an extended rule-based coupling model.

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.

Rule-based modeling and simulations of the inner kinetochore structure

BACKGROUND

Combinatorial complexity is a central problem when modeling biochemical reaction networks, since the association of a few components can give rise to a large variation of protein complexes. Available classical modeling approaches are often insufficient for the analysis of very large and complex networks in detail. Recently, we developed a new rule-based modeling approach that facilitates the analysis of spatial and combinatorially complex problems. Here, we explore for the first time how this approach can be applied to a specific biological system, the human kinetochore, which is a multi-protein complex involving over 100 proteins.

RESULTS

Applying our freely available SRSim software to a large data set on kinetochore proteins in human cells, we construct a spatial rule-based simulation model of the human inner kinetochore. The model generates an estimation of the probability distribution of the inner kinetochore 3D architecture and we show how to analyze this distribution using information theory. In our model, the formation of a bridge between CenpA and an H3 containing nucleosome only occurs efficiently for higher protein concentration realized during S-phase but may be not in G1. Above a certain nucleosome distance the protein bridge barely formed pointing towards the importance of chromatin structure for kinetochore complex formation. We define a metric for the distance between structures that allow us to identify structural clusters. Using this modeling technique, we explore different hypothetical chromatin layouts.

CONCLUSIONS

Applying a rule-based network analysis to the spatial kinetochore complex geometry allowed us to integrate experimental data on kinetochore proteins, suggesting a 3D model of the human inner kinetochore architecture that is governed by a combinatorial algebraic reaction network. This reaction network can serve as bridge between multiple scales of modeling. Our approach can be applied to other systems beyond kinetochores.