Computational modeling of mammalian signaling networks

Dynamic Modeling of Cell Cycle Arrest Through Integrated Single-Cell and Mathematical Modelling Approaches

Highly multiplexed imaging assays allow simultaneous quantification of multiple protein and phosphorylation markers, providing a static snapshots of cell types and states. Pseudo-time techniques can transform these static snapshots of unsynchronized cells into dynamic trajectories, enabling the study of dynamic processes such as development trajectories and the cell cycle. Such ordering also enables training of mathematical models on these data, but technical challenges have hitherto made it difficult to integrate multiple experimental conditions, limiting the predictive power and insights these models can generate. In this work, we propose data processing and model training approaches for integrating multiplexed, multi-condition immunofluorescence data with mathematical modelling. We devise training strategies that are applicable to datasets where cells exhibit oscillatory as well as arrested dynamics and use them to train a cell cycle model on a dataset of MCF-10A mammary epithelial exposed to cell-cycle arresting small molecules. We validate the model by investigating predicted growth factor sensitivities and responses to inhibitors of cells at different initial conditions. We anticipate that our framework will generalise to other highly multiplexed measurement techniques such as mass-cytometry, rendering larger bodies of data accessible to dynamic modelling and paving the way to deeper biological insights.

Mechanistic modeling of cell viability assays with in silico lineage tracing

Data from cell viability assays, which measure cumulative division and death events in a population and reflect substantial cellular heterogeneity, are widely available. However, interpreting such data with mechanistic computational models is hindered because direct model/data comparison is often muddled. We developed an algorithm that tracks simulated division and death events in mechanistically detailed single-cell lineages to enable such a model/data comparison and suggest causes of cell-cell drug response variability. Using our previously developed model of mammalian single-cell proliferation and death signaling, we simulated drug dose response experiments for four targeted anti-cancer drugs (alpelisib, neratinib, trametinib and palbociclib) and compared them to experimental data. Simulations are consistent with data for strong growth inhibition by trametinib (MEK inhibitor) and overall lack of efficacy for alpelisib (PI-3K inhibitor), but are inconsistent with data for palbociclib (CDK4/6 inhibitor) and neratinib (EGFR inhibitor). Model/data inconsistencies suggest (i) the importance of CDK4/6 for driving the cell cycle may be overestimated, and (ii) that the cellular balance between basal (tonic) and ligand-induced signaling is a critical determinant of receptor inhibitor response. Simulations show subpopulations of rapidly and slowly dividing cells in both control and drug-treated conditions. Variations in mother cells prior to drug treatment all impinging on ERK pathway activity are associated with the rapidly dividing phenotype and trametinib resistance. This work lays a foundation for the application of mechanistic modeling to large-scale cell viability assay datasets and better understanding determinants of cellular heterogeneity in drug response.

Self-consistent signal transduction analysis for modeling context-specific signaling cascades and perturbations

Biological signal transduction networks are central to information processing and regulation of gene expression across all domains of life. Dysregulation is known to cause a wide array of diseases, including cancers. Here I introduce self-consistent signal transduction analysis, which utilizes genome-scale -omics data (specifically transcriptomics and/or proteomics) in order to predict the flow of information through these networks in an individualized manner. I apply the method to the study of endocrine therapy in breast cancer patients, and show that drugs that inhibit estrogen receptor α elicit a wide array of antitumoral effects, and that their most clinically-impactful ones are through the modulation of proliferative signals that control the genes GREB1, HK1, AKT1, MAPK1, AKT2, and NQO1. This method offers researchers a valuable tool in understanding how and why dysregulation occurs, and how perturbations to the network (such as targeted therapies) effect the network itself, and ultimately patient outcomes.

Investigating In silico and In vitro Therapeutic Potential of Diosmetin as the Anti-Parkinson Agent

AIM

This study aimed to investigate how diosmetin interacts with seven target receptors associated with oxidative stress (OS) and validate its antioxidant properties for the potential management of Parkinson's disease (PD).

BACKGROUND

In PD, the degeneration of dopaminergic cells is strongly influenced by OS. This stressor is intricately connected to various mechanisms involved in neurodegeneration, such as mitochondrial dysfunction, neuroinflammation, and excitotoxicity induced by nitric oxide.

OBJECTIVE

The aim of this research was to establish a molecular connection between diosmetin and OS-associated target receptors was the goal, and it investigated how this interaction can lessen PD.

METHODS

Seven molecular targets - Adenosine A2A (AA2A), Peroxisome Proliferator-Activated Receptor Gamma (PPARγ), Protein Kinase AKT1, Nucleolar Receptor NURR1, Liver - X Receptor Beta (LXRβ), Monoamine Oxidase - B (MAO-B) and Tropomyosin receptor kinase B (TrkB) were obtained from RCSB. Molecular docking software was employed to determine molecular interactions, while antioxidant activity was assessed through in vitro assays against various free radicals.

RESULTS

Diosmetin exhibited interactions with all seven target receptors at their binding sites. Notably, it showed superior interaction with AA2A and NURR1 compared to native ligands, with binding energies of -7.55, and -6.34 kcal/mol, respectively. Additionally, significant interactions were observed with PPARγ, AKT1, LXRβ, MAO-B, and TrkB with binding energies of -8.34, -5.42, -7.66, -8.82, -8.45 kcal/mol, respectively. Diosmetin also demonstrated antioxidant activity against various free radicals, particularly against hypochlorous acid (HOCl) and nitric oxide (NO) free radicals.

CONCLUSION

Diosmetin possibly acts on several target receptors linked to the pathophysiology of PD, demonstrating promise as an OS inhibitor and scavenger.

Computational modeling of signal transduction networks without kinetic parameters: Petri net approaches

