The interdependent contributions of gravitational and structural features to perfusion distribution in a multiscale model of the pulmonary circulation

Physics-based in silico modelling of microvascular pulmonary perfusion in COVID-19

Due to its ability to induce heterogenous, patient-specific damage in pulmonary alveoli and capillaries, COVID-19 poses challenges in defining a uniform profile to elucidate infection across all patients. Computational models that integrate changes in ventilation and perfusion with heterogeneous damage profiles offer valuable insights into the impact of COVID-19 on pulmonary health. This study aims to develop an in silico hypothesis-testing platform specifically focused on studying microvascular pulmonary perfusion in COVID-19-infected lungs. Through this platform, we explore the effects of various acinar-level pulmonary perfusion abnormalities on global lung function. Our modelling approach simulates changes in pulmonary perfusion and the resulting mismatch of ventilation and perfusion in COVID-19-afflicted lungs. Using this coupled modelling platform, we conducted multiple simulations to assess different scenarios of perfusion abnormalities in COVID-19-infected lungs. The simulation results showed an overall decrease in ventilation-perfusion (V/Q) ratio with inclusion of various types of perfusion abnormalities such as hypoperfusion with and without microangiopathy. This model serves as a foundation for comprehending and comparing the spectrum of findings associated with COVID-19 in the lung, paving the way for patient-specific modelling of microscale lung damage in emerging pulmonary pathologies like COVID-19.

Predicting Patient Status in Chronic Thromboembolic Pulmonary Hypertension Using a Biophysical Model

Chronic thromboembolic pulmonary hypertension (CTEPH) involves abnormally high blood pressure in the pulmonary vessels and is associated with small vessel vasculopathy and pre-capillary proximal occlusions. Management of CTEPH disease is challenging, therefore accurate diagnosis is crucial in ensuring effective treatment and improved patient outcomes. The treatment of choice for CTEPH is pulmonary endarterectomy, which is an invasive surgical intervention to remove thrombi. Following PEA, a number of patients experience poor outcomes or worse-than-expected improvements, which may indicate that they have significant small vessel disease. A method that can predict the extent of distal remodelling may provide useful clinical information to plan appropriate CTEPH patient treatment. Here, a novel biophysical modelling approach has been developed to estimate and quantify the extent of distal remodelling. This method includes a combination of mathematical modelling and computed tomography pulmonary angiography to first model the geometry of the pulmonary arteries and to identify the under-perfused regions in CTEPH. The geometric model is then used alongside haemodynamic measurements from right heart catheterisation to predict distal remodelling. In this study, the method is tested and validated using synthetically generated remodelling data. Then, a preliminary application of this technique to patient data is shown to demonstrate the potential of the approach for use in the clinical setting.Clinical relevance- Patient-specific modelling can help provide useful information regarding the extent of distal vasculopathy on a per-patient basis, which remains challenging. Physicians can be unsure of outcomes following pulmonary endarterectomy. Therefore, the predictive aspect of the patient's response to surgery can help with clinical decision-making.

Roadmap for an imaging and modelling paediatric study in rural NZ

Our study methodology is motivated from three disparate needs: one, imaging studies have existed in silo and study organs but not across organ systems; two, there are gaps in our understanding of paediatric structure and function; three, lack of representative data in New Zealand. Our research aims to address these issues in part, through the combination of magnetic resonance imaging, advanced image processing algorithms and computational modelling. Our study demonstrated the need to take an organ-system approach and scan multiple organs on the same child. We have pilot tested an imaging protocol to be minimally disruptive to the children and demonstrated state-of-the-art image processing and personalized computational models using the imaging data. Our imaging protocol spans brain, lungs, heart, muscle, bones, abdominal and vascular systems. Our initial set of results demonstrated child-specific measurements on one dataset. This work is novel and interesting as we have run multiple computational physiology workflows to generate personalized computational models. Our proposed work is the first step towards achieving the integration of imaging and modelling improving our understanding of the human body in paediatric health and disease.

An in silico approach to understanding the interaction between cardiovascular and pulmonary lymphatic dysfunction

The lung is extremely sensitive to interstitial fluid balance, yet the role of pulmonary lymphatics in lung fluid homeostasis and its interaction with cardiovascular pressures is poorly understood. In health, there is a fine balance between fluid extravasated from the pulmonary capillaries into the interstitium and the return of fluid to the circulation via the lymphatic vessels. This balance is maintained by an extremely interdependent system governed by pressures in the fluids (air and blood) and tissue (interstitium), lung motion during breathing, and the permeability of the tissues. Chronic elevation in left atrial pressure (LAP) due to left heart disease increases the capillary blood pressure. The consequent fluid accumulation in the delicate lung tissue increases its weight, decreases its compliance, and impairs gas exchange. This interdependent system is difficult, if not impossible, to study experimentally. Computational modeling provides a unique perspective to analyze fluid movement in the cardiopulmonary vasculature in health and disease. We have developed an initial in silico model of pulmonary lymphatic function using an anatomically structured model to represent ventilation and perfusion and underlying biophysical laws governing fluid transfer at the interstitium. This novel model was tested against increased LAP and noncardiogenic effects (increased permeability). The model returned physiologically reasonable values for all applications, predicting pulmonary edema when LAP reached 25 mmHg and with increased permeability.NEW & NOTEWORTHY This model presents a novel approach to understanding the interaction between cardiac dysfunction and pulmonary lymphatic function, using anatomically structured models and biophysical equations to estimate regional variation in fluid transport from blood to interstitial and lymphatic flux. This fluid transport model brings together advanced models of ventilation, perfusion, and lung mechanics to produce a detailed model of fluid transport in health and various altered pathological conditions.

Aerosol Transport Modeling: The Key Link Between Lung Infections of Individuals and Populations

The recent COVID-19 pandemic has propelled the field of aerosol science to the forefront, particularly the central role of virus-laden respiratory droplets and aerosols. The pandemic has also highlighted the critical need, and value for, an information bridge between epidemiological models (that inform policymakers to develop public health responses) and within-host models (that inform the public and health care providers how individuals develop respiratory infections). Here, we review existing data and models of generation of respiratory droplets and aerosols, their exhalation and inhalation, and the fate of infectious droplet transport and deposition throughout the respiratory tract. We then articulate how aerosol transport modeling can serve as a bridge between and guide calibration of within-host and epidemiological models, forming a comprehensive tool to formulate and test hypotheses about respiratory tract exposure and infection within and between individuals.

A century of exercise physiology: key concepts on coupling respiratory oxygen flow to muscle energy demand during exercise

