Search

Time- and state-discrete dynamical systems are frequently used to model molecular networks. This paper provides a collection of mathematical and computational tools for the study of robustness in Boolean network models. The focus is on networks governed by k -canalizing functions, a recently introduced class of Boolean functions that contains the well-studied class of nested canalizing functions. The variable activities and sensitivity of a function quantify the impact of input changes on the function output. This paper generalizes the latter concept to c -sensitivity and provides formulas for the activities and c -sensitivity of general k -canalizing functions as well as canalizing functions with more precisely defined structure. A popular measure for the robustness of a network, the Derrida value, can be expressed as a weighted sum of the c -sensitivities of the governing canalizing functions, and can also be calculated for a stochastic extension of Boolean networks. These findings provide a computationally efficient way to obtain Derrida values of Boolean networks, deterministic or stochastic, that does not involve simulation. [ABSTRACT FROM AUTHOR]

This paper presents a modular software design for the construction of computational modeling technology that will help implement precision medicine. In analogy to a common industrial strategy used for preventive maintenance of engineered products, medical digital twins are computational models of disease processes calibrated to individual patients using multiple heterogeneous data streams. They have the potential to help improve diagnosis, prognosis, and personalized treatment for a wide range of medical conditions. Their large-scale development relies on both mechanistic and datadriven techniques and requires the integration and ongoing update of multiple component models developed across many different laboratories. Distributed model building and integration requires an open-source modular software platform for the integration and simulation of models that is scalable and supports a decentralized, community-based model building process. This paper presents such a platform, including a case study in an animal model of a respiratory fungal infection. [ABSTRACT FROM AUTHOR]

Agent-based models (ABMs) have become an increasingly important mode of inquiry for the life sciences. They are particularly valuable for systems that are not understood well enough to build an equation-based model. These advantages, however, are counterbalanced by the difficulty of analyzing and using ABMs, due to the lack of the type of mathematical tools available for more traditional models, which leaves simulation as the primary approach. As models become large, simulation becomes challenging. This paper proposes a novel approach to two mathematical aspects of ABMs, optimization and control, and it presents a few first steps outlining how one might carry out this approach. Rather than viewing the ABM as a model, it is to be viewed as a surrogate for the actual system. For a given optimization or control problem (which may change over time), the surrogate system is modeled instead, using data from the ABM and a modeling framework for which ready-made mathematical tools exist, such as differential equations, or for which control strategies can explored more easily. Once the optimization problem is solved for the model of the surrogate, it is then lifted to the surrogate and tested. The final step is to lift the optimization solution from the surrogate system to the actual system. This program is illustrated with published work, using two relatively simple ABMs as a demonstration, Sugarscape and a consumer-resource ABM. Specific techniques discussed include dimension reduction and approximation of an ABM by difference equations as well systems of PDEs, related to certain specific control objectives. This demonstration illustrates the very challenging mathematical problems that need to be solved before this approach can be realistically applied to complex and large ABMs, current and future. The paper outlines a research program to address them. [ABSTRACT FROM AUTHOR]

The global dynamics of gene regulatory networks are known to show robustness to perturbations in the form of intrinsic and extrinsic noise, as well as mutations of individual genes. One molecular mechanism underlying this robustness has been identified as the action of so-called microRNAs that operate via feedforward loops. We present results of a computational study, using the modeling framework of stochastic Boolean networks, which explores the role that such network motifs play in stabilizing global dynamics. The paper introduces a new measure for the stability of stochastic networks. The results show that certain types of feedforward loops do indeed buffer the network against stochastic effects. [ABSTRACT FROM AUTHOR]

Abstract: Iron is a metal essential for cellular metabolism. However, excess iron available for reactions contributes to the formation of dangerous reactive oxygen species, such as the hydroxyl radical, via the Fenton reaction. Therefore, intracellular iron levels are tightly constrained by a control system of proteins. This paper contains a mathematical model, in the form of a system of five ordinary differential equations, of the core of this control system, including the labile iron pool as well as proteins that regulate uptake, storage, and export and are connected through negative feedback loops. The model is validated using data from an overexpression experiment with cultured human breast epithelial cells. The parameters in the mathematical model are not known for this particular cell culture system, so the analysis of the model was done for a generic choice of parameters. Through a mixture of analytical arguments and extensive simulations it is shown that for any choice of parameters the model reaches a unique stable steady state, thereby ruling out oscillatory behavior. It is shown furthermore that the model parameters are identifiable through suitable experiments. [Copyright &y& Elsevier]

Abstract: Epstein–Barr virus (EBV) is an important human pathogen that establishes a life-long persistent infection and for which no precise animal model exists. In this paper, we describe in detail an agent-based model and computer simulation of EBV infection. Agents representing EBV and sets of B and T lymphocytes move and interact on a three-dimensional grid approximating Waldeyer''s ring, together with abstract compartments for lymph and blood. The simulation allows us to explore the development and resolution of virtual infections in a manner not possible in actual human experiments. Specifically, we identify parameters capable of inducing clearance, persistent infection, or death. [Copyright &y& Elsevier]

Microarray technologies, which can measure tens of thousands of gene expression values simultaneously in a single experiment, have become a common research method for biomedical researchers. Computational tools to analyze microarray data for biological discovery are needed. In this paper, we investigate the feasibility of using formal concept analysis (FCA) as a tool for microarray data analysis. The method of FCA builds a (concept) lattice from the experimental data together with additional biological information. For microarray data, each vertex of the lattice corresponds to a subset of genes that are grouped together according to their expression values and some biological information related to gene function. The lattice structure of these gene sets might reflect biological relationships in the dataset. Similarities and differences between experiments can then be investigated by comparing their corresponding lattices according to various graph measures. We apply our method to microarray data derived from influenza-infected mouse lung tissue and healthy controls. Our preliminary results show the promise of our method as a tool for microarray data analysis. [ABSTRACT FROM AUTHOR]

We consider the problem of reverse-engineering dynamic models of biochemical networks from experimental data using polynomial dynamic systems. In earlier work, we developed an algorithm to identify minimal wiring diagrams, that is, directed graphs that represent the causal relationships between network variables. Here we extend this algorithm to identify a most likely dynamic model from the set of all possible dynamic models that fit the data over a fixed wiring diagram. To illustrate its performance, the method is applied to simulated time-course data from a published gene regulatory network in the fruitfly Drosophila melanogaster. [ABSTRACT FROM AUTHOR]

Background: At the molecular level, nonlinear networks of heterogeneous molecules control many biological processes, so that systems biology provides a valuable approach in this field, building on the integration of experimental biology with mathematical modeling. One of the biggest challenges to making this integration a reality is that many life scientists do not possess the mathematical expertise needed to build and manipulate mathematical models well enough to use them as tools for hypothesis generation. Available modeling software packages often assume some modeling expertise. There is a need for software tools that are easy to use and intuitive for experimentalists. Results: This paper introduces PlantSimLab, a web-based application developed to allow plant biologists to construct dynamic mathematical models of molecular networks, interrogate them in a manner similar to what is done in the laboratory, and use them as a tool for biological hypothesis generation. It is designed to be used by experimentalists, without direct assistance from mathematical modelers. Conclusions: Mathematical modeling techniques are a useful tool for analyzing complex biological systems, and there is a need for accessible, efficient analysis tools within the biological community. PlantSimLab enables users to build, validate, and use intuitive qualitative dynamic computer models, with a graphical user interface that does not require mathematical modeling expertise. It makes analysis of complex models accessible to a larger community, as it is platform-independent and does not require extensive mathematical expertise. [ABSTRACT FROM AUTHOR]