Roeland M. H. Merks

Contact-inhibited chemotaxis in de novo and sprouting blood-vessel growth.

Blood vessels form either when dispersed endothelial cells (the cells lining the inner walls of fully formed blood vessels) organize into a vessel network (vasculogenesis), or by sprouting or splitting of existing blood vessels (angiogenesis). Although they are closely related biologically, no current model explains both phenomena with a single biophysical mechanism. Most computational models describe sprouting at the level of the blood vessel, ignoring how cell behavior drives branch splitting during sprouting. We present a cell-based, Glazier–Graner–Hogeweg model (also called Cellular Potts Model) simulation of the initial patterning before the vascular cords form lumens, based on plausible behaviors of endothelial cells. The endothelial cells secrete a chemoattractant, which attracts other endothelial cells. As in the classic Keller–Segel model, chemotaxis by itself causes cells to aggregate into isolated clusters. However, including experimentally observed VE-cadherin–mediated contact inhibition of chemotaxis in the simulation causes randomly distributed cells to organize into networks and cell aggregates to sprout, reproducing aspects of both de novo and sprouting blood-vessel growth. We discuss two branching instabilities responsible for our results. Cells at the surfaces of cell clusters attempting to migrate to the centers of the clusters produce a buckling instability. In a model variant that eliminates the surface–normal force, a dissipative mechanism drives sprouting, with the secreted chemical acting both as a chemoattractant and as an inhibitor of pseudopod extension. Both mechanisms would also apply if force transmission through the extracellular matrix rather than chemical signaling mediated cell–cell interactions. The branching instabilities responsible for our results, which result from contact inhibition of chemotaxis, are both generic developmental mechanisms and interesting examples of unusual patterning instabilities.

Emergence of tissue polarization from synergy of intracellular and extracellular auxin signaling

Plant development is exceptionally flexible as manifested by its potential for organogenesis and regeneration, which are processes involving rearrangements of tissue polarities. Fundamental questions concern how individual cells can polarize in a coordinated manner to integrate into the multicellular context. In canalization models, the signaling molecule auxin acts as a polarizing cue, and feedback on the intercellular auxin flow is key for synchronized polarity rearrangements. We provide a novel mechanistic framework for canalization, based on up‐to‐date experimental data and minimal, biologically plausible assumptions. Our model combines the intracellular auxin signaling for expression of PINFORMED (PIN) auxin transporters and the theoretical postulation of extracellular auxin signaling for modulation of PIN subcellular dynamics. Computer simulations faithfully and robustly recapitulated the experimentally observed patterns of tissue polarity and asymmetric auxin distribution during formation and regeneration of vascular systems and during the competitive regulation of shoot branching by apical dominance. Additionally, our model generated new predictions that could be experimentally validated, highlighting a mechanistically conceivable explanation for the PIN polarization and canalization of the auxin flow in plants.

VirtualLeaf: An Open-Source Framework for Cell-Based Modeling of Plant Tissue Growth and Development

