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.

Lung Computational Models and the Role of the Small Airways in Asthma.

Rationale: Asthma is characterized by disease within the small airways. Several studies have suggested that forced oscillation technique-derived resistance at 5 Hz (R5) - resistance at 20 Hz (R20) is a measure of small airway disease; however, there has been limited validation of this measurement to date.Objectives: To validate the use of forced oscillation R5 - R20 as a measure of small airway narrowing in asthma, and to investigate the role that small airway narrowing plays in asthma.Methods: Patient-based complete conducting airway models were generated from computed tomography scans to simulate the impact of different degrees of airway narrowing at different levels of the airway tree on forced oscillation R5 - R20 (n = 31). The computational models were coupled with regression models in an asthmatic cohort (n = 177) to simulate the impact of small airway narrowing on asthma control and quality of life. The computational models were used to predict the impact on small airway narrowing of type-2 targeting biologics using pooled data from two similarly design randomized, placebo-controlled biologic trials (n = 137).Measurements and Main Results: Simulations demonstrated that narrowing of the small airways had a greater impact on R5 - R20 than narrowing of the larger airways and was associated (above a threshold of approximately 40% narrowing) with marked deterioration in both asthma control and asthma quality of life, above the minimal clinical important difference. The observed treatment effect on R5 - R20 in the pooled trials equated to a predicted small airway narrowing reversal of approximately 40%.Conclusions: We have demonstrated, using computational modeling, that forced oscillation R5 - R20 is a direct measure of anatomical narrowing in the small airways and that small airway narrowing has a marked impact on both asthma control and quality of life and may be modified by biologics.

Computational modeling of the obstructive lung diseases asthma and COPD.

Asthma and chronic obstructive pulmonary disease (COPD) are characterized by airway obstruction and airflow imitation 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 currentstate-of-the-art in techniques developed for pulmonary image analysis, development of structural models of therespiratory 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 toinvestigate structure-function relationships in asthma and COPD.

Systems Medicine: from molecular features and models to the clinic in COPD.

BACKGROUND AND HYPOTHESIS: Chronic Obstructive Pulmonary Disease (COPD) patients are characterized by heterogeneous clinical manifestations and patterns of disease progression. Two major factors that can be used to identify COPD subtypes are muscle dysfunction/wasting and co-morbidity patterns. We hypothesized that COPD heterogeneity is in part the result of complex interactions between several genes and pathways. We explored the possibility of using a Systems Medicine approach to identify such pathways, as well as to generate predictive computational models that may be used in clinic practice. OBJECTIVE AND METHOD: Our overarching goal is to generate clinically applicable predictive models that characterize COPD heterogeneity through a Systems Medicine approach. To this end we have developed a general framework, consisting of three steps/objectives: (1) feature identification, (2) model generation and statistical validation, and (3) application and validation of the predictive models in the clinical scenario. We used muscle dysfunction and co-morbidity as test cases for this framework. RESULTS: In the study of muscle wasting we identified relevant features (genes) by a network analysis and generated predictive models that integrate mechanistic and probabilistic models. This allowed us to characterize muscle wasting as a general de-regulation of pathway interactions. In the co-morbidity analysis we identified relevant features (genes/pathways) by the integration of gene-disease and disease-disease associations. We further present a detailed characterization of co-morbidities in COPD patients that was implemented into a predictive model. In both use cases we were able to achieve predictive modeling but we also identified several key challenges, the most pressing being the validation and implementation into actual clinical practice. CONCLUSIONS: The results confirm the potential of the Systems Medicine approach to study complex diseases and generate clinically relevant predictive models. Our study also highlights important obstacles and bottlenecks for such approaches (e.g. data availability and normalization of frameworks among others) and suggests specific proposals to overcome them.

Evaluating the role of sex-related structure-function differences on airway aerosol transport and deposition.

Several experimental studies have found that females have higher deposition of particles in the airways compared with males. This has implications for the delivery of aerosolized therapeutics and for understanding sex differences in respiratory system response to environmental exposures. This study evaluates several factors that potentially contribute to sex differences in particle deposition, using scale-specific structure-function models of 1D ventilation distribution, particle transport, and deposition. The impact of gravity, inhalation flow rate, and dead space are evaluated in 12 structure-based models (seven female; five male). Females were found to have significantly higher total, bronchial, and alveolar deposition than males across a particle size range from 0.01 to 10 . Results suggest that higher deposition fraction in females is due to higher alveolar deposition for smaller particle sizes, and higher bronchial deposition for larger particles. Females had higher alveolar deposition in the lower lobes, and slightly lower particle concentration in the left upper lobe. Males were found to be more sensitive to changes due to gravity, showing greater reduction in bronchial deposition fraction. Males were also more sensitive to change in inhalation flow rate, and to scaling of dead space due to the larger male baseline airway size. Predictions of sex differences in particle deposition - that are consistent with the literature - suggest that sex-based characteristics of lung and airway size interacting with particle size gives rise to differences in regional deposition.

