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.     

Hypothesis testing via integrated computer modeling and digital fluorescence microscopy

Computational modeling has the potential to add an entirely new approach to hypothesis testing in yeast cell biology. Here, we present a method for seamless integration of computational modeling with quantitative digital fluorescence microscopy. This integration is accomplished by developing computational models based on hypotheses for underlying cellular processes that may give rise to experimentally observed fluorescent protein localization patterns. Simulated fluorescence images are generated from the computational models of underlying cellular processes via a "model-convolution" process. These simulated images can then be directly compared to experimental fluorescence images in order to test the model. This method provides a framework for rigorous hypothesis testing in yeast cell biology via integrated mathematical modeling and digital fluorescence microscopy.     

Mathematical modeling as a tool for investigating cell cycle control networks

Although not a traditional experimental "method," mathematical modeling can provide a powerful approach for investigating complex cell signaling networks, such as those that regulate the eukaryotic cell division cycle. We describe here one modeling approach based on expressing the rates of biochemical reactions in terms of nonlinear ordinary differential equations. We discuss the steps and challenges in assigning numerical values to model parameters and the importance of experimental testing of a mathematical model. We illustrate this approach throughout with the simple and well-characterized example of mitotic cell cycles in frog egg extracts. To facilitate new modeling efforts, we describe several publicly available modeling environments, each with a collection of integrated programs for mathematical modeling. This review is intended to justify the place of mathematical modeling as a standard method for studying molecular regulatory networks and to guide the non-expert to initiate modeling projects in order to gain a systems-level perspective for complex control systems.     

Distributed sensor networks: a cellular nonlinear network perspective

Large-scale networks of integrated wireless sensors become increasingly tractable. Advances in hardware technology and engineering design have led to dramatic reductions in size, power consumption, and cost for digital circuitry, and wireless communications. Networking, self-organization, and distributed operation are crucial ingredients to harness the sensing, computing, and computational capabilities of the nodes into a complete system. This article shows that those networks can be considered as cellular nonlinear networks (CNNs), and that their analysis and design may greatly benefit from the rich theoretical results available for CNNs.     

基于约束的建模方法在代谢网络中的应用研究

代谢网络在各种细胞功能和生命过程中发挥着至关重要的作用。随着细胞网络重建工程的迅速发展,可用的基因组水平代谢网络越来越多,因而计算方法在这些网络的结构功能分析中越来越重要。基于约束的建模方法不像图论方法那样仅考虑代谢模型的纯拓扑结构,也不像各种动力学建模方法那样需求详尽的热力学参数,因而极具优势。采用基于约束的建模方法对一个含619个基因,655个代谢物和743个代谢反应的金黄色葡萄球菌(Staphylococcusaureus)代谢网络进行了分析,主要研究了该模型的网络结构特征,以及其最优生长率、动态生长情况和基因删除学习等。本研究提供了一个对金黄色葡萄球菌代谢网络进行约束建模分析的初步框架。     

A Population Density Approach That Facilitates Large-Scale Modeling of Neural Networks: Analysis and an Application to Orientation Tuning

We explore a computationally efficient method of simulating realistic networks of neurons introduced by Knight, Manin, and Sirovich (1996) in which integrate-and-fire neurons are grouped into large populations of similar neurons. For each population, we form a probability density that represents the distribution of neurons over all possible states. The populations are coupled via stochastic synapses in which the conductance of a neuron is modulated according to the firing rates of its presynaptic populations. The evolution equation for each of these probability densities is a partial differential-integral equation, which we solve numerically. Results obtained for several example networks are tested against conventional computations for groups of individual neurons.We apply this approach to modeling orientation tuning in the visual cortex. Our population density model is based on the recurrent feedback model of a hypercolumn in cat visual cortex of Somers et al. (1995). We simulate the response to oriented flashed bars. As in the Somers model, a weak orientation bias provided by feed-forward lateral geniculate input is transformed by intracortical circuitry into sharper orientation tuning that is independent of stimulus contrast.The population density approach appears to be a viable method for simulating large neural networks. Its computational efficiency overcomes some of the restrictions imposed by computation time in individual neuron simulations, allowing one to build more complex networks and to explore parameter space more easily. The method produces smooth rate functions with one pass of the stimulus and does not require signal averaging. At the same time, this model captures the dynamics of single-neuron activity that are missed in simple firing-rate models.     

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.     

Specification, annotation, visualization and simulation of a large rule-based model for ERBB receptor signaling

