Vascular compromise and hemodynamics in pulmonary arterial hypertension: Model predictions

computer model of the based on (dog) The model

Where A a(0,35) is the cross sectional area of a given arterial vessel segment defined at a transpulmonary pressure (P tp ) of 0 and P tm of 35 cmH 2 O (a reflection of vessel dimensions under standardized conditions). The characteristic behaviour of A v as a function of P tm was also derived from data obtained by Smith and Mitzner (S3) and described in detail by Bshouty and Younes (S1). It is expressed by the following curvilinear relationship, where A v(0,35) is the cross sectional area of a given venous vessel segment defined at a P tp of 0 and P tm of 35 cmH 2 O (again, a reflection of vessel dimensions under standardized conditions). Perivascular pressures of the first three pulmonary arterial and venous generations (extra-alveolar extraparenchymal vessels) were assumed equal to pleural pressure (P pl ). Perivascular pressures of extra-alveolar intraparenchymal vessels are affected by both lung volume (and hence, P tp ) and P v . The relationship between P x of extra-alveolar intraparenchymal vessels and P v was also derived from data obtained by Smith and Mitzner (S4). The following relationship was used for both arteries and veins, The parameters a, b, and c were dependent on whether the vessel segment was arterial or venous and all three parameters varied linearly with P tp (for details see reference S1).
Vessel segmental lengths (l) of the first three arterial and venous generations were assumed constant and independent of lung volume (extraparenchymal vessels). The segmental length of extra-alveolar intraparenchymal vessels was exclusively a function of lung volume and hence, of P tp and lung compliance (C L ). This relationship was also derived from data obtained by Smith and Mitzner (S3). The following relationship was used for both arteries and veins, where l 0,35 is the vessel segment length at a P tp of 0 (again a reflection of vessel characteristics under standardized conditions), V M is minimal lung volume (at P tp =0) as a fraction of total lung capacity (for details see reference S1).
In the current simulations, P pl at the LA level was assumed equal to −5 cm H 2 O reflecting normal lungs at end-expiration. Pleural pressure for vessels above the LA was decreased by 0.55 cmH 2 O per cm height and increased by the same amount for vessels below the LA (S5). From dimensional changes, resistance of each segment was calculated assuming laminar flow. Total resistance of each vessel was calculated using numerical integration over the length of the vessel (1000 segments per vessel in the model).
The characteristic behaviour of the pulmonary capillary bed was modelled after the sheet model of Fung and Sobin (S6,S7). P x for the capillary bed was assumed equal to alveolar pressure (P al ). Details of capillary resistance calculations, under different zone conditions, and recruitment and derecruitment at the capillary level are beyond the scope of this paper. These were described in detail by Bshouty and Younes (S2).
The effect of gravity was accounted for by assigning the model a certain height (in the current simulations 20 cm) representing an average adult lung height in the supine position. All pressures were referenced to the LA. Flow in the model was assumed laminar and nonpulsatile.
The model was developed to represent 15 orders of pulmonary arteries between the main pulmonary artery and the capillaries and 15 orders of pulmonary veins between the capillaries and the left atrium according to Huang et al (S8). Each resistance in the model may be viewed not as a single vessel but representative of many vessels that reside at the same hydrostatic level. The first-generation vessel (arterial and venous in the current model) represents the cumulative resistances of extra-alveolar extraparenchymal vessels generations one to three. The second generation vessels in the model represent the cumulative resistances of extra-alveolar intraparenchymal vessels generations four to six. Each generation in the model represents three generations in the human lung so that the fifth generation in the model represents lung vessel generations thirteen to fifteen.
Although the wall structure of the pulmonary arteries varies considerably between the main pulmonary artery all the way down to the precapillary terminal arteries (S9), the distensibility of the vessels is independent of vessel diameter (S9-S12). Hence, the distensibilities of the arterial vessels were assumed equal across the five arterial generations. Since a similar observation was seen in pulmonary veins (S13), the distensibilities of the venous vessels across the five venous generations were also assumed equal, under normal conditions.

