Transmural remodeling of right ventricular myocardium in response to pulmonary arterial hypertension

Simulated Effects of Acute Left Ventricular Myocardial Infarction on Mitral Regurgitation in an Ovine Model

Ischemic mitral regurgitation (IMR) occurs from incomplete coaptation of the mitral valve (MV) after myocardial infarction (MI), typically worsened by continued remodeling of the left ventricular (LV). The importance of LV remodeling is clear as IMR is induced by the post-MI dual mechanisms of mitral annular dilation and leaflet tethering from papillary muscle (PM) distension via the MV chordae tendineae (MVCT). However, the detailed etiology of IMR remains poorly understood, in large part due to the complex interactions of the MV and the post-MI LV remodeling processes. Given the patient-specific anatomical complexities of the IMR disease processes, simulation-based approaches represent an ideal approach to improve our understanding of this deadly disease. However, development of patient-specific models of left ventricle-mitral valve (LV-MV) interactions in IMR are complicated by the substantial variability and complexity of the MR etiology itself, making it difficult to extract underlying mechanisms from clinical data alone. To address these shortcomings, we developed a detailed ovine LV-MV finite element (FE) model based on extant comprehensive ovine experimental data. First, an extant ovine LV FE model (Sci. Rep. 2021 Jun 29;11(1):13466) was extended to incorporate the MV using a high fidelity ovine in vivo derived MV leaflet geometry. As it is not currently possible to image the MVCT in vivo, a functionally equivalent MVCT network was developed to create the final LV-MV model. Interestingly, in pilot studies, the MV leaflet strains did not agree well with known in vivo MV leaflet strain fields. We then incorporated previously reported MV leaflet prestrains (J. Biomech. Eng. 2023 Nov 1;145(11):111002) in the simulations. The resulting LV-MV model produced excellent agreement with the known in vivo ovine MV leaflet strains and deformed shapes in the normal state. We then simulated the effects of regional acute infarctions of varying sizes and anatomical locations by shutting down the local myocardial contractility. The remaining healthy (noninfarcted) myocardium mechanical behaviors were maintained, but allowed to adjust their active contractile patterns to maintain the prescribed pressure-volume loop behaviors in the acute post-MI state. For all cases studied, the LV-MV simulation demonstrated excellent agreement with known LV and MV in vivo strains and MV regurgitation orifice areas. Infarct location was shown to play a critical role in resultant MV leaflet strain fields. Specifically, extensional deformations of the posterior leaflets occurred in the posterobasal and laterobasal infarcts, while compressive deformations of the anterior leaflet were observed in the anterobasal infarct. Moreover, the simulated posterobasal infarct induced the largest MV regurgitation orifice area, consistent with experimental observations. The present study is the first detailed LV-MV simulation that reveals the important role of MV leaflet prestrain and functionally equivalent MVCT for accurate predictions of LV-MV interactions. Importantly, the current study further underscored simulation-based methods in understanding MV function as an integral part of the LV.

Dissecting contributions of pulmonary arterial remodeling to right ventricular afterload in pulmonary hypertension

Pulmonary hypertension (PH) is defined as an elevation in the right ventricle (RV) afterload, characterized by increased hemodynamic pressure in the main pulmonary artery (PA). Elevations in RV afterload increase RV wall stress, resulting in RV remodeling and potentially RV failure. From a biomechanical standpoint, the primary drivers for RV afterload elevations include increases in pulmonary vascular resistance (PVR) in the distal vasculature and decreases in vessel compliance in the proximal PA. However, the individual contributions of the various vascular remodeling events toward the progression of PA pressure elevations and altered vascular hemodynamics remain elusive. In this study, we used a subject-specific one-dimensional (1D) fluid-structure interaction (FSI) model to investigate the alteration of pulmonary hemodynamics in PH and to quantify the contributions of vascular stiffening and increased resistance towards increased main pulmonary artery (MPA) pressure. We used a combination of subject-specific hemodynamic measurements, ex-vivo mechanical testing of arterial tissue specimens, and ex-vivo X-ray micro-tomography imaging to develop the 1D-FSI model and dissect the contribution of PA remodeling events towards alterations in the MPA pressure waveform. Both the amplitude and pulsatility of the MPA pressure waveform were analyzed. Our results indicated that increased distal resistance has the greatest effect on the increase in maximum MPA pressure, while increased stiffness caused significant elevations in the characteristic impedance. The method presented in this study will serve as an essential step toward understanding the complex interplay between PA remodeling events that leads to the most severe adverse effect on RV dysfunction.

Role of Microtubule Network in the Passive Anisotropic Viscoelasticity of Healthy Right Ventricle

Cardiomyocytes are viscoelastic and key determinants of right ventricle (RV) mechanics. Intracellularly, microtubules are found to impact the viscoelasticity of isolated cardiomyocytes or trabeculae; whether they contribute to the tissue-level viscoelasticity is unknown. Our goal was to reveal the role of the microtubule network in the passive anisotropic viscoelasticity of the healthy RV. Equibiaxial stress relaxation tests were conducted in healthy RV free wall (RVFW) under early (6%) and end (15%) diastolic strain levels, and at sub- and physiological stretch rates. The viscoelasticity was assessed at baseline and after the removal of microtubule network. Furthermore, a quasi-linear viscoelastic (QLV) model was applied to delineate the contribution of microtubules to the relaxation behavior of RVFW. After removing the microtubule network, RVFW elasticity and viscosity were reduced at the early diastolic strain level and in both directions. The reduction in elasticity was stronger in the longitudinal direction, whereas the degree of changes in viscosity were equivalent between directions. There was insignificant change in RVFW viscoelasticity at late diastolic strain level. Finally, the modeling showed that the tissue's relaxation strength was reduced by the removal of the microtubule network, but the change was present only at a later time scale. These new findings suggest a critical role of cytoskeleton filaments in RVFW passive mechanics in physiological conditions.

Right ventricular stiffening and anisotropy alterations in pulmonary hypertension: Mechanisms and relations to function

Aims

Pulmonary hypertension (PH) results in an increase in RV afterload, leading to RV dysfunction and failure. The mechanisms underlying maladaptive RV remodeling are poorly understood. In this study, we investigated the multiscale and mechanistic nature of RV free wall (RVFW) biomechanical remodeling and its correlations with RV function adaptations.

Methods and Results

