Colorectal Cancer Through Simulation and Experiment

Colorectal cancer has continued to generate a huge amount of research interest over several decades, forming a canonical example of tumourigenesis since its use in Fearon and Vogelstein’s linear model of genetic mutation. Over time, the field has witnessed a transition from solely experimental work to the inclusion of mathematical biology and computer-based modelling. The fusion of these disciplines has the potential to provide valuable insights into oncologic processes, but also presents the challenge of uniting many diverse perspectives. Furthermore, the cancer cell phenotype defined by the ‘Hallmarks of Cancer’ has been extended in recent times and provides an excellent basis for future research. We present a timely summary of the literature relating to colorectal cancer, addressing the traditional experimental findings, summarising the key mathematical and computational approaches, and emphasising the role of the Hallmarks in current and future developments. We conclude with a discussion of interdisciplinary work, outlining areas of experimental interest which would benefit from the insight that mathematical and computational modelling can provide

Coupled Reversible and Irreversible Bistable Switches Underlying TGF-\beta-induced Epithelial to Mesenchymal Transition

Epithelial to mesenchymal transition (EMT) plays important roles in embryonic development, tissue regeneration and cancer metastasis. While several feedback loops have been shown to regulate EMT, it remains elusive how they coordinately modulate EMT response to TGF-\beta treatment. We construct a mathematical model for the core regulatory network controlling TGF-\beta-induced EMT. Through deterministic analyses and stochastic simulations, we show that EMT is a sequential two-step program that an epithelial cell first transits to partial EMT then to the mesenchymal state, depending on the strength and duration of TGF-\beta stimulation. Mechanistically the system is governed by coupled reversible and irreversible bistable switches. The SNAIL1/miR-34 double negative feedback loop is responsible for the reversible switch and regulates the initiation of EMT, while the ZEB/miR-200 feedback loop is accountable for the irreversible switch and controls the establishment of the mesenchymal state. Furthermore, an autocrine TGF-\beta/miR-200 feedback loop makes the second switch irreversible, modulating the maintenance of EMT. Such coupled bistable switches are robust to parameter variation and molecular noise. We provide a mechanistic explanation on multiple experimental observations. The model makes several explicit predictions on hysteretic dynamic behaviors, system response to pulsed stimulation and various perturbations, which can be straightforwardly tested.Comment: 32 pages, 8 figures, accepted by Biophysical Journa

Systems biology in inflammatory bowel diseases

Purpose of review: Ulcerative colitis (UC) and Crohn’s Disease (CD) are the two predominant types of inflammatory bowel disease (IBD), affecting over 1.4 million individuals in the US. IBD results from complex interactions between pathogenic components, including genetic and epigenetic factors, the immune response and the microbiome through an unknown sequence of events. The purpose of this review is to describe a system biology approach to IBD as a novel and exciting methodology aiming at developing novel IBD therapeutics based on the integration of molecular and cellular "omics" data. Recent Findings: Recent evidence suggested the presence of genetic, epigenetic, transcriptomic, proteomic and metabolomic alterations in IBD patients. Furthermore, several studies have shown that different cell types, including fibroblasts, epithelial, immune and endothelial cells together with the intestinal microbiota are involved in IBD pathogenesis. Novel computational methodologies have been developed aiming to integrate high - throughput molecular data. Summary: A systems biology approach could potentially identify the central regulators (hubs) in the IBD interactome and improve our understanding of the molecular mechanisms involved in IBD pathogenesis. The future IBD therapeutics should be developed on the basis of targeting the central hubs in the IBD network

A mathematical insight in the epithelial-mesenchymal-like transition in cancer cells and its effect in the invasion of the extracellular matrix

Current biological knowledge supports the existence of a secondary group of cancer cells within the body of the tumour that exhibits stem cell-like properties. These cells are termed Cancer Stem Cells (CSCs}, and as opposed to the more usual Differentiated Cancer Cells (DCCs), they exhibit higher motility, they are more resilient to therapy, and are able to metastasize to secondary locations within the organism and produce new tumours. The origin of the CSCs is not completely clear; they seem to stem from the DCCs via a transition process related to the Epithelial-Mesenchymal Transition (EMT) that can also be found in normal tissue. In the current work we model and numerically study the transition between these two types of cancer cells, and the resulting "ensemble" invasion of the extracellular matrix. This leads to the derivation and numerical simulation of two systems: an algebraic-elliptic system for the transition and an advection-reaction-diffusion system of Keller-Segel taxis type for the invasion

