Towards building the silicon cell: a modular approach

Systems Biology aims to understand quantitatively how properties of biological systems can be understood as functions of the characteristics of, and interactions between their macromolecular components. Whereas, traditional biochemistry focused on isolation and characterization of cellular components, the challenge for Systems Biology lies in integration of this knowledge and the knowledge about molecular interactions. Computer models play an important role in this integration. We here discuss an approach with which we aim to link kinetic models on small parts of metabolism together, so as to form detailed kinetic models of larger chunks of metabolism, and ultimately of the entire living cell. Specifically, we will discuss techniques that can be used to model a sub-network in isolation of a larger network of which it is a part, while still maintaining the dynamics of the larger complete network. We will start by outlining the JWS online system, the silicon cell project, and the type of models we propose. JWS online is a model repository, which can be used for the storage, simulation and analysis of kinetic models. We advocate to integrate a top-down approach, where measurements on the complete system are used to derive fluxes in a detailed structural model, with a bottom-up approach, consisting of the integration of molecular mechanism-based detailed kinetic models into the structural model.     

Dynamic simulation of regulatory networks using SQUAD

Background  

The ambition of most molecular biologists is the understanding of the intricate network of molecular interactions that control biological systems. As scientists uncover the components and the connectivity of these networks, it becomes possible to study their dynamical behavior as a whole and discover what is the specific role of each of their components. Since the behavior of a network is by no means intuitive, it becomes necessary to use computational models to understand its behavior and to be able to make predictions about it. Unfortunately, most current computational models describe small networks due to the scarcity of kinetic data available. To overcome this problem, we previously published a methodology to convert a signaling network into a dynamical system, even in the total absence of kinetic information. In this paper we present a software implementation of such methodology.     

A framework for integrating the songbird brain

Biological systems by default involve complex components with complex relationships. To decipher how biological systems work, we assume that one needs to integrate information over multiple levels of complexity. The songbird vocal communication system is ideal for such integration due to many years of ethological investigation and a discreet dedicated brain network. Here we announce the beginnings of a songbird brain integrative project that involves high-throughput, molecular, anatomical, electrophysiological and behavioral levels of analysis. We first formed a rationale for inclusion of specific biological levels of analysis, then developed high-throughput molecular technologies on songbird brains, developed technologies for combined analysis of electrophysiological activity and gene regulation in awake behaving animals, and developed bioinformatic tools that predict causal interactions within and between biological levels of organization. This integrative brain project is fitting for the interdisciplinary approaches taken in the current songbird issue of the Journal of Comparative Physiology A and is expected to be conducive to deciphering how brains generate and perceive complex behaviors.     

合成生物系统的组合优化

合成生物学所面临的一项重要挑战是构建具有全新功能的生物系统.由于生物系统固有的复杂性,仅通过理性设计,通常难以使合成基因线路发挥出最优的功能.组合工程的兴起和发展为获得组合优化性状提供了有利条件,并大大促进了具有全新功能的生物系统的构建.文中主要从单个元件的微调、代谢通路的优化以及基因组范围内靶点的识别和组合修饰三个方面入手,总结和评述了近些年表现突出的合成生物系统的组合优化方法.     

Interactome networks: the state of the science

A report on the joint Cold Spring Harbor/Wellcome Trust Meeting 'Interactome Networks', Hinxton, UK, 31 August-4 September 2005.

High-throughput analyses are identifying the DNA, RNA, proteins and metabolites within a biological system with increasing accuracy and speed. As a result, we now have a relatively detailed understanding of the components that make up the dynamic and temporal characteristics of the cell. In most cases, however, we know very little about how the individual components work together to carry out specific biological functions. To get over this hurdle, it will be necessary to map how individual biomolecules interact with one another within a larger network of molecular interactions (the so-called 'interactome') in the cell as a whole.

Mapping this network is the shared goal of an increasing number of researchers from the UK, Europe, US, and Japan, who gathered at the first annual Cold Spring Harbor/Wellcome Trust meeting on interactome networks in Hinxton. This meeting provided an opportunity to review the recent experimental and computational advances that have been applied to uncover biomolecular interactions. Here we report a few of the key advances in the areas of new interaction-mapping techniques, new experimental reagents and resources, and new computational tools for understanding interaction networks that were presented at the meeting.

     

Identification of potential targets in biological signalling systems through network perturbation analysis