After a short historical account, and a discussion of Hill and Meyerhof’s theory of the energetics of muscular exercise, we analyse steady-state rest and exercise as the condition wherein coupling of respiration to metabolism is most perfect. The quantitative relationships show that the homeostatic equilibrium, centred around arterial pH of 7.4 and arterial carbon dioxide partial pressure of 40 mmHg, is attained when the ratio of alveolar ventilation to carbon dioxide flow (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\dot{V}}_{A}/{\dot{V}}_{R}{CO}_{2}$$\end{document}V˙A/V˙RCO2) is − 21.6. Several combinations, exploited during exercise, of pertinent respiratory variables are compatible with this equilibrium, allowing adjustment of oxygen flow to oxygen demand without its alteration. During exercise transients, the balance is broken, but the coupling of respiration to metabolism is preserved when, as during moderate exercise, the respiratory system responds faster than the metabolic pathways. At higher exercise intensities, early blood lactate accumulation suggests that the coupling of respiration to metabolism is transiently broken, to be re-established when, at steady state, blood lactate stabilizes at higher levels than resting. In the severe exercise domain, coupling cannot be re-established, so that anaerobic lactic metabolism also contributes to sustain energy demand, lactate concentration goes up and arterial pH falls continuously. The \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\dot{V}}_{A}/{\dot{V}}_{R}{CO}_{2}$$\end{document}V˙A/V˙RCO2 decreases below − 21.6, because of ensuing hyperventilation, while lactate keeps being accumulated, so that exercise is rapidly interrupted. The most extreme rupture of the homeostatic equilibrium occurs during breath-holding, because oxygen flow from ambient air to mitochondria is interrupted. No coupling at all is possible between respiration and metabolism in this case.

Simulating Multi-Scale Pulmonary Vascular Function by Coupling Computational Fluid Dynamics With an Anatomic Network Model

The function of the pulmonary circulation is truly multi-scale, with blood transported through vessels from centimeter to micron scale. There are scale-dependent mechanisms that govern the flow in the pulmonary vascular system. However, very few computational models of pulmonary hemodynamics capture the physics of pulmonary perfusion across the spatial scales of functional importance in the lung. Here we present a multi-scale model that incorporates the 3-dimensional (3D) complexities of pulmonary blood flow in the major vessels, coupled to an anatomically-based vascular network model incorporating the multiple contributing factors to capillary perfusion, including gravity. Using the model we demonstrate how we can predict the impact of vascular remodeling and occlusion on both macro-scale functional drivers (flow distribution between lungs, and wall shear stress) and micro-scale contributors to gas exchange. The model predicts interactions between 3D and 1D models that lead to a redistribution of blood between postures, both on a macro- and a micro-scale. This allows us to estimate the effect of posture on left and right pulmonary artery wall shear stress, with predictions varying by 0.75-1.35 dyne/cm² between postures.

A computational model of contributors to pulmonary hypertensive disease: impacts of whole lung and focal disease distributions

Pulmonary hypertension has multiple etiologies and so can be difficult to diagnose, prognose, and treat. Diagnosis is typically made via invasive hemodynamic measurements in the main pulmonary artery and is based on observed elevation of mean pulmonary artery pressure. This static mean pressure enables diagnosis, but does not easily allow assessment of the severity of pulmonary hypertension, nor the etiology of the disease, which may impact treatment. Assessment of the dynamic properties of pressure and flow data obtained from catheterization potentially allows more meaningful assessment of the strain on the right heart and may help to distinguish between disease phenotypes. However, mechanistic understanding of how the distribution of disease in the lung leading to pulmonary hypertension impacts the dynamics of blood flow in the main pulmonary artery and/or the pulmonary capillaries is lacking. We present a computational model of the pulmonary vasculature, parameterized to characteristic features of pulmonary arterial hypertension and chronic thromboembolic pulmonary hypertension to help understand how the two conditions differ in terms of pulmonary vascular response to disease. Our model incorporates key features known to contribute to pulmonary vascular function in health and disease, including anatomical structure and multiple contributions from gravity. The model suggests that dynamic measurements obtained from catheterization potentially distinguish between distal and proximal vasculopathy typical of pulmonary arterial hypertension and chronic thromboembolic pulmonary hypertension. However, the model suggests a non-linear relationship between these data and vascular structural changes typical of pulmonary arterial hypertension and chronic thromboembolic pulmonary hypertension which may impede analysis of these metrics to distinguish between cohorts.

Integrative Computational Models of Lung Structure-Function Interactions

Anatomically based integrative models of the lung and their interaction with other key components of the respiratory system provide unique capabilities for investigating both normal and abnormal lung function. There is substantial regional variability in both structure and function within the normal lung, yet it remains capable of relatively efficient gas exchange by providing close matching of air delivery (ventilation) and blood delivery (perfusion) to regions of gas exchange tissue from the scale of the whole organ to the smallest continuous gas exchange units. This is despite remarkably different mechanisms of air and blood delivery, different fluid properties, and unique scale-dependent anatomical structures through which the blood and air are transported. This inherent heterogeneity can be exacerbated in the presence of disease or when the body is under stress. Current computational power and data availability allow for the construction of sophisticated data-driven integrative models that can mimic respiratory system structure, function, and response to intervention. Computational models do not have the same technical and ethical issues that can limit experimental studies and biomedical imaging, and if they are solidly grounded in physiology and physics they facilitate investigation of the underlying interaction between mechanisms that determine respiratory function and dysfunction, and to estimate otherwise difficult-to-access measures. © 2021 American Physiological Society. Compr Physiol 11:1501-1530, 2021.

In Silico Ventilation Within the Dose-Volume is Predictive of Lung Function Post-radiation Therapy in Patients with Lung Cancer

Lung cancer is a leading cause of death worldwide. Radiation therapy (RT) is one method to treat this disease. A common side effect of RT for lung cancer is radiation-induced lung damage (RILD) which leads to loss of lung function. RILD often compounds pre-existing smoking-related regional lung function impairment. It is difficult to predict patient outcomes due to large variability in individual response to RT. In this study, the capability of image-based modelling of regional ventilation in lung cancer patients to predict lung function post-RT was investigated. Twenty-five patient-based models were created using CT images to define the airway geometry, size and location of tumour, and distribution of emphysema. Simulated ventilation within the 20 Gy isodose volume showed a statistically significant negative correlation with the change in forced expiratory volume in 1 s 12-months post-RT (p = 0.001, R = - 0.61). Patients with higher simulated ventilation within the 20 Gy isodose volume had a greater loss in lung function post-RT and vice versa. This relationship was only evident with the combined impact of tumour and emphysema, with the location of the emphysema relative to the dose-volume being important. Our results suggest that model-based ventilation measures can be used in the prediction of patient lung function post-RT.

Wave reflection in an anatomical model of the pulmonary circulation in local and global hypertensive disease

