Static and dynamic stress heterogeneity in a multiscale model of the asthmatic airway wall.

Airway hyperresponsiveness (AHR) is a key characteristic of asthma that remains poorly understood. Tidal breathing and deep inspiration ordinarily cause rapid relaxation of airway smooth muscle (ASM) (as demonstrated via application of length fluctuations to tissue strips) and are therefore implicated in modulation of AHR, but in some cases (such as application of transmural pressure oscillations to isolated intact airways) this mechanism fails. Here we use a multiscale biomechanical model for intact airways that incorporates strain stiffening due to collagen recruitment and dynamic force generation by ASM cells to show that the geometry of the airway, together with interplay between dynamic active and passive forces, gives rise to large stress and compliance heterogeneities across the airway wall that are absent in tissue strips. We show further that these stress heterogeneities result in auxotonic loading conditions that are currently not replicated in tissue-strip experiments; stresses in the strip are similar to hoop stress only at the outer airway wall and are under- or overestimates of stresses at the lumen. Taken together these results suggest that a previously underappreciated factor, stress heterogeneities within the airway wall and consequent ASM cellular response to this micromechanical environment, could contribute to AHR and should be explored further both theoretically and experimentally.

The role of contractile unit reorganization in force generation in airway smooth muscle.

Airway smooth muscle (ASM) cells undergo remodelling and reside in a tissue structure that is subject to heterogenous stress distributions that change dynamically during the breathing cycle. In this paper, we develop a structural model of an ASM cell that consists of contractile units (actin and myosin filaments) in series and parallel, anchored to a nonlinearly elastic cytoskeleton. We mimic a typical experimental protocol that involves isometric force generation through triggering of the contractile machinery, followed by oscillatory length fluctuation of the cell. We use the model to predict the effect of a single instance of rearrangement of the contractile machinery, combined with strain-stiffening of the cytoskeleton, on the force generated by the sarcomeres, and the total force generated by the cell. By linking intra-cellular events to whole-cell behaviour, the model reveals mechanistic relationships between structural properties and cell-level force-length loops. We show how contractile force, shortening velocity and sarcomere operating lengths vary as the internal cell architecture is altered. Additionally, we show how interactions between the internal contractile machinery and cytoskeletal structure play a role in the regulation of force generation and hysteresis of the cell.

Buckling of a growing tissue and the emergence of two-dimensional patterns.

The process of biological growth and the associated generation of residual stress has previously been considered as a driving mechanism for tissue buckling and pattern selection in numerous areas of biology. Here, we develop a two-dimensional thin plate theory to simulate the growth of cultured intestinal epithelial cells on a deformable substrate, with the goal of elucidating how a tissue engineer might best recreate the regular array of invaginations (crypts of Lieberkühn) found in the wall of the mammalian intestine. We extend the standard von Kármán equations to incorporate inhomogeneity in the plate's mechanical properties and surface stresses applied to the substrate by cell proliferation. We determine numerically the configurations of a homogeneous plate under uniform cell growth, and show how tethering to an underlying elastic foundation can be used to promote higher-order buckled configurations. We then examine the independent effects of localised softening of the substrate and spatial patterning of cellular growth, demonstrating that (within a two-dimensional framework, and contrary to the predictions of one-dimensional models) growth patterning constitutes a more viable mechanism for control of crypt distribution than does material inhomogeneity.

A biomechanical model of anther opening reveals the roles of dehydration and secondary thickening.

Understanding the processes that underlie pollen release is a prime target for controlling fertility to enable selective breeding and the efficient production of hybrid crops. Pollen release requires anther opening, which involves changes in the biomechanical properties of the anther wall. In this research, we develop and use a mathematical model to understand how these biomechanical processes lead to anther opening. Our mathematical model describing the biomechanics of anther opening incorporates the bilayer structure of the mature anther wall, which comprises the outer epidermal cell layer, whose turgor pressure is related to its hydration, and the endothecial layer, whose walls contain helical secondary thickening, which resists stretching and bending. The model describes how epidermal dehydration, in association with the thickened endothecial layer, creates forces within the anther wall causing it to bend outwards, resulting in anther opening and pollen release. The model demonstrates that epidermal dehydration can drive anther opening, and suggests why endothecial secondary thickening is essential for this process (explaining the phenotypes presented in the myb26 and nst1nst2 mutants). The research hypothesizes and demonstrates a biomechanical mechanism for anther opening, which appears to be conserved in many other biological situations where tissue movement occurs.

