Search

Right ventricular (RV) function is an important prognostic indicator for pulmonary arterial hypertension (PAH), a vasculopathy that primarily and disproportionally affects women with distinct pre- and postmenopausal clinical outcomes. However, most animal studies have overlooked the impact of sex and ovarian hormones on RV remodeling in PAH. Here, we combined invasive measurements of RV hemodynamics and morphology with computational models of RV biomechanics in sugen-hypoxia (SuHx)-treated male, ovary-intact female, and ovariectomized female rats. Despite similar pressure overload levels, SuHx induced increases in end-diastolic elastance and passive myocardial stiffening, notably in male SuHx animals, corresponding to elevated diastolic intracellular calcium. Increases in end-systolic chamber elastance were largely explained by myocardial hypertrophy in male and ovary-intact female rats, whereas ovariectomized females exhibited contractility recruitment via calcium transient augmentation. Ovary-intact female rats primarily responded with hypertrophy, showing fewer myocardial mechanical alterations and less stiffening. These findings highlight sex-related RV remodeling differences in rats, affecting systolic and diastolic RV function in PAH. NEW & NOTEWORTHY Combining hemodynamic and morphological measurements from male, female, and ovariectomized female pulmonary arterial hypertension (PAH) rats revealed distinct adaptation mechanisms despite similar pressure overload. Males showed the most diastolic stiffening. Ovariectomized females had enhanced myocyte contractility and calcium transient upregulation. Ovary-intact females primarily responded with hypertrophy, experiencing milder passive myocardial stiffening and no changes in myocyte shortening. These findings suggest potential sex-specific pathways in right ventricular (RV) adaptation to PAH, with implications for targeted interventions.

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.

Determining whether an ectopic depolarization will lead to a self-perpetuating arrhythmia is of critical importance in determining arrhythmia risk, so it is necessary to understand what factors impact substrate vulnerability. This study sought to explore the impact of cell-to-cell heterogeneity in ion channel conductance on substrate vulnerability to arrhythmia by measuring the duration of the vulnerable window in computational models of one-dimensional cables of ventricular cardiomyocytes. We began by using a population of uniform cable models to determine the mechanisms underlying the vulnerable window phenomenon. We found that in addition to the known importance of G Na , the conductances G Ca,L and G Kr also play a minor role in determining the vulnerable window duration. We also found that a steeper slope of the repolarizing action potential during the vulnerable window correlated with a shorter vulnerable window duration in uniform cables. We applied our understanding from these initial simulations to an investigation of the vulnerable window in heterogeneous cable models. The heterogeneous cables displayed a great deal of intra-cable variation in vulnerable window duration, highly sensitive to the cardiomyocytes in the local environment of the ectopic stimulus. Coupling strength modulated not only the magnitude of the vulnerable window duration but also the extent of intra-tissue variability in vulnerable window duration. NEW & NOTEWORTHY We investigate the impact of cell-to-cell heterogeneity in ion channel conductance on substrate vulnerability to arrhythmia by measuring the vulnerable window duration in computational cardiomyocyte cable models. We demonstrate a wide range of intra-cable variability in vulnerable window duration (VWD) and show how this is changed by ion channel block and coupling strength perturbations.

Murray's law has been viewed as a fundamental law of physiology. Relating blood flow ([Formula: see text]) to vessel diameter ( D ) ([Formula: see text]·∝· D 3 ), it dictates minimum lumen area (MLA) targets for coronary bifurcation percutaneous coronary intervention (PCI). The cubic exponent (3.0), however, has long been disputed, with alternative theoretical derivations, arguing this should be closer to 2.33 (7/3). The aim of this meta-analysis was to quantify the optimum flow-diameter exponent in human and mammalian coronary arteries. We conducted a systematic review and meta-analysis of all articles quantifying an optimum flow-diameter exponent for mammalian coronary arteries within the Cochrane library, PubMed Medline, Scopus, and Embase databases on 20 March 2023. A random-effects meta-analysis was used to determine a pooled flow-diameter exponent. Risk of bias was assessed with the National Institutes of Health (NIH) quality assessment tool, funnel plots, and Egger regression. From a total of 4,772 articles, 18 were suitable for meta-analysis. Studies included data from 1,070 unique coronary trees, taken from 372 humans and 112 animals. The pooled flow diameter exponent across both epicardial and transmural arteries was 2.39 (95% confidence interval: 2.24-2.54; I 2 = 99%). The pooled exponent of 2.39 showed very close agreement with the theoretical exponent of 2.33 (7/3) reported by Kassab and colleagues. This exponent may provide a more accurate description of coronary morphometric scaling in human and mammalian coronary arteries, as compared with Murray's original law. This has important implications for the assessment, diagnosis, and interventional treatment of coronary artery disease.

Premature ventricular complexes (PVCs) are spontaneous excitations occurring in the ventricles of the heart that are associated with ventricular arrhythmias and sudden cardiac death. Under long QT conditions, PVCs can be mediated by repolarization gradient (RG) and early afterdepolarizations (EADs), yet the effects of heterogeneities or geometry of the RG or EAD regions on PVC genesis remain incompletely understood. In this study, we use computer simulation to systematically investigate the effects of the curvature of the RG border region on PVC genesis under long QT conditions. We show that PVCs can be either promoted or suppressed by negative or positive RG border curvature depending on the source and sink conditions. When the origin of oscillation is in the source region and the source is too strong, a positive RG border curvature can promote PVCs by causing the source area to oscillate. When the origin of oscillation is in the sink region, a negative RG border curvature can promote PVCs by causing the sink area to oscillate. Furthermore, EAD-mediated PVCs are also promoted by negative border curvature. We also investigate the effects of wavefront curvature and show that PVCs are promoted by convex but suppressed by concave wavefronts; however, the effect of wavefront curvature is much smaller than that of RG border curvature. In conclusion, besides the increase of RG and occurrence of EADs caused by QT prolongation, the geometry of the RG border plays important roles in PVC genesis, which can greatly increase the risk of arrhythmias in cardiac diseases. NEW & NOTEWORTHY The effects of the curvature or geometry of the repolarization gradient region and wavefront curvature on the genesis of premature ventricular complexes are systematically investigated using computer modeling and simulation. Premature ventricular complexes can be promoted by either positive or negative curvature of the gradient region depending on the source and sink conditions. The underlying mechanisms of the curvature effects are revealed, which provides mechanistic insights into arrhythmogenesis in cardiac diseases.

Heart growth in the pregnant patient helps maintain cardiovascular function while supporting the growing fetus. However, in some cases, the cardiovascular demand of pregnancy can trigger life-threatening conditions, including hypertensive disorders of pregnancy and peripartum cardiomyopathy. The mechanisms that control heart growth throughout pregnancy are unclear, and treating these diseases remains elusive. We previously developed a computational model that accounts for hormonal and hemodynamic interactions throughout pregnancy and demonstrated its ability to capture realistic cardiac growth in normal rat pregnancy. In this study, we evaluated whether this model could capture heart growth beyond normal pregnancy. After further validation of our normal pregnancy predictions, we tested our model predictions of three rat studies of hypertensive pregnancies. We next simulated the postpartum period and examined the impact of lactation on cardiac growth in rats. We demonstrate that our multiscale model can capture cardiac growth associated with new-onset hypertension during pregnancy and lactation status in the postpartum period. We conclude by elaborating on the potential clinical utility of our model in the future. NEW & NOTEWORTHY Our multiscale model predicts appropriate heart growth beyond normal pregnancy, including elevated heart weights in rats with induced hypertension during pregnancy and in lactating mice and decreased heart weight in nonlactating mice. Our model captures distinct mechanisms that result in similar organ-level growth, highlighting its potential to distinguish healthy from diseased pregnancy-induced growth.

Arterial pulse waves (PWs) such as blood pressure and photoplethysmogram (PPG) signals contain a wealth of information on the cardiovascular (CV) system that can be exploited to assess vascular age and identify individuals at elevated CV risk. We review the possibilities, limitations, complementarity, and differences of reduced-order, biophysical models of arterial PW propagation, as well as theoretical and empirical methods for analyzing PW signals and extracting clinically relevant information for vascular age assessment. We provide detailed mathematical derivations of these models and theoretical methods, showing how they are related to each other. Finally, we outline directions for future research to realize the potential of modeling and analysis of PW signals for accurate assessment of vascular age in both the clinic and in daily life.

Cardiac voltage-gated sodium channel gain-of-function prolongs repolarization in the long-QT syndrome type 3 (LQT3). Previous studies suggest that narrowing the perinexus within the intercalated disc, leading to rapid sodium depletion, attenuates LQT3-associated action potential duration (APD) prolongation. However, it remains unknown whether extracellular sodium concentration modulates APD prolongation during sodium channel gain-of-function. We hypothesized that elevated extracellular sodium concentration and widened perinexus synergistically prolong APD in LQT3. LQT3 was induced with sea anemone toxin (ATXII) in Langendorff-perfused guinea pig hearts ( n = 34). Sodium concentration was increased from 145 to 160 mM. Perinexal expansion was induced with mannitol or the sodium channel β1-subunit adhesion domain antagonist (βadp1). Epicardial ventricular action potentials were optically mapped. Individual and combined effects of varying clefts and sodium concentrations were simulated in a computational model. With ATXII, both mannitol and βadp1 significantly widened the perinexus and prolonged APD, respectively. The elevated sodium concentration alone significantly prolonged APD as well. Importantly, the combination of elevated sodium concentration and perinexal widening synergistically prolonged APD. Computational modeling results were consistent with animal experiments. Concurrently elevating extracellular sodium and increasing intercalated disc edema prolongs repolarization more than the individual interventions alone in LQT3. This synergistic effect suggests an important clinical implication that hypernatremia in the presence of cardiac edema can markedly increase LQT3-associated APD prolongation. Therefore, to our knowledge, this is the first study to provide evidence of a tractable and effective strategy to mitigate LQT3 phenotype by means of managing sodium levels and preventing cardiac edema in patients. NEW & NOTEWORTHY This is the first study to demonstrate that the long-QT syndrome type 3 (LQT3) phenotype can be exacerbated or concealed by regulating extracellular sodium concentrations and/or the intercalated disc separation. The animal experiments and computational modeling in the current study reveal a critically important clinical implication: sodium dysregulation in the presence of edema within the intercalated disc can markedly increase the risk of arrhythmia in LQT3. These findings strongly suggest that maintaining extracellular sodium within normal physiological limits may be an effective and inexpensive therapeutic option for patients with congenital or acquired sodium channel gain-of-function diseases.

