Systems biology by the rules: hybrid intelligent systems for pathway modeling and discovery

Fuzzy logic based approaches for gene regulatory network inference

The rapid advancements in high-throughput techniques have fueled large-scale production of biological data at very affordable costs. Some of these techniques are microarrays and next-generation sequencing that provide genome level insight of living cells. As a result, the size of most of the biological databases, such as NCBI-GEO, NCBI-SRA, etc., is growing exponentially. These biological data are analyzed using various computational techniques for knowledge discovery - which is also one of the objectives of bioinformatics research. Gene regulatory network (GRN) is a gene-gene interaction network which plays a pivotal role in understanding gene regulation processes and disease mechanism at the molecular level. From last couple of decades, researchers are interested in developing computational algorithms for GRN inference (GRNI) from high-throughput experimental data. Several computational approaches have been proposed for inferring GRN from gene expression data including statistical techniques (correlation coefficient), information theory (mutual information), regression-based approaches, probabilistic approaches (Bayesian networks, naïve byes), artificial neural networks and fuzzy logic. The fuzzy logic, along with its hybridization with other intelligent approaches, is a well-studied technique in GRNI due to its several advantages. In this paper, we present a consolidated review on fuzzy logic and its hybrid approaches developed during last two decades for GRNI.

Computational Systems Biology in Cancer: Modeling Methods and Applications

In recent years it has become clear that carcinogenesis is a complex process, both at the molecular and cellular levels. Understanding the origins, growth and spread of cancer, therefore requires an integrated or system-wide approach. Computational systems biology is an emerging sub-discipline in systems biology that utilizes the wealth of data from genomic, proteomic and metabolomic studies to build computer simulations of intra and intercellular processes. Several useful descriptive and predictive models of the origin, growth and spread of cancers have been developed in an effort to better understand the disease and potential therapeutic approaches. In this review we describe and assess the practical and theoretical underpinnings of commonly-used modeling approaches, including ordinary and partial differential equations, petri nets, cellular automata, agent based models and hybrid systems. A number of computer-based formalisms have been implemented to improve the accessibility of the various approaches to researchers whose primary interest lies outside of model development. We discuss several of these and describe how they have led to novel insights into tumor genesis, growth, apoptosis, vascularization and therapy.

Mathematical modelling and a systems science approach to describe the role of cytokines in the evolution of severe dengue

Background

Dengue causes considerable morbidity and mortality in Sri Lanka. Inflammatory mediators such as cytokines, contribute to its evolution from an asymptotic infection to severe forms of dengue. The majority of previous studies have analysed the association of individual cytokines with clinical disease severity. In contrast, we view evolution to Dengue Haemorrhagic Fever as the behaviour of a complex dynamic system. We therefore, analyse the combined effect of multiple cytokines that interact dynamically with each other in order to generate a mathematical model to predict occurrence of Dengue Haemorrhagic Fever. We expect this to have predictive value in detecting severe cases and improve outcomes. Platelet activating factor (PAF), Sphingosine 1- Phosphate (S1P), IL-1β, TNFα and IL-10 are used as the parameters for the model. Hierarchical clustering is used to detect factors that correlated with each other. Their interactions are mapped using Fuzzy Logic mechanisms with the combination of modified Hamacher and OWA operators. Trapezoidal membership functions are developed for each of the cytokine parameters and the degree of unfavourability to attain Dengue Haemorrhagic Fever is measured.

Results

The accuracy of this model in predicting severity level of dengue is 71.43% at 96 h from the onset of illness, 85.00% at 108 h and 76.92% at 120 h. A region of ambiguity is detected in the model for the value range 0.36 to 0.51. Sensitivity analysis indicates that this is a robust mathematical model.

Conclusions

The results show a robust mathematical model that explains the evolution from dengue to its serious forms in individual patients with high accuracy. However, this model would have to be further improved by including additional parameters and should be validated on other data sets.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-017-0415-3) contains supplementary material, which is available to authorized users.

Computational modeling in melanoma for novel drug discovery

INTRODUCTION

There is a growing body of evidence highlighting the applications of computational modeling in the field of biomedicine. It has recently been applied to the in silico analysis of cancer dynamics. In the era of precision medicine, this analysis may allow the discovery of new molecular targets useful for the design of novel therapies and for overcoming resistance to anticancer drugs. According to its molecular behavior, melanoma represents an interesting tumor model in which computational modeling can be applied. Melanoma is an aggressive tumor of the skin with a poor prognosis for patients with advanced disease as it is resistant to current therapeutic approaches.

AREAS COVERED

This review discusses the basics of computational modeling in melanoma drug discovery and development. Discussion includes the in silico discovery of novel molecular drug targets, the optimization of immunotherapies and personalized medicine trials.

EXPERT OPINION

Mathematical and computational models are gradually being used to help understand biomedical data produced by high-throughput analysis. The use of advanced computer models allowing the simulation of complex biological processes provides hypotheses and supports experimental design. The research in fighting aggressive cancers, such as melanoma, is making great strides. Computational models represent the key component to complement these efforts. Due to the combinatorial complexity of new drug discovery, a systematic approach based only on experimentation is not possible. Computational and mathematical models are necessary for bringing cancer drug discovery into the era of omics, big data and personalized medicine.

