Heterogeneous mechanics of the mouse pulmonary arterial network

A computational growth and remodeling framework for adaptive and maladaptive pulmonary arterial hemodynamics

Hemodynamic loading is known to contribute to the development and progression of pulmonary arterial hypertension (PAH). This loading drives changes in mechanobiological stimuli that affect cellular phenotypes and lead to pulmonary vascular remodeling. Computational models have been used to simulate mechanobiological metrics of interest, such as wall shear stress, at single time points for PAH patients. However, there is a need for new approaches that simulate disease evolution to allow for prediction of long-term outcomes. In this work, we develop a framework that models the pulmonary arterial tree through adaptive and maladaptive responses to mechanical and biological perturbations. We coupled a constrained mixture theory-based growth and remodeling framework for the vessel wall with a morphometric tree representation of the pulmonary arterial vasculature. We show that non-uniform mechanical behavior is important to establish the homeostatic state of the pulmonary arterial tree, and that hemodynamic feedback is essential for simulating disease time courses. We also employed a series of maladaptive constitutive models, such as smooth muscle hyperproliferation and stiffening, to identify critical contributors to development of PAH phenotypes. Together, these simulations demonstrate an important step toward predicting changes in metrics of clinical interest for PAH patients and simulating potential treatment approaches.

Sensitivity analysis and uncertainty quantification of 1-D models of pulmonary hemodynamics in mice under control and hypertensive conditions

Pulmonary hypertension (PH), defined as an elevated mean blood pressure in the main pulmonary artery (MPA) at rest, is associated with vascular remodeling of both large and small arteries. PH has several sub-types that are all linked to high mortality rates. In this study, we use a one-dimensional (1-D) fluid dynamics model driven by in vivo measurements of MPA flow to understand how model parameters and network size influence MPA pressure predictions in the presence of PH. We compare model predictions with in vivo MPA pressure measurements from a control and a hypertensive mouse and analyze results in three networks of increasing complexity, extracted from micro-computed tomography (micro-CT) images. We introduce global scaling factors for boundary condition parameters and perform local and global sensitivity analysis to calculate parameter influence on model predictions of MPA pressure and correlation analysis to determine a subset of identifiable parameters. These are inferred using frequentist optimization and Bayesian inference via the Delayed Rejection Adaptive Metropolis (DRAM) algorithm. Frequentist and Bayesian uncertainty is computed for model parameters and MPA pressure predictions. Results show that MPA pressure predictions are most sensitive to distal vascular resistance and that parameter influence changes with increasing network complexity. Our outcomes suggest that PH leads to increased vascular stiffness and decreased peripheral compliance, congruent with clinical observations.

Assessing model mismatch and model selection in a Bayesian uncertainty quantification analysis of a fluid-dynamics model of pulmonary blood circulation

This study uses Bayesian inference to quantify the uncertainty of model parameters and haemodynamic predictions in a one-dimensional pulmonary circulation model based on an integration of mouse haemodynamic and micro-computed tomography imaging data. We emphasize an often neglected, though important source of uncertainty: in the mathematical model form due to the discrepancy between the model and the reality, and in the measurements due to the wrong noise model (jointly called 'model mismatch'). We demonstrate that minimizing the mean squared error between the measured and the predicted data (the conventional method) in the presence of model mismatch leads to biased and overly confident parameter estimates and haemodynamic predictions. We show that our proposed method allowing for model mismatch, which we represent with Gaussian processes, corrects the bias. Additionally, we compare a linear and a nonlinear wall model, as well as models with different vessel stiffness relations. We use formal model selection analysis based on the Watanabe Akaike information criterion to select the model that best predicts the pulmonary haemodynamics. Results show that the nonlinear pressure-area relationship with stiffness dependent on the unstressed radius predicts best the data measured in a control mouse.

Image-based scaling laws for somatic growth and pulmonary artery morphometry from infancy to adulthood