Mild and severe models of PH, consisting of hypoxia (Hx) model in Sprague-Dawley (SD) rats (n=6 each, Control and PH) and Sugen-hypoxia (SuHx) model in Fischer (CDF) rats (n=6 each, Control and PH), were used. Organ-level function and tissue-level stiffness and microstructure were quantified through in-vivo and ex-vivo measures, respectively. Multiscale analysis was used to determine the association between fiber-level remodeling, tissue-level stiffening, and organ-level dysfunction. Animal models with different PH severity provided a wide range of RVFW stiffening and anisotropy alterations in PH. Decreased RV-pulmonary artery (PA) coupling correlated strongly with stiffening but showed a weaker association with the loss of RVFW anisotropy. Machine learning classification identified the range of adaptive and maladaptive RVFW stiffening. Multiscale modeling revealed that increased collagen fiber tautness was a key remodeling mechanism that differentiated severe from mild stiffening. Myofiber orientation analysis indicated a shift away from the predominantly circumferential fibers observed in healthy RVFW specimens, leading to a significant loss of tissue anisotropy.

Conclusion

Multiscale biomechanical analysis indicated that although hypertrophy and fibrosis occur in both mild and severe PH, certain fiber-level remodeling events, including increased tautness in the newly deposited collagen fibers and significant reorientations of myofibers, contributed to excessive biomechanical maladaptation of the RVFW leading to severe RV-PA uncoupling. Collagen fiber remodeling and the loss of tissue anisotropy can provide an improved understanding of the transition from adaptive to maladaptive remodeling.

Translational perspective

Right ventricular (RV) failure is a leading cause of mortality in patients with pulmonary hypertension (PH). RV diastolic and systolic impairments are evident in PH patients. Stiffening of the RV wall tissue and changes in the wall anisotropy are expected to be major contributors to both impairments. Global assessments of the RV function remain inadequate in identifying patients with maladaptive RV wall remodeling primarily due to their confounded and weak representation of RV fiber and tissue remodeling events. This study provides novel insights into the underlying mechanisms of RV biomechanical remodeling and identifies the adaptive-to-maladaptive transition across the RV biomechanics-function spectrum. Our analysis dissecting the contribution of different RV wall remodeling events to RV dysfunction determines the most adverse fiber-level remodeling to RV dysfunction as new therapeutic targets to curtail RV maladaptation and, in turn, RV failure in PH.

A computational study of right ventricular mechanics in a rat model of pulmonary arterial hypertension

Pulmonary arterial hypertension (PAH) presents a significant challenge to right ventricular (RV) function due to progressive pressure overload, necessitating adaptive remodeling in the form of increased wall thickness, enhanced myocardial contractility and stiffness to maintain cardiac performance. However, the impact of these remodeling mechanisms on RV mechanics in not clearly understood. In addition, there is a lack of quantitative understanding of how each mechanism individually influences RV mechanics. Utilizing experimental data from a rat model of PAH at three distinct time points, we developed biventricular finite element models to investigate how RV stress and strain evolved with PAH progression. The finite element models were fitted to hemodynamic and morphological data to represent different disease stages and used to analyze the impact of RV remodeling as well as the altered RV pressure. Furthermore, we performed a number of theoretical simulation studies with different combinations of morphological and physiological remodeling, to assess and quantify their individual impact on overall RV load and function. Our findings revealed a substantial 4-fold increase in RV stiffness and a transient 2-fold rise in contractility, which returned to baseline by week 12. These changes in RV material properties in addition to the 2-fold increase in wall thickness significantly mitigated the increase in wall stress and strain caused by the progressive increase in RV afterload. Despite the PAH-induced cases showing increased wall stress and strain at end-diastole and end-systole compared to the control, our simulations suggest that without the observed remodeling mechanisms, the increase in stress and strain would have been much more pronounced. Our model analysis also indicated that while changes in the RV's material properties-particularly increased RV stiffness - have a notable effect on its mechanics, the primary compensatory factor limiting the stress and strain increase in the early stages of PAH was the significant increase in wall thickness. These findings underscore the importance of RV remodeling in managing the mechanical burden on the right ventricle due to pressure overload.

Establishment of a Biaxial Testing System for Characterization of Right Ventricle Viscoelasticity Under Physiological Loadings

PURPOSE

Prior studies have indicated an impact of cardiac muscle viscoelasticity on systolic and diastolic functions. However, the studies of ventricular free wall viscoelasticity, particularly for that of right ventricles (RV), are limited. Moreover, investigations on ventricular passive viscoelasticity have been restricted to large animals and there is a lack of data on rodent species. To fill this knowledge gap, this study aims to develop a biaxial tester that induces high-speed physiological deformations to characterize the passive viscoelasticity of rat RVs.

METHODS

The biaxial testing system was fabricated so that planar deformation of rat ventricle tissues at physiological strain rates was possible. The testing system was validated using isotropic polydimethylsiloxane (PDMS) sheets. Next, viscoelastic measurements were performed in healthy rat RV free walls by equibiaxial cyclic sinusoidal loadings and stress relaxation.

RESULTS

The biaxial tester's consistency, accuracy, and stability was confirmed from the PDMS samples measurements. Moreover, significant viscoelastic alterations of the RV were found between sub-physiological (0.1 Hz) and physiological frequencies (1-8 Hz). From hysteresis loop analysis, we found as the frequency increased, the elasticity and viscosity were increased in both directions. Interestingly, the ratio of storage energy to dissipated energy (Wd/Ws) remained constant at 0.1-5 Hz. We did not observe marked differences in healthy RV viscoelasticity between longitudinal and circumferential directions.

CONCLUSION

This work provides a new experimental tool to quantify the passive, biaxial viscoelasticity of ventricle free walls in both small and large animals. The dynamic mechanical tests showed frequency-dependent elastic and viscous behaviors of healthy rat RVs. But the ratio of dissipated energy to stored energy was maintained between frequencies. These findings offer novel baseline information on the passive viscoelasticity of healthy RVs in adult rats.

Using LDDMM and a kinematic cardiac growth model to quantify growth and remodelling in rat hearts under PAH

Pulmonary arterial hypertension (PAH) is a rapidly progressive and fatal disease, with right ventricular failure being the primary cause of death in patients with PAH. This study aims to determine the mechanical stimuli that may initiate heart growth and remodelling (G&R). To achieve this, two bi-ventricular models were constructed: one for a control rat heart and another for a rat heart with PAH. The growth of the diseased heart was estimated by warping it to the control heart using an improved large deformation diffeomorphic metric mapping (LDDMM) framework. Correlation analysis was then performed between mechanical cues (stress and strain) and growth tensors, which revealed that principal strains may serve as a triggering stimulus for myocardial growth and remodelling under PAH. The growth tensors, estimated from in vivo images, could explain 84.3% of the observed geometrical changes in the diseased heart with PAH by using a kinematic cardiac growth model. Our approach has the potential to quantify G&R using sparse in vivo images and to provide insights into the underlying mechanism of triggering right heart failure from a biomechanical perspective.

An image-driven micromechanical approach to characterize multiscale remodeling in infarcted myocardium