The network-based representation and analysis of biological systems contributes to a greater understanding of their structures and functions at different levels of complexity. These techniques can also be used to identify potential novel therapeutic targets based on the characterisation of vulnerable or highly influential network components. There is a need to investigate methods for estimating the impact of molecular perturbations. The prediction of high-impact or critical targets can aid in the identification of novel strategies for controlling the level of activation of specific, therapeutically relevant genes or proteins. Here, we report a new computational strategy for the analysis of the vulnerability of cellular signalling networks based on the quantitative assessment of the impact of large-scale, dynamic perturbations. To show the usefulness of this methodology, two complex signalling networks were analysed: the caspase-3 and the adenosine-regulated calcium signalling systems. This allowed us to estimate and rank the perturbation impact of the components defining these networks. Testable hypotheses about how these targets could modify the dynamic operation of the systems are provided. In the case of the caspase-3 system, the predictions and rankings were in line with results obtained from previous experimental validations of computational predictions generated by a relatively more computationally complex technique. In the case of the adenosine-regulated calcium system, we offer new testable predictions on the potential effect of different targets on the control of calcium flux. Unlike previous methods, the proposed approach provides perturbation-specific scores for each network component. The proposed perturbation assessment methodology may be applied to other systems to gain a deeper understanding of their dynamic operation and to assist the discovery of new therapeutic targets and strategies.     

A novel knowledge-driven systems biology approach for phenotype prediction upon genetic intervention

Deciphering the biological networks underlying complex phenotypic traits, e.g., human disease is undoubtedly crucial to understand the underlying molecular mechanisms and to develop effective therapeutics. Due to the network complexity and the relatively small number of available experiments, data-driven modeling is a great challenge for deducing the functions of genes/proteins in the network and in phenotype formation. We propose a novel knowledge-driven systems biology method that utilizes qualitative knowledge to construct a Dynamic Bayesian network (DBN) to represent the biological network underlying a specific phenotype. Edges in this network depict physical interactions between genes and/or proteins. A qualitative knowledge model first translates typical molecular interactions into constraints when resolving the DBN structure and parameters. Therefore, the uncertainty of the network is restricted to a subset of models which are consistent with the qualitative knowledge. All models satisfying the constraints are considered as candidates for the underlying network. These consistent models are used to perform quantitative inference. By in silico inference, we can predict phenotypic traits upon genetic interventions and perturbing in the network. We applied our method to analyze the puzzling mechanism of breast cancer cell proliferation network and we accurately predicted cancer cell growth rate upon manipulating (anti)cancerous marker genes/proteins.     

Competing to coordinate cell fate decisions: the MST2-Raf-1 signaling device

How do biochemical signaling pathways generate biological specificity? This question is fundamental to modern biology, and its enigma has been accentuated by the discovery that most proteins in signaling networks serve multifunctional roles. An answer to this question may lie in analyzing network properties rather than individual traits of proteins in order to elucidate design principles of biochemical networks that enable biological decision-making. We discuss how this is achieved in the MST2/Hippo-Raf-1 signaling network with the help of mathematical modeling and model-based analysis, which showed that competing protein interactions with affinities controlled by dynamic protein modifications can function as Boolean computing devices that determine cell fate decisions. In addition, we discuss areas of interest for future research and highlight how systems approaches would be of benefit.     

Dupuytren's: a systems biology disease

Dupuytren's disease (DD) is an ill-defined fibroproliferative disorder of the palm of the hands leading to digital contracture. DD commonly occurs in individuals of northern European extraction. Cellular components and processes associated with DD pathogenesis include altered gene and protein expression of cytokines, growth factors, adhesion molecules, and extracellular matrix components. Histology has shown increased but varying levels of particular types of collagen, myofibroblasts and myoglobin proteins in DD tissue. Free radicals and localised ischaemia have been suggested to trigger the proliferation of DD tissue. Although the existing available biological information on DD may contain potentially valuable (though largely uninterpreted) information, the precise aetiology of DD remains unknown. Systems biology combines mechanistic modelling with quantitative experimentation in studies of networks and better understanding of the interaction of multiple components in disease processes. Adopting systems biology may be the ideal approach for future research in order to improve understanding of complex diseases of multifactorial origin. In this review, we propose that DD is a disease of several networks rather than of a single gene, and show that this accounts for the experimental observations obtained to date from a variety of sources. We outline how DD may be investigated more effectively by employing a systems biology approach that considers the disease network as a whole rather than focusing on any specific single molecule.     

