Julio M. Belmonte

Multi-scale modeling of tissues using CompuCell3D.

The study of how cells interact to produce tissue development, homeostasis, or diseases was, until recently, almost purely experimental. Now, multi-cell computer simulation methods, ranging from relatively simple cellular automata to complex immersed-boundary and finite-element mechanistic models, allow in silico study of multi-cell phenomena at the tissue scale based on biologically observed cell behaviors and interactions such as movement, adhesion, growth, death, mitosis, secretion of chemicals, chemotaxis, etc. This tutorial introduces the lattice-based Glazier-Graner-Hogeweg (GGH) Monte Carlo multi-cell modeling and the open-source GGH-based CompuCell3D simulation environment that allows rapid and intuitive modeling and simulation of cellular and multi-cellular behaviors in the context of tissue formation and subsequent dynamics. We also present a walkthrough of four biological models and their associated simulations that demonstrate the capabilities of the GGH and CompuCell3D.

A Multi-cell, Multi-scale Model of Vertebrate Segmentation and Somite Formation

Somitogenesis, the formation of the bodys primary segmental structure common to all vertebrate development, requires coordination between biological mechanisms at several scales. Explaining how these mechanisms interact across scales and how events are coordinated in space and time is necessary for a complete understanding of somitogenesis and its evolutionary flexibility. So far, mechanisms of somitogenesis have been studied independently. To test the consistency, integrability and combined explanatory power of current prevailing hypotheses, we built an integrated clock-and-wavefront model including submodels of the intracellular segmentation clock, intercellular segmentation-clock coupling via Delta/Notch signaling, an FGF8 determination front, delayed differentiation, clock-wavefront readout, and differential-cell-cell-adhesion-driven cell sorting. We identify inconsistencies between existing submodels and gaps in the current understanding of somitogenesis mechanisms, and propose novel submodels and extensions of existing submodels where necessary. For reasonable initial conditions, 2D simulations of our model robustly generate spatially and temporally regular somites, realistic dynamic morphologies and spontaneous emergence of anterior-traveling stripes of Lfng. We show that these traveling stripes are pseudo-waves rather than true propagating waves. Our model is flexible enough to generate interspecies-like variation in somite size in response to changes in the PSM growth rate and segmentation-clock period, and in the number and width of Lfng stripes in response to changes in the PSM growth rate, segmentation-clock period and PSM length.

Somites Without a Clock

The formation of body segments in vertebrate embryos involves local cell interactions independent of cyclic gene expression. [Also see Perspective by Kondo] The formation of body segments (somites) in vertebrate embryos is accompanied by molecular oscillations (segmentation clock). Interaction of this oscillator with a wave traveling along the body axis (the clock-and-wavefront model) is generally believed to control somite number, size, and axial identity. Here we show that a clock-and-wavefront mechanism is unnecessary for somite formation. Non-somite mesoderm treated with Noggin generates many somites that form simultaneously, without cyclic expression of Notch-pathway genes, yet have normal size, shape, and fate. These somites have axial identity: The Hox code is fixed independently of somite fate. However, these somites are not subdivided into rostral and caudal halves, which is necessary for neural segmentation. We propose that somites are self-organizing structures whose size and shape is controlled by local cell-cell interactions. Cell-Cell Interactions in Development In vertebrate embryos, the number, size, and positional identity of mesodermal segments (somites) located bilaterally along the anteriorposterior axis is widely believed to be controlled by a molecular clock of oscillating gene expression interacting with a traveling wave of signals to determine how many cells make up a somite. Dias et al. (p. 791, published online 9 January; see the Perspective by Kondo) reveal that it is possible to generate somites of normal size, shape, and identity without either a clock or a wavefront. Instead, the findings suggest that somite size and shape are regulated by local cell-cell interactions.

A theory that predicts behaviors of disordered cytoskeletal networks

Morphogenesis in animal tissues is largely driven by actomyosin networks, through tensions generated by an active contractile process. Although the network components and their properties are known, and networks can be reconstituted in vitro, the requirements for contractility are still poorly understood. Here, we describe a theory that predicts whether an isotropic network will contract, expand, or conserve its dimensions. This analytical theory correctly predicts the behavior of simulated networks, consisting of filaments with varying combinations of connectors, and reveals conditions under which networks of rigid filaments are either contractile or expansile. Our results suggest that pulsatility is an intrinsic behavior of contractile networks if the filaments are not stable but turn over. The theory offers a unifying framework to think about mechanisms of contractions or expansion. It provides the foundation for studying a broad range of processes involving cytoskeletal networks and a basis for designing synthetic networks.

A Liver-Centric Multiscale Modeling Framework for Xenobiotics.