Myocardial function and exercise reserve are important determinants of outcome in pulmonary arterial hypertension (PAH) but are incompletely understood. For this study, we performed subject-specific computer simulations, based on invasive measurements and cardiac magnetic resonance imaging (CMR), to investigate whole circulation properties in PAH at rest and exercise and determinants of exercise reserve. CMR and right heart catheterization were performed in nine patients with idiopathic PAH, and CMR in 10 healthy controls. CMR during exercise was performed in seven patients with PAH. A full-circulation computer model was developed, and model parameters were optimized at the individual level. Patient-specific simulations were used to analyze the effect of right ventricular (RV) inotropic reserve on exercise performance. Simulations achieved a high consistency with observed data. RV contractile force was increased in patients with PAH (127.1 ± 28.7 kPa vs. 70.5 ± 14.5 kPa, P < 0.001), whereas left ventricular contractile force was reduced (107.5 ± 17.5 kPa vs. 133.9 ± 10.3 kPa, P = 0.002). During exercise, RV contractile force increased by 1.56 ± 0.17, P = 0.001. In silico experiments confirmed RV inotropic reserve as the important limiting factor for cardiac output. Subject-specific computer simulation of myocardial mechanics in PAH is feasible and can be used to evaluate myocardial performance. With this method, we demonstrate marked functional myocardial adaptation to PAH in the resting state, primarily composed of increased contractile force development by RV myofibers, and we show the negative impact of reduced RV inotropic reserve on cardiac output during exercise. NEW & NOTEWORTHY Computer simulations of the myocardial mechanics and hemodynamics of rest and exercise were performed in nine patients with pulmonary arterial hypertension and 10 control subjects, with the use of data from invasive catheterization and from cardiac magnetic resonance. This approach allowed a detailed analysis of myocardial adaptation to pulmonary arterial hypertension and showed how reduction in right ventricular inotropic reserve is the important limiting factor for an increase in cardiac output during exercise.

Although pulmonary arterial hypertension (PAH) leads to right ventricle (RV) hypertrophy and structural remodeling, the relative contributions of changes in myocardial geometric and mechanical properties to systolic and diastolic chamber dysfunction and their time courses remain unknown. Using measurements of RV hemodynamic and morphological changes over 10 wk in a male rat model of PAH and a mathematical model of RV mechanics, we discriminated the contributions of RV geometric remodeling and alterations of myocardial material properties to changes in systolic and diastolic chamber function. Significant and rapid RV hypertrophic wall thickening was sufficient to stabilize ejection fraction in response to increased pulmonary arterial pressure by week 4 without significant changes in systolic myofilament activation. After week 4 , RV end-diastolic pressure increased significantly with no corresponding changes in end-diastolic volume. Significant RV diastolic chamber stiffening by week 5 was not explained by RV hypertrophy. Instead, model analysis showed that the increases in RV end-diastolic chamber stiffness were entirely attributable to increased resting myocardial material stiffness that was not associated with significant myocardial fibrosis or changes in myocardial collagen content or type. These findings suggest that whereas systolic volume in this model of RV pressure overload is stabilized by early RV hypertrophy, diastolic dilation is prevented by subsequent resting myocardial stiffening. NEW & NOTEWORTHY Using a novel combination of hemodynamic and morphological measurements over 10 wk in a male rat model of PAH and a mathematical model of RV mechanics, we found that compensated systolic function was almost entirely explained by RV hypertrophy, but subsequently altered RV end-diastolic mechanics were primarily explained by passive myocardial stiffening that was not associated with significant collagen extracellular matrix accumulation.

Ventricular arrhythmias, consisting of single ectopic beats (sEB), multiple EB (mEB), and torsades de pointes (TdP, defined as ≥5 beats with QRS vector twisting around isoelectric line) can be induced in the anesthetized chronic atrioventricular block (CAVB) dog by dofetilide ( I Kr blocker). The interplay between temporal dispersion of repolarization, quantified as short-term variability (STV), and spatial dispersion of repolarization (SDR) in the initiation and perpetuation of these arrhythmias remains unclear. Five inducible (≥3 TdPs/10 min) CAVB dogs underwent one mapping experiment and were observed for 10 min from the start of dofetilide infusion (0.025 mg/kg, 5 min). An intracardiac decapolar electrogram (EGM) catheter and 30 intramural cardiac needles in the left ventricle (LV) were introduced. STV ARI was derived from 31 consecutive activation recovery intervals (ARIs) on the intracardiac EGM, using the formula: [Formula: see text]. The mean SDR 3D in the LV was determined as the three-dimensional repolarization time differences between the intramural cardiac needles. Moments of measurement included baseline (BL) and after dofetilide infusion before first 1 ) sEB (occurrence at 100 ± 35 s), 2 ) mEB (224 ± 96 s), and 3 ) non-self-terminating TdP (454 ± 298 s). STV ARI increased from 2.15 ± 0.32 ms at BL to 3.73 ± 0.99 ms* before the first sEB and remained increased without further significant progression to mEB (4.41 ± 0.45 ms*) and TdP (5.07 ± 0.84 ms*) (* P < 0.05 compared with BL). SDR 3D did not change from 31 ± 11 ms at BL to 43 ± 13 ms before sEB but increased significantly before mEB (68 ± 7 ms*) and to TdP (86 ± 9 ms* + ) ( + P < 0.05 compared with sEB). An increase in STV contributes to the initiation of sEB, whereas an increase in SDR is important for the perpetuation of non-self-terminating TdPs. NEW & NOTEWORTHY This study compared two well-established electrophysiological parameters, being temporal and spatial dispersion of repolarization, and provided new insights into their interplay in the arrhythmogenesis of torsades de pointes arrhythmias. Although it confirmed that an increase in temporal dispersion of repolarization contributes to the initiation of single ectopic beats, it showed that an increase in spatial dispersion of repolarization is important for the perpetuation of non-self-terminating torsades de pointes arrhythmias.

Chronic thromboembolic pulmonary hypertension (CTEPH) is caused by recurrent or unresolved pulmonary thromboemboli, leading to perfusion defects and increased arterial wave reflections. CTEPH treatment aims to reduce pulmonary arterial pressure and reestablish adequate lung perfusion, yet patients with distal lesions are inoperable by standard surgical intervention. Instead, these patients undergo balloon pulmonary angioplasty (BPA), a multisession, minimally invasive surgery that disrupts the thromboembolic material within the vessel lumen using a catheter balloon. However, there still lacks an integrative, holistic tool for identifying optimal target lesions for treatment. To address this insufficiency, we simulate CTEPH hemodynamics and BPA therapy using a multiscale fluid dynamics model. The large pulmonary arterial geometry is derived from a computed tomography (CT) image, whereas a fractal tree represents the small vessels. We model ring- and web-like lesions, common in CTEPH, and simulate normotensive conditions and four CTEPH disease scenarios; the latter includes both large artery lesions and vascular remodeling. BPA therapy is simulated by simultaneously reducing lesion severity in three locations. Our predictions mimic severe CTEPH, manifested by an increase in mean proximal pulmonary arterial pressure above 20 mmHg and prominent wave reflections. Both flow and pressure decrease in vessels distal to the lesions and increase in unobstructed vascular regions. We use the main pulmonary artery (MPA) pressure, a wave reflection index, and a measure of flow heterogeneity to select optimal target lesions for BPA. In summary, this study provides a multiscale, image-to-hemodynamics pipeline for BPA therapy planning for patients with inoperable CTEPH. NEW & NOTEWORTHY This article presents novel computational framework for predicting pulmonary hemodynamics in chronic thromboembolic pulmonary hypertension. The mathematical model is used to identify the optimal target lesions for balloon pulmonary angioplasty, combining simulated pulmonary artery pressure, wave intensity analysis, and a new quantitative metric of flow heterogeneity.

Radial applanation tonometry is a well-established method for clinical hemodynamic assessment and is also becoming popular in wrist-worn fitness trackers. The time difference between the foot and the dicrotic notch of the arterial pressure waveform is a well-accepted approximation for the left ventricular ejection time (ET). However, several clinical studies have shown that ET measured from the radial pressure waveform deviates from that measured centrally. In this work, we consider the systolic wave and the dicrotic wave as two independent traveling waves and hypothesize that their wave speed difference leads to the intersite differences of measured ET (ΔET). Accordingly, we derived a mathematical dicrotic wave decomposition model and identified the most influential factors on ΔET via global sensitivity analysis. In our clinical validation on a heterogeneous cohort ( N = 5,742) from the Framingham Heart Study (FHS), the local sensitivity analysis results resembled the sensitivity variation patterns of ΔET from model simulations. A regression analysis on FHS data, using morphological features of radial pressure waveforms to estimate the carotid ET, produced a root mean square error of 3.76 ms and R 2 of 0.91. The proposed dicrotic wave decomposition model can explain the intersite ET measurement discrepancies observed in the clinical data of FHS and can facilitate the precise identification of ET with radial pressure waveforms. Therefore, the proposed model will improve various physics-based pulse wave analysis methods as well as prospective artificial intelligence methods for tackling the subsequent big data produced from widespread wearable radial pressure monitoring. NEW & NOTEWORTHY Based on a new understanding of pressure wave propagation, we propose a novel dicrotic wave decomposition model considering the dicrotic wave as an independent traveling component. The proposed model can explain the mechanism underlying the intersite discrepancies in ejection time measurement from arterial waveforms and then, in principle, enhance the accuracy of both classical physics-based as well as more contemporary artificial intelligence-based pulse wave analysis methods in clinical and wearable radial blood pressure monitoring applications.

Another type of natural wave, traced from longitudinal wall motion and propagation along the artery, is observed in our in vivo human carotid artery experiments. We coin it as extension wave (EW) and hypothesize that EW velocity (EWV) is associated with arterial longitudinal stiffness. The EW is thus assumed to complement the pulse wave (PW), whose velocity (PWV) is tracked from the radial wall displacement and linked to arterial circumferential stiffness through the Moens-Korteweg equation, as indicators for arterial mechanical anisotropy quantification by noninvasive high-frame-rate ultrasound. The relationship between directional arterial stiffnesses and the two natural wave speeds was investigated in wave theory, finite-element simulations based on isotropic and anisotropic arterial models, and in vivo human common carotid artery ( n = 10) experiments. Excellent agreement between the theory and simulations showed that EWV was 2.57 and 1.03 times higher than PWV in an isotropic and an anisotropic carotid artery model, respectively, whereas in vivo EWV was consistently lower than PWV in all 10 healthy human subjects. A strong linear correlation was substantiated in vivo between EWV and arterial longitudinal stiffness quantified by a well-validated vascular-guided wave imaging technique (VGWI). We thereby proposed a novel index calculated as EWV 2 /PWV 2 as an alternative to assess arterial mechanical anisotropy. Simulations and in vivo results corroborated the effect of mechanical anisotropy on the propagation of spontaneous waves along the arterial wall. The proposed anisotropy index demonstrated the feasibility of the concurrent EW and PW imaged by high frame-rate ultrasound in grading of arterial wall anisotropy. NEW & NOTEWORTHY An extension wave formed by longitudinal wall displacements was observed by high-frame-rate ultrasound in the human common carotid artery in vivo. A strong correlation between extension wave velocity and arterial longitudinal stiffness complements the well-established pulse wave, which is linked to circumferential stiffness, to noninvasively assess direction-dependent wall elasticity of the major artery. The proposed anisotropy index, which directly reflects arterial wall microstructure and function, might be a potential risk factor for screening (sub-) clinical cardiovascular diseases.

