Investigation of the relative effects of vascular branching structure and gravity on pulmonary arterial blood flow heterogeneity via an image-based computational model

The age and gender-dependent changes in pulmonary arterial flow distribution: A cardiac magnetic resonance flow study

OBJECTIVE

To investigate whether the pulmonary artery blood flow distribution in the normal population changes depending on age and gender using the Cardiac Magnetic Resonance Phase Contrast Imaging method.

MATERIALS AND METHODS

A total of 62 healthy volunteers aged 20-60 years were included in the study between September 2020 and November 2021. Participants were first divided into two groups according to their age (group 120-40 years old, group 2: 40-60 years old). Then, each group was divided into two according to gender (groups 1A, 1B and groups 2A, 2B). After routine MR sequences of mediastinum and heart (True Fast Imaging with Steady State Free Precession, short axis for heart (SA), 4-chamber, 2-chamber and 3-chamber CINE sequence). Phase Contrast imaging was performed to Main Pulmonary Artery (MPA), right PA and left PA. On the phase contrast (PC) images obtained, flow velocity, forward flow, backward flow and net forward flow volume and gradient were calculated for each artery. With these data, the distribution ratio for both pulmonary arteries was calculated. Measurements were made simultaneously by two radiologists.

RESULTS

In phase contrast sequence evaluation; No statistically significant correlation was found between RPA/LPA distribution ratio and Group 1 (p > 0.05). A moderate and positive correlation was found statistically between RPA/LPA distribution ratio and Group 2 (r:0.42, p:0.04). No statistically significant correlation was found between the RPA/LPA distribution ratio and the gender variable (p:0.59).

CONCLUSION

The distribution ratio of pulmonary blood flow to both lungs in healthy volunteers can be measured using PC-CMR without the need for invasive procedures. Knowing the normal distribution ratios according to age in the population can be used for the early diagnosis of pathological conditions.

Pathophysiology and pathogenic mechanisms of pulmonary hypertension: role of membrane receptors, ion channels, and Ca2+ signaling

The pulmonary circulation is a low-resistance, low-pressure, and high-compliance system that allows the lungs to receive the entire cardiac output. Pulmonary arterial pressure is a function of cardiac output and pulmonary vascular resistance, and pulmonary vascular resistance is inversely proportional to the fourth power of the intraluminal radius of the pulmonary artery. Therefore, a very small decrease of the pulmonary vascular lumen diameter results in a significant increase in pulmonary vascular resistance and pulmonary arterial pressure. Pulmonary arterial hypertension is a fatal and progressive disease with poor prognosis. Regardless of the initial pathogenic triggers, sustained pulmonary vasoconstriction, concentric vascular remodeling, occlusive intimal lesions, in situ thrombosis, and vascular wall stiffening are the major and direct causes for elevated pulmonary vascular resistance in patients with pulmonary arterial hypertension and other forms of precapillary pulmonary hypertension. In this review, we aim to discuss the basic principles and physiological mechanisms involved in the regulation of lung vascular hemodynamics and pulmonary vascular function, the changes in the pulmonary vasculature that contribute to the increased vascular resistance and arterial pressure, and the pathogenic mechanisms involved in the development and progression of pulmonary hypertension. We focus on reviewing the pathogenic roles of membrane receptors, ion channels, and intracellular Ca2+ signaling in pulmonary vascular smooth muscle cells in the development and progression of pulmonary hypertension.

Lobar pulmonary perfusion quantification with dual-energy CT angiography: Interlobar variability and relationship with regional clot burden in pulmonary embolism

Purpose

Semi-automated lobar segmentation tools enable an anatomical assessment of regional pulmonary perfusion with Dual-Energy CTA (DE-CTA). We aimed to quantify lobar pulmonary perfusion with DE-CTA, analyze the perfusion distribution among the pulmonary lobes in subjects without cardiopulmonary diseases and assess the correlation between lobar perfusion and regional endoluminal clots in patients with acute pulmonary embolism (PE).

Methods

We evaluated 151 consecutive subjects with suspected PE and without cardiopulmonary comorbidities. DE-CTA derived perfused blood volume (PBV) of each pulmonary lobe was measured applying a semi-automated lobar segmentation technique. In patients with PE, blood clot location was assessed, and CT-based vascular obstruction index of each lobe (CTOI_lobe) was calculated and classified into three groups: CTOI_lobe= 0, low CTOI_lobe (1–50%) and high CTOI_lobe (>50%).