More and more computational techniques have been applied to model biological systems, especially signaling pathways in medical systems. Due to the large number of experimental data driven by high-throughput technologies, new computational concepts have been developed. Nevertheless, often the necessary kinetic data cannot be determined in sufficient number and quality because of experimental complexity or ethical reasons. At the same time, the number of qualitative data drastically increased, for example, gene expression data, protein-protein interaction data, and imaging data. Especially for large-scale models, the application of kinetic modeling techniques can fail. On the other hand, many large-scale models have been constructed applying qualitative and semiquantitative techniques, for example, logical models or Petri net models. These techniques make it possible to explore system's dynamics without knowing kinetic parameters. Here, we summarize the work of the last 10 years for modeling signal transduction pathways in medical applications applying Petri net formalism. We focus on analysis techniques based on system's invariants without any kinetic parameters and show predictions of all signaling pathways of the system. We start with an intuitive introduction into Petri nets and system's invariants. We illustrate the main concepts using the tumor necrosis factor receptor 1 (TNFR1)-induced nuclear factor κ-light-chain-enhancer of activated B cells (NF-κB) pathway as a case study. Summarizing recent models, we discuss the advantages and challenges of Petri net applications to medical signaling systems. In addition, we provide exemplarily interesting Petri net applications to model signaling in medical systems of the last years that use the well-known stochastic and kinetic concepts developed about 50 years ago.

A forecast for large-scale, predictive biology: Lessons from meteorology

Quantitative systems biology, in which predictive mathematical models are constructed to guide the design of experiments and predict experimental outcomes, is at an exciting transition point, where the foundational scientific principles are becoming established, but the impact is not yet global. The next steps necessary for mathematical modeling to transform biological research and applications, in the same way it has already transformed other fields, is not completely clear. The purpose of this perspective is to forecast possible answers to this question-what needs to happen next-by drawing on the experience gained in another field, specifically meteorology. We review here a number of lessons learned in weather prediction that are directly relevant to biological systems modeling, and that we believe can enable the same kinds of global impact in our field as atmospheric modeling makes today.

Identification of major depressive disorder disease-related genes and functional pathways based on system dynamic changes of network connectivity

BACKGROUND

Major depressive disorder (MDD) is a leading psychiatric disorder that involves complex abnormal biological functions and neural networks. This study aimed to compare the changes in the network connectivity of different brain tissues under different pathological conditions, analyzed the biological pathways and genes that are significantly related to disease progression, and further predicted the potential therapeutic drug targets.

METHODS