Cardiac myosin binding protein-C (cMyBP-C) is a thick filament protein that influences sarcomere stiffness and modulates cardiac contraction-relaxation through its phosphorylation. Phosphorylation of cMyBP-C and ablation of cMyBP-C have been shown to increase the rate of MgADP release in the acto-myosin cross-bridge cycle in the intact sarcomere. The influence of cMyBP-C on Pi-dependent myosin kinetics has not yet been examined. We investigated the effect of cMyBP-C, and its phosphorylation, on myosin kinetics in demembranated papillary muscle strips bearing the β-cardiac myosin isoform from nontransgenic and homozygous transgenic mice lacking cMyBP-C. We used quick stretch and stochastic length-perturbation analysis to characterize rates of myosin detachment and force development over 0-12 mM Pi and at maximal (pCa 4.8) and near-half maximal (pCa 5.75) Ca 2+ activation. Protein kinase A (PKA) treatment was applied to half the strips to probe the effect of cMyBP-C phosphorylation on Pi sensitivity of myosin kinetics. Increasing Pi increased myosin cross-bridge detachment rate similarly for muscles with and without cMyBP-C, although these rates were higher in muscle without cMyBP-C. Treating myocardial strips with PKA accelerated detachment rate when cMyBP-C was present over all Pi, but not when cMyBP-C was absent. The rate of force development increased with Pi in all muscles. However, Pi sensitivity of the rate force development was reduced when cMyBP-C was present versus absent, suggesting that cMyBP-C inhibits Pi-dependent reversal of the power stroke or stabilizes cross-bridge attachment to enhance the probability of completing the power stroke. These results support a functional role for cMyBP-C in slowing myosin detachment rate, possibly through a direct interaction with myosin or by altering strain-dependent myosin detachment via cMyBP-C-dependent stiffness of the thick filament and myofilament lattice. PKA treatment reduces the role for cMyBP-C to slow myosin detachment and thus effectively accelerates β-myosin detachment in the intact myofilament lattice. NEW & NOTEWORTHY Length perturbation analysis was used to demonstrate that β-cardiac myosin characteristic rates of detachment and recruitment in the intact myofilament lattice are accelerated by Pi, phosphorylation of cMyBP-C, and the absence of cMyBP-C. The results suggest that cMyBP-C normally slows myosin detachment, including Pi-dependent detachment, and that this inhibition is released with phosphorylation or absence of cMyBP-C.

Atherosclerosis is a dynamic process starting with endothelial dysfunction and inflammation and eventually leading to life-threatening arterial plaques. Exercise generally improves endothelial function in a dose-dependent manner by altering hemodynamics, specifically by increased arterial pressure, pulsatility, and shear stress. However, athletes who regularly participate in high-intensity training can develop arterial plaques, suggesting alternative mechanisms through which excessive exercise promotes vascular disease. Understanding the mechanisms that drive atherosclerosis in sedentary versus exercise states may lead to novel rehabilitative methods aimed at improving exercise compliance and physical activity. Preclinical tools, including in vitro cell assays, in vivo animal models, and in silico computational methods, broaden our capabilities to study the mechanisms through which exercise impacts atherogenesis, from molecular maladaptation to vascular remodeling. Here, we describe how preclinical research tools have and can be used to study exercise effects on atherosclerosis. We then propose how advanced bioengineering techniques can be used to address gaps in our current understanding of vascular pathophysiology, including integrating in vitro, in vivo, and in silico studies across multiple tissue systems and size scales. Improving our understanding of the antiatherogenic exercise effects will enable engaging, targeted, and individualized exercise recommendations to promote cardiovascular health rather than treating cardiovascular disease that results from a sedentary lifestyle.

Central systolic blood pressure (cSBP, the peak of the central waveform) is usually regarded as the determinant of peripheral systolic blood pressure with amplification of peripheral systolic blood pressure (pSBP) measured with reference to cSBP. However, the earlier portion of the central waveform up to the first systolic shoulder (P1) may be the major determinant of pSBP. We performed in silico simulation studies and examined previously acquired experimental data ( n = 131) in which peripheral and central blood pressure waveforms had been acquired both invasively and noninvasively to examine the determinants of pSBP. Measurements were made at baseline and during perturbation of hemodynamics by inotropic and vasoactive drugs. In silico simulations using a central-to-peripheral transfer function demonstrated that pSBP is dependent on P1 and the rate of change (dP/d t ) of central pressure up to the time of P1 but not cSBP. In computational simulations, peripheral reflection in the radial artery was closely related to dP/d t , and 97% of the variability in amplification as measured with reference to P1 was explained by dP/d t . In vivo, amplification of pSBP over P1 was correlated with dP/d t ( R  > 0.75, P < 0.0001 for all data sets), and P1 and dP/d t were independently correlated with pSBP, explaining 90% of the variability in pSBP. We conclude that P1 and dP/d t are major determinants of pSBP and that pSBP and cSBP are, in part, determined by different cardiac, central, and peripheral vascular properties. NEW & NOTEWORTHY Peripheral systolic BP is determined mainly by the first shoulder and the rate of rise of the central systolic blood pressure waveform rather than the peak of this waveform (central systolic BP). Peripheral and central systolic blood pressure are determined by different cardiac and vascular properties.

Accurate assessment of the left ventricular (LV) systolic function is indispensable in the clinic. However, estimation of a precise index of cardiac contractility, i.e., the end-systolic elastance ( E es ), is invasive and cannot be established as clinical routine. The aim of this work was to present and validate a methodology that allows for the estimation of E es from simple and readily available noninvasive measurements. The method is based on a validated model of the cardiovascular system and noninvasive data from arm-cuff pressure and routine echocardiography to render the model patient-specific. Briefly, the algorithm first uses the measured aortic flow as model input and optimizes the properties of the arterial system model to achieve correct prediction of the patient's peripheral pressure. In a second step, the personalized arterial system is coupled with the cardiac model (time-varying elastance model) and the LV systolic properties, including E es , are tuned to predict accurately the aortic flow waveform. The algorithm was validated against invasive measurements of E es (multiple pressure-volume loop analysis) taken from n = 10 patients with heart failure with preserved ejection fraction and n = 9 patients without heart failure. Invasive measurements of E es (median = 2.4 mmHg/mL, range = [1.0, 5.0] mmHg/mL) agreed well with method predictions (normalized root mean square error = 9%, ρ = 0.89, bias = -0.1 mmHg/mL, and limits of agreement = [-0.9, 0.6] mmHg/mL). This is a promising first step toward the development of a valuable tool that can be used by clinicians to assess systolic performance of the LV in the critically ill. NEW & NOTEWORTHY In this study, we present a novel model-based method to estimate the left ventricular (LV) end-systolic elastance ( E es ) according to measurement of the patient's arm-cuff pressure and a routine echocardiography examination. The proposed method was validated in vivo against invasive multiple-loop measurements of E es , achieving high correlation and low bias. This tool could be most valuable for clinicians to assess the cardiovascular health of critically ill patients.

The use of physiological models in medicine allows the evaluation of new hypotheses, development of diagnosis and clinical treatment applications, and development of training and medical education tools, as well as medical device design. Although several mathematical models of physiological systems have been presented in the literature, few of them are able to predict the human cardiorespiratory response under physical exercise stimulus adequately. This paper aims to present the building and comparison of an integrated cardiorespiratory model focused on the prediction of the healthy human response under rest and aerobic exercise. The model comprises cardiovascular circulation, respiratory mechanics, and gas exchange system, as well as cardiovascular and respiratory controllers. Every system is based on previously reported physiological models and incorporates reported mechanisms related to the aerobic exercise dynamics. Experimental data of 30 healthy male volunteers undergoing a cardiopulmonary exercise test and simulated data from two of the most current and complete cardiorespiratory models were used to evaluate the performance of the presented model. Experimental design, processing, and exploratory analysis are described in detail. The simulation results were compared against the experimental data in steady state and in transient regime. The predictions of the proposed model closely mimic the experimental data, showing in overall the lowest prediction error (10.35%), the lowest settling times for cardiovascular and respiratory variables, and in general the fastest and similar responses in transient regime. These results suggest that the proposed model is suitable to predict the cardiorespiratory response of healthy adult humans under rest and aerobic exercise conditions. NEW & NOTEWORTHY This paper presents a new cardiorespiratory model focused on the prediction of the healthy human response under rest and aerobic dynamic exercise conditions. Simulation results of cardiorespiratory variables are compared against experimental data and two of the most current and complete cardiorespiratory models.

Mechanical dyssynchrony (MD) affects left ventricular (LV) mechanics and coronary perfusion. To understand the multifactorial effects of MD, we developed a computational model that bidirectionally couples the systemic circulation with the LV and coronary perfusion with flow regulation. In the model, coronary flow in the left anterior descending (LAD) and left circumflex (LCX) arteries affects the corresponding regional contractility based on a prescribed linear LV contractility-coronary flow relationship. The model is calibrated with experimental measurements of LV pressure and volume, as well as LAD and LCX flow rate waveforms acquired under regulated and fully dilated conditions from a swine under right atrial (RA) pacing. The calibrated model is applied to simulate MD. The model can simultaneously reproduce the reduction in mean LV pressure (39.3%), regulated flow (LAD: 7.9%; LCX: 1.9%), LAD passive flow (21.6%), and increase in LCX passive flow (15.9%). These changes are associated with right ventricular pacing compared with RA pacing measured in the same swine only when LV contractility is affected by flow alterations with a slope of 1.4 mmHg/mL 2 in a contractility-flow relationship. In sensitivity analyses, the model predicts that coronary flow reserve (CFR) decreases and increases in the LAD and LCX with increasing delay in LV free wall contraction. These findings suggest that asynchronous activation associated with MD impacts 1 ) the loading conditions that further affect the coronary flow, which may explain some of the changes in CFR, and 2 ) the coronary flow that reduces global contractility, which contributes to the reduction in LV pressure. NEW & NOTEWORTHY A computational model that couples the systemic circulation of the left ventricular (LV) and coronary perfusion with flow regulation is developed to study the effects of mechanical dyssynchrony. The delayed contraction in the LV free wall with respect to the septum has a significant effect on LV function and coronary flow reserve.

