A Systems Biology Approach to Understand the Racial Disparities in Colorectal Cancer

Racial disparities between Black/African Americans (AA) and White patients in colorectal cancer are an ever-growing area of concern. Black/AA show the highest incidence and have the highest mortality among major U.S. racial groups. There is no definite cause other than possible sociodemographic, socioeconomic, education, nutrition, delivery of healthcare, screening, and cultural factors. A primary limitation in this field is the lack of and small sample size of Black/AA studies. Thus, this study aimed to investigate whether differences in gene expression contribute to this ongoing unanswered racial disparity issue. In this study, we examined transcriptomic data of Black/AA and White patient cohorts using a bioinformatic and systems biology approach. We performed a Kaplan-Meier overall survival analysis between both patient cohorts across critical colorectal cancer signal transduction networks (STN), to determine the differences in significant genes across each cohort. Other bioinformatic analyses performed included PROGENy (pathway responsive genes for activity inference), RNA sequencing differential expression using DESeq2, multivariable-adjusted regression, and other associated Kaplan-Meier analyses. These analyses identified novel prognostic genes independent from each cohort, 176 differentially expressed genes, and specific patient cohort STN survival associations. Despite the overarching limitation, the results revealed several novel differences in gene expression between the colorectal cancer Black/AA and White patient cohorts, which allows one to dive deeper into and understand the behavior on a systems level of what could be driving this racial difference across colorectal cancer. Concretely, this information can guide precision medicine approaches tailored specifically for colorectal cancer racial disparities.

SIGNIFICANCE

The purpose of this work is to investigate the racial disparities in colorectal cancer between Black/AA and White patient cohorts using a systems biology and bioinformatic approach. Our study investigates the underlying biology of each patient cohort. Concretely, the findings of this study include disparity-associated genes and pathways, which provide a tangible starting point to guide precision medicine approaches tailored specifically for colorectal cancer racial disparities.

A Negative Feedback Loop and Transcription Factor Cooperation Regulate Zonal Gene Induction by 2, 3, 7, 8-Tetrachlorodibenzo-p-Dioxin in the Mouse Liver

The cytochrome P450 (Cyp) proteins Cyp1A1 and Cyp1A2 are strongly induced in the mouse liver by the potent environmental toxicant 2, 3, 7, 8‐tetrachlorodibenzo‐p‐dioxin (TCDD), acting through the aryl hydrocarbon receptor (AHR). The induction of Cyp1A1 is localized within the centrilobular regions of the mouse liver at low doses of TCDD, progressing to pan‐lobular induction at higher doses. Even without chemical perturbation, metabolic functions and associated genes are basally zonated in the liver lobule along the central‐to‐portal axis. To investigate the mechanistic basis of spatially restricted gene induction by TCDD, we have developed a multiscale computational model of the mouse liver lobule with single‐cell resolution. The spatial location of individual hepatocytes in the model was calibrated from previously published high‐resolution images. A systems biology model of the network of biochemical signaling pathways underlying Cyp1A1 and Cyp1A2 induction was then incorporated into each hepatocyte in the model. Model simulations showed that a negative feedback loop formed by binding of the induced Cyp1A2 protein to TCDD, together with cooperative gene induction by the β‐catenin/AHR/TCDD transcription factor complex and β‐catenin, help produce the spatially localized induction pattern of Cyp1A1. Although endogenous WNT regulates the metabolic zonation of many genes, it was not a driver of zonal Cyp1A1 induction in our model. Conclusion: In this work, we used data‐driven computational modeling to identify the mechanistic basis of zonally restricted gene expression induced by the potent and persistent environmental pollutant TCDD. The multiscale model and derived results clarify the mechanisms of dose‐dependent hepatic gene induction responses to TCDD. Additionally, this work contributes to our broader understanding of spatial gene regulation along the liver lobule.

Personalised Medicine for Colorectal Cancer Using Mechanism-Based Machine Learning Models