DYNAMICS OF GLI REGULATION AND A STRATEGY TO CONTROL CANCEROUS SITUATION: HEDGEHOG SIGNALING PATHWAY REVISITED

The hedgehog signaling cascade generates highly diverse, fine-tuned responses in response to the external stimulus by the sonic hedgehog (SHH) protein. This is required for the flawless functioning of the cell, its development, survival and proliferation; maintained through production of Glioma protein (GLI) and transcriptional activation of its target genes. Any change in the behavior of GLI response by ectopic expression of SHH or mutations in the core pathway components may cause serious consequences in the cell fate through rapid, uncontrolled and elevated production of GLI. Here, we present a simple but extensive computational model that considers the detailed reaction mechanisms involved in the hedgehog signal transduction and provides a detailed insight into regulation of GLI. For the first time, by explicit involvement of suppressor of fused (SUFU) and Hedgehog interacting protein (HHIP) reaction kinetics in the model, we try to demonstrate the vital importance of HHIP and SUFU in maintaining the graded response of GLI in response to SHH. By performing parameter variations, we capture the conversion of a graded response of GLI to an ultrasensitive switch under SUFU-deficient conditions that might predispose abnormal embryonic development and the irreversible switching response of GLI that corresponds to signal-independent pathway activation observed in cancers.

Network reconstruction based on proteomic data and prior knowledge of protein connectivity using graph theory

Modeling of signal transduction pathways is instrumental for understanding cells’ function. People have been tackling modeling of signaling pathways in order to accurately represent the signaling events inside cells’ biochemical microenvironment in a way meaningful for scientists in a biological field. In this article, we propose a method to interrogate such pathways in order to produce cell-specific signaling models. We integrate available prior knowledge of protein connectivity, in a form of a Prior Knowledge Network (PKN) with phosphoproteomic data to construct predictive models of the protein connectivity of the interrogated cell type. Several computational methodologies focusing on pathways’ logic modeling using optimization formulations or machine learning algorithms have been published on this front over the past few years. Here, we introduce a light and fast approach that uses a breadth-first traversal of the graph to identify the shortest pathways and score proteins in the PKN, fitting the dependencies extracted from the experimental design. The pathways are then combined through a heuristic formulation to produce a final topology handling inconsistencies between the PKN and the experimental scenarios. Our results show that the algorithm we developed is efficient and accurate for the construction of medium and large scale signaling networks. We demonstrate the applicability of the proposed approach by interrogating a manually curated interaction graph model of EGF/TNFA stimulation against made up experimental data. To avoid the possibility of erroneous predictions, we performed a cross-validation analysis. Finally, we validate that the introduced approach generates predictive topologies, comparable to the ILP formulation. Overall, an efficient approach based on graph theory is presented herein to interrogate protein–protein interaction networks and to provide meaningful biological insights.

Data-derived modeling characterizes plasticity of MAPK signaling in melanoma

The majority of melanomas have been shown to harbor somatic mutations in the RAS-RAF-MEK-MAPK and PI3K-AKT pathways, which play a major role in regulation of proliferation and survival. The prevalence of these mutations makes these kinase signal transduction pathways an attractive target for cancer therapy. However, tumors have generally shown adaptive resistance to treatment. This adaptation is achieved in melanoma through its ability to undergo neovascularization, migration and rearrangement of signaling pathways. To understand the dynamic, nonlinear behavior of signaling pathways in cancer, several computational modeling approaches have been suggested. Most of those models require that the pathway topology remains constant over the entire observation period. However, changes in topology might underlie adaptive behavior to drug treatment. To study signaling rearrangements, here we present a new approach based on Fuzzy Logic (FL) that predicts changes in network architecture over time. This adaptive modeling approach was used to investigate pathway dynamics in a newly acquired experimental dataset describing total and phosphorylated protein signaling over four days in A375 melanoma cell line exposed to different kinase inhibitors. First, a generalized strategy was established to implement a parameter-reduced FL model encoding non-linear activity of a signaling network in response to perturbation. Next, a literature-based topology was generated and parameters of the FL model were derived from the full experimental dataset. Subsequently, the temporal evolution of model performance was evaluated by leaving time-defined data points out of training. Emerging discrepancies between model predictions and experimental data at specific time points allowed the characterization of potential network rearrangement. We demonstrate that this adaptive FL modeling approach helps to enhance our mechanistic understanding of the molecular plasticity of melanoma.

Signal transduction pathways can be described as static routes, transmitting extrinsic signals to the nucleus to induce a transcriptional response. In contrast to this reductionist view, the emerging paradigm is that signaling networks undergo dynamic crosstalk, both in disease and physiological conditions. To understand complex pathway behavior, it is necessary to develop methods to identify pathway interactions that are active as a consequence of stimuli and, importantly, to describe their evolution in time. To that end, we developed a method relying on prior knowledge networks in order to predict signaling crosstalk evolution, in response to perturbation and over time. The challenge we addressed was to establish a method dependent on information related to the topology of reported interactions, and not their mechanistic characteristics, and at the same time complex enough to reproduce the behavior of the signaling intermediates. The work presented here demonstrates that such an approach can be used to predict mechanisms that melanoma uses to rearrange its signaling and maintain its abnormal proliferation upon treatment.

Seizure detection, seizure prediction, and closed-loop warning systems in epilepsy