The TRPV4 channel is a calcium-permeable channel ( P Ca / P Na ∼ 10). Its expression has been reported in ventricular myocytes, where it is involved in several cardiac pathological mechanisms. In this study, we investigated the implication of TRPV4 in ventricular electrical activity. Left ventricular myocytes were isolated from trpv4 +/+ and trpv4 -/- mice. TRPV4 membrane expression and its colocalization with L-type calcium channels (Ca v 1.2) was confirmed using Western blot biotinylation, immunoprecipitation, and immunostaining experiments. Then, electrocardiograms (ECGs) and patch-clamp recordings showed shortened QTc and action potential (AP) duration in trpv4 -/- compared with trpv4 +/+ mice. Thus, TRPV4 activator GSK1016790A produced a transient and dose-dependent increase in AP duration at 90% of repolarization (APD 90 ) in trpv4 +/+ but not in trpv4 -/- myocytes or when combined with TRPV4 inhibitor GSK2193874 (100 nM). Hence, GSK1016790A increased calcium transient (CaT) amplitude in trpv4 +/+ but not in trpv4 -/- myocytes, suggesting that TRPV4 carries an inward Ca 2+ current in myocytes. Conversely, TRPV4 inhibitor GSK2193874 (100 nM) alone reduced APD 90 in trpv4 +/+ but not in trpv4 -/- myocytes, suggesting that TRPV4 prolongs AP duration in basal condition. Finally, introducing TRPV4 parameters in a mathematical model predicted the development of an inward TRPV4 current during repolarization that increases AP duration and CaT amplitude, in accord with what was found experimentally. This study shows for the first time that TRPV4 modulates AP and QTc durations. It would be interesting to evaluate whether TRPV4 could be involved in long QT-mediated ventricular arrhythmias. NEW & NOTEWORTHY Transient receptor potential vanilloid 4 (TRPV4) is expressed at the membrane of mouse ventricular myocytes and colocalizes with non-T-tubular L-type calcium channels. Deletion of trpv4 gene in mice results in shortened QT interval on electrocardiogram and reduced action potential duration of ventricular myocytes. Pharmacological activation of TRPV4 channel leads to increased action potential duration and increased calcium transient amplitude in trpv4 -/- but not in trpv4 -/- ventricular myocytes. To the contrary, TRPV4 channel pharmacological inhibition reduces action potential duration in trpv4 +/+ but not in trpv4 -/- myocytes. Integration of TRPV4 channel in a computational model of mouse action potential shows that the channel carries an inward current contributing to slowing down action potential repolarization and to increase calcium transient amplitude, similarly to what is observed experimentally. This study highlights for the first time the involvement of TRPV4 channel in ventricular electrical activity.

Comprehensive functional characterization of cardiac tissue includes investigation of length and load dependence. Such measurements have been slow to develop in engineered heart tissues (EHTs), whose mechanical characterizations have been limited primarily to isometric and near-isometric behaviors. A more realistic assessment of myocardial function would include force-velocity curves to characterize power output and force-length loops mimicking the cardiac cycle to characterize work output. We developed a system that produces force-velocity curves and work loops in human EHTs using an adaptive iterative control scheme. We used human EHTs in this system to perform a detailed characterization of the cardiac β-myosin specific inhibitor, mavacamten. Consistent with the clinically proposed application of this drug to treat hypertrophic cardiomyopathy, our data support the premise that mavacamten improves diastolic function through reduction of diastolic stiffness and isometric relaxation time. Meanwhile, the effects of mavacamten on length- and load-dependent muscle performance were mixed. The drug attenuated the length-dependent response at small stretch values but showed normal length dependency at longer lengths. Peak power output of mavacamten-treated EHTs showed reduced power output as expected but also shifted peak power output to a lower load. Here, we demonstrate a robust method for the generation of isotonic contraction series and work loops in engineered heart tissues using an adaptive-iterative method. This approach reveals new features of mavacamten pharmacology, including previously unappreciated effects on intrinsic myosin dynamics and preservation of Frank-Starling behavior at longer muscle lengths. NEW & NOTEWORTHY We applied innovative methods to comprehensively characterize the length and load-dependent behaviors of engineered human cardiac muscle when treated with the cardiac β-myosin specific inhibitor mavacamten, a drug on the verge of clinical implementation for hypertrophic cardiomyopathy. We find mechanistic support for the role of mavacamten in improving diastolic function of cardiac tissue and note novel effects on work and power.

Ever since their origin more than one half-century ago, microneurographic recordings of sympathetic nerve activity have significantly advanced our understanding of the generation and regulation of central sympathetic outflow in human health and disease. For example, it is now appreciated that a myriad of disease states exhibit chronic sympathetic overactivity, a significant predictor of cardiovascular morbidity and mortality. Although microneurographic recordings allow for the direct quantification of sympathetic outflow, they alone do not provide information with respect to the ensuing sympathetically mediated vasoconstriction and blood pressure (BP) response. Therefore, the study of vascular and/or BP responses to sympathetic outflow (i.e., sympathetic transduction) has now emerged as an area of growing interest within the field of neural cardiovascular control in human health and disease. To date, studies have primarily examined sympathetic transduction under two distinct paradigms: when reflexively evoking sympatho-excitation through the induction of a laboratory stressor (i.e., sympathetic transduction during stress) and/or following spontaneous bursts of sympathetic outflow occurring under resting conditions (i.e., sympathetic transduction at rest). The purpose of this brief review is to highlight how our physiological understanding of sympathetic transduction has been advanced by these studies and to evaluate the primary analytical techniques developed to study sympathetic transduction in humans. We also discuss the framework by which the assessment of sympathetic transduction during stress reflects a fundamentally different process relative to sympathetic transduction at rest and why findings from investigations using these different techniques should be interpreted as such and not necessarily be considered one and the same.

With each heartbeat, the right ventricle (RV) inputs blood into the pulmonary vascular (PV) compartment, which conducts blood through the lungs at low pressure and concurrently fills the left atrium (LA) for output to the systemic circulation. This overall hemodynamic function of the integrated RV-PV-LA unit is determined by complex interactions between the components that vary over the cardiac cycle but are often assessed in terms of mean pressure and flow. Exercise challenges these hemodynamic interactions as cardiac filling increases, stroke volume augments, and cycle length decreases, with PV pressures ultimately increasing in association with cardiac output. Recent cardiopulmonary exercise hemodynamic studies have enriched the available data from healthy adults, yielded insight into the underlying mechanisms that modify the PV pressure-flow relationship, and better delineated the normal limits of healthy responses to exercise. This review will examine hemodynamic function of the RV-PV-LA unit using the two-element Windkessel model for the pulmonary circulation. It will focus on acute PV and LA responses that accommodate increased RV output during exercise, including PV recruitment and distension and LA reservoir expansion, and the integrated mean pressure-flow response to exercise in healthy adults. Finally, it will consider how these responses may be impacted by age-related remodeling and modified by sex-related cardiopulmonary differences. Studying the determinants and recognizing the normal limits of PV pressure-flow relations during exercise will improve our understanding of cardiopulmonary mechanisms that facilitate or limit exercise.

Spiral wave reentry as a mechanism of lethal ventricular arrhythmias has been widely demonstrated in animal experiments and recordings from human hearts. It has been shown that in structurally normal hearts spiral waves are unstable, breaking up into multiple wavelets via dynamical instabilities. However, many of the second-generation action potential models give rise only to stable spiral waves, raising issues regarding the underlying mechanisms of spiral wave breakup. In this study, we carried out computer simulations of two-dimensional homogeneous tissues using five ventricular action potential models. We show that the transient outward potassium current ( I to ), although it is not required, plays a key role in promoting spiral wave breakup in all five models. As the maximum conductance of I to increases, it first promotes spiral wave breakup and then stabilizes the spiral waves. In the absence of I to , speeding up the L-type calcium kinetics can prevent spiral wave breakup, however, with the same speedup kinetics, spiral wave breakup can be promoted by increasing I to . Increasing I to promotes single-cell dynamical instabilities, including action potential duration alternans and chaos, and increasing I to further suppresses these action potential dynamics. These cellular properties agree with the observation that increasing I to first promotes spiral wave breakup and then stabilizes spiral waves in tissue. Implications of our observations to spiral wave dynamics in the real hearts and action potential model improvements are discussed. NEW & NOTEWORTHY Spiral wave breakup manifesting as multiple wavelets is a mechanism of ventricular fibrillation. It has been known that spiral wave breakup in cardiac tissue can be caused by a steeply sloped action potential duration restitution curve, a property mainly determined by the recovery of L-type calcium current. Here, we show that the transient outward potassium current ( I to ) is another current that plays a key role in spiral wave breakup, that is, spiral waves can be stable for low and high maximum I to conductance but breakup occurs for intermediate maximum I to conductance. Since I to is present in normal hearts of many species and required for Brugada syndrome, it may play an important role in the spiral wave stability and arrhythmogenesis under both normal condition and Brugada syndrome.

The creatine kinase system facilitates energy transfer between mitochondria and the major ATPases in the heart. Creatine-deficient mice, which lack arginine-glycine amidinotransferase (AGAT) to synthesize creatine and homoarginine, exhibit reduced cardiac contractility. We studied how the absence of a functional CK system influences calcium handling in isolated cardiomyocytes from AGAT-knockouts and wild-type littermates as well as in AGAT-knockout mice receiving lifelong creatine supplementation via the food. Using a combination of whole cell patch clamp and fluorescence microscopy, we demonstrate that the L-type calcium channel (LTCC) current amplitude and voltage range of activation were significantly lower in AGAT-knockout compared with wild-type littermates. Additionally, the inactivation of LTCC and the calcium transient decay were significantly slower. According to our modeling results, these changes can be reproduced by reducing three parameters in knockout mice when compared with wild-type: LTCC conductance, the exchange constant of Ca 2+ transfer between subspace and cytosol, and SERCA activity. Because tissue expression of LTCC and SERCA protein were not significantly different between genotypes, this suggests the involvement of posttranslational regulatory mechanisms or structural reorganization. The AGAT-knockout phenotype of calcium handling was fully reversed by dietary creatine supplementation throughout life. Our results indicate reduced calcium cycling in cardiomyocytes from AGAT-knockouts and suggest that the creatine kinase system is important for the development of calcium handling in the heart. NEW & NOTEWORTHY Creatine-deficient mice lacking arginine-glycine amidinotransferase exhibit compromised cardiac function. Here, we show that this is at least partially due to an overall slowing of calcium dynamics. Calcium influx into the cytosol via the L-type calcium current (LTCC) is diminished, and the rate of the sarcoendoplasmic reticulum calcium ATPase (SERCA) pumping calcium back into the sarcoplasmic reticulum is slower. The expression of LTCC and SERCA did not change, suggesting that the changes are regulatory.

