Surfactant amplifies yield-stress effects in the capillary instability of a film coating a tube.

To assess how the presence of surfactant in lung airways alters the flow of mucus that leads to plug formation and airway closure, we investigate the effect of insoluble surfactant on the instability of a viscoplastic liquid coating the interior of a cylindrical tube. Evolution equations for the layer thickness using thin-film and long-wave approximations are derived that incorporate yield-stress effects and capillary and Marangoni forces. Using numerical simulations and asymptotic analysis of the thin-film system, we quantify how the presence of surfactant slows growth of the Rayleigh-Plateau instability, increases the size of initial perturbation required to trigger instability and decreases the final peak height of the layer. When the surfactant strength is large, the thin-film dynamics coincide with the dynamics of a surfactant-free layer but with time slowed by a factor of four and the capillary Bingham number, a parameter proportional to the yield stress, exactly doubled. By solving the long-wave equations numerically, we quantify how increasing surfactant strength can increase the critical layer thickness for plug formation to occur and delay plugging. The previously established effect of the yield stress in suppressing plug formation [Shemilt et al., J. Fluid Mech., 2022, vol. 944, A22] is shown to be amplified by introducing surfactant. We discuss the implications of these results for understanding the impact of surfactant deficiency and increased mucus yield stress in obstructive lung diseases.

Modeling collagen fibril self-assembly from extracellular medium in embryonic tendon.

Collagen is a key structural component of multicellular organisms and is arranged in a highly organized manner. In structural tissues such as tendons, collagen forms bundles of parallel fibers between cells, which appear within a 24-h window between embryonic day 13.5 (E13.5) and E14.5 during mouse embryonic development. Current models assume that the organized structure of collagen requires direct cellular control, whereby cells actively lay down collagen fibrils from cell surfaces. However, such models appear incompatible with the time and length scales of fibril formation. We propose a phase-transition model to account for the rapid development of ordered fibrils in embryonic tendon, reducing reliance on active cellular processes. We develop phase-field crystal simulations of collagen fibrillogenesis in domains derived from electron micrographs of inter-cellular spaces in embryonic tendon and compare results qualitatively and quantitatively to observed patterns of fibril formation. To test the prediction of this phase-transition model that free protomeric collagen should exist in the inter-cellular spaces before the formation of observable fibrils, we use laser-capture microdissection, coupled with mass spectrometry, which demonstrates steadily increasing free collagen in inter-cellular spaces up to E13.5, followed by a rapid reduction of free collagen that coincides with the appearance of less-soluble collagen fibrils. The model and measurements together provide evidence for extracellular self-assembly of collagen fibrils in embryonic mouse tendon, supporting an additional mechanism for rapid collagen fibril formation during embryonic development.

Couple stresses and discrete potentials in the vertex model of cellular monolayers.

The vertex model is widely used to simulate the mechanical properties of confluent epithelia and other multicellular tissues. This inherently discrete framework allows a Cauchy stress to be attributed to each cell, and its symmetric component has been widely reported, at least for planar monolayers. Here, we consider the stress attributed to the neighbourhood of each tricellular junction, evaluating in particular its leading-order antisymmetric component and the associated couple stresses, which characterise the degree to which individual cells experience (and resist) in-plane bending deformations. We develop discrete potential theory for localised monolayers having disordered internal structure and use this to derive the analogues of Airy and Mindlin stress functions. These scalar potentials typically have broad-banded spectra, highlighting the contributions of small-scale defects and boundary layers to global stress patterns. An affine approximation attributes couple stresses to pressure differences between cells sharing a trijunction, but simulations indicate an additional role for non-affine deformations.

Red blood cell dynamics in extravascular biological tissues modelled as canonical disordered porous media.