Model adaptation to humans
The model was originally developed based on data obtained in open chested dogs. Lung compliance was changed to equal what is normally observed in humans (S14) after correcting static compliance to liter (L) per cmH 2 O according to the following relationship, where Ht is patient height in cm. Because the vascular structure of the pulmonary circulation as seen on pathology is similar in both humans and dogs, the characteristic behaviour (cross-sectional area as a function of P tm and vessel length as a function of lung volume) of extra-alveolar arteries and veins were kept the same. A similar assumption was used for the characteristic behaviour of the capillary bed.
The longitudinal distribution of pulmonary vascular resistance along the vascular bed was based on data obtained by Huang et al (S8). An average diameter ratio of parent to daughter vessel of 1.56 was assumed for all pulmonary arteries and 1.58 for all pulmonary veins under normal conditions. An average length ratio of parent to daughter vessel of 1.49 was assumed for all pulmonary arteries and 1.50 for all pulmonary veins. An average branching ratio (daughter vessels to parent vessel) of 3.36 was used for all pulmonary arteries and 3.33 for veins. These ratios were obtained at a perfusion pressure of 3 cmH 2 O while alveolar gas pressure was maintained at 10 cmH 2 O and pleural pressure was zero (S8). Based on the above, baseline vessel resistance of each vessel in each generation was calculated. Using the characteristic behaviour of arteries and veins (equations 1, 2, and 4) baseline resistances were adjusted to the current model standardized conditions of P tp of 0 cmH 2 O and P tm of 35 cmH 2 O. Transition from the above dimensional characteristics to provide actual resistances along the pulmonary vascular bed was achieved by obtaining a PAPm of 15 mmHg at a CO of 6.45 L/min (Cardiac index 3.38 L/min/m 2 ) and a LAP of 5 mmHg reflecting the normal pulmonary circulation in an adult subject 175 cm tall and weighing 75 kg (S15). The contribution of capillary resistance to total resistance under these conditions was approximately 16% with the arterial and venous systems contributing equally to the remaining resistance (S16). Although in the current simulations a height of 175 cm and weight of 75 kg were used, the model accepts a height range of 30 cm to 225 cm and a weight range of 2 kg to 300 kg. Baseline resistances in the model are adjusted automatically to accommodate different body sizes. This is achieved by assuming that at a normal cardiac index of 3.38, all subjects (irrespective of height and weight) with a normal pulmonary circulation have a PAPm of 15 mmHg at a LAP of 5 mmHg (S15).

Model adaptation to PAH
The predominant cause of the increase in pulmonary vascular resistance in PAH is loss of vascular luminal cross-sectional area due to vascular remodeling produced by excessive hypertrophy and proliferation of the intima, media, and adventitia of small pulmonary arteries (usually <2000 µm in diameter) (S17). Vasoconstriction plays a role in only a minority of cases as less than 13% of the patients are vasoreactive at the time of investigation (S18). Hence, when simulating PAH, vascular compromise in the model was confined to the fifth arterial generation (representing generations thirteen to fifteen in the human lung). Because the same vessels are involved in regulating the regional distribution of blood flow in the lung, they were also assumed the site of vasoconstriction (S19). Although, smooth muscle tone, representing vasoconstriction, could be added to any vessel in the model (For a complete description of modeling smooth muscle tone see reference [S1]).
Initially, simulations were performed without the addition of vasoconstriction. Simulation of PAH was achieved by assigning the precapillary (fifth arterial generation) vessels increasing levels of vascular compromise. Vascular compromise was imposed by assuming a smaller cross-sectional area for the fifth arterial generation at P tp =0 and P tm =35 cmH 2 O. In the current simulation, vascular compromise was applied equally to all sixteen precapillary vessels. The characteristics of all other vessels were assumed normal.