Under normal conditions, coronary blood flow (CBF) provides critical blood supply to the myocardium so that it can appropriately meet the metabolic demands of the body. Dogmatically, there exist several known regulators and modulators of CBF that include local metabolites and neurohormonal factors that can influence the function of the coronary circulation. In disease states such as diabetes and myocardial ischemia, these regulators are impaired or shifted such that CBF is reduced. Although functional considerations have been and continued to be well studied, more recent evidence builds upon established studies that collectively suggest that the relative roles of coronary structure, biomechanics, and the influence of cardiac biomechanics via extravascular compression may also play a significant role in dictating CBF. In this mini review, we discuss these regulators of CBF under normal and pathophysiological conditions and their potential influence on the control of CBF.

Central blood pressure (cBP) is a highly prognostic cardiovascular (CV) risk factor whose accurate, invasive assessment is costly and carries risks to patients. We developed and assessed novel algorithms for estimating cBP from noninvasive aortic hemodynamic data and a peripheral blood pressure measurement. These algorithms were created using three blood flow models: the two- and three-element Windkessel (0-D) models and a one-dimensional (1-D) model of the thoracic aorta. We tested new and existing methods for estimating CV parameters (left ventricular ejection time, outflow BP, arterial resistance and compliance, pulse wave velocity, and characteristic impedance) required for the cBP algorithms, using virtual (simulated) subjects ( n = 19,646) for which reference CV parameters were known exactly. We then tested the cBP algorithms using virtual subjects ( n = 4,064), for which reference cBP were available free of measurement error, and clinical datasets containing invasive ( n = 10) and noninvasive ( n = 171) reference cBP waves across a wide range of CV conditions. The 1-D algorithm outperformed the 0-D algorithms when the aortic vascular geometry was available, achieving central systolic blood pressure (cSBP) errors ≤ 2.1 ± 9.7 mmHg and root-mean-square errors (RMSEs) ≤ 6.4 ± 2.8 mmHg against invasive reference cBP waves ( n = 10). When the aortic geometry was unavailable, the three-element 0-D algorithm achieved cSBP errors ≤ 6.0 ± 4.7 mmHg and RMSEs ≤ 5.9 ± 2.4 mmHg against noninvasive reference cBP waves ( n = 171), outperforming the two-element 0-D algorithm. All CV parameters were estimated with mean percentage errors ≤ 8.2%, except for the aortic characteristic impedance (≤13.4%), which affected the three-element 0-D algorithm's performance. The freely available algorithms developed in this work enable fast and accurate calculation of the cBP wave and CV parameters in datasets containing noninvasive ultrasound or magnetic resonance imaging data. NEW & NOTEWORTHY First, our proposed methods for CV parameter estimation and a comprehensive set of methods from the literature were tested using in silico and clinical datasets. Second, optimized algorithms for estimating cBP from aortic flow were developed and tested for a wide range of cBP morphologies, including catheter cBP data. Third, a dataset of simulated cBP waves was created using a three-element Windkessel model. Fourth, the Windkessel model dataset and optimized algorithms are freely available.

Load, chamber stiffness, and relaxation are the three established determinants of global diastolic function (DF). Coupling of systolic stiffness and isovolumic relaxation has been hypothesized; however, diastolic stiffness-relaxation coupling (DSRC) remains unknown. The parametrized diastolic filling (PDF) formalism, a validated DF model incorporates DSRC. PDF model-predicted DSRC was validated by analysis of 159 Doppler E-waves from a published data set (22 healthy volunteers undergoing bicycle exercise). E-waves at varying (46-120 bpm) heart rates (HR) demonstrated variation in acceleration time (AT), deceleration time (DT), and E-wave peak velocity. AT, DT, and E peak were converted into PDF parameters: stiffness ([Formula: see text]), relaxation ([Formula: see text]), and load ( x o ) using published numerical methods. Univariate linear regression showed that over a twofold increase in HR, AT, and DT decrease ([Formula: see text] = -0.44; P < 0.001 and r = -0.42; P < 0.001, respectively), while, DT/AT remains constant ( r = -0.04; P = 0.67). Similarly, [Formula: see text] increases with HR ( r = 0.55; P < 0.001), while [Formula: see text] has no significant correlation with HR ( r = 0.08; P = 0.32). However, the dimensionless DSRC parameter ψ = c 2 /4 k shows no significant correlation with HR (r = -0.03; P = 0.7). Furthermore, ψ is uniquely determined by DT/AT rather than AT or DT independently. Constancy of ψ in spite of a twofold increase in HR establishes that stiffness ( k ) and relaxation ( c ) are coupled and manifest via a HR-invariant parameter of E-wave asymmetry and should not be considered independent of each other. The manifestation of DSRC through E-wave asymmetry via ψ underscores the value of DT/AT as a physiological, mechanism-derived index of DF. NEW & NOTEWORTHY: Although diastolic stiffness and relaxation are considered independent chamber properties, the cardio-hemic inertial oscillation that generates E-waves obeys Newton's law. E-waves vary with heart rate requiring simultaneous change in stiffness and relaxation. By retrospective analysis of human heart-rate varying transmitral Doppler-data, we show that diastolic stiffness and relaxation are coupled and that the coupling manifests through E-wave asymmetry, quantified through a parametrized diastolic filling model-derived dimensionless parameter, which only depends on deceleration time and acceleration time, readily obtainable via standard echocardiography.

Vasomotor responses conduct among resistance arteries to coordinate blood flow delivery pursuant to energetic demand. Conduction is set by the electrical and mechanical properties of vascular cells, the former tied to how gap junctions and ion channels distribute and dissipate charge, respectively. These membrane proteins are subject to modulation; thus, conduction could be viewed as "pliant" to the current regulatory state. This study used in silico approaches to conceptualize electrical pliancy and to illustrate how gap junctional and ion channel properties distinctly impact conduction along a single skeletal muscle artery or a branching cerebrovascular network. Initial simulations revealed how vascular cells encoded with electrotonic properties best reproduced spreading behavior; the endothelium's importance as a charge source and a longitudinal conduit was readily observed. Alterations in gap junctional conductance produced unique electrical fingerprints: 1 ) decreased endothelial coupling impaired longitudinal but enhanced radial spread, and 2 ) reduced myoendothelial coupling limited radial but enhanced longitudinal spread. Subsequent simulations illustrated how tuning ion channel activity, e.g., inward rectifying- and voltage-gated K + channels, modified charge dissipation, resting membrane potential, and the spread of the electrical phenomenon. Restricting ion channel tuning to a network subregion then revealed how electrical spread could be locally shaped in accordance with the aggregate changes in membrane resistance. In summary, our analysis frames and reimagines electrical conduction as a pliable process, with subtle regulatory changes to membrane proteins shaping network spread and tissue perfusion. NEW & NOTEWORTHY Conducted vasomotor responses depend on initiation and spread of electrical phenomena along arterial walls and their translation into contractile responses. Using computational approaches, we show how subtle but widespread regulation of gap junctions and ion channels can modulate the range and amplitude of electrical spread. Ion channels are regulated by endocrine and mechanical signals and may differ regionally in networks. Subregional electrical changes are not spatially confined but may affect electrical conduction in neighboring regions.

In this study, we mathematically predict retinal vascular resistance (RVR) and retinal blood flow (RBF), we test predictions using laser speckle flowgraphy (LSFG), we estimate the range of vascular autoregulation, and we examine the relationship of RBF with the retinal nerve fiber layer (RNFL) and ganglion cell complex (GCC). Fundus, optical coherence tomography (OCT), and OCT-angiography images, systolic/diastolic blood pressure (SBP/DBP), and intraocular pressure (IOP) measurements were obtained from 36 human subjects. We modeled two circulation markers (RVR and RBF) and estimated individualized lower/higher autoregulation limits (LARL/HARL), using retinal vessel calibers, fractal dimension, perfusion pressure, and population-based hematocrit values. Quantitative LSFG waveforms were extracted from vessels of the same eyes, before and during IOP elevation. LSFG metrics explained most variance in RVR ( R 2  = 0.77/ P = 6.9·10 -9 ) and RBF ( R 2  = 0.65/ P = 1.0·10 -6 ), suggesting that the markers strongly reflect blood flow physiology. Higher RBF was associated with thicker RNFL ( P = 4.0·10 -4 ) and GCC ( P = 0.003), thus also verifying agreement with structural measurements. LARL was at SBP/DBP of 105/65 mmHg for the average subject without arterial hypertension and at 115/75 mmHg for the average hypertensive subject. Moreover, during IOP elevation, changes in RBF were more pronounced than changes in RVR. These observations physiologically imply that healthy subjects are already close to LARL, thus prone to hypoperfusion. In conclusion, we modeled two clinical markers and described a novel method to predict individualized autoregulation limits. These findings could improve understanding of retinal perfusion and pave the way for personalized intervention decisions, when treating patients with coexisting ophthalmic and cardiovascular pathologies. NEW & NOTEWORTHY We describe and test a new approach to quantify retinal blood flow, based on standard clinical examinations and imaging techniques, linked together with a physiological model. We use these findings to generate individualized estimates of the autoregulation range. We provide evidence that healthy subjects are closer to the lower autoregulation limit than thought before. This suggests that some retinas are less prepared to withstand hypoperfusion, even after small intraocular pressure rises or blood pressure drops.

Transcatheter aortic valve replacement (TAVR) is increasingly used to treat severe aortic stenosis (AS) patients. However, little is known regarding the direct effect of TAVR on the ventricular-aortic interaction. In the present study, we aimed to investigate changes in central hemodynamics after successful TAVR. We retrospectively examined 33 cases of severe AS patients (84 ± 6 yr) who underwent TAVR. Invasive measurements of left ventricular and aortic pressures as well as echocardiographic aortic flow were acquired before and after TAVR (maximum within 5 days). We examined alterations in key features of central pressure and flow waveforms, including the aortic augmentation index (AIx), and performed wave separation analysis. Arterial parameters were determined via parameter-fitting on a two-element Windkessel model. Resolution of AS resulted in direct increase in the aortic systolic pressure and maximal aortic flow (131 ± 22 vs. 157 ± 25 mmHg and 237 ± 49 vs. 302 ± 69 mL/s, P < 0.001 for all), whereas the ejection duration decreased ( P < 0.001). We noted a significant decrease in the AIx (from 42 ± 12 to 19 ± 11%, P < 0.001). Of note, the arterial properties remained unchanged. There was a comparable increase in both forward (61 ± 20 vs. 77 ± 20 mmHg, P < 0.001) and backward (35 ± 14 vs. 42 ± 10 mmHg, P = 0.013) pressure wave amplitudes, while their ratio, i.e., the reflection coefficient, was preserved. Our results highlight the impact of TAVR on the ventricular-aortic interaction by affecting the amplitude, shape, and related attributes of the aortic pressure and flow pulse and challenge the interpretation of AIx as a solely vascular measure in AS patients. NEW & NOTEWORTHY Transcatheter aortic valve replacement (TAVR) is linked with an immediate increase in aortic systolic blood pressure and maximal flow, as well as steeper aortic pressure and flow wave upstrokes. After TAVR, the forward wave pumped by the heart is enhanced. Although the arterial properties remain unchanged, the central augmentation index (AIx) is markedly decreased after TAVR. This challenges the interpretation of AIx as a solely vascular measure in patients with aortic valve stenosis.

