Computational modelling of cancerous mutations in the EGFR/ERK signalling pathway

A Phosphoproteomics Data Resource for Systems-level Modeling of Kinase Signaling Networks

Building mechanistic models of kinase-driven signaling pathways requires quantitative measurements of protein phosphorylation across physiologically relevant conditions, but this is rarely done because of the insensitivity of traditional technologies. By using a multiplexed deep phosphoproteome profiling workflow, we were able to generate a deep phosphoproteomics dataset of the EGFR-MAPK pathway in non-transformed MCF10A cells across physiological ligand concentrations with a time resolution of <12 min and in the presence and absence of multiple kinase inhibitors. An improved phosphosite mapping technique allowed us to reliably identify >46,000 phosphorylation sites on >6600 proteins, of which >4500 sites from 2110 proteins displayed a >2-fold increase in phosphorylation in response to EGF. This data was then placed into a cellular context by linking it to 15 previously published protein databases. We found that our results were consistent with much, but not all previously reported data regarding the activation and negative feedback phosphorylation of core EGFR-ERK pathway proteins. We also found that EGFR signaling is biphasic with substrates downstream of RAS/MAPK activation showing a maximum response at <3ng/ml EGF while direct substrates, such as HGS and STAT5B, showing no saturation. We found that RAS activation is mediated by at least 3 parallel pathways, two of which depend on PTPN11. There appears to be an approximately 4-minute delay in pathway activation at the step between RAS and RAF, but subsequent pathway phosphorylation was extremely rapid. Approximately 80 proteins showed a >2-fold increase in phosphorylation across all experiments and these proteins had a significantly higher median number of phosphorylation sites (~18) relative to total cellular phosphoproteins (~4). Over 60% of EGF-stimulated phosphoproteins were downstream of MAPK and included mediators of cellular processes such as gene transcription, transport, signal transduction and cytoskeletal arrangement. Their phosphorylation was either linear with respect to MAPK activation or biphasic, corresponding to the biphasic signaling seen at the level of the EGFR. This deep, integrated phosphoproteomics data resource should be useful in building mechanistic models of EGFR and MAPK signaling and for understanding how downstream responses are regulated.

NF-kB in Signaling Patterns and Its Temporal Dynamics Encode/Decode Human Diseases

Defects in signaling pathways are the root cause of many disorders. These malformations come in a wide variety of types, and their causes are also very diverse. Some of these flaws can be brought on by pathogenic organisms and viruses, many of which can obstruct signaling processes. Other illnesses are linked to malfunctions in the way that cell signaling pathways work. When thinking about how errors in signaling pathways might cause disease, the idea of signalosome remodeling is helpful. The signalosome may be conveniently divided into two types of defects: phenotypic remodeling and genotypic remodeling. The majority of significant illnesses that affect people, including high blood pressure, heart disease, diabetes, and many types of mental illness, appear to be caused by minute phenotypic changes in signaling pathways. Such phenotypic remodeling modifies cell behavior and subverts normal cellular processes, resulting in illness. There has not been much progress in creating efficient therapies since it has been challenging to definitively confirm this connection between signalosome remodeling and illness. The considerable redundancy included into cell signaling systems presents several potential for developing novel treatments for various disease conditions. One of the most important pathways, NF-κB, controls several aspects of innate and adaptive immune responses, is a key modulator of inflammatory reactions, and has been widely studied both from experimental and theoretical perspectives. NF-κB contributes to the control of inflammasomes and stimulates the expression of a number of pro-inflammatory genes, including those that produce cytokines and chemokines. Additionally, NF-κB is essential for controlling innate immune cells and inflammatory T cells' survival, activation, and differentiation. As a result, aberrant NF-κB activation plays a role in the pathogenesis of several inflammatory illnesses. The activation and function of NF-κB in relation to inflammatory illnesses was covered here, and the advancement of treatment approaches based on NF-κB inhibition will be highlighted. This review presents the temporal behavior of NF-κB and its potential relevance in different human diseases which will be helpful not only for theoretical but also for experimental perspectives.

Navigating Multi-Scale Cancer Systems Biology Towards Model-Driven Clinical Oncology and Its Applications in Personalized Therapeutics

Rapid advancements in high-throughput omics technologies and experimental protocols have led to the generation of vast amounts of scale-specific biomolecular data on cancer that now populates several online databases and resources. Cancer systems biology models built using this data have the potential to provide specific insights into complex multifactorial aberrations underpinning tumor initiation, development, and metastasis. Furthermore, the annotation of these single- and multi-scale models with patient data can additionally assist in designing personalized therapeutic interventions as well as aid in clinical decision-making. Here, we have systematically reviewed the emergence and evolution of (i) repositories with scale-specific and multi-scale biomolecular cancer data, (ii) systems biology models developed using this data, (iii) associated simulation software for the development of personalized cancer therapeutics, and (iv) translational attempts to pipeline multi-scale panomics data for data-driven in silico clinical oncology. The review concludes that the absence of a generic, zero-code, panomics-based multi-scale modeling pipeline and associated software framework, impedes the development and seamless deployment of personalized in silico multi-scale models in clinical settings.

MicroRNA‑137 regulates hypoxia‑mediated migration and epithelial‑mesenchymal transition in prostate cancer by targeting LGR4 via the EGFR/ERK signaling pathway

