A dynamical model of TGF-ß activation in asthmatic airways.

Excessive activation of the regulatory cytokine transforming growth factor $\beta $ (TGF-$\beta $) via contraction of airway smooth muscle (ASM) is associated with the development of asthma. In this study, we develop an ordinary differential equation model that describes the change in density of the key airway wall constituents, ASM and extracellular matrix (ECM), and their interplay with subcellular signalling pathways leading to the activation of TGF-$\beta $. We identify bistable parameter regimes where there are two positive steady states, corresponding to either reduced or elevated TGF-$\beta $ concentration, with the latter leading additionally to increased ASM and ECM density. We associate the former with a healthy homeostatic state and the latter with a diseased (asthmatic) state. We demonstrate that external stimuli, inducing TGF-$\beta $ activation via ASM contraction (mimicking an asthmatic exacerbation), can perturb the system irreversibly from the healthy state to the diseased one. We show that the properties of the stimuli, such as their frequency or strength, and the clearance of surplus active TGF-$\beta $, are important in determining the long-term dynamics and the development of disease. Finally, we demonstrate the utility of this model in investigating temporal responses to bronchial thermoplasty, a therapeutic intervention in which ASM is ablated by applying thermal energy to the airway wall. The model predicts the parameter-dependent threshold damage required to obtain irreversible reduction in ASM content, suggesting that certain asthma phenotypes are more likely to benefit from this intervention.

Differential remodeling in small and large murine airways revealed by novel whole lung airway analysis.

Airway remodeling occurs in chronic asthma leading to increased airway smooth muscle (ASM) mass and extracellular matrix (ECM) deposition. Although extensively studied in murine airways, studies report only selected larger airways at one time-point meaning the spatial distribution and resolution of remodeling are poorly understood. Here we use a new method allowing comprehensive assessment of the spatial and temporal changes in ASM, ECM, and epithelium in large numbers of murine airways after allergen challenge. Using image processing to analyze 20-50 airways per mouse from a whole lung section revealed increases in ASM and ECM after allergen challenge were greater in small and large rather than intermediate airways. ASM predominantly accumulated adjacent to the basement membrane, whereas ECM was distributed across the airway wall. Epithelial hyperplasia was most marked in small and intermediate airways. After challenge, ASM changes resolved over 7 days, whereas ECM and epithelial changes persisted. The new method suggests large and small airways remodel differently, and the long-term consequences of airway inflammation may depend more on ECM and epithelial changes than ASM. The improved quantity and quality of unbiased data provided by the method reveals important spatial differences in remodeling and could set new analysis standards for murine asthma models.

A computational cardiopulmonary physiology simulator accurately predicts individual patient responses to changes in mechanical ventilator settings.

We present new results validating the capability of a high-fidelity computational simulator to accurately predict the responses of individual patients with acute respiratory distress syndrome to changes in mechanical ventilator settings. 26 pairs of data-points comprising arterial blood gasses collected before and after changes in inspiratory pressure, PEEP, FiO2, and I:E ratio from six mechanically ventilated patients were used for this study. Parallelized global optimization algorithms running on a high-performance computing cluster were used to match the simulator to each initial data point. Mean absolute percentage errors between the simulator predicted values of PaO2 and PaCO2 and the patient data after changing ventilator parameters were 10.3% and 12.6%, respectively. Decreasing the complexity of the simulator by reducing the number of independent alveolar compartments reduced the accuracy of its predictions. Clinical Relevance- These results provide further evidence that our computational simulator can accurately reproduce patient responses to mechanical ventilation, highlighting its usefulness as a clinical research tool.

Modeling Mechanical Ventilation In Silico-Potential and Pitfalls.

