Strategies for dealing with incomplete information in the modeling of molecular interaction networks

Modelers of molecular interaction networks encounter the paradoxical situation that while large amounts of data are available, these are often insufficient for the formulation and analysis of mathematical models describing the network dynamics. In particular, information on the reaction mechanisms and numerical values of kinetic parameters are usually not available for all but a few well-studied model systems. In this article we review two strategies that have been proposed for dealing with incomplete information in the study of molecular interaction networks: parameter sensitivity analysis and model simplification. These strategies are based on the biologically justified intuition that essential properties of the system dynamics are robust against moderate changes in the value of kinetic parameters or even in the rate laws describing the interactions. Although advanced measurement techniques can be expected to relieve the problem of incomplete information to some extent, the strategies discussed in this article will retain their interest as tools providing an initial characterization of essential properties of the network dynamics.     

A survey on methods for modeling and analyzing integrated biological networks

Understanding how cellular systems build up integrated responses to their dynamically changing environment is one of the open questions in Systems Biology. Despite their intertwinement, signaling networks, gene regulation and metabolism have been frequently modeled independently in the context of well-defined subsystems. For this purpose, several mathematical formalisms have been developed according to the features of each particular network under study. Nonetheless, a deeper understanding of cellular behavior requires the integration of these various systems into a model capable of capturing how they operate as an ensemble. With the recent advances in the "omics" technologies, more data is becoming available and, thus, recent efforts have been driven toward this integrated modeling approach. We herein review and discuss methodological frameworks currently available for modeling and analyzing integrated biological networks, in particular metabolic, gene regulatory and signaling networks. These include network-based methods and Chemical Organization Theory, Flux-Balance Analysis and its extensions, logical discrete modeling, Petri Nets, traditional kinetic modeling, Hybrid Systems and stochastic models. Comparisons are also established regarding data requirements, scalability with network size and computational burden. The methods are illustrated with successful case studies in large-scale genome models and in particular subsystems of various organisms.     

Towards kinetic modeling of genome-scale metabolic networks without sacrificing stoichiometric,thermodynamic and physiological constraints

Mathematical modeling is an essential tool for the comprehensive understanding of cell metabolism and its interactions with the environmental and process conditions. Recent developments in the construction and analysis of stoichiometric models made it possible to define limits on steady-state metabolic behavior using flux balance analysis. However, detailed information on enzyme kinetics and enzyme regulation is needed to formulate kinetic models that can accurately capture the dynamic metabolic responses. The use of mechanistic enzyme kinetics is a difficult task due to uncertainty in the kinetic properties of enzymes. Therefore, the majority of recent works considered only mass action kinetics for reactions in metabolic networks. Herein, we applied the optimization and risk analysis of complex living entities (ORACLE) framework and constructed a large-scale mechanistic kinetic model of optimally grown Escherichia coli. We investigated the complex interplay between stoichiometry, thermodynamics, and kinetics in determining the flexibility and capabilities of metabolism. Our results indicate that enzyme saturation is a necessary consideration in modeling metabolic networks and it extends the feasible ranges of metabolic fluxes and metabolite concentrations. Our results further suggest that enzymes in metabolic networks have evolved to function at different saturation states to ensure greater flexibility and robustness of cellular metabolism.     

RMBNToolbox: random models for biochemical networks

Background  

There is an increasing interest to model biochemical and cell biological networks, as well as to the computational analysis of these models. The development of analysis methodologies and related software is rapid in the field. However, the number of available models is still relatively small and the model sizes remain limited. The lack of kinetic information is usually the limiting factor for the construction of detailed simulation models.     

Mathematical optimization applications in metabolic networks

Genome-scale metabolic models are increasingly becoming available for a variety of microorganisms. This has spurred the development of a wide array of computational tools, and in particular, mathematical optimization approaches, to assist in fundamental metabolic network analyses and redesign efforts. This review highlights a number of optimization-based frameworks developed towards addressing challenges in the analysis and engineering of metabolic networks. In particular, three major types of studies are covered here including exploring model predictions, correction and improvement of models of metabolism, and redesign of metabolic networks for the targeted overproduction of a desired compound. Overall, the methods reviewed in this paper highlight the diversity of queries, breadth of questions and complexity of redesigns that are amenable to mathematical optimization strategies.     

Current Applications of Models of Genetic Effects with Interactions Across the Genome

Models of genetic effects integrate the action of genes, regulatory regions and interactions among alleles across the genome. Such theoretical frameworks are critical for applied studies in at least two ways. First, discovering genetic networks with specific effects underlying traits in populations requires the development of models that implement those effects as parameters—adjusting the implementation of epistasis parameters in genetic models has for instance been a requirement for properly testing for epistasis in gene-mapping studies. Second, studying the properties and implications of models of genetic effects that involve complex genetic networks has proven to be valuable, whether those networks have been revealed for particular organisms or inferred to be of interest from theoretical works and simulations. Here I review the current state of development and recent applications of models of genetic effects. I focus on general models aiming to depict complex genotype-to-phenotype maps and on applications of them to networks of interacting loci.     

Shaping dots and lines: adding modularity into protein interaction networks using structural information

Determining protein interaction networks and generating models to simulate network changes in time and space are crucial for understanding a biological system and for predicting the effect of mutants found in diseases. In this review we discuss the great potential of using structural information together with computational tools towards reaching this goal: the prediction of new protein interactions, the estimation of affinities and kinetic rate constants between protein complexes, and finally the determination of which interactions are compatible with each other and which interactions are exclusive. The latter one will be important to reorganize large scale networks into functional modular networks.     

