Sharon R. Lubkin

Multiphase Mechanics of Capsule Formation in Tumors

The presence of a capsule around a tumor is known to be correlated with benign status, and the absence of a capsule often has negative implications for patient prognosis. A mechanical description is presented of the growth of a tumor and the resulting deformations of surrounding normal tissue. A mathematical model of the mechanics is analyzed using physical parameters measured in vivo and in vitro. The model has only three dimensionless parameters, and its results are very robust with respect to parameter variation. We show that the presence of contractility in the surrounding tissue, corresponding to a host defense, can make an existing capsule denser and constrain the tumor better, but cannot be responsible for the observed pressure gradients in encapsulated tumors. Some implications for treatment are discussed.

A mechanism for early branching in lung morphogenesis

The lung is a highly branched fluid-filled structure, that develops by repeated dichotomous branching of a single bud off the foregut, of epithelium invaginating into mesenchyme. Incorporating the known stress response of developing lung tissues, we model the developing embryonic lung in fluid mechanical terms. We suggest that the repeated branching of the early embryonic lung can be understood as the natural physical consequence of the interactions of two or more plastic substances with surface tension between them. The model makes qualitative and quantitative predictions, as well as suggesting an explanation for such observed phenomena as the asymmetric second branching of the embryonic bronchi.

Spatial pattern formation in biology: I. Dermal wound healing. II. Bacterial patterns

Abstract Although the development of spatial pattern and form is a central issue in biology the mechanisms which generate them are generally unknown. The interdisciplinary modelling challenge is to construct realistic mechanisms which capture the key biological processes and show how they are orchestrated to create the observed pattern. We discuss two specific patterning problems of current widespread interest in biomedicine. In the first, possible mechanisms of dermal wound healing are reviewed with a discussion of what is needed of realistic models for studying wound healing. We then list a series of open problems. In the second problem we describe a model for the complex patterns formed by bacterial colonies, specifically Escherichia coli, and derive and analyse a model firmly based on experimental data. The results from the model compare well with experiment. Mathematically, the class of models discussed gives rise to novel systems of partial differential equations which pose challenging problems, both analytical and numerical. The models have provided the experimentalist with insight as to how such patterns might be formed and have suggested possible experiments to elucidate the underlying biological processes.

Branched organs: mechanics of morphogenesis by multiple mechanisms.

Branching morphogenesis is ubiquitous and important in creating bulk transport systems. Branched ducts can be generated by several different mechanisms including growth, cell rearrangements, contractility, adhesion changes, and other mechanisms. We have developed several models of the mechanics of cleft formation, which we review. We discuss the implications of several candidate mechanisms and review what has been found in models and in experiments.

Morphogenetic implications of peristaltic fluid-tissue dynamics in the embryonic lung.

Peristalsis begins in the lung as soon as the smooth muscle forms, and persists until birth. Since the prenatal lung is liquid-filled, smooth muscle action can deform tissues and transport fluid far from the immediately adjacent tissues. Stretching of embryonic tissues and sensation of internal fluid flows have been shown to have potent morphogenetic effects. We hypothesize that these effects are at work in lung morphogenesis. To place that hypothesis in a quantitative framework, we analyze a model of the fluid-structure interactions between embryonic tissues and lumen fluid resulting from peristaltic waves that partially occlude the airway. We find that if the airway is closed, deformations are synchronized; by contrast, if the trachea is open, maximal occlusion precedes maximal pressure. We perform a parametric analysis of how occlusion, stretch, and flow depend on tissue stiffnesses, smooth muscle force, tissue shape and size, and fluid viscosity. We find that most of these relationships are governed by simple ratios.

Quantifying stretch and secretion in the embryonic lung: Implications for morphogenesis.

Branching in the embryonic lung is controlled by a variety of morphogens. Mechanics is also believed to play a significant role in lung branching. The relative roles and interactions of these two broad factors are challenging to determine. We considered three hypotheses for explaining why tracheal occlusion triples branching with no overall increase in size. Both hypotheses are based on tracheal occlusion blocking the exit of secretions. (H1) Increased lumen pressure stretches tissues; stretch receptors at shoulders of growing tips increase local rate of branching. (H2) Blocking exit of secretions blocks advective transport of morphogens, leading to (H2a) increased overall concentration of morphogens or (H2b) increased flux of morphogens at specific locations. We constructed and analyzed computational models of tissue stretch and solute transport in a 3D lung geometry. Observed tissue stresses and stretches were predominantly in locations unrelated to subsequent branch locations, suggesting that tissue stretch (H1) is not the mechanism of enhancement of branching. Morphogen concentration in the mesenchyme (H2a) increased with tracheal occlusion, consistent with previously reported results. Morphogen flux at the epithelial surface (H2b) completely changed its distribution pattern when the trachea was occluded, tripling the number of locations at which it was elevated. Our results are consistent with the hypothesis that tracheal occlusion blocks outflow of secretions, leading to a higher number of high-flux locations at branching tips, in turn leading to a large increase in number of branching locations.