Myocardial infarction (MI) is accompanied by the formation of a fibrotic scar in the left ventricle (LV) and initiates significant alterations in the architecture and constituents of the LV free wall (LVFW). Previous studies have shown that LV adaptation is highly individual, indicating that the identification of remodeling mechanisms post-MI demands a fully subject-specific approach that can integrate a host of structural alterations at the fiber-level to changes in bulk biomechanical adaptation at the tissue-level. We present an image-driven micromechanical approach to characterize remodeling, assimilating new biaxial mechanical data, histological studies, and digital image correlation data within an in-silico framework to elucidate the fiber-level remodeling mechanisms that drive tissue-level adaptation for each subject. We found that a progressively diffused collagen fiber structure combined with similarly disorganized myofiber architecture in the healthy region leads to the loss of LVFW anisotropy post-MI, offering an important tissue-level hallmark for LV maladaptation. In contrast, our results suggest that reductions in collagen undulation are an adaptive mechanism competing against LVFW thinning. Additionally, we show that the inclusion of subject-specific geometry when modeling myocardial tissue is essential for accurate prediction of tissue kinematics. Our approach serves as an essential step toward identifying fiber-level remodeling indices that govern the transition of MI to systolic heart failure. These indices complement the traditional, organ-level measures of LV anatomy and function that often fall short of early prognostication of heart failure in MI. In addition, our approach offers an integrated methodology to advance the design of personalized interventions, such as hydrogel injection, to reinforce and suppress native adaptive and maladaptive mechanisms, respectively, to prevent the transition of MI to heart failure. STATEMENT OF SIGNIFICANCE: Biomechanical and architectural adaptation of the LVFW remains a central, yet overlooked, remodeling process post-MI. Our study indicates the biomechanical adaptation of the LVFW post-MI is highly individual and driven by altered fiber network architecture and collective changes in collagen fiber content, undulation, and stiffness. Our findings demonstrate the possibility of using cardiac strains to infer such fiber-level remodeling events through in-silico modeling, paving the way for in-vivo characterization of multiscale biomechanical indices in humans. Such indices will complement the traditional, organ-level measures of LV anatomy and function that often fall short of early prognostication of heart failure in MI.

Three-Dimensional, Right Ventricular Surface Strain Computation From Three-Dimensional Echocardiographic Images From Patients With Pediatric Pulmonary Hypertension

Right Ventricular (RV) dysfunction is routinely assessed with echocardiographic-derived global longitudinal strain (GLS). GLS is measured from a two-dimensional echo image and is increasingly accepted as a means for assessing RV function. However, any two-dimensional (2D) analysis cannot visualize the asymmetrical deformation of the RV nor visualize strain over the entire RV surface. We believe three-dimensional surface (3DS) strain, obtained from 3D echo will better evaluate myocardial mechanics. Components of 3DS strain (longitudinal, LS; circumferential, CS; longitudinal-circumferential shear, ɣCL; principal strains PSMax and PSMin; max shear, ɣMax; and principal angle θMax) were computed from RV surface meshes obtained with 3D echo from 50 children with associated pulmonary arterial hypertension (PAH), 43 children with idiopathic PAH, and 50 healthy children by computing strains from a discretized displacement field. All 3DS freewall (FW) normal strain (LS, CS, PSMax, and PSMin) showed significant decline at end-systole in PH groups (p < 0.0001 for all), as did FW-ɣMax (p = 0.0012). FW-θMax also changed in disease (p < 0.0001). Limits of agreement analysis suggest that 3DS LS, PSMax, and PSMin are related to GLS. 3DS strains showed significant heterogeneity over the 3D surface of the RV. Components of 3DS strain agree with existing clinical strain measures, well classify normal -versus- PAH subjects, and suggest that strains change direction on the myocardial surface due to disease. This last finding is similar to that of myocardial fiber realignment in disease, but further work is needed to establish true associations.

Untangling the mechanisms of pulmonary arterial hypertension-induced right ventricular stiffening in a large animal model

Pulmonary hypertension (PHT) is a devastating disease with low survival rates. In PHT, chronic pressure overload leads to right ventricle (RV) stiffening; thus, impeding diastolic filling. Multiple mechanisms may contribute to RV stiffening, including wall thickening, microstructural disorganization, and myocardial stiffening. The relative importance of each mechanism is unclear. Our objective is to use a large animal model to untangle these mechanisms. Thus, we induced pulmonary arterial hypertension (PAH) in sheep via pulmonary artery banding. After eight weeks, the hearts underwent anatomic and diffusion tensor MRI to characterize wall thickening and microstructural disorganization. Additionally, myocardial samples underwent histological and gene expression analyses to quantify compositional changes and mechanical testing to quantify myocardial stiffening. Finally, we used finite element modeling to disentangle the relative importance of each stiffening mechanism. We found that the RVs of PAH animals thickened most at the base and the free wall and that PAH induced excessive collagen synthesis, increased cardiomyocyte cross-sectional area, and led to microstructural disorganization, consistent with increased expression of fibrotic genes. We also found that the myocardium itself stiffened significantly. Importantly, myocardial stiffening correlated significantly with collagen synthesis. Finally, our computational models predicted that myocardial stiffness contributes to RV stiffening significantly more than other mechanisms. Thus, myocardial stiffening may be the most important predictor for PAH progression. Given the correlation between myocardial stiffness and collagen synthesis, collagen-sensitive imaging modalities may be useful for estimating myocardial stiffness and predicting PAH outcomes. STATEMENT OF SIGNIFICANCE: Ventricular stiffening is a significant contributor to pulmonary hypertension-induced right heart failure. However, the mechanisms that lead to ventricular stiffening are not fully understood. The novelty of our work lies in answering this question through the use of a large animal model in combination with spatially- and directionally sensitive experimental techniques. We find that myocardial stiffness is the primary mechanism that leads to ventricular stiffening. Clinically, this knowledge may be used to improve diagnostic, prognostic, and therapeutic strategies for patients with pulmonary hypertension.

Characterizing the Spatiotemporal Transcriptomic Response of the Right Ventricle to Acute Pressure Overload

This study analyzed microarray data of right ventricular (RV) tissue from rats exposed to pulmonary embolism to understand the initial dynamic transcriptional response to mechanical stress and compare it with experimental pulmonary hypertension (PH) models. The dataset included samples harvested from 55 rats at 11 different time points or RV locations. We performed principal component analysis (PCA) to explore clusters based on spatiotemporal gene expression. Relevant pathways were identified from fast gene set enrichment analysis using PCA coefficients. The RV transcriptomic signature was measured over several time points, ranging from hours to weeks after an acute increase in mechanical stress, and was found to be highly dependent on the severity of the initial insult. Pathways enriched in the RV outflow tracts of rats at 6 weeks after severe PE share many commonalities with experimental PH models, but the transcriptomic signature at the RV apex resembles control tissue. The severity of the initial pressure overload determines the trajectory of the transcriptomic response independent of the final afterload, but this depends on the location where the tissue is biopsied. Chronic RV pressure overload due to PH appears to progress toward similar transcriptomic endpoints.

