The influence of hydrostatic pressure on tissue engineered bone development.

The hydrostatic pressure stimulation of an appropriately cell-seeded porous scaffold within a bioreactor is a promising method for engineering bone tissue external to the body. We propose a mathematical model, and employ a suite of candidate constitutive laws, to qualitatively describe the effect of applied hydrostatic pressure on the quantity of minerals deposited in such an experimental setup. By comparing data from numerical simulations with experimental observations under a number of stimulation protocols, we suggest that the response of bone cells to an applied pressure requires consideration of two components; (i) a component describing the cell memory of the applied stimulation, and (ii) a recovery component, capturing the time cells require to recover from high rates of mineralisation.

An investigation of the influence of extracellular matrix anisotropy and cell-matrix interactions on tissue architecture.

Mechanical interactions between cells and the fibrous extracellular matrix (ECM) in which they reside play a key role in tissue development. Mechanical cues from the environment (such as stress, strain and fibre orientation) regulate a range of cell behaviours, including proliferation, differentiation and motility. In turn, the ECM structure is affected by cells exerting forces on the matrix which result in deformation and fibre realignment. In this paper we develop a mathematical model to investigate this mechanical feedback between cells and the ECM. We consider a three-phase mixture of collagen, culture medium and cells, and formulate a system of partial differential equations which represents conservation of mass and momentum for each phase. This modelling framework takes into account the anisotropic mechanical properties of the collagen gel arising from its fibrous microstructure. We also propose a cell-collagen interaction force which depends upon fibre orientation and collagen density. We use a combination of numerical and analytical techniques to study the influence of cell-ECM interactions on pattern formation in tissues. Our results illustrate the wide range of structures which may be formed, and how those that emerge depend upon the importance of cell-ECM interactions.

On a poroviscoelastic model for cell crawling.

In this paper a minimal, one-dimensional, two-phase, viscoelastic, reactive, flow model for a crawling cell is presented. Two-phase models are used with a variety of constitutive assumptions in the literature to model cell motility. We use an upper-convected Maxwell model and demonstrate that even the simplest of two-phase, viscoelastic models displays features relevant to cell motility. We also show care must be exercised in choosing parameters for such models as a poor choice can lead to an ill-posed problem. A stability analysis reveals that the initially stationary, spatially uniform strip of cytoplasm starts to crawl in response to a perturbation which breaks the symmetry of the network volume fraction or network stress. We also demonstrate numerically that there is a steady travelling-wave solution in which the crawling velocity has a bell-shaped dependence on adhesion strength, in agreement with biological observation.

Non-local models for the formation of hepatocyte-stellate cell aggregates.

Liver cell aggregates may be grown in vitro by co-culturing hepatocytes with stellate cells. This method results in more rapid aggregation than hepatocyte-only culture, and appears to enhance cell viability and the expression of markers of liver-specific functions. We consider the early stages of aggregate formation, and develop a new mathematical model to investigate two alternative hypotheses (based on evidence in the experimental literature) for the role of stellate cells in promoting aggregate formation. Under Hypothesis 1, each population produces a chemical signal which affects the other, and enhanced aggregation is due to chemotaxis. Hypothesis 2 asserts that the interaction between the two cell types is by direct physical contact: the stellates extend long cellular processes which pull the hepatocytes into the aggregates. Under both hypotheses, hepatocytes are attracted to a chemical they themselves produce, and the cells can experience repulsive forces due to overcrowding. We formulate non-local (integro-partial differential) equations to describe the densities of cells, which are coupled to reaction-diffusion equations for the chemical concentrations. The behaviour of the model under each hypothesis is studied using a combination of linear stability analysis and numerical simulations. Our results show how the initial rate of aggregation depends upon the cell seeding ratio, and how the distribution of cells within aggregates depends on the relative strengths of attraction and repulsion between the cell types. Guided by our results, we suggest experiments which could be performed to distinguish between the two hypotheses.

The influence of bioreactor geometry and the mechanical environment on engineered tissues.