Function and Compensatory Mechanisms Among the Components of the Chloroplastic Redox Network

Life on earth depends on the presence of photoautotrophic organisms that are able to input carbon into the ecosystems through the process of photosynthesis which, with a few specialized exceptions, takes place within the chloroplast. This organelle contains the most complex redox system in plants being composed of numerous players including thiol reductases, peroxidases, and glutathione-related enzymes. It seems likely that these proteins act together to adjust redox metabolism enabling plants to grow efficiently under both normal and stressed conditions. However, our knowledge concerning how these proteins interact and if they can compensate one another is relatively limited. This is in part due to the failure of considering these components from a systemic perspective. Here, we provide a systemic view of the chloroplastic-redox network highlighting how it operates and how its components co-operate to maintain efficient chloroplastic function. We further explore the cross-talk between chloroplastic-redox metabolism and that of other subcellular compartments. Given the complexity of plant redox metabolism and the compensatory role played by different redox systems, we argue that a unique possibility to understand this system is afforded by systems biology approaches and by characterizing mutants for multiple genes. Taking this into account, we highlight how gene co-expression and protein–protein network analyses coupled with different reverse genetic strategies could be used to reveal the function, potential redundancies, and complementarities among the components of the chloroplastic redox network.     

An integer optimization algorithm for robust identification of non-linear gene regulatory networks

ABSTRACT: BACKGROUND: Reverse engineering gene networks and identifying regulatory interactions are integral to understanding cellular decision making processes. Advancement in high throughput experimental techniques has initiated innovative data driven analysis of gene regulatory networks. However, inherent noise associated with biological systems requires numerous experimental replicates for reliable conclusions. Furthermore, evidence of robust algorithms directly exploiting basic biological traits are few. Such algorithms are expected to be efficient in their performance and robust in their prediction. RESULTS: We have developed a network identification algorithm to accurately infer both the topology and strength of regulatory interactions from time series gene expression data in the presence of significant experimental noise and non-linear behavior. In this novel formulism, we have addressed data variability in biological systems by integrating network identification with the bootstrap resampling technique, hence predicting robust interactions from limited experimental replicates subjected to noise. Furthermore, we have incorporated non-linearity in gene dynamics using the S-system formulation. The basic network identification formulation exploits the trait of sparsity of biological interactions. Towards that, the identification algorithm is formulated as an integer-programming problem by introducing binary variables for each network component. The objective function is targeted to minimize the network connections subjected to the constraint of maximal agreement between the experimental and predicted gene dynamics. The developed algorithm is validated using both in-silico and experimental data-sets. These studies show that the algorithm can accurately predict the topology and connection strength of the in silico networks, as quantified by high precision and recall, and small discrepancy between the actual and predicted kinetic parameters. Furthermore, in both the in silico and experimental case studies, the predicted gene expression profiles are in very close agreement with the dynamics of the input data. CONCLUSIONS: Our integer programming algorithm effectively utilizes bootstrapping to identify robust gene regulatory networks from noisy, non-linear time-series gene expression data. With significant noise and non-linearities being inherent to biological systems, the present formulism, with the incorporation of network sparsity, is extremely relevant to gene regulatory networks, and while the formulation has been validated against in silico and E. Coli data, it can be applied to any biological system.     

Efficient modeling, simulation and coarse-graining of biological complexity with NFsim

Managing the overwhelming numbers of molecular states and interactions is a fundamental obstacle to building predictive models of biological systems. Here we introduce the Network-Free Stochastic Simulator (NFsim), a general-purpose modeling platform that overcomes the combinatorial nature of molecular interactions. Unlike standard simulators that represent molecular species as variables in equations, NFsim uses a biologically intuitive representation: objects with binding and modification sites acted on by reaction rules. During simulations, rules operate directly on molecular objects to produce exact stochastic results with performance that scales independently of the reaction network size. Reaction rates can be defined as arbitrary functions of molecular states to provide powerful coarse-graining capabilities, for example to merge Boolean and kinetic representations of biological networks. NFsim enables researchers to simulate many biological systems that were previously inaccessible to general-purpose software, as we illustrate with models of immune system signaling, microbial signaling, cytoskeletal assembly and oscillating gene expression.     