MicroRNAs (miRs) serve an integral role in prostate cancer. The present study aimed to investigate the effects and mechanisms of miR‑137 in hypoxia‑mediated migration and epithelial‑mesenchymal transition (EMT). PC3 and DU145 prostate cancer cells were exposed to hypoxia for 24 h, after which the expression of miR‑137 was determined by reverse transcription‑quantitative PCR (RT‑qPCR). The cells were transfected with a miR‑137 mimic or inhibitor, followed by hypoxia exposure. The results demonstrated that hypoxia reduced miR‑137 expression. Further results from the Cell Counting Kit‑8, Cell Death Detection ELISA plus kit, Transwell assay, RT‑qPCR and western blotting assays revealed that the miR‑137 mimic prevented cell proliferation, facilitated apoptosis and repressed cell migration, invasiveness, and expression of N‑cadherin, vimentin and matrix metalloproteinase 2; the miR‑137 inhibitor exerted the opposite effects. A dual‑-luciferase reporter assay determined that miR‑137 directly targeted leucine‑rich repeat‑containing G protein‑coupled receptor 4 (LGR4). Additionally, miR‑137 negatively regulated the epidermal growth factor receptor/extracellular signal‑-regulated kinase (EGFR/ERK) signaling pathway by targeting LGR4. LGR4 silencing or EGFR/ERK inhibition abolished the effects of miR‑137 inhibitor on cell migration and EMT. In conclusion, by targeting LGR4 via the EGFR/ERK signaling pathway, miR‑137 inhibited prostate cancer cell migration and EMT.

Applications of personalised signalling network models in precision oncology

As our ability to provide in-depth, patient-specific characterisation of the molecular alterations within tumours rapidly improves, it is becoming apparent that new approaches will be required to leverage the power of this data and derive the full benefit for each individual patient. Systems biology approaches are beginning to emerge within this field as a potential method of incorporating large volumes of network level data and distilling a coherent, clinically-relevant prediction of drug response. However, the initial promise of this developing field is yet to be realised. Here we argue that in order to develop these precise models of individual drug response and tailor treatment accordingly, we will need to develop mathematical models capable of capturing both the dynamic nature of drug-response signalling networks and key patient-specific information such as mutation status or expression profiles. We also review the modelling approaches commonly utilised within this field, and outline recent examples of their use in furthering the application of systems biology for a precision medicine approach to cancer treatment.

The road to ERK activation: Do neurons take alternate routes?

The ERK cascade is a central signaling pathway that regulates a wide variety of cellular processes including proliferation, differentiation, learning and memory, development, and synaptic plasticity. A wide range of inputs travel from the membrane through different signaling pathway routes to reach activation of one set of output kinases, ERK1&2. The classical ERK activation pathway beings with growth factor activation of receptor tyrosine kinases. Numerous G-protein coupled receptors and ionotropic receptors also lead to ERK through increases in the second messengers calcium and cAMP. Though both types of pathways are present in diverse cell types, a key difference is that most stimuli to neurons, e.g. synaptic inputs, are transient, on the order of milliseconds to seconds, whereas many stimuli acting on non-neural tissue, e.g. growth factors, are longer duration. The ability to consolidate these inputs to regulate the activation of ERK in response to diverse signals raises the question of which factors influence the difference in ERK activation pathways. This review presents both experimental studies and computational models aimed at understanding the control of ERK activation and whether there are fundamental differences between neurons and other cells. Our main conclusion is that differences between cell types are quite subtle, often related to differences in expression pattern and quantity of some molecules such as Raf isoforms. In addition, the spatial location of ERK is critical, with regulation by scaffolding proteins producing differences due to colocalization of upstream molecules that may differ between neurons and other cells.

Comparative Network Reconstruction using mixed integer programming

Motivation

Signal-transduction networks are often aberrated in cancer cells, and new anti-cancer drugs that specifically target oncogenes involved in signaling show great clinical promise. However, the effectiveness of such targeted treatments is often hampered by innate or acquired resistance due to feedbacks, crosstalks or network adaptations in response to drug treatment. A quantitative understanding of these signaling networks and how they differ between cells with different oncogenic mutations or between sensitive and resistant cells can help in addressing this problem.

Results

Here, we present Comparative Network Reconstruction (CNR), a computational method to reconstruct signaling networks based on possibly incomplete perturbation data, and to identify which edges differ quantitatively between two or more signaling networks. Prior knowledge about network topology is not required but can straightforwardly be incorporated. We extensively tested our approach using simulated data and applied it to perturbation data from a BRAF mutant, PTPN11 KO cell line that developed resistance to BRAF inhibition. Comparing the reconstructed networks of sensitive and resistant cells suggests that the resistance mechanism involves re-establishing wild-type MAPK signaling, possibly through an alternative RAF-isoform.

Availability and implementation

CNR is available as a python module at https://github.com/NKI-CCB/cnr. Additionally, code to reproduce all figures is available at https://github.com/NKI-CCB/CNR-analyses.

Supplementary information

Supplementary data are available at Bioinformatics online.

Estrogen receptors promote NSCLC progression by modulating the membrane receptor signaling network: a systems biology perspective

Background

Estrogen receptors (ERs) are thought to play an important role in non-small cell lung cancer (NSCLC). However, the effect of ERs in NSCLC is still controversial and needs further investigation. A new consideration is that ERs may affect NSCLC progression through complicated molecular signaling networks rather than individual targets. Therefore, this study aims to explore the effect of ERs in NSCLC from the perspective of cancer systems biology.

Methods

The gene expression profile of NSCLC samples in TCGA dataset was analyzed by bioinformatics method. Variations of cell behaviors and protein expression were detected in vitro. The kinetic process of molecular signaling network was illustrated by a systemic computational model. At last, immunohistochemical (IHC) and survival analysis was applied to evaluate the clinical relevance and prognostic effect of key receptors in NSCLC.

Results

Bioinformatics analysis revealed that ERs might affect many cancer-related molecular events and pathways in NSCLC, particularly membrane receptor activation and signal transduction, which might ultimately lead to changes in cell behaviors. Experimental results confirmed that ERs could regulate cell behaviors including cell proliferation, apoptosis, invasion and migration; ERs also regulated the expression or activation of key members in membrane receptor signaling pathways such as epidermal growth factor receptor (EGFR), Notch1 and Glycogen synthase kinase-3β/β-Catenin (GSK3β/β-Catenin) pathways. Modeling results illustrated that the promotive effect of ERs in NSCLC was implemented by modulating the signaling network composed of EGFR, Notch1 and GSK3β/β-Catenin pathways; ERs maintained and enhanced the output of oncogenic signals by adding redundant and positive-feedback paths into the network. IHC results echoed that high expression of ERs, EGFR and Notch1 had a synergistic effect on poor prognosis of advanced NSCLC.

