Abstract
Right-ventricular function is a good indicator of pulmonary arterial hypertension (PAH) prognosis; however, how the right ventricle (RV) adapts to the pressure overload is not well understood. Here, we aimed at characterizing the time course of RV early remodeling and discriminate the contribution of ventricular geometric remodeling and intrinsic changes in myocardial mechanical properties in a monocrotaline (MCT) animal model. In a longitudinal study of PAH, ventricular morphology and function were assessed weekly during the first four weeks after MCT exposure. Using invasive measurements of RV pressure and volume, heart performance was evaluated at end of systole and diastole to quantify contractility (end-systolic elastance) and chamber stiffness (end-diastolic elastance). To distinguish between morphological and intrinsic mechanisms, a computational model of the RV was developed and used to determine the level of prediction when accounting for wall masses and unloaded volume measurements changes. By four weeks, mean pulmonary arterial pressure and elastance rose significantly. RV pressures rose significantly after the second week accompanied by significant RV hypertrophy, but RV stroke volume and cardiac output were maintained. The model analysis suggested that, after two weeks, this compensation was only possible due to a significant increase in the intrinsic inotropy of RV myocardium. We conclude that this MCT-PAH rat is a model of RV compensation during the first month after treatment, where geometric remodeling on EDPVR and increased myocardial contractility on ESPVR are the major mechanisms by which stroke volume is preserved in the setting of elevated pulmonary arterial pressure. The mediators of this compensation might themselves promote longer-term adverse remodeling and decompensation in this animal model.
Introduction
Pulmonary arterial hypertension (PAH) is a severe disease of the pulmonary vasculature that is characterized by a resting mean pulmonary arterial pressure (mPAP) ≥ 25 mmHg with a pulmonary arterial wedge pressure ≤ 15 mmHg. 1 During the development of PAH, the right ventricle (RV) of the heart becomes subjected to pressure overload, which results in chamber remodeling that initially compensates for the increase in pulmonary vascular resistance (PVR). Eventually this response becomes maladaptive and, if left untreated, PAH results in cardiac decompensation and heart failure. 2 Prior investigations in PAH patients have indicated that increased RV diastolic stiffness is related to adverse clinical outcomes 3 and reduced survival. 4 However, the process of remodeling and adaptation that the RV undergoes in response to pulmonary pressure overload is not clear, making it harder to determine when the RV goes from adapting to failure.
To relate pulmonary hemodynamics to right ventricular mechanics and thereby better characterize the time course of RV early remodeling during PAH, we carried out a longitudinal study in an animal model of PAH. Using invasive measurements of RV pressure and volume in a monocrotaline (MCT) animal model of PAH, ventricular morphology and function were assessed weekly during the first four weeks after MCT exposure. To test whether the observed changes in right ventricular diastolic and systolic function could be explained by the observed changes in right ventricular preload, afterload, and wall geometry alone, we compared measured RV pressure-volume (PV) end-diastolic and end-systolic relations with those predicted by a computational model that included measured geometrical chamber remodeling.
RV systolic and diastolic pressures rose significantly after the second week of PAH induction accompanied by significant RV hypertrophy, but RV stroke volume (SV) and cardiac output (CO) were maintained. The model analysis suggested that, after two weeks, this compensation was only possible due to a significant increase in the intrinsic inotropy of RV myocardium.
Methods
Animal preparation
All experiments were performed with the approval of the Animal Care and Use Committee at the University of Illinois, Chicago. The study included 45 male Sprague Dawley rats weighing 275 ± 12 g (Charles River, Chicago, IL, USA). PAH was induced in the treatment group by a subcutaneous injection of 60 mg/kg of monocrotaline (crotaline, Sigma-Aldrich). The longitudinal study of PAH was carried out by weekly evaluations of the RV-PA hemodynamics for up to four weeks post injections. The treated rats studied after one week of the injection (PAH1) represent the early stage whereas the group studied after four weeks of receiving the injection (PAH4) corresponds to the late stage of PAH. A control group was injected with saline matching the volumetric dose of the treated group. To account for the aging effects, the control group (NCTL = 7) consisted of animals injected after one, two, three, or four weeks (NCTLweek1 = 1, NCTLweek2 = 2, NCTLweek3 = 3, NCTLweek4 = 1). When possible, a treated PAH animal was paired with a control (CTL) animal and housed together for up to four weeks. The group sizes for each disease stage were NPAH1 = 8, NPAH2 = 8, NPAH3 = 6, and NPAH4 = 16.
In-vivo hemodynamic studies
The paired rats underwent surgery the same day, one after another. The animals were anesthetized using isoflurane (Henry Schein, NY, USA). Anesthesia was initially administered by placing the animal in an induction chamber supplied with oxygen mixed with 4% isoflurane. When the animal was confirmed to be unconscious, the anesthesia was decreased to 2.5% and supplied via a nasal cone. Following a tracheotomy, the rat was intubated (E-Z Anesthesia Ventilation System, Palmer, PA, USA) and supplied with a respiratory rate of 50 breaths per minute and a tidal volume based on the animal's weight. 5 The heart and great vessels were exposed via a thoracotomy. Once the chest cavity was opened and the RV was identified, a 21-gauge needle was used to puncture the ventricular wall apically. The needle was immediately replaced by a 1.9F PV admittance catheter (Transonic Scisense, ON, Canada) introduced along the long axis of the chamber. Once PV waveforms recorded were reproducible in steady state (steady-state conditions), the inferior vena cava (IVC) was slowly occluded with a suture line (Fine Science Tools Braided Silk Suture 5-0, Foster City, CA, USA) for a few seconds to change the venous return (occlusion conditions). On average, five occlusions were performed on each animal. Following PV measurements, a second puncture was made on the RV base next to the outflow track using a 23-gauge needle. A 1.6-F dual-pressure sensor catheter (Transonic Scisense, ON, Canada) was inserted into the RV and advanced into the main pulmonary artery (MPA) via the pulmonic valve for simultaneous measurements of the RV and MPA pressures. The MPA pressure time series was used to calculate mPAP and thus confirm PAH. At the conclusion of the hemodynamic measurements, the animals were exsanguinated and the heart and lung vessels were flushed with phosphate buffered saline (PBS; 0.14M NaCl; 0.0027M KCl; 0.0017M KH2PO4; 0.01M Na2HPO4) and heparin (Fresenius Kabi, USA). Data were collected in LabChart Software (AD Instruments Inc., Version 8.0 Pro, Colorado Springs, CO, USA).
Hemodynamic metrics such as SV, ejection fraction (EF), heart rate (HR), CO, total PVR, and ventricular relaxation time constant τ were computed through a custom-written MATLAB (Mathworks Inc, Version 8.6.0, Natick, MA, USA) code. The end-systole (ES) point was determined by finding the maximum value of the pressure to volume ratio in a given pressure-volume loop. The end-diastole (ED) point was the first point in the pressure-volume loop with highest volume and lowest volume rate of change (i.e. dV/dt ∼0). Based on the difference in volume at ED and ES, SV was calculated. EF was determined by the ratio of SV and ED volume. Using the calculated SV and HR, CO was computed by taking the product of the two. Total PVR was computed as the ratio between RV ES pressure and CO.
6
The RV relaxation time τ was calculated by fitting an exponential decay function to the pressure waveform
Right-ventricular and pulmonary vascular elastances
Ventricular contractility or end-systolic elastance Ees corresponds to the slope of the linear regression to the end-systolic points during occlusion of the IVC (approximately 40% drop in blood pressure). Based on the steady-state PV measurements, the effective pulmonary arterial elastance Ea or the slope of the line between the ES point and ED volume at zero pressure was obtained. Ea is a parameter that relates to the ability of the pulmonary vasculature to buffer the incoming blood pressure and volume influx. By taking the ratio of Ees and Ea, the ventriculo-vascular coupling η can be assessed.
10
Finally, passive ventricular chamber stiffness was calculated by fitting the ED PV data points during occlusion of the IVC using the following exponential function
Relating ventricular pressure-volume to ventricular wall mechanics at ES and ED
To compute RV end-diastolic and end-systolic pressure-volume relations from RV geometry and myocardial resting and active material properties, a spherical model of the RV was developed. The RV model comprehends the RV cavity, RV free wall, and septum. Assuming that the RV is a fraction of a sphere with an outer radius R, the volume of the sphere is the sum of the RV free wall volume
Summary statistics of hemodynamics, elastances, and RV structural changes observed during PAH progression.
P < 0.05 vs. CTL.
P < 0.05 vs. PAH1.
P < 0.05 vs. PAH2.
At the unloaded volume, V
RV
= V0 and R = R0, so that
The radius R and wall thickness h of the spherical model in terms of V0 and R0 are computed as:
To avoid over-restricting the model geometry, the experimentally measured wall thickness was not directly incorporated in the computational modeling. Instead, Eq. 3 for wall thickness was used. A schematic of the model is illustrated in Fig. 1.
Schematic of the RV computational model modeled as a fraction of a sphere of outer radius R and mid-wall radius r. In a cross-section of the sphere, the RV cavity VRV, RV free wall Vw, and septum Vsw volumes illustrate the different compartments of the RV model.
For sarcomere mechanics,
12
the myocardium passive stress was defined as
Given that the RV wall is thin compared to its cavity radius,
15
the thin-walled LaPlace's law was used to estimate mean mid-wall stresses. Here, the RV cavity blood pressure P is computed as
To predict ESPVR and EDPVR, PV data from each individual animal obtained during IVC occlusions were fitted using the linear and exponential functions described in the previous section. Using the fitted functions corresponding to each individual animal over a common pressure and volume range, a mean ± standard error curve was generated to represent the measured data per group. While keeping the parameters L0, L
R
, Ca0, and Ca
max
constant, the material property parameters
Statistical analysis
Descriptive statistics were computed for all the hemodynamic measurements, RV tissue measurements, and elastances. To detect significant differences among groups, non-parametric Wilcoxon/Kruskal–Wallis (Chi-square approximation) tests were conducted in cases where observations displayed departures from a normal distribution; otherwise, parametric analysis of variance was carried out. When a difference in treatment was found to be significant, multiple group comparisons were then conducted to identify where the differences occurred between the groups (CTL, PAH1, PAH2, PAH3, PAH4). Multiple comparisons Tukey HSD were used for the normally distributed data; otherwise the Dunn Method was implemented. All the statistical analyses were carried out in JMP Pro Statistical Software (version 13.0.0, © SAS Institute Inc., Cary, NC, USA). For all statistical tests, significance level was set at 0.05. All values are presented as mean ± standard deviation, unless specified.
Results
Animal preparation
Out of 38 animals administered with MCT, 26 developed PAH. Two animals were normotensive, four unexpectedly died at the beginning of the surgical procedure unable to confirm PAH, and six died between weeks 3 and 4 post injection. RV samples and weights were collected from the 19 out of 26 treated and six out of seven control animals.
In-vivo hemodynamics
The PV loops from each group exhibited unique hemodynamic characteristics. The CTL group exhibited a square-shape loop through the four weeks with tight systolic pressure (31.7 ± 6.1 mmHg) and volume ranges (90.8 ± 19.3 to 239.7 ± 19.6 µL). In the case of the treated animals, RV ES pressures rose each week from 27.3 ± 2.4 up to 70.5 ± 6.8 mmHg (P value = 0.0004). ED pressures also increased, significantly in the late PAH group (from 3.5 ± 1.4 to 7.7 ± 1.7 mmHg, P value = 0.028, PAH4 vs. PAH1). While the width of the PV loops or SV remained statistically the same (from 141.1 ± 22.9 to 160.3 ± 6 µL), there was a shift to the right of ES volume (from 79.2 ±19.3 to 140.9 ± 35.6 µL, P value = 0.0190) and ED volume (from 220.3 ± 31.6 to 301.3 ± 37.8 µL, P value = 0.0189) as weeks advanced, indicating the enlargement of the ventricle and an increase in resting volume. CO was also maintained in disease and remained statistically the same for all groups. The increase in mPAP (from 25.2 ± 2 to 42.8 ± 4.9 mmHg, P value = 0.0001) and total PVR (from 0.7 ± 0.2 to 1.6 ± 0.2 mmHg/mL min, P value = 0.0149) was reflected in the PV loop shape by producing the formation of a shoulder during the ejection phase (slanted/triangular shape). The increase in ED volume and preserved SV caused a slight decrease in EF (64.2 ± 6.4% to 53.9 ± 5.5%). Ventricular diastolic function parameters dP/dtmin and τ decreased, with a statistically significant difference only found for dP/dtmin at the late disease stages (vs. normotensive and early PAH, P < 0.0001). Summary statistics of the hemodynamic measurements are shown in Figs. 2 and 3 and Table 1.
Summary statistics of cardiac performance during systole (a) and diastole (b). *P < 0.05. Summary statistics of RV hypertrophy evidenced with Fulton Index, RV thickness, and RV weight. (a) No changes in animal weights are shown. (b) mPAP and total PVR. *P < 0.05.