Pulmonary hypertension is a disease of the pulmonary vasculature which can occur for many different reasons, including pathological remodeling of the pulmonary vessels and occlusion of these vessels (amongst others). Pulmonary hypertension can lead to right heart failure and significantly reduces the quality of life of patients living with the condition. It is difficult to distinguish clinically between different classifications of pulmonary hypertension, and doing so accurately is critical for the management of an individual's condition. In addition, different presentations of the disease (e.g. occlusion versus remodeling) can put different strains on the right heart, despite patients having very similar elevations in pulmonary artery pressure. In this study we use an anatomically based model of the pulmonary circulation to predict pressure and flow wave transmission and reflection in two different kinds of pulmonary hypertension - primary pulmonary hypertension, and chronic thrombotic pulmonary embolism (CTEPH), to enable analysis of the impact of disease type on impedance spectra in the main pulmonary artery.

The Lung Physiome and virtual patient models: From morphometry to clinical translation

The understanding or prediction of specific functions of the lung can be made using compact models that have identifiable parameters and that are custom designed to the problem of interest. However, when structure contributes to function - as is the case with most lung pathologies - structure-based, biophysical models become essential. Here we describe the application of structure-based models within the lung Physiome framework to identifying and explaining patient risk in 12patients diagnosed with acute pulmonary embolism. The model integrates perfusion, ventilation, and gas exchange to predict arterial blood gases and pulmonary artery pressure in individual patient models in response to patient-specific blood clot distribution, with full or partial arterial occlusion. The necessity for a patient-specific approach with biophysical models that account for scale-specific structure and function is demonstrated.

Ventilation/Perfusion Matching: Of Myths, Mice, and Men

Despite a huge range in lung size between species, there is little measured difference in the ability of the lung to provide a well-matched air flow (ventilation) to blood flow (perfusion) at the gas exchange tissue. Here, we consider the remarkable similarities in ventilation/perfusion matching between species through a biophysical lens and consider evidence that matching in large animals is dominated by gravity but in small animals by structure.

Do pulmonary cavity shapes influence lung function?

Distribution of lung tissue within the chest cavity is a key contributor to delivery of both blood and air to the gas exchange regions of the lung. This distribution is multifactorial with influences from parenchyma, gravity and level of inflation. We hypothesize that the manner in which lung inflates, for example, the primarily diaphragmatic nature of normal breathing, is an important contributor to regional lung tissue distribution. To investigate this hypothesis, we present an organ-level model of lung tissue mechanics which incorporates pleural cavity change due to change in lung volume or posture. We quantify the changes using shape and density metrics in ten healthy subjects scanned supine at end-inspiratory and end-expiratory volumes and ten subjects scanned at both supine and prone end-inspiratory volumes. Comparing end-expiratory to end-inspiratory volume, we see primarily a change in the cranial-caudal dimension of the lung, reflective of movement of diaphragm. In the diaphragmatic region there is greater regional lung expansion than in the cranial aspect, which is restricted by the chest wall. When moving from supine to prone, a restriction of lung was observed anteriorly, resulting in a generally reduced lung volume and a redistribution of air volume posteriorly. In general, we see the highest in lung tissue density heterogeneity in regions of the lung that are most inflated. Using our computational model, we quantify the impact of pleural cavity shape change on regional lung distribution, and predict the impact on regional elastic recoil pressure.

Pulmonary Vascular Dynamics

The pulmonary circulation carries deoxygenated blood from the systemic veins through the pulmonary arteries to be oxygenated in the capillaries that line the walls of the pulmonary alveoli. The pulmonary circulation carries the cardiac output with a relatively low driving pressure, and so differs considerably in structure and function from the systemic circulation to maintain a low-resistance vascular system. The pulmonary circulation is often considered to be a quasi-static system in both experimental and computational studies of pulmonary perfusion and its matching to ventilation (air flow) for exchange. However, the system is highly dynamic, with cardiac output and regional perfusion changing with posture, exercise, and over time. Here we review this dynamic system, with a focus on understanding the physiology of pulmonary vascular dynamics across spatial and temporal scales, and the changes to these dynamics that are reflective of disease. © 2019 American Physiological Society. Compr Physiol 9:1081-1100, 2019.

Optimizing human pulmonary perfusion measurement using an in silico model of arterial spin labeling magnetic resonance imaging

Arterial spin labeling (ASL) magnetic resonance imaging (MRI) is an imaging methodology that uses blood as an endogenous contrast agent to quantify flow. One limitation of this method of capillary blood quantification when applied in the lung is the contribution of signals from non-capillary blood. Intensity thresholding is one approach that has been proposed for minimizing the non-capillary blood signal. This method has been tested in previous in silico modeling studies; however, it has only been tested under a restricted set of physiological conditions (supine posture and a cardiac output of 5 L/min). This study presents an in silico approach that extends previous intensity thresholding analysis to estimate the optimal "per-slice" intensity threshold value using the individual components of the simulated ASL signal (signal arising independently from capillary blood as well as pulmonary arterial and pulmonary venous blood). The aim of this study was to assess whether the threshold value should vary with slice location, posture, or cardiac output. We applied an in silico modeling approach to predict the blood flow distribution and the corresponding ASL quantification of pulmonary perfusion in multiple sagittal imaging slices. There was a significant increase in ASL signal and heterogeneity (COV = 0.90 to COV = 1.65) of ASL signals when slice location changed from lateral to medial. Heterogeneity of the ASL signal within a slice was significantly lower (P = 0.03) in prone (COV = 1.08) compared to in the supine posture (COV = 1.17). Increasing stroke volume resulted in an increase in ASL signal and conversely an increase in heart rate resulted in a decrease in ASL signal. However, when cardiac output was increased via an increase in both stroke volume and heart rate, ASL signal remained relatively constant. Despite these differences, we conclude that a threshold value of 35% provides optimal removal of large vessel signal independent of slice location, posture, and cardiac output.

Healthy Lung Vessel Morphology Derived From Thoracic Computed Tomography

Knowledge of the lung vessel morphology in healthy subjects is necessary to improve our understanding about the functional network of the lung and to recognize pathologic deviations beyond the normal inter-subject variation. Established values of normal lung morphology have been derived from necropsy material of only very few subjects. In order to determine morphologic readouts from a large number of healthy subjects, computed tomography pulmonary angiography (CTPA) datasets, negative for pulmonary embolism, and other thoracic pathologies, were analyzed using a fully-automatic, in-house developed artery/vein separation algorithm. The number, volume, and tortuosity of the vessels in a diameter range between 2 and 10 mm were determined. Visual inspection of all datasets was used to exclude subjects with poor image quality or inadequate artery/vein separation from the analysis. Validation of the algorithm was performed manually by a radiologist on randomly selected subjects. In 123 subjects (men/women: 55/68), aged 59 ± 17 years, the median overlap between visual inspection and fully-automatic segmentation was 94.6% (69.2–99.9%). The median number of vessel segments in the ranges of 8–10, 6–8, 4–6, and 2–4 mm diameter was 9, 34, 134, and 797, respectively. Number of vessel segments divided by the subject's lung volume was 206 vessels/L with arteries and veins contributing almost equally. In women this vessel density was about 15% higher than in men. Median arterial and venous volumes were 1.52 and 1.54% of the lung volume, respectively. Tortuosity was best described with the sum-of-angles metric and was 142.1 rad/m (138.3–144.5 rad/m). In conclusion, our fully-automatic artery/vein separation algorithm provided reliable measures of pulmonary arteries and veins with respect to age and gender. There was a large variation between subjects in all readouts. No relevant dependence on age, gender, or vessel type was observed. These data may provide reference values for morphometric analysis of lung vessels.