ABSTRACT: BACKGROUND: Mathematical/computational models are needed to understand cell signaling networks, which are complex. Signaling proteins contain multiple functional components and multiple sites of post-translational modification. The multiplicity of components and sites of modification ensures that interactions among signaling proteins have the potential to generate myriad protein complexes and post-translational modification states. As a result, the number of chemical species that can be populated in a cell signaling network, and hence the number of equations in an ordinary differential equation model required to capture the dynamics of these species, is prohibitively large. To overcome this problem, the rule-based modeling approach has been developed for representing interactions within signaling networks efficiently and compactly through coarse-graining of the chemical kinetics of molecular interactions. RESULTS: Here, we provide a demonstration that the rule-based modeling approach can be used to specify and simulate a large model for ERBB receptor signaling that accounts for site-specific details of protein-protein interactions. The model is considered large because it corresponds to a reaction network containing more reactions than can be practically enumerated. The model encompasses activation of ERK and Akt, and it can be simulated using a network-free simulator, such as NFsim, to generate time courses of phosphorylation for 55 individual serine, threonine, and tyrosine residues. The model is annotated and visualized in the form of an extended contact map. CONCLUSIONS: With the development of software that implements novel computational methods for calculating the dynamics of large-scale rule-based representations of cellular signaling networks, it is now possible to build and analyze models that include a significant fraction of the protein interactions that comprise a signaling network, with incorporation of the site-specific details of the interactions. Modeling at this level of detail is important for understanding cellular signaling.     

Prospects for Biotechnology in Oil Palm (Elaeis guineensis) and Coconut (Cocos nucifera) Improvement

Abstract

Research into human metabolism is expanding rapidly due to the emergence of metabolism as a key factor in common diseases. Mathematical modeling of human cellular metabolism has traditionally been performed via kinetic approaches whose applicability for large-scale systems is limited by lack of kinetic constants data. An alternative computational approach bypassing this hurdle called constraint-based modeling (CBM) serves to analyze the function of large-scale metabolic networks by solely relying on simple physical-chemical constraints. However, while extensive research has been performed on constraint-based modeling of microbial metabolism, large-scale modeling of human metabolism is still in its infancy. Utilizing constraint-based modeling to model human cellular metabolism is significantly more complicated than modeling microbial metabolism as in multi-cellular organisms the metabolic behavior varies across cell-types and tissues. It is further complicated due to lack of data on cell type- and tissue-specific metabolite uptake from the surrounding microenvironments and tissue-specific metabolic objective functions. To overcome these problems, several studies suggested CBM methods that integrate metabolic networks with gene expression data that is easily measurable under various conditions. This specific objective functions are expected to improve the prediction accuracy of the presented methods. Such objective functions may be derived based on computational learning that would give optimal correspondence between predicted and measured metabolic phenotypes (Burgard, 2003).

The CBM methods presented here open the way for future computational investigations of metabolic disorders given the relevant expression data. A first attempt to visualize and interpret changes in gene expression data measured following gastric bypass surgery via a genome-scale metabolic network was done by Duarte et al (Duarte, 2007). Another potential application would be the prediction of diagnostic biomarkers for metabolic diseases that could be identified via biofluid metabolomics (Kell, 2007). Towards this goal, we have recently developed a CBM method for predicting metabolic biomarkers for in-born errors of metabolism by searching for changes in metabolite uptake and secretion rate due to genetic alterations (Shlomi, 2009). Incorporating cell type- and tissue-specific gene expression data within this framework can potentially improve the identification of diagnostic biomarkers. Overall, the methods presented here lay the foundation for studying normal and abnormal human cellular metabolism in tissue-specific manner based on commonly measured gene expression data.     

A link between neuroscience and informatics: large-scale modeling of memory processes

Utilizing advances in functional neuroimaging and computational neural modeling, neuroscientists have increasingly sought to investigate how distributed networks, composed of functionally defined subregions, combine to produce cognition. Large-scale, biologically realistic neural models, which integrate data from cellular, regional, whole brain, and behavioral sources, delineate specific hypotheses about how these interacting neural populations might carry out high-level cognitive tasks. In this review, we discuss neuroimaging, neural modeling, and the utility of large-scale biologically realistic models using modeling of short-term memory as an example. We present a sketch of the data regarding the neural basis of short-term memory from non-human electrophysiological, computational and neuroimaging perspectives, highlighting the multiple interacting brain regions believed to be involved. Through a review of several efforts, including our own, to combine neural modeling and neuroimaging data, we argue that large scale neural models provide specific advantages in understanding the distributed networks underlying cognition and behavior.     

Simulating Metabolism with Statistical Thermodynamics

New methods are needed for large scale modeling of metabolism that predict metabolite levels and characterize the thermodynamics of individual reactions and pathways. Current approaches use either kinetic simulations, which are difficult to extend to large networks of reactions because of the need for rate constants, or flux-based methods, which have a large number of feasible solutions because they are unconstrained by the law of mass action. This report presents an alternative modeling approach based on statistical thermodynamics. The principles of this approach are demonstrated using a simple set of coupled reactions, and then the system is characterized with respect to the changes in energy, entropy, free energy, and entropy production. Finally, the physical and biochemical insights that this approach can provide for metabolism are demonstrated by application to the tricarboxylic acid (TCA) cycle of Escherichia coli. The reaction and pathway thermodynamics are evaluated and predictions are made regarding changes in concentration of TCA cycle intermediates due to 10- and 100-fold changes in the ratio of NAD⁺:NADH concentrations. Finally, the assumptions and caveats regarding the use of statistical thermodynamics to model non-equilibrium reactions are discussed.     

