- 1Stanford Cardiovascular Institute, Stanford University, Stanford, CA, United States
- 2School of Biomedical Engineering, Colorado State University, Fort Collins, CO, United States
- 3Department of Mechanical Engineering, Colorado State University, Fort Collins, CO, United States
- 4Laboratory Animal Resources, Colorado State University, Fort Collins, CO, United States
- 5Department of Biomedical Sciences, Colorado State University, Fort Collins, CO, United States
Introduction: The right ventricle (RV) mechanical property is an important determinant of its function. However, compared to its elasticity, RV viscoelasticity is much less studied, and it remains unclear how pulmonary hypertension (PH) alters RV viscoelasticity. Our goal was to characterize the changes in RV free wall (RVFW) anisotropic viscoelastic properties with PH development and at varied heart rates.
Methods: PH was induced in rats by monocrotaline treatment, and the RV function was quantified by echocardiography. After euthanasia, equibiaxial stress relaxation tests were performed on RVFWs from healthy and PH rats at various strain-rates and strain levels, which recapitulate physiological deformations at varied heart rates (at rest and under acute stress) and diastole phases (at early and late filling), respectively.
Results and Discussion: We observed that PH increased RVFW viscoelasticity in both longitudinal (outflow tract) and circumferential directions. The tissue anisotropy was pronounced for the diseased RVs, not healthy RVs. We also examined the relative change of viscosity to elasticity by the damping capacity (ratio of dissipated energy to total energy), and we found that PH decreased RVFW damping capacity in both directions. The RV viscoelasticity was also differently altered from resting to acute stress conditions between the groups—the damping capacity was decreased only in the circumferential direction for healthy RVs, but it was reduced in both directions for diseased RVs. Lastly, we found some correlations between the damping capacity and RV function indices and there was no correlation between elasticity or viscosity and RV function. Thus, the RV damping capacity may be a better indicator of RV function than elasticity or viscosity alone. These novel findings on RV dynamic mechanical properties offer deeper insights into the role of RV biomechanics in the adaptation of RV to chronic pressure overload and acute stress.
1 Introduction
Right ventricle (RV) failure is a lethal condition that contributes significantly to the mortality and morbidity in a variety of cardiovascular diseases including pulmonary hypertension (PH), congenital heart disease, and left heart failure with preserved ejection fraction (HFpEF) (Voelkel et al., 2006; Haddad et al., 2008; Köhler et al., 2013; Konstam et al., 2018; Lahm et al., 2018). Unfortunately, there has been a lack of effective treatment for these patients, which is due in part to an inadequate understanding of the structure-function relationship of the right ventricle (RV) in physiological and pathological conditions (Konstam et al., 2018; Lahm et al., 2018). The biomechanical properties of the RV free wall (RVFW) are considered to impact its organ function (Jang et al., 2017; Liu, 2020; Liu et al., 2020), and a recent rodent study has reported a correlation between the passive elastic modulus of the RV and the end-diastolic volume (Jang et al., 2017), which is often used as an indicator of heart failure (Diebel et al., 1992; Durham et al., 1995; Cheatham et al., 1998; Zhang et al., 2014; Vonk Noordegraaf et al., 2017). Furthermore, our own ovine study showed that the RV elasticity is correlated with RV geometry and hemodynamic properties (Liu et al., 2020). These studies directly relate the RV tissue biomechanics to clinically relevant parameters, leading to a unique perspective to link the tissue biomechanics with organ function to advance the understanding of biomechanical mechanism of RV failure.
It is well known that the ventricular wall is an anisotropic and viscoelastic material (Dokos et al., 2002; Wang et al., 2016; Ahmad et al., 2018; Nguyen-Truong and Wang, 2018; Liu et al., 2021). Such behavior indicates that the tissue’s mechanical behavior is strain and strain-rate dependent, and there are energies stored and dissipated during the deformation due to the material’s elasticity and viscosity, respectively. These energy expenditure events closely affect the use of metabolic energy from cardiomyocytes during ventricular dilation and contraction in physiological conditions (e.g., at rest or under acute stress). However, most previous studies on RV biomechanics have focused on the elasticity of the tissue (Humphrey et al., 1990; Sacks and Chuong, 1993; Ghaemi et al., 2009; Holzapfel and Ogden, 2009; Valdez-Jasso et al., 2012; Hill et al., 2014; Sommer et al., 2015; Sirry et al., 2016). A few studies that reported the viscoelastic behavior of the RV have used an indentation test of the RVFW (Rubiano et al., 2016), the tensile mechanical test of neonatal RVFW (Ahmad et al., 2018) or the length-tension test of RV papillary muscles (Yamamoto et al., 1998; Stroud et al., 2002). Thus, the biaxial viscoelastic behavior of adult RVFW remains little investigated. Our group originally examined the biaxial viscoelasticity of RVFW in healthy adult sheep (Liu et al., 2021; 2022), but the changes in the viscoelastic behavior in pathological conditions are unclear. Furthermore, at the cell level or in papillary muscles, it has been shown that the myocardial viscoelasticity is increased in hypertrophy/diseased hearts and significantly weakens the muscle’s contractile function (Cooper, 2006; Caporizzo et al., 2018). Thus, the viscoelastic behavior of myocardium is critical to cardiac function. But how exactly the tissue-level viscoelasticity is altered in pathological remodeling and affects the myocardial function remains a knowledge gap.
The strain-rate dependent character of viscoelasticity determines that the ventricle wall’s mechanical behavior is different at different heart rates (HR), such as those at rest (normal HR) and under acute stress or exercise conditions (increased HR) (Schubert et al., 2009). Unlike the chronic stress that reduces the HR, acute stress caused by a psychological or physical stressor typically results in an elevated HR. The adaptation to acute stress involves changes in multiple cardiovascular parameters (e.g., cardiac output, blood pressure), which is often known as cardiovascular reactivity (Huang et al., 2013). But little is studied on the acute mechanical changes of the ventricle wall. The investigation of the RVFW viscoelastic changes under increased strain-rate (HR) will provide novel biomechanical insights for the RV adaptation to acute stress. Additionally, the healthy myocardium adapts to the increased workload during exercise by increasing stroke volume and heart rate, thus increasing the cardiac output—the blood supply to downstream organs. However, an impaired cardiac reserve has been noted in pulmonary arterial hypertension (PAH) and HFpEF patients, and the inability of the RV to increase cardiac output during exercise is one of key contributors (Borlaug et al., 2016; Guazzi et al., 2016; Lin et al., 2016; Jaijee et al., 2018; Malenfant et al., 2021). Whether the diseased (failing) RV has a distinct response to increased HR in its viscoelastic properties than the healthy RV is completely unknown, and the investigation of the viscoelastic changes will bring additional information to the capacity of failing RVs to adapt to the increased myocardial demand.
Therefore, the goal of this study was to investigate the alterations of RVFW biaxial viscoelastic behavior with PH development and at normal and increased heart rates simulating resting and acute stress conditions. We hypothesize that the remodeling in pulmonary hypertensive RVs leads to increased viscoelastic properties and tissue anisotropy and decreased damping capacity, and the increased HR causes a stronger reduction in damping capacity in the diseased RV than the healthy RV. This is the first study to investigate and compare the biaxial dynamic mechanical properties of the RVFW in various physiological and pathological states. The novel findings will improve our understanding of the RV tissue biomechanics in response to PH progression and increased heart rates and offer new insights into the adaptation of RV to chronic pressure overload and acute stress.
2 Materials and methods
2.1 Animal model and in vivo measurement
All procedures were approved by Institutional Animal Care and Use Committee (IACUC #1438) at Colorado State University. All animal experiments were performed in accordance with the guidelines and regulations of the Colorado State University Institutional Animal Care and Use Committee. All procedures were performed in accordance with the ARRIVE guidelines. Briefly, the monocrotaline (MCT) (60 mg/kg, Sigma) was injected once subcutaneously into 6-week-old male Sprague Dawley rats (Charles River) and the animals were housed in normal conditions for 3 weeks to develop PH. Similar healthy rats were used as controls (CTL). One or 2 days before euthanasia, the RV function was obtained by 2D echocardiography. Briefly, rats were anesthetized using inhaled isoflurane (0.5%–2%) and then subjected to transthoracic echocardiography using a Phillips HD11 ultrasound with a 15 mHz sector array probe. Right ventricular dimensions and contractile function including RV end-diastolic area (EDA), end-systolic area (ESA), and fractional shortening (FS) were measured in the apical four-chamber view in 2D B-mode. Pulmonary arterial (PA) outflow dynamics including peak outflow velocity (VPA peak), outflow acceleration time (AT), ejection time (ET), and peak late (atrial) RV filling velocity (VRV max) were recorded using pulse wave Doppler echocardiography in the parasternal long axis (PA outflow) or apical 4-chamber view (RV filling). Prior to the mechanical tests, rats were euthanized by exsanguination after an i.p. injection of urethane (1.2 kg/mg) or CO2 inhalation. After tissue harvest, RV hypertrophy was measured by the thickness of RVFW and the Fulton index, which is the ratio of wet tissue weights calculated as RV/(LV + S), where LV is the left ventricle free wall and S is the septum (Wang and Chesler, 2013; Wang et al., 2013).
2.2 Specimen preparation
Fresh hearts (N = 12 for CTL and N = 10 for MCT) were obtained from CTL and MCT rats. The hearts were placed in cardioplegic solution on ice to maintain tissue viability (Witzenburg et al., 2012) and mechanically tested within 3 h after tissue harvest. The RV was dissected, and the outflow tract direction was marked as the longitudinal direction and its perpendicular direction was defined as the circumferential direction. The samples were then placed in cardioplegic solution combined with 30 mM of 2, 3-butanedione monoxime (BDM) at body temperature for at least 30 min to ensure muscle relaxation. The tissues were immersed in the same solution during the mechanical tests.
2.3 Ex vivo biaxial stress relaxation tests
After mounting, the passive equibiaxial stress relaxation tests were performed after the tissues were preloaded and preconditioned by using a newly established in-house biaxial testing system (Roth et al., 2022). Two sets of stress relaxation tests were included in this study: 1) using different ramp speeds mimicking the sub-physiological and physiological heart rates of rats (0.1, 1, 2, 5 and 8 Hz) at the fixed strain level (20%). In this study, we treated the 0.1–2 Hz testing data as sub-physiological range and the 5 and 8 Hz testing data to represent the heart rates under resting and acute stress conditions, respectively. 2) using different strain levels (6, 9, 12% and 15%) at two fixed ramp speeds (5 and 8 Hz) to reveal the type of viscoelastic behavior (linear, quasilinear or nonlinear viscoelasticity) at rest and under acute stress. Between each test, a recovery time, which was ten times that of the previous testing period, was included to ensure that the tissue was fully recovered from the previous test. Biaxial stretch forces were obtained by 5 lb load cells (Honeywell) at a sampling frequency of 200 Hz, and the engineering stress (
The relaxation modulus was derived as the stress divided by the input strain at a fixed time of the relaxation (Liu et al., 2021) and used to indicate the elastic behavior. Furthermore, from the raw data, we derived the stored energy (Ws) as the product of the input strain and the area beneath the stress relaxation curve and dissipated energy (Wd) as the product of the input strain and the area between the stress relaxation curve and the theoretical curve of a purely elastic material at a fixed time of the relaxation (Figure 1). We then further derived the damping capacity as Wd/(Wd + Ws) to evaluate the relative change of viscosity to elasticity during tissue deformation. The stored and dissipated energies, damping capacity, and relaxation modulus were measured at five different time points (0.01, 0.1, 1, 10, and 100 s) to indicate the viscoelastic relaxation at different stages.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1182703/fbioe-11-1182703-HTML/image_m/fbioe-11-1182703-g001.jpg)
FIGURE 1. The derivations of stored energy and dissipated energy at any fixed time point use different areas calculated from the stress relaxation curve.
Furthermore, the logarithmic scale was plotted for the stress relaxation curves obtained from the Set 2 tests. The relaxation rates and the initial stress were then derived as the slopes and the ordinate intercept of the linear fitting to the logarithmic plots at each input strain level using least square analysis. The goodness of the fit was examined by the R2 value. The calculated relaxation rates from all strain levels were then used to determine the type of tissue viscoelastic behavior, where a fully nonlinear viscoelastic behavior is indicated if a strain-dependent behavior is observed (Troyer and Puttlitz, 2012; Liu et al., 2021; Liu et al., 2022). Lastly, the dependence of the relaxation rate or initial stress on the input strain was examined by curve fitting, and linear and nonlinear (quadratic) relations between these parameters and the input strain were determined as described previously (Shetye et al., 2014; Liu et al., 2022). The fit with a quadratic polynomial formular generated better fitting outcomes, which indicates that the RV exhibits a unique nonlinear viscoelastic behavior with a quadratic relation between relaxation rate or initial stress and the input strain.
2.4 Statistical and correlation analyses
Comparisons between groups (CTL and MCT), directions (longitudinal and circumferential), strain levels (6, 9, 12, and 15%), and ramp speed/frequencies (0.1, 1, 2, 5, and 8 Hz) were performed with the paired Student’s t-test or one-way ANOVA with Tukey’s post hoc analysis. Pearson correlation analysis was used to investigate the correlations between RV damping capacity and RV function measurement. The analyses were performed by Microsoft Excel and GraphPad Prism 9. Data are presented as mean ± SEM and p < 0.05 was considered statistically significant.
3 Results
3.1 RV failure establishment and structural changes in MCT RV
Transthoracic echocardiography revealed significant enlargement of RV chamber size (increased EDA and ESA) and reduced RV FS characteristic of RV dilation and systolic dysfunction, as well as marked changes in PA outflow and RV filling velocities that are characteristic of PH (lower AT/ET) and RV diastolic dysfunction (higher late VRV max) (Table 1). These results demonstrate the expected establishment of pathological RV remodeling in MCT rats. From the tissue measurement, we observed significant increases in the RVFW thickness and Fulton index in the MCT group compared to the CTL group, indicating RV hypertrophy in the MCT group (Table 1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1182703/fbioe-11-1182703-HTML/image_m/fbioe-11-1182703-t001.jpg)
TABLE 1. Summary of the functional and structural measurements of the RV for both control (CTL) and monocrotaline (MCT) groups. Data are present as mean ± SEM.
3.2 Stretch-rate (frequency) dependent viscoelastic behavior of RVFW altered by PH
We firstly investigated the frequency-dependent viscoelastic behavior obtained from the Set 1 tests. Figure 2 shows significant frequency-dependent changes in the relaxation modulus measured at 0.01 s after the peak force in both groups. Firstly, the relaxation modulus curves from the MCT group were above those from the CTL group, indicating increased RV viscoelasticity in the MCT group. We further present detailed changes in the physiologically relevant conditions in the next section. Second, the MCT group presented monotonically decreased relaxation modulus with the increasing frequency in both directions, whereas the CTL group presented weaker frequency-dependent behavior than the MCT group. Therefore, PH altered the frequency-dependent viscoelastic behavior of the RV. Similar frequency-dependent behavior was observed at other time points.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1182703/fbioe-11-1182703-HTML/image_m/fbioe-11-1182703-g002.jpg)
FIGURE 2. Frequency-dependent changes in the relaxation modulus at 0.01 s after peak force with an input strain of 20%. (A): relaxation modulus in the longitudinal (L) direction for both groups; (B): relaxation modulus in the circumferential direction for both groups. # <.05 vs. 0.1 Hz, $ <.05 vs. 1 Hz, % <.05 vs. 2 Hz, ^ <.05 vs. 5 Hz, respectively.
3.3 PH increased RVFW biaxial viscoelasticity and enhanced the tissue anisotropy
Next, we focused on the RV biaxial viscoelasticity obtained from 5 and 8 Hz with an input strain of 20% since these data are physiologically relevant. As shown in Figure 3A, we observed that the MCT RV had increased relaxation modulus than the CTL RV in both directions, and the MCT RV became anisotropic with a larger stiffness in the longitudinal than circumferential direction during PH development. Similar observation was found in the stored energy WS, another index of tissue elasticity (Figure 3B), and the dissipated energy Wd, an index of tissue viscosity (Figure 3C). The same changes in these viscoelastic parameters were also observed under acute stress. Moreover, like the elastic property, the viscous property was similar between directions in the CTL RV, suggesting an isotropic type of tissue. In the MCT RV, the viscous anisotropy was pronounced in the acute stress condition (p < 0.05 at 0.01, 10, and 100 s, p = 0.058 at 1 s) and was absent in the resting condition across all relaxation stages (p = 0.052 at 100 s, Figure 3C). We further derived various anisotropic indices (AI) as the ratio of elastic (or viscous) parameters between longitudinal and circumferential directions (Table 2). We observed that most of the AI values were larger than 1 (except for those derived from Wd in CTL group), indicating a stronger elastic or viscous behavior in the longitudinal direction, and that the AI values were also larger in the elastic than viscous parameters. Moreover, the increases in AI from healthy to diseased RVs was evident only in the viscous behavior. These data indicate that PH increases RVFW elasticity and viscosity and enhances tissue anisotropy. In the diseased RV, while the elastic anisotropy was significant at both resting and acute stress conditions, the viscous anisotropy was only evident under acute stress.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1182703/fbioe-11-1182703-HTML/image_m/fbioe-11-1182703-g003.jpg)
FIGURE 3. RVFW viscoelasticity before and after PH development and at resting and acute stress states. (A): RV elasticity measured by relaxation modulus at 0.01 s after peak force; (B): RV elasticity measured by WS at 100 s after peak force; (C) RV viscosity measured by Wd at 100 s after peak force. *p <.05.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1182703/fbioe-11-1182703-HTML/image_m/fbioe-11-1182703-t002.jpg)
TABLE 2. Summary of various RVFW anisotropic indices (AI) obtained from different elastic or viscous parameters in both groups and testing conditions. Data are present as mean ± SEM.
3.4 PH decreased RVFW damping capacity in both directions
We have shown previously that both elasticity and viscosity of the RV were increased with PH development, but it is unclear if these increases were equal or not. To evaluate the relative change of viscosity to elasticity, we examined the damping capacity as the percentage of dissipated (wasted) energy to the total energy during the relaxation. We found that in the longitudinal direction, the MCT RV had decreased damping capacity than the CTL RV in early relaxation (up to 1 s), and in the circumferential direction the MCT RV had reduced damping capacity in almost all stages of the relaxation (p <0.05 for all time points except for p = 0.09 for 0.01 s) (Table 3). We presented the results at 0.01 s and 100 s (equilibrium) of the relaxation obtained in normal heart rate condition in Figure 4. Similar changes were found in and the acute stress condition as well (Table 3). Therefore, this is the first report that the PH-induced RV remodeling led to reduced damping capacity of the RVFW. Lastly, we found that except for the early stage of the relaxation (at 0.01 s and 0.1 s), the RVFW damping capacity exhibited an anisotropic behavior in both rest and acute stress conditions (Figure 4B; Table 3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1182703/fbioe-11-1182703-HTML/image_m/fbioe-11-1182703-t003.jpg)
TABLE 3. Summary of the RVFW damping capacity obtained from different groups, directions, and testing conditions. Data are present as mean ± SEM.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1182703/fbioe-11-1182703-HTML/image_m/fbioe-11-1182703-g004.jpg)
FIGURE 4. The damping capacity of the RVFW before and after PH development quantified (A) at 0.01 s and (B) 100 s after peak force. **p <.01.
3.5 Acute stress reduced damping capacity differently between healthy and diseased groups
We further investigated the changes in the damping capacity from normal to high heart rates to reveal the viscoelastic alteration from resting to acute stress conditions. We observed that in the early stage of the relaxation (up to 1 s) there was no significant change in damping capacity in either group. At later stage of the relaxation (at 10 s or 100 s), however, the damping capacity was reduced in the tissue for both groups (Figure 5; Table 3). Moreover, the alteration was different between the CTL and MCT groups. The damping capacity was decreased only in the circumferential direction in the CTL group, whereas in the MCT group the decrease was observed in both directions and with small p values (Figure 5B). These data suggest that PH-induced remodeling causes a stronger reduction in damping capacity under acute stress condition than the healthy tissues.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1182703/fbioe-11-1182703-HTML/image_m/fbioe-11-1182703-g005.jpg)
FIGURE 5. Changes in RVFW damping capacity from the resting to acute stress states at (A) 0.01 s after peak force and (B) 100 s after peak force. * <.05, ** <.01, respectively.
3.6 Correlations between the RVFW viscoelastic properties and in vivo measurements
We further investigated the relationships between the RVFW viscoelastic parameters and the echocardiography measurements to explore potential implications of RVFW viscoelasticity. In both directions, we observed that the RVFW damping capacity tended to be strongly or was significantly correlated with the RV EDA, AT and AT/ET (Figure 6). More details of the correlations are shown in Supplementary Table S1. It is worthy of note that neither the elastic nor viscous parameter was correlated with any of these functional parameters. Our findings suggest that the damping capacity is a better indicator of RV function than the tissue elasticity or viscosity alone.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1182703/fbioe-11-1182703-HTML/image_m/fbioe-11-1182703-g006.jpg)
FIGURE 6. Correlations between the RVFW damping capacity at 10 s after peak force and (A) RV EDA, or (B) AT, or (C) ratio of AT and ET. The CTL data were included as the group average measurement (n = 8) and the error bar shows the standard error. The MCT data were included as individual data points (n = 5). L: longitudinal direction; (C) circumferential direction.
3.7 Altered strain-dependent viscoelastic behavior of the RVFW by PH or heart rate
From the second set of mechanical tests (i.e., stress relaxation test at varied strains), we investigated the RVFW viscoelastic relaxation altered by PH or heart rate. Specifically, we compared the rate of relaxation derived from the logarithmic plots of the stress relaxation curves. We observed significant strain-dependent behavior of the relaxation rate in both groups, and in both resting and acute stress conditions (Figure 7). This indicates that the RVFW exhibits a fully nonlinear viscoelastic (NLV) behavior in all groups and conditions. Moreover, the shape of the strain-dependent curve was different between healthy and diseased groups, between rest and acute stress conditions, as well as between directions at certain strain levels. This indicates that the RVFW has varied anisotropic NLV behaviors in different states (healthy vs. diseased, rest vs. acute stress). Therefore, it is necessary to measure the mechanical behavior of the tissue in the corresponding physiological condition.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1182703/fbioe-11-1182703-HTML/image_m/fbioe-11-1182703-g007.jpg)
FIGURE 7. Different strain-dependent and anisotropic behaviors measured by the relaxation rate in the (A,B) healthy (CTL) and (C,D) diseased (MCT) groups in both rest and acute stress conditions. # <.05 vs. 6%, $ <.05 vs. 9%, % <.05 vs. 12%, * <.05, respectively.
Furthermore, we investigated the dependence of the relaxation rate or initial stress on the input strain level using either linear or quadratic curve fitting. From the relations with a good fitting (R2 > 0.8), we found that the CTL and MCT groups had different fitting equations, indicating different NLV behaviors. Similarly, the fitting equations were different from resting to acute stress conditions, indicating that the acute stress altered the NLV behavior as well. Finally, different fitting equations were observed between directions in certain groups (e.g., the linear fit of the initial stress-strain relations at rest in the MCT group), which indicates an anisotropic NLV behavior in this group. Therefore, these data showed that the RVFW viscoelasticity exhibited different types of strain dependence in PH or acute stress conditions.
Lastly, we found that at rest, the linear relation between the relaxation rate and input strain provided good fitting to capture the NLV behavior for the CTL group (R2 > 0.9), but it led to a poor fitting for the MCT group (R2 = 0.29–0.42). Instead, a quadratic fitting was good to describe the relaxation rate-input strain relationship for the MCT group (R2 > 0.9) (Table 4). These data showed that at rest, the healthy RVFW viscoelasticity had a linear dependence on the strain, whereas the diseased RVFW viscoelasticity had a nonlinear (quadratic) dependence on strain. In acute stress condition, both linear and nonlinear relations provided sufficient fitting for the CTL and MCT groups (R2 > 0.8) (Table 4).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1182703/fbioe-11-1182703-HTML/image_m/fbioe-11-1182703-t004.jpg)
TABLE 4. The viscoelastic parameters (relaxation rate and initial stress) as a function of input strain level (
4 Discussion
To best of our knowledge, this is the first study that investigated the changes of RVFW biaxial viscoelasticity during PH progression. The main findings are: 1) PH increased RVFW viscoelasticity in both longitudinal (outflow tract) and circumferential directions, and the tissue anisotropy was pronounced for the diseased RV, not healthy RV. 2) PH decreased RVFW damping capacity (ratio of dissipated energy to total energy) in both directions. 3) The RVFW viscoelasticity was differently altered from resting to acute stress conditions between these groups—the damping capacity was decreased only in the circumferential direction for the healthy RV, but it was reduced in both directions for the diseased RV. 4) Some correlations were observed between the damping capacity and RV function indices. 5) the RVFW exhibited a fully NLV behavior, and this behavior was altered by disease progression (PH) or heart rate. These findings will improve the understanding of RV tissue dynamic mechanical behavior with PH development and under acute stress and shed light onto the biomechanical mechanism of RV failure and the reduced exercise capacity of failing RVs.
4.1 PH increased RVFW viscoelasticity and enhanced RVFW anisotropy
The hemodynamic measurement clearly demonstrated that RV failure was established in the MCT rats (Table 1). Signs of RV failure was further evidenced by the dilation and hypertrophy in the RVs of the MCT group (Table 1). With PH development, we found significant changes in the biaxial viscoelasticity of RVFW. In terms of the elasticity, we observed increases in elasticity in both directions, with more increase in the longitudinal than circumferential direction (Figures 3A, B). Thus, it led to an enhancement in the anisotropic, elastic behavior of the RV (p <0.05). In terms of the viscosity, we observed similar increases (Figure 3C) and an altered anisotropy (from more viscosity in circumferential direction to more viscosity in longitudinal direction, Table 2) during PH development. PH also led to a stronger anisotropy in elasticity from resting to acute stress condition (Figure 3; Table 2).
The increase in myocardium viscoelasticity during heart failure progression (including the RV failure) has been reported previously. Typically, 1D force/tension-length measurements have been performed in isolated cardiomyocytes or papillary muscles, and the viscoelastic property has been quantified from stress relaxation, cyclic stress-strain hysteresis, or the loading stress-strain curve (Dokos et al., 2002; Linke and Fernandez, 2002; Stroud et al., 2002; Cooper, 2006; Caporizzo et al., 2018; Caporizzo et al., 2020). Increases in elasticity and viscosity of the diseased myocardium were noted in these studies. But these methodologies assume an isotropic material of the myocardium and are mainly restrained at the cell level. Rubiano et al. (2016) induced systemic hypertension (not pulmonary hypertension) in adult rats and both left and right ventricles (LV and RV) were hypertrophied. They measured the ventricular tissue viscoelasticity by indentation tests (assuming isotropic material) and a standard linear viscoelasticity model. Interestingly, the RVFW showed increased elasticity and viscosity in the hypertensive rats, whereas the LVFW had only increased viscosity. Thus, our observation of increased elasticity and viscosity in the RVFW is consistent with prior findings in failing hearts. More importantly, the present study provides for the first time the tissue level, anisotropic viscoelastic findings of the RVFW and fills a critical gap of knowledge.
Not surprisingly, we observed increases in stiffness in both longitudinal and circumferential directions, and the increase was more significant in the longitudinal direction. This led to a change of tissue anisotropy and the RVFW became anisotropic (Figure 3). Similar observations have been reported previously when hyperelastic properties of the RVFW is examined by other groups. The enhanced stiffening in the longitudinal direction has been attributed to the re-orientation of the myofibers toward the long axis during PH progression, and the fiber re-orientation contributes to the changes in both RV diastolic and contractile function (Hill et al., 2014; Avazmohammadi et al., 2017a; Avazmohammadi et al., 2017b; Avazmohammadi et al., 2019). Therefore, our observation in RVFW elasticity agrees with prior findings. Moreover, we observed similar increases in viscosity (Wd) in the diseased RVFWs, and the viscous anisotropy was stronger in the acute stress than resting condition. These elastic and viscous changes in the diseased RVFW may be related to several different molecular mechanisms contributing to these mechanical properties. In the PH-induced remodeling, the RV hypertrophies (Kuroha et al., 1991; Heitmeier et al., 2020; Nguyen-Truong et al., 2020) and collagen accumulation occurs (Lahm et al., 2018; Bogaard and Voelkel, 2019; Poble et al., 2019; Nguyen-Truong et al., 2020). Both increases in intracellular microtubules (a cytoskeleton that confers myocyte mechanical strength) and extracellular collagen fibers are expected to stiffen the material property of the tissue in all directions. In addition, it is well accepted that the tissue is the stiffest in the main fiber direction. When the myofibers and collagen fibers get re-aligned and closer to the longitudinal (outflow tract) direction, the RVFW becomes more anisotropic and stiffer in this direction as well. But the molecular mechanism of viscosity is much less known. The myocardial tissue viscosity depends on the sources of frictional forces, which include diverse intra-and inter-molecular frictions—from cytoskeletons (e.g., microtubules) of cardiomyocytes to extracellular matrix proteins (e.g., collagen, proteoglycans) and even interstitial fluids. To the best of our knowledge, there is no known dominant “direction” of frictional sliding in aligned fibrous materials. Indeed, compared to the elastic anisotropy, viscous anisotropy is rarely reported in aligned biomaterials. Our tissue data here seem to suggest that there was weaker anisotropy in viscosity in the remodeled RV in the resting condition, and the increased heart rate (acute stress condition) enhanced the viscous anisotropy (Figure 3B). Since viscosity is related to the energy dissipation due to weak molecular bonds, we speculate that the dynamic distribution of external bonds between collagen fibers (e.g., inter-fiber crosslinking) and/or myofibers is responsible for the observations in the diseased RV. Future work should elaborate the exact molecular mechanisms of the viscosity from various sources such as fiber and non-fiber components.
Lastly, since both elasticity and viscosity were increased in diseased RV, we examined the relative change of viscosity to elasticity by quantifying the damping capacity (Wd/(Wd + WS)). We observed the reduction of damping capacity in both directions (Figure 4). But the reduction was significant only in the early relaxation (up to 0.1 s) in the longitudinal direction and was significant in the most stages of relaxation in the circumferential direction (Table 3). This is observed both in the resting and acute stress conditions. These data suggest that the changes of viscosity and elasticity are not “matched” and the increase in elasticity is more pronounced than that in viscosity in PH development. It is accepted that the stiffer myocardium results in elevated wall stress and it is detrimental to the cardiomyocytes. But the role of viscosity in myocardium physiology is less clear. On one hand, the higher viscosity of the diseased tissue could result in more energy waste; but on the other hand, it may protect the cells by damping the stress. Moreover, it has been shown from cell mechanotransduction studies that the viscous substrate reduces the cell’s focal adhesion area due to a dissipation (loss) of tractional force and cytoskeletal tension, thus producing enhanced cell spreading (morphology) and the associated mechanosensing pathway (Curtis and Russell, 2009; Ryan and O’Brien, 2015). Therefore, the concomitant changes in myocardial elasticity and viscosity may be driven to maintain a mechanical homeostasis for the cardiomyocytes. However, in the failing RVs, the homeostasis (unmatched elasticity and viscosity) cannot be preserved and such disturbance may be responsible for the cellular dysfunction and organ failure.
4.2 Altered response of the diseased RV to acute stress
In the present study, we exposed the RVFW to acute stress condition by increasing the stretch rate and examined the viscoelastic response of the tissue. Our results showed that when changed from the resting to acute stress conditions, the healthy RVFW had decreased damping capacity only in the circumferential direction, whereas the failing RVFW had decreased damping capacity in both directions (Figure 5). Therefore, the healthy RV may be better in handling acute stress in terms of maintaining the mechanical homeostasis. Moreover, the unchanged damping capacity in the longitudinal direction in the healthy RV suggests that the mechanical energy waste in the main blood flow direction is preserved. In contrast, the diseased RV had more significant reduction in damping capacity in both directions, suggesting the tissue has poor response to manage the acute stress. It has been noted recently that the heart failure patients, including the RV failure patients, have reduced exercise tolerance (Borlaug et al., 2016). That is, the exercise challenge imposes a negative impact on the RV function and the dysfunctional RV has poor cardiac reserve for exercise (Borlaug et al., 2016). In response to the acute stress such as exercise, the ventricle typically contracts faster and stronger, thus leading to higher heart rate and stroke volume to increase the blood supply and meet the elevated metabolic demand. However, a clinical study showed that in the RV dysfunction patients with reduced exercise capacity, despite an increased heart rate, the increase in cardiac output was limited due to a weaker contractility (Guazzi et al., 2016). Therefore, our finding offers a mechanical explanation of the impaired exercise capacity in failing RVs: the healthy RVs adapt to the acutely increased heart rate by preserving the mechanical homeostasis along the main blood flow direction, whereas the diseased RVs failed to do so. Our data suggest that the damping capacity may be a critical viscoelastic parameter of the RV biomechanics, and we should include this viscoelastic property of the myocardium into consideration to fully understand the role of RV biomechanics in maladaptive remodeling.
4.3 Altered strain-dependent viscoelastic behavior of RV by PH or acute stress
Linear, quasi-linear and fully nonlinear are the three main types of viscoelastic behavior. While all these three behaviors are time-dependent, the linear viscoelastic material exhibits linear elastic and linear viscous behavior, whereas the quasi-linear viscoelastic material exhibits nonlinear elastic and linear viscous behavior. The “linear viscous behavior” can be revealed by showing the same shape of relaxation or creep curves at different input strains. That is, the time-dependent behavior is strain-independent or the viscosity and elasticity are separable (Lakes, 1998; Troyer and Puttlitz, 2011). In contrast, in a fully nonlinear viscoelastic (NLV) material, both elasticity and viscosity are strain-dependent and they are not separable. Our results showed that the RV viscoelastic behavior was nonlinear viscoelastic (NLV) in different states (physiological or pathological) and at different heart rates (at rest and acute stress). This means that for both healthy and diseased RVs, the tissue’s elasticity and viscosity are not separable. Moreover, we observed different shapes of the strain-dependent relaxation rate curves between healthy and diseased groups, between the two axes and between resting and acute stress conditions (Figure 7). These differences indicate that PH or acute stress (increased heart rate) altered the RV anisotropic NLV behavior, and it is important to obtain the tissue’s dynamic mechanical behavior in the corresponding physiological condition.
4.4 Implications of RV viscoelasticity in RV function
In the present work, we performed some preliminary work to unravel the implications of the RV viscoelasticity in its physiological performance. The correlation between the passive elastic modulus of the RV in the longitudinal direction and the end-diastolic volume has been reported in a rodent study (Jang et al., 2017). Similarly, correlations between the RV longitudinal elasticity and RV end-diastolic diameter/area were found by our group in a recent large animal study (Liu et al., 2020). Compared to the RV elasticity, the RV viscoelasticity is much less studied. As a result, whether and how the RV viscoelasticity is associated with the RV function is still unknown. Our results showed some preliminary correlations between the RV damping capacity and the in vivo parameters obtained from echocardiography (Figure 6); however, these correlations were absent for the RV elasticity indices (relaxation modulus and stored energy) and nearly absent for the RV viscosity indices (dissipated energy and relaxation rate). The discrepant results may be because of the different types of mechanical tests and/or testing conditions. The previous studies derived the tissue elasticity from the loading curve of cyclic tensile mechanical tests and using a sub-physiological stretch rate (quasi-static) (Jang et al., 2017; Liu et al., 2020), whereas the RV elasticity within this study was derived from the stress relaxation with a ramp speed similar to the diastolic stretch rate in rat hearts. Nevertheless, our data highlight the necessity to obtain the dynamic mechanical behavior under physiological conditions when assessing RV performance.
Furthermore, we reported a strong trend of correlation between RV damping capacity and RV EDA. Increase in RV end-diastolic size is closely associated with disease severity and mortality of the RV failure patients (Prihadi et al., 2019; Lejeune et al., 2020; Ji et al., 2022). Thus, our result indicates that the RV damping capacity may potentially indicate the severity of the RV dysfunction. Our results also showed that the damping capacity was strongly correlated with the AT of pulmonary outflow. Preclinical studies of PH subjects has reported a correlation between the ratio of systolic pulmonary artery pressure to AT and pulmonary vascular resistance (afterload) (Thibault et al., 2010). Thus, AT and AT/ET are used to non-invasively estimate RV systolic pressure to indicate the severity of RV dysfunction. The additional correlations between AT or AT/ET and damping capacity found in this study suggest again the role of damping capacity in predicting the severity of RV dysfunction.
4.5 Limitations
All the data were obtained from male rat RVs, and there might be a sex difference in the RV viscoelastic behavior that awaits future research. In this study, we mainly used increased heart rate (thus higher stretch speed at diastole) to mimic the blood filling in the acute stress condition and characterize the passive mechanical properties. However, it is possible that the maximal strain is also increased during acute stress as a result of larger venous blood return. We assume the change in the maximal strain is not significantly as the prior study showed an increase of heart rate from resting to acute stress conditions but the same stroke volume in pulmonary arterial hypertension patients (Groepenhoff et al., 2010). Furthermore, the sample size for our correlation analysis is small, and we used the group average value for the CTL group to estimate the correlations. Future work should confirm the correlations with more samples and using individual values.
4.6 Conclusions
To the best of our knowledge, this is the first study on the changes of RV biaxial viscoelastic behavior during PH development. Our results from the ex vivo stress relaxation tests showed that PH increased RV viscoelasticity in both directions. These changes resulted in an anisotropic viscoelastic behavior of the diseased RV with marked elastic anisotropy in both resting and acute stress conditions and the viscous anisotropy only in acute stress condition. PH development also reduced the damping capacity (ratio of dissipated energy to total energy) of the tissue in both directions. In addition, the RV viscoelasticity was acutely altered from resting to acute stress conditions—the damping capacity was decreased only in the circumferential direction for healthy RVs, but it was decreased in both directions for diseased RVs, indicating an impaired response to acute stress. Lastly, the RVFW exhibited a fully NLV behavior and this behavior was altered by disease progression (PH) or heart rate. These novel findings improve our understanding of RV biomechanics in response to chronically elevated pulsatile mechanical loadings or acutely increased heart rates, and they will shed light on the mechanical mechanism of RV failure and the reduced exercise capacity.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Ethics statement
The animal study was reviewed and approved by Institutional Animal Care and Use Committee (IACUC #1438) at Colorado State University. All animal experiments were performed in accordance with the guidelines and regulations of the Colorado State University IACUC. All procedures were performed in accordance with the ARRIVE guidelines.
Author contributions
WL and ZW contributed to conception and design of the study. KL, KR, JP, JA, AC, and CP assisted with data collection, analysis and/or interpretation. WL and KL organized the database. WL and ZW did data curation. WL performed the statistical analysis. WL wrote the first draft of the manuscript. WL and ZW did supervision. All authors contributed to the article and approved the submitted version.
Acknowledgments
We thank Matt Ahern and Erith Evans from the Colorado State University for the insightful discussions and animal-related work.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fbioe.2023.1182703/full#supplementary-material
References
Ahmad, F., Prabhu, R. j., Liao, J., Soe, S., Jones, M. D., Miller, J., et al. (2018). Biomechanical properties and microstructure of neonatal porcine ventricles. J. Mech. Behav. Biomed. Mat. 88, 18–28. doi:10.1016/j.jmbbm.2018.07.038
Avazmohammadi, R., Hill, M. R., Simon, M. A., Zhang, W., and Sacks, M. S. (2017a). A novel constitutive model for passive right ventricular myocardium: Evidence for myofiber-collagen fiber mechanical coupling. Biomech. Model. Mechanobiol. 16, 561–581. doi:10.1007/s10237-016-0837-7
Avazmohammadi, R., Hill, M., Simon, M., and Sacks, M. (2017b). Transmural remodeling of right ventricular myocardium in response to pulmonary arterial hypertension. Apl. Bioeng. 1, 16105. doi:10.1063/1.5011639
Avazmohammadi, R., Mendiola, E. A., Soares, J. S., Li, D. S., Chen, Z., Merchant, S., et al. (2019). A computational cardiac model for the adaptation to pulmonary arterial hypertension in the rat. Ann. Biomed. Eng. 47, 138–153. doi:10.1007/s10439-018-02130-y
Bogaard, H. J., and Voelkel, N. F. (2019). Is myocardial fibrosis impairing right heart function? Am. J. Respir. Crit. Care Med. 199, 1458–1459. doi:10.1164/rccm.201812-2307ED
Borlaug, B. A., Kane, G. C., Melenovsky, V., and Olson, T. P. (2016). Abnormal right ventricular-pulmonary artery coupling with exercise in heart failure with preserved ejection fraction. Eur. Heart J. 37, 3294–3302. doi:10.1093/eurheartj/ehw241
Caporizzo, M. A., Chen, C. Y., Bedi, K., Margulies, K. B., and Prosser, B. L. (2020). Microtubules increase diastolic stiffness in failing human cardiomyocytes and myocardium. Circulation 2020, 902–915. doi:10.1161/CIRCULATIONAHA.119.043930
Caporizzo, M. A., Chen, C. Y., Salomon, A. K., Margulies, K. B., and Prosser, B. L. (2018). Microtubules provide a viscoelastic resistance to myocyte motion. Biophys. J. 115, 1796–1807. doi:10.1016/j.bpj.2018.09.019
Cheatham, M. L., Nelson, L. D., Chang, M. C., and Safcsak, K. (1998). Right ventricular end-diastolic volume index as a predictor of preload status in patients on positive end-expiratory pressure. Crit. Care Med. 26, 1801–1806. doi:10.1097/00003246-199811000-00017
Cooper, G. (2006). Cytoskeletal networks and the regulation of cardiac contractility: Microtubules, hypertrophy, and cardiac dysfunction. Am. J. Physiol. Heart Circ. Physiol. 291, H1003–H1014. doi:10.1152/ajpheart.00132.2006
Curtis, M. W., and Russell, B. (2009). Cardiac tissue engineering. J. Cardiovasc. Nurs. 24, 87–92. doi:10.1097/01.JCN.0000343562.06614.49
Diebel, L. N., Wilson, R. F., Tagett, M. G., and Kline, R. A. (1992). End-diastolic volume: A better indicator of preload in the critically ill. JAMA Surg. 127, 817–822. doi:10.1001/archsurg.1992.01420070081015
Dokos, S., Smaill, B. H., Young, A. A., Legrice, I. J., and Legrice, I. J. (2002). Shear properties of passive ventricular myocardium. Am. J. Physiol. Hear. Circ. Physiol. 283, 2650–2659. doi:10.1152/ajpheart.00111.2002
Durham, R., Neunaber, K., Vogler, G., Shapiro, M., and Mazuski, J. (1995). Right ventricular end-diastolic volume as a measure of preload. J. Trauma Acute Care Surg. 39, 218–223. doi:10.1097/00005373-199508000-00006
Ghaemi, H., Behdinan, K., and Spence, A. D. (2009). In vitro technique in estimation of passive mechanical properties of bovine heart: Part I. Experimental techniques and data. Med. Eng. Phys. 31, 76–82. doi:10.1016/j.medengphy.2008.04.008
Groepenhoff, H., Westerhof, N., Jacobs, W., Boonstra, A., Postmus, P. E., and Vonk-Noordegraaf, A. (2010). Exercise stroke volume and heart rate response differ in right and left heart failure. Eur. J. Heart Fail. 12, 716–720. doi:10.1093/eurjhf/hfq062
Guazzi, M., Villani, S., Generati, G., Ferraro, O. E., Pellegrino, M., Alfonzetti, E., et al. (2016). Right ventricular contractile reserve and pulmonary circulation uncoupling during exercise challenge in heart failure: Pathophysiology and clinical phenotypes. JACC Hear. Fail. 4, 625–635. doi:10.1016/j.jchf.2016.03.007
Haddad, F., Hunt, S. A., Rosenthal, D. N., and Murphy, D. J. (2008). Right ventricular function in cardiovascular disease, part I: Anatomy, physiology, aging, and functional assessment of the right ventricle. Circulation 117, 1436–1448. doi:10.1161/CIRCULATIONAHA.107.653576
Heitmeier, T., Sydykov, A., Lukas, C., Vroom, C., Korfei, M., Petrovic, A., et al. (2020). Altered proteasome function in right ventricular hypertrophy. Cardiovasc. Res. 116, 406–415. doi:10.1093/cvr/cvz103
Hill, M. R., Simon, M. A., Valdez-Jasso, D., Zhang, W., Champion, H. C., and Sacks, M. S. (2014). Structural and mechanical adaptations of right ventricle free wall myocardium to pressure overload. Ann. Biomed. Eng. 42, 2451–2465. doi:10.1007/s10439-014-1096-3
Holzapfel, G. A., and Ogden, R. W. (2009). Constitutive modelling of passive myocardium: A structurally based framework for material characterization. Philos. Trans. R. Soc. A Math. Phys. Eng. Sci. 367, 3445–3475. doi:10.1098/rsta.2009.0091
Huang, C.-J., Webb, H., Zourdos, M., and Acevedo, E. (2013). Cardiovascular reactivity, stress, and physical activity. Front. Physiol. 4. doi:10.3389/fphys.2013.00314
Humphrey, J. D., Strumpf, R. K., and Yin, F. C. (1990). Biaxial mechanical behavior of excised ventricular epicardium. Am. J. Physiol. Circ. Physiol. 259, H101–H108. doi:10.1152/ajpheart.1990.259.1.H101
Jaijee, S., Quinlan, M., Tokarczuk, P., Clemence, M., Howard, L. S. G. E., Gibbs, J. S. R., et al. (2018). Exercise cardiac MRI unmasks right ventricular dysfunction in acute hypoxia and chronic pulmonary arterial hypertension. Am. J. Physiol. Heart Circ. Physiol. 315, H950–H957. doi:10.1152/ajpheart.00146.2018
Jang, S., Vanderpool, R. R., Avazmohammadi, R., Lapshin, E., Bachman, T. N., Sacks, M., et al. (2017). Biomechanical and hemodynamic measures of right ventricular diastolic function: Translating tissue biomechanics to clinical relevance. J. Am. Hear. Assoc. Cardiovasc. Cerebrovasc. Dis. 6, e006084. doi:10.1161/JAHA.117.006084
Ji, M., Zhang, Y., Lin, Y., Gao, L., Xie, M., and Li, Y. (2022). Research progress of right ventricular longitudinal strain in patients with left heart failure. Chin. J. Ultrason. 31, 365–369. doi:10.3760/cma.j.cn131148-20210821-00580
Köhler, D., Arnold, R., Loukanov, T., and Gorenflo, M. (2013). Right ventricular failure and pathobiology in patients with congenital heart disease – implications for long-term follow-up. Front. Pediatr. 1, 37. doi:10.3389/fped.2013.00037
Konstam, M. A., Kiernan, M. S., Bernstein, D., Bozkurt, B., Jacob, M., Kapur, N. K., et al. (2018). On behalf of the American heart association coun-cil on clinical cardiology; council on cardiovas-cular disease in the young; and council on cardiovascular surgery and anesthesia. Circulation 137, 578–622. doi:10.1161/CIR.0000000000000560
Kuroha, M., Isoyama, S., Ito, N., and Takishima, T. (1991). Effects of age on right ventricular hypertrophic response to pressure-overload in rats. J. Mol. Cell. Cardiol. 23, 1177–1190. doi:10.1016/0022-2828(91)90206-2
Lahm, T., Douglas, I. S., Archer, S. L., Bogaard, H. J., Chesler, N. C., Haddad, F., et al. (2018). Assessment of right ventricular function in the research setting: Knowledge gaps and pathways forward an official American thoracic society research statement. Am. J. Respir. Crit. Care Med. 198, e15–e43. doi:10.1164/rccm.201806-1160ST
Lejeune, S., Roy, C., Ciocea, V., Slimani, A., de Meester, C., Amzulescu, M., et al. (2020). Right ventricular global longitudinal strain and outcomes in heart failure with preserved ejection fraction. J. Am. Soc. Echocardiogr. Off. Publ. Am. Soc. Echocardiogr. 33, 973–984.e2. doi:10.1016/j.echo.2020.02.016
Lin, A. C. W., Strugnell, W. E., Seale, H., Schmitt, B., Schmidt, M., Rourke, R., et al. (2016). Exercise cardiac MRI-derived right ventriculo-arterial coupling ratio detects early right ventricular maladaptation in PAH. Eur. Respir. J. 48, 1797–1800. doi:10.1183/13993003.01145-2016
Linke, W. A., and Fernandez, J. M. (2002). Cardiac titin: Molecular basis of elasticity and cellular contribution to elastic and viscous stiffness components in myocardium. J. Muscle Res. Cell Motil. 23, 483–497. doi:10.1023/a:1023462507254
Liu, W. (2020). Current understanding of the biomechanics of ventricular tissues in heart failure. Bioengineering 7, 2. doi:10.3390/bioengineering7010002
Liu, W., Labus, K. M., Ahern, M., LeBar, K., Avazmohammadi, R., Puttlitz, C. M., et al. (2022). Strain-dependent stress relaxation behavior of healthy right ventricular free wall. Acta Biomater. 152, 290–299. doi:10.1016/j.actbio.2022.08.043
Liu, W., Nguyen-Truong, M., Ahern, M., Labus, K., Puttlitz, C., and Wang, Z. (2021). Different passive viscoelastic properties between the left and right ventricles in healthy adult ovine. J. Biomech. Eng. 143, 1–9. doi:10.1115/1.4052004
Liu, W., Nguyen-Truong, M., Labus, K., Boon, J., Easley, J., Monnet, E., et al. (2020). Correlations between the right ventricular passive elasticity and organ function in adult ovine. J. Integr. Cardiol. 6, 1–6. doi:10.15761/JIC.1000294
Malenfant, S., Lebret, M., Breton-Gagnon, É., Potus, F., Paulin, R., Bonnet, S., et al. (2021). Exercise intolerance in pulmonary arterial hypertension: Insight into central and peripheral pathophysiological mechanisms. Eur. Respir. Rev. 30, 284. doi:10.1183/16000617.0284-2020
Nguyen-Truong, M., and Wang, Z. (2018). Biomechanical properties and mechanobiology of cardiac ECM. Adv. Exp. Med. Biol. 1098, 1–19. doi:10.1007/978-3-319-97421-7_1
Nguyen-Truong, M., Liu, W., Boon, J., Nelson, B., Easley, J., Monnet, E., et al. (2020). Establishment of adult right ventricle failure in ovine using a graded, animal-specific pulmonary artery constriction model. Anim. Model. Exp. Med. 2020, 1–11. doi:10.1002/ame2.12124
Poble, P.-B., Phan, C., Quatremare, T., Bordenave, J., Thuillet, R., Cumont, A., et al. (2019). Therapeutic effect of pirfenidone in the sugen/hypoxia rat model of severe pulmonary hypertension. FASEB J. Off. Publ. Fed. Am. Soc. Exp. Biol. 33, 3670–3679. doi:10.1096/fj.201801659R
Prihadi, E. A., van der Bijl, P., Dietz, M., Abou, R., Vollema, E. M., Marsan, N. A., et al. (2019). Prognostic implications of right ventricular free wall longitudinal strain in patients with significant functional tricuspid regurgitation. Circ. Cardiovasc. Imaging 12, e008666. doi:10.1161/CIRCIMAGING.118.008666
Roth, K., Liu, W., LeBar, K., Ahern, M., and Wang, Z. (2022). A new high-speed biaxial testing system for rodent myocardium viscoelastic measurement under dynamic loadings. FASEB J. 36, 5183. doi:10.1096/fasebj.2022.36.S1.R5183
Rubiano, A., Qi, Y., Guzzo, D., Rowe, K., Pepine, C., and Simmons, C. (2016). Stem cell therapy restores viscoelastic properties of myocardium in rat model of hypertension. J. Mech. Behav. Biomed. Mat. 59, 71–77. doi:10.1016/j.jmbbm.2015.11.041
Ryan, A. J., and O’Brien, F. J. (2015). Insoluble elastin reduces collagen scaffold stiffness, improves viscoelastic properties, and induces a contractile phenotype in smooth muscle cells. Biomaterials 73, 296–307. doi:10.1016/j.biomaterials.2015.09.003
Sacks, M. S., and Chuong, C. J. (1993). Biaxial mechanical properties of passive right ventricular free wall myocardium. J. Biomech. Eng. 115, 202–205. doi:10.1115/1.2894122
Schubert, C., Lambertz, M., Nelesen, R. A., Bardwell, W., Choi, J.-B., and Dimsdale, J. E. (2009). Effects of stress on heart rate complexity--a comparison between short-term and chronic stress. Biol. Psychol. 80, 325–332. doi:10.1016/j.biopsycho.2008.11.005
Shetye, S. S., Troyer, K. L., Streijger, F., Lee, J. H. T., Kwon, B. K., Cripton, P. A., et al. (2014). Nonlinear viscoelastic characterization of the porcine spinal cord. Acta Biomater. 10, 792–797. doi:10.1016/j.actbio.2013.10.038
Sirry, M. S., Butler, J. R., Patnaik, S. S., Brazile, B., Bertucci, R., Claude, A., et al. (2016). Characterisation of the mechanical properties of infarcted myocardium in the rat under biaxial tension and uniaxial compression. J. Mech. Behav. Biomed. Mat. 63, 252–264. doi:10.1016/j.jmbbm.2016.06.029
Sommer, G., Schriefl, A. J., Andrä, M., Sacherer, M., Viertler, C., Wolinski, H., et al. (2015). Biomechanical properties and microstructure of human ventricular myocardium. Acta Biomater. 24, 172–192. doi:10.1016/j.actbio.2015.06.031
Stroud, J. D., Baicu, C. F., Barnes, M. A., Spinale, F. G., and Zile, M. R. (2002). Viscoelastic properties of pressure overload hypertrophied myocardium: Effect of serine protease treatment. Am. J. Physiol. - Hear. Circ. Physiol. 282, 2324–2335. doi:10.1152/ajpheart.00711.2001
Thibault, H. B., Kurtz, B., Raher, M. J., Shaik, R. S., Waxman, A., Derumeaux, G., et al. (2010). Noninvasive assessment of murine pulmonary arterial pressure: Validation and application to models of pulmonary hypertension. Circ. Cardiovasc. Imaging 3, 157–163. doi:10.1161/CIRCIMAGING.109.887109
Troyer, K. L., and Puttlitz, C. M. (2011). Human cervical spine ligaments exhibit fully nonlinear viscoelastic behavior. Acta Biomater. 7, 700–709. doi:10.1016/j.actbio.2010.09.003
Troyer, K. L., and Puttlitz, C. M. (2012). Nonlinear viscoelasticty plays an essential role in the functional behavior of spinal ligaments. J. Biomech. 45, 684–691. doi:10.1016/j.jbiomech.2011.12.009
Valdez-Jasso, D., Simon, M. A., Champion, H. C., and Sacks, M. S. (2012). A murine experimental model for the mechanical behaviour of viable right-ventricular myocardium. J. Physiol. 590, 4571–4584. doi:10.1113/jphysiol.2012.233015
Voelkel, N. F., Quaife, R. A., Leinwand, L. A., Barst, R. J., Mcgoon, M. D., Meldrum, D. R., et al. (2006). Right ventricular function and failure report of a national heart, lung, and blood institute working group on cellular and molecular mechanisms of right heart failure the normal right ventricle the right ventricle in pulmonary hypertension special report. Circulation 114, 1883–1891. doi:10.1161/CIRCULATIONAHA.106.632208
Vonk Noordegraaf, A., Westerhof, B. E., and Westerhof, N. (2017). The relationship between the right ventricle and its load in pulmonary hypertension. J. Am. Coll. Cardiol. 69, 236–243. doi:10.1016/j.jacc.2016.10.047
Wang, Z., and Chesler, N. C. (2013). Pulmonary vascular mechanics: Important contributors to the increased right ventricular afterload of pulmonary hypertension. Exp. Physiol. 98, 1267–1273. doi:10.1113/expphysiol.2012.069096
Wang, Z., Golob, M. J., and Chesler, N. C. (2016). “Viscoelastic properties of cardiovascular tissues,” in Viscoelastic and viscoplastic materials. Editor M. J. Golob (ijeka: IntechOpen). doi:10.5772/64169
Wang, Z., Schreier, D. A., Hacker, T. A., and Chesler, N. C. (2013). Progressive right ventricular functional and structural changes in a mouse model of pulmonary arterial hypertension. Physiol. Rep. 1, e00184. doi:10.1002/phy2.184
Witzenburg, C., Raghupathy, R., Kren, S. M., Taylor, D. A., and Barocas, V. H. (2012). Mechanical changes in the rat right ventricle with decellularization. J. Biomech. 45, 842–849. doi:10.1016/j.jbiomech.2011.11.025
Yamamoto, S., Tsutsui, H., Takahashi, M., Ishibashi, Y., Tagawa, H., Imanaka-Yoshida, K., et al. (1998). Role of microtubules in the viscoelastic properties of isolated cardiac muscle. J. Mol. Cell. Cardiol. 30, 1841–1853. doi:10.1006/jmcc.1998.0747
Zhang, S., Yang, Z.-G., Sun, J.-Y., Wen, L.-Y., Xu, H.-Y., Zhang, G., et al. (2014). Assessing right ventricular function in patients with hypertrophic cardiomyopathy with cardiac MRI: Correlation with the New York heart function assessment (NYHA) classification. PLoS One 9, e104312. doi:10.1371/journal.pone.0104312
Keywords: pressure overload, RV failure, anisotropy, stored energy, dissipated energy, exercise capacity
Citation: Liu W, LeBar K, Roth K, Pang J, Ayers J, Chicco AJ, Puttlitz CM and Wang Z (2023) Alterations of biaxial viscoelastic properties of the right ventricle in pulmonary hypertension development in rest and acute stress conditions. Front. Bioeng. Biotechnol. 11:1182703. doi: 10.3389/fbioe.2023.1182703
Received: 09 March 2023; Accepted: 17 May 2023;
Published: 30 May 2023.
Edited by:
Seungik Baek, Michigan State University, United StatesReviewed by:
Lei Fan, Michigan State University, United StatesChung-Hao Lee, University of Oklahoma, United States
Copyright © 2023 Liu, LeBar, Roth, Pang, Ayers, Chicco, Puttlitz and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zhijie Wang, zhijie.wang@colostate.edu