Results

Among patients without PE (103/151, 68.2%), median lobar PBV was 13.7% (IQR 10.2–18.0%); the right middle lobe presented lower PBV when compared to all the other lobes (p < .001). In patients with PE (48/151, 31.8%), lobar PBV was 12.6% (IQR 9.6–15.7%), 13.7% (IQR 10.1–16.7%) and 6.5% (IQR 5.1–10.2%) in the lobes with CTOI_lobe= 0, low CTOI_lobe and high CTOI_lobe scores, respectively, with a significantly decreased PBV in the lobes with high CTOI_lobe score (p < .001). ROC analysis of lobar PBV for prediction of high CTOI_lobe score revealed AUC of 0.847 (95%CI 0.785–0.908).

Conclusion

Pulmonary perfusion was heterogeneously distributed along the pulmonary lobes in patients without cardiopulmonary diseases. In patients with PE, the lobes with high vascular obstruction score (CTOI_lobe> 50%) presented a decreased lobar perfusion.

•

Semi-automated tools enable assessment of lobar perfusion with Dual-Energy CTA.

•

The pulmonary perfusion is heterogeneously distributed along the pulmonary lobes.

•

Lobar perfusion was decreased only in the lobes with high vascular obstruction index.

Perfusion imaging heterogeneity during NO inhalation distinguishes pulmonary arterial hypertension (PAH) from healthy subjects and has potential as an imaging biomarker

BACKGROUND

Without aggressive treatment, pulmonary arterial hypertension (PAH) has a 5-year mortality of approximately 40%. A patient's response to vasodilators at diagnosis impacts the therapeutic options and prognosis. We hypothesized that analyzing perfusion images acquired before and during vasodilation could identify characteristic differences between PAH and control subjects.

METHODS

We studied 5 controls and 4 subjects with PAH using HRCT and ¹³NN PET imaging of pulmonary perfusion and ventilation. The total spatial heterogeneity of perfusion (CV²_Qtotal) and its components in the vertical (CV²_Qvgrad) and cranio-caudal (CV²_Qzgrad) directions, and the residual heterogeneity (CV²_Qr), were assessed at baseline and while breathing oxygen and nitric oxide (O₂ + iNO). The length scale spectrum of CV²_Qr was determined from 10 to 110 mm, and the response of regional perfusion to O₂ + iNO was calculated as the mean of absolute differences. Vertical gradients in perfusion (Q_vgrad) were derived from perfusion images, and ventilation-perfusion distributions from images of ¹³NN washout kinetics.

RESULTS

O₂ + iNO significantly enhanced perfusion distribution differences between PAH and controls, allowing differentiation of PAH subjects from controls. During O₂ + iNO, CV²_Qvgrad was significantly higher in controls than in PAH (0.08 (0.055-0.10) vs. 6.7 × 10^-3 (2 × 10^-4-0.02), p < 0.001) with a considerable gap between groups. Q_vgrad and CV²_Qtotal showed smaller differences: - 7.3 vs. - 2.5, p = 0.002, and 0.12 vs. 0.06, p = 0.01. CV²_Qvgrad had the largest effect size among the primary parameters during O₂ + iNO. CV²_Qr, and its length scale spectrum were similar in PAH and controls. Ventilation-perfusion distributions showed a trend towards a difference between PAH and controls at baseline, but it was not statistically significant.

CONCLUSIONS

Perfusion imaging during O2 + iNO showed a significant difference in the heterogeneity associated with the vertical gradient in perfusion, distinguishing in this small cohort study PAH subjects from controls.

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.

Abnormal pulmonary perfusion heterogeneity in patients with Fontan circulation and pulmonary arterial hypertension

KEY POINTS

The distribution of pulmonary perfusion is affected by gravity, vascular branching structure and active regulatory mechanisms, which may be disrupted by cardiopulmonary disease, but this is not well studied, particularly in rare conditions. We evaluated pulmonary perfusion in patients who had undergone Fontan procedure, patients with pulmonary arterial hypertension (PAH) and two groups of controls using a proton magnetic resonance imaging technique, arterial spin labelling to measure perfusion. Heterogeneity was assessed by the relative dispersion (SD/mean) and gravitational gradients. Gravitational gradients were similar between all groups, but heterogeneity was significantly increased in both patient groups compared to controls and persisted after removing contributions from large blood vessels and gravitational gradients. Patients with Fontan physiology and patients with PAH have increased pulmonary perfusion heterogeneity that is not explainable by differences in mean perfusion, gravitational gradients, or large vessel anatomy. This probably reflects vascular remodelling in PAH and possibly in Fontan physiology.