A multi-scale computational model for the passive mechanical behavior of right ventricular myocardium

We have previously demonstrated the importance of myofiber-collagen mechanical interactions in modeling the passive mechanical behavior of right ventricle free wall (RVFW) myocardium. To gain deeper insights into these coupling mechanisms, we developed a high-fidelity, micro-anatomically realistic 3D finite element model of right ventricle free wall (RVFW) myocardium by combining high-resolution imaging and supercomputer-based simulations. We first developed a representative tissue element (RTE) model at the sub-tissue scale by specializing the hyperelastic anisotropic structurally-based constitutive relations for myofibers and ECM collagen, and equi-biaxial and non-equibiaxial loading conditions were simulated using the open-source software FEniCS to compute the effective stress-strain response of the RTE. To estimate the model parameters of the RTE model, we first fitted a 'top-down' biaxial stress-strain behavior with our previous structurally based (tissue-scale) model, informed by the measured myofiber and collagen fiber composition and orientation distributions. Next, we employed a multi-scale approach to determine the tissue-level (5 x 5 x 0.7 mm specimen size) RVFW biaxial behavior via 'bottom-up' homogenization of the fitted RTE model, recapitulating the histologically measured myofiber and collagen orientation to the biaxial mechanical data. Our homogenization approach successfully reproduced the tissue-level mechanical behavior of our previous studies in all biaxial deformation modes, suggesting that the 3D micro-anatomical arrangement of myofibers and ECM collagen is indeed a primary mechanism driving myofiber-collagen interactions.

On the possibility of estimating myocardial fiber architecture from cardiac strains

The myocardium is composed of a complex network of contractile myofibers that are organized in such a way as to produce efficient contraction and relaxation of the heart. The myofiber architecture in the myocardium is a key determinant of cardiac motion and the global or organ-level function of the heart. Reports of architectural remodeling in cardiac diseases, such as pulmonary hypertension and myocardial infarction, potentially contributing to cardiac dysfunction call for the inclusion of an architectural marker for an improved assessment of cardiac function. However, the in-vivo quantification of three-dimensional myo-architecture has proven challenging. In this work, we examine the sensitivity of cardiac strains to varying myofiber orientation using a multiscale finite-element model of the LV. Additionally, we present an inverse modeling approach to predict the myocardium fiber structure from cardiac strains. Our results indicate a strong correlation between fiber orientation and LV kinematics, corroborating that the fiber structure is a principal determinant of LV contractile behavior. Our inverse model was capable of accurately predicting the myocardial fiber range and regional fiber angles from strain measures. A concrete understanding of the link between LV myofiber structure and motion, and the development of non-invasive and feasible means of characterizing the myocardium architecture is expected to lead to advanced LV functional metrics and improved prognostic assessment of structural heart disease.

Alterations of biaxial viscoelastic properties of the right ventricle in pulmonary hypertension development in rest and acute stress conditions

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.

Recent Advances in Tissue-Engineered Cardiac Scaffolds-The Progress and Gap in Mimicking Native Myocardium Mechanical Behaviors

Heart failure is the leading cause of death in the US and worldwide. Despite modern therapy, challenges remain to rescue the damaged organ that contains cells with a very low proliferation rate after birth. Developments in tissue engineering and regeneration offer new tools to investigate the pathology of cardiac diseases and develop therapeutic strategies for heart failure patients. Tissue -engineered cardiac scaffolds should be designed to provide structural, biochemical, mechanical, and/or electrical properties similar to native myocardium tissues. This review primarily focuses on the mechanical behaviors of cardiac scaffolds and their significance in cardiac research. Specifically, we summarize the recent development of synthetic (including hydrogel) scaffolds that have achieved various types of mechanical behavior-nonlinear elasticity, anisotropy, and viscoelasticity-all of which are characteristic of the myocardium and heart valves. For each type of mechanical behavior, we review the current fabrication methods to enable the biomimetic mechanical behavior, the advantages and limitations of the existing scaffolds, and how the mechanical environment affects biological responses and/or treatment outcomes for cardiac diseases. Lastly, we discuss the remaining challenges in this field and suggestions for future directions to improve our understanding of mechanical control over cardiac function and inspire better regenerative therapies for myocardial restoration.

Untangling the mechanisms of pulmonary hypertension-induced right ventricular stiffening in a large animal model

Background

Pulmonary arterial hypertension (PHT) is a devastating disease with low survival rates. In PHT, chronic pressure overload leads to right ventricle (RV) remodeling and stiffening; thus, impeding diastolic filling and ventricular function. Multiple mechanisms contribute to RV stiffening, including wall thickening, microstructural disorganization, and myocardial stiffening. The relative importance of each mechanism is unclear. Our objective is to use a large animal model as well as imaging, experimental, and computational approaches to untangle these mechanisms.

Methods

We induced PHT in eight sheep via pulmonary artery banding. After eight weeks, the hearts underwent anatomic and diffusion tensor MRI to characterize wall thickening and microstructural disorganization. Additionally, myocardial samples underwent histological and gene expression analyses to quantify compositional changes and mechanical testing to quantify myocardial stiffening. All findings were compared to 12 control animals. Finally, we used computational modeling to disentangle the relative importance of each stiffening mechanism.

Results

First, we found that the RVs of PHT animals thickened most at the base and the free wall. Additionally, we found that PHT induced excessive collagen synthesis and microstructural disorganization, consistent with increased expression of fibrotic genes. We also found that the myocardium itself stiffened significantly. Importantly, myocardial stiffening correlated significantly with excess collagen synthesis. Finally, our model of normalized RV pressure-volume relationships predicted that myocardial stiffness contributes to RV stiffening significantly more than other mechanisms.

Conclusions

In summary, we found that PHT induces wall thickening, microstructural disorganization, and myocardial stiffening. These remodeling mechanisms were both spatially and directionally dependent. Using modeling, we show that myocardial stiffness is the primary contributor to RV stiffening. Thus, myocardial stiffening may be an important predictor for PHT progression. Given the significant correlation between myocardial stiffness and collagen synthesis, collagen-sensitive imaging modalities may be useful for non-invasively estimating myocardial stiffness and predicting PHT outcomes.

Right Ventricular Architectural Remodeling and Functional Adaptation in Pulmonary Hypertension

BACKGROUND

Global indices of right ventricle (RV) function provide limited insights into mechanisms underlying RV remodeling in pulmonary hypertension (PH). While RV myocardial architectural remodeling has been observed in PH, its effect on RV adaptation is poorly understood.

METHODS

Hemodynamic assessments were performed in 2 rodent models of PH. RV free wall myoarchitecture was quantified using generalized Q-space imaging and tractography analyses. Computational models were developed to predict RV wall strains. Data from animal studies were analyzed to determine the correlations between hemodynamic measurements, RV strains, and structural measures.