Expression of differentially expressed genes (DEGs) were analyzed with postmortem cingulate cortex (ACC) and prefrontal cortex (PFC) mRNA expression profile datasets downloaded from the Gene Expression Omnibus (GEO) database, including 76 MDD patients and 76 healthy subjects in ACC and 63 MDD patients and 63 healthy subjects in PFC. The co-expression network construction was based on system network analysis. The function of the genes was annotated by Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. Human Protein Reference Database (HPRD, http://www.hprd.org/ ) was used for gene interaction relationship mapping.

RESULTS

We filtered 586 DEGs in ACC and 616 DEGs in PFC for further analysis. By constructing the co-expression network, we found that the gene connectivity was significantly reduced under disease conditions (P = 0.04 in PFC and P = 1.227e-09 in ACC). Crosstalk analysis showed that CD19, PTDSS2 and NDST2 were significantly differentially expressed in ACC and PFC of MDD patients. Among them, CD19 and PTDSS2 have been targeted by several drugs in the Drugbank database. KEGG pathway analysis demonstrated that the function of CD19 and PTDSS2 were enriched with the pathway of Glycerophospholipid metabolism and T cell receptor signaling pathway.

CONCLUSION

Co-expression network and tissue comparing analysis can identify signaling pathways and cross talk genes related to MDD, which may provide novel insight for understanding the molecular mechanisms of MDD.

Transcriptomic and Morphological Analysis of Cells Derived from Porcine Buccal Mucosa-Studies on an In Vitro Model

Simple Summary

Domestic pigs express high phylogenetic similarity to humans and are often used as a compatible model in biomedical research. Porcine tissues are used as an accessible biomaterial in human skin transplants and tissue architecture reconstruction. We used transcriptional analysis to investigate the dynamics of complex biological system of the mucosa. Additionally, we performed computer analysis of microscopic images of cultured cells in vitro. Computer analysis of images identified epithelial cells and connective tissue cells in in vitro culture.

Abstract

Transcriptional analysis and live-cell imaging are a powerful tool to investigate the dynamics of complex biological systems. In vitro expanded porcine oral mucosal cells, consisting of populations of epithelial and connective lineages, are interesting and complex systems for study via microarray transcriptomic assays to analyze gene expression profile. The transcriptomic analysis included 56 ontological groups with particular focus on 7 gene ontology groups that are related to the processes of differentiation and development. Most analyzed genes were upregulated after 7 days and downregulated after 15 and 30 days of in vitro culture. The performed transcriptomic analysis was then extended to include automated analysis of differential interference contrast microscopy (DIC) images obtained during in vitro culture. The analysis of DIC imaging allowed to identify the different populations of keratinocytes and fibroblasts during seven days of in vitro culture, and it was possible to evaluate the proportion of these two populations of cells. Porcine mucosa may be a suitable model for reference research on human tissues. In addition, it can provide a reference point for research on the use of cells, scaffolds, or tissues derived from transgenic animals for applications in human tissues reconstruction.

Histamine-induced biphasic activation of RhoA allows for persistent RhoA signaling

The small GTPase RhoA is a central signaling enzyme that is involved in various cellular processes such as cytoskeletal dynamics, transcription, and cell cycle progression. Many signal transduction pathways activate RhoA—for instance, Gα_q-coupled Histamine 1 Receptor signaling via Gα_q-dependent activation of RhoGEFs such as p63. Although multiple upstream regulators of RhoA have been identified, the temporal regulation of RhoA and the coordination of different upstream components in its regulation have not been well characterized. In this study, live-cell measurement of RhoA activation revealed a biphasic increase of RhoA activity upon histamine stimulation. We showed that the first and second phase of RhoA activity are dependent on p63 and Ca²⁺/PKC, respectively, and further identified phosphorylation of serine 240 on p115 RhoGEF by PKC to be the mechanistic link between PKC and RhoA. Combined approaches of computational modeling and quantitative measurement revealed that the second phase of RhoA activation is insensitive to rapid turning off of the receptor and is required for maintaining RhoA-mediated transcription after the termination of the receptor signaling. Thus, two divergent pathways enable both rapid activation and persistent signaling in receptor-mediated RhoA signaling via intricate temporal regulation.

The small GTPase RhoA is a central signaling enzyme that is involved in various cellular processes such as cytoskeletal dynamics, transcription, and cell cycle progression. This study shows that histamine induces biphasic activation of RhoA via two divergent signaling pathways, allowing for intricate regulation of cellular processes.

Profiling Cell Signaling Networks at Single-cell Resolution

Signaling networks process intra- and extracellular information to modulate the functions of a cell. Deregulation of signaling networks results in abnormal cellular physiological states and often drives diseases. Network responses to a stimulus or a drug treatment can be highly heterogeneous across cells in a tissue because of many sources of cellular genetic and non-genetic variance. Signaling network heterogeneity is the key to many biological processes, such as cell differentiation and drug resistance. Only recently, the emergence of multiplexed single-cell measurement technologies has made it possible to evaluate this heterogeneity. In this review, we categorize currently established single-cell signaling network profiling approaches by their methodology, coverage, and application, and we discuss the advantages and limitations of each type of technology. We also describe the available computational tools for network characterization using single-cell data and discuss potential confounding factors that need to be considered in single-cell signaling network analyses.

Connecting signaling and metabolic pathways in EGF receptor-mediated oncogenesis of glioblastoma

As malignant transformation requires synchronization of growth-driving signaling (S) and metabolic (M) pathways, defining cancer-specific S-M interconnected networks (SMINs) could lead to better understanding of oncogenic processes. In a systems-biology approach, we developed a mathematical model for SMINs in mutated EGF receptor (EGFRvIII) compared to wild-type EGF receptor (EGFRwt) expressing glioblastoma multiforme (GBM). Starting with experimentally validated human protein-protein interactome data for S-M pathways, and incorporating proteomic data for EGFRvIII and EGFRwt GBM cells and patient transcriptomic data, we designed a dynamic model for EGFR-driven GBM-specific information flow. Key nodes and paths identified by in silico perturbation were validated experimentally when inhibition of signaling pathway proteins altered expression of metabolic proteins as predicted by the model. This demonstrated capacity of the model to identify unknown connections between signaling and metabolic pathways, explain the robustness of oncogenic SMINs, predict drug escape, and assist identification of drug targets and the development of combination therapies.

Complex and highly dynamic interconnected networks allow cancer to take different routes and circumvent chemotherapy. Therefore, understanding these context-specific networks and their dynamics of molecular interactions driven by different oncogenic signaling and metabolic pathways is very much needed to predict drug targets and the effect of therapeutics. We incorporated high-throughput transcriptome and proteome data into mathematical models to deduce properties of cancer cells through systems biology approach. Here we report the development, testing and validation of an integrated systems biology model of information flow between signaling and metabolic pathways to understand the regulation of the interconnection between them in cancer. Our model efficiently identified unique connections and key nodes important in signaling-metabolic information flow. We predicted some potential novel targets before performing actual drug tests. We have successfully applied this model to identify the interconnections altered in the constitutive signaling of the mutated EGFR by comparing EGF-dependent and wild-type EGFR signaling in glioblastoma multiforme.

Analysis of the Human Kinome and Phosphatome by Mass Cytometry Reveals Overexpression-Induced Effects on Cancer-Related Signaling

Kinase and phosphatase overexpression drives tumorigenesis and drug resistance. We previously developed a mass-cytometry-based single-cell proteomics approach that enables quantitative assessment of overexpression effects on cell signaling. Here, we applied this approach in a human kinome- and phosphatome-wide study to assess how 649 individually overexpressed proteins modulated cancer-related signaling in HEK293T cells in an abundance-dependent manner. Based on these data, we expanded the functional classification of human kinases and phosphatases and showed that the overexpression effects include non-catalytic roles. We detected 208 previously unreported signaling relationships. The signaling dynamics analysis indicated that the overexpression of ERK-specific phosphatases sustains proliferative signaling. This suggests a phosphatase-driven mechanism of cancer progression. Moreover, our analysis revealed a drug-resistant mechanism through which overexpression of tyrosine kinases, including SRC, FES, YES1, and BLK, induced MEK-independent ERK activation in melanoma A375 cells. These proteins could predict drug sensitivity to BRAF-MEK concurrent inhibition in cells carrying BRAF mutations.

Tamoxifen-Induced Apoptosis of MCF-7 Cells via GPR30/PI3K/MAPKs Interactions: Verification by ODE Modeling and RNA Sequencing

Tamoxifen (Nolvadex) is one of the most widely used and effective therapeutic agent for breast cancer. It benefits nearly 75% of patients with estrogen receptor (ER)-positive breast cancer that receive this drug. Its effectiveness is mainly attributed to its capacity to function as an ER antagonist, blocking estrogen binding sites on the receptor, and inhibiting the proliferative action of the receptor-hormone complex. Although, tamoxifen can induce apoptosis in breast cancer cells via upregulation of pro-apoptotic factors, it can also promote uterine hyperplasia in some women. Thus, tamoxifen as a multi-functional drug could have different effects on cells based on the utilization of effective concentrations or availability of specific co-factors. Evidence that tamoxifen functions as a GPR30 (G-Protein Coupled Receptor 30) agonist activating adenylyl cyclase and EGFR (Epidermal Growth Factor Receptor) intracellular signaling networks, provides yet another means of explaining the multi-functionality of tamoxifen. Here ordinary differential equation (ODE) modeling, RNA sequencing and real time qPCR analysis were utilized to establish the necessary data for gene network mapping of tamoxifen-stimulated MCF-7 cells, which express the endogenous ER and GPR30. The gene set enrichment analysis and pathway analysis approaches were used to categorize transcriptionally upregulated genes in biological processes. Of the 2,713 genes that were significantly upregulated following a 48 h incubation with 250 μM tamoxifen, most were categorized as either growth-related or pro-apoptotic intermediates that fit into the Tp53 and/or MAPK signaling pathways. Collectively, our results display that the effects of tamoxifen on the breast cancer MCF-7 cell line are mediated by the activation of important signaling pathways including Tp53 and MAPKs to induce apoptosis.

A blueprint for human whole-cell modeling

Whole-cell dynamical models of human cells are a central goal of systems biology. Such models could help researchers understand cell biology and help physicians treat disease. Despite significant challenges, we believe that human whole-cell models are rapidly becoming feasible. To develop a plan for achieving human whole-cell models, we analyzed the existing models of individual cellular pathways, surveyed the biomodeling community, and reflected on our experience developing whole-cell models of bacteria. Based on these analyses, we propose a plan for a project, termed the Human Whole-Cell Modeling Project, to achieve human whole-cell models. The foundations of the plan include technology development, standards development, and interdisciplinary collaboration.

Late-phase synthesis of IκBα insulates the TLR4-activated canonical NF-κB pathway from noncanonical NF-κB signaling in macrophages

The nuclear factor κB (NF-κB) transcription factors coordinate the inflammatory immune response during microbial infection. Pathogenic substances engage canonical NF-κB signaling through the heterodimer RelA:p50, which is subjected to rapid negative feedback by inhibitor of κBα (IκBα). The noncanonical NF-κB pathway is required for the differentiation of immune cells; however, cross-talk between both pathways can occur. Concomitantly activated noncanonical signaling generates p52 from the p100 precursor. The synthesis of p100 is induced by canonical signaling, leading to the formation of the late-acting RelA:p52 heterodimer. This cross-talk prolongs inflammatory RelA activity in epithelial cells to ensure pathogen clearance. We found that the Toll-like receptor 4 (TLR4)-activated canonical NF-κB signaling pathway is insulated from lymphotoxin β receptor (LTβR)-induced noncanonical signaling in mouse macrophage cell lines. Combined computational and biochemical studies indicated that the extent of NF-κB-responsive expression of Nfkbia, which encodes IκBα, inversely correlated with cross-talk. The Nfkbia promoter showed enhanced responsiveness to NF-κB activation in macrophages compared to that in fibroblasts. We found that this hyperresponsive promoter engaged the RelA:p52 dimer generated during costimulation of macrophages through TLR4 and LTβR to trigger synthesis of IκBα at late time points, which prevented the late-acting RelA cross-talk response. Together, these data suggest that, despite the presence of identical signaling networks in cells of diverse lineages, emergent cross-talk between signaling pathways is subject to cell type-specific regulation. We propose that the insulation of canonical and noncanonical NF-κB pathways limits the deleterious effects of macrophage-mediated inflammation.

Deep Learning Automates the Quantitative Analysis of Individual Cells in Live-Cell Imaging Experiments

Live-cell imaging has opened an exciting window into the role cellular heterogeneity plays in dynamic, living systems. A major critical challenge for this class of experiments is the problem of image segmentation, or determining which parts of a microscope image correspond to which individual cells. Current approaches require many hours of manual curation and depend on approaches that are difficult to share between labs. They are also unable to robustly segment the cytoplasms of mammalian cells. Here, we show that deep convolutional neural networks, a supervised machine learning method, can solve this challenge for multiple cell types across the domains of life. We demonstrate that this approach can robustly segment fluorescent images of cell nuclei as well as phase images of the cytoplasms of individual bacterial and mammalian cells from phase contrast images without the need for a fluorescent cytoplasmic marker. These networks also enable the simultaneous segmentation and identification of different mammalian cell types grown in co-culture. A quantitative comparison with prior methods demonstrates that convolutional neural networks have improved accuracy and lead to a significant reduction in curation time. We relay our experience in designing and optimizing deep convolutional neural networks for this task and outline several design rules that we found led to robust performance. We conclude that deep convolutional neural networks are an accurate method that require less curation time, are generalizable to a multiplicity of cell types, from bacteria to mammalian cells, and expand live-cell imaging capabilities to include multi-cell type systems.

Knowledge-guided fuzzy logic modeling to infer cellular signaling networks from proteomic data

Modeling of signaling pathways is crucial for understanding and predicting cellular responses to drug treatments. However, canonical signaling pathways curated from literature are seldom context-specific and thus can hardly predict cell type-specific response to external perturbations; purely data-driven methods also have drawbacks such as limited biological interpretability. Therefore, hybrid methods that can integrate prior knowledge and real data for network inference are highly desirable. In this paper, we propose a knowledge-guided fuzzy logic network model to infer signaling pathways by exploiting both prior knowledge and time-series data. In particular, the dynamic time warping algorithm is employed to measure the goodness of fit between experimental and predicted data, so that our method can model temporally-ordered experimental observations. We evaluated the proposed method on a synthetic dataset and two real phosphoproteomic datasets. The experimental results demonstrate that our model can uncover drug-induced alterations in signaling pathways in cancer cells. Compared with existing hybrid models, our method can model feedback loops so that the dynamical mechanisms of signaling networks can be uncovered from time-series data. By calibrating generic models of signaling pathways against real data, our method supports precise predictions of context-specific anticancer drug effects, which is an important step towards precision medicine.

60 YEARS OF POMC: Adrenal and extra-adrenal functions of ACTH

The pituitary adrenocorticotropic hormone (ACTH) plays a pivotal role in homeostasis and stress response and is thus the major component of the hypothalamo-pituitary-adrenal axis. After a brief summary of ACTH production from proopiomelanocortin (POMC) and on ACTH receptor properties, the first part of the review covers the role of ACTH in steroidogenesis and steroid secretion. We highlight the mechanisms explaining the differential acute vs chronic effects of ACTH on aldosterone and glucocorticoid secretion. The second part summarizes the effects of ACTH on adrenal growth, addressing its role as either a mitogenic or a differentiating factor. We then review the mechanisms involved in steroid secretion, from the classical Cyclic adenosine monophosphate second messenger system to various signaling cascades. We also consider how the interaction between the extracellular matrix and the cytoskeleton may trigger activation of signaling platforms potentially stimulating or repressing the steroidogenic potency of ACTH. Finally, we consider the extra-adrenal actions of ACTH, in particular its role in differentiation in a variety of cell types, in addition to its known lipolytic effects on adipocytes. In each section, we endeavor to correlate basic mechanisms of ACTH function with the pathological consequences of ACTH signaling deficiency and of overproduction of ACTH.

Multi-label multi-instance transfer learning for simultaneous reconstruction and cross-talk modeling of multiple human signaling pathways

BACKGROUND

Signaling pathways play important roles in the life processes of cell growth, cell apoptosis and organism development. At present the signal transduction networks are far from complete. As an effective complement to experimental methods, computational modeling is suited to rapidly reconstruct the signaling pathways at low cost. To our knowledge, the existing computational methods seldom simultaneously exploit more than three signaling pathways into one predictive model for the discovery of novel signaling components and the cross-talk modeling between signaling pathways.

RESULTS

In this work, we propose a multi-label multi-instance transfer learning method to simultaneously reconstruct 27 human signaling pathways and model their cross-talks. Computational results show that the proposed method demonstrates satisfactory multi-label learning performance and rational proteome-wide predictions. Some predicted signaling components or pathway targeted proteins have been validated by recent literature. The predicted signaling components are further linked to pathways using the experimentally derived PPIs (protein-protein interactions) to reconstruct the human signaling pathways. Thus the map of the cross-talks via common signaling components and common signaling PPIs is conveniently inferred to provide valuable insights into the regulatory and cooperative relationships between signaling pathways. Lastly, gene ontology enrichment analysis is conducted to gain statistical knowledge about the reconstructed human signaling pathways.

CONCLUSIONS

Multi-label learning framework has been demonstrated effective in this work to model the phenomena that a signaling protein belongs to more than one signaling pathway. As results, novel signaling components and pathways targeted proteins are predicted to simultaneously reconstruct multiple human signaling pathways and the static map of their cross-talks for further biomedical research.

A simple feature construction method for predicting upstream/downstream signal flow in human protein-protein interaction networks

Signaling pathways play important roles in understanding the underlying mechanism of cell growth, cell apoptosis, organismal development and pathways-aberrant diseases. Protein-protein interaction (PPI) networks are commonly-used infrastructure to infer signaling pathways. However, PPI networks generally carry no information of upstream/downstream relationship between interacting proteins, which retards our inferring the signal flow of signaling pathways. In this work, we propose a simple feature construction method to train a SVM (support vector machine) classifier to predict PPI upstream/downstream relations. The domain based asymmetric feature representation naturally embodies domain-domain upstream/downstream relations, providing an unconventional avenue to predict the directionality between two objects. Moreover, we propose a semantically interpretable decision function and a macro bag-level performance metric to satisfy the need of two-instance depiction of an interacting protein pair. Experimental results show that the proposed method achieves satisfactory cross validation performance and independent test performance. Lastly, we use the trained model to predict the PPIs in HPRD, Reactome and IntAct. Some predictions have been validated against recent literature.

Quantitative analysis of receptor tyrosine kinase-effector coupling at functionally relevant stimulus levels

A major goal of current signaling research is to develop a quantitative understanding of how receptor activation is coupled to downstream signaling events and to functional cellular responses. Here, we measure how activation of the RET receptor tyrosine kinase on mouse neuroblastoma cells by the neurotrophin artemin (ART) is quantitatively coupled to key downstream effectors. We show that the efficiency of RET coupling to ERK and Akt depends strongly on ART concentration, and it is highest at the low (∼100 pM) ART levels required for neurite outgrowth. Quantitative discrimination between ERK and Akt pathway signaling similarly is highest at this low ART concentration. Stimulation of the cells with 100 pM ART activated RET at the rate of ∼10 molecules/cell/min, leading at 5-10 min to a transient peak of ∼150 phospho-ERK (pERK) molecules and ∼50 pAkt molecules per pRET, after which time the levels of these two signaling effectors fell by 25-50% while the pRET levels continued to slowly rise. Kinetic experiments showed that signaling effectors in different pathways respond to RET activation with different lag times, such that the balance of signal flux among the different pathways evolves over time. Our results illustrate that measurements using high, super-physiological growth factor levels can be misleading about quantitative features of receptor signaling. We propose a quantitative model describing how receptor-effector coupling efficiency links signal amplification to signal sensitization between receptor and effector, thereby providing insight into design principles underlying how receptors and their associated signaling machinery decode an extracellular signal to trigger a functional cellular outcome.

Methodologies for the modeling and simulation of biochemical networks, illustrated for signal transduction pathways: a primer

Biochemical networks depict the chemical interactions that take place among elements of living cells. They aim to elucidate how cellular behavior and functional properties of the cell emerge from the relationships between its components, i.e. molecules. Biochemical networks are largely characterized by dynamic behavior, and exhibit high degrees of complexity. Hence, the interest in such networks is growing and they have been the target of several recent modeling efforts. Signal transduction pathways (STPs) constitute a class of biochemical networks that receive, process, and respond to stimuli from the environment, as well as stimuli that are internal to the organism. An STP consists of a chain of intracellular signaling processes that ultimately result in generating different cellular responses. This primer presents the methodologies used for the modeling and simulation of biochemical networks, illustrated for STPs. These methodologies range from qualitative to quantitative, and include structural as well as dynamic analysis techniques. We describe the different methodologies, outline their underlying assumptions, and provide an assessment of their advantages and disadvantages. Moreover, publicly and/or commercially available implementations of these methodologies are listed as appropriate. In particular, this primer aims to provide a clear introduction and comprehensive coverage of biochemical modeling and simulation methodologies for the non-expert, with specific focus on relevant literature of STPs.

Visualization of the spatial and temporal dynamics of MAPK signaling using fluorescence imaging techniques

Conserved mitogen-activated protein kinase (MAPK) signaling pathways are major mechanisms through which cells perceive and respond properly to their surrounding environment. Such homeostatic responses maintain the life of the organism. Since errors in MAPK signaling pathways can lead to cancers and to defects in immune responses, in the nervous system and metabolism, these pathways have been extensively studied as potential therapeutic targets. Although much has been studied about the roles of MAPKs in various cellular functions, less is known regarding regulation of MAPK in living organisms. This review will focus on the latest understanding of the dynamic regulation of MAPK signaling in intact cells that was revealed by using novel fluorescence imaging techniques and advanced systems-analytical methods. These techniques allowed quantitative analyses of signal transduction in situ with high spatio-temporal resolution and have revealed the nature of the molecular dynamics that determine cellular responses and fates.

Experimental design schemes for learning Boolean network models

MOTIVATION

A holy grail of biological research is a working model of the cell. Current modeling frameworks, especially in the protein-protein interaction domain, are mostly topological in nature, calling for stronger and more expressive network models. One promising alternative is logic-based or Boolean network modeling, which was successfully applied to model signaling regulatory circuits in human. Learning such models requires observing the system under a sufficient number of different conditions. To date, the amount of measured data is the main bottleneck in learning informative Boolean models, underscoring the need for efficient experimental design strategies.

RESULTS

We developed novel design approaches that greedily select an experiment to be performed so as to maximize the difference or the entropy in the results it induces with respect to current best-fit models. Unique to our maximum difference approach is the ability to account for all (possibly exponential number of) Boolean models displaying high fit to the available data. We applied both approaches to simulated and real data from the EFGR and IL1 signaling systems in human. We demonstrate the utility of the developed strategies in substantially improving on a random selection approach. Our design schemes highlight the redundancy in these datasets, leading up to 11-fold savings in the number of experiments to be performed.

AVAILABILITY AND IMPLEMENTATION

Source code will be made available upon acceptance of the manuscript.

Genome scale modeling in systems biology: algorithms and resources

In recent years, in silico studies and trial simulations have complemented experimental procedures. A model is a description of a system, and a system is any collection of interrelated objects; an object, moreover, is some elemental unit upon which observations can be made but whose internal structure either does not exist or is ignored. Therefore, any network analysis approach is critical for successful quantitative modeling of biological systems. This review highlights some of most popular and important modeling algorithms, tools, and emerging standards for representing, simulating and analyzing cellular networks in five sections. Also, we try to show these concepts by means of simple example and proper images and graphs. Overall, systems biology aims for a holistic description and understanding of biological processes by an integration of analytical experimental approaches along with synthetic computational models. In fact, biological networks have been developed as a platform for integrating information from high to low-throughput experiments for the analysis of biological systems. We provide an overview of all processes used in modeling and simulating biological networks in such a way that they can become easily understandable for researchers with both biological and mathematical backgrounds. Consequently, given the complexity of generated experimental data and cellular networks, it is no surprise that researchers have turned to computer simulation and the development of more theory-based approaches to augment and assist in the development of a fully quantitative understanding of cellular dynamics.

Coherent functional modules improve transcription factor target identification, cooperativity prediction, and disease association

Transcription factors (TFs) are fundamental controllers of cellular regulation that function in a complex and combinatorial manner. Accurate identification of a transcription factor's targets is essential to understanding the role that factors play in disease biology. However, due to a high false positive rate, identifying coherent functional target sets is difficult. We have created an improved mapping of targets by integrating ChIP-Seq data with 423 functional modules derived from 9,395 human expression experiments. We identified 5,002 TF-module relationships, significantly improved TF target prediction, and found 30 high-confidence TF-TF associations, of which 14 are known. Importantly, we also connected TFs to diseases through these functional modules and identified 3,859 significant TF-disease relationships. As an example, we found a link between MEF2A and Crohn's disease, which we validated in an independent expression dataset. These results show the power of combining expression data and ChIP-Seq data to remove noise and better extract the associations between TFs, functional modules, and disease.

Transcription factors (TFs) are crucial to the precise regulation of many cellular processes and thus, are responsible for many human phenotypes and diseases. Now that the ENCODE project has mapped hundreds of TFs to their genomic binding locations, extracting functional biological signals is the next step in understanding their role in disease. In this paper, we present a novel approach to identifying TF targets and use these targets to find regulatory relationships between TFs and diseases. We present a large open dataset of putative TF-TF interactions and TF-disease associations which includes known connections as well as novel ones. We validate the association of one of our novel TF-disease associations, MEF2A and Crohn's disease, suggesting that our approach generates testable disease association hypotheses. Integrating these datasets will be crucial for understanding phenotypes and complex diseases.

Signaling networks in Leishmania macrophages deciphered through integrated systems biology: a mathematical modeling approach

Network of signaling proteins and functional interaction between the infected cell and the leishmanial parasite, though are not well understood, may be deciphered computationally by reconstructing the immune signaling network. As we all know signaling pathways are well-known abstractions that explain the mechanisms whereby cells respond to signals, collections of pathways form networks, and interactions between pathways in a network, known as cross-talk, enables further complex signaling behaviours. In silico perturbations can help identify sensitive crosstalk points in the network which can be pharmacologically tested. In this study, we have developed a model for immune signaling cascade in leishmaniasis and based upon the interaction analysis obtained through simulation, we have developed a model network, between four signaling pathways i.e., CD14, epidermal growth factor (EGF), tumor necrotic factor (TNF) and PI3 K mediated signaling. Principal component analysis of the signaling network showed that EGF and TNF pathways can be potent pharmacological targets to curb leishmaniasis. The approach is illustrated with a proposed workable model of epidermal growth factor receptor (EGFR) that modulates the immune response. EGFR signaling represents a critical junction between inflammation related signal and potent cell regulation machinery that modulates the expression of cytokines.

Transfer functions for protein signal transduction: application to a model of striatal neural plasticity

We present a novel formulation for biochemical reaction networks in the context of protein signal transduction. The model consists of input-output transfer functions, which are derived from differential equations, using stable equilibria. We select a set of “source” species, which are interpreted as input signals. Signals are transmitted to all other species in the system (the “target” species) with a specific delay and with a specific transmission strength. The delay is computed as the maximal reaction time until a stable equilibrium for the target species is reached, in the context of all other reactions in the system. The transmission strength is the concentration change of the target species. The computed input-output transfer functions can be stored in a matrix, fitted with parameters, and even recalled to build dynamical models on the basis of state changes. By separating the temporal and the magnitudinal domain we can greatly simplify the computational model, circumventing typical problems of complex dynamical systems. The transfer function transformation of biochemical reaction systems can be applied to mass-action kinetic models of signal transduction. The paper shows that this approach yields significant novel insights while remaining a fully testable and executable dynamical model for signal transduction. In particular we can deconstruct the complex system into local transfer functions between individual species. As an example, we examine modularity and signal integration using a published model of striatal neural plasticity. The modularizations that emerge correspond to a known biological distinction between calcium-dependent and cAMP-dependent pathways. Remarkably, we found that overall interconnectedness depends on the magnitude of inputs, with higher connectivity at low input concentrations and significant modularization at moderate to high input concentrations. This general result, which directly follows from the properties of individual transfer functions, contradicts notions of ubiquitous complexity by showing input-dependent signal transmission inactivation.

Signaling network prediction by the Ontology Fingerprint enhanced Bayesian network

BACKGROUND

Despite large amounts of available genomic and proteomic data, predicting the structure and response of signaling networks is still a significant challenge. While statistical method such as Bayesian network has been explored to meet this challenge, employing existing biological knowledge for network prediction is difficult. The objective of this study is to develop a novel approach that integrates prior biological knowledge in the form of the Ontology Fingerprint to infer cell-type-specific signaling networks via data-driven Bayesian network learning; and to further use the trained model to predict cellular responses.

RESULTS

We applied our novel approach to address the Predictive Signaling Network Modeling challenge of the fourth (2009) Dialog for Reverse Engineering Assessment's and Methods (DREAM4) competition. The challenge results showed that our method accurately captured signal transduction of a network of protein kinases and phosphoproteins in that the predicted protein phosphorylation levels under all experimental conditions were highly correlated (R2 = 0.93) with the observed results. Based on the evaluation of the DREAM4 organizer, our team was ranked as one of the top five best performers in predicting network structure and protein phosphorylation activity under test conditions.

CONCLUSIONS

Bayesian network can be used to simulate the propagation of signals in cellular systems. Incorporating the Ontology Fingerprint as prior biological knowledge allows us to efficiently infer concise signaling network structure and to accurately predict cellular responses.

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.

Integrated experimental and model-based analysis reveals the spatial aspects of EGFR activation dynamics

The epidermal growth factor receptor (EGFR) belongs to the ErbB family of receptor tyrosine kinases, and controls a diverse set of cellular responses relevant to development and tumorigenesis. ErbB activation is a complex process involving receptor-ligand binding, receptor dimerization, phosphorylation, and trafficking (internalization, recycling and degradation), which together dictate the spatio-temporal distribution of active receptors within the cell. The ability to predict this distribution, and elucidation of the factors regulating it, would help to establish a mechanistic link between ErbB expression levels and the cellular response. Towards this end, we constructed mathematical models to determine the contributions of receptor dimerization and phosphorylation to EGFR activation, and to examine the dependence of these processes on sub-cellular location. We collected experimental datasets for EGFR activation dynamics in human mammary epithelial cells, with the specific goal of model parameterization, and used the data to estimate parameters for several alternate models. Model-based analysis indicated that: (1) signal termination via receptor dephosphorylation in late endosomes, prior to degradation, is an important component of the response, (2) less than 40% of the receptors in the cell are phosphorylated at any given time, even at saturating ligand doses, and (3) receptor phosphorylation kinetics at the cell surface and early endosomes are comparable. We validated the last finding by measuring the EGFR dephosphorylation rates at various times following ligand addition both in whole cells and in endosomes using ELISAs and fluorescent imaging. Overall, our results provide important information on how EGFR phosphorylation levels are regulated within cells. This study demonstrates that an iterative cycle of experiments and modeling can be used to gain mechanistic insight regarding complex cell signaling networks.

Static network structure can be used to model the phenotypic effects of perturbations in regulatory networks

MOTIVATION

Biological processes are dynamic, whereas the networks that depict them are typically static. Quantitative modeling using differential equations or logic-based functions can offer quantitative predictions of the behavior of biological systems, but they require detailed experimental characterization of interaction kinetics, which is typically unavailable. To determine to what extent complex biological processes can be modeled and analyzed using only the static structure of the network (i.e. the direction and sign of the edges), we attempt to predict the phenotypic effect of perturbations in biological networks from the static network structure.

RESULTS

We analyzed three networks from different sources: The EGFR/MAPK and PI3K/AKT network from a detailed experimental study, the TNF regulatory network from the STRING database and a large network of all NCI-curated pathways from the Protein Interaction Database. Altogether, we predicted the effect of 39 perturbations (e.g. by one or two drugs) on 433 target proteins/genes. In up to 82% of the cases, an algorithm that used only the static structure of the network correctly predicted whether any given protein/gene is upregulated or downregulated as a result of perturbations of other proteins/genes.

CONCLUSION

While quantitative modeling requires detailed experimental data and heavy computations, which limit its scalability for large networks, a wiring-based approach can use available data from pathway and interaction databases and may be scalable. These results lay the foundations for a large-scale approach of predicting phenotypes based on the schematic structure of networks.

Combinatorial drug therapy for cancer in the post-genomic era

Over the past decade, whole genome sequencing and other 'omics' technologies have defined pathogenic driver mutations to which tumor cells are addicted. Such addictions, synthetic lethalities and other tumor vulnerabilities have yielded novel targets for a new generation of cancer drugs to treat discrete, genetically defined patient subgroups. This personalized cancer medicine strategy could eventually replace the conventional one-size-fits-all cytotoxic chemotherapy approach. However, the extraordinary intratumor genetic heterogeneity in cancers revealed by deep sequencing explains why de novo and acquired resistance arise with molecularly targeted drugs and cytotoxic chemotherapy, limiting their utility. One solution to the enduring challenge of polygenic cancer drug resistance is rational combinatorial targeted therapy.

A quick guide to biomolecular network studies: construction, analysis, applications, and resources

Over the past decade, a rapid increase in network data including signaling, transcription regulation, metabolic reaction, protein-protein interaction and genetic interaction has been observed. Many biology issues have been investigated by analyzing these diverse networks, providing new insights into biology. Networks also play an important role in disease studies including disease gene screening and clinical diagnosis. Large amounts of databases and software have been developed to facilitate the storage, exchange, integration, and analysis of network data and network analysis is becoming a routine procedure for biologists to infer biological information. In this review, several main aspects of network studies are discussed, including network construction, analysis, application, and resources.

Computational modelling of maternal interactions with spermatozoa: potentials and prospects

Understanding the complex interactions between gametes, embryos and the maternal tract is required knowledge for combating infertility and developing new methods of contraception. Here we present some main aspects of spermatozoa interactions with the mammalian oviduct before fertilisation and discuss how computational modelling can be used as an invaluable aid to experimental investigation in this field. A complete predictive computational model of gamete and embryo interactions with the female reproductive tract is a long way off. However, the enormity of this task should not discourage us from working towards it. Computational modelling allows us to investigate aspects of maternal communication with gametes and embryos, which are financially, ethically or practically difficult to look at experimentally. In silico models of maternal communication with gametes and embryos can be used as tools to complement in vivo experiments, in the same way as in vitro and in situ models.

Construction of large signaling pathways using an adaptive perturbation approach with phosphoproteomic data

Construction of large and cell-specific signaling pathways is essential to understand information processing under normal and pathological conditions. On this front, gene-based approaches offer the advantage of large pathway exploration whereas phosphoproteomic approaches offer a more reliable view of pathway activities but are applicable to small pathway sizes. In this paper, we demonstrate an experimentally adaptive approach to construct large signaling pathways from phosphoproteomic data within a 3-day time frame. Our approach--taking advantage of the fast turnaround time of the xMAP technology--is carried out in four steps: (i) screen optimal pathway inducers, (ii) select the responsive ones, (iii) combine them in a combinatorial fashion to construct a phosphoproteomic dataset, and (iv) optimize a reduced generic pathway via an Integer Linear Programming formulation. As a case study, we uncover novel players and their corresponding pathways in primary human hepatocytes by interrogating the signal transduction downstream of 81 receptors of interest and constructing a detailed model for the responsive part of the network comprising 177 species (of which 14 are measured) and 365 interactions.

Systems vaccinology: learning to compute the behavior of vaccine induced immunity

The goal of systems biology is to access and integrate information about the parts (e.g., genes, proteins, cells) of a biological system with a view to computing and predicting the behavior of the system. The past decade has witnessed technological revolutions in the capacity to make high throughput measurements about the behavior of genes, proteins, and cells. Such technologies are widely used in biological research and in medicine, such as toward prognosis and therapy response prediction in cancer patients. More recently, systems biology is being applied to vaccinology, with the goal of: (1) understanding the mechanisms by which vaccines stimulate protective immunity, and (2) predicting the immunogenicity or efficacy of vaccines. Here, we review the recent advances in this area, and highlight the biological and computational challenges posed.

Mechanistic modeling to investigate signaling by oncogenic Ras mutants

Mathematical models based on biochemical reaction mechanisms can be a powerful complement to experimental investigations of cell signaling networks. In principle, such models have the potential to find the behaviors that result from well-understood component interactions and their measurable properties, such as concentrations and rate constants. As cancer results from the acquisition of mutations that alter the expression level and/or the biochemistry of proteins encoded by mutated genes, mathematical models of cell signaling networks would also seem to have the potential to predict how these changes alter cell signaling to produce a cancer phenotype. Ras is commonly found in cancer and has been extensively characterized at the level of detail needed to develop such models. Here, we consider how biochemical mechanism-based models have been used to study mutant Ras signaling. These models demonstrate that it is clearly possible to use observable properties of individual reactions to predict how the entire system behaves to produce the high levels of signal that drive the cancer phenotype. These models also demonstrate differences in how models are developed and studied. Their evaluation suggests which approaches are most promising for future work.

The dynamic state of protein turnover: It's about time

The continual destruction and renewal of proteins that maintain cellular homeostasis has been rigorously studied since the late 1930s. Experimental techniques for measuring protein turnover have evolved to measure the dynamic regulation of key proteins and now, entire proteomes. In the past decade, the proteomics field has aimed to discover how cells adjust their proteomes to execute numerous regulatory programs in response to specific cellular and environmental cues. By combining classical biochemical techniques with modern, high-throughput technologies, researchers have begun to reveal the synthesis and degradation mechanisms that shape protein turnover on a global scale. This review examines several recent developments in protein turnover research, emphasizing the combination of metabolic labeling and mass spectrometry.

Systems biology approaches to dissect mammalian innate immunity

Advances in experimental tools have allowed for the systematic identification of components and biological processes as well as quantification of their activities over time. Together with computational analysis, these measurement and perturbation technologies have given rise to the field of systems biology, which seeks to discover, analyze and model the interactions of physical components in a biological system. Although in its infancy, recent application of this approach has resulted in novel insights into the machinery that regulates and modifies innate immune cell functions. Here, we summarize contributions that have been made through the unbiased interrogation of the mammalian innate immune system, emphasizing the importance of integrating orthogonal datasets into models. To enable application of approaches more broadly, however, a concerted effort across the immunology community to develop reagent and tool platforms will be required.

Systems biology and medicine--meta-issues and frameworks

For further resources related to this article, please visit the WIREs website.