In silico modeling of oxygen-enhanced MRI of specific ventilation

Specific ventilation imaging (SVI) proposes that using oxygen-enhanced 1H MRI to capture signal change as subjects alternatively breathe room air and 100% O2 provides an estimate of specific ventilation distribution in the lung. How well this technique measures SV and the effect of currently adopted approaches of the technique on resulting SV measurement is open for further exploration. We investigated (1) How well does imaging a single sagittal lung slice represent whole lung SV? (2) What is the influence of pulmonary venous blood on the measured MRI signal and resultant SVI measure? and (3) How does inclusion of misaligned images affect SVI measurement? In this study, we utilized two patient-based in silico models of ventilation, perfusion, and gas exchange to address these questions for normal healthy lungs. Simulation results from the two healthy young subjects show that imaging a single slice is generally representative of whole lung SV distribution, with a calculated SV gradient within 90% of that calculated for whole lung distributions. Contribution of O2 from the venous circulation results in overestimation of SV at a regional level where major pulmonary veins cross the imaging plane, resulting in a 10% increase in SV gradient for the imaging slice. A worst-case scenario simulation of image misalignment increased the SV gradient by 11.4% for the imaged slice.

Capturing complexity in pulmonary system modelling

Respiratory disease is a significant problem worldwide, and it is a problem with increasing prevalence. Pathology in the upper airways and lung is very difficult to diagnose and treat, as response to disease is often heterogeneous across patients. Computational models have long been used to help understand respiratory function, and these models have evolved alongside increases in the resolution of medical imaging and increased capability of functional imaging, advances in biological knowledge, mathematical techniques and computational power. The benefits of increasingly complex and realistic geometric and biophysical models of the respiratory system are that they are able to capture heterogeneity in patient response to disease and predict emergent function across spatial scales from the delicate alveolar structures to the whole organ level. However, with increasing complexity, models become harder to solve and in some cases harder to validate, which can reduce their impact clinically. Here, we review the evolution of complexity in computational models of the respiratory system, including successes in translation of models into the clinical arena. We also highlight major challenges in modelling the respiratory system, while making use of the evolving functional data that are available for model parameterisation and testing.

Gravity outweighs the contribution of structure to passive ventilation-perfusion matching in the supine adult human lung

Gravity and matched airway/vascular tree geometries are both hypothesized to be key contributors to ventilation-perfusion (V̇/Q̇) matching in the lung, but their relative contributions are challenging to quantify experimentally. We used a structure-based model to conduct an analysis of the relative contributions of tissue deformation (the "Slinky" effect), other gravitational mechanisms (weight of blood and gravitational gradient in tissue elastic recoil), and matched airway and arterial tree geometry to V̇/Q̇ matching and therefore to total lung oxygen exchange. Our results showed that the heterogeneity in V̇ and Q̇ were lowest and the correlation between V̇ and Q̇ was highest when the only mechanism for V̇/Q̇ matching was either tissue deformation or matched geometry. Heterogeneity in V̇ and Q̇ was highest and their correlation was poorest when all mechanisms were active (that is, at baseline). Eliminating the contribution of matched geometry did not change the correlation between V̇ and Q̇ at baseline. Despite the much larger heterogeneities in V̇ and Q̇ at baseline, the contribution of in-common (to V̇ and Q̇) gravitational mechanisms provided sufficient compensatory V̇/Q̇ matching to minimize the impact on oxygen transfer. In summary, this model predicts that during supine normal breathing under gravitational loading, passive V̇/Q̇ matching is predominantly determined by shared gravitationally induced tissue deformation, compliance distribution, and the effect of the hydrostatic pressure gradient on vessel and capillary size and blood pressures. Contribution from the matching airway and arterial tree geometries in this model is minor under normal gravity in the supine adult human lung. NEW & NOTEWORTHY We use a computational model to systematically analyze contributors to ventilation-perfusion matching in the lung. The model predicts that the multiple effects of gravity are the predominant mechanism in providing passive ventilation-perfusion matching in the supine adult human lung under normal gravitational loads, while geometric matching of airway and arterial trees plays a minor role.

The physiology of submaximal exercise: The steady state concept

The steady state concept implies that the oxygen flow is invariant and equal at each level along the respiratory system. The same is the case with the carbon dioxide flow. This condition has several physiological consequences, which are analysed. First, we briefly discuss the mechanical efficiency of exercise and the energy cost of human locomotion, as well as the roles played by aerodynamic work and frictional work. Then we analyse the equations describing the oxygen flow in lungs and in blood, the effects of ventilation and of the ventilation - perfusion inequality, and the interaction between diffusion and perfusion in the lungs. The cardiovascular responses sustaining gas flow increase in blood are finally presented. An equation linking ventilation, circulation and metabolism is developed, on the hypothesis of constant oxygen flow in mixed venous blood. This equation tells that, if the pulmonary respiratory quotient stays invariant, any increase in metabolic rate is matched by a proportional increase in ventilation, but by a less than proportional increase in cardiac output.

Image-based computational fluid dynamics in the lung: virtual reality or new clinical practice?

The development and implementation of personalized medicine is paramount to improving the efficiency and efficacy of patient care. In the respiratory system, function is largely dictated by the choreographed movement of air and blood to the gas exchange surface. The passage of air begins in the upper airways, either via the mouth or nose, and terminates at the alveolar interface, while blood flows from the heart to the alveoli and back again. Computational fluid dynamics (CFD) is a well-established tool for predicting fluid flows and pressure distributions within complex systems. Traditionally CFD has been used to aid in the effective or improved design of a system or device; however, it has become increasingly exploited in biological and medical-based applications further broadening the scope of this computational technique. In this review, we discuss the advancement in application of CFD to the respiratory system and the contributions CFD is currently making toward improving precision medicine. The key areas CFD has been applied to in the pulmonary system are in predicting fluid transport and aerosol distribution within the airways. Here we focus our discussion on fluid flows and in particular on image-based clinically focused CFD in the ventilatory system. We discuss studies spanning from the paranasal sinuses through the conducting airways down to the level of the alveolar airways. The combination of imaging and CFD is enabling improved device design in aerosol transport, improved biomarkers of lung function in clinical trials, and improved predictions and assessment of surgical interventions in the nasal sinuses. WIREs Syst Biol Med 2017, 9:e1392. doi: 10.1002/wsbm.1392 For further resources related to this article, please visit the WIREs website.

