Evaluation of vascular aging on measures of cardiac function and mechanical efficiency: insights from in-silico modeling.

Introduction: This study evaluated the hypothesis that vascular aging (VA) reduces ventricular contractile function and mechanical efficiency (ME) using the left ventricular pressure-volume (PV) construct. Methods: A previously published in-silico computational model (CM) was modified to evaluate the hypothesis in two phases. In phase I, the CM included five settings of aortic compliance (CA) from normal to stiff, studied at a heart rate of 80 bpm, and phase II included the normal to stiff CA settings evaluated at 60, 100, and 140 bpm. The PV construct provided steady-state and transient data through a simulated vena caval occlusion (VCO). The steady-state data included left ventricular volumes (EDV and ESV), stroke work (SW), and VCO provided the PV area (PVA) data in addition to the three measures of contractile state (CS): end-systolic pressure-volume relationship (ESPVR), dP/dtmax-EDV and preload recruitable stroke work (PRSW). Finally, ME was calculated with the SW/PVA parameter. Results: In phase I, EDV and ESV increased, as did SW and PVA. The impact on the CS parameters demonstrated a small decrease in ESPVR, no change in dP/dtmax-EDV, and a large increase in PRSW. ME decreased from 71.5 to 60.8%, respectively. In phase II, at the normal and stiff CA settings, across the heart rates studied, EDV and ESV decreased, ESPVR and dP/dtmax-EDV increased and PRSW decreased. ME decreased from 76.4 to 62.6% at the normal CA and 65.8 to 53.2% at the stiff CA. Discussion: The CM generated new insights regarding how the VA process impacts the contractile state of the myocardium and ME.

Evaluation of ventricular-vascular coupling with critical care metrics: An in silico approach.

Mean arterial pressure and cardiac output provide insufficient guidance for the management of intraoperative hypotension (IOH). In silico models offer additional insights into acute changes in hemodynamic parameters that may be encountered during IOH. A computational model (CM) generated parameters quantifying ventricular-vascular coupling, and pressure-volume construct across levels of aortic compliance (CA ). We studied how a loss from normal-to-stiff CA impacts critical care metrics of hemodynamics during vascular occlusion. Pulse pressure (PP), end-systolic pressure (Pes ), arterial compliance (Art-ca), arterial elastance (Art-ea), and dynamic arterial elastance (Eadyn), along mechanical efficiency (ME) were measured at five levels of CA . A loss in CA impacted all variables. During steady-state conditions, PP, Pes , and stroke work increased significantly as CA decreased. Art-ca decreased and Art-ea increased similarly; Eadyn increased and ME decreased. During a decrease in preload across all CA levels, arterial dynamics measures remained linear. The CM demonstrated that a loss in CA impacts measures of arterial dynamics during steady-state and transient conditions and the model demonstrates that critical care metrics are sensitive to changes in CA . While Art-ca and Art-ea were sensitive to changes in preload, Eadyn did not change.

Impact of myocardial viability and left ventricular lead location on clinical outcome in cardiac resynchronization therapy recipients with ischemic cardiomyopathy.

INTRODUCTION: Cardiac resynchronization therapy (CRT) recipients with ischemic cardiomyopathy (ICM) have scar segments that may limit ventricular resynchronization and clinical response. The impact of myocardial viability at the left ventricular (LV) pacing site on CRT response is poorly elucidated. METHODS AND RESULTS: A retrospective cohort of 160 ICM patients with single photon emission computed tomography-myocardial perfusion imaging before device implantation were included. Coronary venous angiography and chest radiographs helped classify segmental location of LV lead (LVL). The primary outcome was a composite of heart failure (HF) hospitalization and mortality at 3 years, and secondary outcome was change in systolic function at 6 months. The patients were divided into groups based on the myocardial substrate at the site of LVL: LVL on or adjacent to (1) normal myocardium (LVL-N, n = 64), (2) segmental scar (LVL-S, n = 62), and (3) scar and ischemia (LVL-SI, n = 34). Upon follow-up, 75 (47%) patients reached primary endpoint with a higher incidence noted in LVL-S (60%), and LVL-SI (53%), compared to 31% in LVL-N (P = 0.004). Kaplan Meier method demonstrated poor event free survival for primary outcome in LVL-S (P = 0.002), and LVL-SI (P = 0.03). In Cox proportional hazard model, LVL-S (HR: 2.26, P = 0.004), and LVL-SI (1.9, P = 0.047) were independent predictors of primary outcome. CONCLUSION: In CRT recipients with ICM, scar and reversible ischemia in or adjacent to LV pacing site were independent predictors of HF hospitalization and death.