A model of crosslink kinetics in the expanding plant cell wall: yield stress and enzyme action.

The plant primary cell wall is a composite material containing stiff cellulose microfibrils that are embedded within a pectin matrix and crosslinked through a network of hemicellulose polymers. This microstructure endows the wall with nonlinear anisotropic mechanical properties and allows enzymatic regulation of expansive cell growth. We present a mathematical model of hemicellulose crosslink dynamics in an expanding cell wall incorporating strain-enhanced breakage and enzyme-mediated crosslink kinetics. The model predicts the characteristic yielding behaviour in the relationship between stress and strain-rate seen experimentally, and suggests how the effective yield and extensibility of the wall depend on microstructural parameters and on the action of enzymes of the XTH and expansin families. The model suggests that the yielding behaviour encapsulated in the classical Lockhart equation can be explained by the strongly nonlinear dependence of crosslink breakage rate on crosslink elongation. The model also demonstrates how enzymes that target crosslink binding can be effective in softening the wall in its pre-yield state, whereas its post-yield extensibility is determined primarily by the pectin matrix.

Growth-induced buckling of an epithelial layer.

We use a proof-of-concept experiment and two mathematical models to explore growth-induced tissue buckling, as may occur in colorectal crypt formation. Our experiment reveals how growth of a cultured epithelial monolayer on a thin flexible substrate can cause out-of-plane substrate deflections. We describe this system theoretically using a 'bilayer' model in which a growing cell layer adheres to a thin compressible elastic beam. We compare this with the 'supported-monolayer' model due to Edwards and Chapman (Bull Math Biol 69:1927-1942, 2007) for an incompressible expanding beam (representing crypt epithelium), which incorporates viscoelastic tethering to underlying stroma. We show that the bilayer model can exhibit buckling via parametric growth (in which the system passes through a sequence of equilibrium states, parameterised by the total beam length); in this case, non-uniformities in cell growth and variations in cell-substrate adhesion are predicted to have minimal effect on the shape of resulting buckled states. The supported-monolayer model reveals how competition between lateral supports and stromal adhesion influences the wavelength of buckled states (in parametric growth), and how non-equilibrium relaxation of tethering forces influences post-buckled shapes. This model also predicts that non-uniformities in growth patterns have a much weaker influence on buckled shapes than non-uniformities in material properties. Together, the experiment and models support the concept of patterning by growth-induced buckling and suggest that targeted softening of a growing cell layer provides greater control in shaping tissues than non-uniform growth.

A hybrid approach to multi-scale modelling of cancer.

In this paper, we review multi-scale models of solid tumour growth and discuss a middle-out framework that tracks individual cells. By focusing on the cellular dynamics of a healthy colorectal crypt and its invasion by mutant, cancerous cells, we compare a cell-centre, a cell-vertex and a continuum model of cell proliferation and movement. All models reproduce the basic features of a healthy crypt: cells proliferate near the crypt base, they migrate upwards and are sloughed off near the top. The models are used to establish conditions under which mutant cells are able to colonize the crypt either by top-down or by bottom-up invasion. While the continuum model is quicker and easier to implement, it can be difficult to relate system parameters to measurable biophysical quantities. Conversely, the greater detail inherent in the multi-scale models means that experimentally derived parameters can be incorporated and, therefore, these models offer greater scope for understanding normal and diseased crypts, for testing and identifying new therapeutic targets and for predicting their impacts.

Continuum approximations of individual-based models for epithelial monolayers.

This work examines a 1D individual-based model (IBM) for a system of tightly adherent cells, such as an epithelial monolayer. Each cell occupies a bounded region, defined by the location of its endpoints, has both elastic and viscous mechanical properties and is subject to drag generated by adhesion to the substrate. Differential-algebraic equations governing the evolution of the system are obtained from energy considerations. This IBM is then approximated by continuum models (systems of partial differential equations) in the limit of a large number of cells, N, when the cell parameters vary slowly in space or are spatially periodic (and so may be heterogeneous, with substantial variation between adjacent cells). For spatially periodic cell properties with significant cell viscosity, the relationship between the mean cell pressure and length for the continuum model is found to be history dependent. Terms involving convective derivatives, not normally included in continuum tissue models, are identified. The specific problem of the expansion of an aggregate of cells through cell growth (but without division) is considered in detail, including the long-time and slow-growth-rate limits. When the parameters of neighbouring cells vary slowly in space, the O(1/N(2)) error in the continuum approximation enables this approach to be used even for modest values of N. In the spatially periodic case, the neglected terms are found to be O(1/N). The model is also used to examine the acceleration of a wound edge observed in wound-healing assays.