RESULTS

In contrast to the PH rats with severe RV maladaptation, PH rats with mild RV maladaptation showed a decrease in helical range of fiber orientation in the RV free wall (139º versus 97º; P=0.029), preserved global circumferential strain, and exhibited less reduction in right ventricular-pulmonary arterial coupling (0.029 versus 0.017 mm/mm Hg; P=0.037). Helical range correlated positively with coupling (P=0.036) and stroke volume index (P<0.01). Coupling correlated with global circumferential strain (P<0.01) and global radial strain (P<0.01) but not global longitudinal strain.

CONCLUSIONS

Data analysis suggests that adaptive RV architectural remodeling could improve RV function in PH. Our findings suggest the need to assess RV architecture within routine screenings of PH patients to improve our understanding of its prognostic and therapeutic significance in PH.

Computational models of ventricular mechanics and adaptation in response to right-ventricular pressure overload

Pulmonary arterial hypertension (PAH) is associated with substantial remodeling of the right ventricle (RV), which may at first be compensatory but at a later stage becomes detrimental to RV function and patient survival. Unlike the left ventricle (LV), the RV remains understudied, and with its thin-walled crescent shape, it is often modeled simply as an appendage of the LV. Furthermore, PAH diagnosis is challenging because it often leaves the LV and systemic circulation largely unaffected. Several treatment strategies such as atrial septostomy, right ventricular assist devices (RVADs) or RV resynchronization therapy have been shown to improve RV function and the quality of life in patients with PAH. However, evidence of their long-term efficacy is limited and lung transplantation is still the most effective and curative treatment option. As such, the clinical need for improved diagnosis and treatment of PAH drives a strong need for increased understanding of drivers and mechanisms of RV growth and remodeling (G&R), and more generally for targeted research into RV mechanics pathology. Computational models stand out as a valuable supplement to experimental research, offering detailed analysis of the drivers and consequences of G&R, as well as a virtual test bench for exploring and refining hypotheses of growth mechanisms. In this review we summarize the current efforts towards understanding RV G&R processes using computational approaches such as reduced-order models, three dimensional (3D) finite element (FE) models, and G&R models. In addition to an overview of the relevant literature of RV computational models, we discuss how the models have contributed to increased scientific understanding and to potential clinical treatment of PAH patients.

Mechanical Interaction of the Pericardium and Cardiac Function in the Normal and Hypertensive Rat Heart

The pericardium is a thin connective tissue membrane that surrounds the heart and is an integral regulatory component of cardiopulmonary performance. Pathological growth and remodeling of the right ventricle (RV) stemming from structural heart diseases are thought to include a significant role of the pericardium, but its exact role remains unclear. The objective of this study was to investigate potential biomechanical adaptations of the pericardium in response to pulmonary hypertension and their effects on heart behavior. Integrated computational-experimental modeling of the heart offers a robust platform to achieve this objective. We built upon our recently developed high-fidelity finite-element models of healthy and hypertensive rodent hearts via addition of the pericardial sac. In-silico experiments were performed to investigate changes in pericardium reserve elasticity and their effects on cardiac function in hypertensive hearts. Our results suggest that contractile forces would need to increase in the RV and decrease in the left ventricle (LV) in the hypertensive heart to compensate for reductions in pericardium reserve elasticity. The discrepancies between chamber responses to pericardium addition result, in part, from differences in the impact of pericardium on the RV and LV preload. We further demonstrated the capability of our platform to predict the effect of pericardiectomy on heart function. Consistent with previous results, the effect of pericardiectomy on the chamber pressure-volume loop was the largest in the hypertensive RV. These insights are expected to motivate further computational investigations of the effect of pericardiectomy on cardiac function which remains an important factor in surgical planning of constrictive pericarditis and coronary artery bypass grafting.

Simulation of the 3D Hyperelastic Behavior of Ventricular Myocardium using a Finite-Element Based Neural-Network Approach

High-fidelity cardiac models using attribute-rich finite element based models have been developed to a very mature stage. However, such finite-element based approaches remain time consuming, which have limited their clinical use. There remains a need for alternative methods for novel cardiac simulation methods of capable of high fidelity simulations in clinically relevant time frames. Surrogate models are one approach, which traditionally use a data-driven approach for training, requiring the generation of a sufficiently large number of simulation results as the training dataset. Alternatively, a physics-informed neural network can be trained by minimizing the PDE residuals or energy potentials. However, this approach does not provide for a general method to easily using existing finite element models. To address these challenges, we developed a hybrid approach that seamlessly bridged a neural network surrogate model with a differentiable finite element domain representation (NNFE). Given its importance in cardiac simulations, we applied this approach to simulations of the hyperelastic mechanical behavior of ventricular myocardium from recent 3D kinematic constitutive model (J Mech Behav Biomed Mater, 2020 doi: 10.1016/j.jmbbm.2019.103508). We utilized cuboidal domain and conducted numerical studies of individual myocardium specimens discretized by a finite element mesh and assigned with experimentally obtained myofiber architectures. Both parameterized Dirichlet and Neumann boundary conditions were studied. We developed a second-order Newton optimization method, instead of using stochastic gradient descent method, to train the neural network efficiently. The resulting trained neural network surrogate model demonstrated excellent agreement with the corresponding 'ground truth' finite element solutions over the entire physiological deformation range. More importantly, the NNFE approach provided a significantly decreased computational time for a range of finite element mesh sizes for online predictions. For example, as the finite element mesh sized increased from 2744 to 175615 elements the NNFE computational time increased from 0.1108 s to 0.1393 s, while the 'ground truth' FE model increased from 4.541 s to 719.9 s. These results suggests that NNFE run times can be significantly reduced compared with the traditional large-deformation based finite element solution methods. The trade off is to train the NNFE off-line within a range of anticipated physiological responses. However, training time would only have to be performed once before any number of application uses. Moreover, since the NNFE is an analytical function its computational performance will be amplified when the corresponding problem becomes more complex.

Multiscale Contrasts Between the Right and Left Ventricle Biomechanics in Healthy Adult Sheep and Translational Implications

