Systems analysis of the mechanisms governing the cardiovascular response to changes in posture and in peripheral demand during exercise

The importance of incorporating ventricular-ventricular interaction (VVI) in the study of pulmonary hypertension

Ventricular ventricular interaction (VVI) affects blood volume and pressure in the right and left ventricles of the heart due to the location and balance of forces on the septal wall separating the ventricles. In healthy patients, the pressure of the left ventricle is considerably higher than the right, resulting in a septal wall that bows into the right ventricle. However, in patients with pulmonary hypertension, the pressure in the right ventricle increases significantly to a point where the pressure is similar to or surpasses that of the left ventricle during portions of the cardiac cycle. For these patients, the septal wall deviates towards the left ventricle, impacting its function. It is possible to study this effect using mathematical modeling, but existing models are nonlinear, leading to a system of algebraic differential equations that can be challenging to solve in patient-specific optimizations of clinical data. This study demonstrates that a simplified linearized model is sufficient to account for the effect of VVI and that, as expected, the impact is significantly more pronounced in patients with pulmonary hypertension.

Guidelines for mechanistic modeling and analysis in cardiovascular research

Computational, or in silico, models are an effective, noninvasive tool for investigating cardiovascular function. These models can be used in the analysis of experimental and clinical data to identify possible mechanisms of (ab)normal cardiovascular physiology. Recent advances in computing power and data management have led to innovative and complex modeling frameworks that simulate cardiovascular function across multiple scales. While commonly used in multiple disciplines, there is a lack of concise guidelines for the implementation of computer models in cardiovascular research. In line with recent calls for more reproducible research, it is imperative that scientists adhere to credible practices when developing and applying computational models to their research. The goal of this manuscript is to provide a consensus document that identifies best practices for in silico computational modeling in cardiovascular research. These guidelines provide the necessary methods for mechanistic model development, model analysis, and formal model calibration using fundamentals from statistics. We outline rigorous practices for computational, mechanistic modeling in cardiovascular research and discuss its synergistic value to experimental and clinical data.

Integrated modeling and simulation of recruitment of myocardial perfusion and oxygen delivery in exercise

While exercise-mediated vasoregulation in the myocardium is understood to be governed by autonomic, myogenic, and metabolic-mediated mechanisms, we do not yet understand the spatial heterogeneity of vasodilation or its effects on microvascular flow patterns and oxygen delivery. This study uses a simulation and modeling approach to explore the mechanisms underlying the recruitment of myocardial perfusion and oxygen delivery in exercise. The simulation approach integrates model components representing: whole-body cardiovascular hemodynamics, cardiac mechanics and myocardial work; myocardial perfusion; and myocardial oxygen transport. Integrating these systems together, model simulations reveal: (1.) To match expected flow and transmural flow ratios at increasing levels of exercise, a greater degree of vasodilation must occur in the subendocardium compared to the subepicardium. (2.) Oxygen extraction and venous oxygenation are predicted to substantially decrease with increasing exercise level preferentially in the subendocardium, suggesting that an oxygen-dependent error signal driving metabolic mediated recruitment of flow would be operative only in the subendocardium. (3.) Under baseline physiological conditions approximately 4% of the oxygen delivered to the subendocardium may be supplied via retrograde flow from coronary veins.

A multiscale finite element model of left ventricular mechanics incorporating baroreflex regulation

Cardiovascular function is regulated by a short-term hemodynamic baroreflex loop, which tries to maintain arterial pressure at a normal level. In this study, we present a new multiscale model of the cardiovascular system named MyoFE. This framework integrates a mechanistic model of contraction at the myosin level into a finite-element-based model of the left ventricle pumping blood through the systemic circulation. The model is coupled with a closed-loop feedback control of arterial pressure inspired by a baroreflex algorithm previously published by our team. The reflex loop mimics the afferent neuron pathway via a normalized signal derived from arterial pressure. The efferent pathway is represented by a kinetic model that simulates the net result of neural processing in the medulla and cell-level responses to autonomic drive. The baroreflex control algorithm modulates parameters such as heart rate and vascular tone of vessels in the lumped-parameter model of systemic circulation. In addition, it spatially modulates intracellular Ca2+ dynamics and molecular-level function of both the thick and the thin myofilaments in the left ventricle. Our study demonstrates that the baroreflex algorithm can maintain arterial pressure in the presence of perturbations such as acute cases of altered aortic resistance, mitral regurgitation, and myocardial infarction. The capabilities of this new multiscale model will be utilized in future research related to computational investigations of growth and remodeling.

