Computational model of cell positioning: directed and collective migration in the intestinal crypt epithelium

Chitosan nanocomposite for tissue engineering and regenerative medicine: A review

Regenerative medicine and tissue engineering have emerged as a multidisciplinary promising field in the quest to address the limitations of traditional medical approaches. One of the key aspects of these fields is the development of such types of biomaterials that can mimic the extracellular matrix and provide a conducive environment for tissue regeneration. In this regard, chitosan has played a vital role which is a naturally derived linear bi-poly-aminosaccharide, and has gained significant attention due to its biocompatibility and unique properties. Chitosan possesses many unique physicochemical properties, making it a significant polysaccharide for different applications such as agriculture, nutraceutical, biomedical, food, nutraceutical, packaging, etc. as well as significant material for developing next-generation hydrogel and bio-scaffolds for regenerative medicinal applications. Moreover, chitosan can be easily modified to incorporate desirable properties, such as improved mechanical strength, enhanced biodegradability, and controlled release of bioactive molecules. Blending chitosan with other polymers or incorporating nanoparticles into its matrix further expands its potential in tissue engineering applications. This review summarizes the most recent studies of the last 10 years based on chitosan, blends, and nanocomposites and their application in bone tissue engineering, hard tissue engineering, dental implants, dental tissue engineering, dental fillers, and cartilage tissue engineering.

Mechano-regulated cell-cell signaling in the context of cardiovascular tissue engineering

Cardiovascular tissue engineering (CVTE) aims to create living tissues, with the ability to grow and remodel, as replacements for diseased blood vessels and heart valves. Despite promising results, the (long-term) functionality of these engineered tissues still needs improvement to reach broad clinical application. The functionality of native tissues is ensured by their specific mechanical properties directly arising from tissue organization. We therefore hypothesize that establishing a native-like tissue organization is vital to overcome the limitations of current CVTE approaches. To achieve this aim, a better understanding of the growth and remodeling (G&R) mechanisms of cardiovascular tissues is necessary. Cells are the main mediators of tissue G&R, and their behavior is strongly influenced by both mechanical stimuli and cell-cell signaling. An increasing number of signaling pathways has also been identified as mechanosensitive. As such, they may have a key underlying role in regulating the G&R of tissues in response to mechanical stimuli. A more detailed understanding of mechano-regulated cell-cell signaling may thus be crucial to advance CVTE, as it could inspire new methods to control tissue G&R and improve the organization and functionality of engineered tissues, thereby accelerating clinical translation. In this review, we discuss the organization and biomechanics of native cardiovascular tissues; recent CVTE studies emphasizing the obtained engineered tissue organization; and the interplay between mechanical stimuli, cell behavior, and cell-cell signaling. In addition, we review past contributions of computational models in understanding and predicting mechano-regulated tissue G&R and cell-cell signaling to highlight their potential role in future CVTE strategies.

Intestinal villus structure contributes to even shedding of epithelial cells

In the intestinal epithelium, proliferated epithelial cells ascend the crypts and villi and shed at the villus tips into the gut lumen. In this study, we theoretically investigate the roles of the villi on cell turnover. We present a stochastic model that focuses on the duration over which cells migrate the shortest paths between the crypt orifices and the villus tips, where shedding cells are randomly chosen from among those older than the shortest-path cell migration times. By extending the length of the shortest path to delay cell shedding, the finger-like shape of the villus would tightly regulate shedding-cell ages compared with flat surfaces and shorter projections; the villus allows epithelial cells to shed at around the same age, which limits them from shedding early or staying in the epithelium for long periods. Computational simulations of cell dynamics agreed well with the predictions. We also examine various mechanical conditions of cells and confirm that coordinated collective cell migration supports the predictions. These results suggest the important roles of the villi in homeostatic maintenance of the small intestine, and we discuss the applicability of our approach to other tissues with collective cell movement.

Multivalent Binding of a Ligand-Coated Particle: Role of Shape, Size, and Ligand Heterogeneity

We utilize a multiscale modeling framework to study the effect of shape, size, and ligand composition on the efficacy of binding of a ligand-coated particle to a substrate functionalized with the target receptors. First, we show how molecular dynamics along with steered molecular dynamics calculations can be used to accurately parameterize the molecular-binding free energy and the effective spring constant for a receptor-ligand pair. We demonstrate this for two ligands that bind to the α₅β₁-domain of integrin. Next, we show how these effective potentials can be used to build computational models at the meso- and continuum-scales. These models incorporate the molecular nature of the receptor-ligand interactions and yet provide an inexpensive route to study the multivalent interaction of receptors and ligands through the construction of Bell potentials customized to the molecular identities. We quantify the binding efficacy of the ligand-coated-particle in terms of its multivalency, binding free-energy landscape, and the losses in the configurational entropies. We show that 1) the binding avidity for particle sizes less than 350 nm is set by the competition between the enthalpic and entropic contributions, whereas that for sizes above 350 nm is dominated by the enthalpy of binding; 2) anisotropic particles display higher levels of multivalent binding compared to those of spherical particles; and 3) variations in ligand composition can alter binding avidity without altering the average multivalency. The methods and results presented here have wide applications in the rational design of functionalized carriers and also in understanding cell adhesion.