We describe a multi-scale, liver-centric in silico modeling framework for acetaminophen pharmacology and metabolism. We focus on a computational model to characterize whole body uptake and clearance, liver transport and phase I and phase II metabolism. We do this by incorporating sub-models that span three scales; Physiologically Based Pharmacokinetic (PBPK) modeling of acetaminophen uptake and distribution at the whole body level, cell and blood flow modeling at the tissue/organ level and metabolism at the sub-cellular level. We have used standard modeling modalities at each of the three scales. In particular, we have used the Systems Biology Markup Language (SBML) to create both the whole-body and sub-cellular scales. Our modeling approach allows us to run the individual sub-models separately and allows us to easily exchange models at a particular scale without the need to extensively rework the sub-models at other scales. In addition, the use of SBML greatly facilitates the inclusion of biological annotations directly in the model code. The model was calibrated using human in vivo data for acetaminophen and its sulfate and glucuronate metabolites. We then carried out extensive parameter sensitivity studies including the pairwise interaction of parameters. We also simulated population variation of exposure and sensitivity to acetaminophen. Our modeling framework can be extended to the prediction of liver toxicity following acetaminophen overdose, or used as a general purpose pharmacokinetic model for xenobiotics.

A disassembly-driven mechanism explains F-actin-mediated chromosome transport in starfish oocytes

While contraction of sarcomeric actomyosin assemblies is well understood, this is not the case for disordered networks of actin filaments (F-actin) driving diverse essential processes in animal cells. For example, at the onset of meiosis in starfish oocytes a contractile F-actin network forms in the nuclear region transporting embedded chromosomes to the assembling microtubule spindle. Here, we addressed the mechanism driving contraction of this 3D disordered F-actin network by comparing quantitative observations to computational models. We analyzed 3D chromosome trajectories and imaged filament dynamics to monitor network behavior under various physical and chemical perturbations. We found no evidence of myosin activity driving network contractility. Instead, our observations are well explained by models based on a disassembly-driven contractile mechanism. We reconstitute this disassembly-based contractile system in silico revealing a simple architecture that robustly drives chromosome transport to prevent aneuploidy in the large oocyte, a prerequisite for normal embryonic development.

A Notch positive feedback in the intestinal stem cell niche is essential for stem cell self‐renewal

The intestinal epithelium is the fastest regenerative tissue in the body, fueled by fast‐cycling stem cells. The number and identity of these dividing and migrating stem cells are maintained by a mosaic pattern at the base of the crypt. How the underlying regulatory scheme manages this dynamic stem cell niche is not entirely clear. We stimulated intestinal organoids with Notch ligands and inhibitors and discovered that intestinal stem cells employ a positive feedback mechanism via direct Notch binding to the second intron of the Notch1 gene. Inactivation of the positive feedback by CRISPR/Cas9 mutation of the binding sequence alters the mosaic stem cell niche pattern and hinders regeneration in organoids. Dynamical system analysis and agent‐based multiscale stochastic modeling suggest that the positive feedback enhances the robustness of Notch‐mediated niche patterning. This study highlights the importance of feedback mechanisms in spatiotemporal control of the stem cell niche.

Large-scale microtubule networks contract quite well

The quantitative investigation of how networks of microtubules contract can boost our understanding of actin biology.

Filopodial-Tension Model of Convergent-Extension of Tissues

In convergent-extension (CE), a planar-polarized epithelial tissue elongates (extends) in-plane in one direction while shortening (converging) in the perpendicular in-plane direction, with the cells both elongating and intercalating along the converging axis. CE occurs during the development of most multicellular organisms. Current CE models assume cell or tissue asymmetry, but neglect the preferential filopodial activity along the convergent axis observed in many tissues. We propose a cell-based CE model based on asymmetric filopodial tension forces between cells and investigate how cell-level filopodial interactions drive tissue-level CE. The final tissue geometry depends on the balance between external rounding forces and cell-intercalation traction. Filopodial-tension CE is robust to relatively high levels of planar cell polarity misalignment and to the presence of non-active cells. Addition of a simple mechanical feedback between cells fully rescues and even improves CE of tissues with high levels of polarity misalignments. Our model extends easily to three dimensions, with either one converging and two extending axes, or two converging and one extending axes, producing distinct tissue morphologies, as observed in vivo.

Fibroblast state switching orchestrates dermal maturation and wound healing

Murine dermis contains functionally and spatially distinct fibroblast lineages that cease to proliferate in early postnatal life. Here, we propose a model in which a negative feedback loop between extracellular matrix (ECM) deposition and fibroblast proliferation determines dermal architecture. Virtual‐tissue simulations of our model faithfully recapitulate dermal maturation, predicting a loss of spatial segregation of fibroblast lineages and dictating that fibroblast migration is only required for wound healing. To test this, we performed in vivo live imaging of dermal fibroblasts, which revealed that homeostatic tissue architecture is achieved without active cell migration. In contrast, both fibroblast proliferation and migration are key determinants of tissue repair following wounding. The results show that tissue‐scale coordination is driven by the interdependence of cell proliferation and ECM deposition, paving the way for identifying new therapeutic strategies to enhance skin regeneration.