Myofiber prestretch magnitude determines regional systolic function during ectopic activation in the tachycardia-induced failing canine heart.

Electrical dyssynchrony leads to prestretch in late-activated regions and alters the sequence of mechanical contraction, although prestretch and its mechanisms are not well defined in the failing heart. We hypothesized that in heart failure, fiber prestretch magnitude increases with the amount of early-activated tissue and results in increased end-systolic strains, possibly due to length-dependent muscle properties. In five failing dog hearts with scars, three-dimensional strains were measured at the anterolateral left ventricle (LV). Prestretch magnitude was varied via ventricular pacing at increasing distances from the measurement site and was found to increase with activation time at various wall depths. At the subepicardium, prestretch magnitude positively correlated with the amount of early-activated tissue. At the subendocardium, local end-systolic strains (fiber shortening, radial wall thickening) increased proportionally to prestretch magnitude, resulting in greater mean strain values in late-activated compared with early-activated tissue. Increased fiber strains at end systole were accompanied by increases in preejection fiber strain, shortening duration, and the onset of fiber relengthening, which were all positively correlated with local activation time. In a dog-specific computational failing heart model, removal of length and velocity dependence on active fiber stress generation, both separately and together, alter the correlations between local electrical activation time and timing of fiber strains but do not primarily account for these relationships.

Left ventricular midwall fibrosis as a predictor of mortality and morbidity after cardiac resynchronization therapy in patients with nonischemic cardiomyopathy.

OBJECTIVES: The aim of this study was to determine whether left ventricular (LV) midwall fibrosis, detected by midwall hyperenhancement (MWHE) on late gadolinium enhancement cardiovascular magnetic resonance (CMR) imaging, predicts mortality and morbidity in patients with dilated cardiomyopathy (DCM) undergoing cardiac resynchronization therapy (CRT). BACKGROUND: Midwall fibrosis predicts mortality and morbidity in patients with DCM. METHODS: Patients with DCM with (+) or without (-) MWHE (n = 20 and n = 77, respectively) as well as 161 patients with ischemic cardiomyopathy (ICM) undergoing CRT (n = 258) were followed up for a maximum of 8.7 years. RESULTS: Among patients with DCM, +MWHE predicted cardiovascular mortality (hazard ratio [HR]: 18.6; 95% confidence intervals [CI]: 3.51 to 98.5; p = 0.0008), total mortality or hospitalization for major adverse cardiovascular events (HR: 7.57; 95% CI: 2.71 to 21.2; p < 0.0001), and cardiovascular mortality or heart failure hospitalizations (HR: 9.56; 95% CI: 2.72 to 33.6; p = 0.0004), independent of New York Heart Association class, QRS duration, atrial fibrillation, LV volumes, LV ejection fraction, and a CMR-derived measure of dyssynchrony. Among patients with DCM and ICM, the risk of cardiovascular mortality for DCM +MWHE (adjusted HR: 18.5; 95% CI: 3.93 to 87.3; p = 0.0002) was similar to that for ICM (adjusted HR: 21.0; 95% CI: 5.06 to 87.2; p < 0.0001). Both DCM +MWHE and ICM were predictors of pump failure death as well as sudden cardiac death. LV reverse remodeling was observed in DCM -MWHE and in ICM but not in DCM +MWHE. CONCLUSIONS: Midwall fibrosis is an independent predictor of mortality and morbidity in patients with DCM undergoing CRT. The outcome of DCM with midwall fibrosis is similar to that of ICM. This relationship is mediated by both pump failure and sudden cardiac death.

Efficient characterization of inhomogeneity in contraction strain pattern.