Cell-Cell Adhesion and Cortical Actin Bending Govern Cell Elongation on Negatively Curved Substrates

Physiologically, cells experience and respond to a variety of mechanical stimuli such as rigidity and topography of the extracellular matrix. However, little is known about the effects of substrate curvature on cell behavior. We developed a novel, to our knowledge, method to fabricate cell culture substrates with semicylindrical grooves of negative curvatures (radius of curvature, R_c = 20-100 μm). We found that negative substrate curvatures induced elongation of mesenchymal and epithelial cells along the cylinder axis. As R_c decreases, mesenchymal National Institutes of Health 3T3 fibroblasts increasingly elongate along the long axis of the grooves, whereas elongation of epithelial Madin-Darby Canine Kidney (MDCK) cells is biphasic with maximal cell elongation when R_c = 40 μm. Addition of blebbistatin to MDCK cells to reduce cortical actin rigidity resulted in a decrease in cell elongation across all curvatures while preserving the biphasic trend. However, addition of calyculin A or ethylene glycol-bis(2-aminoethylether)-N,N,N',N'-tetraacetic acid, to increase cortical rigidity or reduce intercellular adhesion, respectively, resulted in a monotonic increase in MDCK cell elongation with decreasing R_c. Using an energy minimization model, we showed that cell elongation in epithelial cell sheet is governed by the competition between two energies as R_c decreases: curvature-dependent intercellular adhesion that prevents elongation; and intracellular cortical actin bending that enhances elongation. Therefore, our results of cellular elongation induced by negatively curved substrates offer insights into how tubule elongation or growth of tubular structures such as kidney tubules can be controlled by the substrate curvature in vivo.

Stem cell competition in the gut: insights from multi-scale computational modelling

Three-dimensional (3D) computational tissue models can provide a comprehensive description of tissue dynamics at the molecular, cellular and tissue level. Moreover, they can support the development of hypotheses about cellular interactions and about synergies between major signalling pathways. We exemplify these capabilities by simulation of a 3D single-cell-based model of mouse small intestinal crypts. We analyse the impact of lineage specification, distribution and cellular lifespan on clonal competition and study effects of Notch- and Wnt activation on fixation of mutations within the tissue. Based on these results, we predict that experimentally observed synergistic effects between autonomous Notch- and Wnt signalling in triggering intestinal tumourigenesis originate in the suppression of Wnt-dependent secretory lineage specification by Notch, giving rise to an increased fixation probability of Wnt-activating mutations. Our study demonstrates that 3D computational tissue models can support a mechanistic understanding of long-term tissue dynamics under homeostasis and during transformation.

The role of apical contractility in determining cell morphology in multilayered epithelial sheets and tubes

A multilayered epithelium is made up of individual cells that are stratified in an orderly fashion, layer by layer. In such tissues, individual cells can adopt a wide range of shapes ranging from columnar to squamous. From histological images, we observe that, in flat epithelia such as the skin, the cells in the top layer are squamous while those in the middle and bottom layers are columnar, whereas in tubular epithelia, the cells in all layers are columnar. We develop a computational model to understand how individual cell shape is governed by the mechanical forces within multilayered flat and curved epithelia. We derive the energy function for an epithelial sheet of cells considering intercellular adhesive and intracellular contractile forces. We determine computationally the cell morphologies that minimize the energy function for a wide range of cellular parameters. Depending on the dominant adhesive and contractile forces, we find four dominant cell morphologies for the multilayered-layered flat sheet and three dominant cell morphologies for the two-layered curved sheet. We study the transitions between the dominant cell morphologies for the two-layered flat sheet and find both continuous and discontinuous transitions and also the presence of multistable states. Matching our computational results with histological images, we conclude that apical contractile forces from the actomyosin belt in the epithelial cells is the dominant force determining cell shape in multilayered epithelia. Our computational model can guide tissue engineers in designing artificial multilayered epithelia, in terms of figuring out the cellular parameters needed to achieve realistic epithelial morphologies.

Deciphering the Complex Signaling Systems That Regulate Intestinal Epithelial Cell Death Processes and Shedding