Cardiac biomechanics play a significant role in the progression of structural heart diseases (SHDs). SHDs alter baseline myocardial biomechanics leading to single or bi-ventricular dysfunction. But therapies for left ventricle (LV) failure patients do not always work well for right ventricle (RV) failure patients. This is partly because the basic knowledge of baseline contrasts between the RV and LV biomechanics remains elusive with limited discrepant findings. The aim of the study was to investigate the multiscale contrasts between LV and RV biomechanics in large animal species. We hypothesize that the adult healthy LV and RV have distinct passive anisotropic biomechanical properties. Ex vivo biaxial tests were performed in fresh sheep hearts. Histology and immunohistochemistry were performed to measure tissue collagen. The experimental data were then fitted to a Fung type model and a structurally informed model, separately. We found that the LV was stiffer in the longitudinal (outflow tract) than circumferential direction, whereas the RV showed the opposite anisotropic behavior. The anisotropic parameter K from the Fung type model accurately captured contrasting anisotropic behaviors in the LV and RV. When comparing the elasticity in the same direction, the LV was stiffer than the RV longitudinally and the RV was stiffer than the LV circumferentially, suggesting different filling patterns of these ventricles during diastole. Results from the structurally informed model suggest potentially stiffer collagen fibers in the LV than RV, demanding further investigation. Finally, type III collagen content was correlated with the low-strain elastic moduli in both ventricles. In summary, our findings provide fundamental biomechanical differences between the chambers. These results provide valuable insights for guiding cardiac tissue engineering and regenerative studies to implement chamber-specific matrix mechanics, which is particularly critical for identifying biomechanical mechanisms of diseases or mechanical regulation of therapeutic responses. In addition, our results serve as a benchmark for image-based inverse modeling technologies to non-invasively estimate myocardial properties in the RV and LV.

The Effects of Healthy Aging on Right Ventricular Structure and Biomechanical Properties: A Pilot Study

Healthy aging has been associated with alterations in pulmonary vascular and right ventricular (RV) hemodynamics, potentially leading to RV remodeling. Despite the current evidence suggesting an association between aging and alterations in RV function and higher prevalence of pulmonary hypertension in the elderly, limited data exist on age-related differences in RV structure and biomechanics. In this work, we report our preliminary findings on the effects of healthy aging on RV structure, function, and biomechanical properties. Hemodynamic measurements, biaxial mechanical testing, constitutive modeling, and quantitative transmural histological analysis were employed to study two groups of male Sprague-Dawley rats: control (11 weeks) and aging (80 weeks). Aging was associated with increases in RV peak pressures (+17%, p = 0.017), RV contractility (+52%, p = 0.004), and RV wall thickness (+38%, p = 0.001). Longitudinal realignment of RV collagen (16.4°, p = 0.013) and myofibers (14.6°, p = 0.017) were observed with aging, accompanied by transmural cardiomyocyte loss and fibrosis. Aging led to increased RV myofiber stiffness (+141%, p = 0.003), in addition to a bimodal alteration in the biaxial biomechanical properties of the RV free wall, resulting in increased tissue-level stiffness in the low-strain region, while progressing into decreased stiffness at higher strains. Our results demonstrate that healthy aging may modulate RV remodeling via increased peak pressures, cardiomyocyte loss, fibrosis, fiber reorientation, and altered mechanical properties in male Sprague-Dawley rats. Similarities were observed between aging-induced remodeling patterns and those of RV remodeling in pressure overload. These findings may help our understanding of age-related changes in the cardiovascular fitness and response to disease.

Current Understanding of the Right Ventricle Structure and Function in Pulmonary Arterial Hypertension

Pulmonary arterial hypertension (PAH) is a disease resulting in increased right ventricular (RV) afterload and RV remodeling. PAH results in altered RV structure and function at different scales from organ-level hemodynamics to tissue-level biomechanical properties, fiber-level architecture, and cardiomyocyte-level contractility. Biomechanical analysis of RV pathophysiology has drawn significant attention over the past years and recent work has found a close link between RV biomechanics and physiological function. Building upon previously developed techniques, biomechanical studies have employed multi-scale analysis frameworks to investigate the underlying mechanisms of RV remodeling in PAH and effects of potential therapeutic interventions on these mechanisms. In this review, we discuss the current understanding of RV structure and function in PAH, highlighting the findings from recent studies on the biomechanics of RV remodeling at organ, tissue, fiber, and cellular levels. Recent progress in understanding the underlying mechanisms of RV remodeling in PAH, and effects of potential therapeutics, will be highlighted from a biomechanical perspective. The clinical relevance of RV biomechanics in PAH will be discussed, followed by addressing the current knowledge gaps and providing suggested directions for future research.

Right ventricular myocardial mechanics: Multi-modal deformation, microstructure, modeling, and comparison to the left ventricle

The right ventricular myocardium, much like the rest of the right side of the heart, has been consistently understudied. Presently, little is known about its mechanics, its microstructure, and its constitutive behavior. In this work, we set out to provide the first data on the mechanics of the mature right ventricular myocardium in both simple shear and uniaxial loading and to compare these data to the mechanics of the left ventricular myocardium. To this end, we tested ovine tissue samples of the right and left ventricle under a comprehensive mechanical testing protocol that consisted of six simple shear modes and three tension/compression modes. After mechanical testing, we conducted a histology-based microstructural analysis on each right ventricular sample that yielded high resolution fiber distribution maps across the entire samples. Equipped with this detailed mechanical and histological data, we employed an inverse finite element framework to determine the optimal form and parameters for microstructure-based constitutive models. The results of our study show that right ventricular myocardium is less stiff then the left ventricular myocardium in the fiber direction, but similarly exhibits non-linear, anisotropic, and tension/compression asymmetric behavior with direction-dependent Poynting effect. In addition, we found that right ventricular myocardial fibers change angles transmurally and are dispersed within the sheet plane and normal to it. Through our inverse finite element analysis, we found that the Holzapfel model successfully fits these data, even when selectively informed by rudimentary microstructural information. And, we found that the inclusion of higher-fidelity microstructural data improved the Holzapfel model's predictive ability. Looking forward, this investigation is a critical step towards understanding the fundamental mechanical behavior of right ventricular myocardium and lays the groundwork for future whole-organ mechanical simulations.

An exploratory assessment of stretch-induced transmural myocardial fiber kinematics in right ventricular pressure overload

Right ventricular (RV) remodeling and longitudinal fiber reorientation in the setting of pulmonary hypertension (PH) affects ventricular structure and function, eventually leading to RV failure. Characterizing the kinematics of myocardial fibers helps better understanding the underlying mechanisms of fiber realignment in PH. In the current work, high-frequency ultrasound imaging and structurally-informed finite element (FE) models were employed for an exploratory evaluation of the stretch-induced kinematics of RV fibers. Image-based experimental evaluation of fiber kinematics in porcine myocardium revealed the capability of affine assumptions to effectively approximate myofiber realignment in the RV free wall. The developed imaging framework provides a noninvasive modality to quantify transmural RV myofiber kinematics in large animal models. FE modeling results demonstrated that chronic pressure overload, but not solely an acute rise in pressures, results in kinematic shift of RV fibers towards the longitudinal direction. Additionally, FE simulations suggest a potential protective role for concentric hypertrophy (increased wall thickness) against fiber reorientation, while eccentric hypertrophy (RV dilation) resulted in longitudinal fiber realignment. Our study improves the current understanding of the role of different remodeling events involved in transmural myofiber reorientation in PH. Future experimentations are warranted to test the model-generated hypotheses.