ABSTRACT

Many factors affect the distribution of pulmonary perfusion, which may be disrupted by cardiopulmonary disease, but this is not well studied, particularly in rare conditions. An example is following the Fontan procedure, where pulmonary perfusion is passive, and heterogeneity may be increased because of the underlying pathophysiology leading to Fontan palliation, remodelling, or increased gravitational gradients from low flow. Another is pulmonary arterial hypertension (PAH), where gravitational gradients may be reduced secondary to high pressures, but remodelling may increase perfusion heterogeneity. We evaluated regional pulmonary perfusion in Fontan patients (n = 5), healthy young controls (Fontan control, n = 5), patients with PAH (n = 6) and healthy older controls (PAH control) using proton magnetic resonance imaging. Regional perfusion was measured using arterial spin labelling. Heterogeneity was assessed by the relative dispersion (SD/mean) and gravitational gradients. Mean perfusion was similar (Fontan = 2.50 ± 1.02 ml min^-1  ml^-1 ; Fontan control = 3.09 ± 0.58, PAH = 3.63 ± 1.95; PAH control = 3.98 ± 0.91, P = 0.26), and the slopes of gravitational gradients were not different (Fontan = -0.23 ± 0.09 ml min^-1  ml^-1  cm^-1 ; Fontan control = -0.29 ± 0.23, PAH = -0.27 ± 0.09, PAH control = -0.25 ± 0.18, P = 0.91) between groups. Perfusion relative dispersion was greater in both Fontan and PAH than controls (Fontan = 1.46 ± 0.18; Fontan control = 0.99 ± 0.21, P = 0.005; PAH = 1.22 ± 0.27, PAH control = 0.91 ± 0.12, P = 0.02) but similar between patient groups (P = 0.13). These findings persisted after removing contributions from large blood vessels and gravitational gradients (all P < 0.05). We conclude that patients with Fontan physiology and PAH have increased pulmonary perfusion heterogeneity that is not explained by differences in mean perfusion, gravitational gradients, or large vessel anatomy. This probably reflects the effects of remodelling in PAH and possibly in Fontan physiology.

Invasive Pulmonary Adenocarcinomas Versus Preinvasive Lesions Appearing as Pure Ground-Glass Nodules: Differentiation Using Enhanced Dual-Source Dual-Energy CT

OBJECTIVE. The objective of our study was to investigate the potentials of enhanced dual-source dual-energy CT (DECT) and three-planar measurements for differentiating invasive pulmonary adenocarcinomas (IPAs) from preinvasive lesions appearing as pure ground-glass nodules (pGGNs). MATERIALS AND METHODS. Thirty-nine patients with 53 pGGNs who underwent enhanced dual-source DECT were included in this retrospective study. All pGGNs were pathologically confirmed and categorized into two groups: preinvasive lesions or IPAs. The traditional CT features of the pGGNs were evaluated on unenhanced images. Quantitative parameters were measured on iodine-enhanced images of dual-source DECT in three planes, and both intra- and interobserver reproducibility analyses were performed to assess the measurement reproducibility of quantitative parameters. To identify significant factors for differentiating IPAs from preinvasive lesions, we performed logistic regression analysis and ROC curve analysis. RESULTS. For traditional CT features, only lesion size and unenhanced CT attenuation value showed significant differences between preinvasive lesions and IPAs (p < 0.05). Preinvasive lesions and IPAs exhibited significant differences in attenuation on virtual images, so-called "virtual HU" or "VHU," and the modified normalized iodine concentration (NIC) (p < 0.05), and both intra- and interobserver agreement for the quantitative measurements were excellent. Multivariate logistic regression analysis revealed that larger lesion size (adjusted odds ratio [OR], 3.65) and higher modified NIC (adjusted OR, 19.01) were significant differentiators of IPAs from preinvasive lesions (p < 0.05). ROC curve analysis revealed that modified NIC showed excellent performance (AUC, 0.924) and significantly higher performance than lesion size (AUC, 0.711) for differentiating IPAs from preinvasive lesions. CONCLUSION. In pGGNs, a lesion with a modified NIC value of more than 0.29 can be a very specific discriminator of IPAs from preinvasive lesions, and IPAs can be accurately and reliably differentiated from preinvasive lesions using enhanced dual-source DECT and three-planar measurements.