Intestinal epithelial cells play a fundamental role in maintaining homeostasis. Shedding of intestinal cells in a controlled manner is critical to maintenance of barrier function. Barrier function is maintained during this shedding process by a redistribution of tight junctional proteins to facilitate closure of the gap left by the shedding cell. However, despite the obvious importance of epithelial cell shedding to gut health, a central question is how the extrusion of epithelial cells is achieved, enabling barrier integrity to be maintained in the healthy gut and restored during inflammation remains largely unanswered. Recent studies have provided evidence that excessive epithelial cell shedding and loss of epithelial barrier integrity is triggered by exposure to lipopolysaccharide or tumor necrosis factor alpha. Subsequent studies have provided evidence of the involvement of specific cellular components and signaling mechanisms as well as the functionality of microbiota that can be either detrimental or beneficial for intestinal barrier integrity. This review will focus on the evidence and decipher how the signaling systems through which the mucosal immune system and microbiota can regulate epithelial cell shedding and how these mechanisms interact to preserve the viability of the epithelium.

Spatio-temporal morphology changes in and quenching effects on the 2D spreading dynamics of cell colonies in both plain and methylcellulose-containing culture media

To deal with complex systems, microscopic and global approaches become of particular interest. Our previous results from the dynamics of large cell colonies indicated that their 2D front roughness dynamics is compatible with the standard Kardar-Parisi-Zhang (KPZ) or the quenched KPZ equations either in plain or methylcellulose (MC)-containing gel culture media, respectively. In both cases, the influence of a non-uniform distribution of the colony constituents was significant. These results encouraged us to investigate the overall dynamics of those systems considering the morphology and size, the duplication rate, and the motility of single cells. For this purpose, colonies with different cell populations (N) exhibiting quasi-circular and quasi-linear growth fronts in plain and MC-containing culture media are investigated. For small N, the average radial front velocity and its change with time depend on MC concentration. MC in the medium interferes with cell mitosis, contributes to the local enlargement of cells, and increases the distribution of spatio-temporal cell density heterogeneities. Colony spreading in MC-containing media proceeds under two main quenching effects, I and II; the former mainly depending on the culture medium composition and structure and the latter caused by the distribution of enlarged local cell domains. For large N, colony spreading occurs at constant velocity. The characteristics of cell motility, assessed by measuring their trajectories and the corresponding velocity field, reflect the effect of enlarged, slow-moving cells and the structure of the medium. Local average cell size distribution and individual cell motility data from plain and MC-containing media are qualitatively consistent with the predictions of both the extended cellular Potts models and the observed transition of the front roughness dynamics from a standard KPZ to a quenched KPZ. In this case, quenching effects I and II cooperate and give rise to the quenched-KPZ equation. Seemingly, these results show a possible way of linking the cellular Potts models and the 2D colony front roughness dynamics.

Multiscale Model of Colorectal Cancer Using the Cellular Potts Framework

Colorectal cancer (CRC) is one of the major causes of death in the developed world and forms a canonical example of tumorigenesis. CRC arises from a string of mutations of individual cells in the colorectal crypt, making it particularly suited for multiscale multicellular modeling, where mutations of individual cells can be clearly represented and their effects readily tracked. In this paper, we present a multicellular model of the onset of colorectal cancer, utilizing the cellular Potts model (CPM). We use the model to investigate how, through the modification of their mechanical properties, mutant cells colonize the crypt. Moreover, we study the influence of mutations on the shape of cells in the crypt, suggesting possible cell- and tissue-level indicators for identifying early-stage cancerous crypts. Crucially, we discuss the effect that the motility parameters of the model (key factors in the behavior of the CPM) have on the distribution of cells within a homeostatic crypt, resulting in an optimal parameter regime that accurately reflects biological assumptions. In summary, the key results of this paper are 1) how to couple the CPM with processes occurring on other spatial scales, using the example of the crypt to motivate suitable motility parameters; 2) modeling mutant cells with the CPM; 3) and investigating how mutations influence the shape of cells in the crypt.

CoGNaC: A Chaste Plugin for the Multiscale Simulation of Gene Regulatory Networks Driving the Spatial Dynamics of Tissues and Cancer

We introduce a Chaste plugin for the generation and the simulation of Gene Regulatory Networks (GRNs) in multiscale models of multicellular systems. Chaste is a widely used and versatile computational framework for the multiscale modeling and simulation of multicellular biological systems. The plugin, named CoGNaC (Chaste and Gene Networks for Cancer), allows the linking of the regulatory dynamics to key properties of the cell cycle and of the differentiation process in populations of cells, which can subsequently be modeled using different spatial modeling scenarios. The approach of CoGNaC focuses on the emergent dynamical behavior of gene networks, in terms of gene activation patterns characterizing the different cellular phenotypes of real cells and, especially, on the overall robustness to perturbations and biological noise. The integration of this approach within Chaste’s modular simulation framework provides a powerful tool to model multicellular systems, possibly allowing for the formulation of novel hypotheses on gene regulation, cell differentiation, and, in particular, cancer emergence and development. In order to demonstrate the usefulness of CoGNaC over a range of modeling paradigms, two example applications are presented. The first of these concerns the characterization of the gene activation patterns of human T-helper cells. The second example is a multiscale simulation of a simplified intestinal crypt, in which, given certain conditions, tumor cells can emerge and colonize the tissue.