Computer simulation offers a fresh approach to traditional medical research that is particularly well suited to investigating issues related to mechanical ventilation. Patients receiving mechanical ventilation are routinely monitored in great detail, providing extensive high-quality data-streams for model design and configuration. Models based on such data can incorporate very complex system dynamics that can be validated against patient responses for use as investigational surrogates. Crucially, simulation offers the potential to "look inside" the patient, allowing unimpeded access to all variables of interest. In contrast to trials on both animal models and human patients, in silico models are completely configurable and reproducible; for example, different ventilator settings can be applied to an identical virtual patient, or the same settings applied to different patients, to understand their mode of action and quantitatively compare their effectiveness. Here, we review progress on the mathematical modeling and computer simulation of human anatomy, physiology, and pathophysiology in the context of mechanical ventilation, with an emphasis on the clinical applications of this approach in various disease states. We present new results highlighting the link between model complexity and predictive capability, using data on the responses of individual patients with acute respiratory distress syndrome to changes in multiple ventilator settings. The current limitations and potential of in silico modeling are discussed from a clinical perspective, and future challenges and research directions highlighted.

The Critical Importance of Spatial and Temporal Scales in Designing and Interpreting Immune Cell Migration Assays.

Intravital microscopy and other direct-imaging techniques have allowed for a characterisation of leukocyte migration that has revolutionised the field of immunology, resulting in an unprecedented understanding of the mechanisms of immune response and adaptive immunity. However, there is an assumption within the field that modern imaging techniques permit imaging parameters where the resulting cell track accurately captures a cell's motion. This notion is almost entirely untested, and the relationship between what could be observed at a given scale and the underlying cell behaviour is undefined. Insufficient spatial and temporal resolutions within migration assays can result in misrepresentation of important physiologic processes or cause subtle changes in critical cell behaviour to be missed. In this review, we contextualise how scale can affect the perceived migratory behaviour of cells, summarise the limited approaches to mitigate this effect, and establish the need for a widely implemented framework to account for scale and correct observations of cell motion. We then extend the concept of scale to new approaches that seek to bridge the current "black box" between single-cell behaviour and systemic response.

Reduced biomechanical models for precision-cut lung-slice stretching experiments.

Precision-cut lung-slices (PCLS), in which viable airways embedded within lung parenchyma are stretched or induced to contract, are a widely used ex vivo assay to investigate bronchoconstriction and, more recently, mechanical activation of pro-remodelling cytokines in asthmatic airways. We develop a nonlinear fibre-reinforced biomechanical model accounting for smooth muscle contraction and extracellular matrix strain-stiffening. Through numerical simulation, we describe the stresses and contractile responses of an airway within a PCLS of finite thickness, exposing the importance of smooth muscle contraction on the local stress state within the airway. We then consider two simplifying limits of the model (a membrane representation and an asymptotic reduction in the thin-PCLS-limit), that permit analytical progress. Comparison against numerical solution of the full problem shows that the asymptotic reduction successfully captures the key elements of the full model behaviour. The more tractable reduced model that we develop is suitable to be employed in investigations to elucidate the time-dependent feedback mechanisms linking airway mechanics and cytokine activation in asthma.

Switching behaviour in vascular smooth muscle cell-matrix adhesion during oscillatory loading.

Integrins regulate mechanotransduction between smooth muscle cells (SMCs) and the extracellular matrix (ECM). SMCs resident in the walls of airways or blood vessels are continuously exposed to dynamic mechanical forces due to breathing or pulsatile blood flow. However, the resulting effects of these forces on integrin dynamics and associated cell-matrix adhesion are not well understood. Here we present experimental results from atomic force microscopy (AFM) experiments, designed to study the integrin response to external oscillatory loading of varying amplitudes applied to live aortic SMCs, together with theoretical results from a mathematical model. In the AFM experiments, a fibronectin-coated probe was used cyclically to indent and retract from the surface of the cell. We observed a transition between states of firm adhesion and of complete detachment as the amplitude of oscillatory loading increased, revealed by qualitative changes in the force timecourses. Interestingly, for some of the SMCs in the experiments, switching behaviour between the two adhesion states is observed during single timecourses at intermediate amplitudes. We obtain two qualitatively similar adhesion states in the mathematical model, where we simulate the cell, integrins and ECM as an evolving system of springs, incorporating local integrin binding dynamics. In the mathematical model, we observe a region of bistability where both the firm adhesion and detachment states can occur depending on the initial adhesion state. The differences are seen to be a result of mechanical cooperativity of integrins and cell deformation. Switching behaviour is a phenomenon associated with bistability in a stochastic system, and bistability in our deterministic mathematical model provides a potential physical explanation for the experimental results. Physiologically, bistability provides a means for transient mechanical stimuli to induce long-term changes in adhesion dynamics-and thereby the cells' ability to transmit force-and we propose further experiments for testing this hypothesis.