Conclusions

This study indicated that ERs were likely to promote NSCLC progression by modulating the integrated membrane receptor signaling network composed of EGFR, Notch1 and GSK3β/β-Catenin pathways and then affecting tumor cell behaviors. It also complemented the molecular mechanisms underlying the progression of NSCLC and provided new opportunities for optimizing therapeutic scheme of NSCLC.

A 3D Multiscale Model to Explore the Role of EGFR Overexpression in Tumourigenesis

The epidermal growth factor receptor (EGFR) signalling cascade is one of the main pathways that regulate the survival and division of mammalian cells. It is also one of the most altered transduction pathways in cancer. Acquired mutations in the EGFR/ERK pathway can cause the overexpression of EGFR on the surface of the cell, while others downregulate the inactivation of switched on intracellular proteins such as Ras and Raf. This upregulates the activity of ERK and promotes cell division. We develop a 3D multiscale model to explore the role of EGFR overexpression on tumour initiation. In this model, cells are described as individual objects that move, interact, divide, proliferate, and die by apoptosis. We use Brownian Dynamics to describe the extracellular and intracellular regulations of cells as well as the spatial and stochastic effects influencing them. The fate of each cell depends on the number of active transcription factors in the nucleus. We use numerical simulations to investigate the individual and combined effects of mutations on the intracellular regulation of individual cells. Next, we show that the distance between active receptors increase the level of EGFR/ERK signalling. We demonstrate the usefulness of the model by quantifying the impact of mutational alterations in the EGFR/ERK pathway on the growth rate of in silico tumours.

Ribophorin II potentiates P-glycoprotein- and ABCG2-mediated multidrug resistance via activating ERK pathway in gastric cancer

Multidrug resistance (MDR) is a critical reason of cancer chemotherapy failure. Ribophorin II (RPN2) has emerged as a vital regulator of MDR in multiple cancers including gastric cancer (GC). However, the roles and molecular mechanisms of RPN2 in MDR have not been well featured till now. The present study aimed to explore the roles and molecular mechanisms of RPN2 in MDR of drug-resistant GC cells. Results showed that the expressions of RPN2, multidrug resistance 1 (MDR1), and ATP binding cassette subfamily G member 2 (ABCG2) were upregulated in SGC7901/DDP and SGC7901/VCR cells. Knockdown of RPN2 alleviated MDR through downregulating MDR1 and ABCG2 expressions in SGC7901/DDP and SGC7901/VCR cells. RPN2 depletion inhibited the activation of MEK/ERK pathway. RPN2 overexpression enhanced MDR by upregulating P-glycoprotein (P-gp) and ABCG2 protein expressions in SGC7901/DDP or SGC7901/VCR cells, while this effect of RPN2 was abrogated by ERK knockdown or treatment with ERK inhibitor PD98059. Our findings suggested that RPN2 potentiated P-gp- and ABCG2-mediated MDR via activating MEK/ERK pathway in GC, hinting the critical values of RPN2 in ameliorating MDR and providing a promising target for GC therapy.

The extracellular matrix as a key regulator of intracellular signalling networks

The extracellular matrix (ECM) is a salient feature of all solid tissues within the body. This complex, acellular entity is composed of hundreds of individual molecules whose assembly, architecture and biomechanical properties are critical to controlling the behaviour and phenotype of the different cell types residing within tissues. Cells are the basic unit of life and the core building block of tissues and organs. At their simplest, they follow a set of rules, governed by their genetic code and effected through the complex protein signalling networks that these genes encode. These signalling networks assimilate and process the information received by the cell to control cellular decisions that govern cell fate. The ECM is the biggest provider of external stimuli to cells and as such is responsible for influencing intracellular signalling dynamics. In this review, we discuss the inclusion of ECM as a central regulatory signalling sub-network in computational models of cellular decision making, with a focus on its role in diseases such as cancer. LINKED ARTICLES: This article is part of a themed section on Translating the Matrix. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v176.1/issuetoc.

New insights into RAS biology reinvigorate interest in mathematical modeling of RAS signaling

RAS is the most frequently mutated gene across human cancers, but developing inhibitors of mutant RAS has proven to be challenging. Given the difficulties of targeting RAS directly, drugs that impact the other components of pathways where mutant RAS operates may potentially be effective. However, the system-level features, including different localizations of RAS isoforms, competition between downstream effectors, and interlocking feedback and feed-forward loops, must be understood to fully grasp the opportunities and limitations of inhibiting specific targets. Mathematical modeling can help us discern the system-level impacts of these features in normal and cancer cells. New technologies enable the acquisition of experimental data that will facilitate development of realistic models of oncogenic RAS behavior. In light of the wealth of empirical data accumulated over decades of study and the advancement of experimental methods for gathering new data, modelers now have the opportunity to advance progress toward realization of targeted treatment for mutant RAS-driven cancers.

Cancerous perturbations within the ERK, PI3K/Akt, and Wnt/β-catenin signaling network constitutively activate inter-pathway positive feedback loops

Perturbations in molecular signaling pathways are a result of genetic or epigenetic alterations, which may lead to malignant transformation of cells. Despite cellular robustness, specific genetic or epigenetic changes of any gene can trigger a cascade of failures, which result in the malfunctioning of cell signaling pathways and lead to cancer phenotypes. The extent of cellular robustness has a link with the architecture of the network such as feedback and feedforward loops. Perturbation in components within feedback loops causes a transition from a regulated to a persistently activated state and results in uncontrolled cell growth. This work represents the mathematical and quantitative modeling of ERK, PI3K/Akt, and Wnt/β-catenin signaling crosstalk to show the dynamics of signaling responses during genetic and epigenetic changes in cancer. ERK, PI3K/Akt, and Wnt/β-catenin signaling crosstalk networks include both intra and inter-pathway feedback loops which function in a controlled fashion in a healthy cell. Our results show that cancerous perturbations of components such as EGFR, Ras, B-Raf, PTEN, and components of the destruction complex cause extreme fragility in the network and constitutively activate inter-pathway positive feedback loops. We observed that the aberrant signaling response due to the failure of specific network components is transmitted throughout the network via crosstalk, generating an additive effect on cancer growth and proliferation.