2•A a(0,35) . Notice that as vascular compromise increases, vessel compliance decreases
There are no human data comparing pulmonary arterial distensibility of healthy subjects to that in patients with pulmonary arterial hypertension. In a study by Molthen RC et al (S12), arterial distensibility (over a pressure range of 5 mmHg to 30 mmHg) of Sprague-Dawley rats exposed to 21 days of hypoxia was compared with controls. The relationship between vessel diameter and pressure remained linear in arteries ranging in size from 70 µm to 2800 µm, and distensibility (the slope of the relationship) dropped from 0.028 to 0.015. This fixed drop in distensibility is likely the result of the group (Sprague-Dawley rats) and study conditions (degree of exposure) homogeneity. In patients suffering from pulmonary arterial hypertension, the decline in distensibility is likely gradual and correlates with disease progression. As shown in equation 1, the relationship between A a and P tm is linear. Because A a(0,35) is multiplied by both slope and intercept, reducing A a(0,35) does not only decrease the intercept (ie, parallel shift) of the relationship between A a and P tm , but also the slope of this relationship reflecting a change in vascular compliance. As shown in Figure S2, compared with baseline (solid line) a progressive decline in vascular cross-sectional area in steps of 0.2 (from 0.8•A a[0,35] to 0.2•A a [0,35] ) is associated with a progressive and proportional decline in both intercepts and slopes thus describing a change in the passive characteristics of these vessels (hence, the use of vascular compromise in this article rather than simple narrowing). This type of behaviour is anticipated even when merely based on simple proliferation of vascular wall layers (media and intima).
Simulation of vasoconstriction without vascular compromise is achieved by keeping A a(0,35) at baseline. In this case, compromise is active and occurs because of an effective reduction in P tm (equivalent to the amount of applied vasoconstriction) along the original curve (again, for details on the simulation of smooth muscle tone see reference S1).

Determination of PAPm and total PVR
As mentioned above, baseline resistances of each vessel in each generation are assigned under standardized conditions (P tp =0 and P tm =35 cmH 2 O) and adjusted automatically to accommodate different patient heights and weights. The operator enters three groups (general, patient specific, and hemodynamic) of data into the model. General data includes the reference point for pressure measurement (top or bottom of the lung, or LA) and the tolerance level for a given run (see below). Patient specific data includes patient height, weight, lung height, lung compliance, and lung volume at a P tp =0. Hemodynamic data includes the point of measurement during the respiratory cycle (as determined by P al and P pl ) and LAP and CO. The operator is given the opportunity to modify any baseline characteristics of any given vessel (including the addition of vasoconstriction to any vessel) in the model before proceeding with the calculations. PAPm and total PVR are determined by iteration. Initially, flow is assumed to be equally distributed across all vascular channels. Given the initial distribution of resistances, vascular pressures across the pulmonary circulation are calculated. After accounting for the effects of gravity and P x changes, transmural pressures at each segment of each vessel are calculated and using numerical integration total resistance of each vessel is recalculated. Using the new distribution of resistances, a new distribution of flows is recalculated with a new distribution of pressures. The computer continues cycling (iterating) through the above process until a stable solution is reached when flows, pressures and resistances stabilize (stop changing) across the whole vascular bed (in the current simulations less than 0.1%). The output includes the distribution of individual pulmonary vascular flows, arterial, capillary, and venous resistances across the whole pulmonary circulation including summary data of PAPm, PVR, upstream (arterial), middle (capillary), and downstream (venous) resistances.

Effect of branching on model predictions
To assess whether model predictions are affected by the levels of branching, results obtained with the four (original) and five generations models, over the whole range of cardiac outputs and left atrial pressures reported here, were compared. In the four generations model, adaptations to humans and PAH (see above) were similar to the five generations model except for the following differences. In the four generations model, the first generation represented the first three generations in human lung and the rest of the generations (two, three and four), each represented four generations in the human lung so that the total number of generations represented was the same between the two models (15). Simulation of PAH was achieved by compromising the fourth generation vessels in the four generations model, and similarly, vasoconstriction was applied to the fourth-generation vessels. Simulation results obtained with the four and five generations models were identical to the second decimal point and hence, were not clinically significant. All simulations and results presented here were generated using the five-generation model.

FOOTnOTE:
The latest version of the simulation program includes both a five-and a 15-generation (16,384 parallel vessels) model, with each vessel divided into 10,000 segments. The program is available, free of charge and upon request, to all PH centres.