Pulmonary artery (PA) morphometry has been extensively explored in adults, with particular focus on intra-acinar arteries. However, scaling law relationships for length and diameter of extensive preacinar PAs by age have not been previously reported for in vivo human data. To understand preacinar PA growth spanning children to adults, we performed morphometric analyses of all PAs visible in the computed tomography (CT) and magnetic resonance (MR) images from a healthy subject cohort [n = 16; age: 1-51 yr; body surface area (BSA): 0.49-2.01 m²]. Subject-specific anatomic PA models were constructed from CT and MR images, and morphometric information-diameter, length, tortuosity, bifurcation angle, and connectivity-was extracted and sorted into diameter-defined Strahler orders. Validation of Murray's law, describing optimal scaling exponents of radii for branching vessels, was performed to determine how closely PAs conform to this classical relationship. Using regression analyses of vessel diameters and lengths against orders and patient metrics (BSA, age, height), we found that diameters increased exponentially with order and allometrically with patient metrics. Length increased allometrically with patient metrics, albeit weakly. The average tortuosity index of all vessels was 0.026 ± 0.024, average bifurcation angle was 28.2 ± 15.1°, and average Murray's law exponent was 2.92 ± 1.07. We report a set of scaling laws for vessel diameter and length, along with other morphometric information. These provide an initial understanding of healthy structural preacinar PA development with age, which can be used for computational modeling studies and comparison with diseased PA anatomy.NEW & NOTEWORTHY Pulmonary artery (PA) morphometry studies to date have focused primarily on large arteries and intra-acinar arteries in either adults or children, neglecting preacinar arteries in both populations. Our study is the first to quantify in vivo preacinar PA morphometry changes spanning infants to adults. For preacinar arteries > 1 mm in diameter, we identify scaling laws for vessel diameters and lengths with patient metrics of growth and establish a healthy PA morphometry baseline for most preacinar PAs.

Hemodynamics and Wall Mechanics after Surgical Repair of Aortic Arch: Implication for Better Clinical Decisions

Graft repair of aortic coarctation is commonly used to mimic the physiological aortic arch shape and function. Various graft materials and shapes have been adopted for the surgery. The goal of this work is to quantitatively assess the impact of graft materials and shapes in the hemodynamics and wall mechanics of the restored aortic arch and its correlation with clinical outcomes. A three-dimensional aortic arch model was reconstructed from magnetic resonance images. The fluid–structure interaction (FSI) analysis was performed to characterize the hemodynamics and solid wall mechanics of the repaired aortic arch. Two graft shapes (i.e., a half-moon shape and a crescent one) were considered. Material choices of the aortic arch repair included three commonly used graft materials (i.e., polytetrafluoroethylene (PTFE) synthetic graft, CorMatrix extracellular matrix, and pulmonary homograft) as well as one native tissue serving as a control. The pathological hemodynamic parameters, in terms of the percentage area of low wall shear stress (WSS), high oscillatory shear index (OSI), and high relative residence time (RRT), were quantified to be associated with potential clinical outcomes. Results have shown that the peak von Mises stress for the aortic arch repaired by the crescent graft was 76% less than that of the half-moon graft. Flow disturbance and recirculation were also minimized with the crescent graft. Moreover, pathological hemodynamic parameters were significantly reduced with the crescent graft. The graft material mismatch with the surrounding tissue aggregated the stress concentration on the aortic wall, but had minimal impact on flow dynamics. The present work demonstrated the role and importance of the graft geometry and materials on hemodynamics and wall mechanics, which could guide optimal graft decisions towards better clinical outcomes.

Uncertainty Quantification in a Patient-Specific One-Dimensional Arterial Network Model: EnKF-Based Inflow Estimator

Successful clinical use of patient-specific models for cardiovascular dynamics depends on the reliability of the model output in the presence of input uncertainties. For 1D fluid dynamics models of arterial networks, input uncertainties associated with the model output are related to the specification of vessel and network geometry, parameters within the fluid and wall equations, and parameters used to specify inlet and outlet boundary conditions. This study investigates how uncertainty in the flow profile applied at the inlet boundary of a 1D model affects area and pressure predictions at the center of a single vessel. More specifically, this study develops an iterative scheme based on the ensemble Kalman filter (EnKF) to estimate the temporal inflow profile from a prior distribution of curves. The EnKF-based inflow estimator provides a measure of uncertainty in the size and shape of the estimated inflow, which is propagated through the model to determine the corresponding uncertainty in model predictions of area and pressure. Model predictions are compared to ex vivo area and blood pressure measurements in the ascending aorta, the carotid artery, and the femoral artery of a healthy male Merino sheep. Results discuss dynamics obtained using a linear and a nonlinear viscoelastic wall model.