The Interplay between Wnt Mediated Expansion and Negative Regulation of Growth Promotes Robust Intestinal Crypt Structure and Homeostasis

The epithelium of the small intestinal crypt, which has a vital role in protecting the underlying tissue from the harsh intestinal environment, is completely renewed every 4–5 days by a small pool of stem cells at the base of each crypt. How is this renewal controlled and homeostasis maintained, particularly given the rapid nature of this process? Here, based on the recent observations from in vitro “mini gut” studies, we use a hybrid stochastic model of the crypt to investigate how exogenous niche signaling (from Wnt and BMP) combines with auto-regulation to promote homeostasis. This model builds on the sub-cellular element method to account for the three-dimensional structure of the crypt, external regulation by Wnt and BMP, internal regulation by Notch signaling, as well as regulation by internally generated diffusible signals. Results show that Paneth cell derived Wnt signals, which have been observed experimentally to sustain crypts in cultured organs, have a dramatically different influence on niche dynamics than does mesenchyme derived Wnt. While this signaling can indeed act as a redundant backup to the exogenous gradient, it introduces a positive feedback that destabilizes the niche and causes its uncontrolled expansion. We find that in this setting, BMP has a critical role in constraining this expansion, consistent with observations that its removal leads to crypt fission. Further results also point to a new hypothesis for the role of Ephrin mediated motility of Paneth cells, specifically that it is required to constrain niche expansion and maintain the crypt’s spatial structure. Combined, these provide an alternative view of crypt homeostasis where the niche is in a constant state of expansion and the spatial structure of the crypt arises as a balance between this expansion and the action of various sources of negative regulation that hold it in check.

The small intestinal epithelium, like our skin, is constantly being renewed. In the intestine however, this epithelium is exposed to the harsh digestive environment, necessitating much more rapid renewal. Remarkably, the entire epithelium is renewed every 4–5 days. This raises the question, how can the size and structure of this tissue be maintained given this pace. Motivated by recent experimental observations, we construct a three-dimensional, hybrid stochastic model to investigate the mechanisms responsible for homeostasis of this tissue. We find that there are redundant signals created by both the epithelium itself and surrounding tissues that act in parallel to maintain epithelial structure. This redundancy comes at a price however: it introduces the possibility of explosive stem cell population growth. Additional results suggest that other signals along with choreographed motion of cells are responsible for repressing this expansion. Taken together, our results provide a novel hypothesis for how robust but fast renewal of the crypt is achieved: as a balance between expansion, which drives fast renewal and repression, which holds that expansion in check to maintain the crypt’s structure.

A mathematical model for eph/ephrin-directed segregation of intermingled cells

Eph receptors, the largest family of receptor tyrosine kinases, control cell-cell adhesion/de-adhesion, cell morphology and cell positioning through interaction with cell surface ephrin ligands. Bi-directional signalling from the Eph and ephrin complexes on interacting cells have a significant role in controlling normal tissue development and oncogenic tissue patterning. Eph-mediated tissue patterning is based on the fine-tuned balance of adhesion and de-adhesion reactions between distinct Eph- and ephrin-expressing cell populations, and adhesion within like populations (expressing either Eph or ephrin). Here we develop a stochastic, Lagrangian model that is based on Eph/ephrin biology: incorporating independent Brownian motion to describe cell movement and a deterministic term (the drift term) to represent repulsive and adhesive interactions between neighbouring cells. Comparison between the experimental and computer simulated Eph/ephrin cell patterning events shows that the model recapitulates the dynamics of cell-cell segregation and cell cluster formation. Moreover, by modulating the term for Eph/ephrin-mediated repulsion, the model can be tuned to match the actual behaviour of cells with different levels of Eph expression or activity. Together the results of our experiments and modelling suggest that the complexity of Eph/ephrin signalling mechanisms that control cell-cell interactions can be described well by a mathematical model with a single term balancing adhesion and de-adhesion between interacting cells. This model allows reliable prediction of Eph/ephrin-dependent control of cell patterning behaviour.

Computational modeling reveals that a combination of chemotaxis and differential adhesion leads to robust cell sorting during tissue patterning