Gaining insight into the mechanisms of signal transduction networks (STNs) by using critical features from patient-specific mathematical models can improve patient stratification and help to identify potential drug targets. To achieve this, these models should focus on the critical STNs for each cancer, include prognostic genes and proteins, and correctly predict patient-specific differences in STN activity. Focussing on colorectal cancer and the WNT STN, we used mechanism-based machine learning models to identify genes and proteins with significant associations to event-free patient survival and predictive power for explaining patient-specific differences of STN activity. First, we identified the WNT pathway as the most significant pathway associated with event-free survival. Second, we built linear-regression models that incorporated both genes and proteins from established mechanistic models in the literature and novel genes with significant associations to event-free patient survival. Data from The Cancer Genome Atlas and Clinical Proteomic Tumour Analysis Consortium were used, and patient-specific STN activity scores were computed using PROGENy. Three linear regression models were built, based on; (1) the gene-set of a state-of-the-art mechanistic model in the literature, (2) novel genes identified, and (3) novel proteins identified. The novel genes and proteins were genes and proteins of the extant WNT pathway whose expression was significantly associated with event-free survival. The results show that the predictive power of a model that incorporated novel event-free associated genes is better compared to a model focussing on the genes of a current state-of-the-art mechanistic model. Several significant genes that should be integrated into future mechanistic models of the WNT pathway are DVL3, FZD5, RAC1, ROCK2, GSK3B, CTB2, CBT1, and PRKCA. Thus, the study demonstrates that using mechanistic information in combination with machine learning can identify novel features (genes and proteins) that are important for explaining the STN heterogeneity between patients and their association to clinical outcomes.

Identification of the Role of Wnt/β-Catenin Pathway Through Integrated Analyses and in vivo Experiments in Vitiligo

Purpose

Vitiligo is an acquired depigmentation skin disease, which affects an average of 1% of the world’s population. The purpose of this study is to identify the key genes and pathways responsible for vitiligo and find new therapeutic targets.

Methods

The datasets GSE65127, GSE53146, and GSE75819 were downloaded from the Gene Expression Omnibus (GEO) database. R language was used to identify the differentially expressed genes (DEGs) between lesional skin of vitiligo and non-lesional skin. Next, the key pathways were obtained by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. The protein–protein interaction (PPI) networks were conducted by STRING database and Cytoscape software. Subsequently, module analysis was performed by Cytoscape. Among these results, the Wnt/β-catenin pathway and melanogenesis pathway caught our attention. The expression level of β-catenin, microphthalmia-associated transcription factor (MITF) and tyrosinase (TYR) was detected by immunofluorescence in vitiligo lesions and healthy skin. Moreover, zebrafish was treated with XAV-939, an inhibitor of the Wnt/β-catenin pathway. After that, the area of melanin granules as a percentage of the head area was measured. The mRNA expression of β-catenin, lymphoid-enhancing factor 1(lef1), tyr and mitf were detected by q-PCR (quantitative polymerase chain reaction) in zebrafish (Danio rerio).

Results

A total of 2442 DEGs were identified, including 1068 upregulated and 1374 downregulated DEGs. The key pathways were identified by GO and KEGG analyses, such as “NOD-like receptor signaling pathway”, “Wnt signaling pathway”, “Melanogenesis”, “mTOR signaling pathway”, “PI3K-Akt signaling pathway”, “Calcium signaling pathway” and “Rap1 signaling pathway”. The immunofluorescence results showed that the level of β-catenin, MITF and TYR was significantly downregulated in vitiligo lesional skin. In zebrafish, the mean percentage area of melanin granules and the expression of β-catenin, lef1, tyr and mitf were decreased after treated with XAV-939.

Conclusion

The present study identified key genes and signaling pathways associated with the pathophysiology of vitiligo. Among them, the Wnt/β-catenin pathway played an essential role in pigmentation and could be a breakthrough point in vitiligo treatment.

Relating simulation studies by provenance-Developing a family of Wnt signaling models