The effect of lung deformation on the spatial distribution of pulmonary blood flow

KEY POINTS

Pulmonary perfusion measurement using magnetic resonance imaging combined with deformable image registration enabled us to quantify the change in the spatial distribution of pulmonary perfusion at different lung volumes. The current study elucidated the effects of tidal volume lung inflation [functional residual capacity (FRC) + 500 ml and FRC + 1 litre] on the change in pulmonary perfusion distribution. Changes in hydrostatic pressure distribution as well as transmural pressure distribution due to the change in lung height with tidal volume inflation are probably bigger contributors to the redistribution of pulmonary perfusion than the changes in pulmonary vasculature resistance caused by lung tissue stretch.

ABSTRACT

Tidal volume lung inflation results in structural changes in the pulmonary circulation, potentially affecting pulmonary perfusion. We hypothesized that perfusion is recruited to regions receiving the greatest deformation from a tidal breath, thus ensuring ventilation-perfusion matching. Density-normalized perfusion (DNP) magnetic resonance imaging data were obtained in healthy subjects (n = 7) in the right lung at functional residual capacity (FRC), FRC+500 ml, and FRC+1.0 l. Using deformable image registration, the displacement of a sagittal lung slice acquired at FRC to the larger volumes was calculated. Registered DNP images were normalized by the mean to estimate perfusion redistribution (nDNP). Data were evaluated across gravitational regions (dependent, middle, non-dependent) and by lobes (upper, RUL; middle, RML; lower, RLL). Lung inflation did not alter mean DNP within the slice (P = 0.10). The greatest expansion was seen in the dependent region (P < 0.0001: dependent vs non-dependent, P < 0.0001: dependent vs middle) and RLL (P = 0.0015: RLL vs RUL, P < 0.0001: RLL vs RML). Neither nDNP recruitment to RLL [+500 ml = -0.047(0.145), +1 litre = 0.018(0.096)] nor to dependent lung [+500 ml = -0.058(0.126), +1 litre = -0.023(0.106)] were found. Instead, redistribution was seen in decreased nDNP in the non-dependent [+500 ml = -0.075(0.152), +1 litre = -0.137(0.167)) and increased nDNP in the gravitational middle lung [+500 ml = 0.098(0.058), +1 litre = 0.093(0.081)] (P = 0.01). However, there was no significant lobar redistribution (P < 0.89). Contrary to our hypothesis, based on the comparison between gravitational and lobar perfusion data, perfusion was not redistributed to the regions of the most inflation. This suggests that either changes in hydrostatic pressure or transmural pressure distribution in the gravitational direction are implicated in the redistribution of perfusion away from the non-dependent lung.

Multi-scale Modeling of the Cardiovascular System: Disease Development, Progression, and Clinical Intervention

Cardiovascular diseases (CVDs) are the leading cause of death in the western world. With the current development of clinical diagnostics to more accurately measure the extent and specifics of CVDs, a laudable goal is a better understanding of the structure-function relation in the cardiovascular system. Much of this fundamental understanding comes from the development and study of models that integrate biology, medicine, imaging, and biomechanics. Information from these models provides guidance for developing diagnostics, and implementation of these diagnostics to the clinical setting, in turn, provides data for refining the models. In this review, we introduce multi-scale and multi-physical models for understanding disease development, progression, and designing clinical interventions. We begin with multi-scale models of cardiac electrophysiology and mechanics for diagnosis, clinical decision support, personalized and precision medicine in cardiology with examples in arrhythmia and heart failure. We then introduce computational models of vasculature mechanics and associated mechanical forces for understanding vascular disease progression, designing clinical interventions, and elucidating mechanisms that underlie diverse vascular conditions. We conclude with a discussion of barriers that must be overcome to provide enhanced insights, predictions, and decisions in pre-clinical and clinical applications.

Multiscale modelling of the feto-placental vasculature

The placenta provides all the nutrients required for the fetus through pregnancy. It develops dynamically, and, to avoid rejection of the fetus, there is no mixing of fetal and maternal blood; rather, the branched placental villi 'bathe' in blood supplied from the uterine arteries. Within the villi, the feto-placental vasculature also develops a complex branching structure in order to maximize exchange between the placental and maternal circulations. To understand the development of the placenta, we must translate functional information across spatial scales including the interaction between macro- and micro-scale haemodynamics and account for the effects of a dynamically and rapidly changing structure through the time course of pregnancy. Here, we present steps towards an anatomically based and multiscale approach to modelling the feto-placental circulation. We assess the effect of the location of cord insertion on feto-placental blood flow resistance and flow heterogeneity and show that, although cord insertion does not appear to directly influence feto-placental resistance, the heterogeneity of flow in the placenta is predicted to increase from a 19.4% coefficient of variation with central cord insertion to 23.3% when the cord is inserted 2 cm from the edge of the placenta. Model geometries with spheroidal and ellipsoidal shapes, but the same volume, showed no significant differences in flow resistance or heterogeneity, implying that normal asymmetry in shape does not affect placental efficiency. However, the size and number of small capillary vessels is predicted to have a large effect on feto-placental resistance and flow heterogeneity. Using this new model as an example, we highlight the importance of taking an integrated multi-disciplinary and multiscale approach to understand development of the placenta.

Comparison of generic and subject-specific models for simulation of pulmonary perfusion and forced expiration

The goal of translating multiscale model analysis of pulmonary function into population studies is challenging because of the need to derive a geometric model for each subject. This could be addressed by using a generic model with appropriate customization to subject-specific data. Here, we present a quantitative comparison of simulating two fundamental behaviours of the lung-its haemodynamic response to vascular occlusion, and the forced expiration in 1 s (FEV1) following bronchoconstriction-in subject-specific and generic models. When the subjects are considered as a group, there is no significant difference between predictions of mean pulmonary artery pressure (mPAP), pulmonary vascular resistance or forced expiration; however, significant differences are apparent in the prediction of arterial oxygen, for both baseline and post-occlusion. Despite the apparent consistency of the generic and subject-specific models, a third of subjects had generic model under-prediction of the increase in mPAP following occlusion, and half had the decrease in arterial oxygen over-predicted; two subjects had considerable differences in the percentage reduction of FEV1 following bronchoconstriction. The generic model approach can be useful for physiologically directed studies but is not appropriate for simulating pathophysiological function that is strongly dependent on interaction with lung structure.

Computational modeling of the obstructive lung diseases asthma and COPD

Asthma and chronic obstructive pulmonary disease (COPD) are characterized by airway obstruction and airflow limitation and pose a huge burden to society. These obstructive lung diseases impact the lung physiology across multiple biological scales. Environmental stimuli are introduced via inhalation at the organ scale, and consequently impact upon the tissue, cellular and sub-cellular scale by triggering signaling pathways. These changes are propagated upwards to the organ level again and vice versa. In order to understand the pathophysiology behind these diseases we need to integrate and understand changes occurring across these scales and this is the driving force for multiscale computational modeling.