Robust tissue patterning is crucial to many processes during development. The "French Flag" model of patterning, whereby naïve cells in a gradient of diffusible morphogen signal adopt different fates due to exposure to different amounts of morphogen concentration, has been the most widely proposed model for tissue patterning. However, recently, using time-lapse experiments, cell sorting has been found to be an alternative model for tissue patterning in the zebrafish neural tube. But it remains unclear what the sorting mechanism is. In this article, we used computational modeling to show that two mechanisms, chemotaxis and differential adhesion, are needed for robust cell sorting. We assessed the performance of each of the two mechanisms by quantifying the fraction of correct sorting, the fraction of stable clusters formed after correct sorting, the time needed to achieve correct sorting, and the size variations of the cells having different fates. We found that chemotaxis and differential adhesion confer different advantages to the sorting process. Chemotaxis leads to high fraction of correct sorting as individual cells will either migrate towards or away from the source depending on its cell type. However after the cells have sorted correctly, there is no interaction among cells of the same type to stabilize the sorted boundaries, leading to cell clusters that are unstable. On the other hand, differential adhesion results in low fraction of correct clusters that are more stable. In the absence of morphogen gradient noise, a combination of both chemotaxis and differential adhesion yields cell sorting that is both accurate and robust. However, in the presence of gradient noise, the simple combination of chemotaxis and differential adhesion is insufficient for cell sorting; instead, chemotaxis coupled with delayed differential adhesion is required to yield optimal sorting.

Investigating the relation between stochastic differentiation, homeostasis and clonal expansion in intestinal crypts via multiscale modeling

Colorectal tumors originate and develop within intestinal crypts. Even though some of the essential phenomena that characterize crypt structure and dynamics have been effectively described in the past, the relation between the differentiation process and the overall crypt homeostasis is still only partially understood. We here investigate this relation and other important biological phenomena by introducing a novel multiscale model that combines a morphological description of the crypt with a gene regulation model: the emergent dynamical behavior of the underlying gene regulatory network drives cell growth and differentiation processes, linking the two distinct spatio-temporal levels. The model relies on a few a priori assumptions, yet accounting for several key processes related to crypt functioning, such as: dynamic gene activation patterns, stochastic differentiation, signaling pathways ruling cell adhesion properties, cell displacement, cell growth, mitosis, apoptosis and the presence of biological noise. We show that this modeling approach captures the major dynamical phenomena that characterize the regular physiology of crypts, such as cell sorting, coordinate migration, dynamic turnover, stem cell niche correct positioning and clonal expansion. All in all, the model suggests that the process of stochastic differentiation might be sufficient to drive the crypt to homeostasis, under certain crypt configurations. Besides, our approach allows to make precise quantitative inferences that, when possible, were matched to the current biological knowledge and it permits to investigate the role of gene-level perturbations, with reference to cancer development. We also remark the theoretical framework is general and may be applied to different tissues, organs or organisms.

Investigation of inflammation and tissue patterning in the gut using a Spatially Explicit General-purpose Model of Enteric Tissue (SEGMEnT)

The mucosa of the intestinal tract represents a finely tuned system where tissue structure strongly influences, and is turn influenced by, its function as both an absorptive surface and a defensive barrier. Mucosal architecture and histology plays a key role in the diagnosis, characterization and pathophysiology of a host of gastrointestinal diseases. Inflammation is a significant factor in the pathogenesis in many gastrointestinal diseases, and is perhaps the most clinically significant control factor governing the maintenance of the mucosal architecture by morphogenic pathways. We propose that appropriate characterization of the role of inflammation as a controller of enteric mucosal tissue patterning requires understanding the underlying cellular and molecular dynamics that determine the epithelial crypt-villus architecture across a range of conditions from health to disease. Towards this end we have developed the Spatially Explicit General-purpose Model of Enteric Tissue (SEGMEnT) to dynamically represent existing knowledge of the behavior of enteric epithelial tissue as influenced by inflammation with the ability to generate a variety of pathophysiological processes within a common platform and from a common knowledge base. In addition to reproducing healthy ileal mucosal dynamics as well as a series of morphogen knock-out/inhibition experiments, SEGMEnT provides insight into a range of clinically relevant cellular-molecular mechanisms, such as a putative role for Phosphotase and tensin homolog/phosphoinositide 3-kinase (PTEN/PI3K) as a key point of crosstalk between inflammation and morphogenesis, the protective role of enterocyte sloughing in enteric ischemia-reperfusion and chronic low level inflammation as a driver for colonic metaplasia. These results suggest that SEGMEnT can serve as an integrating platform for the study of inflammation in gastrointestinal disease.

Unraveling intestinal stem cell behavior with models of crypt dynamics