For many biological systems, a variety of simulation models exist. A new simulation model is rarely developed from scratch, but rather revises and extends an existing one. A key challenge, however, is to decide which model might be an appropriate starting point for a particular problem and why. To answer this question, we need to identify entities and activities that contributed to the development of a simulation model. Therefore, we exploit the provenance data model, PROV-DM, of the World Wide Web Consortium and, building on previous work, continue developing a PROV ontology for simulation studies. Based on a case study of 19 Wnt/β-catenin signaling models, we identify crucial entities and activities as well as useful metadata to both capture the provenance information from individual simulation studies and relate these forming a family of models. The approach is implemented in WebProv, a web application for inserting and querying provenance information. Our specialization of PROV-DM contains the entities Research Question, Assumption, Requirement, Qualitative Model, Simulation Model, Simulation Experiment, Simulation Data, and Wet-lab Data as well as activities referring to building, calibrating, validating, and analyzing a simulation model. We show that most Wnt simulation models are connected to other Wnt models by using (parts of) these models. However, the overlap, especially regarding the Wet-lab Data used for calibration or validation of the models is small. Making these aspects of developing a model explicit and queryable is an important step for assessing and reusing simulation models more effectively. Exposing this information helps to integrate a new simulation model within a family of existing ones and may lead to the development of more robust and valid simulation models. We hope that our approach becomes part of a standardization effort and that modelers adopt the benefits of provenance when considering or creating simulation models.

We revise a provenance ontology for simulation studies of cellular biochemical models. Provenance information is useful for understanding the creation of a simulation model because it not only contains information about the entities and activities that have led to a simulation model but also their relations, all of which can be visualized. It provides additional structure by explicitly recording research questions, assumptions, and requirements and relating them along with data, qualitative models, simulation models, and simulation experiments through a small set of predefined but extensible activities.

We have applied our concept to a family of 19 Wnt signaling models and implemented a web-based tool (WebProv) to store the provenance information from these studies. The resulting provenance graph visualizes the story line of simulation studies and demonstrates the creation and calibration of simulation models, the successive attempts of validation and extension, and shows, beyond an individual simulation study, how the Wnt models are related. Thereby, the steps and sources that contributed to a simulation model are made explicit.

Our approach complements other approaches aimed at facilitating the reuse and assessment of simulation products in systems biology such as model repositories as well as annotation and documentation guidelines.

Paracrine and autocrine regulation of gene expression by Wnt-inhibitor Dickkopf in wild-type and mutant hepatocytes

Background

Cells are able to communicate and coordinate their function within tissues via secreted factors. Aberrant secretion by cancer cells can modulate this intercellular communication, in particular in highly organised tissues such as the liver. Hepatocytes, the major cell type of the liver, secrete Dickkopf (Dkk), which inhibits Wnt/ β-catenin signalling in an autocrine and paracrine manner. Consequently, Dkk modulates the expression of Wnt/ β-catenin target genes. We present a mathematical model that describes the autocrine and paracrine regulation of hepatic gene expression by Dkk under wild-type conditions as well as in the presence of mutant cells.

Results

Our spatial model describes the competition of Dkk and Wnt at receptor level, intra-cellular Wnt/ β-catenin signalling, and the regulation of target gene expression for 21 individual hepatocytes. Autocrine and paracrine regulation is mediated through a feedback mechanism via Dkk and Dkk diffusion along the porto-central axis. Along this axis an APC concentration gradient is modelled as experimentally detected in liver. Simulations of mutant cells demonstrate that already a single mutant cell increases overall Dkk concentration. The influence of the mutant cell on gene expression of surrounding wild-type hepatocytes is limited in magnitude and restricted to hepatocytes in close proximity. To explore the underlying molecular mechanisms, we perform a comprehensive analysis of the model parameters such as diffusion coefficient, mutation strength and feedback strength.

Conclusions

Our simulations show that Dkk concentration is elevated in the presence of a mutant cell. However, the impact of these elevated Dkk levels on wild-type hepatocytes is confined in space and magnitude. The combination of inter- and intracellular processes, such as Dkk feedback, diffusion and Wnt/ β-catenin signal transduction, allow wild-type hepatocytes to largely maintain their gene expression.

Electronic supplementary material

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

Extracellular inhibitors can attenuate tumorigenic Wnt pathway activity in adenomatous polyposis coli mutants: Predictions of a validated mathematical model

Background