Nearly one-third of patients with epilepsy continue to have seizures despite optimal medication management. Systems employed to detect seizures may have the potential to improve outcomes in these patients by allowing more tailored therapies and might, additionally, have a role in accident and SUDEP prevention. Automated seizure detection and prediction require algorithms which employ feature computation and subsequent classification. Over the last few decades, methods have been developed to detect seizures utilizing scalp and intracranial EEG, electrocardiography, accelerometry and motion sensors, electrodermal activity, and audio/video captures. To date, it is unclear which combination of detection technologies yields the best results, and approaches may ultimately need to be individualized. This review presents an overview of seizure detection and related prediction methods and discusses their potential uses in closed-loop warning systems in epilepsy.

Introductory review of computational cell cycle modeling

Recent advances in the modeling of the cell cycle through computer simulation demonstrate the power of systems biology. By definition, systems biology has the goal to connect a parts list, prioritized through experimental observation or high-throughput screens, by the topology of interactions defining intracellular networks to predict system function. Computer modeling of biological systems is often compared to a process of reverse engineering. Indeed, designed or engineered technical systems share many systems-level properties with biological systems; thus studying biological systems within an engineering framework has proven successful. Here we review some aspects of this process as it pertains to cell cycle modeling.

Detecting and removing inconsistencies between experimental data and signaling network topologies using integer linear programming on interaction graphs

Cross-referencing experimental data with our current knowledge of signaling network topologies is one central goal of mathematical modeling of cellular signal transduction networks. We present a new methodology for data-driven interrogation and training of signaling networks. While most published methods for signaling network inference operate on Bayesian, Boolean, or ODE models, our approach uses integer linear programming (ILP) on interaction graphs to encode constraints on the qualitative behavior of the nodes. These constraints are posed by the network topology and their formulation as ILP allows us to predict the possible qualitative changes (up, down, no effect) of the activation levels of the nodes for a given stimulus. We provide four basic operations to detect and remove inconsistencies between measurements and predicted behavior: (i) find a topology-consistent explanation for responses of signaling nodes measured in a stimulus-response experiment (if none exists, find the closest explanation); (ii) determine a minimal set of nodes that need to be corrected to make an inconsistent scenario consistent; (iii) determine the optimal subgraph of the given network topology which can best reflect measurements from a set of experimental scenarios; (iv) find possibly missing edges that would improve the consistency of the graph with respect to a set of experimental scenarios the most. We demonstrate the applicability of the proposed approach by interrogating a manually curated interaction graph model of EGFR/ErbB signaling against a library of high-throughput phosphoproteomic data measured in primary hepatocytes. Our methods detect interactions that are likely to be inactive in hepatocytes and provide suggestions for new interactions that, if included, would significantly improve the goodness of fit. Our framework is highly flexible and the underlying model requires only easily accessible biological knowledge. All related algorithms were implemented in a freely available toolbox SigNetTrainer making it an appealing approach for various applications.

Cellular signal transduction is orchestrated by communication networks of signaling proteins commonly depicted on signaling pathway maps. However, each cell type may have distinct variants of signaling pathways, and wiring diagrams are often altered in disease states. The identification of truly active signaling topologies based on experimental data is therefore one key challenge in systems biology of cellular signaling. We present a new framework for training signaling networks based on interaction graphs (IG). In contrast to complex modeling formalisms, IG capture merely the known positive and negative edges between the components. This basic information, however, already sets hard constraints on the possible qualitative behaviors of the nodes when perturbing the network. Our approach uses Integer Linear Programming to encode these constraints and to predict the possible changes (down, neutral, up) of the activation levels of the involved players for a given experiment. Based on this formulation we developed several algorithms for detecting and removing inconsistencies between measurements and network topology. Demonstrated by EGFR/ErbB signaling in hepatocytes, our approach delivers direct conclusions on edges that are likely inactive or missing relative to canonical pathway maps. Such information drives the further elucidation of signaling network topologies under normal and pathological phenotypes.

A role of TGFß1 dependent 14-3-3σ phosphorylation at Ser69 and Ser74 in the regulation of gene transcription, stemness and radioresistance

Transforming growth factor-β (TGFβ) is a potent regulator of tumorigenesis, although mechanisms defining its tumor suppressing and tumor promoting activities are not understood. Here we describe phosphoproteome profiling of TGFβ signaling in mammary epithelial cells, and show that 60 identified TGFβ-regulated phosphoproteins form a network with scale-free characteristics. The network highlighted interactions, which may distribute signaling inputs to regulation of cell proliferation, metabolism, differentiation and cell organization. In this report, we identified two novel and TGFβ-dependent phosphorylation sites of 14-3-3σ, i.e. Ser69 and Ser74. We observed that 14-3-3σ phosphorylation is a feed-forward mechanism in TGFβ/Smad3-dependent transcription. TGFβ-dependent 14-3-3σ phosphorylation may provide a scaffold for the formation of the protein complexes which include Smad3 and p53 at the Smad3-specific CAGA element. Furthermore, breast tumor xenograft studies in mice and radiobiological assays showed that phosphorylation of 14-3-3σ at Ser69 and Ser74 is involved in regulation of cancer progenitor population and radioresistance in breast cancer MCF7 cells. Our data suggest that TGFβ-dependent phosphorylation of 14-3-3σ orchestrates a functional interaction of TGFβ/Smad3 with p53, plays a role in the maintenance of cancer stem cells and could provide a new potential target for intervention in breast cancer.