Right-ventricular morphology
Right ventricular morphology changed with the progression of the disease (Fig. 2). The ventricle progressively became hypertrophic (Fulton index from 0.31 ± 0.05 to 0.52 ± 0.03, P value = 0.0011, RV weight from 0.23 ± 0.04 to 0.47 ± 0.05 g, P value = 0.0316) and thickened significantly (from 0.6 ± 0.1 mm in the normotensive group to 1.1 ± 0.1 mm in advanced PAH, P value < 0.0001) by week 2 post treatment. The weight of the animals, however, remained the same.
Right-ventricular and pulmonary vascular elastances
Effective arterial elastance Ea significantly increased at the late PAH stages (PAH3 and PAH4 groups). Ees initially did not change (CTL, PAH1, and PAH2 groups) but then increased significantly in the PAH3 group (0.25 ± 0.08 to 0.34 ± 0.1 mmHg/µL, P value = 0.0435) and further increased in PAH4 (0.45 ± 0.13 mmHg/µL, P value = 0.0460 vs. CTL). End-diastolic elastance Eed dropped initially (CTL: 0.030 ± 0.01 to PAH1: 0.015 ± 0.006 mmHg/µL), but then progressively increased reaching baseline value at week 4 (0.033 ± 0.005 mmHg/µL). γ and β were not statistically different across the groups. γ was unaffected at the early PAH stages and displayed an increase after three weeks of MCT treatment. Diastolic chamber stiffness β, on the other hand, decayed in the first week of MCT treatment and remained overall constant until late PAH. The RV and pulmonary system remained coupled during the progression of PAH. This was reflected in the ventricular-vascular coupling efficiency η, which was not statistically different among the treatment groups. Ees, Ea, Eed, and η are summarized in Fig. 4.
(a) Representative occlusion PV loops of control and after 4 weeks of MCT administration (chronic PAH). Ees: slope of dashed black line (ESPVR); Ea: slope with dashed dark gray line; Eed: slope with light gray line. (b) Summary statistics of Ees, Ea, Eed, and η are also depicted. *P < 0.05.
Ventricular hemodynamics related to ventricular wall mechanics at ES and ED
The RV model was able to recapitulate the experimentally measured mean RV pressure-volume relations at ES and ED (Fig. 5a, left). Using geometric measurements from PAH groups (solid lines), the RV model with the control-group material properties overestimated EDPVR at week 1, closely matched the mean EDPVRs measured in the PAH groups at weeks 2 and 3, and underestimated EDPVR at week 4. The ESPVR model predictions were highly underestimated for the week 2, 3, and 4 PAH groups. However, when accounting for the geometrical changes and adjusting the material properties k1 and T
max
to increase intrinsic myocardial passive and active tension development, predictions from the model (dashed line) all showed good agreement with experimentally measured mean ESPVRs across the ranges of volumes measured in the PAH groups. Similar findings and the same conclusions were obtained when the entire analysis was repeated with altered assumed fraction of septal wall in the LV wall mass. Model parameters are shown in Table 2, with the adjusted k1 and T
max
values to account for a given PAH group material properties.
(a) Comparison of measured raw (open circles) and averaged (blue dots and bars) right-ventricular end-diastolic and end-systolic pressure-volume relations with computational model predictions. Control: Model-predicted end-diastolic and end-systolic PV relations (squares and gray lines) using control-group geometry and control-matched material properties. PAH1-4: Model-predicted end-diastolic and end-systolic PV relations (squares and gray lines) using control-group geometry and control material properties. Model-predicted end-diastolic and end-systolic PV relations (solid lines) using PAH geometry and control material properties. Model-predicted end-diastolic and end-systolic PV relations (dashed lines) using PAH geometry and PAH-matched material properties. (b) Model-predicted isometric end-systolic and end-diastolic sarcomere length-tension relations for all 5 animal groups determined by optimizing material parameters to match measured RV pressure-volume relations. Note the increase in predicted systolic sarcomere elastance 2–4 weeks after MCT. End-diastolic sarcomere stiffness tension decreases at week 1 of PAH and increases by. Parameters used in the computational model of the RV to predict PV relations at end-systole and end-diastole. The morphology parameters V0, Vw, Vsw are the unloaded volumes, RV free wall volume, and volume of septum wall, respectively. The material property parameters θC = {k1C, k2C, TCmax, BC} are based on the model predictions of ES PV and ED PV of the control animals. When changes in morphology did not match ESPVR and EDPVR, this gave rise to a new set of adjusted material properties – in this case, k1 and Tmax.
Based on the RV model predictions of ESPVR and EDPVR with the adjusted parameters (Table 2), active and passive sarcomere length-stress curves were generated for each treatment group. As shown in Fig. 5b, the control, PAH2 and PAH3 passive length-stress curves closely overlap. However, there was a decrease in muscle stiffness in PAH1 (from k1 = 2.3 to 1.4 mmHg) and a large increase in muscle stiffness at PAH4 (from k1 = 2.3 to 3.5 mmHg). The active sarcomere length-stress relation in the control and PAH1 groups were the same. However, this relation changed for the PAH2 group, where the contractility of the muscle increased—from T max = 640 to 1300 mmHg—and remained high (T max = 1100 and 1200 mmHg) in the following weeks.
To discriminate the contributions of RV geometric remodeling from changes in diastolic and systolic myocardial material properties on the preservation of RV SV, the model was used to compute RV end-diastolic and end-systolic pressure-volume relations using three combinations of model geometry and material properties (Fig. 6): control-group geometry and material properties (dotted lines in Fig. 6b); PAH geometries and control material properties (solid lines in Fig. 6b); PAH geometries and PAH material properties (dashed lines in Fig. 6b). From these computed ED and ES PV relations, SVs were computed to quantify how the changes in hemodynamics, wall geometry and passive and active myocardial material properties combined to maintain CO in the presence of altered RV hemodynamics.
(a) Model-predicted changes in RV SVs from weeks 1–4 due to: (SV1) changes in preload and afterload only; (SV2–SV1) the additional contribution of geometric remodeling on the EDPVR; (SV3–SV2) the further contribution of altered diastolic compliance on the EDPVR; (SV4–SV3) the further contribution of geometric remodeling on the ESPVR; (SV5–SV4) the further contribution of altered contractility on the ESPVR. (b) Graphical definitions of model-predicted SVs SV1–SV5 obtained from model-predicted EDPVR and ESPVR computed using: control-group geometry and control-group material properties (dotted lines); PAH geometries and control material properties (solid lines); PAH geometries and PAH material properties (dashed lines).
The SV labeled SV1 in Fig. 6b represents the SV predicted under PAH loading conditions if there had been no geometric remodeling or changes in myocardial diastolic or systolic material properties. Without compensatory changes in ventricular geometry, elastic and contractile properties, the steep decrease in SV1 indicates that right ventricular function would have been reduced by 75% after four weeks, and that elevated RV EDP had not been sufficient to compensate after week 1. SV2 is the SV computed when the contribution of geometric remodeling on the EDPVR was accounted for in the models. It shows that RV dilation after week 1 was a major contributor to the preservation of RV SV. The contribution of altered geometry and diastolic compliance on the EDPVR is represented by SV3. The transient decrease in RV myocardial stiffness was the major compensatory mechanism at week 1, but after week 1, the contribution of altered RV passive stiffness was minor and actually negative as the myocardial stiffness increased by week 2. SV4 represents the contribution to maintenance of SV of diastolic material and chamber remodeling and the effects of geometric remodeling on the ESPVR. Subtracting SV3 shows that the effects of changes in ventricular geometry on the ESPVR actually worked slightly against compensation. The increase in ventricular radius that increased ED volume was not fully overcome by the increase in RV wall thickness. Finally, when the effects on SV of chamber remodeling and diastolic and systolic material properties on both the EDPV and ESPV relations were all included in the model (SV5), the additional effect of increased contractility was also shown to be an important compensatory mechanism. In summary, the effects of geometric remodeling on the EDPVR and of increased myocardial contractility on the ESPVR were the major mechanisms by which SV was preserved in the setting of elevated PAP, and both mechanisms contributed equally during weeks 2–4. In week 1, the minor loss of SV that the modest hemodynamic changes would have caused was fully compensated by a transient decrease in diastolic myocardial stiffness before any geometric remodeling or contractile compensation had begun.
Discussion
In the present study, RV hemodynamic and morphological measurements were taken to investigate the compensated phase of RV hypertrophy in the MCT-animal model of PAH. Our goal was to discriminate the contributions of ventricular geometric remodeling and intrinsic changes in myocardial mechanical properties at ES and ED to the maintenance of CO using a computational model. In this longitudinal study, the myocardium was found to initially decrease its stiffness—an acute adaptation that preserved SV at week 1 of PAH induction and reversing by week 2. This transient early decrease in diastolic stiffness allowed the RV to soften initially in order to preserve diastolic filling and therefore maintain SV. By week 2, the PAP, PVR, and end-systolic volume all increased while maintaining SV, CO, and EF. This compensated response could be explained in the model almost entirely by the changes in chamber volume and mass. However, after week 2, the changes in RV volume and wall thickness and increased filling were not sufficient to maintain SV and CO as PAP continued to increase. This led us to conclude that there was a significant increase in RV contractility during this period as shown by the increased T max values required to match the experiment-derived ESPVR. Furthermore, this MCT-PAH rat is a model of RV compensation during the first month after treatment, where geometric remodeling on EDPVR and increased myocardial contractility on ESPVR are the major mechanisms by which SV is preserved in the setting of elevated PAP. The mediators of this compensation might themselves promote longer-term adverse remodeling and decompensation in this animal model.
MCT-model of PAH shows features of preserved cardiac function
In our study, RV pressures rose and volumes increased while maintaining CO, SV, and EF. Our results are in agreement with Handoko et al., 17 where CO was preserved in rats with no significant difference in total PVR by day 14 after MCT injections. In their study, like ours, by day 28 after MCT injections, there was a significant increase in total PVR while CO continued to be preserved. Yin et al. also reported maintenance of EF and CO in their four-week of MCT-rat model. 18 In addition to CO being preserved, Werchan et al. also measured significant decrease (more negative values) in dP/dtmin in four-week MCT-treated rats. 19 Like Werchan et al., magnitude of the rate of pressure change during relaxation (dP/dtmin) became significantly more negative, but relaxation time τ only slightly decreased in advanced PAH. This discrepancy has been previously reported by Umar et al., where no change in the relaxation time τ was found in the presence of an increasing trend in dP/dtmin magnitude 20 or vice versa by Correia-Pinto et al. in MCT-induced PAH rats. 21 These findings suggest that in the four-week MCT model, the RV maintains diastolic functional elasticity and is still in a compensatory phase. The work by Handoko et al. showed that after these four weeks, there is a statistically significant decline in CO, indicating the end of the RV compensation phase.
Sudden deaths between week 3 and week 4 after MCT injections occurred in our study. MCT is a toxic alkaloid that when ingested produces a syndrome characterized by pulmonary hypertension, pulmonary vasculitis, and RV hypertrophy. 22 To date, the toxicological mechanisms of the MCT in the lung are unknown. Abe et al. have suggested that MCT-injected rats probably die of cardiac and renal dysfunction. 23 In a study by Umar et al. they reported sudden death in their 30% death rate 23 days after MCT injections and 100% by day 32 post injections, attributing this to spontaneous ventricular fibrillation. 24 Several groups have reported varying rates after three weeks of treatment.25–28 Some prominent studies indicate the presence of a proarrhythmic substrate in hearts of MCT-treated rats, and among the causes of sudden death are the changes in Ca2+ and K+ channel activity and electrical restitution.29–32
Pulmonary arterial stiffness increases in late PAH
Pulmonary arterial stiffening markedly increased and became significantly different from the control group to the MCT-injected rats after three weeks. Pulmonary arterial elastance Ea values increased as PAH advanced, which has also been reported by other groups when using the rat MCT model: Yin et al. from 0.5 ± 0.14 to 0.79 ± 0.3 mmHg/µL; 18 and de Man et al. from 0.15 to 1.0 mmHg/µL. 33 Total PVR increased by twofold in our studies, which suggests the pruning of pulmonary vessels. While we did not carry out a histological analysis, Yin et al. did in MCT-treated animals after four weeks of injections. In their histological analysis of the lungs, the pulmonary arterioles showed marked lesions (i.e. thickening, narrowed, and inflamed arterioles 18 ), all consistent with characteristics found in PAH patients.
Ventricular contractility and stiffness in PAH cannot be explained by morphological changes alone
SV could only be maintained after two weeks of MCT treatment by the significant increase in end-systolic elastance. Increase in Ees has been previously reported to increase by Yin et al. (from 0.67 ± 0.12 to 0.95 ± 0.42 mmHg/µL), 18 de Man et al. (from 0.2 to 0.5 mmHg/µL), 33 and Blaudszun et al. (from 0.24 ± 0.12 to 0.56 ± 0.33mmHg/µL) 34 from normo- to hypertensive states in MCT-models of PAH in rats. The model analysis showed that wall thickening and geometric remodeling alone were unable to explain all of the increased chamber elastance. The RV model used to predict ESPVR underestimates the data measured at ES of the PAH2, PAH3, and PAH4 groups unless the material parameters were adjusted and the maximal isometric tension Tmax is increased by twofold. Predicted peak systolic stress-strain relations showed that intrinsic myocardial contractility increased significantly by week 2 and remained significantly higher than the control group. Rain et al. 1 found that active tension of RV myocytes from PAH patients was higher than those of normotensive subjects. Interestingly, CO in the PAH group was lower than the control group but without reaching statistical significance. Sabourin et al. studied RV cardiomyocytes from a three-week rat MCT model of PAH and in this group of hypertrophied RV myocytes found an increase in the expression of SERCA2 and the amplitude of calcium transients, which resulted in an increase in contractility. This would explain the results from our model analysis and Rain et al.'s measurements in samples from PAH patients. Sabourin et al. concluded that enhanced SR calcium load was mediated by STIM1L and TRPCs/Orail. The parallel in these studies indicates that MCT is a valuable model of PAH for mimicking the initial RV adaptation found in PAH patients. 2
In our study, we observed a transient decrease in Eed after one week of PAH, followed by a recovery to control values during PAH progression. The RV model analysis suggests that this initial decrease was due to a temporary increase in myocardial compliance. This myocardium intrinsic change, along with RV wall thickening as hypertrophy progressed, contributed to the restoration of end-diastolic chamber stiffness subsequently. The temporary increase in diastolic myocardial compliance followed by the subsequent increase in ED pressure contributed significantly to the maintenance of SV via the Frank-Starling mechanism by increasing end-diastolic volume. In humans with PAH, RV diastolic stiffness Eed has been shown to be sensitive to the stage of the disease. Trip et al. reported an increase in Eed from control (∼0.25 mmHg/mL) to PAH patients. 3 In their study, they determined that diastolic stiffening of the RV is independent of the systolic adaptation process. They also suggested that Eed increment in PAH patients with a survival > 5 years may be explained by the hypertrophy of the RV, while in patients with < 5 years such increase may be related to additional intrinsic alterations of the myocardium. One possible mechanism of the acute loss of myocardial stiffness may be an increase in the fraction of the more compliant N2BA isoform of titin as a fraction of total titin, as was reported by Rain et al. 35 in rats after seven weeks of mild PA banding. While Eed was not statistically different from the control and PAH4 groups, the model shows an increase in EDPVR stiffness and points towards stiffening of the passive sarcomere length-tension relation, which rose in PAH4.
Ventriculo-vascular coupling in PAH remains unaffected
Several studies have used η to evaluate the adaptation and performance of the RV during PAH. In a study carried out by de Man et al. in a MCT-induced PAH rat model, they reported a decrease in η from the control to severe disease state (from 1.0 to 0.5). 33 Even though Yin et al. used the same animal model, they reported no change in η when comparing normo- to hypertensive MCT rats (control: η = 1.44 ± 0.08, PAH: η = 1.42 ± 0.33). 18 Our outcomes showed that the RV and the pulmonary arteries are coupled through the disease progression, indicating that the heart is still keeping up with the high-pressure overload. Here, we found a hypertrophic heart, with an increase in ES volume and ED volume but a preserved SV, no change in CO, EF nor HR. CO is preserved even though the resting volume of the RV increased. This supports the claim that the RV is still in the adapting phase and working along the increase in PA pressures.
Limitations
Hemodynamic measurements were carried out using isoflurane as an anesthetic during open-chest surgery. Some studies have shown that isoflurane can affect hemodynamics measurements by reducing PAP, total PVR, and contractility.36,37 Ewalenko et al. demonstrated that while true that pulmonary vascular tone decreased with isoflurane compared to propofol, this decrease was negligible and did not affect impedance. Given that all animals underwent the same in-vivo measurements, and there was a statistically significant difference in mPAP and total PVR across the groups, we think that the effects of isoflurane can be neglected. Open-chest surgeries may produce an artifactual reduction of right ventricular pressure. 38 Pacher et al. have shown that it is feasible to obtain similar results in open- vs. closed-chest surgery, if body temperature, fluid balance, and anesthesia are adequately regulated and monitored during the surgical procedure, 5 which we carefully monitored in our experiments. Given that the study focused on RV mechanics and hemodynamics and relied on the literature for expected biological changes, no biological data are presented.
Although the RV wall is not a membrane and the thin-walled assumption is less well justified when the wall thickens in PAH, the wall thickness to radius ratio did not exceed 0.14, and the model was only used to estimate global pressure-volume relations not non-homogeneous myocardial stress distributions. To compute RV end-diastolic and end-systolic pressure-volume relations from RV geometry and myocardial resting and active material properties, a spherical model of the RV was adapted from the more comprehensive biventricular TriSeg model. 37 Like previous works,12,13 our model includes length-dependent calcium-sensitivity mechanics. However, we simplified the TriSeg model to exclude the LV, ventricular interactions, and circulatory dynamics because the current study was only focused on the effects of RV geometry and material properties on PV relations and because all the hemodynamic and structural data necessary to parameterize the more comprehensive model were not measured. In spite of these simplifications, the model was able to recapitulate the experimentally measured mean RV pressure-volume relations, using sarcomere length-tension relations similar to those reported in rat cardiac trabeculae by ter Keurs et al. 16
Conclusion
In conclusion, these results suggest the MCT-animal model at four weeks at 60 mg/kg is a model of RV compensation. At week 1, SV was maintained by the RV decreasing its diastolic chamber stiffness, not via changes in systolic function. However, by week 2 that mechanism was not enough: diastolic function had reversed closer to baseline and increased contractility was needed as the RV pressures started to rise. Changes in the EDPVR (i.e. decrease in passive tension) could also be an effect of monocrotaline as there was no change/increase in pressure. By week 2, however, the RV contractility increases (increase in active sarcomere tension) in order to maintain SV and CO. Interestingly, by weeks 3 and 4, the RV does not appear to need to recruit as much muscle fibers because the ventricle started to remodel. Hence, the effects of geometric remodeling on the EDPVR and of increased myocardial contractility on the ESPVR were the major mechanisms by which SV was preserved in the setting of elevated PAP, and both mechanisms contributed equally during weeks 2–4. Therefore, compensated right ventricular function depends on chamber remodeling and contractile augmentation in the onset of PAH in this rat animal model.
Footnotes
Conflict of interest
The author(s) declare that there is no conflict of interest.
Funding
This work was supported, in part, by grants from the US American Heart Association Scientist Development Grant 16SDG29670010 (DVJ), UIC College of Engineering Seed Fund (DVJ) and UIC Honors College Undergraduate Research Grant (JG).