Quantifying cellular and subcellular stretches in embryonic lung epithelia under peristalsis: where to look for mechanosensing

Peristalsis begins in the lung as soon as the smooth muscle (SM) forms, and persists until birth. As the prenatal lung is filled with liquid, SM action can, through lumen pressure, deform tissues far from the immediately adjacent tissues. Stretching of embryonic tissues has been shown to have potent morphogenetic effects. We hypothesize that these effects are at work in lung morphogenesis. In order to refine that broad hypothesis in a quantitative framework, we geometrically analyse cell shapes in an epithelial tissue, and individual cell deformations resulting from peristaltic waves that completely occlude the airway. Typical distortions can be very large, with opposite orientations in the stalk and tip regions. Apical distortions are always greater than basal distortions. We give a quantitative estimate of the relationship between length of occluded airway and the resulting tissue stretch in the distal tip. We refine our analysis of cell stresses and strains from peristalsis with a simple mechanical model of deformation of cells within an epithelium, which accounts for basic subcellular geometry and material properties. The model identifies likely stress concentrations near the nucleus and at the apical cell–cell junction. The surprisingly large strains of airway peristalsis may serve to rearrange cells and stimulate other mechanosensitive processes by repeatedly aligning cytoskeletal components and/or breaking and reforming lateral cell–cell adhesions. Stress concentrations between nuclei of adjacent cells may serve as a mechanical control mechanism guiding the alignment of nuclei as an epithelium matures.

Convergent extension by intercalation without mediolaterally fixed cell motion

We construct and implement a stochastic model of convergent extension, using a minimal set of assumptions on cell behavior. In addition to the basic assumptions of volume conservation, random cell motion, and cell-cell and cell-ECM adhesion, and a non-standard assumption that cytoskeletal polymerization generates an internal pressure tending to keep cells convex, we find that we need only two conditions for convergent extension. (1) Each cell type has a particular aspect ratio towards which it regulates its geometry. We do not require that cells align in a specific orientation, e.g. to be oriented mediolaterally. (2) The elongating tissue is composed of cells that prefer to be elongated, and these cells must be accompanied by cells which prefer to be round. The latter effectively provide a boundary to capture. In simulations, our model tissue extends and converges to a stacked arrangement of elongated cells one cell wide, an arrangement which is seen in ascidian notochords, but which has not been observed in other models. This arrangement is achieved without any direct mediolateral bias other than that which is provided by the physical edge of the adjacent tissue.

Modelling the spatial patterning of teeth primordia in the alligator

We propose a model mechanism for the initiation and spatial positioning of teeth primordia in the alligator, Alligator mississippiensis. Detailed embryological studies by Westergaard and Ferguson (1986, 1987, 1990) have shown that jaw growth plays a crucial role in the developmental patterning of the tooth initiation process. Based on biological data we develop a dynamic patterning mechanism, which crucially includes domain growth. The mechanism can reproduce the spatial pattern development of the first seven teeth primordia in each half jaw of A. mississippiensis. The results for the precise spatio-temporal sequence compare well with experiment. Simulation of the model also predicts that certain transplantations can alter the spatial sequence of teeth primordia initiation.

Tissue geometry may govern lung branching mode selection

Lung branching morphogenesis proceeds in three stereotyped modes (domain, planar, and orthogonal branching). Much is known about the molecular players, including growth factors such as fibroblast growth factor 10 but it is unknown how these signals could actuate the different branching patterns. With the aim of identifying mechanisms that may determine the different branching modes, we developed a computational model of the epithelial lung bud and its surrounding mesenchyme. We studied transport of morphogens and localization of morphogen flux at lobe surfaces and lobe edges. We find that a single simple mechanism is theoretically capable of directing an epithelial tubule to elongate, bend, flatten, or bifurcate, depending solely on geometric ratios of the tissues in the vicinity of a growing tubule tip. Furthermore, the same simple mechanism is capable of generating orthogonal or planar branching, depending only on the same geometric ratios.