Multiscale Computational Analysis of Right Ventricular Mechanoenergetics

An in-silico analysis of experimental designs to study ventricular function: A focus on the right ventricle

In-vivo studies of pulmonary vascular disease and pulmonary hypertension (PH) have provided key insight into the progression of right ventricular (RV) dysfunction. Additional in-silico experiments using multiscale computational models have provided further details into biventricular mechanics and hemodynamic function in the presence of PH, yet few have assessed whether model parameters are practically identifiable prior to data collection. Moreover, none have used modeling to devise synergistic experimental designs. To address this knowledge gap, we conduct a practical identifiability analysis of a multiscale cardiovascular model across four simulated experimental designs. We determine a set of parameters using a combination of Morris screening and local sensitivity analysis, and test for practical identifiability using profile likelihood-based confidence intervals. We employ Markov chain Monte Carlo (MCMC) techniques to quantify parameter and model forecast uncertainty in the presence of noise corrupted data. Our results show that model calibration to only RV pressure suffers from practical identifiability issues and suffers from large forecast uncertainty in output space. In contrast, parameter and model forecast uncertainty is substantially reduced once additional left ventricular (LV) pressure and volume data is included. A comparison between single point systolic and diastolic LV data and continuous, time-dependent LV pressure-volume data reveals that at least some quantitative data from both ventricles should be included for future experimental studies.

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.

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.

Reversal of right ventricular failure by chronic α1A-subtype adrenergic agonist therapy

Right ventricular (RV) failure (RVF) is a serious disease with no effective treatment available. We recently reported a disease prevention study showing that chronic stimulation of α1A-adrenergic receptors (α1A-ARs), started at the time of RV injury, prevented the development of RVF. The present study used a clinically relevant disease reversal design to test if chronic α1A-AR stimulation, started after RVF was established, could reverse RVF. RVF was induced surgically by pulmonary artery constriction in mice. Two weeks after pulmonary artery constriction, in vivo RV fractional shortening as assessed by MRI was reduced by half relative to sham-operated controls (25 ± 2%, n = 27, vs. 52 ± 2%, n = 13, P < 10-11). Subsequent chronic treatment with the α1A-AR agonist A61603 for a further 2 wk resulted in a substantial recovery of RV fractional shortening (to 41 ± 2%, n = 17, P < 10-7 by a paired t-test) along with recovery of voluntary exercise capacity. Mechanistically, chronic A61603 treatment resulted in increased activation of the prosurvival kinase ERK, increased abundance of the antiapoptosis factor Bcl-2, and decreased myocyte necrosis evidenced by a decreased serum level of cardiac troponin. Moreover, A61603 treatment caused increased abundance of the antioxidant glutathione peroxidase-1, decreased level of reactive oxygen species, and decreased oxidative modification (carbonylation) of myofilament proteins. Consistent with these effects, A61603 treatment resulted in increased force development by cardiac myofilaments, which might have contributed to increased RV function. These findings suggest that the α1A-AR is a therapeutic target to reverse established RVF. NEW & NOTEWORTHY Currently, there are no effective therapies for right ventricular (RV) failure (RVF). This project evaluated a novel therapy for RVF. In a mouse model of RVF, chronic stimulation of α1A-adrenergic receptors with the agonist A61603 resulted in recovery of in vivo RV function, improved exercise capacity, reduced oxidative stress-related carbonylation of contractile proteins, and increased myofilament force generation. These results suggest that the α1A-adrenergic receptor is a therapeutic target to treat RVF.

Impaired Myofilament Contraction Drives Right Ventricular Failure Secondary to Pressure Overload: Model Simulations, Experimental Validation, and Treatment Predictions

Introduction: Pulmonary hypertension (PH) causes pressure overload leading to right ventricular failure (RVF). Myocardial structure and myocyte mechanics are altered in RVF but the direct impact of these cellular level factors on organ level function remain unclear. A computational model of the cardiovascular system that integrates cellular function into whole organ function has recently been developed. This model is a useful tool for investigating how changes in myocyte structure and mechanics contribute to organ function. We use this model to determine how measured changes in myocyte and myocardial mechanics contribute to RVF at the organ level and predict the impact of myocyte-targeted therapy. Methods: A multiscale computational framework was tuned to model PH due to bleomycin exposure in mice. Pressure overload was modeled by increasing the pulmonary vascular resistance (PVR) and decreasing pulmonary artery compliance (CPA). Myocardial fibrosis and the impairment of myocyte maximum force generation (Fmax) were simulated by increasing the collagen content (↑PVR + ↓CPA + fibrosis) and decreasing Fmax (↑PVR + ↓CPA + fibrosis + ↓Fmax). A61603 (A6), a selective α1A-subtype adrenergic receptor agonist, shown to improve Fmax was simulated to explore targeting myocyte generated Fmax in PH. Results: Increased afterload (RV systolic pressure and arterial elastance) in simulations matched experimental results for bleomycin exposure. Pressure overload alone (↑PVR + ↓CPA) caused decreased RV ejection fraction (EF) similar to experimental findings but preservation of cardiac output (CO). Myocardial fibrosis in the setting of pressure overload (↑PVR + ↓PAC + fibrosis) had minimal impact compared to pressure overload alone. Including impaired myocyte function (↑PVR + ↓PAC + fibrosis + ↓Fmax) reduced CO, similar to experiment, and impaired EF. Simulations predicted that A6 treatment preserves EF and CO despite maintained RV pressure overload. Conclusion: Multiscale computational modeling enabled prediction of the contribution of cellular level changes to whole organ function. Impaired Fmax is a key feature that directly contributes to RVF. Simulations further demonstrate the therapeutic benefit of targeting Fmax, which warrants additional study. Future work should incorporate growth and remodeling into the computational model to enable prediction of the multiscale drivers of the transition from dysfunction to failure.