There is an urgent need for improved diagnosis and assessment of obstructive lung diseases. Standard clinical measures are based on global function tests which ignore the highly heterogeneous regional changes that are characteristic of obstructive lung disease pathophysiology. Advances in scanning technology such as hyperpolarized gas MRI has led to new regional measurements of ventilation, perfusion and gas diffusion in the lungs, while new image processing techniques allow these measures to be combined with information from structural imaging such as Computed Tomography (CT). However, it is not yet known how to derive clinical measures for obstructive diseases from this wealth of new data. Computational modeling offers a powerful approach for investigating this relationship between imaging measurements and disease severity, and understanding the effects of different disease subtypes, which is key to developing improved diagnostic methods.

Gaining an understanding of a system as complex as the respiratory system is difficult if not impossible via experimental methods alone. Computational models offer a complementary method to unravel the structure-function relationships occurring within a multiscale, multiphysics system such as this. Here we review the current state-of-the-art in techniques developed for pulmonary image analysis, development of structural models of the respiratory system and predictions of function within these models. We discuss application of modeling techniques to obstructive lung diseases, namely asthma and emphysema and the use of models to predict response to therapy. Finally we introduce a large European project, AirPROM that is developing multiscale models to investigate structure-function relationships in asthma and COPD.

High temporal versus high spatial resolution in MR quantitative pulmonary perfusion imaging of two-year old children after congenital diaphragmatic hernia repair

OBJECTIVES

Congenital diaphragmatic hernia (CDH) leads to lung hypoplasia. Using dynamic contrast-enhanced (DCE) MR imaging, lung perfusion can be quantified. As MR perfusion values depend on temporal resolution, we compared two protocols to investigate whether ipsilateral lung perfusion is impaired after CDH, whether there are protocol-dependent differences, and which protocol is preferred.

METHODS

DCE-MRI was performed in 36 2-year old children after CDH on a 3 T MRI system; protocol A (n = 18) based on a high spatial (3.0 s; voxel: 1.25 mm(3)) and protocol B (n = 18) on a high temporal resolution (1.5 s; voxel: 2 mm(3)). Pulmonary blood flow (PBF), pulmonary blood volume (PBV), mean transit time (MTT), and peak-contrast-to-noise-ratio (PCNR) were quantified.

RESULTS

PBF was reduced ipsilaterally, with ipsilateral PBF of 45 ± 26 ml/100 ml/min to contralateral PBF of 63 ± 28 ml/100 ml/min (p = 0.0016) for protocol A; and for protocol B, side differences were equivalent (ipsilateral PBF = 62 ± 24 vs. contralateral PBF = 85 ± 30 ml/100 ml/min; p = 0.0034). PCNR was higher for protocol B (30 ± 18 vs. 20 ± 9; p = 0.0294). Protocol B showed higher values of PBF in comparison to protocol A (p always <0.05).

CONCLUSIONS

Ipsilateral lung perfusion is reduced in 2-year old children following CDH repair. Higher temporal resolution and increased voxel size show a gain in PCNR and lead to higher perfusion values. Protocol B is therefore preferred.

KEY POINTS

• Quantitative lung perfusion parameters depend on temporal and spatial resolution. • Reduction of lung perfusion in CDH can be measured with different MR protocols. • Temporal resolution of 1.5 s with spatial resolution of 2 mm (3) is suitable.

Maximal oxygen consumption in healthy humans: theories and facts

This article reviews the concept of maximal oxygen consumption ([Formula: see text]) from the perspective of multifactorial models of [Formula: see text] limitation. First, I discuss procedural aspects of [Formula: see text] measurement: the implications of ramp protocols are analysed within the theoretical work of Morton. Then I analyse the descriptive physiology of [Formula: see text], evidencing the path that led to the view of monofactorial cardiovascular or muscular [Formula: see text] limitation. Multifactorial models, generated by the theoretical work of di Prampero and Wagner around the oxygen conductance equation, represented a radical change of perspective. These models are presented in detail and criticized with respect to the ensuing experimental work. A synthesis between them is proposed, demonstrating how much these models coincide and converge on the same conclusions. Finally, I discuss the cases of hypoxia and bed rest, the former as an example of the pervasive effects of the shape of the oxygen equilibrium curve, the latter as a neat example of adaptive changes concerning the entire respiratory system. The conclusion is that the concept of cardiovascular [Formula: see text] limitation is reinforced by multifactorial models, since cardiovascular oxygen transport provides most of the [Formula: see text] limitation, at least in normoxia. However, the same models show that the role of peripheral resistances is significant and cannot be neglected. The role of peripheral factors is greater the smaller is the active muscle mass. In hypoxia, the intervention of lung resistances as limiting factors restricts the role played by cardiovascular and peripheral factors.

Hypoxic pulmonary vasoconstriction as a contributor to response in acute pulmonary embolism

Hypoxic pulmonary vasoconstriction (HPV) is an adaptive response unique to the lung whereby blood flow is diverted away from areas of low alveolar oxygen to improve ventilation-perfusion matching and resultant gas exchange. Some previous experimental studies have suggested that the HPV response to hypoxia is blunted in acute pulmonary embolism (APE), while others have concluded that HPV contributes to elevated pulmonary blood pressures in APE. To understand these contradictory observations, we have used a structure-based computational model of integrated lung function in 10 subjects to study the impact of HPV on pulmonary hemodynamics and gas exchange in the presence of regional arterial occlusion. The integrated model includes an experimentally-derived model for HPV. Its function is validated against measurements of pulmonary vascular resistance in normal subjects at four levels of inspired oxygen. Our results show that the apparently disparate observations of previous studies can be explained within a single model: the model predicts that HPV increases mean pulmonary artery pressure in APE (by 8.2 ± 7.0% in these subjects), and concurrently shows a reduction in response to hypoxia in the subjects who have high levels of occlusion and therefore maximal HPV in normoxia.

Multi-scale computational models of the airways to unravel the pathophysiological mechanisms in asthma and chronic obstructive pulmonary disease (AirPROM)

THE RESPIRATORY SYSTEM COMPRISES SEVERAL SCALES OF BIOLOGICAL COMPLEXITY: the genes, cells and tissues that work in concert to generate resultant function. Malfunctions of the structure or function of components at any spatial scale can result in diseases, to the detriment of gas exchange, right heart function and patient quality of life. Vast amounts of data emerge from studies across each of the biological scales; however, the question remains: how can we integrate and interpret these data in a meaningful way? Respiratory disease presents a huge health and economic burden, with the diseases asthma and chronic obstructive pulmonary disease (COPD) affecting over 500 million people worldwide. Current therapies are inadequate owing to our incomplete understanding of the disease pathophysiology and our lack of recognition of the enormous disease heterogeneity: we need to characterize this heterogeneity on a patient-specific basis to advance healthcare. In an effort to achieve this goal, the AirPROM consortium (Airway disease Predicting Outcomes through patient-specific computational Modelling) brings together a multi-disciplinary team and a wealth of clinical data. Together we are developing an integrated multi-scale model of the airways in order to unravel the complex pathophysiological mechanisms occurring in the diseases asthma and COPD.