Activated Oncogenic Pathway Modifies Iron Network in Breast Epithelial Cells: A Dynamic Modeling Perspective

Dysregulation of iron metabolism in cancer is well documented and it has been suggested that there is interdependence between excess iron and increased cancer incidence and progression. In an effort to better understand the linkages between iron metabolism and breast cancer, a predictive mathematical model of an expanded iron homeostasis pathway was constructed that includes species involved in iron utilization, oxidative stress response and oncogenic pathways. The model leads to three predictions. The first is that overexpression of iron regulatory protein 2 (IRP2) recapitulates many aspects of the alterations in free iron and iron-related proteins in cancer cells without affecting the oxidative stress response or the oncogenic pathways included in the model. This prediction was validated by experimentation. The second prediction is that iron-related proteins are dramatically affected by mitochondrial ferritin overexpression. This prediction was validated by results in the pertinent literature not used for model construction. The third prediction is that oncogenic Ras pathways contribute to altered iron homeostasis in cancer cells. This prediction was validated by a combination of simulation experiments of Ras overexpression and catalase knockout in conjunction with the literature. The model successfully captures key aspects of iron metabolism in breast cancer cells and provides a framework upon which more detailed models can be built.

Iron is required for cellular metabolism and growth, but can be toxic due to its ability to cause high oxidative stress and consequently DNA damage. To prevent damage, all organisms that require iron have developed mechanisms to tightly control iron levels. Dysregulation of iron metabolism is detrimental and can contribute to a wide range of diseases, including cancer. This paper presents a predictive mathematical model of iron regulation linked to iron utilization, oxidative stress, and the oncogenic response specific to normal breast epithelial cells. The model uses a discrete modeling framework to generate novel biological hypotheses for an investigation of how normal breast cells become malignant cells, capturing a breast cancer phenotype of iron homeostasis through overexpression and knockout simulations. The new biology discovered is (1) IRP2 overexpression alters the iron homeostasis pathway in breast cells, without affecting the oxidative stress response or oncogenic pathways, (2) an activated oncogenic pathway disrupts iron regulation in breast cancer cells.

Probing the Heterogeneity of Protein Kinase Activation in Cells by Super-resolution Microscopy

Heterogeneity of mitogen-activated protein kinase (MAPK) activation in genetically identical cells, which occurs in response to epidermal growth factor receptor (EGFR) signaling, remains poorly understood. MAPK cascades integrate signals emanating from different EGFR spatial locations, including the plasma membrane and endocytic compartment. We previously hypothesized that in EGF-stimulated cells the MAPK phosphorylation (pMAPK) level and activity are largely determined by the spatial organization of the EGFR clusters within the cell. For experimental testing of this hypothesis, we used super-resolution microscopy to define EGFR clusters by receptor numbers (N) and average intracluster distances (d). From these data, we predicted the extent of pMAPK with 85% accuracy on a cell-to-cell basis with control data returning 54% accuracy (P < 0.001). For comparison, the prediction accuracy was only 61% (P = 0.382) when the diffraction-limited averaged fluorescence intensity/cluster was used. Large clusters (N ≥ 3) with d > 50 nm were most predictive for pMAPK level in cells. Electron microscopy revealed that these large clusters were primarily localized to the limiting membrane of multivesicular bodies (MVB). Many tighter packed dimers/multimers (d < 50 nm) were found on intraluminal vesicles within MVBs, where they were unlikely to activate MAPK because of the physical separation. Our results suggest that cell-to-cell differences in N and d contain crucial information to predict EGFR-activated cellular pMAPK levels and explain pMAPK heterogeneity in isogenic cells.

Integrating network reconstruction with mechanistic modeling to predict cancer therapies

Signal transduction networks are often rewired in cancer cells. Identifying these alterations will enable more effective cancer treatment. We developed a computational framework that can identify, reconstruct, and mechanistically model these rewired networks from noisy and incomplete perturbation response data and then predict potential targets for intervention. As a proof of principle, we analyzed a perturbation data set targeting epidermal growth factor receptor (EGFR) and insulin-like growth factor 1 receptor (IGF1R) pathways in a panel of colorectal cancer cells. Our computational approach predicted cell line-specific network rewiring. In particular, feedback inhibition of insulin receptor substrate 1 (IRS1) by the kinase p70S6K was predicted to confer resistance to EGFR inhibition, suggesting that disrupting this feedback may restore sensitivity to EGFR inhibitors in colorectal cancer cells. We experimentally validated this prediction with colorectal cancer cell lines in culture and in a zebrafish (Danio rerio) xenograft model.

A rule-based model of insulin signalling pathway

BACKGROUND

The insulin signalling pathway (ISP) is an important biochemical pathway, which regulates some fundamental biological functions such as glucose and lipid metabolism, protein synthesis, cell proliferation, cell differentiation and apoptosis. In the last years, different mathematical models based on ordinary differential equations have been proposed in the literature to describe specific features of the ISP, thus providing a description of the behaviour of the system and its emerging properties. However, protein-protein interactions potentially generate a multiplicity of distinct chemical species, an issue referred to as "combinatorial complexity", which results in defining a high number of state variables equal to the number of possible protein modifications. This often leads to complex, error prone and difficult to handle model definitions.

RESULTS