A biomechanical model of agonist-initiated contraction in the asthmatic airway.

This paper presents a modelling framework in which the local stress environment of airway smooth muscle (ASM) cells may be predicted and cellular responses to local stress may be investigated. We consider an elastic axisymmetric model of a layer of connective tissue and circumferential ASM fibres embedded in parenchymal tissue and model the active contractile force generated by ASM via a stress acting along the fibres. A constitutive law is proposed that accounts for active and passive material properties as well as the proportion of muscle to connective tissue. The model predicts significantly different contractile responses depending on the proportion of muscle to connective tissue in the remodelled airway. We find that radial and hoop-stress distributions in remodelled muscle layers are highly heterogenous with distinct regions of compression and tension. Such patterns of stress are likely to have important implications, from a mechano-transduction perspective, on contractility, short-term cytoskeletal adaptation and long-term airway remodelling in asthma.

A mathematical model of intervillous blood flow in the human placentone.

We present a mathematical model for maternal blood flow in a placental circulatory unit (a placentone), describing flow of maternal blood via Darcy's law and steady advective transport of a dissolved nutrient. The method of images and computational integration along streamlines are employed to find flow and solute concentration distributions, which are illustrated for a range of governing system parameters. The model shows how the calibre of the basal vessels can be a dominant determinant of the maternal blood flow rate through the placentone, given a driving pressure difference between the spiral arteries and decidual veins. The model supports the hypothesis that basal veins are located on the periphery of the placentone in order to optimise delivery of nutrients and suggests the existence of an optimal volume fraction of villous tissue.

An integrative computational model for intestinal tissue renewal.

OBJECTIVES: The luminal surface of the gut is lined with a monolayer of epithelial cells that acts as a nutrient absorptive engine and protective barrier. To maintain its integrity and functionality, the epithelium is renewed every few days. Theoretical models are powerful tools that can be used to test hypotheses concerning the regulation of this renewal process, to investigate how its dysfunction can lead to loss of homeostasis and neoplasia, and to identify potential therapeutic interventions. Here we propose a new multiscale model for crypt dynamics that links phenomena occurring at the subcellular, cellular and tissue levels of organisation. METHODS: At the subcellular level, deterministic models characterise molecular networks, such as cell-cycle control and Wnt signalling. The output of these models determines the behaviour of each epithelial cell in response to intra-, inter- and extracellular cues. The modular nature of the model enables us to easily modify individual assumptions and analyse their effects on the system as a whole. RESULTS: We perform virtual microdissection and labelling-index experiments, evaluate the impact of various model extensions, obtain new insight into clonal expansion in the crypt, and compare our predictions with recent mitochondrial DNA mutation data. CONCLUSIONS: We demonstrate that relaxing the assumption that stem-cell positions are fixed enables clonal expansion and niche succession to occur. We also predict that the presence of extracellular factors near the base of the crypt alone suffices to explain the observed spatial variation in nuclear beta-catenin levels along the crypt axis.

Three-dimensional elastohydrodynamics of a thin plate oscillating above a wall.

We consider deflections of a thin rectangular elastic plate that is submerged within a Newtonian fluid. The plate is clamped along one edge and supported horizontally over a plane horizontal wall. We consider both external driving, where the clamped edge is vibrated vertically at high frequencies, and thermal driving, where the plate fluctuates under Brownian motion. In both cases, the amplitude of oscillation is assumed sufficiently small that the resulting flow has little convective inertia, although the oscillation frequency is sufficiently high to generate substantial unsteady inertia in the flow, a common scenario in many nano- and microdevices. We exploit the plate's thinness to develop an integral-equation representation for the three-dimensional flow (a so-called thin-plate theory) which offers considerable computational savings over a full boundary-integral formulation. Limiting cases of high oscillation frequencies and small wall-plate separation distances are studied separately, leading to further simplified descriptions for the hydrodynamics. We validate these reduced integral representations against full boundary-integral computations, and identify the parameter ranges over which these simplified formulations are valid. Addressing the full flow-structure interaction, we also examine the limits of simpler two-dimensional hydrodynamic models. We compare the responses of a narrow plate under two- and three-dimensional hydrodynamic loading, and report differences in the frequency response curves that occur when the plate operates in water, in contrast to the excellent agreement observed in air.