A Wearable Open-Source electrical impedance tomography device.

Electrical impedance tomography (EIT) is medical imaging technique in which small electrical signals are used to map the electrical impedance distribution within the body. It is safe and non-invasive, which make it attractive for use in continuous monitoring or outpatient applications, but the high cost of commercial devices is an impediment to its adoption. Over the last 10 years, many research groups have developed their own EIT devices, but few designs for open-source EIT hardware are available. In this work, we present a complete open-source EIT system that is designed to be suitable for monitoring the lungs of free breathing subjects. The device is low-cost, wearable, and is designed to comply with the industry accepted safety standard for EIT. The device has been tested in two regimes: Firstly in terms of measurement uncertainty as a voltage measurement system, and secondly against a set of measures that have been proposed specifically for EIT hardware. The voltage measurement uncertainty of the device was measured to be - 0.7 % ± 0.36 mV. The EIT specific performance was measured in a phantom test designed to be as physiologically representative as practicable, and the device performed similarly to other published devices. This work will contribute to increased accessibility of EIT for study and will contribute to consensus on testing methodology for EIT devices.

Evaluating Tissue Heterogeneity in the Radiologically Normal-Appearing Tissue in IPF Compared to Healthy Controls.

RATIONALE AND OBJECTIVES: Idiopathic Pulmonary Fibrosis (IPF) is a progressive interstitial lung disease characterised by heterogeneously distributed fibrotic lesions. The inter- and intra-patient heterogeneity of the disease has meant that useful biomarkers of severity and progression have been elusive. Previous quantitative computed tomography (CT) based studies have focussed on characterising the pathological tissue. However, we hypothesised that the remaining lung tissue, which appears radiologically normal, may show important differences from controls in tissue characteristics. MATERIALS AND METHODS: Quantitative metrics were derived from CT scans in IPF patients (N = 20) and healthy controls with a similar age (N = 59). An automated quantitative software (CALIPER, Computer-Aided Lung Informatics for Pathology Evaluation and Rating) was used to classify tissue as normal-appearing, fibrosis, or low attenuation area. Densitometry metrics were calculated for all lung tissue and for only the normal-appearing tissue. Heterogeneity of lung tissue density was quantified as coefficient of variation and by quadtree. Associations between measured lung function and quantitative metrics were assessed and compared between the two cohorts. RESULTS: All metrics were significantly different between controls and IPF (p < 0.05), including when only the normal tissue was evaluated (p < 0.04). Density in the normal tissue was 14% higher in the IPF participants than controls (p < 0.001). The normal-appearing tissue in IPF had heterogeneity metrics that exhibited significant positive relationships with the percent predicted diffusion capacity for carbon monoxide. CONCLUSION: We provide quantitative assessment of IPF lung tissue characteristics compared to a healthy control group of similar age. Tissue that appears visually normal in IPF exhibits subtle but quantifiable differences that are associated with lung function and gas exchange.

Quantitative Analysis of Lung Shape in Idiopathic Pulmonary Fibrosis: Insights Into Disease- and Age-Associated Patterns.

RATIONALE AND OBJECTIVES: Fibrotic scarring in idiopathic pulmonary fibrosis (IPF) typically develops first in the posterior-basal lung tissue before advancing to involve more of the lung. The complexity of lung shape in the costo-diaphragmatic region has been proposed as a potential factor in this regional development. Intrinsic and disease-related shape could therefore be important for understanding IPF risk and its staging. We hypothesized that lung and lobe shape in IPF would have important differences from controls. MATERIALS AND METHODS: A principal component (PC) analysis was used to derive a statistical shape model (SSM) of the lung for a control cohort aged > 50 years (N = 39), using segmented lung and fissure surface data from CT imaging. Individual patient shape models derived for baseline (N = 18) and follow-up (N = 16) CT scans in patients with IPF were projected to the SSM to describe shape as the sum of the SSM average and weighted PC modes. Associations between the first four PC shape modes, lung function, percentage of fibrosis (fibrosis%) and pulmonary vessel-related structures (PVRS%), and other tissue metrics were assessed and compared between the two cohorts. RESULTS: Shape was different between IPF and controls (P < 0.05 for all shape modes), with IPF shape forming a distinct shape cluster. Shape had a negative relationship with age in controls (P = 0.013), but a positive relationship with age in IPF (P = 0.026). Some features of shape changed on follow-up. Shape in IPF was associated with fibrosis% (P < 0.05) and PVRS% (P < 0.05). CONCLUSION: Quantitative comparison of lung and lobe shape in IPF with controls of a similar age reveals shape differences that are strongly associated with age and percent fibrosis. The clustering of IPF cohort shape suggests that it could be an important feature to describe disease.