Modeling Lung Architecture in the XCAT Series of Phantoms: Physiologically Based Airways, Arteries and Veins

The purpose of this paper was to extend the extended cardiac-torso (XCAT) series of computational phantoms to include a detailed lung architecture including airways and pulmonary vasculature. Eleven XCAT phantoms of varying anatomy were used in this paper. The lung lobes and initial branches of the airways, pulmonary arteries, and veins were previously defined in each XCAT model. These models were extended from the initial branches of the airways and vessels to the level of terminal branches using an anatomically-based volume-filling branching algorithm. This algorithm grew the airway and vasculature branches separately and iteratively without intersecting each other using cylindrical models with diameters estimated by order-based anatomical measurements. Geometrical features of the extended branches were compared with the literature anatomy values to quantitatively evaluate the models. These features include branching angle, length to diameter ratio, daughter to parent diameter ratio, asymmetrical branching pattern, diameter, and length ratios. The XCAT phantoms were then used to simulate CT images to qualitatively compare them with the original phantom images. The proposed growth model produced 46369 ± 12521 airways, 44737 ± 11773 arteries, and 39819 ± 9988 veins to the XCAT phantoms. Furthermore, the growth model was shown to produce asymmetrical airway, artery, and vein networks with geometrical attributes close to morphometry and model based studies. The simulated CT images of the phantoms were judged to be more realistic, including more airways and pulmonary vessels compared with the original phantoms. Future work will seek to add a heterogeneous parenchymal background into the XCAT lungs to make the phantoms even more representative of human anatomy, paving the way towards the use of XCAT models as a tool to virtually evaluate the current and emerging medical imaging technologies.

Regional Distribution of Pulmonary Blood Volume with Dual-Energy Computed Tomography: Results in 42 Subjects

RATIONALE AND OBJECTIVES

The noninvasive approach of lung perfusion generated from dual-energy computed tomography acquisitions has entered clinical practice. The purpose of this study was to analyze the regional distribution of iodine within distal portions of the pulmonary arterial bed on dual-source, dual-energy computed tomography examinations in a cohort of subjects without cardiopulmonary pathologies.

MATERIALS AND METHODS

The study population included 42 patients without cardiorespiratory disease, enabling quantitative and qualitative analysis of pulmonary blood volume after administration of a 40% contrast agent. Qualitative analysis was based on visual assessment. Quantitative analysis was obtained after semiautomatic division of each lung into 18 areas.

RESULTS

The iodine concentration did not significantly differ between the right (R) and left (L) lungs (P = .49), with a mean attenuation of 41.35 Hounsfield units (HU) and 41.14 HU, respectively. Three regional gradients of attenuation were observed between: (a) lung bases and apices (P < .001), linked to the conditions of examination (mean Δ: 6.23 in the R lung; 5.96 in the L lung); (b) posterior and anterior parts of the lung (P < .001) due to gravity (mean Δ: 11.92 in the R lung ; 15.93 in the L lung); and (c) medullary and cortical lung zones (P < .001) (mean Δ: 9.35 in the R lung ; 8.37 in the L lung). The intensity of dependent-nondependent (r = 0.42; P < .001) and corticomedullary (r = 0.58; P < .0001) gradients was correlated to the overall iodine concentration.

CONCLUSION

Distribution of pulmonary blood volume is influenced by physiological gradients and scanning conditions.

Nanoparticle transport and delivery in a heterogeneous pulmonary vasculature

Quantitative understanding of nanoparticles delivery in a complex vascular networks is very challenging because it involves interplay of transport, hydrodynamic force, and multivalent interactions across different scales. Heterogeneous pulmonary network includes up to 16 generations of vessels in its arterial tree. Modeling the complete pulmonary vascular system in 3D is computationally unrealistic. To save computational cost, a model reconstructed from MRI scanned images is cut into an arbitrary pathway consisting of the upper 4-generations. The remaining generations are represented by an artificially rebuilt pathway. Physiological data such as branch information and connectivity matrix are used for geometry reconstruction. A lumped model is used to model the flow resistance of the branches that are cut off from the truncated pathway. Moreover, since the nanoparticle binding process is stochastic in nature, a binding probability function is used to simplify the carrier attachment and detachment processes. The stitched realistic and artificial geometries coupled with the lumped model at the unresolved outlets are used to resolve the flow field within the truncated arterial tree. Then, the biodistribution of 200nm, 700nm and 2µm particles at different vessel generations is studied. At the end, 0.2-0.5% nanocarrier deposition is predicted during one time passage of drug carriers through pulmonary vascular tree. Our truncated approach enabled us to efficiently model hemodynamics and accordingly particle distribution in a complex 3D vasculature providing a simple, yet efficient predictive tool to study drug delivery at organ level.

