A theoretical investigation of the effect of proliferation & adhesion on monoclonal conversion in the colonic crypt

The surface epithelium lining the intestinal tract renews itself rapidly by a coordinated programme of cell proliferation, migration and differentiation events that is initiated in the crypts of Lieberkühn. It is generally believed that colorectal cancer arises due to mutations that disrupt the normal cellular dynamics of the crypts. Using a spatially structured cell-based model of a colonic crypt, we investigate the likelihood that the progeny of a mutated cell will dominate, or be sloughed out of, a crypt. Our approach is to perform multiple simulations, varying the spatial location of the initial mutation, and the proliferative and adhesive properties of the mutant cells, to obtain statistical distributions for the probability of their domination. Our simulations lead us to make a number of predictions. The process of monoclonal conversion always occurs, and does not require that the cell which initially gave rise to the population remains in the crypt. Mutations occurring more than one to two cells from the base of the crypt are unlikely to become the dominant clone. The probability of a mutant clone persisting in the crypt is sensitive to dysregulation of adhesion. By comparing simulation results with those from a simple one-dimensional stochastic model of population dynamics at the base of the crypt, we infer that this sensitivity is due to direct competition between wild-type and mutant cells at the base of the crypt. We also predict that increases in the extent of the spatial domain in which the mutant cells proliferate can give rise to counter-intuitive, non-linear changes to the probability of their fixation, due to effects that cannot be captured in simpler models

A theoretical investigation of the effect of proliferation and\ud adhesion on monoclonal conversion in the colonic crypt

Colorectal cancers are initiated by the accumulation of mutations in the colonic epithelium. Using a spatially structured cell-based model of a colonic crypt, we investigate the likelihood that the progeny of a mutated cell will dominate, or be sloughed out of, a crypt. Our approach is to perform multiple simulations, varying the spatial location of the initial mutation, and its proliferative and adhesive properties, to obtain statistical distributions for the probability of domination. Our simulations lead us to make a number of predictions. The process of monoclonal conversion always occurs, and does not require that the cell which initially gave rise to the population remains in the crypt. Mutations occurring more than one to two cells from the base of the crypt are unlikely to become the dominant clone. The probability of a mutant clone persisting in the crypt is sensitive to dysregulation of adhesion, and comparison with a one-dimensional model suggests that this is caused by competition directly at the base of the crypt.\ud We also predict that increases in the extent of the spatial domain in which the mutant cells proliferate cause counter-intuitive non-linear changes to the probability of its fixation, due to effects that cannot be captured in simpler models

Modelling spatially regulated B-catenin dynamics & invasion in intestinal crypts

Experimental data (e.g., genetic lineage and cell population studies) on intestinal crypts reveal that regulatory features of crypt behavior, such as control via morphogen gradients, are remarkably well conserved among numerous organisms (e.g., from mouse and rat to human) and throughout the different regions of the small and large intestines. In this article, we construct a partial differential equation model of a single colonic crypt that describes the spatial distribution of Wnt pathway proteins along the crypt axis. The novelty of our continuum model is that it is based upon assumptions that can be directly related to processes at the cellular and subcellular scales. We use the model to predict how the distributions of Wnt pathway proteins are affected by mutations. The model is then extended to investigate how mutant cell populations can invade neighboring crypts. The model simulations suggest that cell crowding caused by increased proliferation and decreased cell loss may be sufficient for a mutant cell population to colonize a neighboring healthy crypt

Chaste: an open source C++ library for computational physiology and biology

Chaste - Cancer, Heart And Soft Tissue Environment - is an open source C++ library for the computational simulation of mathematical models developed for physiology and biology. Code development has been driven by two initial applications: cardiac electrophysiology and cancer development. A large number of cardiac electrophysiology studies have been enabled and performed, including high performance computational investigations of defibrillation on realistic human cardiac geometries. New models for the initiation and growth of tumours have been developed. In particular, cell-based simulations have provided novel insight into the role of stem cells in the colorectal crypt. Chaste is constantly evolving and is now being applied to a far wider range of problems. The code provides modules for handling common scientific computing components, such as meshes and solvers for ordinary and partial differential equations (ODEs/PDEs). Re-use of these components avoids the need for researchers to "re-invent the wheel" with each new project, accelerating the rate of progress in new applications. Chaste is developed using industrially-derived techniques, in particular test-driven development, to ensure code quality, re-use and reliability. In this article we provide examples that illustrate the types of problems Chaste can be used to solve, which can be run on a desktop computer. We highlight some scientific studies that have used or are using Chaste, and the insights they have provided. The source code, both for specific releases and the development version, is available to download under an open source Berkeley Software Distribution (BSD) licence at http://www.cs.ox.ac.uk/chaste, together with details of a mailing list and links to documentation and tutorials

Accounting for variability in ion current recordings using a mathematical model of artefacts in voltage-clamp experiments

