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.

Optimal control of diffuser shapes for non-uniform flow

A simplified model is used to identify the diffuser shape that maximises pressure recovery for several classes of non-uniform inflow. We find that optimal diffuser shapes strike a balance between not widening too soon, as this accentuates the non-uniform flow, and not staying narrow for too long, which is detrimental for wall drag. Three classes of non-uniform inflow are considered, with the axial velocity varying across the width of the diffuser entrance. The first case has inner and outer streams of different speeds, with a velocity jump between them that evolves into a shear layer downstream. The second case is a limiting case when these streams are of similar speed. The third case is a pure shear profile with linear velocity variation between the centre and outer edge of the diffuser. We describe the evolution of the time-averaged flow profile using a reduced mathematical model that has been previously tested against experiments and computational fluid dynamics models. The model consists of integrated mass and momentum equations, where wall drag is treated with a friction factor parameterisation. The governing equations of this model form the dynamics of an optimal control problem where the control is the diffuser channel shape. A numerical optimisation approach is used to solve the optimal control problem and Pontryagin's maximum principle is used to find analytical solutions in the second and third cases. We show that some of the optimal diffuser shapes can be well approximated by piecewise linear sections. This suggests a low-dimensional parameterisation of the shapes, providing a structure in which more detailed and computationally expensive turbulence models can be used to find optimal shapes for more realistic flow behaviour.

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.

Stochastic modelling of membrane filtration

Membrane fouling during particle filtration occurs through a variety of mechanisms, including internal pore clogging by contaminants, coverage of pore entrances and deposition on the membrane surface. In this paper, we present an efficient method for modelling the behaviour of a filter, which accounts for different retention mechanisms, particle sizes and membrane geometries. The membrane is assumed to be composed of a series of, possibly interconnected, pores. The central feature is a conductivity function, which describes the blockage of each individual pore as particles arrive, which is coupled with a mechanism to account for the stochastic nature of the arrival times of particles at the pore. The result is a system of ordinary differential equations based on the pore-level interactions. We demonstrate how our model can accurately describe a wide range of filtration scenarios. Specifically, we consider a case where blocking via multiple mechanisms can occur simultaneously, which have previously required the study through individual models; the filtration of a combination of small and large particles by a track-etched membrane and particle separation using interconnected pore networks. The model is significantly faster than comparable stochastic simulations for small networks, enabling its use as a tool for efficient future simulations.

Computational modelling of placental amino acid transfer as an integrated system

Placental amino acid transfer is essential for fetal development and its impairment is associated with poor fetal growth. Amino acid transfer is mediated by a broad array of specific plasma membrane transporters with overlapping substrate specificity. However, it is not fully understood how these different transporters work together to mediate net flux across the placenta. Therefore the aim of this study was to develop a new computational model to describe how human placental amino acid transfer functions as an integrated system. Amino acid transfer from mother to fetus requires transport across the two plasma membranes of the placental syncytiotrophoblast, each of which contains a distinct complement of transporter proteins. A compartmental modelling approach was combined with a carrier based modelling framework to represent the kinetics of the individual accumulative, exchange and facilitative classes of transporters on each plasma membrane. The model successfully captured the principal features of transplacental transfer. Modelling results clearly demonstrate how modulating transporter activity and conditions such as phenylketonuria, can increase the transfer of certain groups of amino acids, but that this comes at the cost of decreasing the transfer of others, which has implications for developing clinical treatment options in the placenta and other transporting epithelia.

•

First computational model of placental amino acid transfer as an integrated system

•

Increased activity of a transporter does not mean increased transfer to the fetus.

•

Increasing transfer of certain amino acids can reduce the transfer of others.

•

Amino acid composition as well as concentration determines transfer to the fetus.

•

Modelling of phenylketonuria suggests inhibition by excess maternal phenylalanine.

Mathematical modelling of tissue formation in chondrocyte filter cultures

In the field of cartilage tissue engineering, filter cultures are a frequently used three-dimensional differentiation model. However, understanding of the governing processes of in vitro growth and development of tissue in these models is limited. Therefore, this study aimed to further characterise these processes by means of an approach combining both experimental and applied mathematical methods. A mathematical model was constructed, consisting of partial differential equations predicting the distribution of cells and glycosaminoglycans (GAGs), as well as the overall thickness of the tissue. Experimental data was collected to allow comparison with the predictions of the simulation and refinement of the initial models. Healthy mature equine chondrocytes were expanded and subsequently seeded on collagen-coated filters and cultured for up to 7 weeks. Resulting samples were characterised biochemically, as well as histologically. The simulations showed a good representation of the experimentally obtained cell and matrix distribution within the cultures. The mathematical results indicate that the experimental GAG and cell distribution is critically dependent on the rate at which the cell differentiation process takes place, which has important implications for interpreting experimental results. This study demonstrates that large regions of the tissue are inactive in terms of proliferation and growth of the layer. In particular, this would imply that higher seeding densities will not significantly affect the growth rate. A simple mathematical model was developed to predict the observed experimental data and enable interpretation of the principal underlying mechanisms controlling growth-related changes in tissue composition.

