Linking Changes in Epithelial Morphogenesis to Cancer Mutations Using Computational Modeling

Most tumors arise from epithelial tissues, such as mammary glands and lobules, and their initiation is associated with the disruption of a finely defined epithelial architecture. Progression from intraductal to invasive tumors is related to genetic mutations that occur at a subcellular level but manifest themselves as functional and morphological changes at the cellular and tissue scales, respectively. Elevated proliferation and loss of epithelial polarization are the two most noticeable changes in cell phenotypes during this process. As a result, many three-dimensional cultures of tumorigenic clones show highly aberrant morphologies when compared to regular epithelial monolayers enclosing the hollow lumen (acini). In order to shed light on phenotypic changes associated with tumor cells, we applied the bio-mechanical IBCell model of normal epithelial morphogenesis quantitatively matched to data acquired from the non-tumorigenic human mammary cell line, MCF10A. We then used a high-throughput simulation study to reveal how modifications in model parameters influence changes in the simulated architecture. Three parameters have been considered in our study, which define cell sensitivity to proliferative, apoptotic and cell-ECM adhesive cues. By mapping experimental morphologies of four MCF10A-derived cell lines carrying different oncogenic mutations onto the model parameter space, we identified changes in cellular processes potentially underlying structural modifications of these mutants. As a case study, we focused on MCF10A cells expressing an oncogenic mutant HER2-YVMA to quantitatively assess changes in cell doubling time, cell apoptotic rate, and cell sensitivity to ECM accumulation when compared to the parental non-tumorigenic cell line. By mapping in vitro mutant morphologies onto in silico ones we have generated a means of linking the morphological and molecular scales via computational modeling. Thus, IBCell in combination with 3D acini cultures can form a computational/experimental platform for suggesting the relationship between the histopathology of neoplastic lesions and their underlying molecular defects

Models, measurement and inference in epithelial tissue dynamics

The majority of solid tumours arise in epithelia and therefore much research effort has gone into investigating the growth, renewal and regulation of these tissues. Here we review different mathematical and computational approaches that have been used to model epithelia. We compare different models and describe future challenges that need to be overcome in order to fully exploit new data which present, for the first time, the real possibility for detailed model validation and comparison

Advances in modelling of epithelial to mesenchymal transition

Epithelial to Mesenchymal Transition (EMT) is a cellular transformation process that is employed repeatedly and ubiquitously during vertebrate morphogenesis to build complex tissues and organs. Cellular transformations that occur during cancer cell invasion are phenotypically similar to developmental EMT, and involve the same molecular signalling pathways. EMT processes are diverse, but are characterised by: a loss of cell-cell adhesion; a gain in cell-matrix adhesion; an increase in cell motility; the secretion of proteases that degrade basement membrane proteins; an increased resistance to apoptosis; a loss of polarisation; increased production of extracellular matrix components; a change from a rounded to a fibroblastic morphology; and an invasive phenotype. This thesis focuses explicitly on endocardial EMT, which is the EMT that occurs during vertebrate embryonic heart development. The embryonic heart initially forms as a tube, with myocardium externally, endocardium internally, with these tissue layers separated by a thick extracellular matrix termed the cardiac jelly. Some of the endocardial cells in specific regions of the embryonic heart tube undergo EMT and invade the cardiac jelly. This causes cellularised swellings inside the embryonic heart tube termed the endocardial cushions. The emergence of the four chambered double pump heart of mammals involves a complex remodelling that the endocardial cushions play an active role in. Even while heart remodelling is taking place, the heart tube is operating as a single-circulation pump, and the endocardial cushions are performing a valve-like function that is critical to the survival of the embryo (Nomura-Kitabayashi et al. 2009). As the endocardial cushions grow and remodel, they become the valve leaflets of the foetal heart. The endocardial cushions also contribute tissue to the septa (walls) of the heart. Their correct formation is thus essential to the development of a fully functional, fully divided, double-pump system. It has been shown that genetic mutations that cause impaired endocardial EMT lead to the development of a range of congenital heart defects (Fischer et al. 2007). An extensive review is conducted of existing experimental investigations into endocardial EMT. The information extracted from this review is used to develop a multiscale conceptual model of endocardial EMT, including the major protein signalling pathways involved, and the cellular phenotypes that they induce or inhibit. After considering the requirements for computational simulations of EMT, and reviewing the various techniques and simulation packages available for multi-cell modelling, cellular Potts modelling is selected as having the most appropriate combination of features. The open source simulation platform Compucell3D is selected for model development, due to the flexibility, range of features provided and an existing implementation of multiscale models; that include subcellular models of reaction pathways. Based on the conceptual model of endocardial EMT, abstract computational simulations of key aspects are developed, in order to investigate qualitative behaviour under different simulated conditions. The abstract simulations include a 2D multiscale model of Notch signalling lateral induction, which is the mechanism by which the embryonic heart tube is patterned into cushion and non-cushion forming regions. Additionally, a 3D simulation is used to investigate the possible role of contact-inhibited mitosis, upregulated by the VEGF protein, in maintaining an epithelial phenotype. One particular in vitro investigation of endocardial EMT (Luna-Zurita et al. 2010) is used to develop quantitative simulations. The quantitative data used for fitting the simulations consist of cell shape metrics that are derived from simple processing of the imaging results. Single cell simulations are used to investigate the relationship between cell motility and cell shape in the cellular Potts model. The findings are then implemented in multi-cell models, in order to investigate the relationship between cell-cell adhesion, cell-matrix adhesion, cell motility and cell shape during EMT