Biventricular Interaction During Acute Left Ventricular Ischemia in Mice: A Combined In-Vivo and In-Silico Approach

Computational models provide an efficient paradigm for integrating and linking multiple spatial and temporal scales. However, these models are difficult to parameterize and match to experimental data. Recent advances in both data collection and model analyses have helped overcome this limitation. Here, we combine a multiscale, biventricular interaction model with mouse data before and after left ventricular (LV) ischemia. Sensitivity analyses are used to identify the most influential parameters on pressure and volume predictions. The subset of influential model parameters are calibrated to biventricular pressure-volume loop data (n = 3) at baseline. Each mouse underwent left anterior descending coronary artery ligation, during which changes in fractional shortening and RV pressure-volume dynamics were recorded. Using the calibrated model, we simulate acute LV ischemia and contrast outputs at baseline and in simulated ischemia. Our baseline simulations align with the LV and RV data, and our predictions during ischemia complement recorded RV data and prior studies on LV function during myocardial infarction. We show that a model with both biventricular mechanical interaction and systems-level cardiovascular dynamics can quantitatively reproduce in-vivo data and qualitatively match prior findings from animal studies on LV ischemia.

Computational modeling of ventricular-ventricular interactions suggest a role in clinical conditions involving heart failure

Introduction: The left (LV) and right (RV) ventricles are linked biologically, hemodynamically, and mechanically, a phenomenon known as ventricular interdependence. While LV function has long been known to impact RV function, the reverse is increasingly being realized to have clinical importance. Investigating ventricular interdependence clinically is challenging given the invasive measurements required, including biventricular catheterization, and confounding factors such as comorbidities, volume status, and other aspects of subject variability. Methods: Computational modeling allows investigation of mechanical and hemodynamic interactions in the absence of these confounding factors. Here, we use a threesegment biventricular heart model and simple circulatory system to investigate ventricular interdependence under conditions of systolic and diastolic dysfunction of the LV and RV in the presence of compensatory volume loading. We use the end-diastolic pressure-volume relationship, end-systolic pressure-volume relationship, Frank Starling curves, and cardiac power output as metrics. Results: The results demonstrate that LV systolic and diastolic dysfunction lead to RV compensation as indicated by increases in RV power. Additionally, RV systolic and diastolic dysfunction lead to impaired LV filling, interpretable as LV stiffening especially with volume loading to maintain systemic pressure. Discussion: These results suggest that a subset of patients with intact LV systolic function and diagnosed to have impaired LV diastolic function, categorized as heart failure with preserved ejection fraction (HFpEF), may in fact have primary RV failure. Application of this computational approach to clinical data sets, especially for HFpEF, may lead to improved diagnosis and treatment strategies and consequently improved outcomes.

Portal Venous Remodeling Determines the Pattern of Cirrhosis Decompensation: A Systems Analysis

INTRODUCTION

As liver disease progresses, scarring results in worsening hemodynamics ultimately culminating in portal hypertension. This process has classically been quantified through the portosystemic pressure gradient (PSG), which is clinically estimated by hepatic venous pressure gradient (HVPG); however, PSG alone does not predict a given patient's clinical trajectory regarding the Baveno stage of cirrhosis. We hypothesize that a patient's PSG sensitivity to venous remodeling could explain disparate disease trajectories.

METHODS

We created a computational model of the portal system in the context of worsening liver disease informed by physiologic measurements from the field of portal hypertension. We simulated progression of clinical complications, HVPG, and transjugular intrahepatic portosystemic shunt placement while only varying a patient's likelihood of portal venous remodeling.

RESULTS

Our results unify hemodynamics, venous remodeling, and the clinical progression of liver disease into a mathematically consistent model of portal hypertension. We find that by varying how sensitive patients are to create venous collaterals with rising PSG we can explain variation in patterns of decompensation for patients with liver disease. Specifically, we find that patients who have higher proportions of portosystemic shunting earlier in disease have an attenuated rise in HVPG, delayed onset of ascites, and less hemodynamic shifting after transjugular intrahepatic portosystemic shunt placement.