In this work, we present a comprehensive model of the ISP, which integrates three models previously available in the literature by using the rule-based modelling (RBM) approach. RBM allows for a simple description of a number of signalling pathway characteristics, such as the phosphorylation of signalling proteins at multiple sites with different effects, the simultaneous interaction of many molecules of the signalling pathways with several binding partners, and the information about subcellular localization where reactions take place. Thanks to its modularity, it also allows an easy integration of different pathways. After RBM specification, we simulated the dynamic behaviour of the ISP model and validated it using experimental data. We the examined the predicted profiles of all the active species and clustered them in four clusters according to their dynamic behaviour. Finally, we used parametric sensitivity analysis to show the role of negative feedback loops in controlling the robustness of the system.

CONCLUSIONS

The presented ISP model is a powerful tool for data simulation and can be used in combination with experimental approaches to guide the experimental design. The model is available at http://sysbiobig.dei.unipd.it/ was submitted to Biomodels Database ( https://www.ebi.ac.uk/biomodels-main/ # MODEL 1604100005).

Network and systems biology: essential steps in virtualising drug discovery and development

The biological processes that keep us healthy or cause disease, as well as the mechanisms of action of possible drugs are inherently complex. In the face of this complexity, attempts at discovering new drugs to treat diseases have alternated between trial-and-error (typically on experimental systems) and grand simplification, usually based on much too little information. We now have the chance to combine these strategies through establishment of 'virtual patient' models, centred on a detailed molecular characterisation of thousands or even, in the future, millions of patients. In doing so, we lay the foundations for truly personalised therapy, as well as a far-reaching virtualisation of drug discovery and development in oncology and other areas of medicine.

Alcohol Disrupts Human Liver Stem/Progenitor Cell Proliferation and Differentiation

OBJECTIVE

Excessive alcohol consumption injures the liver resulting in various liver diseases including liver cirrhosis. Advanced liver disease continues to be a major challenge to human health. Liver stem/progenitor cells (LSPCs) are tissue specific precursors with a distinct capacity of multi-lineage differentiation. These precursor cells may play an important role in the process of tissue injury repair and pathological transition of liver structures. At the present time, knowledge about the effect of alcohol on LSPC function during the development of alcoholic liver disease remains absent. This study was conducted to investigate changes in LSPC activity of proliferation and differentiation following alcohol exposure. The disruption of cell signaling mechanisms underlying alcohol-induced alteration of LSPC activities was also examined.

METHODS

Primary and immortalized human liver stem cells (HL1-1 cells and HL1-hT1 cells, respectively) were cultured in media optimized for cell proliferation and hepatocyte differentiation in the absence and presence of ethanol. Changes in cell morphology, proliferation and differentiation were determined. Functional disruption of cell signaling components following alcohol exposure was examined.

RESULTS

Ethanol exposure suppressed HL1-1 cell growth [as measured by cell 5-bromo-2-deoxyuridine (BrdU) incorporation] mediated by epidermal growth factor (EGF) or EGF plus interleukin-6 (IL-6) in an ethanol dose-dependent manner. Similarly, ethanol inhibited BrdU incorporation into HL1-hT1 cells. Cyclin D1 mRNA expression by HL1-hT1 cells was suppressed when cells were cultured with 50 and 100 mM ethanol. Ethanol exposure induced morphological change of HL1-1 cells toward a myofibroblast-like phenotype. Furthermore, ethanol down-regulated E-cadherin expression while increasing collagen I expression by HL1-1 cells. Ethanol also stimulated Snail transcriptional repressor (Snail) and α-smooth muscle actin (α-SMA) gene expression by HL1-1 cells.

CONCLUSION

These results demonstrate that the direct effect of alcohol on LSPCs is inhibiting their proliferation and promoting mesenchymal transition during their differentiation. Alcohol interrupts LSPC differentiation through interfering Snail signaling.

Leucine modulates dynamic phosphorylation events in insulin signaling pathway and enhances insulin-dependent glycogen synthesis in human skeletal muscle cells

BACKGROUND

Branched-chain amino acids, especially leucine, are known to interact with insulin signaling pathway and glucose metabolism. However, the mechanism by which this is exerted, remain to be clearly defined. In order to examine the effect of leucine on muscle insulin signaling, a set of experiments was carried out to quantitate phosphorylation events along the insulin signaling pathway in human skeletal muscle cell cultures. Cells were exposed to insulin, leucine or both, and phosphorylation events of key insulin signaling molecules were tracked over time so as to monitor time-related responses that characterize the signaling events and could be missed by a single sampling strategy limited to pre/post stimulus events.

RESULTS

Leucine is shown to increase the magnitude of insulin-dependent phosphorylation of protein kinase B (AKT) at Ser473 and glycogen synthase kinase (GSK3β) at Ser21-9. Glycogen synthesis follows the same pattern of GSK3β, with a significant increase at 100 μM leucine plus insulin stimulus. Moreover, data do not show any statistically significant increase of pGSK3β and glycogen synthesis at higher leucine concentrations. Leucine is also shown to increase the magnitude of insulin-mediated extracellularly regulated kinase (ERK) phosphorylation; however, differently from AKT and GSK3β, ERK shows a transient behavior, with an early peak response, followed by a return to the baseline condition.

CONCLUSIONS

These experiments demonstrate a complementary effect of leucine on insulin signaling in a human skeletal muscle cell culture, promoting insulin-activated GSK3β phosphorylation and glycogen synthesis.

Cancer systems biology and modeling: microscopic scale and multiscale approaches

Cancer has become known as a complex and systematic disease on macroscopic, mesoscopic and microscopic scales. Systems biology employs state-of-the-art computational theories and high-throughput experimental data to model and simulate complex biological procedures such as cancer, which involves genetic and epigenetic, in addition to intracellular and extracellular complex interaction networks. In this paper, different systems biology modeling techniques such as systems of differential equations, stochastic methods, Boolean networks, Petri nets, cellular automata methods and agent-based systems are concisely discussed. We have compared the mentioned formalisms and tried to address the span of applicability they can bear on emerging cancer modeling and simulation approaches. Different scales of cancer modeling, namely, microscopic, mesoscopic and macroscopic scales are explained followed by an illustration of angiogenesis in microscopic scale of the cancer modeling. Then, the modeling of cancer cell proliferation and survival are examined on a microscopic scale and the modeling of multiscale tumor growth is explained along with its advantages.