Shock formation and non-linear dispersion in a microvascular capillary network.

Temporal and spatial fluctuations are a common feature of blood flow in microvascular networks. Among many possible causes, previous authors have suggested that the non-linear rheological properties of capillary blood flow (notably the Fåhraeus effect, the Fåhraeus-Lindqvist effect and the phase-separation effect at bifurcations) may be sufficient to generate temporal fluctuations even in very simple networks. We have simulated blood flow driven by a fixed pressure drop through a simple arcade network using coupled hyperbolic partial differential equations (PDEs) that incorporate well-established empirical descriptions of these rheological effects, accounting in particular for spatially varying haematocrit distributions; we solved the PDE system using a characteristic-based method. Our computations indicate that, under physiologically realistic conditions, there is a unique steady flow in an arcade network which is linearly stable and that plasma skimming suppresses the oscillatory decay of perturbations. In addition, we find that non-linear perturbations to haematocrit distributions can develop shocks via the Fåhraeus effect, providing a novel mechanism for non-linear dispersion in microvascular networks.

T-cell motility in the early stages of the immune response modeled as a random walk amongst targets.

The transport process by which a T cell makes high-frequency encounters with antigen-presenting cells following infection is an important element of adaptive immunity. Recent experimental work has allowed in vivo cell motility to be characterized in detail. On the basis of experimental data we develop a quantitative model for encounters between T cells and antigen-presenting cells. We model this as a transport-limited chemical reaction with the dynamics dependent on physical contact between randomly moving reactants. We use asymptotic methods to calculate a time distribution which characterizes the delay before a T cell is activated and use Monte Carlo simulations to verify the analysis. We find that the density of antigen-primed dendritic cells within the lymph node paracortex must be greater than 35 cells/mm3 for a T cell to have a more than 50% chance of encountering a dendritic cell within 24 h. This density is much larger than existing estimates based on calculations which neglect the transport process. We also use simulations to compare a T cell which re-orients isotropically with a T cell which turns according to an experimentally observed distribution and find that the effects of anisotropy on the solution are small.

Crypt dynamics and colorectal cancer: advances in mathematical modelling.

Mathematical modelling forms a key component of systems biology, offering insights that complement and stimulate experimental studies. In this review, we illustrate the role of theoretical models in elucidating the mechanisms involved in normal intestinal crypt dynamics and colorectal cancer. We discuss a range of modelling approaches, including models that describe cell proliferation, migration, differentiation, crypt fission, genetic instability, APC inactivation and tumour heterogeneity. We focus on the model assumptions, limitations and applications, rather than on the technical details. We also present a new stochastic model for stem-cell dynamics, which predicts that, on average, APC inactivation occurs more quickly in the stem-cell pool in the absence of symmetric cell division. This suggests that natural niche succession may protect stem cells against malignant transformation in the gut. Finally, we explain how we aim to gain further understanding of the crypt system and of colorectal carcinogenesis with the aid of multiscale models that cover all levels of organization from the molecular to the whole organ.

Stochastic elastohydrodynamics of a microcantilever oscillating near a wall.

We consider the thermally driven motion of a microcantilever in a fluid environment near a wall, a configuration characteristic of the atomic force microscope. A theoretical model is presented which accounts for hydrodynamic interactions between the cantilever and wall over a wide range of frequencies and which exploits the fluctuation-dissipation theorem to capture the Brownian dynamics of the coupled fluid-cantilever system. Model predictions are tested against experimental thermal spectra for a cantilever in air and water. The model shows how, in a liquid environment, the effects of non-delta-correlated Brownian forcing appear in the power spectrum, particularly at low frequencies. The model also predicts accurately changes in the spectrum in liquid arising through hydrodynamic wall effects, which we show are strongly mediated by the angle at which the cantilever is tilted relative to the wall.