The construction of an amino acid network for understanding protein structure and function

Amino acid networks (AANs) are undirected networks consisting of amino acid residues and their interactions in three-dimensional protein structures. The analysis of AANs provides novel insight into protein science, and several common amino acid network properties have revealed diverse classes of proteins. In this review, we first summarize methods for the construction and characterization of AANs. We then compare software tools for the construction and analysis of AANs. Finally, we review the application of AANs for understanding protein structure and function, including the identification of functional residues, the prediction of protein folding, analyzing protein stability and protein–protein interactions, and for understanding communication within and between proteins.     

Evaluating Communities of Practice and Knowledge Networks: A Systematic Scoping Review of Evaluation Frameworks

Communities of Practice (CoPs) are increasingly considered a part of ecohealth and other sectors such as health care, education, and business. However, there is little agreement on approaches to evaluate the influence and effectiveness of CoPs. The purpose of this review was to understand what frameworks and methods have been proposed or used to evaluate CoPs and/or knowledge networks. The review searched electronic databases in interdisciplinary, health, education, and business fields, and further collected references and forward citations from relevant articles. Nineteen articles with 16 frameworks were included in the synthesis. The purposes of the evaluation frameworks varied; while some focused on assessing the performance of CoPs, several frameworks sought to learn about CoPs and their critical success factors. Nine of the frameworks had been applied or tested in some way, most frequently to guide a case study. With limited applications of the frameworks, strong claims about generalizability could not be made. The review results can inform the development of tailored frameworks. However, there is a need for more detailed and targeted CoP evaluation frameworks, as many imperative CoP evaluation needs would be unmet by the available frameworks.     

Modeling biological systems using Dynetica--a simulator of dynamic networks

We present Dynetica, a user-friendly simulator of dynamic networks for constructing, visualizing, and analyzing kinetic models of biological systems. In addition to generic reaction networks, Dynetica facilitates construction of models of genetic networks, where many reactions are gene expression and interactions among gene products. Further, it integrates the capability of conducting both deterministic and stochastic simulations. AVAILABILITY AND SUPPLEMENTARY INFORMATION: Dynetica 1.0, example models, and the user's guide are available at http://www.its.caltech.edu/~you/Dynetica/Dynetica_page.htm     

Non-linear reduction for kinetic models of metabolic reaction networks

Kinetic models of metabolic networks are essential for predicting and optimizing the transient behavior of cells in culture. However, such models are inherently high dimensional and stiff due to the large number of species and reactions involved and to kinetic rate constants of widely different orders of magnitude. In this paper we address the problem of deriving non-stiff, reduced-order non-linear models of the dominant dynamics of metabolic networks with fast and slow reactions. We present a method, based on singular perturbation analysis, which allows the systematic identification of quasi-steady-state conditions for the fast reactions, and the derivation of explicit non-linear models of the slow dynamics independent of the fast reaction rate expressions. The method is successfully applied to detailed models of metabolism in human erythrocytes and Saccharomyces cerevisiae.     

k-Cone analysis: determining all candidate values for kinetic parameters on a network scale

The absence of comprehensive measured kinetic values and the observed inconsistency in the available in vitro kinetic data has hindered the formulation of network-scale kinetic models of biochemical reaction networks. To meet this challenge we present an approach to construct a convex space, termed the k-cone, which contains all the allowable numerical values of the kinetic constants in large-scale biochemical networks. The definition of the k-cone relies on the incorporation of in vivo concentration data and a simplified approach to represent enzyme kinetics within an established constraint-based modeling approach. The k-cone approach was implemented to define the allowable combination of numerical values for a full kinetic model of human red blood cell metabolism and to study its correlated kinetic parameters. The k-cone approach can be used to determine consistency between in vitro measured kinetic values and in vivo concentration and flux measurements when used in a network-scale kinetic model. k-Cone analysis was successful in determining whether in vitro measured kinetic values used in the reconstruction of a kinetic-based model of Saccharomyces cerevisiae central metabolism could reproduce in vivo measurements. Further, the k-cone can be used to determine which numerical values of in vitro measured parameters are required to be changed in a kinetic model if in vivo measured values are not reproduced. k-Cone analysis could identify what minimum number of in vitro determined kinetic parameters needed to be adjusted in the S. cerevisiae model to be consistent with the in vivo data. Applying the k-cone analysis a priori to kinetic model development may reduce the time and effort involved in model building and parameter adjustment. With the recent developments in high-throughput profiling of metabolite concentrations at a whole-cell scale and advances in metabolomics technologies, the k-cone approach presented here may hold the promise for kinetic characterization of metabolic networks as well as other biological functions at a whole-cell level.     

Modeling and analysis of biopathways dynamics

Cellular processes are governed and coordinated by a multitude of biopathways. A pathway can be viewed as a complex network of biochemical reactions. The dynamics of this network largely determines the functioning of the pathway. Hence the modeling and analysis of biochemical networks dynamics is an important problem and is an active area of research. Here we review quantitative models of biochemical networks based on ordinary differential equations (ODEs). We mainly focus on the parameter estimation and sensitivity analysis problems and survey the current methods for tackling them. In this context we also highlight a recently developed probabilistic approximation technique using which these two problems can be considerably simplified.