DISCUSSION

This article builds a computational model of portal hypertension which supports that patient-level differences in venous remodeling may explain disparate clinical trajectories of disease.

Biventricular interaction during acute left ventricular ischemia in mice: a combined in-vivo and in-silico approach

Computational models provide an efficient paradigm for integrating and linking multiple spatial and temporal scales. However, these models are difficult to parameterize and match to experimental data. Recent advances in both data collection and model analyses have helped overcome this limitation. Here, we combine a multiscale, biventricular interaction model with mouse data before and after left ventricular (LV) ischemia. Sensitivity analyses are used to identify the most influential parameters on pressure and volume predictions. The subset of influential model parameters are calibrated to biventricular pressure-volume loop data (n=3) at baseline. Each mouse underwent left anterior descending coronary artery ligation, during which changes in fractional shortening and RV pressure-volume dynamics were recorded. Using the calibrated model, we simulate acute LV ischemia and contrast outputs at baseline and in simulated ischemia. Our baseline simulations align with the LV and RV data, and our predictions during ischemia complement recorded RV data and prior studies on LV function during myocardial infarction. We show that a model with both biventricular mechanical interaction and systems level cardiovascular dynamics can quantitatively reproduce in-vivo data and qualitatively match prior findings from animal studies on LV ischemia.

A multiscale model of the cardiovascular system that regulates arterial pressure via closed loop baroreflex control of chronotropism, cell-level contractility, and vascular tone

Multiscale models of the cardiovascular system can provide new insights into physiological and pathological processes. PyMyoVent is a computer model that bridges from molecular- to organ-level function and which simulates a left ventricle pumping blood through the systemic circulation. Initial work with PyMyoVent focused on the end-systolic pressure volume relationship and ranked potential therapeutic strategies by their impact on contractility. This manuscript extends the PyMyoVent framework by adding closed-loop feedback control of arterial pressure. The control algorithm mimics important features of the physiological baroreflex and was developed as part of a long-term program that focuses on growth and biological remodeling. Inspired by the underlying biology, the reflex algorithm uses an afferent signal derived from arterial pressure to drive a kinetic model that mimics the net result of neural processing in the medulla and cell-level responses to autonomic drive. The kinetic model outputs control signals that are constrained between limits that represent maximum parasympathetic and maximum sympathetic drive and which modulate heart rate, intracellular Ca2+ dynamics, the molecular-level function of both the thick and the thin myofilaments, and vascular tone. Simulations show that the algorithm can regulate mean arterial pressure at user-defined setpoints as well as maintaining arterial pressure when challenged by changes in blood volume and/or valve resistance. The reflex also regulates arterial pressure when cell-level contractility is modulated to mimic the idealized impact of myotropes. These capabilities will be important for future work that uses computer modeling to investigate clinical conditions and treatments.

Aerobic Training for Healthy Men and Women: Determining Intensities by Different Equations

The aim of this study is to develop equations for aerobic exercise prescription for the intensities of 50, 60, 70 and 80% in healthy subjects of both sexes.

METHOD

This is a cross-sectional study with convenience sampling drawn from a database of 228 healthy subjects who were randomized into the regression group (GR: 197 subjects (male = 143 and female = 54)) and cross-validation group (CVG: 31 individuals [men = 20 and women = 11]). Kohavi's assumptions were followed in relation to cross-validation and bootstrap for precision estimation and model selection. The GR was used to build the estimation equations from the multiple linear regression. The CVG was determined to analyze the validity in the estimation equations. The equations to determine the intensities were constructed by means of multiple regression, the independent variables were determined by the stepwise method, observing the significance level of p < 0.05.

RESULTS

The reliability level of Cronbach's alpha of the multiple linear regression equations was moderate for the intensity of 50% (0.51); for the intensities of 60, 70 and 80%, it corresponded to 0.50, 0.53 and 0.57, respectively.

CONCLUSION

The results show that it is possible to apply the equations in the determination of aerobic exercise intensities for healthy individuals. However, the need for further studies in other populations to prove the reliability of the proposed equations is evident.