Defining systems biology: a brief overview of the term and field

Here we provide a broad overview of the definition of the term "systems biology" as well as pinpoint specific events in biological research and beyond that are consistently cited to have contributed and led to the current science of in silico systems biology. Since there have been many reviews and historical accounts describing the term, it would be impossible to include all single references. However, we do attempt to provide a consensus vision of how the field has evolved and consequently the terminology that followed it. We also highlight the development and general acceptance, and use, of standards for model representations as being crucial to the continued success of the in silico systems biology field.

Non Linear Programming (NLP) formulation for quantitative modeling of protein signal transduction pathways

Modeling of signal transduction pathways plays a major role in understanding cells' function and predicting cellular response. Mathematical formalisms based on a logic formalism are relatively simple but can describe how signals propagate from one protein to the next and have led to the construction of models that simulate the cells response to environmental or other perturbations. Constrained fuzzy logic was recently introduced to train models to cell specific data to result in quantitative pathway models of the specific cellular behavior. There are two major issues in this pathway optimization: i) excessive CPU time requirements and ii) loosely constrained optimization problem due to lack of data with respect to large signaling pathways. Herein, we address both issues: the former by reformulating the pathway optimization as a regular nonlinear optimization problem; and the latter by enhanced algorithms to pre/post-process the signaling network to remove parts that cannot be identified given the experimental conditions. As a case study, we tackle the construction of cell type specific pathways in normal and transformed hepatocytes using medium and large-scale functional phosphoproteomic datasets. The proposed Non Linear Programming (NLP) formulation allows for fast optimization of signaling topologies by combining the versatile nature of logic modeling with state of the art optimization algorithms.

Inference of biological pathway from gene expression profiles by time delay boolean networks

One great challenge of genomic research is to efficiently and accurately identify complex gene regulatory networks. The development of high-throughput technologies provides numerous experimental data such as DNA sequences, protein sequence, and RNA expression profiles makes it possible to study interactions and regulations among genes or other substance in an organism. However, it is crucial to make inference of genetic regulatory networks from gene expression profiles and protein interaction data for systems biology. This study will develop a new approach to reconstruct time delay boolean networks as a tool for exploring biological pathways. In the inference strategy, we will compare all pairs of input genes in those basic relationships by their corresponding p-scores for every output gene. Then, we will combine those consistent relationships to reveal the most probable relationship and reconstruct the genetic network. Specifically, we will prove that O(log n) state transition pairs are sufficient and necessary to reconstruct the time delay boolean network of n nodes with high accuracy if the number of input genes to each gene is bounded. We also have implemented this method on simulated and empirical yeast gene expression data sets. The test results show that this proposed method is extensible for realistic networks.

Multi-parametric analysis and modeling of relationships between mitochondrial morphology and apoptosis

Mitochondria exist as a network of interconnected organelles undergoing constant fission and fusion. Current approaches to study mitochondrial morphology are limited by low data sampling coupled with manual identification and classification of complex morphological phenotypes. Here we propose an integrated mechanistic and data-driven modeling approach to analyze heterogeneous, quantified datasets and infer relations between mitochondrial morphology and apoptotic events. We initially performed high-content, multi-parametric measurements of mitochondrial morphological, apoptotic, and energetic states by high-resolution imaging of human breast carcinoma MCF-7 cells. Subsequently, decision tree-based analysis was used to automatically classify networked, fragmented, and swollen mitochondrial subpopulations, at the single-cell level and within cell populations. Our results revealed subtle but significant differences in morphology class distributions in response to various apoptotic stimuli. Furthermore, key mitochondrial functional parameters including mitochondrial membrane potential and Bax activation, were measured under matched conditions. Data-driven fuzzy logic modeling was used to explore the non-linear relationships between mitochondrial morphology and apoptotic signaling, combining morphological and functional data as a single model. Modeling results are in accordance with previous studies, where Bax regulates mitochondrial fragmentation, and mitochondrial morphology influences mitochondrial membrane potential. In summary, we established and validated a platform for mitochondrial morphological and functional analysis that can be readily extended with additional datasets. We further discuss the benefits of a flexible systematic approach for elucidating specific and general relationships between mitochondrial morphology and apoptosis.

A hybrid model of mammalian cell cycle regulation

The timing of DNA synthesis, mitosis and cell division is regulated by a complex network of biochemical reactions that control the activities of a family of cyclin-dependent kinases. The temporal dynamics of this reaction network is typically modeled by nonlinear differential equations describing the rates of the component reactions. This approach provides exquisite details about molecular regulatory processes but is hampered by the need to estimate realistic values for the many kinetic constants that determine the reaction rates. It is difficult to estimate these kinetic constants from available experimental data. To avoid this problem, modelers often resort to ‘qualitative’ modeling strategies, such as Boolean switching networks, but these models describe only the coarsest features of cell cycle regulation. In this paper we describe a hybrid approach that combines the best features of continuous differential equations and discrete Boolean networks. Cyclin abundances are tracked by piecewise linear differential equations for cyclin synthesis and degradation. Cyclin synthesis is regulated by transcription factors whose activities are represented by discrete variables (0 or 1) and likewise for the activities of the ubiquitin-ligating enzyme complexes that govern cyclin degradation. The discrete variables change according to a predetermined sequence, with the times between transitions determined in part by cyclin accumulation and degradation and as well by exponentially distributed random variables. The model is evaluated in terms of flow cytometry measurements of cyclin proteins in asynchronous populations of human cell lines. The few kinetic constants in the model are easily estimated from the experimental data. Using this hybrid approach, modelers can quickly create quantitatively accurate, computational models of protein regulatory networks in cells.