The dynamics of blood flow in the smallest vessels and passages of the human body, where the cellular character of blood becomes prominent, plays a dominant role in the transport and exchange of solutes. Recent studies have revealed that the microhaemodynamics of a vascular network is underpinned by its interconnected structure, and certain structural alterations such as capillary dilation and blockage can substantially change blood flow patterns. However, for extravascular media with disordered microstructure (e.g. the porous intervillous space in the placenta), it remains unclear how the medium's structure affects the haemodynamics. Here, we simulate cellular blood flow in simple models of canonical porous media representative of extravascular biological tissue, with corroborative microfluidic experiments performed for validation purposes. For the media considered here, we observe three main effects: first, the relative apparent viscosity of blood increases with the structural disorder of the medium; second, the presence of red blood cells (RBCs) dynamically alters the flow distribution in the medium; third, symmetry breaking introduced by moderate structural disorder can promote more homogeneous distribution of RBCs. Our findings contribute to a better understanding of the cell-scale haemodynamics that mediates the relationship linking the function of certain biological tissues to their microstructure.

Advection-dominated transport past isolated disordered sinks: stepping beyond homogenization.

We investigate the transport of a solute past isolated sinks in a bounded domain when advection is dominant over diffusion, evaluating the effectiveness of homogenization approximations when sinks are distributed uniformly randomly in space. Corrections to such approximations can be non-local, non-smooth and non-Gaussian, depending on the physical parameters (a Péclet number Pe, assumed large, and a Damköhler number Da) and the compactness of the sinks. In one spatial dimension, solute distributions develop a staircase structure for large Pe , with corrections being better described with credible intervals than with traditional moments. In two and three dimensions, solute distributions are near-singular at each sink (and regularized by sink size), but their moments can be smooth as a result of ensemble averaging over variable sink locations. We approximate corrections to a homogenization approximation using a moment-expansion method, replacing the Green's function by its free-space form, and test predictions against simulation. We show how, in two or three dimensions, the leading-order impact of disorder can be captured in a homogenization approximation for the ensemble mean concentration through a modification to Da that grows with diminishing sink size.

Model-based Bayesian inference of the ventilation distribution in patients with cystic fibrosis from multiple breath washout, with comparison to ventilation MRI.

BACKGROUND: Indices of ventilation heterogeneity (VH) from multiple breath washout (MBW) have been shown to correlate well with VH indices derived from hyperpolarised gas ventilation MRI. Here we report the prediction of ventilation distributions from MBW data using a mathematical model, and the comparison of these predictions with imaging data. METHODS: We developed computer simulations of the ventilation distribution in the lungs to model MBW measurement with 3 parameters: σV, determining the extent of VH; V0, the lung volume; and VD, the dead-space volume. These were inferred for each individual from supine MBW data recorded from 25 patients with cystic fibrosis (CF) using approximate Bayesian computation. The fitted models were used to predict the distribution of gas imaged by 3He ventilation MRI measurements collected from the same visit. RESULTS: The MRI indices measured (I1/3, the fraction of pixels below one-third of the mean intensity and ICV, the coefficient of variation of pixel intensity) correlated strongly with those predicted by the MBW model fits (r=0.93,0.88 respectively). There was also good agreement between predicted and measured MRI indices (mean bias ± limits of agreement: I1/3:-0.003±0.118 and ICV:-0.004±0.298). Fitted model parameters were robust to truncation of MBW data. CONCLUSION: We have shown that the ventilation distribution in the lung can be inferred from an MBW signal, and verified this using ventilation MRI. The Bayesian method employed extracts this information with fewer breath cycles than required for LCI, reducing acquisition time required, and gives uncertainty bounds, which are important for clinical decision making.

Collagen fibril assembly: New approaches to unanswered questions.

Collagen fibrils are essential for metazoan life. They are the largest, most abundant, and most versatile protein polymers in animals, where they occur in the extracellular matrix to form the structural basis of tissues and organs. Collagen fibrils were first observed at the turn of the 20th century. During the last 40 years, the genes that encode the family of collagens have been identified, the structure of the collagen triple helix has been solved, the many enzymes involved in the post-translational modifications of collagens have been identified, mutations in the genes encoding collagen and collagen-associated proteins have been linked to heritable disorders, and changes in collagen levels have been associated with a wide range of diseases, including cancer. Yet despite extensive research, a full understanding of how cells assemble collagen fibrils remains elusive. Here, we review current models of collagen fibril self-assembly, and how cells might exert control over the self-assembly process to define the number, length and organisation of fibrils in tissues.