In patients with heart failure, atrial septal defect (ASD) closure has a risk of inducing life-threatening acute pulmonary edema. The objective of this study was to develop a novel framework for quantitative prediction of hemodynamics after ASD closure. The generalized circulatory equilibrium comprises right and left cardiac output (CO) curves and pulmonary and systemic venous return surfaces. We incorporated ASD into the framework of circulatory equilibrium by representing ASD shunt flow (Q ASD ) by the difference between pulmonary flow (Q P ) and systemic flow (Q S ). To examine the accuracy of prediction, we created ASD in six dogs. Four weeks after ASD creation, we measured left atrial pressure (P LA ), right atrial pressure (P RA ), Q P , and Qs before and after ASD balloon occlusion. We then predicted postocclusion hemodynamics from measured preocclusion hemodynamics. Finally, we numerically simulated hemodynamics under various ASD diameters while changing left and right ventricular function. Predicted postocclusion P LA , P RA , and Q S from preocclusion hemodynamics matched well with those measured [P LA : coefficient of determination ( r 2 ) = 0.96, standard error of estimate (SEE) = 0.89 mmHg, P RA : r 2  = 0.98, SEE = 0.26 mmHg, Q S : r 2  = 0.97, SEE = 5.6 mL·min -1 ·kg -1 ]. A simulation study demonstrated that ASD closure increases the risk of pulmonary edema in patients with impaired left ventricular function and normal right ventricular function, indicating the importance of evaluation for the balance between right and left ventricular function. ASD shunt incorporated into the generalized circulatory equilibrium accurately predicted hemodynamics after ASD closure, which would facilitate safety management of ASD closure. NEW & NOTEWORTHY We developed a framework to predict the impact of atrial septal defect (ASD) closure on hemodynamics by incorporating ASD shunt flow into the framework of circulatory equilibrium. The proposed framework accurately predicted hemodynamics after ASD closure. Patient-specific prediction of hemodynamics may be useful for safety management of ASD closure.

Diastolic dysfunction (DD) is a major component of heart failure with preserved ejection fraction (HFpEF). Accordingly, a profound understanding of the underlying biomechanical mechanisms involved in DD is needed to elucidate all aspects of HFpEF. In this study, we have developed a computational model of DD by leveraging the power of an advanced one-dimensional arterial network coupled to a four-chambered zero-dimensional cardiac model. The two main pathologies investigated were linked to the active relaxation of the myocardium and the passive stiffness of the left ventricular wall. These pathologies were quantified through two parameters for the biphasic delay of active relaxation, which simulate the early and late-phase relaxation delay, and one parameter for passive stiffness, which simulates the increased nonlinear stiffness of the ventricular wall. A parameter sensitivity analysis was conducted on each of the three parameters to investigate their effect in isolation. The three parameters were then concurrently adjusted to produce the three main phenotypes of DD. It was found that the impaired relaxation phenotype can be replicated by mainly manipulating the active relaxation, the pseudo-normal phenotype was replicated by manipulating both the active relaxation and passive stiffness, and, finally, the restricted phenotype was replicated by mainly changing the passive stiffness. This article presents a simple model producing a holistic and comprehensive replication of the main DD phenotypes and presents novel biomechanical insights on how key parameters defining the relaxation and stiffness properties of the myocardium affect the development and manifestation of DD. NEW & NOTEWORTHY This study uses a complete and validated computational model of the cardiovascular system to simulate the two main pathologies involved in diastolic dysfunction (DD), i.e., abnormal active relaxation and increased ventricular diastolic stiffness. The three phenotypes of DD were successfully replicated according to literature data. We elucidate the biomechanical effect of the relaxation pathologies involved and how these pathologies interact to create the various phenotypes of DD.

Analyses of individual atherosclerotic plaques are mostly descriptive, relying, for example, on histological classification by spectral analysis of ultrasound waves or staining and observing particular cellular components. Such passive methods have proved useful for characterizing the structure and vulnerability of plaques but have little quantitative predictive power. Our aim is to introduce and discuss a computational framework to provide insight to clinicians and help them visualize internal plaque dynamics. We use partial differential equations (PDEs) with macrophages, necrotic cells, oxidized lipids, oxygen concentration, and platelet-derived growth factor (PDGF) as primary variables coupled to a biomechanical model to describe vessel growth. The model is deterministic, providing mechanical, morphological, and histological characteristics of an atherosclerotic vessel at any desired future time point. We use our model to create computer-generated animations of a plaque evolution that are in qualitative agreement with published serial ultrasound images and hypothesize possible atherogenic mechanisms. A systems biology model consisting of five differential equations is able to capture the morphology of necrotic cores residing within vulnerable atherosclerotic plaque. In the context of the model, the distribution of oxidized low-density lipoprotein (Ox-LDL) particles, endothelial inflammation, plaque oxygenation (via the presence of vasa vasora), and intimal oxygenation are four important factors that drive changes in core morphology. NEW & NOTEWORTHY In this article, we propose a quantitative framework to describe the evolution of atherosclerotic plaque. We use partial differential equations (PDEs) with macrophages, necrotic cells, oxidized lipids, oxygen concentration, and PDGF as primary variables coupled to a biomechanical model to describe vessel growth. A feature of our method is that it outputs color-coded vessel sections corresponding to regions of the plaque that are necrotic and fibrous, qualitatively similar to images generated by enhanced intravascular ultrasound.

Cardiac electrophysiology and mechanics are strongly interconnected. Calcium is crucial in this complex interplay through its role in cellular electrophysiology and sarcomere contraction. We aim to differentiate the effects of acute β-adrenergic stimulation (β-ARS) and cardiomyocyte stretch (increased sarcomere length) on calcium-transient dynamics and force generation, using a novel computational model of cardiac electromechanics. We implemented a bidirectional coupling between the O'Hara-Rudy model of human ventricular electrophysiology and the MechChem model of sarcomere mechanics through the buffering of calcium by troponin. The coupled model was validated using experimental data from large mammals or human samples. Calcium transient and force were simulated for various degrees of β-ARS and initial sarcomere lengths. The model reproduced force-frequency, quick-release, and isotonic contraction experiments, validating the bidirectional electromechanical interactions. An increase in β-ARS increased the amplitudes of force (augmented inotropy) and calcium transient, and shortened both force and calcium-transient duration (lusitropy). An increase in sarcomere length increased force amplitude even more, but decreased calcium-transient amplitude and increased both force and calcium-transient duration. Finally, a gradient in relaxation along the thin filament may explain the nonmonotonic decay in cytosolic calcium observed with high tension. Using a novel coupled human electromechanical model, we identified differential effects of β-ARS and stretch on calcium and force. Stretch mostly contributed to increased force amplitude and β-ARS to the reduction of calcium and force duration. We showed that their combination, rather than individual contributions, is key to ensure force generation, rapid relaxation, and low diastolic calcium levels. NEW & NOTEWORTHY This work identifies the contribution of electrical and mechanical alterations to regulation of calcium and force under exercise-like conditions using a novel human electromechanical model integrating ventricular electrophysiology and sarcomere mechanics. By better understanding their individual and combined effects, this can uncover arrhythmogenic mechanisms in exercise-like situations. This publicly available model is a crucial step toward understanding the complex interplay between cardiac electrophysiology and mechanics to improve arrhythmia risk prediction and treatment.

Proper inlet boundary conditions are essential for accurate computational fluid dynamics (CFD) modeling. We developed methodology to derive noninvasive FFR B using CFD and computed tomography coronary angiography (CTCA) images. This study aims to assess the influence of brachial mean blood pressure (MBP) and total coronary inflow on FFR B computation. Twenty-two patients underwent both CTCA and FFR measurements. Total coronary flow was computed from left ventricular mass (LVM) measured from CTCA. A total of 286 CFD simulations were run by varying MBP and LVM at 70, 80, 90, 100, 110, 120, and 130% of the measured values. FFR B increased with incrementally higher input values of MBP: 0.78 ± 0.12, 0.80 ± 0.11, 0.82 ± 0.10, 0.84 ± 0.09, 0.85 ± 0.08, 0.86 ± 0.08, and 0.87 ± 0.07, respectively. Conversely, FFR B decreased with incrementally higher inputs value of LVM: 0.86 ± 0.08, 0.85 ± 0.08, 0.84 ± 0.09, 0.84 ± 0.09, 0.83 ± 0.10, 0.83 ± 0.10, and 0.82 ± 0.10, respectively. Noninvasive FFR B calculated using measured MBP and LVM on a total of 30 vessels was 0.84 ± 0.09 and correlated well with invasive FFR (0.83 ± 0.09) ( r  = 0.92, P < 0.001). Positive association was observed between FFR B and MBP input values (mmHg) and negative association between FFR B and LVM values (g). Respective slopes were 0.0016 and -0.005, respectively, suggesting potential application of FFR B in a clinical setting. Inaccurate MBP and LVM inputs differing from patient-specific values could result in misclassification of borderline ischemic lesions. NEW & NOTEWORTHY While brachial mean blood pressure (MBP) and left ventricular mass (LVM) measured from CTCA are the two CFD simulation input parameters, their effects on noninvasive fractional flow reserve (FFR B ) have not been systematically investigated. We demonstrate that inaccurate MBP and LVM inputs differing from patient-specific values could result in misclassification of borderline ischemic lesions. This is important in the clinical application of noninvasive FFR in coronary artery disease diagnosis.

Pulmonary artery (PA) morphometry has been extensively explored in adults, with particular focus on intra-acinar arteries. However, scaling law relationships for length and diameter of extensive preacinar PAs by age have not been previously reported for in vivo human data. To understand preacinar PA growth spanning children to adults, we performed morphometric analyses of all PAs visible in the computed tomography (CT) and magnetic resonance (MR) images from a healthy subject cohort [ n = 16; age: 1-51 yr; body surface area (BSA): 0.49-2.01 m 2 ]. Subject-specific anatomic PA models were constructed from CT and MR images, and morphometric information-diameter, length, tortuosity, bifurcation angle, and connectivity-was extracted and sorted into diameter-defined Strahler orders. Validation of Murray's law, describing optimal scaling exponents of radii for branching vessels, was performed to determine how closely PAs conform to this classical relationship. Using regression analyses of vessel diameters and lengths against orders and patient metrics (BSA, age, height), we found that diameters increased exponentially with order and allometrically with patient metrics. Length increased allometrically with patient metrics, albeit weakly. The average tortuosity index of all vessels was 0.026 ± 0.024, average bifurcation angle was 28.2 ± 15.1°, and average Murray's law exponent was 2.92 ± 1.07. We report a set of scaling laws for vessel diameter and length, along with other morphometric information. These provide an initial understanding of healthy structural preacinar PA development with age, which can be used for computational modeling studies and comparison with diseased PA anatomy. NEW & NOTEWORTHY Pulmonary artery (PA) morphometry studies to date have focused primarily on large arteries and intra-acinar arteries in either adults or children, neglecting preacinar arteries in both populations. Our study is the first to quantify in vivo preacinar PA morphometry changes spanning infants to adults. For preacinar arteries > 1 mm in diameter, we identify scaling laws for vessel diameters and lengths with patient metrics of growth and establish a healthy PA morphometry baseline for most preacinar PAs.