Plant organs, including leaves and roots, develop by means of a multilevel cross talk between gene regulation, patterned cell division and cell expansion, and tissue mechanics. The multilevel regulatory mechanisms complicate classic molecular genetics or functional genomics approaches to biological development, because these methodologies implicitly assume a direct relation between genes and traits at the level of the whole plant or organ. Instead, understanding gene function requires insight into the roles of gene products in regulatory networks, the conditions of gene expression, etc. This interplay is impossible to understand intuitively. Mathematical and computer modeling allows researchers to design new hypotheses and produce experimentally testable insights. However, the required mathematics and programming experience makes modeling poorly accessible to experimental biologists. Problem-solving environments provide biologically intuitive in silico objects (“cells”, “regulation networks”) required for setting up a simulation and present those to the user in terms of familiar, biological terminology. Here, we introduce the cell-based computer modeling framework VirtualLeaf for plant tissue morphogenesis. The current version defines a set of biologically intuitive C++ objects, including cells, cell walls, and diffusing and reacting chemicals, that provide useful abstractions for building biological simulations of developmental processes. We present a step-by-step introduction to building models with VirtualLeaf, providing basic example models of leaf venation and meristem development. VirtualLeaf-based models provide a means for plant researchers to analyze the function of developmental genes in the context of the biophysics of growth and patterning. VirtualLeaf is an ongoing open-source software project (http://virtualleaf.googlecode.com) that runs on Windows, Mac, and Linux.

Dynamic mechanisms of blood vessel growth

The formation of a polygonal configuration of proto-blood-vessels from initially dispersed cells is the first step in the development of the circulatory system in vertebrates. This initial vascular network later expands to form new blood vessels, primarily via a sprouting mechanism. We review a range of recent results obtained with a Monte Carlo model of chemotactically migrating cells which can explain both de novo blood vessel growth and aspects of blood vessel sprouting. We propose that the initial network forms via a percolation-like instability depending on cell shape, or through an alternative contact-inhibition of motility mechanism which also reproduces aspects of sprouting blood vessel growth.

Mechanical cell-matrix feedback explains pairwise and collective endothelial cell behavior in vitro.

In vitro cultures of endothelial cells are a widely used model system of the collective behavior of endothelial cells during vasculogenesis and angiogenesis. When seeded in an extracellular matrix, endothelial cells can form blood vessel-like structures, including vascular networks and sprouts. Endothelial morphogenesis depends on a large number of chemical and mechanical factors, including the compliancy of the extracellular matrix, the available growth factors, the adhesion of cells to the extracellular matrix, cell-cell signaling, etc. Although various computational models have been proposed to explain the role of each of these biochemical and biomechanical effects, the understanding of the mechanisms underlying in vitro angiogenesis is still incomplete. Most explanations focus on predicting the whole vascular network or sprout from the underlying cell behavior, and do not check if the same model also correctly captures the intermediate scale: the pairwise cell-cell interactions or single cell responses to ECM mechanics. Here we show, using a hybrid cellular Potts and finite element computational model, that a single set of biologically plausible rules describing (a) the contractile forces that endothelial cells exert on the ECM, (b) the resulting strains in the extracellular matrix, and (c) the cellular response to the strains, suffices for reproducing the behavior of individual endothelial cells and the interactions of endothelial cell pairs in compliant matrices. With the same set of rules, the model also reproduces network formation from scattered cells, and sprouting from endothelial spheroids. Combining the present mechanical model with aspects of previously proposed mechanical and chemical models may lead to a more complete understanding of in vitro angiogenesis.

Morphogenesis of the branching reef coral Madracis mirabilis

Understanding external deciding factors in growth and morphology of reef corals is essential to elucidate the role of corals in marine ecosystems, and to explain their susceptibility to pollution and global climate change. Here, we extend on a previously presented model for simulating the growth and form of a branching coral and we compare the simulated morphologies to three–dimensional (3D) images of the coral species Madracis mirabilis. Simulation experiments and isotope analyses of M. mirabilis skeletons indicate that external gradients of dissolved inorganic carbon (DIC) determine the morphogenesis of branching, phototrophic corals. In the simulations we use a first principle model of accretive growth based on local interactions between the polyps. The only species–specific information in the model is the average size of a polyp. From flow tank and simulation studies it is known that a relatively large stagnant and diffusion dominated region develops within a branching colony. We have used this information by assuming in our model that growth is entirely driven by a diffusion–limited process, where DIC supply represents the limiting factor. With such model constraints it is possible to generate morphologies that are virtually indistinguishable from the 3D images of the actual colonies.

SHORT-ROOT and SCARECROW regulate leaf growth in Arabidopsis by stimulating S-phase progression of the cell cycle

SHORT-ROOT (SHR) and SCARECROW (SCR) are required for stem cell maintenance in the Arabidopsis (Arabidopsis thaliana) root meristem, ensuring its indeterminate growth. Mutation of SHR and SCR genes results in disorganization of the quiescent center and loss of stem cell activity, resulting in the cessation of root growth. This paper reports on the role of SHR and SCR in the development of leaves, which, in contrast to the root, have a determinate growth pattern and lack a persistent stem cell niche. Our results demonstrate that inhibition of leaf growth in shr and scr mutants is not a secondary effect of the compromised root development but is caused by an effect on cell division in the leaves: a reduced cell division rate and early exit of the proliferation phase. Consistent with the observed cell division phenotype, the expression of SHR and SCR genes in leaves is closely associated with cell division activity in most cell types. The increased cell cycle duration is due to a prolonged S-phase duration, which is mediated by up-regulation of cell cycle inhibitors known to restrain the activity of the transcription factor, E2Fa. Therefore, we conclude that, in contrast to their specific roles in cortex/endodermis differentiation and stem cell maintenance in the root, SHR and SCR primarily function as general regulators of cell proliferation in leaves.

Cellular Potts Modeling of Tumor Growth, Tumor Invasion, and Tumor Evolution

Despite a growing wealth of available molecular data, the growth of tumors, invasion of tumors into healthy tissue, and response of tumors to therapies are still poorly understood. Although genetic mutations are in general the first step in the development of a cancer, for the mutated cell to persist in a tissue, it must compete against the other, healthy or diseased cells, for example by becoming more motile, adhesive, or multiplying faster. Thus, the cellular phenotype determines the success of a cancer cell in competition with its neighbors, irrespective of the genetic mutations or physiological alterations that gave rise to the altered phenotype. What phenotypes can make a cell “successful” in an environment of healthy and cancerous cells, and how? A widely used tool for getting more insight into that question is cell-based modeling. Cell-based models constitute a class of computational, agent-based models that mimic biophysical and molecular interactions between cells. One of the most widely used cell-based modeling formalisms is the cellular Potts model (CPM), a lattice-based, multi particle cell-based modeling approach. The CPM has become a popular and accessible method for modeling mechanisms of multicellular processes including cell sorting, gastrulation, or angiogenesis. The CPM accounts for biophysical cellular properties, including cell proliferation, cell motility, and cell adhesion, which play a key role in cancer. Multiscale models are constructed by extending the agents with intracellular processes including metabolism, growth, and signaling. Here we review the use of the CPM for modeling tumor growth, tumor invasion, and tumor progression. We argue that the accessibility and flexibility of the CPM, and its accurate, yet coarse-grained and computationally efficient representation of cell and tissue biophysics, make the CPM the method of choice for modeling cellular processes in tumor development.

Models of coral growth: spontaneous branching, compactification and the Laplacian growth assumption.

In stony corals it is often observed that specimens collected from a sheltered growth site have more open and more thinly branched growth forms than specimens of the same species from more exposed growth sites, where stronger water currents are found. This observation was explained using an abiotic computational model inspired by coral growth, in which the growth velocity depended locally on the absorption of a resource dispersed by advection and diffusion (Kaandorp and Sloot, J. Theor. Biol 209 (2001) 257). In that model a morphological range was found; as the Péclet-number (indicating the relative importance of advective and diffusive nutrient transport) was increased, more compact and spherical growth forms were found. Two unsatisfactory items have remained in this model, which we address in the present paper. First, an explicit curvature rule was responsible for branching. In this work we show that the curvature rule is not needed: the model exhibits spontaneous branching, provided that the resource field is computed with enough precision. Second, previously no explanation was given for the morphological range found in the simulations. Here we show that such an explanation is given by the conditions under which spontaneous branching occurs in our model, in which the compactness of the growth forms depends on the ratio of the rates of growth and nutrient transport. We did not find an effect of flow. This suggests that the computational evidence that hydrodynamics influences the compactness of corals in laminar flows may not be conclusive. The applicability of the Laplacian growth paradigm to understand coral growth is discussed.

From Genes to Organisms Via the Cell: A Problem-Solving Environment for Multicellular Development

To gain performance, developers often build scientific applications in procedural languages, such as C or Fortran, which unfortunately reduces flexibility. To address this imbalance, the authors present CompuCell3D, a multitiered, flexible, and scalable problem-solving environment for morphogenesis simulations thats written in C++ using object-oriented design patterns.