Force networks, torque balance and Airy stress in the planar vertex model of a confluent epithelium.

The vertex model is a popular framework for modelling tightly packed biological cells, such as confluent epithelia. Cells are described by convex polygons tiling the plane and their equilibrium is found by minimizing a global mechanical energy, with vertex locations treated as degrees of freedom. Drawing on analogies with granular materials, we describe the force network for a localized monolayer and derive the corresponding discrete Airy stress function, expressed for each N-sided cell as N scalars defined over kites covering the cell. We show how a torque balance (commonly overlooked in implementations of the vertex model) requires each internal vertex to lie at the orthocentre of the triangle formed by neighbouring edge centroids. Torque balance also places a geometric constraint on the stress in the neighbourhood of cellular trijunctions, and requires cell edges to be orthogonal to the links of a dual network that connect neighbouring cell centres and thereby triangulate the monolayer. We show how the Airy stress function depends on cell shape when a standard energy functional is adopted, and discuss implications for computational implementations of the model.

Quantifying the impact of tissue metabolism on solute transport in feto-placental microvascular networks.

The primary exchange units in the human placenta are terminal villi, in which fetal capillary networks are surrounded by a thin layer of villous tissue, separating fetal from maternal blood. To understand how the complex spatial structure of villi influences their function, we use an image-based theoretical model to study the effect of tissue metabolism on the transport of solutes from maternal blood into the fetal circulation. For solute that is taken up under first-order kinetics, we show that the transition between flow-limited and diffusion-limited transport depends on two new dimensionless parameters defined in terms of key geometric quantities, with strong solute uptake promoting flow-limited transport conditions. We present a simple algebraic approximation for solute uptake rate as a function of flow conditions, metabolic rate and villous geometry. For oxygen, accounting for nonlinear kinetics using physiological parameter values, our model predicts that villous metabolism does not significantly impact oxygen transfer to fetal blood, although the partitioning of fluxes between the villous tissue and the capillary network depends strongly on the flow regime.

Physical and geometric determinants of transport in fetoplacental microvascular networks.

Across mammalian species, solute exchange takes place in complex microvascular networks. In the human placenta, the primary exchange units are terminal villi that contain disordered networks of fetal capillaries and are surrounded externally by maternal blood. We show how the irregular internal structure of a terminal villus determines its exchange capacity for diverse solutes. Distilling geometric features into three parameters, obtained from image analysis and computational fluid dynamics, we capture archetypal features of the structure-function relationship of terminal villi using a simple algebraic approximation, revealing transitions between flow- and diffusion-limited transport at vessel and network levels. Our theory accommodates countercurrent effects, incorporates nonlinear blood rheology, and offers an efficient method for testing network robustness. Our results show how physical estimates of solute transport, based on carefully defined geometrical statistics, provide a viable method for linking placental structure and function and offer a framework for assessing transport in other microvascular systems.

Decoupling the Roles of Cell Shape and Mechanical Stress in Orienting and Cueing Epithelial Mitosis.

Distinct mechanisms involving cell shape and mechanical force are known to influence the rate and orientation of division in cultured cells. However, uncoupling the impact of shape and force in tissues remains challenging. Combining stretching of Xenopus tissue with mathematical methods of inferring relative mechanical stress, we find separate roles for cell shape and mechanical stress in orienting and cueing division. We demonstrate that division orientation is best predicted by an axis of cell shape defined by the position of tricellular junctions (TCJs), which align with local cell stress rather than tissue-level stress. The alignment of division to cell shape requires functional cadherin and the localization of the spindle orientation protein, LGN, to TCJs but is not sensitive to relative cell stress magnitude. In contrast, proliferation rate is more directly regulated by mechanical stress, being correlated with relative isotropic stress and decoupled from cell shape when myosin II is depleted.

Blood flow and transport in the human placenta.