Cardiac dyssynchrony often accompanies patients with heart failure (HF) and can lead to an increase in mortality rate. Cardiac resynchronization therapy (CRT) has been shown to provide substantial benefits to the HF population with ventricular dyssynchrony; however, there still exists a group of patients who do not respond to this treatment. In order to better understand patient response to CRT, it is necessary to quantitatively characterize both electrical and mechanical dyssynchrony. The quantification of mechanical dyssynchrony via characterization of contraction strain field inhomogeneity is the focus of this modeling investigation. Raw data from a 3D finite element (FE) model were received from Roy Kerckhoffs et al. and analyzed in MATLAB. The FE model consisted of canine left and right ventricles coupled to a closed circulation with the effects of the pericardium acting as a pressure on the epicardial surface. For each of three simulations (normal synchronous, SYNC, right ventricular apical pacing, RVA, and left ventricular free wall pacing, LVFW) the Gauss point locations and values were used to generate lookup tables (LUTs) with each entry representing a location in the heart. In essence, we employed piecewise cubic interpolation to generate a fine point cloud (LUTs) from a course point cloud (Gauss points). Strain was calculated in the fiber direction and was then displayed in multiple ways to better characterize strain inhomogeneity. By plotting average strain and standard deviation over time, the point of maximum contraction and the point of maximal inhomogeneity were found for each simulation. Strain values were organized into seven strain bins to show operative strain ranges and extent of inhomogeneity throughout the heart wall. In order to visualize strain propagation, magnitude, and inhomogeneity over time, we created 2D area maps displaying strain over the entire cardiac cycle. To visualize spatial strain distribution at the time point of maximum inhomogeneity, a 3D point cloud was created for each simulation, and a CURE index was calculated. We found that both the RVA and LFVW simulations took longer to reach maximum contraction than the SYNC simulation, while also exhibiting larger disparities in strain values during contraction. Strain in the hoop direction was also analyzed and was found to be similar to the fiber strain results. It was found that our method of analyzing contraction strain pattern yielded more detailed spacial and temporal information about fiber strain in the heart over the cardiac cycle than the more conventional CURE index method. We also observed that our method of strain binning aids in visualization of the strain fields, and in particular, the separation of the mass points into separate images associated with each strain bin allows the strain pattern to be explicitly compartmentalized.

Improvement in pump function with endocardial biventricular pacing increases with activation time at the left ventricular pacing site in failing canine hearts.

Recently, attention has been focused on comparing left ventricular (LV) endocardial (ENDO) with epicardial (EPI) pacing for cardiac resynchronization therapy. However, the effects of ENDO and EPI lead placement at multiple sites have not been studied in failing hearts. We hypothesized that differences in the improvement of ventricular function due to ENDO vs. EPI pacing in dyssynchronous (DYSS) heart failure may depend on the position of the LV lead in relation to the original activation pattern. In six nonfailing and six failing dogs, electrical DYSS was created by atrioventricular sequential pacing of the right ventricular apex. ENDO was compared with EPI biventricular pacing at five LV sites. In failing hearts, increases in the maximum rate of LV pressure change (dP/dt; r = 0.64), ejection fraction (r = 0.49), and minimum dP/dt (r = 0.51), relative to DYSS, were positively correlated (P < 0.01) with activation time at the LV pacing site during ENDO but not EPI pacing. ENDO pacing at sites with longer activation delays led to greater improvements in hemodynamic parameters and was associated with an overall reduction in electrical DYSS compared with EPI pacing (P < 0.05). These findings were qualitatively similar for nonfailing hearts. Improvement in hemodynamic function increased with activation time at the LV pacing site during ENDO but not EPI pacing. At the anterolateral wall, end-systolic transmural function was greater with local ENDO compared with EPI pacing. ENDO pacing and intrinsic activation delay may have important implications for management of DYSS heart failure.

Ventricular dilation and electrical dyssynchrony synergistically increase regional mechanical nonuniformity but not mechanical dyssynchrony: a computational model.

BACKGROUND: Heart failure (HF) in combination with mechanical dyssynchrony is associated with a high mortality rate. To quantify contractile dysfunction in patients with HF, investigators have proposed several indices of mechanical dyssynchrony, including percentile range of time to peak shortening (WTpeak), circumferential uniformity ratio estimate (CURE), and internal stretch fraction (ISF). The goal of this study was to compare the sensitivity of these indices to 4 major abnormalities responsible for cardiac dysfunction in dyssynchronous HF: dilation, negative inotropy, negative lusitropy, and dyssynchronous activation. METHODS AND RESULTS: All combinations of these 4 major abnormalities were included in 3D computational models of ventricular electromechanics. Compared with a nonfailing heart model, ventricles were dilated, inotropy was reduced, twitch duration was prolonged, and activation sequence was changed from normal to left bundle branch block. In the nonfailing heart, CURE, ISF, and WTpeak were 0.97+/-0.004, 0.010+/-0.002, and 78+/-1 milliseconds, respectively. With dilation alone, CURE decreased 2.0+/-0.07%, ISF increased 58+/-47%, and WTpeak increased 31+/-3%. With dyssynchronous activation alone, CURE decreased 15+/-0.6%, ISF increased 14-fold (+/-3), and WTpeak increased 121+/-4%. With the combination of dilation and dyssynchronous activation, CURE decreased 23+/-0.8%, ISF increased 20-fold (+/-5), and WTpeak increased 147+/-5%. CONCLUSIONS: Dilation and left bundle branch block combined synergistically decreased regional cardiac function. CURE and ISF were sensitive to this combination, but WTpeak was not. CURE and ISF also reflected the relative nonuniform distribution of regional work better than WTpeak. These findings might explain why CURE and ISF are better predictors of reverse remodeling in cardiac resynchronization therapy.