A three phase model for the growth of a tissue construct within a perfusion bioreactor is examined. The cell population (and attendant extracellular matrix), culture medium, and porous scaffold are treated as distinct phases. The bioreactor system is represented by a two-dimensional channel containing a cell-seeded rigid porous scaffold (tissue construct), which is perfused with a culture medium. Through the prescription of appropriate functional forms for cell proliferation and extracellular matrix deposition rates, the model is used to compare the influence of cell density-, pressure-, and culture medium shear stress-regulated growth on the composition of the engineered tissue. The governing equations are derived in O'Dea et al. "A Three Phase Model for Tissue Construct Growth in a Perfusion Bioreactor," Math. Med. Biol., in which the long-wavelength limit was exploited to aid analysis; here, finite element methods are used to construct two-dimensional solutions to the governing equations and to investigate thoroughly their behavior. Comparison of the total tissue yield and averaged pressures, velocities, and shear stress demonstrates that quantitative agreement between the two-dimensional and long-wavelength approximation solutions is obtained for channel aspect ratios of order 10(-2) and that much of the qualitative behavior of the model is captured in the long-wavelength limit, even for relatively large channel aspect ratios. However, we demonstrate that in order to capture accurately the effect of mechanotransduction mechanisms on tissue construct growth, spatial effects in at least two dimensions must be included due to the inherent spatial variation of mechanical stimuli relevant to perfusion bioreactors, most notably, fluid shear stress, a feature not captured in the long-wavelength limit.

Multiscale modelling and homogenisation of fibre-reinforced hydrogels for tissue engineering.

Tissue engineering aims to grow artificial tissues in vitro to replace those in the body that have been damaged through age, trauma or disease. A recent approach to engineer artificial cartilage involves seeding cells within a scaffold consisting of an interconnected 3D-printed lattice of polymer fibres combined with a cast or printed hydrogel, and subjecting the construct (cell-seeded scaffold) to an applied load in a bioreactor. A key question is to understand how the applied load is distributed throughout the construct. To address this, we employ homogenisation theory to derive equations governing the effective macroscale material properties of a periodic, elastic-poroelastic composite. We treat the fibres as a linear elastic material and the hydrogel as a poroelastic material, and exploit the disparate length scales (small inter-fibre spacing compared with construct dimensions) to derive macroscale equations governing the response of the composite to an applied load. This homogenised description reflects the orthotropic nature of the composite. To validate the model, solutions from finite element simulations of the macroscale, homogenised equations are compared to experimental data describing the unconfined compression of the fibre-reinforced hydrogels. The model is used to derive the bulk mechanical properties of a cylindrical construct of the composite material for a range of fibre spacings and to determine the local mechanical environment experienced by cells embedded within the construct.

Pattern formation in multiphase models of chemotactic cell aggregation.

We develop a continuum model for the aggregation of cells cultured in a nutrient-rich medium in a culture well. We consider a 2D geometry, representing a vertical slice through the culture well, and assume that the cell layer depth is small compared with the typical lengthscale of the culture well. We adopt a continuum mechanics approach, treating the cells and culture medium as a two-phase mixture. Specifically, the cells and culture medium are treated as fluids. Additionally, the cell phase can generate forces in response to environmental cues, which include the concentration of a chemoattractant that is produced by the cells within the culture medium. The model leads to a system of coupled nonlinear partial differential equations for the volume fraction and velocity of the cell phase, the culture medium pressure and the chemoattractant concentration, which must be solved subject to appropriate boundary and initial conditions. To gain insight into the system, we consider two model reductions, appropriate when the cell layer depth is thin compared to the typical length scale of the culture well: a (simple) 1D and a (more involved) thin-film extensional flow reduction. By investigating the resulting systems of equations analytically and numerically, we identify conditions under which small amplitude perturbations to a homogeneous steady state (corresponding to a spatially uniform cell distribution) can lead to a spatially varying steady state (pattern formation). Our analysis reveals that the simpler 1D reduction has the same qualitative features as the thin-film extensional flow reduction in the linear and weakly nonlinear regimes, motivating the use of the simpler 1D modelling approach when a qualitative understanding of the system is required. However, the thin-film extensional flow reduction may be more appropriate when detailed quantitative agreement between modelling predictions and experimental data is desired. Furthermore, full numerical simulations of the two model reductions in regions of parameter space when the system is not close to marginal stability reveal significant differences in the evolution of the volume fraction and velocity of the cell phase, and chemoattractant concentration.

Synchronization between arterial blood pressure and cerebral oxyhaemoglobin concentration investigated by wavelet cross-correlation.

Wavelet cross-correlation (WCC) is used to analyse the relationship between low-frequency oscillations in near-infrared spectroscopy (NIRS) measured cerebral oxyhaemoglobin (O(2)Hb) and mean arterial blood pressure (MAP) in patients suffering from autonomic failure and age-matched controls. Statistically significant differences are found in the wavelet scale of maximum cross-correlation upon posture change in patients, but not in controls. We propose that WCC analysis of the relationship between O(2)Hb and MAP provides a useful method of investigating the dynamics of cerebral autoregulation using the spontaneous low-frequency oscillations that are typically observed in both variables without having to make the assumption of stationarity of the time series. It is suggested that for a short-duration clinical test previous transfer-function-based approaches to analyse this relationship may suffer due to the inherent nonstationarity of low-frequency oscillations that are observed in the resting brain.