The physiological behaviors of cells (growth and division, differentiation, movement, death, etc.) are controlled by complex networks of interacting genes and proteins, and a fundamental goal of computational cell biology is to develop dynamical models of these regulatory networks that are realistic, accurate and predictive. Historically, these models have divided along two basic lines: deterministic or stochastic, and continuous or discrete; with scattered efforts to develop hybrid approaches that bridge these divides. Using the cell cycle control system in eukaryotes as an example, we propose a hybrid approach that combines a continuous representation of slowly changing protein concentrations with a discrete representation of components that switch rapidly between ‘on’ and ‘off’ states, and that combines the deterministic causality of network interactions with the stochastic uncertainty of random events. The hybrid approach can be easily tailored to the available knowledge of control systems, and it provides both qualitative and quantitative results that can be compared to experimental data to test the accuracy and predictive power of the model.

Computational systems biology of aging

Computational systems biology is expected to make major contributions to unravel the complex molecular mechanisms underlying the progression of aging in cells, tissues, and organisms. The development of computational approaches is, however, challenged by a wide spectrum of aging mechanisms participating on different levels of biological organization. The tight connectivity between the molecular constituents, functions, and cell states requires frameworks and strategies that extend beyond current practice to model, simulate, and predict the progression of aging and the emerging aging phenotypes. We provide a general overview of the specific computational tasks and opportunities in aging research, and discuss some illustrative systems level concepts in more detail. One example provided here is the assembly of a conceptual whole cell model that considers the temporal dynamics of the aging process grounded on molecular mechanisms. Another application is the assembly of interactomes, such as protein networks that allow us to analyze changes in network topology and interaction of proteins that have been implicated in aging with other cellular constituents and processes. We introduce the necessary key steps to build these applications and discuss their merits and future extensions for aging research. WIREs Syst Biol Med 2011 3 414-428 DOI: 10.1002/wsbm.126

The role of noise and positive feedback in the onset of autosomal dominant diseases

Background

Autosomal dominant (AD) diseases result when a single mutant or non-functioning gene is present on an autosomal chromosome. These diseases often do not emerge at birth. There are presently two prevailing theories explaining the expression of AD diseases. One explanation originates from the Knudson two-hit theory of hereditary cancers, where loss of heterozygosity or occurrence of somatic mutations impairs the function of the wild-type copy. While these somatic second hits may be sufficient for stable disease states, it is often difficult to determine if their occurrence necessarily marks the initiation of disease progression. A more direct consequence of a heterozygous genetic background is haploinsufficiency, referring to a lack of sufficient gene function due to reduced wild-type gene copy number; however, haploinsufficiency can involve a variety of additional mechanisms, such as noise in gene expression or protein levels, injury and second hit mutations in other genes. In this study, we explore the possible contribution to the onset of autosomal dominant diseases from intrinsic factors, such as those determined by the structure of the molecular networks governing normal cellular physiology.

Results

First, simple models of single gene insufficiency using the positive feedback loops that may be derived from a three-component network were studied by computer simulation using Bionet software. The network structure is shown to affect the dynamics considerably; some networks are relatively stable even when large stochastic variations in are present, while others exhibit switch-like dynamics. In the latter cases, once the network switches over to the disease state it remains in that state permanently. Model pathways for two autosomal dominant diseases, AD polycystic kidney disease and mature onset diabetes of youth (MODY) were simulated and the results are compared to known disease characteristics.

Conclusions

By identifying the intrinsic mechanisms involved in the onset of AD diseases, it may be possible to better assess risk factors as well as lead to potential new drug targets. To illustrate the applicability of this study of pathway dynamics, we simulated the primary pathways involved in two autosomal dominant diseases, Polycystic Kidney Disease (PKD) and mature onset diabetes of youth (MODY). Simulations demonstrate that some of the primary disease characteristics are consistent with the positive feedback - stochastic variation theory presented here. This has implications for new drug targets to control these diseases by blocking the positive feedback loop in the relevant pathways.

Rule-based cell systems model of aging using feedback loop motifs mediated by stress responses

Investigating the complex systems dynamics of the aging process requires integration of a broad range of cellular processes describing damage and functional decline co-existing with adaptive and protective regulatory mechanisms. We evolve an integrated generic cell network to represent the connectivity of key cellular mechanisms structured into positive and negative feedback loop motifs centrally important for aging. The conceptual network is casted into a fuzzy-logic, hybrid-intelligent framework based on interaction rules assembled from a priori knowledge. Based upon a classical homeostatic representation of cellular energy metabolism, we first demonstrate how positive-feedback loops accelerate damage and decline consistent with a vicious cycle. This model is iteratively extended towards an adaptive response model by incorporating protective negative-feedback loop circuits. Time-lapse simulations of the adaptive response model uncover how transcriptional and translational changes, mediated by stress sensors NF-κB and mTOR, counteract accumulating damage and dysfunction by modulating mitochondrial respiration, metabolic fluxes, biosynthesis, and autophagy, crucial for cellular survival. The model allows consideration of lifespan optimization scenarios with respect to fitness criteria using a sensitivity analysis. Our work establishes a novel extendable and scalable computational approach capable to connect tractable molecular mechanisms with cellular network dynamics underlying the emerging aging phenotype.