Dynamic blood flow and wall shear stress in pulmonary hypertensive disease

This study provides new model of pulsatile flow in the pulmonary circulation in health and pulmonary hypertensive disease. Structural vascular remodeling typical of pulmonary hypertensive disease was implemented in the model by progressively altering the mechanical properties of the arterial geometry and progressively increasing the inlet pulse pressure (PP). The transmission of PP throughout the tree was shown to increase in advanced stages of disease, creating the potential for a `vicious-cycle' of damage to vasculature. Wall shear stress (WSS) was shown to be highest in the terminal arteries of the model and increased significantly with disease. A further trend observed in WSS results was that high WSS values began to `climb' the arterial tree towards the proximal vessels as disease progressed. This suggests a link between WSS and distal remodeling in pulmonary hypertensive disease, which initiates in the small muscular arteries and arterioles and spreads into larger arteries as the disease progresses.

Integrative approaches for modeling regulation and function of the respiratory system

Mathematical models have been central to understanding the interaction between neural control and breathing. Models of the entire respiratory system-which comprises the lungs and the neural circuitry that controls their ventilation-have been derived using simplifying assumptions to compartmentalize each component of the system and to define the interactions between components. These full system models often rely-through necessity-on empirically derived relationships or parameters, in addition to physiological values. In parallel with the development of whole respiratory system models are mathematical models that focus on furthering a detailed understanding of the neural control network, or of the several functions that contribute to gas exchange within the lung. These models are biophysically based, and rely on physiological parameters. They include single-unit models for a breathing lung or neural circuit, through to spatially distributed models of ventilation and perfusion, or multicircuit models for neural control. The challenge is to bring together these more recent advances in models of neural control with models of lung function, into a full simulation for the respiratory system that builds upon the more detailed models but remains computationally tractable. This requires first understanding the mathematical models that have been developed for the respiratory system at different levels, and which could be used to study how physiological levels of O2 and CO2 in the blood are maintained.

Blood Flow in Capillaries of the Human Lung

A novel model for the blood system is postulated focusing on the flow rate and pressure distribution inside the arterioles and venules of the pulmonary acinus. Based upon physiological data it is devoid of any ad hoc constants. The model comprises nine generations of arterioles, venules, and capillaries in the acinus, the gas exchange unit of the lung. Blood is assumed incompressible and Newtonian and the blood vessels are assumed inextensible. Unlike previous models of the blood system, the venules and arterioles open up to the capillary network in numerous locations along each generation. The large number of interconnected capillaries is perceived as a porous medium in which the flow is macroscopically unidirectional from arterioles to venules openings. In addition, the large number of capillaries extending from each arteriole and venule allows introduction of a continuum theory and formulation of a novel system of ordinary, nonlinear differential equations which governs the blood flow and pressure fields along the arterioles, venules, and capillaries. The solution of the differential equations is semianalytical and requires the inversion of three diagonal, 9 × 9 matrices only. The results for the total flow rate of blood through the acinus are within the ballpark of physiological observations despite the simplifying assumptions used in our model. The results also manifest that the contribution of the nonlinear convection term of the Navier-Stokes equations has little effect (less than 2%) on the total blood flow entering/leaving the acinus despite the fact that the Reynolds number is not much smaller than unity at the proximal generations. The model makes it possible to examine some pathological cases. Here, centri-acinar and distal emphysema were investigated yielding a reduction in inlet blood flow rate.

Ventilatory and cardiac responses to pulmonary embolism: consequences for gas exchange and blood pressure

Acute thromboembolic pulmonary embolism (PE) is a life threatening condition that can lead to pulmonary hypertension and right ventricular dysfunction or failure. There is typically an increase in ventilation rate and cardiac output as a response to PE prior to cardiac failure, which is at least in part due to systemic hypoxemia. Here we assess the response of the lungs to changes in these parameters using anatomically-based computational models of pulmonary perfusion, ventilation and gas exchange. We show that increases in ventilation and cardiac output improve overall gas exchange in PE. However, this comes at the cost of an increased pulmonary blood pressure, which may contribute to pulmonary hypertension as a result of PE.

Blood flow redistribution and ventilation-perfusion mismatch during embolic pulmonary arterial occlusion

Acute pulmonary embolism causes redistribution of blood in the lung, which impairs ventilation/perfusion matching and gas exchange and can elevate pulmonary arterial pressure (PAP) by increasing pulmonary vascular resistance (PVR). An anatomically-based multi-scale model of the human pulmonary circulation was used to simulate pre- and post-occlusion flow, to study blood flow redistribution in the presence of an embolus, and to evaluate whether reduction in perfused vascular bed is sufficient to increase PAP to hypertensive levels, or whether other vasoconstrictive mechanisms are necessary. A model of oxygen transfer from air to blood was included to assess the impact of vascular occlusion on oxygen exchange. Emboli of 5, 7, and 10 mm radius were introduced to occlude increasing proportions of the vasculature. Blood flow redistribution was calculated after arterial occlusion, giving predictions of PAP, PVR, flow redistribution, and micro-circulatory flow dynamics. Because of the large flow reserve capacity (via both capillary recruitment and distension), approximately 55% of the vasculature was occluded before PAP reached clinically significant levels indicative of hypertension. In contrast, model predictions showed that even relatively low levels of occlusion could cause localized oxygen deficit. Flow preferentially redistributed to gravitationally non-dependent regions regardless of occlusion location, due to the greater potential for capillary recruitment in this region. Red blood cell transit times decreased below the minimum time for oxygen saturation (<0.25 s) and capillary pressures became high enough to initiate cell damage (which may result in edema) only after ~80% of the lung was occluded.

What can imaging tell us about physiology? Lung growth and regional mechanical strain

The interplay of mechanical forces transduces diverse physico-biochemical processes to influence lung morphogenesis, growth, maturation, remodeling and repair. Because tissue stress is difficult to measure in vivo, mechano-sensitive responses are commonly inferred from global changes in lung volume, shape, or compliance and correlated with structural changes in tissue blocks sampled from postmortem-fixed lungs. Recent advances in noninvasive volumetric imaging technology, nonrigid image registration, and deformation analysis provide valuable tools for the quantitative analysis of in vivo regional anatomy and air and tissue-blood distributions and when combined with transpulmonary pressure measurements, allow characterization of regional mechanical function, e.g., displacement, strain, shear, within and among intact lobes, as well as between the lung and the components of its container-rib cage, diaphragm, and mediastinum-thereby yielding new insights into the inter-related metrics of mechanical stress-strain and growth/remodeling. Here, we review the state-of-the-art imaging applications for mapping asymmetric heterogeneous physical interactions within the thorax and how these interactions permit as well as constrain lung growth, remodeling, and compensation during development and following pneumonectomy to illustrate how advanced imaging could facilitate the understanding of physiology and pathophysiology. Functional imaging promises to facilitate the formulation of realistic computational models of lung growth that integrate mechano-sensitive events over multiple spatial and temporal scales to accurately describe in vivo physiology and pathophysiology. Improved computational models in turn could enhance our ability to predict regional as well as global responses to experimental and therapeutic interventions.

