Exploring hypotheses of the actions of TGF-beta 1 in epidermal wound healing using a 3D computational multiscale model of the human epidermis

In vivo and in vitro studies give a paradoxical picture of the actions of the key regulatory factor TGF-beta 1 in epidermal wound healing with it stimulating migration of keratinocytes but also inhibiting their proliferation. To try to reconcile these into an easily visualized 3D model of wound healing amenable for experimentation by cell biologists, a multiscale model of the formation of a 3D skin epithelium was established with TGF-beta 1 literature-derived rule sets and equations embedded within it. At the cellular level, an agent-based bottom-up model that focuses on individual interacting units ( keratinocytes) was used. This was based on literature-derived rules governing keratinocyte behavior and keratinocyte/ECM interactions. The selection of these rule sets is described in detail in this paper. The agent-based model was then linked with a subcellular model of TGF-beta 1 production and its action on keratinocytes simulated with a complex pathway simulator. This multiscale model can be run at a cellular level only or at a combined cellular/subcellular level. It was then initially challenged ( by wounding) to investigate the behavior of keratinocytes in wound healing at the cellular level. To investigate the possible actions of TGF-beta 1, several hypotheses were then explored by deliberately manipulating some of these rule sets at subcellular levels. This exercise readily eliminated some hypotheses and identified a sequence of spatial-temporal actions of TGF-beta 1 for normal successful wound healing in an easy-to-follow 3D model. We suggest this multiscale model offers a valuable, easy-to-visualize aid to our understanding of the actions of this key regulator in wound healing, and provides a model that can now be used to explore pathologies of wound healing

An Adhesion-Dependent Switch between Mechanisms That Determine Motile Cell Shape

Keratocytes are fast-moving cells in which adhesion dynamics are tightly coupled to the actin polymerization motor that drives migration, resulting in highly coordinated cell movement. We have found that modifying the adhesive properties of the underlying substrate has a dramatic effect on keratocyte morphology. Cells crawling at intermediate adhesion strengths resembled stereotypical keratocytes, characterized by a broad, fan-shaped lamellipodium, clearly defined leading and trailing edges, and persistent rates of protrusion and retraction. Cells at low adhesion strength were small and round with highly variable protrusion and retraction rates, and cells at high adhesion strength were large and asymmetrical and, strikingly, exhibited traveling waves of protrusion. To elucidate the mechanisms by which adhesion strength determines cell behavior, we examined the organization of adhesions, myosin II, and the actin network in keratocytes migrating on substrates with different adhesion strengths. On the whole, our results are consistent with a quantitative physical model in which keratocyte shape and migratory behavior emerge from the self-organization of actin, adhesions, and myosin, and quantitative changes in either adhesion strength or myosin contraction can switch keratocytes among qualitatively distinct migration regimes

Model for in vivo progression of tumors based on co-evolving cell population and vasculature

With countless biological details emerging from cancer experiments, there is a growing need for minimal mathematical models which simultaneously advance our understanding of single tumors and metastasis, provide patient-personalized predictions, whilst avoiding excessive hard-to-measure input parameters which complicate simulation, analysis and interpretation. Here we present a model built around a co-evolving resource network and cell population, yielding good agreement with primary tumors in a murine mammary cell line EMT6-HER2 model in BALB/c mice and with clinical metastasis data. Seeding data about the tumor and its vasculature from in vivo images, our model predicts corridors of future tumor growth behavior and intervention response. A scaling relation enables the estimation of a tumor's most likely evolution and pinpoints specific target sites to control growth. Our findings suggest that the clinically separate phenomena of individual tumor growth and metastasis can be viewed as mathematical copies of each other differentiated only by network structure

A Cellular Potts Model simulating cell migration on and in matrix environments

Cell migration on and through extracellular matrix plays a critical role in a wide variety of physiological and pathological phenomena, and in scaffold-based tissue engineering. Migration is regulated by a number of extracellular matrix- or cell-derived biophysical parameters, such as matrix fiber orientation, gap size, and elasticity, or cell deformation, proteolysis, and adhesion. We here present an extended Cellular Potts Model (CPM) able to qualitatively and quantitatively describe cell migratory phenotype on both two-dimensional substrates and within three-dimensional environments, in a close comparison with experimental evidence. As distinct features of our approach, the cells are represented by compartmentalized discrete objects, differentiated in the nucleus and in the cytosolic region, while the extracellular matrix is composed of a fibrous mesh and of a homogeneous fluid. Our model provides a strong correlation of the directionality of migration with the topological ECM distribution and, further, a biphasic dependence of migration on the matrix density, and in part adhesion, in both two-dimensional and three-dimensional settings. Moreover, we demonstrate that the directional component of cell movement is strongly correlated with the topological distribution of the ECM fibrous network. In the three-dimensional networks, we also investigate the effects of the matrix mechanical microstructure, observing that, at a given distribution of fibers, cell motility has a subtle bimodal relation with the elasticity of the scaffold. Finally, cell locomotion requires deformation of the cell's nucleus and/or cell-derived proteolysis of steric fibrillar obstacles within rather rigid matrices characterized by small pores, not, however, for sufficiently large pores. In conclusion, we here propose a mathematical modeling approach that serves to characterize cell migration as a biological phenomen in health, disease and tissue engineering applications. The research that led to the present paper was partially supported by a grant of the group GNFM of INdA