The global process of aging disturbs a broad range of cellular mechanisms in a complex fashion and is not well understood. One important goal of computational approaches in aging is to develop integrated models in terms of a unifying aging theory, predicting progression of aging phenotypes grounded on molecular mechanisms. However, current experimental data incoherently reflects many isolated processes from a large diversity of approaches, biological model systems, and species, which makes such integration a challenging task. In an attempt to close this gap, we iteratively develop a fuzzy-logic cell systems model considering the interplay of damage, metabolism, and signaling by positive and negative feedback-loop motifs using relationships drawn from literature data. Because cellular biodynamics may be considered a complex control system, this approach seems particularly suitable. Here, we demonstrate that rule-based fuzzy-logic models provide semi-quantitative predictions that enhance our understanding of complex and interlocked molecular mechanisms and their implications on the aging physiome.

Cancer systems biology: a network modeling perspective

Cancer is now appreciated as not only a highly heterogenous pathology with respect to cell type and tissue origin but also as a disease involving dysregulation of multiple pathways governing fundamental cell processes such as death, proliferation, differentiation and migration. Thus, the activities of molecular networks that execute metabolic or cytoskeletal processes, or regulate these by signal transduction, are altered in a complex manner by diverse genetic mutations in concert with the environmental context. A major challenge therefore is how to develop actionable understanding of this multivariate dysregulation, with respect both to how it arises from diverse genetic mutations and to how it may be ameliorated by prospective treatments. While high-throughput experimental platform technologies ranging from genomic sequencing to transcriptomic, proteomic and metabolomic profiling are now commonly used for molecular-level characterization of tumor cells and surrounding tissues, the resulting data sets defy straightforward intuitive interpretation with respect to potential therapeutic targets or the effects of perturbation. In this review article, we will discuss how significant advances can be obtained by applying computational modeling approaches to elucidate the pathways most critically involved in tumor formation and progression, impact of particular mutations on pathway operation, consequences of altered cell behavior in tissue environments and effects of molecular therapeutics.

Genes regulating molecular and cellular functions in noninfectious nonallergic rhinitis

BACKGROUND

Chronic noninfectious, nonallergic rhinitis (NINAR) is a complex syndrome with a principally unknown pathophysiology. New technology has made it possible to examine differentially expressed genes and according to network theory, genes connected by their function that might have key roles in the disease.

METHODS

Connectivity analysis was used to identify NINAR key genes. mRNA was extracted from nasal biopsies from 12 NINAR patients and 12 healthy volunteers. Microarrays were performed using Affymetrix chips with 54 613 genes. Data were analysed with the Ingenuity Pathway System for organization of genes into annotated biological functions and, thereafter, linking genes into networks due to their connectivity. The regulation of key genes was confirmed with reverse transcription-polymerase chain reaction (RT-PCR).

RESULTS

In all, 43 genes were differentially expressed. The functional analysis showed that these genes were primarily involved in cellular movement, haematological system development and immune response. Merging these functions, 10 genes were found to be shared. Network analysis generated three networks and two of these 'shared genes' in key positions, c-fos and cell division cycle 42 (Cdc42). These genes were upregulated in both the array and the RT-PCR analysis.

CONCLUSION

Ten genes were found to be of pathophysiological interest for NINAR and of these, c-fos and Cdc42 seemed to be of specific interest due to their ability to interact with other genes of interest within this context. Although the role of c-fos and Cdc42 in upper airway inflammation remains unknown, they might be used as potential disease markers.

Fuzzy logic analysis of kinase pathway crosstalk in TNF/EGF/insulin-induced signaling

When modeling cell signaling networks, a balance must be struck between mechanistic detail and ease of interpretation. In this paper we apply a fuzzy logic framework to the analysis of a large, systematic dataset describing the dynamics of cell signaling downstream of TNF, EGF, and insulin receptors in human colon carcinoma cells. Simulations based on fuzzy logic recapitulate most features of the data and generate several predictions involving pathway crosstalk and regulation. We uncover a relationship between MK2 and ERK pathways that might account for the previously identified pro-survival influence of MK2. We also find unexpected inhibition of IKK following EGF treatment, possibly due to down-regulation of autocrine signaling. More generally, fuzzy logic models are flexible, able to incorporate qualitative and noisy data, and powerful enough to produce quantitative predictions and new biological insights about the operation of signaling networks.

Cells use networks of interacting proteins to interpret intra-cellular state and extra-cellular cues and to execute cell-fate decisions. Even when individual proteins are well understood at a molecular level, the dynamics and behavior of networks as a whole are harder to understand. However, deciphering the operation of such networks is key to understanding disease processes and therapeutic opportunities. As a means to study signaling networks, we have modified and applied a fuzzy logic approach originally developed for industrial control. We use fuzzy logic to model the responses of colon cancer cells in culture to combinations of pro-survival and pro-death cytokines, making it possible to interpret quantitative data in the context of abstract information drawn from the literature. Our work establishes that fuzzy logic can be used to understand complex signaling pathways with respect to multi-factorial activity-based protein data and prior knowledge.