Current thinking suggests that wave reflection in arteries limits pulse pressure and hydraulic energy (HE) transmission to the microvasculature and that this protective effect reduces with advancing age. However, according to transmission line theory, pressure transmission (T p ) and reflection (R) coefficients are proportional (T p  = 1 + R), implying that wave reflection would promote rather than limit pressure transmission. We hypothesized that increasing distal pulse pressure (PPd) with age is instead related to increased proximal pulse pressure (PPp) and its forward component and that these are modulated by arterial compliance. A one-dimensional model of a fractal arterial tree containing 21 generations was constructed. Wave speed in each vessel was prescribed to achieve a uniform R at every junction, with changes in R achieved by progressively stiffening proximal or distal vessels. For both stiffening scenarios, decreasing reflection led to a decrease or no change in PPd when forward pressure or compliance were held constant, respectively, suggesting that wave reflection per se does not limit pressure transmission. Proximal pulse pressure, its forward component, and PPd increased with decreasing compliance; furthermore, proximal and distal pulse pressures were approximately proportional. With fixed compliance but decreasing reflection, HE transmission increased, whereas pressure transmission decreased, consistent with transmission line theory. In conclusion, wave reflection does not protect the microvasculature from high PPd; rather, PPp and PPd are modulated by arterial compliance, which reduces with age. Wave reflection has opposing effects on pressure and HE transmission; hence, the relative importance of pressure versus HE in contributing to microvascular damage warrants investigation. NEW & NOTEWORTHY With aging, a reduction in the stiffness gradient between elastic and muscular arteries is thought to reduce wave reflection in conduit arteries, leading to increased pulsatile pressure transmission into the microvasculature. This assumes that wave reflection limits pressure transmission in arteries. However, using a computational model, we showed that wave reflection promotes pulsatile pressure transmission, although it does limit hydraulic energy transmission. Increased microvascular pulse pressure with aging is instead related to decreasing arterial compliance.

Endothelial dysfunction is thought to underpin atherosclerotic cardiovascular disease. The most widely used in vivo test of endothelial function is flow-mediated dilation (FMD). However, the results of FMD may be subject to some confounding factors that are not fully understood. We investigated potential biophysical confounding factors that could cause a disassociation between FMD and true endothelial cell shear stress response (the release of endothelium-dependent relaxing factors in response to wall shear stress). Arterial hemodynamics during FMD was simulated using a novel computational modeling approach. The model included an endothelial response function relating changes in wall shear stress to changes in local vascular stiffness in the arm arteries and accounted for vascular stiffening with increasing blood pressure. The hemodynamic effects of cuff inflation and deflation were modeled by prescribing intraluminal arterial pressure changes and peripheral vasodilation. Evolution of arterial diameter and flow velocity during FMD was assessed by comparison against in vivo data. Our model revealed that vasoconstriction occurring immediately after cuff deflation is independent of endothelial response function and entirely caused by the change in transmural pressure along conduit arteries. Moreover, for the same endothelial response function model, FMD values increased exponentially with increasing peak flow velocity, decreased linearly with increasing arterial stiffness at a rate of 0.95%/MPa, and increased linearly with increasing central blood pressure at a rate of 0.22%/mmHg. Dependence of FMD on confounding factors, such as arterial stiffness and blood pressure, suggests that the current FMD test may not reflect the true endothelial cell response. NEW & NOTEWORTHY First, a novel computational model simulating arterial hemodynamics during flow-mediated dilation (FMD) was proposed. Second, the model was used to explain why FMD may be influenced by endothelium-independent factors, showing that FMD results are 1 ) partly masked by the vasoconstriction due to the change in transmural pressure and 2 ) affected by physiological factors (i.e., arterial stiffness and arterial blood pressure) that are difficult to eliminate due to their multiple interactions.

The Fontan circulation is characterized as a nonpulsatile flow propagation without a pressure-generating ventricle. However, flow through the Fontan circulation still exhibits oscillatory waves as a result of pressure changes generated by the systemic single ventricle. Identification of discrete flow patterns through the Fontan circuit may be important to understand single ventricle performance. Ninety-seven patients with Fontan circulation underwent phase-contrast MRI of the right pulmonary artery, yielding subject-specific flow waveforms. Principal component (PC) analysis was performed on preprocessed flow waveforms. Principal components were then correlated with standard MRI indices of function, volume, and aortopulmonary collateral flow. The first principal component (PC) described systolic versus diastolic-dominant flow through the Fontan circulation, accounting for 31.3% of the variance in all waveforms. The first PC correlated with end-diastolic volume ( R  = 0.34, P = 0.001), and end-systolic volume ( R  = 0.30, P = 0.003), cardiac index ( R  = 0.51, P < 0.001), and the amount of aortopulmonary collateral flow ( R  = 0.25, P = 0.027)-lower ventricular volumes and a smaller volume of collateral flow-were associated with diastolic-dominant cavopulmonary flow. The second PC accounted for 19.5% of variance and described late diastolic acceleration versus deceleration and correlated with ejection fraction-diastolic deceleration was associated with higher ejection fraction. Principal components describing the diastolic flow variations in pulmonary arteries are related to the single ventricle function and volumes. Particularly, diastolic-dominant flow without late acceleration appears to be related to preserved ventricular volume and function, respectively. NEW & NOTEWORTHY The exact physiological significance of flow oscillations of phasic and temporal flow variations in Fontan circulation is unknown. With the use of principal component analysis, we discovered that flow variations in the right pulmonary artery of Fontan patients are related to the single ventricle function and volumes. Particularly, diastolic-dominant flow without late acceleration appears to be related to more ideal ventricular volume and systolic function, respectively.