Convergence in parameters and predictions using computational experimental design

Typically, biological models fitted to experimental data suffer from significant parameter uncertainty, which can lead to inaccurate or uncertain predictions. One school of thought holds that accurate estimation of the true parameters of a biological system is inherently problematic. Recent work, however, suggests that optimal experimental design techniques can select sets of experiments whose members probe complementary aspects of a biochemical network that together can account for its full behaviour. Here, we implemented an experimental design approach for selecting sets of experiments that constrain parameter uncertainty. We demonstrated with a model of the epidermal growth factor–nerve growth factor pathway that, after synthetically performing a handful of optimal experiments, the uncertainty in all 48 parameters converged below 10 per cent. Furthermore, the fitted parameters converged to their true values with a small error consistent with the residual uncertainty. When untested experimental conditions were simulated with the fitted models, the predicted species concentrations converged to their true values with errors that were consistent with the residual uncertainty. This paper suggests that accurate parameter estimation is achievable with complementary experiments specifically designed for the task, and that the resulting parametrized models are capable of accurate predictions.

Customizing the therapeutic response of signaling networks to promote antitumor responses by drug combinations

Drug resistance, de novo and acquired, pervades cellular signaling networks (SNs) from one signaling motif to another as a result of cancer progression and/or drug intervention. This resistance is one of the key determinants of efficacy in targeted anti-cancer drug therapy. Although poorly understood, drug resistance is already being addressed in combination therapy by selecting drug targets where SN sensitivity increases due to combination components or as a result of de novo or acquired mutations. Additionally, successive drug combinations have shown low resistance potential. To promote a rational, systematic development of combination therapies, it is necessary to establish the underlying mechanisms that drive the advantages of combination therapies, and design methods to determine drug targets for combination regimens. Based on a joint systems analysis of cellular SN response and its sensitivity to drug action and oncogenic mutations, we describe an in silico method to analyze the targets of drug combinations. Our method explores mechanisms of sensitizing the SN through a combination of two drugs targeting vertical signaling pathways. We propose a paradigm of SN response customization by one drug to both maximize the effect of another drug in combination and promote a robust therapeutic response against oncogenic mutations. The method was applied to customize the response of the ErbB/PI3K/PTEN/AKT pathway by combination of drugs targeting HER2 receptors and proteins in the down-stream pathway. The results of a computational experiment showed that the modification of the SN response from hyperbolic to smooth sigmoid response by manipulation of two drugs in combination leads to greater robustness in therapeutic response against oncogenic mutations determining cancer heterogeneity. The application of this method in drug combination co-development suggests a combined evaluation of inhibition effects together with the capability of drug combinations to suppress resistance mechanisms before they become clinically manifest.

In silico identification of potential therapeutic targets in the TGF-β signal transduction pathway

The transforming growth factor-β (TGF-β) superfamily of cytokines controls fundamental cellular processes, such as proliferation, motility, differentiation, and apoptosis. This fundamental role is emphasized by the widespread presence of mutations of the core components of the TGF-β signal transduction pathway in a number of human diseases. Therefore, there is an increasing interest in the development of therapies to specifically target this pathway. Here we develop a computational approach to identify potential intervention points that are capable of restoring the normal signaling dynamics to the mutated system while maintaining the behavior of normal cells substantially unperturbed. We apply this approach explicitly to the TGF-β pathway to study the signaling dynamics of mutated and normal cells treated with inhibitory drugs and identify the processes in the pathway that are most susceptible to therapeutic intervention.

Integrative modelling of the influence of MAPK network on cancer cell fate decision

The Mitogen-Activated Protein Kinase (MAPK) network consists of tightly interconnected signalling pathways involved in diverse cellular processes, such as cell cycle, survival, apoptosis and differentiation. Although several studies reported the involvement of these signalling cascades in cancer deregulations, the precise mechanisms underlying their influence on the balance between cell proliferation and cell death (cell fate decision) in pathological circumstances remain elusive. Based on an extensive analysis of published data, we have built a comprehensive and generic reaction map for the MAPK signalling network, using CellDesigner software. In order to explore the MAPK responses to different stimuli and better understand their contributions to cell fate decision, we have considered the most crucial components and interactions and encoded them into a logical model, using the software GINsim. Our logical model analysis particularly focuses on urinary bladder cancer, where MAPK network deregulations have often been associated with specific phenotypes. To cope with the combinatorial explosion of the number of states, we have applied novel algorithms for model reduction and for the compression of state transition graphs, both implemented into the software GINsim. The results of systematic simulations for different signal combinations and network perturbations were found globally coherent with published data. In silico experiments further enabled us to delineate the roles of specific components, cross-talks and regulatory feedbacks in cell fate decision. Finally, tentative proliferative or anti-proliferative mechanisms can be connected with established bladder cancer deregulations, namely Epidermal Growth Factor Receptor (EGFR) over-expression and Fibroblast Growth Factor Receptor 3 (FGFR3) activating mutations.

Depending on environmental conditions, strongly intertwined cellular signalling pathways are activated, involving activation/inactivation of proteins and genes in response to external and/or internal stimuli. Alterations of some components of these pathways can lead to wrong cell behaviours. For instance, cancer-related deregulations lead to high proliferation of malignant cells enabling sustained tumour growth. Understanding the precise mechanisms underlying these pathways is necessary to delineate efficient therapeutical approaches for each specific tumour type. We particularly focused on the Mitogen-Activated Protein Kinase (MAPK) signalling network, whose involvement in cancer is well established, although the precise conditions leading to its positive or negative influence on cell proliferation are still poorly understood. We tackled this problem by first collecting sparse published biological information into a comprehensive map describing the MAPK network in terms of stylised chemical reactions. This information source was then used to build a dynamical Boolean model recapitulating network responses to characteristic stimuli observed in selected bladder cancers. Systematic model simulations further allowed us to link specific network components and interactions with proliferative/anti-proliferative cell responses.