Despite considerable investigational efforts, no method to overcome the pathogenesis caused by loss of function (LoF) mutations in tumor suppressor genes has been successfully translated to the clinic. The most frequent LoF mutation in human cancers is Adenomatous polyposis coli (APC), causing aberrant activation of the Wnt pathway. In nearly all colon cancer tumors, the APC protein is truncated, but still retains partial binding abilities.

Objective & methods

Here, we tested the hypothesis that extracellular inhibitors of the Wnt pathway, although acting upstream of the APC mutation, can restore normal levels of pathway activity in colon cancer cells. To this end, we developed and simulated a mathematical model for the Wnt pathway in different APC mutants, with or without the effects of the extracellular inhibitors, Secreted Frizzled-Related Protein1 (sFRP1) and Dickhopf1 (Dkk1). We compared our model predictions to experimental data in the literature.

Results

Our model accurately predicts T-cell factor (TCF) activity in mutant cells that vary in APC mutation. Model simulations suggest that both sFRP1 and DKK1 can reduce TCF activity in APC^1638N/1572T and Apc^min/min mutants, but restoration of normal activity levels is possible only in the former. When applied in combination, synergism between the two inhibitors can reduce their effective doses to one-fourth of the doses required under single inhibitor application. Overall, re-establishment of normal Wnt pathway activity is predicted for every APC mutant in whom TCF activity is increased by up to 11 fold.

Conclusions

Our work suggests that extracellular inhibitors can effectively restore normal Wnt pathway activity in APC-truncated cancer cells, even though these LoF mutations occur downstream of the inhibitory action. The insufficient activity of the truncated APC can be quantitatively balanced by the upstream intervention. This new concept of upstream intervention to control the effects of downstream mutations may be considered also for other partial LoF mutations in other signaling pathways.

Βeta-catenin N-terminal domain: An enigmatic region prone to cancer causing mutations

The Wnt/β-catenin is a highly conserved signaling pathway involved in cell fate decisions during various stages of development. Dysregulation of canonical Wnt/β-catenin signaling has been associated with various diseases including cancer. β-Catenin, the central component of canonical Wnt signaling pathway, is a multi-functional protein playing both structural and signaling roles. β-Catenin is composed of three distinct domains: N-terminal domain, C-terminal domain and a central armadillo repeat domain. N-terminal domain of β-catenin harbours almost all of the cancer causing mutations, thus deciphering its critical structural and functional roles offers great potential in cancer detection and therapy. Here, in this review, we have collected information from pharmacological analysis, bio-physical and structural studies, molecular modeling, in-vivo and in-vitro assays, and transgenic animal experiments employing various N-terminal domain variants of β-catenin to discuss the interaction of β-catenin with its binding partners that specifically interact with this domain and the implications of these interactions on signaling, cell fate determination, and in tumorigenesis. A thorough understanding of interactions between β-catenin and its binding partners will enable us to more effectively understand how β-catenin switches between its multiple roles, and will lead to the development of specific assays for the identification of small molecules as chemotherapeutic agents to treat diseases, including cancer and neurological disorders, where Wnt/β-catenin signaling is dysregulated.

An ongoing role for Wnt signaling in differentiating melanocytes in vivo

A role for Wnt signaling in melanocyte specification from neural crest is conserved across vertebrates, but possible ongoing roles in melanocyte differentiation have received little attention. Using a systems biology approach to investigate the gene regulatory network underlying stable melanocyte differentiation in zebrafish highlighted a requirement for a positive-feedback loop involving the melanocyte master regulator Mitfa. Here, we test the hypothesis that Wnt signaling contributes to that positive feedback. We show firstly that Wnt signaling remains active in differentiating melanocytes and secondly that enhanced Wnt signaling drives elevated transcription of mitfa. We show that chemical activation of the Wnt signaling pathway at early stages of melanocyte development enhances melanocyte specification as expected, but importantly that at later (differentiation) stages, it results in altered melanocyte morphology, although melanisation is not obviously affected. Downregulation of Wnt signaling also results in altered melanocyte morphology and organization. We conclude that Wnt signaling plays a role in regulating ongoing aspects of melanocyte differentiation in zebrafish.