An Integrated Pipeline for Combining in vitro Data and Mathematical Models Using a Bayesian Parameter Inference Approach to Characterize Spatio-temporal Chemokine Gradient Formation.

All protective and pathogenic immune and inflammatory responses rely heavily on leukocyte migration and localization. Chemokines are secreted chemoattractants that orchestrate the positioning and migration of leukocytes through concentration gradients. The mechanisms underlying chemokine gradient establishment and control include physical as well as biological phenomena. Mathematical models offer the potential to both understand this complexity and suggest interventions to modulate immune function. Constructing models that have powerful predictive capability relies on experimental data to estimate model parameters accurately, but even with a reductionist approach most experiments include multiple cell types, competing interdependent processes and considerable uncertainty. Therefore, we propose the use of reduced modeling and experimental frameworks in complement, to minimize the number of parameters to be estimated. We present a Bayesian optimization framework that accounts for advection and diffusion of a chemokine surrogate and the chemokine CCL19, transport processes that are known to contribute to the establishment of spatio-temporal chemokine gradients. Three examples are provided that demonstrate the estimation of the governing parameters as well as the underlying uncertainty. This study demonstrates how a synergistic approach between experimental and computational modeling benefits from the Bayesian approach to provide a robust analysis of chemokine transport. It provides a building block for a larger research effort to gain holistic insight and generate novel and testable hypotheses in chemokine biology and leukocyte trafficking.

DP2 antagonism reduces airway smooth muscle mass in asthma by decreasing eosinophilia and myofibroblast recruitment.

Increased airway smooth muscle mass, a feature of airway remodeling in asthma, is the strongest predictor of airflow limitation and contributes to asthma-associated morbidity and mortality. No current drug therapy for asthma is known to affect airway smooth muscle mass. Although there is increasing evidence that prostaglandin D2 type 2 receptor (DP2) is expressed in airway structural and inflammatory cells, few studies have addressed the expression and function of DP2 in airway smooth muscle cells. We report that the DP2 antagonist fevipiprant reduced airway smooth muscle mass in bronchial biopsies from patients with asthma who had participated in a previous randomized placebo-controlled trial. We developed a computational model to capture airway remodeling. Our model predicted that a reduction in airway eosinophilia alone was insufficient to explain the clinically observed decrease in airway smooth muscle mass without a concomitant reduction in the recruitment of airway smooth muscle cells or their precursors to airway smooth muscle bundles that comprise the airway smooth muscle layer. We experimentally confirmed that airway smooth muscle migration could be inhibited in vitro using DP2-specific antagonists in an airway smooth muscle cell culture model. Our analyses suggest that fevipiprant, through antagonism of DP2, reduced airway smooth muscle mass in patients with asthma by decreasing airway eosinophilia in concert with reduced recruitment of myofibroblasts and fibrocytes to the airway smooth muscle bundle. Fevipiprant may thus represent a potential therapy to ameliorate airway remodeling in asthma.

Small airway hyperresponsiveness in COPD: relationship between structure and function in lung slices.

The direct relationship between pulmonary structural changes and airway hyperresponsiveness (AHR) in chronic obstructive pulmonary disease (COPD) is unclear. We investigated AHR in relation to airway and parenchymal structural changes in a guinea pig model of COPD and in COPD patients. Precision-cut lung slices (PCLS) were prepared from guinea pigs challenged with lipopolysaccharide or saline two times weekly for 12 wk. Peripheral PCLS were obtained from patients with mild to moderate COPD and non-COPD controls. AHR to methacholine was measured in large and small airways using video-assisted microscopy. Airway smooth muscle mass and alveolar airspace size were determined in the same slices. A mathematical model was used to identify potential changes in biomechanical properties underlying AHR. In guinea pigs, lipopolysaccharide increased the sensitivity of large (>150 µm) airways toward methacholine by 4.4-fold and the maximal constriction of small airways (<150 µm) by 1.5-fold. Similarly increased small airway responsiveness was found in COPD patients. In both lipopolysaccharide-challenged guinea pigs and patients, airway smooth muscle mass was unaltered, whereas increased alveolar airspace correlated with small airway hyperresponsiveness in guinea pigs. Fitting the parameters of the model indicated that COPD weakens matrix mechanical properties and enhances stiffness differences between the airway and the parenchyma, in both species. In conclusion, this study demonstrates small airway hyperresponsiveness in PCLS from COPD patients. These changes may be related to reduced parenchymal retraction forces and biomechanical changes in the airway wall. PCLS from lipopolysaccharide-exposed guinea pigs may be useful to study mechanisms of small airway hyperresponsiveness in COPD.