Vertical gradients in regional alveolar oxygen tension in supine human lung imaged by hyperpolarized 3He MRI

The purpose of this study was to evaluate whether regional alveolar oxygen tension (P(A)O2) vertical gradients imaged with hyperpolarized (3)He can identify smoking-induced pulmonary alterations. These gradients are compared with common clinical measurements including pulmonary function tests (PFTs), the six minute walk test, and the St. George's Respiratory Questionnaire. 8 healthy non-smokers, 12 asymptomatic smokers, and 7 symptomatic subjects with chronic obstructive pulmonary disease (COPD) underwent two sets of back-to-back P(A)O2 imaging acquisitions in the supine position in two opposite directions (top to bottom and bottom to top), followed by clinically standard pulmonary tests. The whole-lung mean, standard deviation (DP(A)O2) and vertical gradients of P(A)O2 along the slices were extracted, and the results were compared with clinically derived metrics. Statistical tests were performed to analyze the differences between cohorts. The anterior-posterior vertical gradients and DP(A)O2 effectively differentiated all three cohorts (p < 0.05). The average vertical gradient P(A)O2 in healthy subjects was -1.03 ± 0.51 Torr/cm toward lower values in the posterior/dependent regions. The directional gradient was absent in smokers (0.36 ± 1.22 Torr/cm) and was in the opposite direction in COPD subjects (2.18 ± 1.54 Torr/cm). The vertical gradients correlated with smoking history (p = 0.004); body mass index (p = 0.037), PFT metrics (forced expiratory volume in 1 s, p = 0.025; residual volume/total lung capacity percent predicted, p = 0.033) and with distance walked in 6 min (p = 0.009). Regional P(A)O2 data indicate that cigarette smoke induces physiological alterations that are not being detected by the most widely used physiological tests.

Endothelial and smooth muscle cell ion channels in pulmonary vasoconstriction and vascular remodeling

The pulmonary circulation is a low resistance and low pressure system. Sustained pulmonary vasoconstriction and excessive vascular remodeling often occur under pathophysiological conditions such as in patients with pulmonary hypertension. Pulmonary vasoconstriction is a consequence of smooth muscle contraction. Many factors released from the endothelium contribute to regulating pulmonary vascular tone, while the extracellular matrix in the adventitia is the major determinant of vascular wall compliance. Pulmonary vascular remodeling is characterized by adventitial and medial hypertrophy due to fibroblast and smooth muscle cell proliferation, neointimal proliferation, intimal, and plexiform lesions that obliterate the lumen, muscularization of precapillary arterioles, and in situ thrombosis. A rise in cytosolic free Ca(2+) concentration ([Ca(2+)]cyt) in pulmonary artery smooth muscle cells (PASMC) is a major trigger for pulmonary vasoconstriction, while increased release of mitogenic factors, upregulation (or downregulation) of ion channels and transporters, and abnormalities in intracellular signaling cascades are key to the remodeling of the pulmonary vasculature. Changes in the expression, function, and regulation of ion channels in PASMC and pulmonary arterial endothelial cells play an important role in the regulation of vascular tone and development of vascular remodeling. This article will focus on describing the ion channels and transporters that are involved in the regulation of pulmonary vascular function and structure and illustrating the potential pathogenic role of ion channels and transporters in the development of pulmonary vascular disease.

Determinants of pulmonary blood flow distribution

The primary function of the pulmonary circulation is to deliver blood to the alveolar capillaries to exchange gases. Distributing blood over a vast surface area facilitates gas exchange, yet the pulmonary vascular tree must be constrained to fit within the thoracic cavity. In addition, pressures must remain low within the circulatory system to protect the thin alveolar capillary membranes that allow efficient gas exchange. The pulmonary circulation is engineered for these unique requirements and in turn these special attributes affect the spatial distribution of blood flow. As the largest organ in the body, the physical characteristics of the lung vary regionally, influencing the spatial distribution on large-, moderate-, and small-scale levels.