Knowledge-based expert systems and a proof-of-concept case study for multiple sequence alignment construction and analysis

The traditional approach to bioinformatics analyses relies on independent task-specific services and applications, using different input and output formats, often idiosyncratic, and frequently not designed to inter-operate. In general, such analyses were performed by experts who manually verified the results obtained at each step in the process. Today, the amount of bioinformatics information continuously being produced means that handling the various applications used to study this information presents a major data management and analysis challenge to researchers. It is now impossible to manually analyse all this information and new approaches are needed that are capable of processing the large-scale heterogeneous data in order to extract the pertinent information. We review the recent use of integrated expert systems aimed at providing more efficient knowledge extraction for bioinformatics research. A general methodology for building knowledge-based expert systems is described, focusing on the unstructured information management architecture, UIMA, which provides facilities for both data and process management. A case study involving a multiple alignment expert system prototype called AlexSys is also presented.

Robust signaling networks of the adipose secretome

Type 2 diabetes is a prototypical complex systems disease that has a strong hereditary component and etiologic links with a sedentary lifestyle, overeating and obesity. Adipose tissue has been shown to be a central driver of type 2 diabetes progression, establishing and maintaining a chronic state of low-level inflammation. The number and diversity of identified endocrine factors from adipose tissue (adipokines) is growing rapidly. Here, I argue that a systems biology approach to understanding the robust multi-level signaling networks established by the adipose secretome will be crucial for developing efficient type 2 diabetes treatment. Recent advances in whole-genome association studies, global molecular profiling and quantitative modeling are currently fueling the emergence of this novel research strategy.

Development of therapeutic approaches to 'triple negative' phenotype breast cancer

BACKGROUND

Triple negative phenotype (TNP) breast cancers are characterised by the lack of expression of oestrogen and progesterone receptors and of human EGF receptor 2 (HER2) overexpression/amplification. This subgroup of cancers has an aggressive clinical behaviour and is associated with poorer overall survival compared with other subtypes. Given the lack of targets for current tailored therapies in TNP tumours, chemotherapy is the only systemic treatment available; however, overall outcomes remain poor. Therefore, optimal treatment regimens and targeted therapies are urgently needed.

OBJECTIVE

We discuss characteristics of TNP cancers that underpin the rationale of current and novel therapeutic strategies, and an approach for finding and validating new therapeutic targets.

RESULTS/CONCLUSION

The results of large prospective randomised controlled trials are currently awaited. Efforts to unravel the heterogeneity and complexity of TNP cancers using the latest high-throughput molecular techniques and integrating these findings with biology-driven therapeutic strategies in clinical trials will be of paramount importance for the development of treatment approaches for this breast cancer subtype.

Systems biology and cancer prevention: all options on the table

In this paper, we outline the status quo and approaches to further development of the systems biology concepts with focus on applications in cancer prevention science. We discuss the biological aspects of cancer research that are of primary importance in cancer prevention, motivations for their mathematical modeling and some recent advances in computational oncology. We also make an attempt to outline in big conceptual terms the contours of future work aimed at creation of large-scale computational and informational infrastructure for using as a routine tool in cancer prevention science and decision making.

Computer modelling of epilepsy

Epilepsy is a complex set of disorders that can involve many areas of the cortex, as well as underlying deep-brain systems. The myriad manifestations of seizures, which can be as varied as déjà vu and olfactory hallucination, can therefore give researchers insights into regional functions and relations. Epilepsy is also complex genetically and pathophysiologically: it involves microscopic (on the scale of ion channels and synaptic proteins), macroscopic (on the scale of brain trauma and rewiring) and intermediate changes in a complex interplay of causality. It has long been recognized that computer modelling will be required to disentangle causality, to better understand seizure spread and to understand and eventually predict treatment efficacy. Over the past few years, substantial progress has been made in modelling epilepsy at levels ranging from the molecular to the socioeconomic. We review these efforts and connect them to the medical goals of understanding and treating the disorder.

Using genetic markers to orient the edges in quantitative trait networks: the NEO software

Background

Systems genetic studies have been used to identify genetic loci that affect transcript abundances and clinical traits such as body weight. The pairwise correlations between gene expression traits and/or clinical traits can be used to define undirected trait networks. Several authors have argued that genetic markers (e.g expression quantitative trait loci, eQTLs) can serve as causal anchors for orienting the edges of a trait network. The availability of hundreds of thousands of genetic markers poses new challenges: how to relate (anchor) traits to multiple genetic markers, how to score the genetic evidence in favor of an edge orientation, and how to weigh the information from multiple markers.

Results

We develop and implement Network Edge Orienting (NEO) methods and software that address the challenges of inferring unconfounded and directed gene networks from microarray-derived gene expression data by integrating mRNA levels with genetic marker data and Structural Equation Model (SEM) comparisons. The NEO software implements several manual and automatic methods for incorporating genetic information to anchor traits. The networks are oriented by considering each edge separately, thus reducing error propagation. To summarize the genetic evidence in favor of a given edge orientation, we propose Local SEM-based Edge Orienting (LEO) scores that compare the fit of several competing causal graphs. SEM fitting indices allow the user to assess local and overall model fit. The NEO software allows the user to carry out a robustness analysis with regard to genetic marker selection. We demonstrate the utility of NEO by recovering known causal relationships in the sterol homeostasis pathway using liver gene expression data from an F2 mouse cross. Further, we use NEO to study the relationship between a disease gene and a biologically important gene co-expression module in liver tissue.