Angiotensin Receptor-Neprilysin Inhibition Attenuates Right Ventricular Remodeling in Pulmonary Hypertension

Background

Pulmonary hypertension (PH) results in increased right ventricular (RV) afterload and ventricular remodeling. Sacubitril/valsartan (sac/val) is a dual acting drug, composed of the neprilysin inhibitor sacubitril and the angiotensin receptor blocker valsartan, that has shown promising outcomes in reducing the risk of death and hospitalization for chronic systolic left ventricular heart failure. In this study, we aimed to examine if angiotensin receptor‐neprilysin inhibition using sac/val attenuates RV remodeling in PH.

Methods and Results

RV pressure overload was induced in Sprague–Dawley rats via banding the main pulmonary artery. Three different cohorts of controls, placebo‐treated PH, and sac/val‐treated PH were studied in a 21‐day treatment window. Terminal invasive hemodynamic measurements, quantitative histological analysis, biaxial mechanical testing, and constitutive modeling were employed to conduct a multiscale analysis on the effects of sac/val on RV remodeling in PH. Sac/val treatment decreased RV maximum pressures (29% improvement, P=0.002), improved RV contractile (30%, P=0.012) and relaxation (29%, P=0.043) functions, reduced RV afterload (35% improvement, P=0.016), and prevented RV‐pulmonary artery uncoupling. Furthermore, sac/val attenuated RV hypertrophy (16% improvement, P=0.006) and prevented transmural reorientation of RV collagen and myofibers (P=0.011). The combined natriuresis and vasodilation resulting from sac/val led to improved RV biomechanical properties and prevented increased myofiber stiffness in PH (61% improvement, P=0.032).

Conclusions

Sac/val may prevent maladaptive RV remodeling in a pressure overload model via amelioration of RV pressure rise, hypertrophy, collagen, and myofiber reorientation as well as tissue stiffening both at the tissue and myofiber level.

Pathophysiology of Acute and Chronic Right Heart Failure

Right-sided heart failure (RHF) occurs from impaired contractility of the right ventricle caused by pressure, volume overload, or intrinsic myocardial contractile dysfunction. The development of subclinical right ventricle (RV) dysfunction or overt RHF is a negative prognostic indicator. Recent attention has focused on RV-specific inflammatory growth factors and mediators of myocardial fibrosis to elucidate the mechanisms leading to RHF and potentially guide the development of novel therapeutics. This article focuses on the distinct changes in RV structure, mechanics, and function, as well as molecular and inflammatory mediators involved in the pathophysiology of acute and chronic RHF.

Insights into the passive mechanical behavior of left ventricular myocardium using a robust constitutive model based on full 3D kinematics

Myocardium possesses a hierarchical structure that results in complex three-dimensional (3D) mechanical behavior, forming a critical component of ventricular function in health and disease. A wide range of constitutive model forms have been proposed for myocardium since the first planar biaxial studies were performed by Demer and Yin (J. Physiol. 339 (1), 1983). While there have been extensive studies since, none have been based on full 3D kinematic data, nor have they utilized optimal experimental design to estimate constitutive parameters, which may limit their predictive capability. Herein we have applied our novel 3D numerical-experimental methodology (Avazmohammadi et al., Biomechanics Model. Mechanobiol. 2018) to explore the applicability of an orthotropic constitutive model for passive ventricular myocardium (Holzapfel and Ogden, Philos. Trans. R. Soc. Lond.: Math. Phys. Eng. Sci. 367, 2009) by integrating 3D optimal loading paths, spatially varying material structure, and inverse modeling techniques. Our findings indicated that the initial model form was not successful in reproducing all optimal loading paths, due to previously unreported coupling behaviors via shearing of myofibers and extracellular collagen fibers in the myocardium. This observation necessitated extension of the constitutive model by adding two additional terms based on the I8(C) pseudo-invariant in the fiber-normal and sheet-normal directions. The modified model accurately reproduced all optimal loading paths and exhibited improved predictive capabilities. These unique results suggest that more complete constitutive models are required to fully capture the full 3D biomechanical response of left ventricular myocardium. The present approach is thus crucial for improved understanding and performance in cardiac modeling in healthy, diseased, and treatment scenarios.

Oxidative Stress and Its Implications in the Right Ventricular Remodeling Secondary to Pulmonary Hypertension

Pulmonary hypertension (PH) is a pulmonary vascular disease characterized by increased pulmonary artery pressures. Long standing pulmonary arterial pressure overload leads to right ventricular (RV) hypertrophy, RV failure, and death. RV failure is a major determinant of survival in PH. Oxidative stress has been associated with the development of RV failure secondary to PH. Here we summarize the structural and functional changes in the RV in response to sustained pulmonary arterial pressure overload. Furthermore, we review the pre-clinical and clinical studies highlighting the association of oxidative stress with pulmonary vasculature and RV remodeling in chronic PH. Targeting oxidative stress promises to be an effective therapeutic strategy for the treatment of RV failure.

Relative Contributions of Matrix and Myocytes to Biaxial Mechanics of the Right Ventricle in Pulmonary Arterial Hypertension

Pulmonary arterial hypertension (PAH) commonly leads to right ventricular (RV) hypertrophy and fibrosis that affect the mechanical properties of the RV myocardium (MYO). To investigate the effects of PAH on the mechanics of the RV MYO and extracellular matrix (ECM), we compared RV wall samples, isolated from rats in which PAH was induced using the SuHx protocol, with samples from control animals before and after the tissues were decellularized. Planar biaxial mechanical testing, a technique first adapted to living soft biological tissues by Fung, was performed on intact and decellularized samples. Fung's anisotropic exponential strain energy function fitted the full range of biaxial test results with high fidelity in control and PAH samples both before and after they were decellularized. Mean RV myocardial apex-to-outflow tract and circumferential stresses during equibiaxial strain were significantly greater in PAH than control samples. Mean RV ECM circumferential but not apex-to-outflow tract stresses during equibiaxial strain were significantly greater in the PAH than control group. The ratio of ECM to myocardial stresses at matched strains did not change significantly between groups. Circumferential stresses were significantly higher than apex-to-outflow tract stresses for all groups. These findings confirm the predictions of a mathematical model based on changes in RV hemodynamics and morphology in rat PAH, and may provide a foundation for a new constitutive analysis of the contributions of ECM remodeling to changes in RV filling properties during PAH.

On the AIC-based model reduction for the general Holzapfel-Ogden myocardial constitutive law

