Cascading failures in networks of heterogeneous node behavior.

Variability in the dynamical function of nodes comprising a complex network impacts upon cascading failures that can compromise the network's ability to operate. Node types correspond to sources, sinks, or passive conduits of a current flow, applicable to renewable electrical power microgrids containing a variable number of intermittently operating generators and consumers of power. The resilience to cascading failures of ensembles of synthetic networks with different topology is examined as a function of the edge current carrying capacity and mix of node types, together with exemplar real-world networks. While a network with a homogeneous composition of node types can be resilient to failure, onewith an identical topology but with heterogeneous nodes can be strongly susceptible to failure. For networks with similar numbers of sources, sinks, and passive nodes the mean resilience decreases as networks become more disordered. Nevertheless all network topologies have enhanced regions of resilience, accessible by the manipulation of node composition and functionality.

Cellular Uptake and Efflux of Palbociclib In Vitro in Single Cell and Spheroid Models.

Adequate drug distribution through tumors is essential for treatment to be effective. Palbociclib is a cyclin-dependent kinase 4/6 inhibitor approved for use in patients with hormone receptor positive, human epidermal growth factor receptor 2 negative metastatic breast cancer. It has unusual physicochemical properties, which may significantly influence its distribution in tumor tissue. We studied the penetration and distribution of palbociclib in vitro, including the use of multicellular three-dimensional models and mathematical modeling. MCF-7 and DLD-1 cell lines were grown as single cell suspensions (SCS) and spheroids; palbociclib uptake and efflux were studied using liquid chromatography-tandem mass spectrometry. Intracellular concentrations of palbociclib for MCF-7 SCS (C max 3.22 µM) and spheroids (C max 2.91 µM) were 32- and 29-fold higher and in DLD-1, 13- and 7-fold higher, respectively, than the media concentration (0.1 µM). Total palbociclib uptake was lower in DLD-1 cells than MCF-7 cells in both SCS and spheroids. Both uptake and efflux of palbociclib were slower in spheroids than SCS. These data were used to develop a mathematical model of palbociclib transport that quantifies key parameters determining drug penetration and distribution. The model reproduced qualitatively most features of the experimental data and distinguished between SCS and spheroids, providing additional support for hypotheses derived from the experimental data. Mathematical modeling has the potential for translating in vitro data into clinically relevant estimates of tumor drug concentrations. SIGNIFICANCE STATEMENT: This study explores palbociclib uptake and efflux in single cell suspension and spheroid models of cancer. Large intracellular concentrations of palbociclib are found after drug exposure. The data from this study may aid understanding of the intratumoural pharmacokinetics of palbociclib, which is useful in understanding how drug distributes within tumor tissue and optimizing drug efficacy. Biomathematical modelling has the potential to derive intratumoural drug concentrations from plasma pharmacokinetics in patients.

Effective equations governing an active poroelastic medium.

In this work, we consider the spatial homogenization of a coupled transport and fluid-structure interaction model, to the end of deriving a system of effective equations describing the flow, elastic deformation and transport in an active poroelastic medium. The 'active' nature of the material results from a morphoelastic response to a chemical stimulant, in which the growth time scale is strongly separated from other elastic time scales. The resulting effective model is broadly relevant to the study of biological tissue growth, geophysical flows (e.g. swelling in coals and clays) and a wide range of industrial applications (e.g. absorbant hygiene products). The key contribution of this work is the derivation of a system of homogenized partial differential equations describing macroscale growth, coupled to transport of solute, that explicitly incorporates details of the structure and dynamics of the microscopic system, and, moreover, admits finite growth and deformation at the pore scale. The resulting macroscale model comprises a Biot-type system, augmented with additional terms pertaining to growth, coupled to an advection-reaction-diffusion equation. The resultant system of effective equations is then compared with other recent models under a selection of appropriate simplifying asymptotic limits.

A multiscale analysis of nutrient transport and biological tissue growth in vitro.