The placenta is a multi-functional organ that exchanges blood gases and nutrients between a mother and her developing fetus. In humans, fetal blood flows through intricate networks of vessels confined within villous trees, the branches of which are bathed in pools of maternal blood. Fluid mechanics and transport processes play a central role in understanding how these elaborate structures contribute to the function of the placenta, and how their disorganization may lead to disease. Recent advances in imaging and computation have spurred significant advances in simulations of fetal and maternal flows within the placenta, across a range of lengthscales. Models describe jets of maternal blood emerging from spiral arteries into a disordered and deformable porous medium, and solute uptake by fetal blood flowing through elaborate three-dimensional capillary networks. We survey recent developments and emerging challenges in modeling flow and transport in this complex organ.

Modelling structural determinants of ventilation heterogeneity: A perturbative approach.

We have developed a computational model of gas mixing and ventilation in the human lung represented as a bifurcating network. We have simulated multiple-breath washout (MBW), a clinical test for measuring ventilation heterogeneity (VH) in patients with obstructive lung conditions. By applying airway constrictions inter-regionally, we have predicted the response of MBW indices to obstructions and found that they detect a narrow range of severe constrictions that reduce airway radius to 10%-30% of healthy values. These results help to explain the success of the MBW test to distinguish obstructive lung conditions from healthy controls. Further, we have used a perturbative approach to account for intra-regional airway heterogeneity that avoids modelling each airway individually. We have found, for random airway heterogeneity, that the variance in MBW indices is greater when indices are already elevated due to constrictions. By quantifying this effect, we have shown that variability in lung structure and mechanical properties alone can lead to clinically significant variability in MBW indices (specifically the Lung Clearance Index-LCI, and the gradient of phase-III slopes-Scond), but only in cases simulating obstructive lung conditions. This method is a computationally efficient way to probe the lung's sensitivity to structural changes, and to quantify uncertainty in predictions due to random variations in lung mechanical and structural properties.

A Mechanosensitive RhoA Pathway that Protects Epithelia against Acute Tensile Stress.

Adherens junctions are tensile structures that couple epithelial cells together. Junctional tension can arise from cell-intrinsic application of contractility or from the cell-extrinsic forces of tissue movement. Here, we report a mechanosensitive signaling pathway that activates RhoA at adherens junctions to preserve epithelial integrity in response to acute tensile stress. We identify Myosin VI as the force sensor, whose association with E-cadherin is enhanced when junctional tension is increased by mechanical monolayer stress. Myosin VI promotes recruitment of the heterotrimeric Gα12 protein to E-cadherin, where it signals for p114 RhoGEF to activate RhoA. Despite its potential to stimulate junctional actomyosin and further increase contractility, tension-activated RhoA signaling is necessary to preserve epithelial integrity. This is explained by an increase in tensile strength, especially at the multicellular vertices of junctions, that is due to mDia1-mediated actin assembly.

Effect of intermittent inspiratory leaks on measurement of lung clearance index using nitrogen and sulfur hexafluoride.

Even small air leaks can have a significant impact on LCI measured using N2. This is particularly the case for leaks that occur towards the end of washout. In contrast, leaks generally have a much smaller impact on LCI measured by SF6. http://ow.ly/az7b30lG5Ku.

Curvature-Sensitive Kinesin Binding Can Explain Microtubule Ring Formation and Reveals Chaotic Dynamics in a Mathematical Model.

Microtubules are filamentous tubular protein polymers which are essential for a range of cellular behaviour, and are generally straight over micron length scales. However, in some gliding assays, where microtubules move over a carpet of molecular motors, individual microtubules can also form tight arcs or rings, even in the absence of crosslinking proteins. Understanding this phenomenon may provide important explanations for similar highly curved microtubules which can be found in nerve cells undergoing neurodegeneration. We propose a model for gliding assays where the kinesins moving the microtubules over the surface induce ring formation through differential binding, substantiated by recent findings that a mutant version of the motor protein kinesin applied in solution is able to lock-in microtubule curvature. For certain parameter regimes, our model predicts that both straight and curved microtubules can exist simultaneously as stable steady states, as has been seen experimentally. Additionally, unsteady solutions are found, where a wave of differential binding propagates down the microtubule as it glides across the surface, which can lead to chaotic motion. Whilst this model explains two-dimensional microtubule behaviour in an experimental gliding assay, it has the potential to be adapted to explain pathological curling in nerve cells.