Quantifying cancer epithelial-mesenchymal plasticity and its association with stemness and immune response

Cancer cells can acquire a spectrum of stable hybrid epithelial/mesenchymal (E/M) states during epithelial-mesenchymal transition (EMT). Cells in these hybrid E/M phenotypes often combine epithelial and mesenchymal features and tend to migrate collectively commonly as small clusters. Such collectively migrating cancer cells play a pivotal role in seeding metastases and their presence in cancer patients indicates an adverse prognostic factor. Moreover, cancer cells in hybrid E/M phenotypes tend to be more associated with stemness which endows them with tumor-initiation ability and therapy resistance. Most recently, cells undergoing EMT have been shown to promote immune suppression for better survival. A systematic understanding of the emergence of hybrid E/M phenotypes and the connection of EMT with stemness and immune suppression would contribute to more effective therapeutic strategies. In this review, we first discuss recent efforts combining theoretical and experimental approaches to elucidate mechanisms underlying EMT multi-stability (i.e. the existence of multiple stable phenotypes during EMT) and the properties of hybrid E/M phenotypes. Following we discuss non-cell-autonomous regulation of EMT by cell cooperation and extracellular matrix. Afterwards, we discuss various metrics that can be used to quantify EMT spectrum. We further describe possible mechanisms underlying the formation of clusters of circulating tumor cells. Last but not least, we summarize recent systems biology analysis of the role of EMT in the acquisition of stemness and immune suppression.Comment: 50 pages, 6 figure

Integrated signaling pathway and gene expression regulatory model to dissect dynamics of <em>Escherichia coli </em>challenged mammary epithelial cells

AbstractCells transform external stimuli, through the activation of signaling pathways, which in turn activate gene regulatory networks, in gene expression. As more omics data are generated from experiments, eliciting the integrated relationship between the external stimuli, the signaling process in the cell and the subsequent gene expression is a major challenge in systems biology. The complex system of non-linear dynamic protein interactions in signaling pathways and gene networks regulates gene expression.The complexity and non-linear aspects have resulted in the study of the signaling pathway or the gene network regulation in isolation. However, this limits the analysis of the interaction between the two components and the identification of the source of the mechanism differentiating the gene expression profiles. Here, we present a study of a model of the combined signaling pathway and gene network to highlight the importance of integrated modeling.Based on the experimental findings we developed a compartmental model and conducted several simulation experiments. The model simulates the mRNA expression of three different cytokines (RANTES, IL8 and TNFα) regulated by the transcription factor NFκB in mammary epithelial cells challenged with E. coli. The analysis of the gene network regulation identifies a lack of robustness and therefore sensitivity for the transcription factor regulation. However, analysis of the integrated signaling and gene network regulation model reveals distinctly different underlying mechanisms in the signaling pathway responsible for the variation between the three cytokine's mRNA expression levels. Our key findings reveal the importance of integrating the signaling pathway and gene expression dynamics in modeling. Modeling infers valid research questions which need to be verified experimentally and can assist in the design of future biological experiments

Investigating biocomplexity through the agent-based paradigm.

Capturing the dynamism that pervades biological systems requires a computational approach that can accommodate both the continuous features of the system environment as well as the flexible and heterogeneous nature of component interactions. This presents a serious challenge for the more traditional mathematical approaches that assume component homogeneity to relate system observables using mathematical equations. While the homogeneity condition does not lead to loss of accuracy while simulating various continua, it fails to offer detailed solutions when applied to systems with dynamically interacting heterogeneous components. As the functionality and architecture of most biological systems is a product of multi-faceted individual interactions at the sub-system level, continuum models rarely offer much beyond qualitative similarity. Agent-based modelling is a class of algorithmic computational approaches that rely on interactions between Turing-complete finite-state machines--or agents--to simulate, from the bottom-up, macroscopic properties of a system. In recognizing the heterogeneity condition, they offer suitable ontologies to the system components being modelled, thereby succeeding where their continuum counterparts tend to struggle. Furthermore, being inherently hierarchical, they are quite amenable to coupling with other computational paradigms. The integration of any agent-based framework with continuum models is arguably the most elegant and precise way of representing biological systems. Although in its nascence, agent-based modelling has been utilized to model biological complexity across a broad range of biological scales (from cells to societies). In this article, we explore the reasons that make agent-based modelling the most precise approach to model biological systems that tend to be non-linear and complex