In silico prediction of e-cigarette aerosol particle transport and deposition within the airways.

Electronic cigarettes (ECs) generate aerosols by heating up a liquid ('e-liquid') that typically consists of propylene glycol (PG), vegetable glycerol (VG), nicotine and flavouring agents. These aerosols transport through the airway tree, and lung and deposit non-uniformly in the bronchi and alveoli. Studying the transport of aerosols through lung airways is necessary because it provides information about the concentration and deposition of particles in the upper and lower airways. Here, particle transport and deposition were simulated within an anatomically-realistic airway model, which was constructed from computed tomography imaging. Particle transport was simulated using the advection-diffusion equations. Particle deposition was estimated using three different mechanisms; including sedimentation, impaction and Brownian diffusion. Results show that by increasing the particle size (PS) from 50 nm to 500 nm, the total deposition efficiency decreased from 50% to 10%, and then by increasing the PS to 3 µm, it increased to 60%. In addition, Brownian deposition was the dominant mechanism for nanoparticles (PS≪0.5µm), while the sedimentation deposition mechanism was the dominant one for microparticles (PS≫0.5µm).Clinical relevance-There is an urgent need to understand the risk that ECs pose to human health and to determine the safest methods for using these devices to support smoking cessation whilst also minimising harm. The results of this study will be used to simulate the conditions such as aerosol concentration and flow rate in airways and alveoli to use in in vitro studies.

De novo discovery of traits co-occurring with chronic obstructive pulmonary disease.

Chronic obstructive pulmonary disease (COPD) is a heterogeneous group of chronic lung conditions. Genome-wide association studies have identified single-nucleotide polymorphisms (SNPs) associated with COPD and the co-occurring conditions, suggesting common biological mechanisms underlying COPD and these co-occurring conditions. To identify them, we have integrated information across different biological levels (i.e., genetic variants, lung-specific 3D genome structure, gene expression and protein-protein interactions) to build lung-specific gene regulatory and protein-protein interaction networks. We have queried these networks using disease-associated SNPs for COPD, unipolar depression and coronary artery disease. COPD-associated SNPs can control genes involved in the regulation of lung or pulmonary function, asthma, brain region volumes, cortical surface area, depressed affect, neuroticism, Parkinson's disease, white matter microstructure and smoking behaviour. We describe the regulatory connections, genes and biochemical pathways that underlay these co-occurring trait-SNP-gene associations. Collectively, our findings provide new avenues for the investigation of the underlying biology and diverse clinical presentations of COPD. In so doing, we identify a collection of genetic variants and genes that may aid COPD patient stratification and treatment.

Pulmonary vessel volume in idiopathic pulmonary fibrosis compared with healthy controls aged > 50 years.

Idiopathic pulmonary fibrosis (IPF) is characterised by progressive fibrosing interstitial pneumonia with an associated irreversible decline in lung function and quality of life. IPF prevalence increases with age, appearing most frequently in patients aged > 50 years. Pulmonary vessel-like volume (PVV) has been found to be an independent predictor of mortality in IPF and other interstitial lung diseases, however its estimation can be impacted by artefacts associated with image segmentation methods and can be confounded by adjacent fibrosis. This study compares PVV in IPF patients (N = 21) with PVV from a healthy cohort aged > 50 years (N = 59). The analysis includes a connected graph-based approach that aims to minimise artefacts contributing to calculation of PVV. We show that despite a relatively low extent of fibrosis in the IPF cohort (20% of the lung volume), PVV is 2-3 times higher than in controls. This suggests that a standardised method to calculate PVV that accounts for tree connectivity could provide a promising tool to provide early diagnostic or prognostic information in IPF patients and other interstitial lung disease.

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 evaluation of a New Zealand "vape to quit smoking" programme.