The definition, regulation and function of intestinal stem cells (ISCs) has been hotly debated. Recent discoveries have started to clarify the nature of ISCs, but many questions remain. This review discusses the current advances and controversies of ISC biology as well as theoretical compartmental models that have been coupled with in vivo experimentation to investigate the mechanisms of ISC dynamics during homeostasis, tumorigenesis, repair and development. We conclude our review by discussing the key lingering questions in the field and proposing how many of these questions can be addressed using both compartmental models and experimental techniques.

Fluctuations of cell population in a colonic crypt

The number of stem cells in a colonic crypt is often very small, which leads to large intrinsic fluctuations in the cell population. Based on the model of cell population dynamics with linear feedback in a colonic crypt, we present a stochastic dynamics of the cell population [including stem cells (SCs), transit amplifying cells (TACs), and fully differentiated cells (FDCs)]. The Fano factor, covariance, and susceptibility formulas of the cell population around the steady state are derived by using the Langevin theory. In the range of physiologically reasonable parameter values, it is found that the stationary populations of TACs and FDCs exhibit an approximately threshold behavior as a function of the net growth rate of TACs, and the reproductions of TACs and FDCs can be classified into three regimens: controlled, crossover, and uncontrolled. With the increasing of the net growth rate of TACs, there is a maximum of the relative intrinsic fluctuations (i.e., the Fano factors) of TACs and FDCs in the crossover region. For a fixed differentiation rate and the net growth rate of SCs, the covariance of fluctuations between SCs and TACs has a maximum in the crossover region. However, the susceptibilities of both TACs and FDCs to the net growth rate of TACs have a minimum in the crossover region.

Modeling of stem cell dynamics in human colonic crypts in silico

BACKGROUND

Several possible scenarios of cellular dynamics in human colonic crypts have been inferred from transgenic animal experiments. However, because of the discrepancy in size and physiology between humans and animals, quantitative predictions of tissue renewal and cancer development are difficult to execute.

METHODS

A two-dimensional individual based model was developed for the first time to predict cellular dynamics in human colonic crypts. A simple scenario, in which stem cells were not fixed positionally, divide symmetrically and asymmetrically in a stochastic fashion in the lower part of the crypt, was proposed and implemented in the developed model. Numerical simulations of the model were executed in silico.

RESULTS

By comparing the results of computational simulations with available experimental data, the presented scenario was consistent with various experimental evidence. Using this scenario, we simulated and visualized monoclonal conversion in the human colonic crypt. We also predicted that the propensity for monoclonal expansion of a mutant cell was largely dependent on the phenotype, the cell type, the position and the state of the crypt.

CONCLUSIONS

Using the computational framework developed in this study, model users can verify possible scenarios of stem cell dynamics occurring in human colonic crypts and quantitatively predict cell behavior. Its applicability in scenario verification and predictability makes it a valuable tool for elucidation of stem cell dynamics in human colonic crypts.

Computational models reveal a passive mechanism for cell migration in the crypt

Cell migration in the intestinal crypt is essential for the regular renewal of the epithelium, and the continued upward movement of cells is a key characteristic of healthy crypt dynamics. However, the driving force behind this migration is unknown. Possibilities include mitotic pressure, active movement driven by motility cues, or negative pressure arising from cell loss at the crypt collar. It is possible that a combination of factors together coordinate migration. Here, three different computational models are used to provide insight into the mechanisms that underpin cell movement in the crypt, by examining the consequence of eliminating cell division on cell movement. Computational simulations agree with existing experimental results, confirming that migration can continue in the absence of mitosis. Importantly, however, simulations allow us to infer mechanisms that are sufficient to generate cell movement, which is not possible through experimental observation alone. The results produced by the three models agree and suggest that cell loss due to apoptosis and extrusion at the crypt collar relieves cell compression below, allowing cells to expand and move upwards. This finding suggests that future experiments should focus on the role of apoptosis and cell extrusion in controlling cell migration in the crypt.

A calibrated agent-based computer model of stochastic cell dynamics in normal human colon crypts useful for in silico experiments

BACKGROUND

Normal colon crypts consist of stem cells, proliferating cells, and differentiated cells. Abnormal rates of proliferation and differentiation can initiate colon cancer. We have measured the variation in the number of each of these cell types in multiple crypts in normal human biopsy specimens. This has provided the opportunity to produce a calibrated computational model that simulates cell dynamics in normal human crypts, and by changing model parameter values, to simulate the initiation and treatment of colon cancer.

RESULTS