Mechanics of Morphogenesis in Neural Development: in vivo, in vitro, and in silico

Morphogenesis in the central nervous system has received intensive attention as elucidating fundamental mechanisms of morphogenesis will shed light on the physiology and pathophysiology of the developing central nervous system. Morphogenesis of the central nervous system is of a vast topic that includes important morphogenetic events such as neurulation and cortical folding. Here we review three types of methods used to improve our understanding of morphogenesis of the central nervous system: in vivo experiments, organoids (in vitro), and computational models (in silico). The in vivo experiments are used to explore cellular- and tissue-level mechanics and interpret them on the roles of neurulation morphogenesis. Recent advances in human brain organoids have provided new opportunities to study morphogenesis and neurogenesis to compensate for the limitations of in vivo experiments, as organoid models are able to recapitulate some critical neural morphogenetic processes during early human brain development. Due to the complexity and costs of in vivo and in vitro studies, a variety of computational models have been developed and used to explain the formation and morphogenesis of brain structures. We review and discuss the Pros and Cons of these methods and their usage in the studies on morphogenesis of the central nervous system. Notably, none of these methods alone is sufficient to unveil the biophysical mechanisms of morphogenesis, thus calling for the interdisciplinary approaches using a combination of these methods in order to test hypotheses and generate new insights on both normal and abnormal development of the central nervous system

In Vitro Modeling of Mechanics in Cancer Metastasis

In addition to a multitude of genetic and biochemical alterations, abnormal morphological, structural, and mechanical changes in cells and their extracellular environment are key features of tumor invasion and metastasis. Furthermore, it is now evident that mechanical cues alongside biochemical signals contribute to critical steps of cancer initiation, progression, and spread. Despite its importance, it is very challenging to study mechanics of different steps of metastasis in the clinic or even in animal models. While considerable progress has been made in developing advanced in vitro models for studying genetic and biological aspects of cancer, less attention has been paid to models that can capture both biological and mechanical factors realistically. This is mainly due to lack of appropriate models and measurement tools. After introducing the central role of mechanics in cancer metastasis, we provide an outlook on the emergence of novel in vitro assays and their combination with advanced measurement technologies to probe and recapitulate mechanics in conditions more relevant to the metastatic disease

A strategy for tissue self-organization that is robust to cellular heterogeneity and plasticity

Developing tissues contain motile populations of cells that can self-organize into spatially ordered tissues based on differences in their interfacial surface energies. However, it is unclear how self-organization by this mechanism remains robust when interfacial energies become heterogeneous in either time or space. The ducts and acini of the human mammary gland are prototypical heterogeneous and dynamic tissues comprising two concentrically arranged cell types. To investigate the consequences of cellular heterogeneity and plasticity on cell positioning in the mammary gland, we reconstituted its self-organization from aggregates of primary cells in vitro. We find that self-organization is dominated by the interfacial energy of the tissue–ECM boundary, rather than by differential homo- and heterotypic energies of cell–cell interaction. Surprisingly, interactions with the tissue–ECM boundary are binary, in that only one cell type interacts appreciably with the boundary. Using mathematical modeling and cell-type-specific knockdown of key regulators of cell–cell cohesion, we show that this strategy of self-organization is robust to severe perturbations affecting cell–cell contact formation. We also find that this mechanism of self-organization is conserved in the human prostate. Therefore, a binary interfacial interaction with the tissue boundary provides a flexible and generalizable strategy for forming and maintaining the structure of two-component tissues that exhibit abundant heterogeneity and plasticity. Our model also predicts that mutations affecting binary cell–ECM interactions are catastrophic and could contribute to loss of tissue architecture in diseases such as breast cancer