Interlead distance and left ventricular lead electrical delay predict reverse remodeling during cardiac resynchronization therapy.

BACKGROUND: Both anatomic interlead separation and left ventricle lead electrical delay (LVLED) have been associated with outcomes following cardiac resynchronization therapy (CRT). However, the relationship between interlead distance and electrical delay in predicting CRT outcomes has not been defined. METHODS: We studied 61 consecutive patients undergoing CRT for standard clinical indications. All patients underwent intraprocedural measurement of LVLED. Interlead distances in the horizontal (HD), vertical (VD), and direct (DD) dimensions were measured from postprocedure chest radiographs (CXR). Remodeling indices [percent change in left ventricle (LV) ejection fraction, end-diastolic, end-systolic dimensions] were assessed by transthoracic echocardiogram. RESULTS: There was a positive correlation between corrected LVLED and HD on lateral CXR (r = 0.361, P = 0.004) and a negative correlation between LVLED and VD on posteroanterior (PA) CXR (r =-0.281, P = 0.028). To account for this inverse relationship, we developed a composite anatomic distance (defined as: lateral HD-PA VD), which correlated most closely with LVLED (r = 0.404, P = 0.001). Follow-up was available for 48 patients. At a mean of 4.1 +/- 3.2 months, patients with optimal values for both corrected LVLED (>or=75%) and composite anatomic distance (>or=15 cm) demonstrated greater reverse LV remodeling than patients with either one or neither of these optimized values. CONCLUSIONS: We identified a significant correlation between LV-right ventricular interlead distance and LVLED; additionally, both parameters act synergistically in predicting LV anatomic reverse remodeling. Efforts to optimize both interlead distance and electrical delay may improve CRT outcomes.

Decay of postextrasystolic potentiation in the left and right ventricles of intact canine hearts.

Intracellular regulation of myocardial Ca2+ has long been of interest to physiologists. The force-interval relationship provides a phenomenological approach that permits insight into aspects of calcium regulation. The response to an extrasystole is a potentiation in contractile force and the recovery in contractile force is described by the recirculation fraction (RF). The RF provides a gross estimation of calcium uptake by sarcoplasmic reticulum (SR), leading to myocardial relaxation. The current study focused on the relationship of right (RV) and left ventricular (LV) RF in canines under several contractile states. Anesthetized canines (n = 5) were catheterized for RV and LV pressure measurements. dP/dt(max) for the RV and LV was calculated for three baseline beats, one extrasystole and the first five postextrasystolic beats. The relationship between the LV dP/dt(max) and RV dP/dt(max) for all of the mentioned beats was then examined. Contractility was increased with calcium chloride and extrasystoles were delivered. Once cardiac function returned to a baseline level, contractility was reduced by increasing the concentration of isoflurane and the evaluation repeated. All ventricular contractions were controlled by RA pacing to maintain intrinsic conduction. A strong linear relationship between RV and LV dP/dt(max) (r = 0.94 +/- .06) existed for most canine's contractile states. These results build on findings in isolated hearts and demonstrate that biventricular response to extrasystoles and subsequent contractile recovery is both linear and correlated, suggesting that intracellular calcium regulation in a given heart across contractile state is static.

Left ventricular septal and left ventricular apical pacing chronically maintain cardiac contractile coordination, pump function and efficiency.