A Model for the Interplay of Receptor Recycling and Receptor-Mediated Contact in T Cells

Orientation of organelles inside T cells (TC) toward antigen-presenting cells (APC) ensures that the immune response is properly directed, but the orientation mechanisms remain largely unknown. Structural dynamics of TC are coupled to dynamics of T-cell receptor (TCR), which recognizes antigen on the APC surface. Engagement of the TCR triggers its internalization followed by delayed polarized recycling to the plasma membrane through the submembrane recycling compartment (RC), which organelle shares intracellular location with the TC effector apparatus. TCR engagement also triggers TC-APC interface expansion enabling further receptor engagement. To analyze the interplay of the cell-cell contact and receptor dynamics, we constructed a new numerical model. The new model displays the experimentally observed selective stabilization of the contact initiated next to the RC, and only transient formation of contact diametrically opposed to the RC. In the general case wherein the TC-APC contact is initiated in an arbitrary orientation to the RC, the modeling predicts that the contact dynamics and receptor recycling can interact, resulting effectively in migration of the contact to the TC surface domain adjacent to the submembrane RC. Using three-dimensional live-cell confocal microscopy, we obtain data consistent with this unexpected behavior. We conclude that a TC can stabilize its contact with an APC by aligning it with the polarized intracellular traffic of TCR. The results also suggest that the orientation of TC organelles, such as the RC and the effector apparatus, toward the APC can be achieved without any intracellular translocation of the organelles

A Quorum-Sensing Factor in Vegetative Dictyostelium Discoideum Cells Revealed by Quantitative Migration Analysis

Background: Many cells communicate through the production of diffusible signaling molecules that accumulate and once a critical concentration has been reached, can activate or repress a number of target genes in a process termed quorum sensing (QS). In the social amoeba Dictyostelium discoideum, QS plays an important role during development. However little is known about its effect on cell migration especially in the growth phase. Methods and Findings: To investigate the role of cell density on cell migration in the growth phase, we use multisite timelapse microscopy and automated cell tracking. This analysis reveals a high heterogeneity within a given cell population, and the necessity to use large data sets to draw reliable conclusions on cell motion. In average, motion is persistent for short periods of time (tƒ5min), but normal diffusive behavior is recovered over longer time periods. The persistence times are positively correlated with the migrated distances. Interestingly, the migrated distance decreases as well with cell density. The adaptation of cell migration to cell density highlights the role of a secreted quorum sensing factor (QSF) on cell migration. Using a simple model describing the balance between the rate of QSF generation and the rate of QSF dilution, we were able to gather all experimental results into a single master curve, showing a sharp cell transition between high and low motile behaviors with increasing QSF. Conclusion: This study unambiguously demonstrates the central role played by QSF on amoeboid motion in the growt

Interstitial cell migration: integrin-dependent and alternative adhesion mechanisms

Adhesion and migration are integrated cell functions that build, maintain and remodel the multicellular organism. In migrating cells, integrins are the main transmembrane receptors that provide dynamic interactions between extracellular ligands and actin cytoskeleton and signalling machineries. In parallel to integrins, other adhesion systems mediate adhesion and cytoskeletal coupling to the extracellular matrix (ECM). These include multifunctional cell surface receptors (syndecans and CD44) and discoidin domain receptors, which together coordinate ligand binding with direct or indirect cytoskeletal coupling and intracellular signalling. We review the way that the different adhesion systems for ECM components impact cell migration in two- and three-dimensional migration models. We further discuss the hierarchy of these concurrent adhesion systems, their specific tasks in cell migration and their contribution to migration in three-dimensional multi-ligand tissue environments

Acidic Extracellular pH Promotes Activation of Integrin αvβ3