The Integrated Role of Wnt/β-Catenin, N-Glycosylation, and E-Cadherin-Mediated Adhesion in Network Dynamics

The cellular network composed of the evolutionarily conserved metabolic pathways of protein N-glycosylation, Wnt/β-catenin signaling pathway, and E-cadherin-mediated cell-cell adhesion plays pivotal roles in determining the balance between cell proliferation and intercellular adhesion during development and in maintaining homeostasis in differentiated tissues. These pathways share a highly conserved regulatory molecule, β-catenin, which functions as both a structural component of E-cadherin junctions and as a co-transcriptional activator of the Wnt/β-catenin signaling pathway, whose target is the N-glycosylation-regulating gene, DPAGT1. Whereas these pathways have been studied independently, little is known about the dynamics of their interaction. Here we present the first numerical model of this network in MDCK cells. Since the network comprises a large number of molecules with varying cell context and time-dependent levels of expression, it can give rise to a wide range of plausible cellular states that are difficult to track. Using known kinetic parameters for individual reactions in the component pathways, we have developed a theoretical framework and gained new insights into cellular regulation of the network. Specifically, we developed a mathematical model to quantify the fold-change in concentration of any molecule included in the mathematical representation of the network in response to a simulated activation of the Wnt/ β-catenin pathway with Wnt3a under different conditions. We quantified the importance of protein N-glycosylation and synthesis of the DPAGT1 encoded enzyme, GPT, in determining the abundance of cytoplasmic β-catenin. We confirmed the role of axin in β-catenin degradation. Finally, our data suggest that cell-cell adhesion is insensitive to E-cadherin recycling in the cell. We validate the model by inhibiting β-catenin-mediated activation of DPAGT1 expression and predicting changes in cytoplasmic β-catenin concentration and stability of E-cadherin junctions in response to DPAGT1 inhibition. We show the impact of pathway dysregulation through measurements of cell migration in scratch-wound assays. Collectively, our results highlight the importance of numerical analyses of cellular networks dynamics to gain insights into physiological processes and potential design of therapeutic strategies to prevent epithelial cell invasion in cancer.

Construction and Experimental Validation of a Petri Net Model of Wnt/β-Catenin Signaling

The Wnt/β-catenin signaling pathway is important for multiple developmental processes and tissue maintenance in adults. Consequently, deregulated signaling is involved in a range of human diseases including cancer and developmental defects. A better understanding of the intricate regulatory mechanism and effect of physiological (active) and pathophysiological (hyperactive) WNT signaling is important for predicting treatment response and developing novel therapies. The constitutively expressed CTNNB1 (commonly and hereafter referred to as β-catenin) is degraded by a destruction complex, composed of amongst others AXIN1 and GSK3. The destruction complex is inhibited during active WNT signaling, leading to β-catenin stabilization and induction of β-catenin/TCF target genes. In this study we investigated the mechanism and effect of β-catenin stabilization during active and hyperactive WNT signaling in a combined in silico and in vitro approach. We constructed a Petri net model of Wnt/β-catenin signaling including main players from the plasma membrane (WNT ligands and receptors), cytoplasmic effectors and the downstream negative feedback target gene AXIN2. We validated that our model can be used to simulate both active (WNT stimulation) and hyperactive (GSK3 inhibition) signaling by comparing our simulation and experimental data. We used this experimentally validated model to get further insights into the effect of the negative feedback regulator AXIN2 upon WNT stimulation and observed an attenuated β-catenin stabilization. We furthermore simulated the effect of APC inactivating mutations, yielding a stabilization of β-catenin levels comparable to the Wnt-pathway activities observed in colorectal and breast cancer. Our model can be used for further investigation and viable predictions of the role of Wnt/β-catenin signaling in oncogenesis and development.

Dissecting the role of Wnt signaling and its interactions with FGF signaling during midbrain neurogenesis