AIM: To compare the use of smoking cessation aids across different ethnic groups and age groups within a large New Zealand cohort and to assess the uptake and effectiveness of e-cigarettes for smoking cessation via a "vape to quit" initiative. METHODS: Retrospective analysis of Te Ha - Waitaha smoking cessation service, including a telephone interview of a subgroup, who opted into the "vape to quit" programme. The uptake of different smoking cessation aids, including the use of medications and other products, was evaluated and the self-reported quit rate in a "vape to quit" cohort was evaluated. RESULTS: The final dataset analysed consisted of 1,118 participants: 66.6% NZ European; 28.1% Maori; 3.1% Pacific; and 2.2% Asian. Maori participants were younger on average and had increasing vaping use. Maori were less likely to receive varenicline to assist with smoking cessation. Vaping use increased over time in all groups. Nicotine containing e-cigarettes were the most common smoking cessation products used, with >65% of each ethnic cohort utilising these products. Of the 100 participants in the "vape to quit" cohort 16% were smokefree and vapefree, 31% were smokefree and vaping, 31% were smoking and not vaping, and 22% were smoking and vaping. CONCLUSIONS: The Te Ha - Waitaha service was successful in engaging Maori in their smoking cessation programme. Nicotine containing e-cigarette products were popular in all cohorts. Nicotine containing e-cigarettes are showing potential in smoking cessation programmes in support of the Smokefree Aotearoa 2025; however, 22% of those in the "vape to quit" programme became dual users.

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/cm2 between postures.

Bridging the gap between respiratory research and health literacy: an interactive web-based platform.

Many patients with respiratory disease lack an understanding of basic respiratory physiology and the changes occurring in their lungs due to disease. Describing how the lungs work using realistic 3D visualisation of lung structure and function will improve communication of complicated concepts, resulting in improved health literacy. We developed a web-based platform, using anatomically realistic 3D lung models, to create an interactive visualisation tool to improve health literacy for patients with respiratory disease. A small amount of non-identifying personal information including gender, age, weight, height and smoking history can be used to customise the visualisation to an individual user. 3D computer modelling was used to create a web-based application that helps people understand how their lungs work in health and disease. The web-based application includes pages describing and visualising how the lungs work and the changes that occur during asthma and damage that smoking may be doing to their lungs. The application is freely available and located at https://sites.bioeng.auckland.ac.nz/silo6/lung_new/. This application bridges the gap between computational modelling and patient education, giving a visually compelling view into the patient's body that cannot be provided with any existing tools, hence providing a novel platform for enhancing patient-clinician interaction.

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.

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.

Human lungs are created to breathe clean air: the questionable quantification of vaping safety "95% less harmful".

The New Zealand government is aiming for Smokefree Aotorea, equivalent to a reduction in smoking prevalence to 5% or less by 2025. E-cigarettes may be one tool to meet this target, but how safe are they? Little is known about their long-term health implications in humans. In 2015, Public Health England commissioned a report summarising the available literature on e-cigarettes and coined the now well-known quantification that "e-cigarettes are 95% less harmful to your health than normal cigarettes". In this article, we argue that this is an unfounded quantification because the data required to make this quantification are not yet available. The value of '95% safer' was based on a study estimating the relative harms of nicotine-containing products that utilised scoring from a selected panel of experts. One of the key limitations of this quantification is that while the scores provided by the panellists were informed by knowledge, they are fundamentally value judgements and are not an exact science. E-cigarettes are probably safer than conventional cigarettes, however, there is mounting evidence that they are not without harm and the long-term health impacts are not yet known.

Integrated lung tissue mechanics one piece at a time: Computational modeling across the scales of biology.

The lung is a delicately balanced and highly integrated mechanical system. Lung tissue is continuously exposed to the environment via the air we breathe, making it susceptible to damage. As a consequence, respiratory diseases present a huge burden on society and their prevalence continues to rise. Emergent function is produced not only by the sum of the function of its individual components but also by the complex feedback and interactions occurring across the biological scales - from genes to proteins, cells, tissue and whole organ - and back again. Computational modeling provides the necessary framework for pulling apart and putting back together the pieces of the body and organ systems so that we can fully understand how they function in both health and disease. In this review, we discuss models of lung tissue mechanics spanning from the protein level (the extracellular matrix) through to the level of cells, tissue and whole organ, many of which have been developed in isolation. This is a vital step in the process but to understand the emergent behavior of the lung, we must work towards integrating these component parts and accounting for feedback across the scales, such as mechanotransduction. These interactions will be key to unlocking the mechanisms occurring in disease and in seeking new pharmacological targets and improving personalized healthcare.

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.