Acidic extracellular pH is characteristic of the cell microenvironment in several important physiological and pathological contexts. Although it is well established that acidic extracellular pH can have profound effects on processes such as cell adhesion and migration, the underlying molecular mechanisms are largely unknown. Integrin receptors physically connect cells to the extracellular matrix, and are thus likely to modulate cell responses to extracellular conditions. Here, we examine the role of acidic extracellular pH in regulating activation of integrin [alpha]v[beta]3. Through computational molecular dynamics simulations, we find that acidic extracellular pH promotes opening of the [alpha]v[beta]3 headpiece, indicating that acidic pH can thereby facilitate integrin activation. This prediction is consistent with our flow cytometry and atomic force microscope-mediated force spectroscopy assays of integrin [alpha]v[beta]3 on live cells, which both demonstrate that acidic pH promotes activation at the intact cell surface. Finally, quantification of cell morphology and migration measurements shows that acidic extracellular pH affects cell behavior in a manner that is consistent with increased integrin activation. Taken together, these computational and experimental results suggest a new and complementary mechanism of integrin activation regulation, with associated implications for cell adhesion and migration in regions of altered pH that are relevant to wound healing and cancer.National Institute of Biomedical Imaging and Bioengineering (U.S.) (Award Number T32EB006348)Massachusetts Institute of Technology (Collamore-Rogers Fellowship)National Institutes of Health (U.S.) (NIH Cell Migration Consortium Grant U54-GM069668)National Science Foundation (U.S.) (CAREER Award)Singapore-MIT Alliance for Research and Technology (BioSystem and Micromechanics (BioSyM) Interdisciplinary Research Group

Development of a Three Dimensional Multiscale Computational Model of the Human Epidermis

Transforming Growth Factor (TGF-β1) is a member of the TGF-beta superfamily ligand-receptor network. and plays a crucial role in tissue regeneration. The extensive in vitro and in vivo experimental literature describing its actions nevertheless describe an apparent paradox in that during re-epithelialisation it acts as proliferation inhibitor for keratinocytes. The majority of biological models focus on certain aspects of TGF-β1 behaviour and no one model provides a comprehensive story of this regulatory factor's action. Accordingly our aim was to develop a computational model to act as a complementary approach to improve our understanding of TGF-β1. In our previous study, an agent-based model of keratinocyte colony formation in 2D culture was developed. In this study this model was extensively developed into a three dimensional multiscale model of the human epidermis which is comprised of three interacting and integrated layers: (1) an agent-based model which captures the biological rules governing the cells in the human epidermis at the cellular level and includes the rules for injury induced emergent behaviours, (2) a COmplex PAthway SImulator (COPASI) model which simulates the expression and signalling of TGF-β1 at the sub-cellular level and (3) a mechanical layer embodied by a numerical physical solver responsible for resolving the forces exerted between cells at the multi-cellular level. The integrated model was initially validated by using it to grow a piece of virtual epidermis in 3D and comparing the in virtuo simulations of keratinocyte behaviour and of TGF-β1 signalling with the extensive research literature describing this key regulatory protein. This research reinforces the idea that computational modelling can be an effective additional tool to aid our understanding of complex systems. In the accompanying paper the model is used to explore hypotheses of the functions of TGF-β1 at the cellular and subcellular level on different keratinocyte populations during epidermal wound healing

An agent-based model of the response to angioplasty and bare-metal stent deployment in an atherosclerotic blood vessel

Purpose: While animal models are widely used to investigate the development of restenosis in blood vessels following an intervention, computational models offer another means for investigating this phenomenon. A computational model of the response of a treated vessel would allow investigators to assess the effects of altering certain vessel- and stent-related variables. The authors aimed to develop a novel computational model of restenosis development following an angioplasty and bare-metal stent implantation in an atherosclerotic vessel using agent-based modeling techniques. The presented model is intended to demonstrate the body's response to the intervention and to explore how different vessel geometries or stent arrangements may affect restenosis development. Methods: The model was created on a two-dimensional grid space. It utilizes the post-procedural vessel lumen diameter and stent information as its input parameters. The simulation starting point of the model is an atherosclerotic vessel after an angioplasty and stent implantation procedure. The model subsequently generates the final lumen diameter, percent change in lumen cross-sectional area, time to lumen diameter stabilization, and local concentrations of inflammatory cytokines upon simulation completion. Simulation results were directly compared with the results from serial imaging studies and cytokine levels studies in atherosclerotic patients from the relevant literature. Results: The final lumen diameter results were all within one standard deviation of the mean lumen diameters reported in the comparison studies. The overlapping-stent simulations yielded results that matched published trends. The cytokine levels remained within the range of physiological levels throughout the simulations. Conclusion: We developed a novel computational model that successfully simulated the development of restenosis in a blood vessel following an angioplasty and bare-metal stent deployment based on the characteristics of the vessel crosssection and stent. A further development of this model could ultimately be used as a predictive tool to depict patient outcomes and inform treatment options. © 2014 Curtin, Zhou