Kinetic effects regularize the mass-flux singularity at the contact line of a thin evaporating drop.

We consider the transport of vapour caused by the evaporation of a thin, axisymmetric, partially wetting drop into an inert gas. We take kinetic effects into account through a linear constitutive law that states that the mass flux through the drop surface is proportional to the difference between the vapour concentration in equilibrium and that at the interface. Provided that the vapour concentration is finite, our model leads to a finite mass flux in contrast to the contact-line singularity in the mass flux that is observed in more standard models that neglect kinetic effects. We perform a local analysis near the contact line to investigate the way in which kinetic effects regularize the mass-flux singularity at the contact line. An explicit expression is derived for the mass flux through the free surface of the drop. A matched-asymptotic analysis is used to further investigate the regularization of the mass-flux singularity in the physically relevant regime in which the kinetic timescale is much smaller than the diffusive one. We find that the effect of kinetics is limited to an inner region near the contact line, in which kinetic effects enter at leading order and regularize the mass-flux singularity. The inner problem is solved explicitly using the Wiener-Hopf method and a uniformly valid composite expansion is derived for the mass flux in this asymptotic limit.

A mathematical model for cell infiltration and proliferation in a chondral defect.

We develop a mathematical model to describe the regeneration of a hydrogel inserted into an ex vivo osteochondral explant. Specifically we use partial differential equations to describe the evolution of two populations of cells that migrate from the tissue surrounding the defect, proliferate, and compete for space and resources within the hydrogel. The two cell populations are chondrocytes and cells that infiltrate from the subchondral bone. Model simulations are used to investigate how different seeding strategies and growth factor placement within the hydrogel affect the spatial distribution of both cell types. Since chondrocyte migration is extremely slow, we conclude that the hydrogel should be seeded with chondrocytes prior to culture in order to obtain zonal chondrocyte distributions typical of those associated with healthy cartilage.

A tidal breathing model for the multiple inert gas elimination technique.

The tidal breathing lung model described for the sine-wave technique (D. J. Gavaghan and C. E. W. Hahn. Respir. Physiol. 106: 209-221, 1996) is generalized to continuous ventilation-perfusion and ventilation-volume distributions. This tidal breathing model is then applied to the multiple inert gas elimination technique (P. D. Wagner, H. A. Saltzman, and J. B. West. J. Appl. Physiol. 36: 588-599, 1974). The conservation of mass equations are solved, and it is shown that 1) retentions vary considerably over the course of a breath, 2) the retentions are dependent on alveolar volume, and 3) the retentions depend only weakly on the width of the ventilation-volume distribution. Simulated experimental data with a unimodal ventilation-perfusion distribution are inserted into the parameter recovery model for a lung with 1 or 2 alveolar compartments and for a lung with 50 compartments. The parameters recovered using both models are dependent on the time interval over which the blood sample is taken. For best results, the blood sample should be drawn over several breath cycles.

The effect of diffusion in the respiratory tree on the alveolar amplitude response technique (AART).

Theoretical data for the alveolar amplitude response technique (AART) (J. Appl. Physiol. 41 (1976) 419-424) for assessing lung function was simulated using a single path lung model. This model takes account of stratified inhomogeneities in gas concentrations within the respiratory tree. The data was inserted into previously published parameter recovery techniques that may be used to estimate dead-space volume, alveolar volume and cardiac output. These parameter recovery techniques are based on much simpler mathematical models that do not allow stratified inhomogeneities in gas concentrations. It was found that: (i) recovered dead-space volume depended significantly on the ventilation pattern and on the distribution of volume within of the conducting airways; (ii) alveolar volume was recovered to a good degree of accuracy; and (iii) the recovered value of cardiac output was highly dependent on both the choice of inert gas and parameter recovery technique.

A tidal ventilation model for oxygenation in respiratory failure.

We develop tidal-ventilation pulmonary gas-exchange equations that allow pulmonary shunt to have different values during expiration and inspiration, in accordance with lung collapse and recruitment during lung dysfunction (Am. J. Respir. Crit. Care Med. 158 (1998) 1636). Their solutions are tested against published animal data from intravascular oxygen tension and saturation sensors. These equations provide one explanation for (i) observed physiological phenomena, such as within-breath fluctuations in arterial oxygen saturation and blood-gas tension; and (ii) conventional (time averaged) blood-gas sample oxygen tensions. We suggest that tidal-ventilation models are needed to describe within-breath fluctuations in arterial oxygen saturation and blood-gas tension in acute respiratory distress syndrome (ARDS) subjects. Both the amplitude of these oxygen saturation and tension fluctuations, and the mean oxygen blood-gas values, are affected by physiological variables such as inspired oxygen concentration, lung volume, and the inspiratory:expiratory (I:E) ratio, as well as by changes in pulmonary shunt during the respiratory cycle.