An agent-based model of stochastic cell dynamics in human colon crypts was developed in the multi-platform open-source application NetLogo. It was assumed that each cell's probability of proliferation and probability of death is determined by its position in two gradients along the crypt axis, a divide gradient and in a die gradient. A cell's type is not intrinsic, but rather is determined by its position in the divide gradient. Cell types are dynamic, plastic, and inter-convertible. Parameter values were determined for the shape of each of the gradients, and for a cell's response to the gradients. This was done by parameter sweeps that indicated the values that reproduced the measured number and variation of each cell type, and produced quasi-stationary stochastic dynamics. The behavior of the model was verified by its ability to reproduce the experimentally observed monocolonal conversion by neutral drift, the formation of adenomas resulting from mutations either at the top or bottom of the crypt, and by the robust ability of crypts to recover from perturbation by cytotoxic agents. One use of the virtual crypt model was demonstrated by evaluating different cancer chemotherapy and radiation scheduling protocols.

CONCLUSIONS

A virtual crypt has been developed that simulates the quasi-stationary stochastic cell dynamics of normal human colon crypts. It is unique in that it has been calibrated with measurements of human biopsy specimens, and it can simulate the variation of cell types in addition to the average number of each cell type. The utility of the model was demonstrated with in silico experiments that evaluated cancer therapy protocols. The model is available for others to conduct additional experiments.

Cell organisation in the colonic crypt: a theoretical comparison of the pedigree and niche concepts

The intestinal mucosa is a monolayer of rapidly self-renewing epithelial cells which is not only responsible for absorption of water and nutrients into the bloodstream but also acts as a protective barrier against harmful microbes entering the body. New functional epithelial cells are produced from stem cells, and their proliferating progeny. These stem cells are found within millions of crypts (tubular pits) spaced along the intestinal tract. The entire intestinal epithelium is replaced every 2–3 days in mice (3–5 days in humans) and hence cell production, differentiation, migration and turnover need to be tightly regulated. Malfunctions in this regulation are strongly linked to inflammatory bowel diseases and to the formation of adenomas and ultimately cancerous tumours. Despite a great deal of biological experimentation and observation, precisely how colonic crypts are regulated to produce mature colonocytes remains unclear. To assist in understanding how cell organisation in crypts is achieved, two very different conceptual models of cell behaviour are developed here, referred to as the ‘pedigree’ and the ‘niche’ models. The pedigree model proposes that crypt cells are largely preprogrammed and receive minimal prompting from the environment as they move through a routine of cell differentiation and proliferation to become mature colonocytes. The niche model proposes that crypt cells are primarily influenced by the local microenvironments along the crypt, and that predetermined cell behaviour plays a negligible role in their development. In this paper we present a computational model of colonic crypts in the mouse, which enables a comparison of the quality and controllability of mature coloncyte production by crypts operating under these two contrasting conceptual models of crypt regulation.

Patterns of proliferative activity in the colonic crypt determine crypt stability and rates of somatic evolution

Epithelial cells in the colon are arranged in cylindrical structures called crypts in which cellular proliferation and migration are tightly regulated. We hypothesized that the proliferation patterns of cells may determine the stability of crypts as well as the rates of somatic evolution towards colorectal tumorigenesis. Here, we propose a linear process model of colonic epithelial cells that explicitly takes into account the proliferation kinetics of cells as a function of cell position within the crypt. Our results indicate that proliferation kinetics has significant influence on the speed of cell movement, kinetics of mutation propagation, and sensitivity of the system to selective effects of mutated cells. We found that, of all proliferation curves tested, those with mitotic activities concentrated near the stem cell, including the actual proliferation kinetics determined in in vivo labeling experiments, have a greater ability of delaying the rate of mutation accumulation in colonic stem cells compared to hypothetical proliferation curves with mitotic activities focused near the top of the crypt column. Our model can be used to investigate the dynamics of proliferation and mutation accumulation in spatially arranged tissues.

Mathematical and computational models have a long and rich history in enhancing our understanding of intestinal epithelial cells. A plethora of models have been proposed to describe different aspects of cellular behavior, including cell proliferation, migration, differentiation, and mutation accumulation. Here, we present a novel approach to examine the effects of proliferation kinetics on the rate of somatic evolution in a spatially arranged model of the colon. Based on our simulation results, we demonstrate that spatially determined proliferation kinetics has the ability to delay the rate of somatic evolution, and changes in proliferation patterns can significantly affect the speed of mutation accumulation. Our work highlights the importance of considering proliferation kinetics as well as the spatial organization of tissues when investigating the dynamics of cancer initiation.

Eph/Ephrin signaling in injury and inflammation