Mathematical models of ion channels, which constitute indispensable components of action potential models, are commonly constructed by fitting to whole-cell patch-clamp data. In a previous study, we fitted cell-specific models to hERG1a (Kv11.1) recordings simultaneously measured using an automated high-throughput system, and studied cell-cell variability by inspecting the resulting model parameters. However, the origin of the observed variability was not identified. Here, we study the source of variability by constructing a model that describes not just ion current dynamics, but the entire voltage-clamp experiment. The experimental artefact components of the model include: series resistance, membrane and pipette capacitance, voltage offsets, imperfect compensations made by the amplifier for these phenomena, and leak current. In this model, variability in the observations can be explained by either cell properties, measurement artefacts, or both. Remarkably, by assuming that variability arises exclusively from measurement artefacts, it is possible to explain a larger amount of the observed variability than when assuming cell-specific ion current kinetics. This assumption also leads to a smaller number of model parameters. This result suggests that most of the observed variability in patch-clamp data measured under the same conditions is caused by experimental artefacts, and hence can be compensated for in post-processing by using our model for the patch-clamp experiment. This study has implications for the question of the extent to which cell-cell variability in ion channel kinetics exists, and opens up routes for better correction of artefacts in patch-clamp data. This article is part of the theme issue 'Uncertainty quantification in cardiac and cardiovascular modelling and simulation'

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

A nonlinear and time-dependent leak current in the presence of calcium fluoride patch-clamp seal enhancer [version 1; peer review: 2 approved with reservations]

Automated patch-clamp platforms are widely used and vital tools in both academia and industry to enable high-throughput studies such as drug screening. A leak current to ground occurs whenever the seal between a pipette and cell (or internal solution and cell in high-throughput machines) is not perfectly insulated from the bath (extracellular) solution. Over 1 GΩ seal resistance between pipette and bath solutions is commonly used as a quality standard for manual patch work. With automated platforms it can be difficult to obtain such a high seal resistance between the intra- and extra-cellular solutions. One suggested method to alleviate this problem is using an F− containing internal solution together with a Ca2+ containing external solution — so that a CaF2 crystal forms when the two solutions meet which ‘plugs the holes’ to enhance the seal resistance. However, we observed an unexpected nonlinear-in-voltage and time-dependent current using these solutions on an automated patch-clamp platform. We performed manual patch-clamp experiments with the automated patch-clamp solutions, but no biological cell, and observed the same nonlinear time-dependent leak current. The current could be completely removed by washing out F− ions to leave a conventional leak current that was linear and not time-dependent. We therefore conclude fluoride ions interacting with the CaF2 crystal are the origin of the nonlinear time-dependent leak current. The consequences of such a nonlinear and time-dependent leak current polluting measurements should be considered carefully if it cannot be isolated and subtracted

An integrative computational model for intestinal tissue renewal

Objectives\ud \ud The luminal surface of the gut is lined with a monolayer of epithelial cells that acts as a nutrient absorptive engine and protective barrier. To maintain its integrity and functionality, the epithelium is renewed every few days. Theoretical models are powerful tools that can be used to test hypotheses concerning the regulation of this renewal process, to investigate how its dysfunction can lead to loss of homeostasis and neoplasia, and to identify potential therapeutic interventions. Here we propose a new multiscale model for crypt dynamics that links phenomena occurring at the subcellular, cellular and tissue levels of organisation.\ud \ud Methods\ud \ud At the subcellular level, deterministic models characterise molecular networks, such as cell-cycle control and Wnt signalling. The output of these models determines the behaviour of each epithelial cell in response to intra-, inter- and extracellular cues. The modular nature of the model enables us to easily modify individual assumptions and analyse their effects on the system as a whole.\ud \ud Results\ud \ud We perform virtual microdissection and labelling-index experiments, evaluate the impact of various model extensions, obtain new insight into clonal expansion in the crypt, and compare our predictions with recent mitochondrial DNA mutation data. \ud \ud Conclusions\ud \ud We demonstrate that relaxing the assumption that stem-cell positions are fixed enables clonal expansion and niche succession to occur. We also predict that the presence of extracellular factors near the base of the crypt alone suffices to explain the observed spatial variation in nuclear beta-catenin levels along the crypt axis

Empirical Quantification of Predictive Uncertainty Due to Model Discrepancy by Training with an Ensemble of Experimental Designs: An Application to Ion Channel Kinetics

When using mathematical models to make quantitative predictions for clinical or industrial use, it is important that predictions come with a reliable estimate of their accuracy (uncertainty quantification). Because models of complex biological systems are always large simplifications, model discrepancy arises—models fail to perfectly recapitulate the true data generating process. This presents a particular challenge for making accurate predictions, and especially for accurately quantifying uncertainty in these predictions. Experimentalists and modellers must choose which experimental procedures (protocols) are used to produce data used to train models. We propose to characterise uncertainty owing to model discrepancy with an ensemble of parameter sets, each of which results from training to data from a different protocol. The variability in predictions from this ensemble provides an empirical estimate of predictive uncertainty owing to model discrepancy, even for unseen protocols. We use the example of electrophysiology experiments that investigate the properties of hERG potassium channels. Here, ‘information-rich’ protocols allow mathematical models to be trained using numerous short experiments performed on the same cell. In this case, we simulate data with one model and fit it with a different (discrepant) one. For any individual experimental protocol, parameter estimates vary little under repeated samples from the assumed additive independent Gaussian noise model. Yet parameter sets arising from the same model applied to different experiments conflict—highlighting model discrepancy. Our methods will help select more suitable ion channel models for future studies, and will be widely applicable to a range of biological modelling problems