Spatio-logical processes in intracellular signalling

Classical models of intracellular signalling describe how small changes in a cell's external environment can bring about major changes in cellular activity. Recent findings from experimental biology indicate that many intracellular signalling systems show a high level of spatial organisation. This permits the modification, by protein kinase or protein phosphatase action, of specific subsets of intracellular proteins - an attribute that is not addressed in classical signalling models. Here we use ideas and concepts from computer science to describe the information processing nature of intracellular signalling pathways and the impact of spatial heterogeneity of their components (e.g. protein kinases and protein phosphatases) on signalling activity. We argue that it is useful to view the signalling ecology as a vast parallel distributed processing network of agents operating in heterogeneous microenvironments, and we conclude with an overview of the mathematical and semantic methodologies that might help clarify this analogy between biological and computational systems.     

Bridging from molecular simulation to biochemical networks

How can we make the connection between the three-dimensional structures of individual proteins and understanding how complex biological systems involving many proteins work? The modelling and simulation of protein structures can help to answer this question for systems ranging from multimacromolecular complexes to organelles and cells. On one hand, multiscale modelling and simulation techniques are advancing to permit the spatial and temporal properties of large systems to be simulated using atomic-detail structures. On the other hand, the estimation of kinetic parameters for the mathematical modelling of biochemical pathways using protein structure information provides a basis for iterative manipulation of biochemical pathways guided by protein structure. Recent advances include the structural modelling of protein complexes on the genomic level, novel coarse-graining strategies to increase the size of the system and the time span that can be simulated, and comparative molecular field analyses to estimate enzyme kinetic parameters.     

Toward a systems level view of the ECM and related proteins: a framework for the systematic definition and analysis of biological systems

Advances in high throughput 'omic technologies are starting to provide unprecedented insights into how components of biological systems are organized and interact. Key to exploiting these datasets is the definition of the components that comprise the system of interest. Although a variety of knowledge bases exist that capture such information, a major challenge is determining how these resources may be best utilized. Here we present a systematic curation strategy to define a systems-level view of the human extracellular matrix (ECM)--a three-dimensional meshwork of proteins and polysaccharides that impart structure and mechanical stability to tissues. Employing our curation strategy we define a set of 357 proteins that represent core components of the ECM, together with an additional 524 genes that mediate related functional roles, and construct a map of their physical interactions. Topological properties help identify modules of functionally related proteins, including those involved in cell adhesion, bone formation and blood clotting. Because of its major role in cell adhesion, proliferation and morphogenesis, defects in the ECM have been implicated in cancer, atherosclerosis, asthma, fibrosis, and arthritis. We use MeSH annotations to identify modules enriched for specific disease terms that aid to strengthen existing as well as predict novel gene-disease associations. Mapping expression and conservation data onto the network reveal modules evolved in parallel to convey tissue-specific functionality on otherwise broadly expressed units. In addition to demonstrating an effective workflow for defining biological systems, this study crystallizes our current knowledge surrounding the organization of the ECM.     

From computational modelling of the intrinsic apoptosis pathway to a systems-based analysis of chemotherapy resistance: achievements,perspectives and challenges in systems medicine

Our understanding of the mitochondrial or intrinsic apoptosis pathway and its role in chemotherapy resistance has increased significantly in recent years by a combination of experimental studies and mathematical modelling. This combined approach enhanced the quantitative and kinetic understanding of apoptosis signal transduction, but also provided new insights that systems-emanating functions (i.e., functions that cannot be attributed to individual network components but that are instead established by multi-component interplay) are crucial determinants of cell fate decisions. Among these features are molecular thresholds, cooperative protein functions, feedback loops and functional redundancies that provide systems robustness, and signalling topologies that allow ultrasensitivity or switch-like responses. The successful development of kinetic systems models that recapitulate biological signal transduction observed in living cells have now led to the first translational studies, which have exploited and validated such models in a clinical context. Bottom-up strategies that use pathway models in combination with higher-level modelling at the tissue, organ and whole body-level therefore carry great potential to eventually deliver a new generation of systems-based diagnostic tools that may contribute to the development of personalised and predictive medicine approaches. Here we review major achievements in the systems biology of intrinsic apoptosis signalling, discuss challenges for further model development, perspectives for higher-level integration of apoptosis models and finally discuss requirements for the development of systems medical solutions in the coming years.