The ion channel inverse problem: neuroinformatics meets biophysics

Ion channels are the building blocks of the information processing capability of neurons: any realistic computational model of a neuron must include reliable and effective ion channel components. Sophisticated statistical and computational tools have been developed to study the ion channel structure–function relationship, but this work is rarely incorporated into the models used for single neurons or small networks. The disjunction is partly a matter of convention. Structure–function studies typically use a single Markov model for the whole channel whereas until recently whole-cell modeling software has focused on serial, independent, two-state subunits that can be represented by the Hodgkin–Huxley equations. More fundamentally, there is a difference in purpose that prevents models being easily reused. Biophysical models are typically developed to study one particular aspect of channel gating in detail, whereas neural modelers require broad coverage of the entire range of channel behavior that is often best achieved with approximate representations that omit structural features that cannot be adequately constrained. To bridge the gap so that more recent channel data can be used in neural models requires new computational infrastructure for bringing together diverse sources of data to arrive at best-fit models for whole-cell modeling. We review the current state of channel modeling and explore the developments needed for its conclusions to be integrated into whole-cell modeling.     

Bridging the gaps in systems biology

Systems biology aims at creating mathematical models, i.e., computational reconstructions of biological systems and processes that will result in a new level of understanding—the elucidation of the basic and presumably conserved “design” and “engineering” principles of biomolecular systems. Thus, systems biology will move biology from a phenomenological to a predictive science. Mathematical modeling of biological networks and processes has already greatly improved our understanding of many cellular processes. However, given the massive amount of qualitative and quantitative data currently produced and number of burning questions in health care and biotechnology needed to be solved is still in its early phases. The field requires novel approaches for abstraction, for modeling bioprocesses that follow different biochemical and biophysical rules, and for combining different modules into larger models that still allow realistic simulation with the computational power available today. We have identified and discussed currently most prominent problems in systems biology: (1) how to bridge different scales of modeling abstraction, (2) how to bridge the gap between topological and mechanistic modeling, and (3) how to bridge the wet and dry laboratory gap. The future success of systems biology largely depends on bridging the recognized gaps.     

Design constraints on a synthetic metabolism

A metabolism is a complex network of chemical reactions that converts sources of energy and chemical elements into biomass and other molecules. To design a metabolism from scratch and to implement it in a synthetic genome is almost within technological reach. Ideally, a synthetic metabolism should be able to synthesize a desired spectrum of molecules at a high rate, from multiple different nutrients, while using few chemical reactions, and producing little or no waste. Not all of these properties are achievable simultaneously. We here use a recently developed technique to create random metabolic networks with pre-specified properties to quantify trade-offs between these and other properties. We find that for every additional molecule to be synthesized a network needs on average three additional reactions. For every additional carbon source to be utilized, it needs on average two additional reactions. Networks able to synthesize 20 biomass molecules from each of 20 alternative sole carbon sources need to have at least 260 reactions. This number increases to 518 reactions for networks that can synthesize more than 60 molecules from each of 80 carbon sources. The maximally achievable rate of biosynthesis decreases by approximately 5 percent for every additional molecule to be synthesized. Biochemically related molecules can be synthesized at higher rates, because their synthesis produces less waste. Overall, the variables we study can explain 87 percent of variation in network size and 84 percent of the variation in synthesis rate. The constraints we identify prescribe broad boundary conditions that can help to guide synthetic metabolism design.     

Thermodynamic-based computational profiling of cellular regulatory control in hepatocyte metabolism

Thermodynamic-based constraints on biochemical fluxes and concentrations are applied in concert with mass balance of fluxes in glycogenesis and glycogenolysis in a model of hepatic cell metabolism. Constraint-based modeling methods that facilitate predictions of reactant concentrations, reaction potentials, and enzyme activities are introduced to identify putative regulatory and control sites in biological networks by computing the minimal control scheme necessary to switch between metabolic modes. Computational predictions of control sites in glycogenic and glycogenolytic operational modes in the hepatocyte network compare favorably with known regulatory mechanisms. The developed hepatic metabolic model is used to computationally analyze the impairment of glucose production in von Gierke's and Hers' diseases, two metabolic diseases impacting glycogen metabolism. The computational methodology introduced here can be generalized to identify downstream targets of agonists, to systematically probe possible drug targets, and to predict the effects of specific inhibitors (or activators) on integrated network function.