An asymptotic model of unsteady airway reopening.

We consider a simple physical model for the reopening of a collapsed lung airway involving the unsteady propagation of a long bubble of air, driven at a prescribed flow-rate, into a liquid-filled channel formed by two flexible membranes that are held under large longitudinal tension and are confined between two parallel rigid plates. This system is described theoretically using an asymptotic approximation, valid for uniformly small membrane slopes, which reduces to a fourth-order nonlinear evolution equation for the channel width ahead of the bubble tip, from which the time-evolution of the bubble pressure pb* and bubble speed may be determined. The model shows that there can be a substantial delay between the time at which the bubble starts to grow in volume and the time at which its tip starts to move. Under certain conditions, the start of the bubble's motion is accompanied by a transient overshoot in pb*, as seen previously in experiment; the model predicts that the overshoot is greatest in narrow channels when the bubble is driven with a large volume flux. It is also shown how the threshold pressure for steady bubble propagation in wide channels has distinct contributions from the capillary pressure drop across the bubble tip and viscous dissipation in the channel ahead of the bubble.

Capillary-elastic instabilities of liquid-lined lung airways.

To model the competition between capillary and elastic forces in controlling the shape of a small lung airway and its interior liquid lining, we compute the equilibrium configurations of a liquid-lined, externally pressurized, buckled elastic tube. We impose axial uniformity and assume that the liquid wets the tube wall with zero contact angle. Non-zero surface tension has a profound effect on the tube's quasi-steady inflation-deflation characteristics. At low liquid volumes, hysteresis arises through two distinct mechanisms, depending on the buckling wavenumber. Sufficient compression always leads to abrupt and irreversible collapse and flooding of the tube; flooding is promoted by increasing liquid volumes or surface tension. The model captures mechanisms whereby capillary-elastic instabilities can lead to airway closure.

Surfactant transport over airway liquid lining of nonuniform depth.

Numerous effects (e.g., airway wall buckling, gravity, airway curvature, capillary instabilities) give rise to nonuniformities in the depth of the liquid lining of peripheral lung airways. The effects of such thickness variations on the unsteady spreading of a surfactant monolayer along an airway are explored theoretically here. Flow-induced film deformations are shown to have only a modest influence on spreading rates, motivating the use of a simplified model in which the liquid-lining depth is prescribed and the monolayer concentration satisfies a spatially inhomogeneous nonlinear diffusion equation. Two generic situations are considered: spreading along a continuous annular liquid lining of nonuniform depth, and spreading along a rivulet that wets the airway wall with zero contact angle. In both cases, transverse averaging at large times yields a one-dimensional approximation of axial spreading that is valid for the majority of the monolayer. However, a localized monolayer remains persistently two dimensional in a region at its leading edge having axial length scales comparable to the length scale of transverse depth variation. It is also shown how the transverse spreading of a monolayer may be arrested as it approaches a static contact line at the edge of a rivulet. Implications for Surfactant Replacement Therapy are discussed.

A theoretical study of surfactant and liquid delivery into the lung.

A computational study is presented for the transport of liquids and insoluble surfactant through the lung airways, delivered from a source at the distal end of the trachea. Four distinct transport regimes are considered: 1) the instilled bolus may create a liquid plug that occludes the large airways but is forced peripherally during mechanical ventilation; 2) the bolus creates a deposited film on the airway walls, either from the liquid plug transport or from direct coating, that drains under the influence of gravity through the first few airway generations; 3) in smaller airways, surfactant species form a surface layer that spreads due to surface-tension gradients, i.e., Marangoni flows; and 4) the surfactant finally reaches the alveolar compartment where it is cleared according to first-order kinetics. The time required for a quasi-steady-state transport process to evolve and for the subsequent delivery of the dose is predicted. Following fairly rapid transients, on the order of seconds, steady-state transport develops and is governed by the interaction of Marangoni flow and alveolar kinetics. Total delivery time is approximately 24 h for a typical first dose. Numerical solutions show that both transit and delivery times are strongly influenced by the strength of the preexisting surfactant and the geometric properties of the airway network. Delivery times for follow-up doses can increase significantly as the level of preexisting surfactant rises.