BACKGROUND: Conventional right ventricular (RV) apex pacing can lead to adverse clinical outcome associated with asynchronous activation and reduced left ventricular (LV) pump function. We investigated to what extent alternate RV (septum) and LV (septum, apex) pacing sites improve LV electric activation, mechanics, hemodynamic performance, and efficiency over 4 months of pacing. METHODS AND RESULTS: After AV nodal ablation, mongrel dogs were randomized to receive 16 weeks of VDD pacing at the RV apex, RV septum, LV apex, or LV septum (transventricular septal approach). Electric activation maps (combined epicardial contact and endocardial noncontact) showed that RV apical and RV septal pacing induced significantly greater electric desynchronization than LV apical and LV septal pacing. RV apex and RV septal pacing also significantly increased mechanical dyssynchrony, discoordination (MRI tagging) and blood flow redistribution (microspheres) and reduced LV contractility, relaxation, and myocardial efficiency (stroke work/myocardial oxygen consumption). In contrast, LV apical and LV septal pacing did not significantly alter these parameters as compared with the values during intrinsic conduction. At 16 weeks, acute intrasubject comparison showed that single-site LV apical and LV septal pacing generally resulted in similar or better contractility, relaxation, and efficiency as compared with acute biventricular pacing. CONCLUSIONS: Acute and chronic LV apical and LV septal pacing maintain regional cardiac mechanics, contractility, relaxation, and efficiency near native levels, whereas RV apical or RV septal pacing diminish these variables. Acute LV apical and LV septal pacing tend to maintain or improve contractility and efficiency compared with biventricular pacing.

Effect of transmurally heterogeneous myocyte excitation-contraction coupling on canine left ventricular electromechanics.

The excitation-contraction coupling properties of cardiac myocytes isolated from different regions of the mammalian left ventricular wall have been shown to vary considerably, with uncertain effects on ventricular function. We embedded a cell-level excitation-contraction coupling model with region-dependent parameters within a simple finite element model of left ventricular geometry to study effects of electromechanical heterogeneity on local myocardial mechanics and global haemodynamics. This model was compared with one in which heterogeneous myocyte parameters were assigned randomly throughout the mesh while preserving the total amount of each cell subtype. The two models displayed nearly identical transmural patterns of fibre and cross-fibre strains at end-systole, but showed clear differences in fibre strains at earlier points during systole. Haemodynamic function, including peak left ventricular pressure, maximal rate of left ventricular pressure development and stroke volume, were essentially identical in the two models. These results suggest that in the intact ventricle heterogeneously distributed myocyte subtypes primarily impact local deformation of the myocardium, and that these effects are greatest during early systole.

Effects of biventricular pacing and scar size in a computational model of the failing heart with left bundle branch block.

OBJECTIVES: To study the impact of biventricular pacing (BiV) and scar size on left ventricular (LV) regional and global function using a detailed finite element model of three-dimensional electromechanics in the failing canine heart. BACKGROUND: Cardiac resynchronization therapy (CRT) clinical trials have demonstrated that up to 30% of patients may be classified as non-responders. The presence of a scar appears to contribute to those that do not respond to CRT. A recent study in patients with myocardial scar showed that LV dyssynchrony was the sole independent predictor of reverse remodeling, and not scar location or size. METHODS: Two activation sequences were simulated: left bundle branch block (LBBB) and acute simultaneous BiV (with leads in the left and right ventricle) in hearts with chronic scars of various sizes. The dependence of regional function (mean fiber ejection strain, variance of fiber isovolumic strain and fraction of tissue stretched during ejection) and global function (left ventricular dP/dt(max), ejection fraction, stroke work) on scar size and pacing protocol was tested. RESULTS: Global function and regional function averaged over the whole LV during LBBB and BiV decreased as a function of scar size. In the non-scarred regions, however, regional function was largely independent of scar size for a fixed pacing site. CONCLUSIONS: The model results suggest that uniformity of mechanical contraction in non-scarred regions in the failing heart during biventricular pacing is independent of scar size for a fixed pacing site.

Computational modeling for bedside application.

Advances in computer power, novel diagnostic and therapeutic medical technologies, and an increasing knowledge of pathophysiology from gene to organ systems make it increasingly feasible to apply multiscale patient-specific modeling based on proven disease mechanisms. Such models may guide and predict the response to therapy in many areas of medicine. This is an exciting and relatively new approach, for which efficient methods and computational tools are of the utmost importance. Investigators have designed patient-specific models in almost all areas of human physiology. Not only will these models be useful in clinical settings to predict and optimize the outcome from surgery and non-interventional therapy, but they will also provide pathophysiologic insights from the cellular level to the organ system level. Models, therefore, will provide insight as to why specific interventions succeed or fail.