Spatial distribution of ventilation and perfusion: mechanisms and regulation

With increasing spatial resolution of regional ventilation and perfusion, it has become more apparent that ventilation and blood flow are quite heterogeneous in the lung. A number of mechanisms contribute to this regional variability, including hydrostatic gradients, pleural pressure gradients, lung compressibility, and the geometry of the airway and vascular trees. Despite this marked heterogeneity in both ventilation and perfusion, efficient gas exchange is possible through the close regional matching of the two. Passive mechanisms, such as the shared effect of gravity and the matched branching of vascular and airway trees, create efficient gas exchange through the strong correlation between ventilation and perfusion. Active mechanisms that match local ventilation and perfusion play little if no role in the normal healthy lung but are important under pathologic conditions.

Distribution of perfusion

Local driving pressures and resistances within the pulmonary vascular tree determine the distribution of perfusion in the lung. Unlike other organs, these local determinants are significantly influenced by regional hydrostatic and alveolar pressures. Those effects on blood flow distribution are further magnified by the large vertical height of the human lung and the relatively low intravascular pressures in the pulmonary circulation. While the distribution of perfusion is largely due to passive determinants such as vascular geometry and hydrostatic pressures, active mechanisms such as vasoconstriction induced by local hypoxia can also redistribute blood flow. This chapter reviews the determinants of regional lung perfusion with a focus on vascular tree geometry, vertical gradients induced by gravity, the interactions between vascular and surrounding alveolar pressures, and hypoxic pulmonary vasoconstriction. While each of these determinants of perfusion distribution can be examined in isolation, the distribution of blood flow is dynamically determined and each component interacts with the others so that a change in one region of the lung influences the distribution of blood flow in other lung regions.

Imaging for lung physiology: what do we wish we could measure?

The role of imaging as a tool for investigating lung physiology is growing at an accelerating pace. Looking forward, we wished to identify unresolved issues in lung physiology that might realistically be addressed by imaging methods in development or imaging approaches that could be considered. The role of imaging is framed in terms of the importance of good spatial and temporal resolution and the types of questions that could be addressed as these technical capabilities improve. Recognizing that physiology is fundamentally a quantitative science, a recurring emphasis is on the need for imaging methods that provide reliable measurements of specific physiological parameters. The topics included necessarily reflect our perspective on what are interesting questions and are not meant to be a comprehensive review. Nevertheless, we hope that this essay will be a spur to physiologists to think about how imaging could usefully be applied in their research and to physical scientists developing new imaging methods to attack challenging questions imaging could potentially answer.

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.

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.

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

Recent experimental and imaging studies suggest that the influence of gravity on the measured distribution of blood flow in the lung is largely through deformation of the parenchymal tissue. To study the contribution of hydrostatic effects to regional perfusion in the presence of tissue deformation, we have developed an anatomically structured computational model of the pulmonary circulation (arteries, capillaries, veins), coupled to a continuum model of tissue deformation, and including scale-appropriate fluid dynamics for blood flow in each vessel type. The model demonstrates that both structural and the multiple effects of gravity on the pulmonary circulation make a distinct contribution to the distribution of blood. It shows that postural differences in perfusion gradients can be explained by the combined effect of tissue deformation and extra-acinar blood vessel resistance to flow in the dependent tissue. However, gravitational perfusion gradients persist when the effect of tissue deformation is eliminated, highlighting the importance of the hydrostatic effects of gravity on blood distribution in the pulmonary circulation. Coupling of large- and small-scale models reveals variation in microcirculatory driving pressures within isogravitational planes due to extra-acinar vessel resistance. Variation in driving pressures is due to heterogeneous large-vessel resistance as a consequence of geometric asymmetry in the vascular trees and is amplified by the complex balance of pressures, distension, and flow at the microcirculatory level.

Image-based modeling of lung structure and function

The current state-of-the-art in image-based modeling allows derivation of patient-specific models of the lung, lobes, airways, and pulmonary vascular trees. The application of traditional engineering analyses of fluid and structural mechanics to image-based subject-specific models has the potential to provide new insight into structure-function relationships in the individual via functional interpretation that complements imaging and experimental studies. Three major issues that are encountered in studies of airflow through the bronchial airways are the representation of airway geometry, the imposition of physiological boundary conditions, and the treatment of turbulence. Here we review some efforts to resolve each of these issues, with particular focus on image-based models that have been developed to simulate airflow from the mouth to the terminal bronchiole, and subjected to physiologically meaningful boundary conditions via image registration and soft-tissue mechanics models.