Chemokine Transport Dynamics and Emerging Recognition of Their Role in Immune Function.

Leukocyte migration is critically important during all protective and pathological immune and inflammatory responses. Chemokines play fundamental roles in this process, and chemokine concentration gradients stimulate the directional migration of leukocytes. The formation and regulation of these gradients is poorly understood. These are complex processes that depend on the specific properties of each chemokine and interactions between physical, biological and biochemical processes, including production, diffusion, advection, scavenging, post-translational modification, and extracellular matrix (ECM) binding. While some of these mechanisms have been investigated in isolation or limited combinations, more integrative research is required to provide a quantitative knowledge base that explains how chemokine gradients are established and maintained, and how cells respond to, and modify, these gradients.

A theoretical model of inflammation- and mechanotransduction-driven asthmatic airway remodelling.

Inflammation, airway hyper-responsiveness and airway remodelling are well-established hallmarks of asthma, but their inter-relationships remain elusive. In order to obtain a better understanding of their inter-dependence, we develop a mechanochemical morphoelastic model of the airway wall accounting for local volume changes in airway smooth muscle (ASM) and extracellular matrix in response to transient inflammatory or contractile agonist challenges. We use constrained mixture theory, together with a multiplicative decomposition of growth from the elastic deformation, to model the airway wall as a nonlinear fibre-reinforced elastic cylinder. Local contractile agonist drives ASM cell contraction, generating mechanical stresses in the tissue that drive further release of mitogenic mediators and contractile agonists via underlying mechanotransductive signalling pathways. Our model predictions are consistent with previously described inflammation-induced remodelling within an axisymmetric airway geometry. Additionally, our simulations reveal novel mechanotransductive feedback by which hyper-responsive airways exhibit increased remodelling, for example, via stress-induced release of pro-mitogenic and pro-contractile cytokines. Simulation results also reveal emergence of a persistent contractile tone observed in asthmatics, via either a pathological mechanotransductive feedback loop, a failure to clear agonists from the tissue, or a combination of both. Furthermore, we identify various parameter combinations that may contribute to the existence of different asthma phenotypes, and we illustrate a combination of factors which may predispose severe asthmatics to fatal bronchospasms.

Effect of Loading History on Airway Smooth Muscle Cell-Matrix Adhesions.

Integrin-mediated adhesions between airway smooth muscle (ASM) cells and the extracellular matrix (ECM) regulate how contractile forces generated within the cell are transmitted to its external environment. Environmental cues are known to influence the formation, size, and survival of cell-matrix adhesions, but it is not yet known how they are affected by dynamic fluctuations associated with tidal breathing in the intact airway. Here, we develop two closely related theoretical models to study adhesion dynamics in response to oscillatory loading of the ECM, representing the dynamic environment of ASM cells in vivo. Using a discrete stochastic-elastic model, we simulate individual integrin binding and rupture events and observe two stable regimes in which either bond formation or bond rupture dominate, depending on the amplitude of the oscillatory loading. These regimes have either a high or low fraction of persistent adhesions, which could affect the level of strain transmission between contracted ASM cells and the airway tissue. For intermediate loading, we observe a region of bistability and hysteresis due to shared loading between existing bonds; the level of adhesion depends on the loading history. These findings are replicated in a related continuum model, which we use to investigate the effect of perturbations mimicking deep inspirations (DIs). Because of the bistability, a DI applied to the high adhesion state could either induce a permanent switch to a lower adhesion state or allow a return of the system to the high adhesion state. Transitions between states are further influenced by the frequency of oscillations, cytoskeletal or ECM stiffnesses, and binding affinities, which modify the magnitudes of the stable adhesion states as well as the region of bistability. These findings could explain (in part) the transient bronchodilatory effect of a DI observed in asthmatics compared to a more sustained effect in normal subjects.