Effect of pacing site on systolic mechanical restitution curves in the in vivo canine model.

INTRODUCTION: Pacing site is known to influence the contractile state of the ventricle. Non-physiologic pacing sites such as the right ventricular apex (RVA) or left ventricular freewall (LVFW) have been shown to decrease the contractile state of normal myocardium, due to abnormal electrical propagation. The impact of pacing at these sites may alter mechanical restitution (MR), a fundamental cardiac property involving the electro-mechanical regulation of contraction. This, in turn, may affect cardiac function. The present study was conducted to determine if pacing site alters the time constant of MR: tau. METHODS AND RESULTS: Anesthetized canines (n = 6) were acutely paced at four sites: right atrium (RA), RVA, right ventricular septum (RVS), and LVFW. MR data was captured by the S1-S2 pacing protocol and used to create MR curves, generating a restitution time constant, tau, at each site. No significant difference in tau was found between pacing sites. A linear regression analysis of MR curves revealed that there was no significant difference in slope between pacing sites. CONCLUSION: Although pacing site has been found to influence the contractile state of the ventricle, this is the first known study to demonstrate no change in tau in an in vivo preparation. This suggests that alteration of electro-mechanical coupling described by MR is not sufficiently robust to provide insight into pacing site and cardiac function in healthy hearts.

Novel means to monitor cardiac performance: the impact of the force-frequency and force-interval relationships on recirculation fraction and potentiation ratio.

Insights into intracellular calcium regulation and contractile state can be accomplished by changing pacing rate. Steady-state increases in heart rate (HR) (force-frequency relationship, FFR), and introduction of extrasystoles (ES) (force-interval relationship, FIR) have been used to investigate this relationship. This study focused on the recirculation fraction (RF) and potentiation ratio (PR), obtained from the recovery of the FFR and FIR. These parameters may provide insight on intracellular Ca(2+) regulation. Left ventricular (LV) pressures and HR were assessed in anesthetized canines (n = 7). Intrinsic data were collected prior to and following HR increases to 150, 180, and 200 bpm, as well as following delivery of an ES at 280 ms. The RF was calculated as the slope of dP/dt(max(n + 1)) vs. dP/dt(max(n)), where n = beat number. The PR was calculated by normalizing dP/dt(max) from the first beat following the ES (or the last paced beat) to the steady-state dP/dt(max). The RF due to an ES was not significantly different than that from a HR of 200 bpm. The PR from an ES was not significantly different than from a HR of 150 bpm. The impact of an ES delivered at an interval of 280 ms produces a PR similar to that from a HR of 150 bpm; yet, it recovers similarly to the termination of pacing at 200 bpm, eliciting a similar RF value. The method of measuring RF by an ES versus an increased HR may provide a safer and more feasible approach to collecting diagnostic information.

Effects of pacing rate on mechanical restitution within the in vivo canine heart: study of the force-frequency relationship.

INTRODUCTION: It is known that as stimulation frequency is increased in a healthy heart, a corresponding increase in LV contractile function (dP/dt(max)) is observed, i.e., force-frequency relationship. The impact of this relationship on systolic and diastolic mechanical restitution in an ejecting, in vivo preparation has yet to be explored. Understanding this relationship may lead to further insight on the cellular processes that govern the contraction and relaxation of the heart, in addition to providing a safer, more feasible clinical diagnostic tool. METHODS AND RESULTS: Anesthetized canines (n = 8) were paced from the RA at rates of 130, 150, and 180 bpm. At each rate, extrasystoles were delivered at varying intervals. The LV dP/dt(max) and dP/dt(min) associated with the extrasystolic beat were expressed as a percentage of steady-state levels and plotted as a function of the extrasystolic interval to obtain mechanical restitution curves. The systolic restitution time constant length decreased significantly with all increases in heart rate, P < 0.05. In the diastolic case, significant decreases in restitution time constants were seen when heart rate was increased from 130 bpm to 180 bpm, and from 150 bpm to 180 bpm, P < 0.05. CONCLUSION: This study was the first to quantify the finding that the time constant of restitution significantly and consistently decreased with a consistent increase in heart rate. The identification of such behavior may be employed to develop stimulation protocols and chronic diagnostic tools to more safely and sensitively identify and optimize the clinical status of patients receiving pacing therapy.