Identification of alterations in the Jacobian of biochemical reaction networks from steady state covariance data at two conditions

Model building of biochemical reaction networks typically involves experiments in which changes in the behavior due to natural or experimental perturbations are observed. Computational models of reaction networks are also used in a systems biology approach to study how transitions from a healthy to a diseased state result from changes in genetic or environmental conditions. In this paper we consider the nonlinear inverse problem of inferring information about the Jacobian of a Langevin type network model from covariance data of steady state concentrations associated to two different experimental conditions. Under idealized assumptions on the Langevin fluctuation matrices we prove that relative alterations in the network Jacobian can be uniquely identified when comparing the two data sets. Based on this result and the premise that alteration is locally confined to separable parts due to network modularity we suggest a computational approach using hybrid stochastic-deterministic optimization for the detection of perturbations in the network Jacobian using the sparsity promoting effect of [Formula: see text]-penalization. Our approach is illustrated by means of published metabolomic and signaling reaction networks.

Chemical kinetic mechanistic models to investigate cancer biology and impact cancer medicine

Traditional experimental biology has provided a mechanistic understanding of cancer in which the malignancy develops through the acquisition of mutations that disrupt cellular processes. Several drugs developed to target such mutations have now demonstrated clinical value. These advances are unequivocal testaments to the value of traditional cellular and molecular biology. However, several features of cancer may limit the pace of progress that can be made with established experimental approaches alone. The mutated genes (and resultant mutant proteins) function within large biochemical networks. Biochemical networks typically have a large number of component molecules and are characterized by a large number of quantitative properties. Responses to a stimulus or perturbation are typically nonlinear and can display qualitative changes that depend upon the specific values of variable system properties. Features such as these can complicate the interpretation of experimental data and the formulation of logical hypotheses that drive further research. Mathematical models based upon the molecular reactions that define these networks combined with computational studies have the potential to deal with these obstacles and to enable currently available information to be more completely utilized. Many of the pressing problems in cancer biology and cancer medicine may benefit from a mathematical treatment. As work in this area advances, one can envision a future where such models may meaningfully contribute to the clinical management of cancer patients.

A mathematical model for microRNA in lung cancer

Lung cancer is the leading cause of cancer-related deaths worldwide. Lack of early detection and limited options for targeted therapies are both contributing factors to the dismal statistics observed in lung cancer. Thus, advances in both of these areas are likely to lead to improved outcomes. MicroRNAs (miRs or miRNAs) represent a class of non-coding RNAs that have the capacity for gene regulation and may serve as both diagnostic and prognostic biomarkers in lung cancer. Abnormal expression patterns for several miRNAs have been identified in lung cancers. Specifically, let-7 and miR-9 are deregulated in both lung cancers and other solid malignancies. In this paper, we construct a mathematical model that integrates let-7 and miR-9 expression into a signaling pathway to generate an in silico model for the process of epithelial mesenchymal transition (EMT). Simulations of the model demonstrate that EGFR and Ras mutations in non-small cell lung cancers (NSCLC), which lead to the process of EMT, result in miR-9 upregulation and let-7 suppression, and this process is somewhat robust against random input into miR-9 and more strongly robust against random input into let-7. We elected to validate our model in vitro by testing the effects of EGFR inhibition on downstream MYC, miR-9 and let-7a expression. Interestingly, in an EGFR mutated lung cancer cell line, treatment with an EGFR inhibitor (Gefitinib) resulted in a concentration specific reduction in c-MYC and miR-9 expression while not changing let-7a expression. Our mathematical model explains the signaling link among EGFR, MYC, and miR-9, but not let-7. However, very little is presently known about factors that regulate let-7. It is quite possible that when such regulating factors become known and integrated into our model, they will further support our mathematical model.

Inverse problems from biomedicine: inference of putative disease mechanisms and robust therapeutic strategies

Many complex diseases that are difficult to treat cannot be mapped onto a single cause, but arise from the interplay of multiple contributing factors. In the study of such diseases, it is becoming apparent that therapeutic strategies targeting a single protein or metabolite are often not efficacious. Rather, a systems perspective describing the interaction of physiological components is needed. In this paper, we demonstrate via examples of disease models the kind of inverse problems that arise from the need to infer disease mechanisms and/or therapeutic strategies. We identify the challenges that arise, in particular the need to devise strategies that are robust against variable physiological states and parametric uncertainties.

Reactive oxygen species and epidermal growth factor are antagonistic cues controlling SHP-2 dimerization

The SHP-2 tyrosine phosphatase plays key regulatory roles in the modulation of the cell response to growth factors and cytokines. Over the past decade, the integration of genetic, biochemical, and structural data has helped in interpreting the pathological consequences of altered SHP-2 function. Using complementary approaches, we provide evidence here that endogenous SHP-2 can dimerize through the formation of disulfide bonds that may also involve the catalytic cysteine. We show that the fraction of dimeric SHP-2 is modulated by growth factor stimulation and by the cell redox state. Comparison of the phosphatase activities of the monomeric self-inhibited and dimeric forms indicated that the latter is 3-fold less active, thus pointing to the dimerization process as an additional mechanism for controlling SHP-2 activity. Remarkably, dimers formed by different SHP-2 mutants displaying diverse biochemical properties were found to respond differently to epidermal growth factor (EGF) stimulation. Although this differential behavior cannot be rationalized mechanistically yet, these findings suggest a possible regulatory role of dimerization in SHP-2 function.

Features of the reversible sensitivity-resistance transition in PI3K/PTEN/AKT signalling network after HER2 inhibition