In vitro, in silico and in vivo study challenges the impact of bronchial thermoplasty on acute airway smooth muscle mass loss.

Bronchial thermoplasty is a treatment for asthma. It is currently unclear whether its histopathological impact is sufficiently explained by the proportion of airway wall that is exposed to temperatures necessary to affect cell survival.Airway smooth muscle and bronchial epithelial cells were exposed to media (37-70°C) for 10 s to mimic thermoplasty. In silico we developed a mathematical model of airway heat distribution post-thermoplasty. In vivo we determined airway smooth muscle mass and epithelial integrity pre- and post-thermoplasty in 14 patients with severe asthma.In vitro airway smooth muscle and epithelial cell number decreased significantly following the addition of media heated to ≥65°C. In silico simulations showed a heterogeneous heat distribution that was amplified in larger airways, with <10% of the airway wall heated to >60°C in airways with an inner radius of ∼4 mm. In vivo at 6 weeks post-thermoplasty, there was an improvement in asthma control (measured via Asthma Control Questionnaire-6; mean difference 0.7, 95% CI 0.1-1.3; p=0.03), airway smooth muscle mass decreased (absolute median reduction 5%, interquartile range (IQR) 0-10; p=0.03) and epithelial integrity increased (14%, IQR 6-29; p=0.007). Neither of the latter two outcomes was related to improved asthma control.Integrated in vitro and in silico modelling suggest that the reduction in airway smooth muscle post-thermoplasty cannot be fully explained by acute heating, and nor did this reduction confer a greater improvement in asthma control.

Repeated airway constrictions in mice do not alter respiratory function.

It is suggested that the frequent strain the airways undergo in asthma because of repeated airway smooth muscle (ASM)-mediated constrictions contributes to airway wall remodeling. However, the effects of repeated constrictions on airway remodeling, as well as the ensuing impact of this presumptive remodeling on respiratory mechanics, have never been investigated in subjects without asthma. In this study, we set out to determine whether repeated constrictions lead to features that are reminiscent of asthma in mice without asthma. BALB/c mice were subjected to a 30-min constriction elicited by aerosolized methacholine every other day over 6 wk. Forty-eight hours after the last constriction, the mechanics of the respiratory system was evaluated at baseline and in response to incremental doses of nebulized methacholine with the flexiVent. The whole-lung lavages, the tracheas, and the lungs were also collected to evaluate inflammation, the contractile capacity of ASM, and the structural components of the airway wall, respectively. The resistance and the compliance of the respiratory system, as well as the Newtonian resistance and the resistive and elastic properties of the lung tissue, were not affected by repeated constrictions, both at baseline and in response to methacholine. All the other examined features also remained unaltered, except the number of goblet cells in the epithelium and the number of macrophages in the whole-lung lavages, which both increased with repeated constrictions. This study demonstrates that, despite causing goblet cell hyperplasia and a mild macrophagic inflammation, repeated constrictions with methacholine do not lead to structural changes that adversely impact the physiology. NEW & NOTEWORTHY Repeated airway constrictions led to signs of remodeling that are typically observed in asthma, which neither altered respiratory mechanics nor the contractile capacity of airway smooth muscle. These findings shed light on a debate between those claiming that constrictions induce remodeling and those convinced that methacholine challenges are harmless. Insofar as our results with mice relate to humans, the findings indicate that repeated challenges with methacholine can be performed safely.

A Novel Computational Model Predicts Key Regulators of Chemokine Gradient Formation in Lymph Nodes and Site-Specific Roles for CCL19 and ACKR4.