In this paper, we consider the derivation of macroscopic equations appropriate to describe the growth of biological tissue, employing a multiple-scale homogenization method to accommodate explicitly the influence of the underlying microscale structure of the material, and its evolution, on the macroscale dynamics. Such methods have been widely used to study porous and poroelastic materials; however, a distinguishing feature of biological tissue is its ability to remodel continuously in response to local environmental cues. Here, we present the derivation of a model broadly applicable to tissue engineering applications, characterized by cell proliferation and extracellular matrix deposition in porous scaffolds used within tissue culture systems, which we use to study coupling between fluid flow, nutrient transport, and microscale tissue growth. Attention is restricted to surface accretion within a rigid porous medium saturated with a Newtonian fluid; coupling between the various dynamics is achieved by specifying the rate of microscale growth to be dependent upon the uptake of a generic diffusible nutrient. The resulting macroscale model comprises a Darcy-type equation governing fluid flow, with flow characteristics dictated by the assumed periodic microstructure and surface growth rate of the porous medium, coupled to an advection-reaction equation specifying the nutrient concentration. Illustrative numerical simulations are presented to indicate the influence of microscale growth on macroscale dynamics, and to highlight the importance of including experimentally relevant microstructural information to correctly determine flow dynamics and nutrient delivery in tissue engineering applications.

The interplay between tissue growth and scaffold degradation in engineered tissue constructs.

In vitro tissue engineering is emerging as a potential tool to meet the high demand for replacement tissue, caused by the increased incidence of tissue degeneration and damage. A key challenge in this field is ensuring that the mechanical properties of the engineered tissue are appropriate for the in vivo environment. Achieving this goal will require detailed understanding of the interplay between cell proliferation, extracellular matrix (ECM) deposition and scaffold degradation. In this paper, we use a mathematical model (based upon a multiphase continuum framework) to investigate the interplay between tissue growth and scaffold degradation during tissue construct evolution in vitro. Our model accommodates a cell population and culture medium, modelled as viscous fluids, together with a porous scaffold and ECM deposited by the cells, represented as rigid porous materials. We focus on tissue growth within a perfusion bioreactor system, and investigate how the predicted tissue composition is altered under the influence of (1) differential interactions between cells and the supporting scaffold and their associated ECM, (2) scaffold degradation, and (3) mechanotransduction-regulated cell proliferation and ECM deposition. Numerical simulation of the model equations reveals that scaffold heterogeneity typical of that obtained from [Formula: see text]CT scans of tissue engineering scaffolds can lead to significant variation in the flow-induced mechanical stimuli experienced by cells seeded in the scaffold. This leads to strong heterogeneity in the deposition of ECM. Furthermore, preferential adherence of cells to the ECM in favour of the artificial scaffold appears to have no significant influence on the eventual construct composition; adherence of cells to these supporting structures does, however, lead to cell and ECM distributions which mimic and exaggerate the heterogeneity of the underlying scaffold. Such phenomena have important ramifications for the mechanical integrity of engineered tissue constructs and their suitability for implantation in vivo.

The isolation of spatial patterning modes in a mathematical model of juxtacrine cell signalling.

Juxtacrine signalling mechanisms are known to be crucial in tissue and organ development, leading to spatial patterns in gene expression. We investigate the patterning behaviour of a discrete model of juxtacrine cell signalling due to Owen & Sherratt (1998, Mathematical modelling of juxtacrine cell signalling. Math. Biosci., 153, 125-150) in which ligand molecules, unoccupied receptors and bound ligand-receptor complexes are modelled. Feedback between the ligand and receptor production and the level of bound receptors is incorporated. By isolating two parameters associated with the feedback strength and employing numerical simulation, linear stability and bifurcation analysis, the pattern-forming behaviour of the model is analysed under regimes corresponding to lateral inhibition and induction. Linear analysis of this model fails to capture the patterning behaviour exhibited in numerical simulations. Via bifurcation analysis, we show that since the majority of periodic patterns fold subcritically from the homogeneous steady state, a wide variety of stable patterns exists at a given parameter set, providing an explanation for this failure. The dominant pattern is isolated via numerical simulation. Additionally, by sampling patterns of non-integer wavelength on a discrete mesh, we highlight a disparity between the continuous and discrete representations of signalling mechanisms: in the continuous case, patterns of arbitrary wavelength are possible, while sampling such patterns on a discrete mesh leads to longer wavelength harmonics being selected where the wavelength is rational; in the irrational case, the resulting aperiodic patterns exhibit 'local periodicity', being constructed from distorted stable shorter wavelength patterns. This feature is consistent with experimentally observed patterns, which typically display approximate short-range periodicity with defects.