The spatiotemporal features of normal in vivo cardiac motion are well established. Longitudinal velocity has become a focus of diastolic function (DF) characterization, particularly the tissue Doppler e' -wave, manifesting in early diastole when the left ventricle (LV) is a mechanical suction pump (dP/dV < 0). To characterize DF and elucidate mechanistic features, several models have been proposed and have been previously compared algebraically, numerically, and in their ability to fit physiological velocity data. We analyze two previously noncompared models of early rapid-filling lengthening velocity (Doppler e' -wave): parametrized diastolic filling (PDF) and force balance model (FBM). Our initial numerical experiments sampled FBM-generated e' ( t ) contours as input to determine PDF model predicted fit. The resulting exact numerical agreement [standard error of regression (SER) = 9.06 × 10 -16 ] was not anticipated. Therefore, we analyzed all published FBM-generated e' ( t ) contours and observed identical agreement. We re-expressed FBM's algebraic expressions for e' ( t ) and observed for the first time that model-based predictions for lengthening velocity by the FBM and the PDF model are mathematically identical: e' ( t ) = γe -α t sinh(β t ), thereby providing exact algebraic relations between the three PDF parameters and the six FBM parameters. Previous pioneering experiments have independently established the unique determinants of e'( t ) to be LV relaxation, restoring forces (stiffness), and load. In light of the exact intermodel agreement, we conclude that the three PDF parameters, relaxation, stiffness (restoring forces), and load, are unique determinants of DF and e' ( t ). Thus, we show that only the PDF formalism can compute the three unique, independent, physiological determinants of long-axis LV myocardial velocity from e' ( t ). NEW & NOTEWORTHY We show that two separate, independently derived physiological (kinematic) models predict mathematically identical expressions for LV-lengthening velocity (Doppler e' -wave), indicating that damped harmonic oscillatory motion is a physiologically accurate model of diastolic function. Although both models predict the same "overdamped" velocity contour, only one model solves the "inverse problem" and generates unique, lumped parameters of relaxation, stiffness (restoring force), and load from the e' -wave.

The oxygen consumption by the heart and its extraction from the coronary arterial blood are the highest among all organs. Any increase in oxygen demand due to a change in heart metabolic activity requires an increase in coronary blood flow. This functional requirement of adjustment of coronary blood flow is mediated by coronary flow regulation to meet the oxygen demand without any discomfort, even under strenuous exercise conditions. The goal of this article is to provide an overview of the theoretical and computational models of coronary flow regulation and to reveal insights into the functioning of a complex physiological system that affects the perfusion requirements of the myocardium. Models for three major control mechanisms of myogenic, flow, and metabolic control are presented. These explain how the flow regulation mechanisms operating over multiple spatial scales from the precapillaries to the large coronary arteries yield the myocardial perfusion characteristics of flow reserve, autoregulation, flow dispersion, and self-similarity. The review not only introduces concepts of coronary blood flow regulation but also presents state-of-the-art advances and their potential to impact the assessment of coronary microvascular dysfunction (CMD), cardiac-coronary coupling in metabolic diseases, and therapies for angina and heart failure. Experimentalists and modelers not trained in these models will have exposure through this review such that the nonintuitive and highly nonlinear behavior of coronary physiology can be understood from a different perspective. This survey highlights knowledge gaps, key challenges, future research directions, and novel paradigms in the modeling of coronary flow regulation.

Various experimental mouse models are extensively used to research human diseases, including atrial fibrillation, the most common cardiac rhythm disorder. Despite this, there are no comprehensive mathematical models that describe the complex behavior of the action potential and [Ca 2+ ] i transients in mouse atrial myocytes. Here, we develop a novel compartmentalized mathematical model of mouse atrial myocytes that combines the action potential, [Ca 2+ ] i dynamics, and β-adrenergic signaling cascade for a subpopulation of right atrial myocytes with developed transverse-axial tubule system. The model consists of three compartments related to β-adrenergic signaling (caveolae, extracaveolae, and cytosol) and employs local control of Ca 2+ release. It also simulates ionic mechanisms of action potential generation and describes atrial-specific Ca 2+ handling as well as frequency dependences of the action potential and [Ca 2+ ] i transients. The model showed that the T-type Ca 2+ current significantly affects the later stage of the action potential, with little effect on [Ca 2+ ] i transients. The block of the small-conductance Ca 2+ -activated K + current leads to a prolongation of the action potential at high intracellular Ca 2+ . Simulation results obtained from the atrial model cells were compared with those from ventricular myocytes. The developed model represents a useful tool to study complex electrical properties in the mouse atria and could be applied to enhance the understanding of atrial physiology and arrhythmogenesis. NEW & NOTEWORTHY A new compartmentalized mathematical model of mouse right atrial myocytes was developed. The model simulated action potential and Ca 2+ dynamics at baseline and after stimulation of the β-adrenergic signaling system. Simulations showed that the T-type Ca 2+ current markedly prolonged the later stage of atrial action potential repolarization, with a minor effect on [Ca 2+ ] i transients. The small-conductance Ca 2+ -activated K + current block resulted in prolongation of the action potential only at the relatively high intracellular Ca 2+ .

Computational modeling based on experimental data remains an important component in cardiac electrophysiological research, especially because clinical data such as human action potential (AP) dynamics are scarce or limited by practical or ethical concerns. Such modeling has been used to develop and test a variety of mechanistic hypotheses, with the majority of these studies involving the rate dependence of AP duration (APD) including APD restitution and conduction velocity (CV). However, there is very little information regarding the complex dynamics at the boundary of repolarization (or refractoriness) and reexcitability. Here, we developed a "minimal" ionic model of the human AP, based on in vivo human monophasic AP (MAP) recordings obtained during clinical programmed electrical stimulation (PES) to address the progressive decrease in AP take-off potential (TOP) and associated CV slowing seen during three tightly spaced extrastimuli. Recent voltage-clamp data demonstrating the effect of intracellular calcium on sodium current availability were incorporated and were required to reproduce large (>15 mV) elevations in take-off potential and progressive encroachment. Introducing clinically observed APD gradients into the model enabled us to replicate the dynamic response to PES in patients leading to conduction block and reentry formation for the positive, but not the negative, APD gradient. Finally, we modeled the dynamics of reentry and show that spiral waves follow a meandering trajectory with a period of ~180 ms. We conclude that our model reproduces a variety of electrophysiological behavior including the response to sequential premature stimuli and provides a basis for studies of the initiation of reentry in human ventricular tissue. NEW & NOTEWORTHY This work presents a new model of the action potential of the human which reproduces the complex dynamics during premature stimulation in patients.

Clinical investigations have demonstrated that mitral regurgitation (MR) quantification using echocardiography (echo) may significantly underestimate or overestimate the regurgitant volume, especially for two-dimensional (2D) echo. Computer modeling and simulated echo were conducted to evaluate the fundamental assumptions in the echo quantification of primary MR that is due to posterior mitral leaflet prolapse. The theoretical flaw of the proximal isovelocity surface area (PISA) method originates from the assumption that the MR flow rate is the product of the isovelocity surface area and aliasing velocity, which is only valid when the velocity vectors are perpendicular to the isovelocity surface. Other factors such as the Doppler angle effect, the view planes of 2D echo, and the single time instant of PISA were also analyzed. We find that the hemielliptic PISA method gives the smallest error for moderate and severe MR cases compared with other PISA methods. Compared with the PISA method, the volumetric technique (VT) is theoretically more robust. By considering correction factors that are caused by nonflat velocity profiles and the closing volume of the aortic valve, the accuracy of the VT method can be significantly improved. The corrected volumetric technique provides more accurate results compared with the PISA methods, especially for mild MR. NEW & NOTEWORTHY We evaluate the accuracy of common echocardiography techniques for the quantification of primary mitral regurgitations using computer modeling. The hemielliptic proximal isovelocity surface area (PISA) method gives the smallest error (within 15%) for moderate and severe mitral regurgitation cases compared with other PISA methods. The volumetric method is theoretically more robust than the PISA method. The accuracy of the volumetric method can be improved by a correction factor around 0.7 because of the nonflat velocity profiles and the closing volume of the aortic valve.

Recognition that coronary blood flow is tightly coupled with myocardial metabolism has been appreciated for well over half a century. However, exactly how coronary microvascular resistance is tightly coupled with myocardial oxygen consumption (MV̇o 2 ) remains one of the most highly contested mysteries of the coronary circulation to this day. Understanding the mechanisms responsible for local metabolic control of coronary blood flow has been confounded by continued debate regarding both anticipated experimental outcomes and data interpretation. For a number of years, coronary venous Po 2 has been generally accepted as a measure of myocardial tissue oxygenation and thus the classically proposed error signal for the generation of vasodilator metabolites in the heart. However, interpretation of changes in coronary venous Po 2 relative to MV̇o 2 are quite nuanced, inherently circular in nature, and subject to confounding influences that remain largely unaccounted for. The purpose of this review is to highlight difficulties in interpreting the complex interrelationship between key coronary outcome variables and the arguments that emerge from prior studies performed during exercise, hemodilution, hypoxemia, and alterations in perfusion pressure. Furthermore, potential paths forward are proposed to help to facilitate further dialogue and study to ultimately unravel what has become the Gordian knot of the coronary circulation.

Myocardial bridging (MB) is linked to angina and myocardial ischemia and may lead to sudden cardiac death in patients with hypertrophic cardiomyopathy (HCM). However, it remains unclear how MB affect the coronary blood flow in HCM patients. The aim of this study was to assess the effects of MB on coronary hemodynamics in HCM patients. Fifteen patients with MB (7 HCM and 8 non-HCM controls) in their left anterior descending (LAD) coronary artery were chosen. Transient computational fluid dynamics (CFD) simulations were conducted in anatomically realistic models of diseased (with MB) and virtually healthy (without MB) LAD from these patients, reconstructed from biplane angiograms. Our CFD simulation results demonstrated that dynamic compression of MB led to diastolic flow disturbances and could significantly reduce the coronary flow in HCM patients as compared with non-HCM group ( P < 0.01). The pressure drop coefficient was remarkably higher ( P < 0.05) in HCM patients. The flow rate change is strongly correlated with both upstream Reynolds number and MB compression ratio, while the MB length has less impact on coronary flow. The hemodynamic results and clinical outcomes revealed that HCM patients with an MB compression ratio higher than 65% required a surgical intervention. In conclusion, the transient MB compression can significantly alter the diastolic flow pattern and wall shear stress distribution in HCM patients. HCM patients with severe MB may need a surgical intervention. NEW & NOTEWORTHY In this study, the hemodynamic significance of myocardial bridging (MB) in patients with hypertrophic cardiomyopathy (HCM) was investigated to provide valuable information for surgical decision-making. Our results illustrated that the transient MB compression led to complex flow patterns, which can significantly alter the diastolic flow and wall shear stress distribution. The hemodynamic results and clinical outcomes demonstrated that patients with HCM and an MB compression ratio higher than 65% required a surgical intervention.

The Fontan procedure for univentricular heart defects creates a nonphysiologic circulation where systemic venous blood drains directly into the pulmonary arteries, leading to multiorgan dysfunction secondary to chronic low-shear nonpulsatile pulmonary blood flow and central venous hypertension. Although blood viscosity increases exponentially in this low-shear environment, the role of shear-dependent ("non-Newtonian") blood viscosity in this pathophysiology is unclear. We studied three-dimensional (3D)-printed Fontan models in an in vitro flow loop with a Philips 3-T magnetic resonance imaging (MRI) scanner. A 4D flow phase-contrast sequence was used to acquire a time-varying 3D velocity field for each experimental condition. On the basis of blood viscosity of a cohort of patients who had undergone the Fontan procedure, it was decided to use 0.04% xanthan gum as a non-Newtonian blood analog; 45% glycerol was used as a Newtonian control fluid. MRI data were analyzed using GTFlow and MATLAB software. The primary outcome, power loss, was significantly higher with the Newtonian fluid [14.8 (13.3, 16.4) vs. 8.1 (6.4, 9.8)%, medians with 95% confidence interval, P < 0.0001]. The Newtonian fluid also demonstrated marginally higher right pulmonary artery flow, marginally lower shear stress, and a trend toward higher caval flow mixing. Outcomes were modulated by Fontan model complexity, cardiac output, and caval flow ratio. Vortexes, helical flow, and stagnant flow were more prevalent with the non-Newtonian fluid. Our data demonstrate that shear-dependent viscosity significantly alters qualitative flow patterns, power loss, pulmonary flow distribution, shear stress, and caval flow mixing in synthetic models of the Fontan circulation. Potential clinical implications include effects on exercise capacity, ventilation-perfusion matching, risk of pulmonary arteriovenous malformations, and risk of thromboembolism. NEW & NOTEWORTHY Although blood viscosity increases exponentially in low-shear environments, the role of shear-dependent ("non-Newtonian") blood viscosity in the pathophysiology of the low-shear Fontan circulation is unclear. We demonstrate that shear-dependent viscosity significantly alters qualitative flow patterns, power loss, pulmonary flow distribution, shear stress, and caval flow mixing in synthetic models of the Fontan circulation. Potential clinical implications include effects on exercise capacity, ventilation-perfusion matching, risk of pulmonary arteriovenous malformations, and risk of thromboembolism.

The deteriorating nature of severe functional tricuspid regurgitation (FTR) has led to the heightened interest in this pathology. However, therapies are heterogeneous and an ideal technique is uncertain. The hemodynamic impact on the cardiac chamber following therapeutic repairs has not been well studied, while its analysis could be used to predict the treatment success. In this study, the hemodynamics of the right ventricle (RV) after 1 ) clover edge-to-edge tricuspid repair, and 2 ) double orifice tricuspid repair was evaluated in three right heart models using an ex vivo pulsatile platform emulating severe FTR with the aid of stereoscopic particle image velocimetry. Although all repairs substantially reduced tricuspid regurgitant area, they resulted in a more than 50% reduction in diastolic tricuspid valve (TV) opening area. Splitting the TV orifice into multiple smaller orifices by both repairs eliminated the ring-shaped vortical structure inside the RV observed in FTR cases. Postrepair RV domain was mostly occupied with irregular vortical features and isolated vortex residuals. Moreover, vortical features varied among repair samples, indicating enhanced sensitivity of RV flow to postrepair TV morphology. Compared with clover repair, double orifice subjected the RV to enhanced swirling motions and exposed more regions to vortical motions, potentially indicating better rinsing and lower risk of mural thrombus formation. Double orifice repair increased the levels of RV mean kinetic energy and viscous energy loss than those observed in clover repair, although the impact of these on the cardiac efficiency remains unclear. These preliminary insights could be used to improve future treatment design and planning. NEW & NOTEWORTHY While clover and double orifice tricuspid repairs markedly improved leaflet coaptation, they substantially reduced diastolic tricuspid opening area. Postrepair right ventricle (RV) exhibited specific hemodynamic traits, including the loss of ring-shaped vortical structure and the enhanced sensitivity of RV flow to postrepair tricuspid valve morphology. Compared with clover technique, double orifice repair led to higher swirling motions in the RV domain, which could indicate lower risk of mural thrombus formation.