The chemokine receptor CCR7 drives leukocyte migration into and within lymph nodes (LNs). It is activated by chemokines CCL19 and CCL21, which are scavenged by the atypical chemokine receptor ACKR4. CCR7-dependent navigation is determined by the distribution of extracellular CCL19 and CCL21, which form concentration gradients at specific microanatomical locations. The mechanisms underpinning the establishment and regulation of these gradients are poorly understood. In this article, we have incorporated multiple biochemical processes describing the CCL19-CCL21-CCR7-ACKR4 network into our model of LN fluid flow to establish a computational model to investigate intranodal chemokine gradients. Importantly, the model recapitulates CCL21 gradients observed experimentally in B cell follicles and interfollicular regions, building confidence in its ability to accurately predict intranodal chemokine distribution. Parameter variation analysis indicates that the directionality of these gradients is robust, but their magnitude is sensitive to these key parameters: chemokine production, diffusivity, matrix binding site availability, and CCR7 abundance. The model indicates that lymph flow shapes intranodal CCL21 gradients, and that CCL19 is functionally important at the boundary between B cell follicles and the T cell area. It also predicts that ACKR4 in LNs prevents CCL19/CCL21 accumulation in efferent lymph, but does not control intranodal gradients. Instead, it attributes the disrupted interfollicular CCL21 gradients observed in Ackr4-deficient LNs to ACKR4 loss upstream. Our novel approach has therefore generated new testable hypotheses and alternative interpretations of experimental data. Moreover, it acts as a framework to investigate gradients at other locations, including those that cannot be visualized experimentally or involve other chemokines.

Corrigendum: Airway and Parenchymal Strains during Bronchoconstriction in the Precision Cut Lung Slice.

[This corrects the article on p. 309 in vol. 7, PMID: 27559314.].

Mathematical modelling of cytokines, MMPs and fibronectin fragments in osteoarthritic cartilage.

Osteoarthritis (OA) is a degenerative disease which causes pain and stiffness in joints. OA progresses through excessive degradation of joint cartilage, eventually leading to significant joint degeneration and loss of function. Cytokines, a group of cell signalling proteins, present in raised concentrations in OA joints, can be classified into pro-inflammatory and anti-inflammatory groups. They mediate cartilage degradation through several mechanisms, primarily the up-regulation of matrix metalloproteinases (MMPs), a group of collagen-degrading enzymes. In this paper we show that the interactions of cytokines within cartilage have a crucial role to play in OA progression and treatment. We develop a four-variable ordinary differential equation model for the interactions between pro- and anti-inflammatory cytokines, MMPs and fibronectin fragments (Fn-fs), a by-product of cartilage degradation and up-regulator of cytokines. We show that the model has four classes of dynamic behaviour: homoeostasis, bistable inflammation, tristable inflammation and persistent inflammation. We show that positive and negative feedbacks controlling cytokine production rates can determine either a pre-disposition to OA or initiation of OA. Further, we show that manipulation of cytokine, MMP and Fn-fs levels can be used to treat OA, but we suggest that multiple treatment targets may be essential to halt or slow disease progression.

Airway and Parenchymal Strains during Bronchoconstriction in the Precision Cut Lung Slice.

The precision-cut lung slice (PCLS) is a powerful tool for studying airway reactivity, but biomechanical measurements to date have largely focused on changes in airway caliber. Here we describe an image processing tool that reveals the associated spatio-temporal changes in airway and parenchymal strains. Displacements of sub-regions within the PCLS are tracked in phase-contrast movies acquired after addition of contractile and relaxing drugs. From displacement maps, strains are determined across the entire PCLS or along user-specified directions. In a representative mouse PCLS challenged with 10(-4)M methacholine, as lumen area decreased, compressive circumferential strains were highest in the 50 µm closest to the airway lumen while expansive radial strains were highest in the region 50-100 µm from the lumen. However, at any given distance from the airway the strain distribution varied substantially in the vicinity of neighboring small airways and blood vessels. Upon challenge with the relaxant agonist chloroquine, although most strains disappeared, residual positive strains remained a long time after addition of chloroquine, predominantly in the radial direction. Taken together, these findings establish strain mapping as a new tool to elucidate local dynamic mechanical events within the constricting airway and its supporting parenchyma.