Computational Simulation of the Pulmonary Arteries and its Role in the Study of Pediatric Pulmonary Hypertension

The hemodynamic state of the pulmonary arteries is challenging to routinely measure in children due to the vascular circuit's position in the lungs. The resulting relative scarcity of quantitative clinical diagnostic and prognostic information impairs management of diseases such as pulmonary hypertension, or high blood pressure of the pulmonary circuit, and invites new techniques of measurement. Here we examine recent applications of macro-scale computational mechanics methods for fluids and solids - traditionally used by engineers in the design and virtual testing of complex metal and composite structures - applied to study the pulmonary vasculature, both in healthy and diseased states. In four subject areas, we briefly outline advances in computational methodology and provide examples of clinical relevance.

Lung volume does not alter the distribution of pulmonary perfusion in dependent lung in supine humans

There is a gravitational influence on pulmonary perfusion, including in the most dependent lung, where perfusion is reduced, termed Zone 4. Studies using xenon-133 show Zone 4 behaviour, present in the dependent 4 cm at total lung capacity (TLC), affects the dependent 11 cm at functional residual capacity (FRC) and almost all the lung at residual volume (RV). These differences were ascribed to increased resistance in extra-alveolar vessels at low lung volumes although other mechanisms have been proposed. To further evaluate the behaviour of perfusion in dependent lung using a technique that directly measures pulmonary perfusion and corrects for tissue distribution by measuring regional proton density, seven healthy subjects (age = 38 ± 6 years, FEV₁ = 104 ± 7% predicted) underwent magnetic resonance imaging in supine posture. Data were acquired in the right lung during breath-holds at RV, FRC and TLC. Arterial spin labelling quantified regional pulmonary perfusion, which was normalized for regional proton density measured using a fast low-angle shot technique. The height of the onset of Zone 4 behaviour was not different between lung volumes (P = 0.23). There were no significant differences in perfusion (expressed as ml min⁻¹ g⁻¹) between lung volumes in the gravitationally intermediate (RV = 8.9 ± 3.1, FRC = 8.1 ± 2.9, TLC = 7.4 ± 3.6; P = 0.26) and dependent lung (RV = 6.6 ± 2.4, FRC = 6.1 ± 2.1, TLC = 6.4 ± 2.6; P = 0.51). However, at TLC perfusion was significantly lower in non-dependent lung than at FRC or RV (3.6 ± 3.3, 7.7 ± 1.5, 7.9 ± 2.0, respectively; P < 0.001). These data suggest that the mechanism of the reduction in perfusion in dependent lung is unlikely to be a result of lung volume related increases in resistance in extra-alveolar vessels. In supine posture, the gravitational influence on perfusion is remarkably similar over most of the lung, irrespective of lung volume.

Species-specific pulmonary arterial asymmetry determines species differences in regional pulmonary perfusion

The functional significance of differences in pulmonary vascular branching and diameter asymmetry between the human and quadruped lung has not previously been addressed. To evaluate the contribution of branching asymmetry to observable species differences in blood flow gradients, computed distributions of blood flow were compared in structure-based models of the human and ovine pulmonary arteries. The models were derived using a combination of computed tomography and a volume-filling algorithm. Pressure, flow, and deformed vessel diameter were calculated in both species models using equations representing conservation of mass and momentum, and a pressure-diameter relationship. The major difference between the human and ovine results was the presence of a large region of "zone 4" flow and higher mean flows in the central region of the ovine lung compared to that in the human. Heterogeneity in tissue perfusion and the contribution of gravity were similar in both species models; however, the gravitationally directed gradients of perfusion in the human and ovine models were different and each consistent with human and quadruped measurements, respectively. The results suggest that measured species differences in pulmonary perfusion gradients are largely determined by differences in branching asymmetry.

Modelling pulmonary blood flow