Constitutive laws that describe the mechanical responses of cardiac tissue under loading hold the key to accurately model the biomechanical behaviour of the heart. There have been ample choices of phenomenological constitutive laws derived from experiments, some of which are quite sophisticated and include effects of microscopic fibre structures of the myocardium. A typical example is the strain-invariant-based Holzapfel-Ogden 2009 model that is excellently fitted to simple shear tests. It has been widely used and regarded as the state-of-the-art constitutive law for myocardium. However, there has been no analysis to show if it has both adequate descriptive and predictive capabilities for other tissue tests of myocardium. Indeed, such an analysis is important for any constitutive laws for clinically useful computational simulations. In this work, we perform such an analysis using combinations of tissue tests, uniaxial tension, biaxial tension and simple shear from three different sets of myocardial tissue studies. Starting from the general 14-parameter myocardial constitutive law developed by Holzapfel and Ogden, denoted as the general HO model, we show that this model has good descriptive and predictive capabilities for all the experimental tests. However, to reliably determine all 14 parameters of the model from experiments remains a great challenge. Our aim is to reduce the constitutive law using Akaike information criterion, to maintain its mechanical integrity whilst achieving minimal computational cost. A competent constitutive law should have descriptive and predictive capabilities for different tissue tests. By competent, we mean the model has least terms but is still able to describe and predict experimental data. We also investigate the optimal combinations of tissue tests for a given constitutive model. For example, our results show that using one of the reduced HO models, one may need just one shear response (along normal-fibre direction) and one biaxial stretch (ratio of 1 mean fibre : 1 cross-fibre) to satisfactorily describe Sommer et al. human myocardial mechanical properties. Our study suggests that single-state tests (i.e. simple shear or stretching only) are insufficient to determine the myocardium responses. We also found it is important to consider transmural fibre rotations within each myocardial sample of tests during the fitting process. This is done by excluding un-stretched fibres using an "effective fibre ratio", which depends on the sample size, shape, local myofibre architecture and loading conditions. We conclude that a competent myocardium material model can be obtained from the general HO model using AIC analysis and a suitable combination of tissue tests.

Interactions between structural remodeling and volumetric growth in right ventricle in response to pulmonary arterial hypertension

Pulmonary arterial hypertension (PAH) exerts substantial pressure overload on the right ventricle (RV). The associated RV free wall (RVFW) adaptation could consist of myocardial hypertrophy, augmented intrinsic contractility, collagen fibrosis, and structural remodeling in an attempt to cope with pressure overload. If RVFW adaptation cannot maintain the RV stroke volume, RV dilation will prevail as an exit mechanism which usually decompensates the RV function leading to RV failure. Our knowledge of the factors determining the transition from the upper limit of RVFW adaptation to RV decompensation and the role of fiber remodeling events in this transition remains very limited. Computational heart models that connect the growth and remodeling (G\&R) events at the fiber and tissue levels with alterations in the organ-level function are essential to predict the temporal order and the compensatory level of the underlying mechanisms. In this work, building upon our recent rodent heart models (RHM) of PAH, we integrated mathematical models that describe time-evolution volumetric growth of the RV and structural remodeling of the RVFW. Results suggest that augmentation of the intrinsic contractility of myofibers accompanied by an increase in passive stiffness of RVFW is among the first remodeling events through which the RV strives to maintain the cardiac output. Interestingly, we found that the observed reorientation of the myofibers towards the longitudinal (apex-to-base) direction was a maladaptive mechanism that impaired the contractile pattern of RVFW and advanced along with RV dilation at later stages of PAH development.

A Contemporary Look at Biomechanical Models of Myocardium

Understanding and predicting the mechanical behavior of myocardium under healthy and pathophysiological conditions are vital to developing novel cardiac therapies and promoting personalized interventions. Within the past 30 years, various constitutive models have been proposed for the passive mechanical behavior of myocardium. These models cover a broad range of mathematical forms, microstructural observations, and specific test conditions to which they are fitted. We present a critical review of these models, covering both phenomenological and structural approaches, and their relations to the underlying structure and function of myocardium. We further explore the experimental and numerical techniques used to identify the model parameters. Next, we provide a brief overview of continuum-level electromechanical models of myocardium, with a focus on the methods used to integrate the active and passive components of myocardial behavior. We conclude by pointing to future directions in the areas of optimal form as well as new approaches for constitutive modeling of myocardium.

In-silico assessment of the effects of right ventricular assist device on pulmonary arterial hypertension using an image based biventricular modeling framework

Pulmonary arterial hypertension (PAH) is a heart disease that is characterized by an abnormally high pressure in the pulmonary artery (PA). While right ventricular assist device (RVAD) has been considered recently as a treatment option for the end-stage PAH patients, its effects on biventricular mechanics are, however, largely unknown. To address this issue, we developed an image-based modeling framework consisting of a biventricular finite element (FE) model that is coupled to a lumped model describing the pulmonary and systemic circulations in a closed-loop system. The biventricular geometry was reconstructed from the magnetic resonance images of two PAH patients showing different degree of RV remodeling and a normal subject. The framework was calibrated to match patient-specific measurements of the left ventricular (LV) and RV volume and pressure waveforms. An RVAD model was incorporated into the calibrated framework and simulations were performed with different pump speeds. Results showed that RVAD unloads the RV, improves cardiac output and increases septum curvature, which are more pronounced in the PAH patient with severe RV remodeling. These improvements, however, are also accompanied by an adverse increase in the PA pressure. These results suggest that the RVAD implantation may need to be optimized depending on disease progression.

A Computational Cardiac Model for the Adaptation to Pulmonary Arterial Hypertension in the Rat

Pulmonary arterial hypertension (PAH) imposes pressure overload on the right ventricle (RV), leading to RV enlargement via the growth of cardiac myocytes and remodeling of the collagen fiber architecture. The effects of these alterations on the functional behavior of the right ventricular free wall (RVFW) and organ-level cardiac function remain largely unexplored. Computational heart models in the rat (RHMs) of the normal and hypertensive states can be quite valuable in simulating the effects of PAH on cardiac function to gain insights into the pathophysiology of underlying myocardium remodeling. We thus developed high-fidelity biventricular finite element RHMs for the normal and post-PAH hypertensive states using extensive experimental data collected from rat hearts. We then applied the RHM to investigate the transmural nature of RVFW remodeling and its connection to wall stress elevation under PAH. We found a strong correlation between the longitudinally-dominated fiber-level adaptation of the RVFW and the transmural alterations of relevant wall stress components. We further conducted several numerical experiments to gain new insights on how the RV responds both normally and in the post-PAH state. We found that the effect of pressure overload alone on the increased contractility of the RV is comparable to the effects of changes in the RV geometry and stiffness. Furthermore, our RHMs provided fresh perspectives on long-standing questions of the functional role of the interventricular septum in RV function. Specifically, we demonstrated that an inaccurate identification of the mechanical adaptation of the septum can lead to a significant underestimation of RVFW contractility in the post-PAH state. These findings show how integrated experimental-computational models can facilitate a more comprehensive understanding of the cardiac remodeling events during PAH.