Continuum limits of pattern formation in hexagonal-cell monolayers.

Intercellular signalling is key in determining cell fate. In closely packed tissues such as epithelia, juxtacrine signalling is thought to be a mechanism for the generation of fine-grained spatial patterns in cell differentiation commonly observed in early development. Theoretical studies of such signalling processes have shown that negative feedback between receptor activation and ligand production is a robust mechanism for fine-grained pattern generation and that cell shape is an important factor in the resulting pattern type. It has previously been assumed that such patterns can be analysed only with discrete models since significant variation occurs over a lengthscale concomitant with an individual cell; however, considering a generic juxtacrine signalling model in square cells, in O'Dea and King (Math Biosci 231(2):172-185 2011), a systematic method for the derivation of a continuum model capturing such phenomena due to variations in a model parameter associated with signalling feedback strength was presented. Here, we extend this work to derive continuum models of the more complex fine-grained patterning in hexagonal cells, constructing individual models for the generation of patterns from the homogeneous state and for the transition between patterning modes. In addition, by considering patterning behaviour under the influence of simultaneous variation of feedback parameters, we construct a more general continuum representation, capturing the emergence of the patterning bifurcation structure. Comparison with the steady-state and dynamic behaviour of the underlying discrete system is made; in particular, we consider pattern-generating travelling waves and the competition between various stable patterning modes, through which we highlight an important deficiency in the ability of continuum representations to accommodate certain dynamics associated with discrete systems.

Multiscale analysis of pattern formation via intercellular signalling.

Lateral inhibition, a juxtacrine signalling mechanism by which a cell adopting a particular fate inhibits neighbouring cells from doing likewise, has been shown to be a robust mechanism for the formation of fine-grained spatial patterns (in which adjacent cells in developing tissues diverge to achieve contrasting states of differentiation), provided that there is sufficiently strong feedback. The fine-grained nature of these patterns poses problems for analysis via traditional continuum methods since these require that significant variation takes place only over lengthscales much larger than an individual cell and such systems have therefore been investigated primarily using discrete methods. Here, however, we apply a multiscale method to derive systematically a continuum model from the discrete Delta-Notch signalling model of Collier et al. (J.R. Collier, N.A.M. Monk, P.K. Maini, J.H. Lewis, Pattern formation by lateral inhibition with feedback: a mathematical model of Delta-Notch intercellular signalling, J. Theor. Biol., 183, 1996, 429-446) under particular assumptions on the parameters, which we use to analyse the generation of fine-grained patterns. We show that, on the macroscale, the contact-dependent juxtacrine signalling interaction manifests itself as linear diffusion, motivating the use of reaction-diffusion-based models for such cell-signalling systems. We also analyse the travelling-wave behaviour of our system, obtaining good quantitative agreement with the discrete system.

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.

A multiphase model for tissue construct growth in a perfusion bioreactor.

The growth of a cell population within a rigid porous scaffold in a perfusion bioreactor is studied, using a three-phase continuum model of the type presented by Lemon et al. (2006, Multiphase modelling of tissue growth using the theory of mixtures. J. Math. Biol., 52, 571-594) to represent the cell population (and attendant extracellular matrix), culture medium and porous scaffold. The bioreactor system is modelled as a 2D channel containing the cell-seeded rigid porous scaffold (tissue construct) which is perfused with culture medium. The study concentrates on (i) the cell-cell and cell-scaffold interactions and (ii) the impact of mechanotransduction mechanisms on construct composition. A numerical and analytical analysis of the model equations is presented and, depending upon the relative importance of cell aggregation and repulsion, markedly different cell movement is revealed. Additionally, mechanotransduction effects due to cell density, pressure and shear stress-mediated tissue growth are shown to generate qualitative differences in the composition of the resulting construct. The results of our simulations indicate that this model formulation (in conjunction with appropriate experimental data) has the potential to provide a means of identifying the dominant regulatory stimuli in a cell population.