The Eph/ephrin receptor-ligand system plays an important role in embryogenesis and adult life, principally by influencing cell behavior through signaling pathways, resulting in modification of the cell cytoskeleton and cell adhesion. There are 10 EphA receptors, and six EphB receptors, distinguished on sequence difference and binding preferences, that interact with the six glycosylphosphatidylinositol-linked ephrin-A ligands and the three transmembrane ephrin-B ligands, respectively. The Eph/ephrin proteins, originally described as developmental regulators that are expressed at low levels postembryonically, are re-expressed after injury to the optic nerve, spinal cord, and brain in fish, amphibians, rodents, and humans. In rodent spinal cord injury, the up-regulation of EphA4 prevents recovery by inhibiting axons from crossing the injury site. Eph/ephrin proteins may be partly responsible for the phenotypic changes to the vascular endothelium in inflammation, which allows fluid and inflammatory cells to pass from the vascular space into the interstitial tissues. Specifically, EphA2/ephrin-A1 signaling in the lung may be responsible for pulmonary inflammation in acute lung injury. A role in T-cell maturation and chronic inflammation (heart failure, inflammatory bowel disease, and rheumatoid arthritis) is also reported. Although there remains much to learn about Eph/ephrin signaling in human disease, and specifically in injury and inflammation, this area of research raises the exciting prospect that novel therapies will be developed that precisely target these pathways.

A review of spatial computational models for multi-cellular systems, with regard to intestinal crypts and colorectal cancer development

Colon rectal cancers (CRC) are the result of sequences of mutations which lead the intestinal tissue to develop in a carcinoma following a "progression" of observable phenotypes. The actual modeling and simulation of the key biological structures involved in this process is of interest to biologists and physicians and, at the same time, it poses significant challenges from the mathematics and computer science viewpoints. In this report we give an overview of some mathematical models for cell sorting (a basic phenomenon that underlies several dynamical processes in an organism), intestinal crypt dynamics and related problems and open questions. In particular, major attention is devoted to the survey of so-called in-lattice (or grid) models and off-lattice (off-grid) models. The current work is the groundwork for future research on semi-automated hypotheses formation and testing about the behavior of the various actors taking part in the adenoma-carcinoma progression, from regulatory processes to cell-cell signaling pathways.

A review of spatial computational models for multi-cellular systems, with regard to intestinal crypts and colorectal cancer development

Colon rectal cancers (CRC) are the result of sequences of mutations which lead the intestinal tissue to develop in a carcinoma following a "progression" of observable phenotypes. The actual modeling and simulation of the key biological structures involved in this process is of interest to biologists and physicians and, at the same time, it poses significant challenges from the mathematics and computer science viewpoints. In this report we give an overview of some mathematical models for cell sorting (a basic phenomenon that underlies several dynamical processes in an organism), intestinal crypt dynamics and related problems and open questions. In particular, major attention is devoted to the survey of so-called in-lattice (or grid) models and off-lattice (off-grid) models. The current work is the groundwork for future research on semi-automated hypotheses formation and testing about the behavior of the various actors taking part in the adenoma-carcinoma progression, from regulatory processes to cell-cell signaling pathways.

Control of cell adhesion and compartmentalization in the intestinal epithelium

Continuous cell renewal in the intestinal mucosa occurs without disrupting the integrity of the epithelial layer. Despite the restrictions imposed by strong cell-to-cell adhesions, epithelial intestinal cells migrate constantly between tissue compartments. Alterations in cell adhesion and compartmentalization play key roles in diseases of the intestine. In particular, decreased E-cadherin-mediated adhesion during inflammatory bowel disease and loss of EphB/ephrin-B-mediated compartmentalization in colorectal cancer have recently emerged as key players of these prevalent pathologies. Here we will review our current knowledge on how cell-to-cell adhesion, migration and cell positioning are coordinated in the intestinal epithelium. We will highlight what the in vivo genetic analysis of intestinal epithelium has taught us about the complex regulation of cell adhesion and migration in homeostasis and disease.

Epithelial and mesenchymal phenotypic switchings modulate cell motility in metastasis

The most ominous stage of cancer progression is metastasis, or the dissemination of carcinoma cells from the primary site into distant organs. Metastases are often resistant to current extirpative therapies and even the newest biological agents cure only a small subset of patients. Therefore a greater understanding of tumor biology that integrates properties intrinsic to carcinomas with tissue environmental modulators of behavior is needed. In no aspect of tumor progression is this more evident than the acquisition of cell motility that is critical for both escape from the primary tumor and colonization. In this overview, we discuss how this behavior is modified by carcinoma cell phenotypic plasticity that is evidenced by reversible switching between epithelial and mesenchymal phenotypes. The presence or absence of intercellular adhesions mediate these switches and dictate the receptivity towards signals from the extracellular milieu. These signals, which include soluble growth factors, cytokines, and extracellular matrix embedded with matrikines and matricryptines will be discussed in depth. Finally, we will describe a new mode of discerning the balance between epithelioid and mesenchymal movement.