Experimental-computational evaluation of human bone marrow stromal cell spreading on trabecular bone structures

The clinical application of macro-porous scaffolds for bone regeneration is significantly affected by the problem of insufficient cell colonization. Given the wide variety of different scaffold structures used for tissue engineering it is essential to derive relationships for cell colonization independent of scaffold architecture. To study cell population spreading on 3D structures decoupled from nutrient limitations, an in vitro culture system was developed consisting of thin slices of human trabecular bone seeded with Human Bone Marrow Stromal Cells, combined with dedicated microCT imaging and computational modeling of cell population spreading. Only the first phase of in vitro scaffold colonization was addressed, in which cells migrate and proliferate up to the stage when the surface of the bone is covered as a monolayer, a critical prerequisite for further tissue formation. The results confirm the model's ability to represent experimentally observed cell population spreading. The key advantage of the computational model was that by incorporating complex 3D structure, cell behavior can be characterized quantitatively in terms of intrinsic migration parameters, which could potentially be used for predictions on different macro-porous scaffolds subject to additional experimental validation. This type of modeling will prove useful in predicting cell colonization and improving strategies for skeletal tissue engineering.

Modelling the formation of necrotic regions in avascular tumours

The mechanisms underlying the formation of necrotic regions within avascular tumours are not well understood. In this paper, we examine the relative roles of nutrient deprivation and of cell death, from both the proliferating phase of the cell cycle via apoptosis and from the quiescent phase via necrosis, in changing the structure within multicellular tumour spheroids and particularly the accumulation of dead cell material in the centre. A mathematical model is presented and studied that accounts for nutrient diffusion, changes in cell cycling rates, the two different routes to cell death as well as active motion of cells and passive motion of the dead cell material. In studying the accumulation of dead cell matter we do not distinguish between the route by which each was formed. The resulting mathematical model is examined for a number of scenarios. Results show that in many cases the size of the necrotic core is closely correlated with low levels in nutrient concentration. However, in certain cases, particularly where the rate of necrosis is large, the resulting necrotic core can lead to regions of non-negligible nutrient concentration-dependent upon the mode of cell death.

Modelling the cell cycle and cell movement in multicellular tumour spheroids

This paper analyses a recent mathematical model of avascular tumour spheroid growth which accounts for both cell cycle dynamics and chemotactic driven cell movement. The model considers cells to exist in one of two compartments: proliferating and quiescent, as well as accounting for necrosis and apoptosis. One particular focus of this paper is the behaviour created when proliferating and quiescent cells have different chemotactic responses to an extracellular nutrient supply. Two very different steady-state behaviours are identified corresponding to those cases where proliferating cells move either more quickly or more slowly than quiescent cells in response to a gradient in the extracellular nutrient supply. The case where proliferating cells move more rapidly leads to the commonly accepted spheroid structure of a thin layer of proliferating cells surrounding an inner quiescent core. In the case where proliferating cells move more slowly than quiescent cells the model predicts an interesting structure of a thin layer of quiescent cells surrounding an inner core of proliferating and quiescent cells. The sensitivity of this tumour structure to the cell cycle model parameters is also discussed. In particular variations in the steady-state size of the tumour and the types of transient behaviour are explored. The model reveals interesting transient behaviour with sharply delineated regions of proliferating and quiescent cells.

Mechanisms to explain the reverse perivascular transport of solutes out of the brain

Experimental studies and observations in the human brain indicate that interstitial fluid and solutes, such as amyloid-beta (Abeta), are eliminated from grey matter of the brain along pericapillary and periarterial pathways. It is unclear, however, what constitutes the motive force for such transport within blood vessel walls, which is in the opposite direction to blood flow. In this paper the potential for global pressure differences to achieve such transport are considered. A mathematical model is constructed in order to test the hypothesis that perivascular drainage of interstitial fluid and solutes out of brain tissue is driven by pulsations of the blood vessel walls. Here it is assumed that drainage occurs through a thin layer between astrocytes and endothelial cells or between smooth muscle cells. The model suggests that, during each pulse cycle, there are periods when fluid and solutes are driven along perivascular spaces in the reverse direction to the flow of blood. It is shown that successful drainage may depend upon some attachment of solutes to the lining of the perivascular space, in order to produce a valve-like effect, although an alternative without this requirement is also postulated. Reduction in pulse amplitude, as in ageing cerebral vessels, would prolong the attachment time, encourage precipitation of Abeta peptides in vessel walls, and impair elimination of Abeta from the brain. These factors may play a role in the pathogenesis of cerebral amyloid angiopathy and in the accumulation of Abeta in the brain in Alzheimer's disease.