Osmotic and electroosmotic fluid transport across the retinal pigment epithelium: A mathematical model.

The retinal pigment epithelium (RPE) is the outermost cell layer of the retina. It has several important physiological functions, among which is removal of excess fluid from the sub-retinal space by pumping it isotonically towards the choroid. Failure of this pumping leads to fluid accumulation, which is closely associated with several pathological conditions, such as age-related macular degeneration, macular oedema and retinal detachment. In the present work we study mechanisms responsible for fluid transport across the RPE with the aim of understanding how fluid accumulation can be prevented. We focus on two possible mechanisms, osmosis and electroosmosis, and develop a spatially resolved mathematical model that couples fluid and ion transport across the epithelium, accounting for the presence of Na+,K+ and Cl- ions. Our model predicts spatial variability of ion concentrations and the electrical potential along the cleft gap between two adjacent cells, which osmotically drives the flow across the lateral membranes. This flow is directed from the sub-retinal space to the choroid and has a magnitude close to measured values. Electroosmosis is subdominant by three orders of magnitude to osmosis and has an opposite direction, suggesting that local osmosis is the main driving mechanism for water transport across the RPE.

Mechanical characterization of disordered and anisotropic cellular monolayers.

We consider a cellular monolayer, described using a vertex-based model, for which cells form a spatially disordered array of convex polygons that tile the plane. Equilibrium cell configurations are assumed to minimize a global energy defined in terms of cell areas and perimeters; energy is dissipated via dynamic area and length changes, as well as cell neighbor exchanges. The model captures our observations of an epithelium from a Xenopus embryo showing that uniaxial stretching induces spatial ordering, with cells under net tension (compression) tending to align with (against) the direction of stretch, but with the stress remaining heterogeneous at the single-cell level. We use the vertex model to derive the linearized relation between tissue-level stress, strain, and strain rate about a deformed base state, which can be used to characterize the tissue's anisotropic mechanical properties; expressions for viscoelastic tissue moduli are given as direct sums over cells. When the base state is isotropic, the model predicts that tissue properties can be tuned to a regime with high elastic shear resistance but low resistance to area changes, or vice versa.

Three-dimensional plant architecture and sunlit-shaded patterns: a stochastic model of light dynamics in canopies.

Background and Aims: Diurnal changes in solar position and intensity combined with the structural complexity of plant architecture result in highly variable and dynamic light patterns within the plant canopy. This affects productivity through the complex ways that photosynthesis responds to changes in light intensity. Current methods to characterize light dynamics, such as ray-tracing, are able to produce data with excellent spatio-temporal resolution but are computationally intensive and the resulting data are complex and high-dimensional. This necessitates development of more economical models for summarizing the data and for simulating realistic light patterns over the course of a day. Methods: High-resolution reconstructions of field-grown plants are assembled in various configurations to form canopies, and a forward ray-tracing algorithm is applied to the canopies to compute light dynamics at high (1 min) temporal resolution. From the ray-tracer output, the sunlit or shaded state for each patch on the plants is determined, and these data are used to develop a novel stochastic model for the sunlit-shaded patterns. The model is designed to be straightforward to fit to data using maximum likelihood estimation, and fast to simulate from. Key Results: For a wide range of contrasting 3-D canopies, the stochastic model is able to summarize, and replicate in simulations, key features of the light dynamics. When light patterns simulated from the stochastic model are used as input to a model of photoinhibition, the predicted reduction in carbon gain is similar to that from calculations based on the (extremely costly) ray-tracer data. Conclusions: The model provides a way to summarize highly complex data in a small number of parameters, and a cost-effective way to simulate realistic light patterns. Simulations from the model will be particularly useful for feeding into larger-scale photosynthesis models for calculating how light dynamics affects the photosynthetic productivity of canopies.