Systems biology approaches that combine experimental data and theoretical modelling to understand cellular signalling network dynamics offer a useful platform to investigate the mechanisms of resistance to drug interventions and to identify combination drug treatments. Extending our work on modelling the PI3K/PTEN/AKT signalling network (SN), we analyse the sensitivity of the SN output signal, phospho-AKT, to inhibition of HER2 receptor. We model typical aberrations in this SN identified in cancer development and drug resistance: loss of PTEN activity, PI3K and AKT mutations, HER2 overexpression, and overproduction of GSK3β and CK2 kinases controlling PTEN phosphorylation. We show that HER2 inhibition by the monoclonal antibody pertuzumab increases SN sensitivity, both to external signals and to changes in kinetic parameters of the proteins and their expression levels induced by mutations in the SN. This increase in sensitivity arises from the transition of SN functioning from saturation to non-saturation mode in response to HER2 inhibition. PTEN loss or PIK3CA mutation causes resistance to anti-HER2 inhibitor and leads to the restoration of saturation mode in SN functioning with a consequent decrease in SN sensitivity. We suggest that a drug-induced increase in SN sensitivity to internal perturbations, and specifically mutations, causes SN fragility. In particular, the SN is vulnerable to mutations that compensate for drug action and this may result in a sensitivity-to-resistance transition. The combination of HER2 and PI3K inhibition does not sensitise the SN to internal perturbations (mutations) in the PI3K/PTEN/AKT pathway: this combination treatment provides both synergetic inhibition and may prevent the SN from acquired mutations causing drug resistance. Through combination inhibition treatments, we studied the impact of upstream and downstream interventions to suppress resistance to the HER2 inhibitor in the SN with PTEN loss. Comparison of experimental results of PI3K inhibition in the PTEN upstream pathway with PDK1 inhibition in the PTEN downstream pathway shows that upstream inhibition abrogates resistance to pertuzumab more effectively than downstream inhibition. This difference in inhibition effect arises from the compensatory mechanism of an activation loop induced in the downstream pathway by PTEN loss. We highlight that drug target identification for combination anti-cancer therapy needs to account for the mutation effects on the upstream and downstream pathways.

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.

Dynamic modeling and analysis of cancer cellular network motifs

With the advent of high-throughput biology, we now routinely scan cells and organisms at practically all levels, from genome to protein, metabolism, signaling and other cellular functions. This methodology allowed biological studies to move from a reductionist approach, such as isolation of specific pathways and mechanisms, to a more integrative approach, where biological systems are seen as a network of interconnected components that provide specific outputs and functions in response to stimuli. Recent literature on biological networks demonstrates two important concepts that we will consider in this review: (i) cellular pathways are highly interconnected and should not be studied separately, but as a network; (ii) simple, recurrent feedback motifs within the network can produce very specific functions that favor their modular use. The first theme differs from the traditional approach in biology because it provides a framework (i.e., the network view) in which large datasets are analyzed with an unbiased view. The second theme (feedback motifs) shows the importance of locally analyzing the dynamic properties of biological networks in order to better understand their functionality. We will review these themes with examples from cell signaling networks, gene regulatory networks and metabolic pathways. The deregulation of cellular networks (metabolism, signaling etc.) is involved in cancer, but the size of the networks and resulting non-linear behavior do not allow for intuitive reasoning. In that context, we argue that the qualitative classification of the 'building blocs' of biological networks (i.e. the motifs) in terms of dynamics and functionality will be critical to improve our understanding of cancer biology and rationalize the wealth of information from high-throughput experiments. From the examples highlighted in this review, it is clear that dynamic feedback motifs can be used to provide a unified view of various cellular processes involved in cancer and this will be critical for future research on personalized and predictive cancer therapies.

Computational model of EGFR and IGF1R pathways in lung cancer: a Systems Biology approach for Translational Oncology

In this paper we propose a Systems Biology approach to understand the molecular biology of the Epidermal Growth Factor Receptor (EGFR, also known as ErbB1/HER1) and type 1 Insulin-like Growth Factor (IGF1R) pathways in non-small cell lung cancer (NSCLC). This approach, combined with Translational Oncology methodologies, is used to address the experimental evidence of a close relationship among EGFR and IGF1R protein expression, by immunohistochemistry (IHC) and gene amplification, by in situ hybridization (FISH) and the corresponding ability to develop a more aggressive behavior. We develop a detailed in silico model, based on ordinary differential equations, of the pathways and study the dynamic implications of receptor alterations on the time behavior of the MAPK cascade down to ERK, which in turn governs proliferation and cell migration. In addition, an extensive sensitivity analysis of the proposed model is carried out and a simplified model is proposed which allows us to infer a similar relationship among EGFR and IGF1R activities and disease outcome.

Computation of restoration of ligand response in the random kinetics of a prostate cancer cell signaling pathway

Mutation and/or dysfunction of signaling proteins in the mitogen activated protein kinase (MAPK) signal transduction pathway are frequently observed in various kinds of human cancer. Consistent with this fact, in the present study, we experimentally observe that the epidermal growth factor (EGF) induced activation profile of MAP kinase signaling is not straightforward dose-dependent in the PC3 prostate cancer cells. To find out what parameters and reactions in the pathway are involved in this departure from the normal dose-dependency, a model-based pathway analysis is performed. The pathway is mathematically modeled with 28 rate equations yielding those many ordinary differential equations (ODE) with kinetic rate constants that have been reported to take random values in the existing literature. This has led to us treating the ODE model of the pathways kinetics as a random differential equations (RDE) system in which the parameters are random variables. We show that our RDE model captures the uncertainty in the kinetic rate constants as seen in the behavior of the experimental data and more importantly, upon simulation, exhibits the abnormal EGF dose-dependency of the activation profile of MAP kinase signaling in PC3 prostate cancer cells. The most likely set of values of the kinetic rate constants obtained from fitting the RDE model into the experimental data is then used in a direct transcription based dynamic optimization method for computing the changes needed in these kinetic rate constant values for the restoration of the normal EGF dose response. The last computation identifies the parameters, i.e., the kinetic rate constants in the RDE model, that are the most sensitive to the change in the EGF dose response behavior in the PC3 prostate cancer cells. The reactions in which these most sensitive parameters participate emerge as candidate drug targets on the signaling pathway.