Dynamics of tear film deposition and draining

This paper investigates the deposition of the tear film on the cornea of the human eye. The tear film is laid down by the motion of the upper eyelid and then subsequently flows and thins. Of particular interest is the stability of the tear layer and the development of dry patches on the cornea. While there has been significant research on the behaviour of tear films between blinks, this paper focuses on understanding the mechanisms which control the shape and thickness of the deposited film and how this affects the subsequent film behaviour. Numerical and analytical methods are applied to a lubrication model which includes the effects of surface tension, viscosity, gravity and evaporation. The model reveals the importance of the eyelid velocity, motion of the surface lipid layer and the storage of tear film between blinks.

Mathematical modelling of skeletal repair

Tissue engineering offers significant promise as a viable alternative to current clinical strategies for replacement of damaged tissue as a consequence of disease or trauma. Since mathematical modelling is a valuable tool in the analysis of complex systems, appropriate use of mathematical models has tremendous potential for advancing the understanding of the physical processes involved in such tissue reconstruction. In this review, the potential benefits, and limitations, of theoretical modelling in tissue engineering applications are examined with specific emphasis on tissue engineering of bone. A central tissue engineering approach is the in vivo implantation of a biomimetic scaffold seeded with an appropriate population of stem or progenitor cells. This review will therefore consider the theory behind a number of key factors affecting the success of such a strategy including: stem cell or progenitor population expansion and differentiation ex vivo; cell adhesion and migration, and the effective design of scaffolds; and delivery of nutrient to avascular structures. The focus will be on current work in this area, as well as on highlighting limitations and suggesting possible directions for future work to advance health-care for all.

Tumour dynamics and necrosis: surface tension and stability

A model is developed for the motion of cells within a multicell spherical tumour. The model allows either for the intercellular forces to be in compression and cells to be compacted to a fixed number density, or for the cell number density to fall and cells to become isolated from each other. The model develops necrotic regions naturally due to force balances rather than being directly attributable to a critical oxygen concentration. These necrotic regions may result in a gradual reduction in local cell density rather than jump to a completely dead region. Numerical and analytical analysis of the spherically symmetric model shows that the long time behaviour of the spheroid depends on any surface tension effects created by cells on the outer surface. For small surface tension the spheroid grows linearly in time developing a large necrotic region, while for larger surface tension the growth can be halted. The linear stability to spherically symmetric perturbations of all the possible resulting steady states is revealed.

The migration of cells in multicell tumor spheroids

A mathematical model is proposed to explain the observed internalization of microspheres and 3H-thymidine labelled cells in steady-state multicellular spheroids. The model uses the conventional ideas of nutrient diffusion and consumption by the cells. In addition, a very simple model of the progress of the cells through the cell cycle is considered. Cells are divided into two classes, those proliferating (being in G1, S, G2 or M phases) and those that are quiescent (being in G0). Furthermore, the two categories are presumed to have different chemotactic responses to the nutrient gradient. The model accounts for the spatial and temporal variations in the cell categories together with mitosis, conversion between categories and cell death. Numerical solutions demonstrate that the model predicts the behavior similar to existing models but has some novel effects. It allows for spheroids to approach a steady-state size in a non-monotonic manner, it predicts self-sorting of the cell classes to produce a thin layer of rapidly proliferating cells near the outer surface and significant numbers of cells within the spheroid stalled in a proliferating state. The model predicts that overall tumor growth is not only determined by proliferation rates but also by the ability of cells to convert readily between the classes. Moreover, the steady-state structure of the spheroid indicates that if the outer layers are removed then the tumor grows quickly by recruiting cells stalled in a proliferating state. Questions are raised about the chemotactic response of cells in differing phases and to the dependency of cell cycle rates to nutrient levels.

Traveling-wave relaxation in elongated liquid crystal cells

We have made a theoretical study of Freedericksz relaxation in a long thin nematic liquid crystal cell subject to strong anchoring on the short ends and weak anchoring on the long sides. On removing an imposed magnetic field, three different types of relaxation behavior may be observed. Two of these are simple generalizations of one-dimensional relaxation channels, and are dominated by either the ends or the sides. The third is a traveling wave, nucleated by the strong anchoring ends of the cell but driven by the weak anchoring sides and is the result of a subtle balance between the two classical mechanisms. A phase diagram is derived, identifying the relaxation regimes as a function of the nondimensional initial field and the anchoring strength in the long cell limit. A comparison is made between numerical results and a simple one-dimensional theory derived from an asymptotic analysis. Surprisingly, the traveling wave behavior occurs for a large region of parameter space.