Conclusion

The NEO software can be used to orient the edges of gene co-expression networks or quantitative trait networks if the edges can be anchored to genetic marker data. R software tutorials, data, and supplementary material can be downloaded from: .

MetaReg: a platform for modeling, analysis and visualization of biological systems using large-scale experimental data

MetaReg http://acgt.cs.tau.ac.il/metareg/application.html is a computational tool that models cellular networks and integrates experimental results with such models. MetaReg represents established knowledge about a biological system, available today mostly in informal form in the literature, as probabilistic network models with underlying combinatorial regulatory logic. MetaReg enables contrasting predictions with measurements, model improvements and studying what-if scenarios. By summarizing prior knowledge and providing visual and computational aids, it helps the expert explore and understand her system better.

Computational systems biology in cancer: modeling methods and applications

In recent years it has become clear that carcinogenesis is a complex process, both at the molecular and cellular levels. Understanding the origins, growth and spread of cancer, therefore requires an integrated or system-wide approach. Computational systems biology is an emerging sub-discipline in systems biology that utilizes the wealth of data from genomic, proteomic and metabolomic studies to build computer simulations of intra and intercellular processes. Several useful descriptive and predictive models of the origin, growth and spread of cancers have been developed in an effort to better understand the disease and potential therapeutic approaches. In this review we describe and assess the practical and theoretical underpinnings of commonly-used modeling approaches, including ordinary and partial differential equations, petri nets, cellular automata, agent based models and hybrid systems. A number of computer-based formalisms have been implemented to improve the accessibility of the various approaches to researchers whose primary interest lies outside of model development. We discuss several of these and describe how they have led to novel insights into tumor genesis, growth, apoptosis, vascularization and therapy.

A multi-layered approach to protein data integration for diabetes research

OBJECTIVE

Recent advances in high-throughput experimental techniques have enabled many protein-protein interactions to be identified and stored in large databases. Understanding protein interactions is fundamental to the advancement of science and medical knowledge, unfortunately the scale of the requires an automated approach to analysis. We describe our graph-mining techniques to identify important structures within protein-protein interaction networks to aid in human comprehension and computerised analysis.

METHODS AND MATERIALS

We describe our techniques for characterizing graph type and associated properties which is constructed from data collated from the Human Protein Reference Database. Using random graph rewiring comparative techniques and cross-validation with other identification methods a further analysis of the identified essential proteins is presented to illustrate the accuracy of these measures. We argue for using techniques based upon graph structure for separating and encapsulating proteins based upon functionality.

RESULTS

We demonstrate how rational Erdos numbers may be used as a method to identify collaborating proteins based solely upon network structure. Further, by using dynamic cut-off limit it demonstrates how collaboration subgraphs can be generated for each protein within the network, and how graph containment can be used as a means of identifying which of many possible graphs are likely to be actual protein complexes. The demonstration protein interaction network built for diabetes is found to be a scale-free, small-world graph with a power-law degree distribution of interactions on nodes. These findings are consistent with many other protein interaction networks.

Getting it right: designing microarray (and not 'microawry') comparative genomic hybridization studies for cancer research

The development of high-resolution microarray-based comparative genomic hybridization (aCGH), using cDNA, bacterial artificial chromosome (BAC) and oligonucleotide probes, is providing tremendous opportunities for translational research by facilitating detailed analysis of entire cancer genomes in a single experiment. However, this technology will only fulfil its promise if studies incorporating aCGH are designed with a full understanding of its current limitations and the strategies available to circumvent them. While there have been several excellent reviews on the current status of this technology, there is currently very little guidance available regarding the appropriate design of experiments incorporating aCGH (including the strengths and weaknesses of each platform), and how best to combine the results obtained from aCGH with other 'omic' technologies, including gene expression. In this review, we present the key design issues that need to be considered in order to optimize aCGH studies, including sample selection, the definition of appropriate experimental objectives, arguments for and against the various microarray platforms that are currently available, and methods for data validation and integration. It is envisaged that future well-designed aCGH studies will enhance our understanding of the genetic basis of cancer, and lead to the identification of novel predictive and prognostic cancer biomarkers, as well as molecular therapeutic targets in cancer.

Cancer epidemiology: study designs and data analysis

Among the scientific interests of cancer epidemiology is the identification of both environmental and genetic factors associated with cancer development. Observational designs requiring sophisticated methodology are applied to control for potential confounding factors. The enormous biotechnological potential developed in the last two decades has allowed the integration of a plethora of new biomarkers in epidemiological studies to better define the exposure and "neoclassic" outcomes, as well as incorporating genetic susceptibility factors in both classical and new epidemiological designs. The integration of scopes, objectives, data and tools coming from different disciplines also benefits epidemiology, thus evolving into "systems epidemiology". In this manuscript, we review the basic concepts of study designs and data analysis and introduce readers to the more innovative aspects that are now being applied in epidemiological studies.