A computational model of the topographic distribution of ventilation in healthy human lungs

The topographic distribution of ventilation in the lungs is determined by the interaction of several factors, including lung shape, airway tree geometry, posture, and tissue deformation. Inter-species differences in lung structure-function and technical difficulty in obtaining high resolution imaging of the upright human lung means that it is not straightforward to experimentally determine the contribution of each of these factors to ventilation distribution. We present a mathematical model for predicting the topological distribution of inhaled air in the upright healthy human lung, based on anatomically structured model geometries and biophysical equations for model function. Gravitational deformation of the lung tissue is predicted using a continuum model. Airflow is simulated in anatomically based conducting airways coupled to geometrically simplified terminal acinar units with varying volume-dependent compliances. The predicted ventilation distribution is hence governed by local tissue density and elastic recoil pressure, airway resistance and acinar compliance. Results suggest that there is significant spatial variation in intrinsic tissue properties in the lungs. The model confirms experimental evidence that in the healthy lungs tissue compliance has a far greater effect than airway resistance on the spatial distribution of ventilation, and hence a realistic description of tissue deformation is essential in models of ventilation.

Assessing potential errors of MRI-based measurements of pulmonary blood flow using a detailed network flow model

MRI images of pulmonary blood flow using arterial spin labeling (ASL) measure the delivery of magnetically tagged blood to an image plane during one systolic ejection period. However, the method potentially suffers from two problems, each of which may depend on the imaging plane location: 1) the inversion plane is thicker than the imaging plane, resulting in a gap that blood must cross to be detected in the image; and 2) ASL includes signal contributions from tagged blood in conduit vessels (arterial and venous). By using an in silico model of the pulmonary circulation we found the gap reduced the ASL signal to 64-74% of that in the absence of a gap in the sagittal plane and 53-84% in the coronal. The contribution of the conduit vessels varied markedly as a function of image plane ranging from ∼90% of the overall signal in image planes that encompass the central hilar vessels to <20% in peripheral image planes. A threshold cutoff removing voxels with intensities >35% of maximum reduced the conduit vessel contribution to the total ASL signal to ∼20% on average; however, planes with large contributions from conduit vessels underestimate acinar flow due to a high proportion of in-plane flow, making ASL measurements of perfusion impractical. In other image planes, perfusion dominated the resulting ASL images with good agreement between ASL and acinar flow. Similarly, heterogeneity of the ASL signal as measured by relative dispersion is a reliable measure of heterogeneity of the acinar flow distribution in the same image planes.

Contralateral decubitus positioning enhances computed tomographic angiographic evaluation of pulmonary vasculature in a patient with a pulmonary arteriovenous malformation

Computed tomographic pulmonary angiography has become the diagnostic procedure of choice in patients suspected of having a pulmonary embolus. However, intrapulmonary shunting of blood in a variety of pathologic conditions can cause suboptimal opacification of the pulmonary arterial circulation and result in a suboptimal or even nondiagnostic study. Radiologists should be aware of these conditions and be familiar with positioning techniques to minimize such shunting. We report a patient suspected of having pulmonary embolism, in whom a preexisting unilateral arteriovenous malformation prevented adequate evaluation of the pulmonary circulation. Positioning the patient in the contralateral decubitus position significantly enhanced image quality.

Computational modeling of airway and pulmonary vascular structure and function: development of a "lung physiome"

Computational models of lung structure and function necessarily span multiple spatial and temporal scales, i.e., dynamic molecular interactions give rise to whole organ function, and the link between these scales cannot be fully understood if only molecular or organ-level function is considered. Here, we review progress in constructing multiscale finite element models of lung structure and function that are aimed at providing a computational framework for bridging the spatial scales from molecular to whole organ. These include structural models of the intact lung, embedded models of the pulmonary airways that couple to model lung tissue, and models of the pulmonary vasculature that account for distinct structural differences at the extra- and intra-acinar levels. Biophysically based functional models for tissue deformation, pulmonary blood flow, and airway bronchoconstriction are also described. The development of these advanced multiscale models has led to a better understanding of complex physiological mechanisms that govern regional lung perfusion and emergent heterogeneity during bronchoconstriction.

Pulmonary embolism: predicting disease severity

Pulmonary embolism (PE) is the most common cause of acute pulmonary hypertension, yet it is commonly undiagnosed, with risk of death if not recognized promptly and managed accordingly. Patients typically present with hypoxemia and hypomania, although the presentation varies greatly, being confounded by co-morbidities such as pre-existing cardio-respiratory disease. Previous studies have demonstrated variable patient outcomes in spite of similar extent and distribution of pulmonary vascular occlusion, but the path physiological determinants of outcome remain unclear. Computational models enable exact control over many of the compounding factors leading to functional outcomes and therefore provide a useful tool to understand and assess these mechanisms. We review the current state of pulmonary blood flow models. We present a pilot study within 10 patients presenting with acute PE, where patient-derived vascular occlusions are imposed onto an existing model of the pulmonary circulation enabling predictions of resultant haemodynamic after embolus occlusion. Results show that mechanical obstruction alone is not sufficient to cause pulmonary arterial hypertension, even when up to 65 per cent of lung tissue is occluded. Blood flow is found to preferentially redistribute to the gravitationally non-dependent regions. The presence of an additional downstream occlusion is found to significantly increase pressures.

Computational models of the pulmonary circulation: Insights and the move towards clinically directed studies

Biophysically-based computational models provide a tool for integrating and explaining experimental data, observations, and hypotheses. Computational models of the pulmonary circulation have evolved from minimal and efficient constructs that have been used to study individual mechanisms that contribute to lung perfusion, to sophisticated multi-scale and -physics structure-based models that predict integrated structure-function relationships within a heterogeneous organ. This review considers the utility of computational models in providing new insights into the function of the pulmonary circulation, and their application in clinically motivated studies. We review mathematical and computational models of the pulmonary circulation based on their application; we begin with models that seek to answer questions in basic science and physiology and progress to models that aim to have clinical application. In looking forward, we discuss the relative merits and clinical relevance of computational models: what important features are still lacking; and how these models may ultimately be applied to further increasing our understanding of the mechanisms occurring in disease of the pulmonary circulation.