Interactions between FGF and Wnt/ bcat signaling control development of the midbrain. The nature of this interaction and how these regulate patterning, growth and differentiation is less clear, as it has not been possible to temporally dissect the effects of one pathway relative to the other. We have employed pharmacological and genetic tools to probe the temporal and spatial roles of FGF and Wnt in controlling the specification of early midbrain neurons. We identify a β-catenin (bcat) independent role for GSK-3 in modulating FGF activity and hence neuronal patterning. This function is complicated by an overlap with bcat-dependent regulation of FGF signaling, through the regulation of sprouty4. Additionally we reveal how attenuation of Axin protein function can promote fluctuating levels of bcat activity that are dependent on FGF activity. This highlights the complex nature of the interactions between FGF and Wnt/ bcat and reveals that they act at multiple levels to control each others activity in the midbrain.

A regulatory receptor network directs the range and output of the Wingless signal

The potent activity of Wnt/Wingless (Wg) signals necessitates sophisticated mechanisms that spatially and temporally regulate their distribution and range of action. The two main receptor components for Wg - Arrow (Arr) and Frizzled 2 (Fz2) - are transcriptionally downregulated by Wg signaling, thus forming gradients that oppose that of Wg. Here, we analyze the relevance of this transcriptional regulation for the formation of the Wg gradient in the Drosophila wing disc by combining in vivo receptor overexpression with an in silico model of Wg receptor interactions. Our experiments show that ubiquitous upregulation of Arr and Fz2 has no significant effects on Wg output, whereas clonal overexpression of these receptors leads to signaling discontinuities that have detrimental phenotypic consequences. These findings are supported by our in silico model for Wg diffusion and signal transduction, which suggests that abrupt changes in receptor levels causes discontinuities in Wg signaling. Furthermore, we identify a 200 bp regulatory element in the arr locus that can account for the Arr gradient, and we show that this is indirectly negatively controlled by Wg activity. Finally, we analyze the role of Frizzled 3 (Fz3) in this system and find that its expression, which is induced by Wg, contributes to the establishment of the Arr and Fz2 gradients through counteracting canonical signaling. Taken together, our results provide a model in which the regulatory network of Wg and the three receptor components account for the range and shape of this prototypical morphogen system.

Analysis of Wnt signaling β-catenin spatial dynamics in HEK293T cells

Background

Wnt/β-catenin signaling is involved in different stages of mammalian development and implicated in various cancers (e.g. colorectal cancer). Recent experimental and computational studies have revealed characteristics of the pathway, however a cell-specific spatial perspective is lacking. In this study, a novel 3D confocal quantitation protocol is developed to acquire spatial (two cellular compartments: nucleus and cytosol-membrane) and temporal quantitative data on target protein (e.g. β-catenin) concentrations in Human Epithelial Kidney cells (HEK293T) during perturbation (with either cycloheximide or Wnt3A). Computational models of the Wnt pathway are constructed and interrogated based on this data.

Results

A single compartment Wnt pathway model is compared with a simple β-catenin two compartment model to investigate Wnt3A signaling in HEK293T cells. When protein synthesis is inhibited, β-catenin decreases at the same rate in both cellular compartments, suggesting diffusional transport is fast compared to β-catenin degradation in the cytosol. With Wnt3A stimulation, the total amount of β-catenin rises throughout the cell, however the increase is initially (~first hour) faster in the nuclear compartment. While both models were able to reproduce the whole cell changes in β-catenin, only the compartment model reproduced the Wnt3A induced changes in β-catenin distribution and it was also the best fit for the data obtained when active transport was included alongside passive diffusion transport.

Conclusions

This integrated 3D quantitation imaging protocol and computational modeling approach allowed cell-specific compartment models of the signaling pathways to be constructed and analyzed. The Wnt models constructed in this study are the first for HEK293T and have suggested potential roles of inter-compartment transport to the dynamics of signaling.

Analysing the impact of nucleo-cytoplasmic shuttling of β-catenin and its antagonists APC, Axin and GSK3 on Wnt/β-catenin signalling