The effects of ventilation pattern on carbon dioxide transfer in three computer models of the airways.

We investigate the effects on arterial P(CO(2)) and on arterial-end tidal P(CO(2)) difference of six different ventilation patterns of equal tidal volume, and also of various combinations of tidal volume and respiratory rate that maintain a constant alveolar ventilation. We use predictions from three different mathematical models. Models 1 (distributed) and 2 (compartmental) include combined convection and diffusion effects. Model 3 incorporates a single well-mixed alveolar compartment and an anatomical dead-space in which plug flow occurs. We found that: (i) breathing patterns with longer inspiratory times yield lower arterial P(CO(2)); (ii) varying tidal volume and respiratory rate so that alveolar ventilation is kept constant may change both PA(CO(2)) and the PA(CO(2))-PET(CO(2)) difference; (iii) the distributed model predicts higher end-tidal and arterial P(CO(2)) than the compartmental models under similar conditions; and (iv) P(CO(2)) capnograms predicted by the distributed model exhibit longer phase I and steeper phase II than other models.

Some factors affecting oxygen uptake by red blood cells in the pulmonary capillaries.

In this study we investigate the equations governing the transport of oxygen in pulmonary capillaries. We use a mathematical model consisting of a red blood cell completely surrounded by plasma within a cylindrical pulmonary capillary. This model takes account of convection and diffusion of oxygen through plasma, diffusion of oxygen through the red blood cell, and the reaction between oxygen and haemoglobin molecules. The velocity field within the plasma is calculated by solving the slow flow equations. We investigate the effect on the solution of the governing equations of: (i) mixed-venous blood oxygen partial pressure (the initial conditions); (ii) alveolar gas oxygen partial pressure (the boundary conditions); (iii) neglecting the convection term; and (iv) assuming an instantaneous reaction between the oxygen and haemoglobin molecules. It is found that: (a) equilibrium is reached much more rapidly for high values of mixed-venous blood and alveolar gas oxygen partial pressure; (b) the convection term has a negligible effect on the time taken to reach a prescribed degree of equilibrium; and (c) an instantaneous reaction may be assumed. Explanations are given for each of these results.

Periodic breathing induced by arterial oxygen partial pressure oscillations.

The Grodins model of respiratory control (Grodins et al., 1967) describes cardio-respiratory control for a lung with homogeneous gas concentrations. In this study we modify the Grodins model to take account of the inhomogeneities in gas concentration within the lung that are seen in many subjects with respiratory illnesses. This modification has the effect of lowering arterial oxygen partial pressure significantly. We investigate the effect on cardio-respiratory control of this low arterial oxygen signal and find that the governing equations may be reduced to a single delay-differential equation. This reduced model is found to be a good approximation to the full model and gives predictions that are similar to reported clinical data.

Multiple travelling-wave solutions in a minimal model for cell motility.

Two-phase flow models have been used previously to model cell motility. In order to reduce the complexity inherent with describing the many physical processes, we formulate a minimal model. Here we demonstrate that even the simplest 1D, two-phase, poroviscous, reactive flow model displays various types of behaviour relevant to cell crawling. We present stability analyses that show that an asymmetric perturbation is required to cause a spatially uniform, stationary strip of cytoplasm to move, which is relevant to cell polarization. Our numerical simulations identify qualitatively distinct families of travelling-wave solutions that coexist at certain parameter values. Within each family, the crawling speed of the strip has a bell-shaped dependence on the adhesion strength. The model captures the experimentally observed behaviour that cells crawl quickest at intermediate adhesion strengths, when the substrate is neither too sticky nor too slippy.

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.

A mathematical model of electron transfer within the mitochondrial respiratory cytochromes.

A simple mathematical model of electron flow along the mitochondrial respiratory cytochrome assembly and the transfer of electrons to molecular oxygen is presented. First, an expression for the current-voltage relationship for a biological oxygen electrode is derived, and from this the relationship between oxygen consumption rate and oxygen partial pressure is determined. An independent relationship between mitochondrial oxygen partial pressure and oxygen supply rate is then derived. By eliminating oxygen partial pressure from these two expressions, we may obtain a relationship between oxygen supply rate and oxygen consumption rate. This model is then used to investigate the effects of tissue dysoxia, uncoupling of oxidative phosphorylation, increased cellular diffusional resistance and inhomogeneities in oxygen supply on oxygen consumption. It is concluded that each of the above contribute in varying degrees to the phenomenon of "pathological oxygen supply dependency".