Chaste: a test-driven approach to software development for biological modelling

Chaste (‘Cancer, heart and soft-tissue environment’) is a software library and a set of test suites for computational simulations in the domain of biology. Current functionality has arisen from modelling in the fields of cancer, cardiac physiology and soft-tissue mechanics. It is released under the LGPL 2.1 licence.\ud \ud Chaste has been developed using agile programming methods. The project began in 2005 when it was reasoned that the modelling of a variety of physiological phenomena required both a generic mathematical modelling framework, and a generic computational/simulation framework. The Chaste project evolved from the Integrative Biology (IB) e-Science Project, an inter-institutional project aimed at developing a suitable IT infrastructure to support physiome-level computational modelling, with a primary focus on cardiac and cancer modelling

How Changes in Extracellular Matrix Mechanics and Gene Expression Variability Might Combine to Drive Cancer Progression

Changes in extracellular matrix (ECM) structure or mechanics can actively drive cancer progression; however, the underlying mechanism remains unknown. Here we explore whether this process could be mediated by changes in cell shape that lead to increases in genetic noise, given that both factors have been independently shown to alter gene expression and induce cell fate switching. We do this using a computer simulation model that explores the impact of physical changes in the tissue microenvironment under conditions in which physical deformation of cells increases gene expression variability among genetically identical cells. The model reveals that cancerous tissue growth can be driven by physical changes in the microenvironment: when increases in cell shape variability due to growth-dependent increases in cell packing density enhance gene expression variation, heterogeneous autonomous growth and further structural disorganization can result, thereby driving cancer progression via positive feedback. The model parameters that led to this prediction are consistent with experimental measurements of mammary tissues that spontaneously undergo cancer progression in transgenic C3(1)-SV40Tag female mice, which exhibit enhanced stiffness of mammary ducts, as well as progressive increases in variability of cell-cell relations and associated cell shape changes. These results demonstrate the potential for physical changes in the tissue microenvironment (e.g., altered ECM mechanics) to induce a cancerous phenotype or accelerate cancer progression in a clonal population through local changes in cell geometry and increased phenotypic variability, even in the absence of gene mutation

The effects of ageing on microenvironment-contextual epithelial cell signalling

Age is the greatest risk factor for developing cancer. Two-thirds of cancer diagnoses occur in people over the age of 65. This increased vulnerability to tumourigenesis is attributed to intrinsic cellular changes, in particular the agerelated telomere shortening and the accumulation of mutations over time and. We propose that extrinsic factors comprising age-related alterations in the tissue microenvironment are also important in cancer development. We sought to elucidate how the microenvironmental affects cells, and how this is related to cancer development. This knowledge can be utilized to improve prevention, diagnosis and treatment of cancer in our ageing population. Cellular function is coordinated by microenvironmental factors such as growth factors, cytokines as well as extracellular matrix proteins. Soluble factormediated signal transduction is strongly influenced by microenvironmental context. To allow single cell level measurement of the microenvironmental contextual effect on cell signaling, we developed a novel flow cytometry method: microsphere cytometry. Single normal or neoplastic cells were adhered to uniform microspheres that display mimetic-microenvironments comprising surface combinations of extracellular matrix (ECM) in the presence of soluble agonists or antagonists. Temporal signaling responses were measured with fluorophore-conjugated antibodies that recognize responsedependent epitopes by multiparametric flow cytometry. Using this approach, we demonstrated that microenvironment-mimetic combinations of growth factors and extracellular matrix proteins generate distinct cellular signaling signatures from normal and patient biopsy-derived neoplastic cells. We asked whether the ageing process affects how human mammary epithelial cells (HMEC) respond to microenvironmental signals, and if this altered response results in increased susceptibility to oncogenic transformation. Using microsphere cytometry we analyzed age-dependent changes in human mammary myoepithelial and luminal epithelial cells exposed to different ECM and growth factors. We found that ECM–mediated MAP kinase and PI3K pathway activation levels in HMEC are attenuated with age, and that the diminished signaling magnitude in HMEC from ageing women correlated with reduced probability of activating oncogene-induced senescence. Our results suggest that attenuated cell signaling response may reduce the likelihood of activating oncogene induced senescence, for cells in ageing women. We hypothesize this is the result of age-related changes to the microenvironment that support age-emergent cellular phenotypes with increased cancer susceptibility