The canonical Wnt signalling pathway plays a critical role in development and disease. The key player of the pathway is β-catenin. Its activity is mainly regulated by the destruction complex consisting of APC, Axin and GSK3. In the nucleus, the complex formation of β-catenin and TCF initiates target gene expression. Our study provides a comprehensive analysis of the role of nucleo-cytoplasmic shuttling of APC, Axin, and GSK3 and the inactivation of β-catenin by the destruction complex in Wnt/β-catenin signalling. We address the following questions: Can nucleo-cytoplasmic shuttling of APC, Axin and GSK3 increase the [β-catenin/TCF] concentration? And, how is the [β-catenin/TCF] concentration influenced by phosphorylation and subsequent degradation of nuclear β-catenin? Based on experimental findings, we develop a compartmental model and conduct several simulation experiments. Our analysis reveals the following key findings: 1) nucleo-cytoplasmic shuttling of β-catenin and its antagonists can yield a spatial separation between the said proteins, which results in a breakdown of β-catenin degradation, followed by an accumulation of β-catenin and hence leads to an increase of the [β-catenin/TCF] concentration. Our results strongly suggest that Wnt signalling can benefit from nucleo-cytoplasmic shuttling of APC, Axin and GSK3, although they are in general β-catenin antagonising proteins. 2) The total robustness of the [β-catenin/TCF] output is closely linked to its absolute concentration levels. We demonstrate that the compartmental separation of β-catenin and the destruction complex does not only lead to a maximization, but additionally to an increased robustness of [β-catenin/TCF] signalling against perturbations in the cellular environment. 3) A nuclear accumulation of the destruction complex renders the pathway robust against fluctuations in Wnt signalling and against changes in the compartmental distribution of β-catenin. 4) Elucidating the impact of destruction complex inhibition, we show that the [β-catenin/TCF] concentration is more effectively enhanced by inhibition of the kinase GSK3 rather than the binding of β-catenin to the destruction complex.

Modeling Wnt/β-Catenin Target Gene Expression in APC and Wnt Gradients Under Wild Type and Mutant Conditions

The Wnt/β-catenin pathway is involved in the regulation of a multitude of physiological processes by controlling the differential expression of target genes. In certain tissues such as the adult liver, the Wnt/β-catenin pathway can attain different levels of activity due to gradients of Wnt ligands and/or intracellular pathway components like APC. How graded pathway activity is converted into regionally distinct patterns of Wnt/β-catenin target gene expression is largely unknown. Here, we apply a mathematical modeling approach to investigate the impact of different regulatory mechanisms on target gene expression within Wnt or APC concentration gradients. We develop a minimal model of Wnt/β-catenin signal transduction and combine it with various mechanisms of target gene regulation. In particular, the effects of activation, inhibition, and an incoherent feedforward loop (iFFL) are compared. To specify activation kinetics, we analyze experimental data that quantify the response of β-catenin/TCF reporter constructs to different Wnt concentrations, and demonstrate that the induction of these constructs occurs in a cooperative manner with Hill coefficients between 2 and 5. In summary, our study shows that the combination of specific gene regulatory mechanisms with a time-independent gradient of Wnt or APC is sufficient to generate distinct target gene expression patterns as have been experimentally observed in liver. We find that cooperative gene activation in combination with a TCF feedback can establish sharp borders of target gene expression in Wnt or APC gradients. In contrast, the iFFL renders gene expression independent of gradients of the upstream signaling components. Our subsequent analysis of carcinogenic pathway mutations reveals that their impact on gene expression is determined by the gene regulatory mechanism and the APC concentration of the cell in which the mutation occurs.

Elucidating the sources of β-catenin dynamics in human neural progenitor cells

Human neural progenitor cells (hNPCs) form a new prospect for replacement therapies in the context of neurodegenerative diseases. The Wnt/β-catenin signaling pathway is known to be involved in the differentiation process of hNPCs. RVM cells form a common cell model of hNPCs for in vitro investigation. Previous observations in RVM cells raise the question of whether observed kinetics of the Wnt/β-catenin pathway in later differentiation phases are subject to self-induced signaling. However, a concern when investigating RVM cells is that experimental results are possibly biased by the asynchrony of cells w.r.t. the cell cycle. In this paper, we present, based on experimental data, a computational modeling study on the Wnt/β-catenin signaling pathway in RVM cell populations asynchronously distributed w.r.t. to their cell cycle phases. Therefore, we derive a stochastic model of the pathway in single cells from the reference model in literature and extend it by means of cell populations and cell cycle asynchrony. Based on this, we show that the impact of the cell cycle asynchrony on wet-lab results that average over cell populations is negligible. We then further extend our model and the thus-obtained simulation results provide additional evidence that self-induced Wnt signaling occurs in RVM cells. We further report on significant stochastic effects that directly result from model parameters provided in literature and contradict experimental observations.