Computational model analysis has been used widely to understand and interpret complexity of interactions in the pulmonary system. Pulmonary blood transport is a multi-scale phenomenon that involves scale-dependent structure and function, therefore requiring different model assumptions for the microcirculation and the arterial or venous flows. The blood transport systems interact with the surrounding lung tissue, and are dependent on hydrostatic pressure gradients, control of vasoconstriction, and the topology and material composition of the vascular trees. This review focuses on computational models that have been developed to study the different mechanisms contributing to regional perfusion of the lung. Different models for the microcirculation and the pulmonary arteries are considered, including fractal approaches and anatomically-based methods. The studies that are reviewed illustrate the different complementary approaches that can be used to address the same physiological question of flow heterogeneity.

Computer simulation analysis of normal and abnormal development of the mammalian diaphragm

BACKGROUND

Congenital diaphragmatic hernia (CDH) is a birth defect with significant morbidity and mortality. Knowledge of diaphragm morphogenesis and the aberrations leading to CDH is limited. Although classical embryologists described the diaphragm as arising from the septum transversum, pleuroperitoneal folds (PPF), esophageal mesentery and body wall, animal studies suggest that the PPF is the major, if not sole, contributor to the muscular diaphragm. Recently, a posterior defect in the PPF has been identified when the teratogen nitrofen is used to induce CDH in fetal rodents. We describe use of a cell-based computer modeling system (Nudge++) to study diaphragm morphogenesis.

METHODS AND RESULTS

Key diaphragmatic structures were digitized from transverse serial sections of paraffin-embedded mouse embryos at embryonic days 11.5 and 13. Structure boundaries and simulated cells were combined in the Nudge++ software. Model cells were assigned putative behavioral programs, and these programs were progressively modified to produce a diaphragm consistent with the observed anatomy in rodents. Homology between our model and recent anatomical observations occurred under the following simulation conditions: (1) cell mitoses are restricted to the edge of growing tissue; (2) cells near the chest wall remain mitotically active; (3) mitotically active non-edge cells migrate toward the chest wall; and (4) movement direction depends on clonal differentiation between anterior and posterior PPF cells.

CONCLUSION

With the PPF as the sole source of mitotic cells, an early defect in the PPF evolves into a posteromedial diaphragm defect, similar to that of the rodent nitrofen CDH model. A posterolateral defect, as occurs in human CDH, would be more readily recreated by invoking other cellular contributions. Our results suggest that recent reports of PPF-dominated diaphragm morphogenesis in the rodent may not be strictly applicable to man. The ability to recreate a CDH defect using a combination of experimental data and testable hypotheses gives impetus to simulation modeling as an adjunct to experimental analysis of diaphragm morphogenesis.

Computational models of structure-function relationships in the pulmonary circulation and their validation

The pulmonary airway, arterial, venous and capillary networks are vast complex branching and converging systems that are mechanically coupled to the surrounding lung tissue. Early studies that examined vascular or airway geometry relied on measurements from casts, but medical imaging now enables measurement of the lung in vivo, at controlled lung volumes. The high-quality data that imaging provides have prompted development of increasingly sophisticated models of the geometry of the airway and pulmonary vascular trees. The accurate spatial relationships between airway, vessel and tissue in these imaging-derived models are necessary for computational analysis that aims to elucidate regional airway-vessel-tissue interactions. Predictions of blood flow through multiscale imaging-derived models of the pulmonary arteries and capillary bed reveal geometry-dependent patterns of perfusion in response to gravity and lung orientation that cannot be predicted with simplified, summary representations of the pulmonary transport trees. Validation of such predictions against measures from functional imaging holds significant potential for explaining and differentiating normal and disease-related heterogeneity in regional blood flow calculated using perfusion imaging.

Computational predictions of pulmonary blood flow gradients: gravity versus structure

A computational model of blood flow through the human pulmonary arterial tree has been developed to investigate the mechanisms contributing to regional pulmonary perfusion in the isolated network when the lung is in different orientations. The arterial geometric model was constructed using a combination of computed tomography and a volume-filling branching algorithm. Equations governing conservation of mass, momentum, and vessel distension, incorporating gravity, were solved to predict pressure, flow, and vessel radius. Analysis of results in the upright posture, with and without gravity, and in the inverted, prone, and supine postures reveals significant flow heterogeneity and a persistent decrease in flow in the cranial and caudal regions for all postures suggesting that vascular geometry makes a major contribution to regional flow with gravity having a lesser role. Results in the isolated arterial tree demonstrate that the vascular path lengths and therefore the positioning of the pulmonary trunk relative to the rest of the network play a significant role in the determination of flow.