A new validated mathematical model of the Wnt signalling pathway predicts effective combinational therapy by sFRP and Dkk

The Wnt signalling pathway controls cell proliferation and differentiation, and its deregulation is implicated in different diseases including cancer. Learning how to manipulate this pathway could substantially contribute to the development of therapies. We developed a mathematical model describing the initial sequence of events in the Wnt pathway, from ligand binding to β-catenin accumulation, and the effects of inhibitors, such as sFRPs (secreted Frizzled-related proteins) and Dkk (Dickkopf). Model parameters were retrieved from experimental data reported previously. The model was retrospectively validated by accurately predicting the effects of Wnt3a and sFRP1 on β-catenin levels in two independent published experiments (R2 between 0.63 and 0.91). Prospective validation was obtained by testing the model's accuracy in predicting the effect of Dkk1 on Wnt-induced β-catenin accumulation (R2≈0.94). Model simulations under different combinations of sFRP1 and Dkk1 predicted a clear synergistic effect of these two inhibitors on β-catenin accumulation, which may point towards a new treatment avenue. Our model allows precise calculation of the effect of inhibitors applied alone or in combination, and provides a flexible framework for identifying potential targets for intervention in the Wnt signalling pathway.

Expression and clinical significance of von Hippel-Lindau downstream genes: Jade-1 and β-catenin related to renal cell carcinoma

OBJECTIVE

To investigate the expression of Jade-1 and β-catenin and their effect in the development of renal cell carcinoma (RCC).

METHODS

The expression of Jade-1 and β-catenin in 11 normal kidney tissue specimens (normal group) and 60 RCC specimens (RCC group) was determined using real-time polymerase chain reaction and immunohistochemistry. Also, their effect on early relapses of RCC was analyzed after 5 years of follow-up.

RESULTS

The expression of Jade-1 protein in the RCC group was significantly lower than that in the normal group (0.1655 vs 0.7438, P < .05), and the expression of β-catenin protein was significantly greater than that in the normal group (0.2756 vs 0.0855, P < .05). The expression of Jade-1 mRNA in the RCC group was 0.202 times that in the normal group, which were lower (P < .05). The expression of β-catenin mRNA was 1.014 times that in the normal group (P > .05). The expression of Jade-1 protein and β-catenin protein in poorly differentiated RCC specimens was significantly lower and higher than the expression in the well-differentiated RCC specimen (P < .05), respectively. Patients with negative Jade-1 expression and positive β-catenin expression were found to have shorter survival on univariate analysis (P < .05).

CONCLUSION

RCC with a low expression of Jade-1 and high expression of β-catenin is associated with a poor outcome and decreased survival.

Modeling mechanistic biological networks: an advanced Boolean approach

MOTIVATION

The understanding of the molecular sources for diseases like cancer can be significantly improved by computational models. Recently, Boolean networks have become very popular for modeling signaling and regulatory networks. However, such models rely on a set of Boolean functions that are in general not known. Unfortunately, while detailed information on the molecular interactions becomes available in large scale through electronic databases, the information on the Boolean functions does not become available simultaneously and has to be included manually into the models, if at all known.

RESULTS

We propose a new Boolean approach which can directly utilize the mechanistic network information available through modern databases. The Boolean function is implicitly defined by the reaction mechanisms. Special care has been taken for the treatment of kinetic features like inhibition. The method has been applied to a signaling model combining the Wnt and MAPK pathway.

AVAILABILITY

A sample C++ implementation of the proposed method is available for Linux and compatible systems through http://code.google.com/p/libscopes/wiki/Paper2011.