AJP - Heart Calcium Transients and Cell-Sarcomere
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Am J Physiol Heart Circ Physiol 273: H2044-H2061, 1997;
0363-6135/97 $5.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Campbell, K. B.
Right arrow Articles by Slinker, B. K.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Campbell, K. B.
Right arrow Articles by Slinker, B. K.
Vol. 273, Issue 4, H2044-H2061, October 1997

MODELING IN PHYSIOLOGY
Left ventricular pressure response to small-amplitude, sinusoidal volume changes in isolated rabbit heart

Kenneth B. Campbell1,2, Yiming Wu1, Robert D. Kirkpatrick1, and Bryan K. Slinker1

1 Departments of Veterinary and Comparative Anatomy, Pharmacology, and Physiology and 2 Department of Biological Systems Engineering, Washington State University, Pullman, Washington 99164-6520

    ABSTRACT
Top
Abstract
Introduction
Methods
Conclusions
Appendix
References

The objective was to determine the dynamics of contractile processes from pressure responses to small-amplitude, sinusoidal volume changes in the left ventricle of the beating heart. Hearts were isolated from 14 anesthetized rabbits and paced at 1 beats/s. Volume was perturbed sinusoidally at four frequencies ( f ) (25, 50, 76.9, and 100 Hz) and five amplitudes (0.50, 0.75, 1.00, 1.25, and 1.50% of baseline volume). A prominent component of the pressure response occurred at the f of perturbation [in-frequency response, &Dgr;P<SUB><IT>f</IT></SUB> (t)]. A model, based on cross-bridge mechanisms and containing both pre- and postpower stroke states, was constructed to interpret &Dgr;P<SUB><IT>f</IT></SUB> (t). Model predictions were that &Dgr;P<SUB><IT>f</IT></SUB> (t) consisted of two parts: a part with an amplitude rising and falling in proportion to the pressure around that which &Dgr;P<SUB><IT>f</IT></SUB> (t) occurred [Pr(t)], and a part with an amplitude rising and falling in proportion to the derivative of Pr(t) with time. Statistical analysis revealed that both parts were significant. Additional model predictions concerning response amplitude and phase were also confirmed statistically. The model was further validated by fitting simultaneously to all &Dgr;P<SUB><IT>f</IT></SUB> (t) over the full range of f and Delta V in a given heart. Residual errors from fitting were small (R2 = 0.978) and were not systematically distributed. Elaborations of the model to include noncontractile series elastance and distortion-dependent cross-bridge detachment did not improve the ability to represent the data. We concluded that the model could be used to identify cross-bridge rate constants in the whole heart from responses to 25- to 100-Hz sinusoidal volume perturbations.

cross-bridge model; cross-bridge detachment; cross-bridge power stroke; heart muscle

    INTRODUCTION
Top
Abstract
Introduction
Methods
Conclusions
Appendix
References

PRESSURE RESPONSES to controlled volume perturbation of the left ventricle (LV) of the beating heart have long been used to characterize ventricular mechanodynamics. These pressure responses have been elicited using a variety of volume perturbation protocols including 1) sustained constant-flow volume withdrawal over periods sufficient to achieve a given volume at a specified time in the cardiac cycle (29-31, 35); 2) rapid small-volume withdrawal at each of several times during the cardiac cycle (1, 13, 25, 26, 38); 3) small-volume withdrawal of varying amplitude and rate at the time of peak systolic pressure (4, 5, 9, 27); and 4) sinusoidal volume change over the entire time course of the cardiac cycle (32, 33). Analyses of the resultant data have related pressure to volume and flow. Conclusions from all these analyses point to an organ with complex mechanodynamics that varies with time over the course of the cardiac cycle. These analyses have emphasized global features such as chamber elastance, viscous resistance, and series elastance and have made only indirect associations between these global features and the underlying muscle properties responsible for them. Recently, there has been an effort to interpret LV pressure responses directly in terms of underlying muscle mechanisms (3, 5, 7, 27). The promise from these early studies is that by using carefully controlled small-amplitude perturbations and appropriate model-based analysis, detailed kinetic behavior of cardiac muscle may be elucidated from observations of pressure responses to volume changes in the whole heart.

In this study, we pursued that promise by using one of the original volume perturbation protocols: a continuous high-frequency, small-amplitude, sinusoidal volume change delivered over the entire cardiac cycle (32, 33). We chose this protocol not only because of the history of its previous use but also because it provides a means for identifying underlying contractile behavior from observations in the whole heart; these observations extend over the entire cycle period and are not confined to a very brief 20-ms interval at the time of peak isovolumic systole as in our previous studies (5, 27, 28). Furthermore, we employed a cross-bridge model with pre- and postpower stroke elastance states to describe, predict, and explain the observed pressure responses. We concluded that these experimental and analytic techniques could be used to extract information about underlying cross-bridge mechanisms from observations made in the whole heart.

Glossary

A1 Amplitude scaling factor for pressure response component that varies proportionately with Pr(t)
A2 Amplitude scaling factor for pressure response component that varies proportionately with time drivative of &Pdot;r(t)
AIC Aikake Information Criterion
Ap Amplitude of dynamic passive component
b, d Rate constants governing formation and dissolution of prepower stroke state
Bn Amplitude of the nth harmonic of Delta Pd(t)
 Delta P(t) Pressure response to volume perturbation
 Delta Pd(t) Depressive component of pressure response
 Delta Pf(t) In-frequency component of pressure response
 Delta Pfa(t) Active part of in-frequency response
 Delta Pp(t) Dynamic passive component of pressure response
 Delta V Measured amplitude of volume perturbation
 Delta V(t) Time-varying volume perturbation
 Delta Vc Computer-commanded amplitude of volume perturbation
Ee0(t) Elastance of pressure generators in the prepower stroke state
Eep(t) Elastance of pressure generators in the postpower stroke state
ESE Elastance of series-coupled noncontractile element
ci regression coefficient
 epsilon Elastance of a single generator
f Frequency of volume perturbation
g Rate constant governing cross-bridge detachment
h Rate constant governing power stroke
K Number of parameters
 Delta lm/2 Change in half-mass wall circumference
N Number of sampled data points
Ne0(t) Number of pressure generators in the prepower stroke state
Nep(t) Number of pressure generators in the postpower stroke state
P(t) Pressure of perturbed beat
Piso(t) Pressure of isovolumic (unperturbed) beat
Pr(t) Pressure around which Delta Pf (t) occurred
Q10 Relative rate of change with a 10°C increase in temperature
 phi1 Phase of pressure response component that varies proportionately with Pr(t)
 phi2 Phase of pressure response component that varies proportionately with &Pdot;r(t)
 phip Phase of dynamic passive component
RSS Residual Sum of Squares
SC Schwartz Criterion
T Period of a heartbeat
 theta n Phase of the nth harmonic of Delta Pd(t)
VBL Baseline volume
VW Wall volume
 sigma Wall stress
 omega Angular frequency
X0 Average isovolumic distortion of postpower stroke generators
Xe0(t) Average distortion of prepower stroke generators
Xep(t) Average distortion of postpower stroke generators
Zep Total distortion among all ep generators

    EXPERIMENTAL METHODS AND PROCEDURES
Top
Abstract
Introduction
Methods
Conclusions
Appendix
References

Experimental Preparation

Hearts were isolated from 14 adult male rabbits (avg wt = 3.1 kg). Procedures for isolating the heart and attaching it to a volume-servo device have been described in detail elsewhere (7, 19). Briefly, the brachiocephalic artery was cannulated, and perfusion was begun with oxygenated relaxing solution (in mM: 121.4 Na+, 35.0 K+, 137.4 Cl-, 0.1 Ca2+, 1.1 Mg2+, 21.0 HCO<SUP>−</SUP><SUB>3</SUB>, 0.36 PO<SUP>3−</SUP><SUB>4</SUB>, and 11.1 glucose and 2.5 U/l insulin) to stop the heart before it was isolated from the animal. The perfusate was oxygenated by vigorously bubbling with 95% O2-5% CO2.

The heart was transferred to a perfusion support system consisting of a gas-exchange chamber, a roller pump, a constant-pressure chamber, and an environmental chamber. The heart was placed within an environmental chamber where the coronary arteries were perfused at 90 mmHg. Temperature was kept constant at 30°C. The heart was submerged in perfusate at all times by allowing the coronary effluent to accumulate in the environmental chamber until it reached the chamber overflow at the level of the base of the heart. The perfusate was not recirculated.

A thin latex balloon, secured to the piston cylinder of a volume-servo system, was drawn into the LV chamber such that its tip was anchored through a puncture in the apex, which also served as a vent for any fluids between the balloon and chamber wall. A draw-string suture in the mitral annulus was tightened around the obturator of a piston-cylinder device, which secured the balloon in the LV chamber. The balloon was filled with degassed distilled water until passive chamber pressure reached 10 mmHg. Balloons were sized to fill the LV without excessive folding and without developing pressure at the volumes encountered in these ventricles. Thus balloons did not contribute to measured pressure.

The perfusing solution was changed from the relaxing solution to one that allowed spontaneous beating (in mM: 148.4 Na+, 7.4 K+, 139.1 Cl-, 1.24 Ca2+, 1.1 Mg2+, 21.0 HCO<SUP>−</SUP><SUB>3</SUB>, 0.36 PO<SUP>3−</SUP><SUB>4</SUB>, and 11.1 glucose and 2.5 U/l insulin). The heart, suspended in the spent perfusate, was paced at 1 Hz with field stimulation using 5-ms pulses of 15 mV and 250 mA from 5 cm × 5 cm copper plates placed 4.5-cm apart on either side of the heart.

The volume-servo system consisted of a linear motor, a piston-cylinder device, and a linear variable differential transformer (LVDT, model 0294-0000, Trans-Tek). The piston-cylinder device was a modified 5-ml glass syringe (East Rutherford Syringes) with two side ports. One side port allowed calibrated infusion of fluid into the LV balloon to establish a baseline volume (VBL). The second port was used to introduce a 5-Fr catheter-tip pressure transducer (Millar, Houston, TX) into the balloon. The piston was driven by the armature shaft of the linear motor. Motions of the piston produced LV volume changes around VBL at a resolution of 0.001 ml. Both the pressure measurement system and the LVDT system had frequency responses of 1 kHz.

Motion of the motor armature, and consequently piston motion, was controlled to achieve specified changes in LV volume by feeding back the position signal from the LVDT transducer, comparing it with a reference position signal from a supervisory-control computer, and passing the difference through an analog proportional-integral-derivative compensator. Output from the compensator was used to drive a high-current amplifier, which delivered electrical current to the motor, causing piston position to match the volume command.

The supervisory-control computer controlled experimental protocols according to programmed instructions and also acquired data for later analysis. Pressure and volume signals were amplified to make maximal use of the 12-bit range of an analog-to-digital converter and were acquired at a 2-kHz sampling rate.

Experimental use of animals was approved by the Animal Care and Use Committee at Washington State University. The investigation conforms with the Guide for the Care and Use of Laboratory Animals published by the National Institutes of Health (NIH publication No. 85-23, Revised 1985).

Protocols

A single-beat Frank-Starling protocol (7) was conducted to establish VBL for each heart. VBL was chosen as the volume equal to 80% of the volume at which maximum pressure was developed. This protocol was also used to establish the passive pressure-volume relationship. A monoexponential equation was fit to points over the range of end-diastolic pressure and volume values generated in this protocol. Thus the contribution to pressure by parallel passive structures at any volume was estimated and removed from all ensuing data records in order to allow us to focus on just active contractile properties.

After VBL was established, a high-frequency volume perturbation protocol was conducted as follows. Twenty pairs of data records consisting of pressure and volume signals were taken. One record in a pair contained a single volume-perturbed beat, and the other record contained an unperturbed beat that served as a reference. Volume perturbation was administered only on a selected single beat. On the perturbed beat, the linear motor was commanded to deliver a sinusoidal volume change at one of four frequencies (100, 76.9, 50, or 25 Hz, corresponding to periods of 10, 13, 20, or 40 ms) and one of five amplitudes (0.5, 0.75, 1.0, 1.25, or 1.5% of VBL). Repeated records of perturbed beats were taken until all combinations of frequencies and amplitudes (20 perturbed beats) were recorded. Pressure responses to the volume perturbation were then analyzed.

Because the volume-servo system was underdamped, the actual volume perturbation did not exactly equal the commanded sinusoid from the supervisory-control computer. The frequency ( f ) of actual and commanded signals was the same, but there were differences between actual and commanded amplitudes (Delta V and Delta Vc, respectively), and there was 1-2 ms delay in the actual signal relative to the commanded signal. Consequently for some analyses (see below), each of the measured volume perturbation signals was fitted with the analytic function
&Dgr;V(<IT>t</IT>) = &Dgr;V sin (2&pgr;<IT>ft</IT> + ϕ) (1)
where phi was a phase relative to the recorded time window. Equation 1 fitted all measured Delta V(t) with a correlation coefficient (R2) > 0.99 and was thus judged to be an adequate representation of the actual perturbation signal for specific analyses. The underdamped character of the volume-servo system produced an actual Delta V that, for a given Delta Vc, increased with frequency over the 25- to 100-Hz frequency range with the result that the actual Delta V at 100 Hz was 145% of that at 25 Hz. Thus Delta V from the fit to the measured signal was used rather than the commanded Delta Vc in all data analyses.

After the high-frequency volume perturbation protocol, a second single-beat Frank-Starling protocol was conducted to generate a Frank-Starling curve that could be compared with the one collected previously. This allowed detection of any deterioration of the preparation during the course of an experiment. No detectable deterioration occurred.

    PRESSURE RESPONSE

Peak isovolumic pressure generated by these 14 hearts (averaged over all 280 observations) was 120.2 ± 13.6 (SD) mmHg at an average VBL of 2.11 ± 0.09 ml. The average LV weight, including the septum plus LV free wall, was 5.96 ± 0.54 g.

The pressure response [Delta P(t)] to Delta V(t) was defined as the difference between active pressure of the reference isovolumic beat [Piso(t)], i.e., the pressure that would have developed had no volume perturbation been administered, and active pressure of the perturbed beat, P(t)
&Dgr;P(<IT>t</IT>) = P(<IT>t</IT>) − P<SUB>iso</SUB>(<IT>t</IT>) (2)
Representative Piso(t), P(t), and Delta P(t) are shown in Fig. 1 (f = 50 Hz, Delta Vc = 1% VBL). All responses [Delta P(t)] contained two components: a depressive response, Delta Pd(t) [called "depressive" because it represented a sustained decrease in pressure below Piso(t) that was not at the perturbation frequency] and an in-frequency response [&Dgr;P<SUB><IT>f</IT></SUB>(t)] (that part of the response at the perturbation frequency). Thus
&Dgr;P(<IT>t</IT>) = &Dgr;P<SUB>d</SUB>(<IT>t</IT>) + &Dgr;P<SUB><IT>f</IT></SUB>(<IT>t</IT>) (3)
The Delta Pd(t) was extracted from Delta P(t) by fitting a curve to Delta P(t) that did not contain frequency content of the perturbation frequency. Delta Pd(t) was taken as the sum of the first 10 harmonics in the Fourier series
&Dgr;P<SUB>d</SUB>(<IT>t</IT>) = <LIM><OP>∑</OP><LL><IT>n</IT>=1</LL><UL>10</UL></LIM> <IT>B</IT><SUB><IT>n</IT></SUB> sin <FENCE><IT>n</IT> <FR><NU>2&pgr;</NU><DE><IT>T</IT></DE></FR> + &thgr;<SUB><IT>n</IT></SUB></FENCE> (4)
where n is the harmonic number, Bn and theta n are harmonic amplitude and phase, respectively, of n, and T is the heart period. Because the shortest heart period used in these studies was 1 s, the 10th harmonic (10 Hz) was well below 25 Hz, the lowest frequency used for volume perturbation. The amplitude and phase parameters for the ith harmonic (Bi and &thgr;<SUB><IT>i</IT></SUB>) had no particular significance other than to give a Delta Pd(t) waveshape, identifiable within Delta P(t), that did not include components of the in-frequency response. Once Delta Pd(t) was identified by fitting with Eq. 4, it was subtracted from Delta P(t) to yield &Dgr;P<SUB><IT>f</IT></SUB> (t). Subtraction of Delta Pd(t) from Piso(t) generated a signal representing the pressure around which &Dgr;P<SUB><IT>f</IT></SUB> (t) took place [Pr(t)]. Pr(t), with corresponding Delta P(t), Delta Pd(t), and &Dgr;P<SUB><IT>f</IT></SUB> (t), is shown in Fig. 1.


View larger version (18K):
[in this window]
[in a new window]
 


View larger version (14K):
[in this window]
[in a new window]
 
Fig. 1.   Method of determining pressure response [volume perturbation = 50 Hz, 1% baseline volume (VBL)]. A: 2 left panels: Piso(t), pressure of an isovolumic beat in which no volume perturbation was applied; P(t), pressure of a beat that received volume perturbation. Middle panel: Delta P(t), pressure response to volume perturbation [= P(t- Piso(t)]. Right panel: Delta Pd(t), depressive component of Delta P(t) obtained by low-pass filtering Delta P(t). Bottom panel: &Dgr;P<SUB><IT>f</IT></SUB> (t), in-frequency component of Delta P(t) [=Delta P(t- Delta Pd(t)]. B: pressure around which in-frequency response occurred [Pr(t)] was obtained by subtracting Delta Pd(t) from Piso(t). Vertical scale is in mmHg, horizontal scale is in s.

This report concerns just &Dgr;P<SUB><IT>f</IT></SUB> (t); the Delta Pd(t) is the subject of another report (unpublished observations). Twelve of twenty &Dgr;P<SUB><IT>f</IT></SUB> (t) responses obtained in one heart are shown in Fig. 2. &Dgr;P<SUB><IT>f</IT></SUB> (t) rose and fell during the course of the heartbeat but also contained a small contribution that was present during diastolic periods when there was no active contraction. This small component was assumed to be due to dynamic features of parallel passive properties that were not included in the static passive pressure-volume relationship, which had already been subtracted from the response. Furthermore, these passive dynamic properties were assumed to be expressed continuously throughout the period of the heartbeat and, for a given Delta V(t) sinusoid, they were represented as
&Dgr;P<SUB>p</SUB>(<IT>t</IT>) = <IT>A</IT><SUB>p</SUB> sin (2&pgr;<IT>ft</IT> + ϕ<SUB>p</SUB>) (5)
This passive dynamic response was subtracted from the response signal using data-fitting procedures described below. The great majority of &Dgr;P<SUB><IT>f</IT></SUB>(t) waxed and waned as pressure rose and fell and was considered to be an expression of active processes. Thus
&Dgr;P<SUB><IT>fa</IT></SUB>(<IT>t</IT>) = &Dgr;P<SUB><IT>f</IT></SUB>(<IT>t</IT>) − &Dgr;P<SUB>p</SUB>(<IT>t</IT>) (6)
where &Dgr;P<SUB><IT>fa</IT></SUB>(t) was the contribution of active process to &Dgr;P<SUB><IT>f</IT></SUB> (t).


View larger version (33K):
[in this window]
[in a new window]
 
Fig. 2.   &Dgr;P<SUB><IT>f</IT></SUB> (t) over single heartbeat for 12 of 20 responses in one heart, showing extremes and midrange of responses to various amplitudes and frequencies of Delta V(t). Note growth in amplitude of response with both amplitude and frequency of Delta V(t). A: 0.5%; B: 1.0%; C: 1.5%.

It is clear from Fig. 1 that, in accordance with results from several studies (1, 13, 25, 29, 32, 38), &Dgr;P<SUB><IT>fa</IT></SUB>(t) waxed and waned during the heart period as Pr(t) rose and fell. One objective of the current work was to predict from model considerations whether other time-varying components were contained within &Dgr;P<SUB><IT>fa</IT></SUB>(t) and, then, to test these predictions.

    MODEL DESCRIPTION

A model for describing and predicting the active part of the pressure response [&Dgr;P<SUB><IT>fa</IT></SUB>(t)] was constructed on the assumption that elements responsible for force generation in cardiac muscle (i.e., cross-bridges between thick and thin myofilaments) were also responsible for pressure generation in the LV chamber. Furthermore, it was assumed that there was a straightforward linear transformation between force-length relationships in the wall of the heart and pressure-volume relationships in the LV chamber. Acceptability of the linear transformation assumption requires small-amplitude perturbations and homogenous myocardium. Criteria for satisfying the small-amplitude requirement are detailed in the APPENDIX. Evidence that the homogenous myocardium requirement is satisfied is given in the findings relative to the unimportance of noncontractile series elasticity in these hearts (see MODEL VALIDATION). However, a strong reason for employing the linear transformation assumption is the success that has been achieved with its application in earlier studies (3, 5, 6, 27).

It was further assumed that during a heartbeat mechanodynamics were from two sources: 1) dynamics of activation as activator Ca2+ comes and goes and numbers of force-bearing cross-bridges rise and fall, and 2) dynamics of cross-bridge cycling as myosin heads cyclically attach to and detach from the actin binding site. In accordance with an earlier hypothesis (6), we argue that the only dynamics expressed within the brief cycle period of frequencies >= 25 Hz were those associated with steps in the cross-bridge cycle and that the dynamics of activation were too slow to contribute to changes within these brief time periods. Such separation of time scales in the study of muscle dynamics is in accordance with analyses conducted by Kawai and co-workers (17, 24, 39) and in accordance with our recent demonstration that cooperativity between force-bearing cross-bridges and activation can cause activation to be slow relative to cross-bridge dynamics (2).

We refer to cross-bridges as force generators. Generators contributing to the pressure response were assumed to be in two states: 1) a state that possessed elastance but did not, under isometric (isovolumic) conditions, generate pressure (state e0), and 2) a state that both possessed elastance and also generated isometric (isovolumic) pressure (state ep). When we assume linear, independent, and parallel generators, the elastance associated with each state is the number of parallel generators in that state (N) times the elastance of a single generator (epsilon ). Assuming that all generators in states e0 and ep possess the same epsilon , we show the net elastance of all parallel generators in each of the two states as
<IT>E</IT><SUB><IT>e0</IT></SUB>(<IT>t</IT>) = &egr; ⋅ <IT>N</IT><SUB><IT>e0</IT></SUB>(<IT>t</IT>) (7)
<IT>E</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>) = &egr; ⋅ <IT>N</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>) (8)
where, because of our linear transformation assumption, Ee0(t) and Eep(t) may be taken as volumetric elastances (with units of mmHg/ml).

Generators are in continual transition as they progress from one state to another in the cross-bridge cycle (Fig. 3). We assumed that state e0 preceded state ep. Furthermore, transitions into, out of, and between states e0 and ep were assumed to be governed by rate constants b, d, g, and h. State 0 in Fig. 3 is without elastance and is a precursor to state e0; state 00 is also without elastance and follows state ep. These nonelastance states represent all other states needed to complete a cross-bridge cycle. Given these relationships between states and assuming that there is no noncontractile series elastance (4), it is shown in the APPENDIX that Ee0(t) and Eep(t) may be calculated from Pr(t) and its first time derivative [&Pdot;r(t)] according to
<IT>E</IT><SUB><IT>e0</IT></SUB>(<IT>t</IT>) = <FR><NU><A><AC>P</AC><AC>˙</AC></A><SUB>r</SUB>(<IT>t</IT>) + <IT>g</IT>P<SUB>r</SUB>(<IT>t</IT>)</NU><DE><IT>hX</IT><SUB>0</SUB></DE></FR> (9)
and
<IT>E</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>) = <FR><NU>P<SUB>r</SUB>(<IT>t</IT>)</NU><DE><IT>X</IT><SUB>0</SUB></DE></FR> (10)
where X0 is a parameter representing average volumetric distortion among generators in state ep during isovolumic conditions. The transitional step between states e0 and ep, which is governed by h, is the cross-bridge power stroke, and this step is responsible for inducing X0 distortion in generators as they go through the power stroke to enter the ep state. The dissolution of the postpower stroke (ep) state, which is governed by g, is the cross-bridge detachment step.


View larger version (17K):
[in this window]
[in a new window]
 
Fig. 3.   Schematic drawing of states in cross-bridge cycle. Ni is the number of generators in ith state. Generators in states e0 and ep possess elastance, whereas generators in states 0 and 00 do not. Generators possessing elastance may be distorted during a volume perturbation such that both contribute to pressure response. Under isovolumic conditions, generators enter state e0 without distortion and do not generate pressure. Transitions between states are governed by rate constants b, d, h, and g. Transition between state e0 and state ep is the power stroke and induces a baseline distortion in postpower stroke (ep) generators, which, as a result of their elastance, causes development of isovolumic pressure. Isovolumic pressure is modified during a volume perturbation by induced distortion in postpower stroke and prepower stroke states, ep and e0, respectively, and by whatever influence volume perturbation has on recruitment of generators into or out of cross-bridge cycle.

Because of their elastic nature, generators within each of these states are distorted during volume perturbation. The net volume-induced distortion is determined by the rate of volume change relative to the rates of formation and dissolution of the respective states. Differential equations for these volume-induced distortions [Delta Xe0(t) and Delta Xep(t) in states e0 and ep, respectively] are derived in the APPENDIX as
&Dgr;<IT><A><AC>X</AC><AC>˙</AC></A></IT><SUB><IT>e0</IT></SUB>(<IT>t</IT>) = − <FENCE><FR><NU><IT><A><AC>E</AC><AC>˙</AC></A></IT><SUB><IT>e0</IT></SUB>(<IT>t</IT>)</NU><DE><IT>E</IT><SUB><IT>e0</IT></SUB>(<IT>t</IT>)</DE></FR> + <IT>h</IT> + <IT>d</IT></FENCE> &Dgr;<IT>X</IT><SUB><IT>e0</IT></SUB>(<IT>t</IT>) + &Dgr;<A><AC>V</AC><AC>˙</AC></A>(<IT>t</IT>) (11)
&Dgr;<IT><A><AC>X</AC><AC>˙</AC></A></IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>) = − <FENCE><FR><NU><IT><A><AC>E</AC><AC>˙</AC></A></IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>)</NU><DE><IT>E</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>)</DE></FR> + <IT>g</IT></FENCE> &Dgr;<IT>X</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>) + &Dgr;<A><AC>V</AC><AC>˙</AC></A>(<IT>t</IT>) (12)
A dot over a variable indicates its first time derivative. In words, Eqs. 11 and 12 state that the time rate of change of generator distortion is negatively related to distortion itself and is positively related to the first time derivative of volume. Thus volume-induced distortion is driven directly by flow (velocity) and dies away at a rate proportional to the distortion.

Given these relationships, the model equation predicting the active part of the in-frequency pressure response [Delta &Pcirc;fa(t)] is
&Dgr;<A><AC>P</AC><AC>ˆ</AC></A><SUB><IT>fa</IT></SUB>(<IT>t</IT>) = <IT>E</IT><SUB><IT>e0</IT></SUB>(<IT>t</IT>)&Dgr;<IT>X</IT><SUB><IT>e0</IT></SUB>(<IT>t</IT>) + <IT>E<SUB>ep</SUB></IT>(<IT>t</IT>)&Dgr;<IT>X</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>) (13)
where the "hat" in Delta &Pcirc;fa(t) indicates that the quantity is model predicted.

Equations 9-13 constitute a 2-state, 4-parameter model. In this model, Delta V(t) is the input [although distortion is driven directly by Delta V(t)]; Ee0(t), Eep(t), and their time derivatives combine to form time-varying parameters; Delta Xe0(t) and Delta Xep(t) are state variables; and Delta &Pcirc;fa(t)is the output variable. Equations 11 and 12 are a set of linear, uncoupled, first-order, time-varying differential equations.

    MODEL PREDICTIONS

Two Dynamic Components of In-Frequency Response

In-frequency response consists of two dynamic components: 1) a component with an amplitude varying with Pr(t) and 2) a component with an amplitude varying with the time derivative of Pr(t). The model output equation Eq. 13 can be rearranged by substituting elastance Eqs. 9 and 10 into Eq. 13 to create an alternative formulation.
&Dgr;<A><AC>P</AC><AC>ˆ</AC></A><SUB><IT>fa</IT></SUB>(<IT>t</IT>) = P<SUB>r</SUB>(<IT>t</IT>) <FENCE><FR><NU>1</NU><DE><IT>X</IT><SUB>0</SUB></DE></FR></FENCE><FENCE><FR><NU><IT>g</IT></NU><DE><IT>h</IT></DE></FR> &Dgr;<IT>X</IT><SUB><IT>e0</IT></SUB>(<IT>t</IT>) + &Dgr;<IT>X</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>)</FENCE>
+ <A><AC>P</AC><AC>˙</AC></A><SUB>r</SUB>(<IT>t</IT>) <FENCE><FR><NU>1</NU><DE><IT>hX</IT><SUB>0</SUB></DE></FR></FENCE> &Dgr;<IT>X</IT><SUB><IT>e0</IT></SUB>(<IT>t</IT>) (14)

The first term on the right-hand side of Eq. 14 is a response component with an amplitude varying proportionately with Pr(t), and the second term on the right-hand side is a component with an amplitude varying proportionately with &Pdot;r(t). This development clearly identifies the contribution of the prepower stroke, e0 state as the sole source of the dynamic response component with an amplitude rising and falling in proportion to the derivative of the pressure around which the response is occurring.

To determine the relative roles of these two components, an approximate sinusoidal solution of the model equations Eqs. 11 and 12 was developed as follows. When we ignore the influence of the time-varying part of the coefficients in Eqs. 11 and 12, a steady-state solution of these equations for Delta Xe0(t) and Delta Xep(t) when Delta V(t) is a volume sinusoid of frequency f results in an expression for the first and second terms on the right-hand side of Eq. 14 as
<FENCE><FR><NU>1</NU><DE><IT>X</IT><SUB>0</SUB></DE></FR></FENCE><FENCE><FR><NU><IT>g</IT></NU><DE><IT>h</IT></DE></FR> &Dgr;<IT>X</IT><SUB><IT>e0</IT></SUB>(<IT>t</IT>) + &Dgr;<IT>X</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>)</FENCE> = <IT>A</IT><SUB>1</SUB> sin (2&pgr;<IT>ft</IT> + ϕ<SUB>1</SUB>)
<FENCE><FR><NU>1</NU><DE><IT>hX</IT><SUB>0</SUB></DE></FR></FENCE> &Dgr;<IT>X<SUB>e0</SUB></IT>(<IT>t</IT>) = <IT>A</IT><SUB>2</SUB> sin (2&pgr;<IT>ft</IT> + ϕ<SUB>2</SUB>) (15)
Such that for a volume sinusoid of frequency f
&Dgr;<A><AC>P</AC><AC>ˆ</AC></A><SUB><IT>fa</IT></SUB>(<IT>t</IT>) = P<SUB>r</SUB>(<IT>t</IT>)[<IT>A</IT><SUB>1</SUB> sin (2&pgr;<IT>ft</IT> + ϕ<SUB>1</SUB>)]
+ <A><AC>P</AC><AC>˙</AC></A><SUB>r</SUB>(<IT>t</IT>)[<IT>A</IT><SUB>2</SUB> sin (2&pgr;<IT>ft</IT> + ϕ<SUB>2</SUB>)] (16)

The relative role of the &Pdot;r(t) component in the observed responses was determined by an incremental approach. Equations 6 and 16 were combined and fit to the observed &Dgr;P<SUB><IT>f</IT></SUB> (t) by first excluding the &Pdot;r(t) component
&Dgr;<A><AC>P</AC><AC>ˆ</AC></A><SUB><IT>f</IT>1</SUB>(<IT>t</IT>) = <IT>A</IT><SUB>p</SUB> sin (2&pgr;<IT>ft</IT> + ϕ<SUB>p</SUB>)
+ P<SUB>r</SUB>(<IT>t</IT>)[<IT>A</IT><SUB>1</SUB> sin (2&pgr;<IT>ft</IT> + ϕ<SUB>1</SUB>)] (17)
and then including the &Pdot;r(t) component
&Dgr;<A><AC>P</AC><AC>ˆ</AC></A><SUB><IT>f</IT>2</SUB>(<IT>t</IT>) = <IT>A</IT><SUB>p</SUB> sin (2&pgr;<IT>ft</IT> + ϕ<SUB>p</SUB>)
+ P<SUB>r</SUB>(<IT>t</IT>)[<IT>A</IT><SUB>1</SUB> sin (2&pgr;<IT>ft</IT> + ϕ<SUB>1</SUB>)]
+ <A><AC>P</AC><AC>˙</AC></A><SUB>r</SUB>(<IT>t</IT>)[<IT>A</IT><SUB>2</SUB> sin (2&pgr;<IT>ft</IT> + ϕ<SUB>2</SUB>)] (18)
Fitting of Eqs. 17 and 18 to &Dgr;P<SUB><IT>f</IT></SUB> (t) was by an heuristic search algorithm (Levenberg-Marquardt algorithm, Argonne National Laboratory) to minimize the residual sum of squares (RSS).

To test for degradation or improvement in the representation of &Dgr;P<SUB><IT>f</IT></SUB> (t) signal with the addition of the two additional parameters in Eq. 18 that relate to the &Pdot;r(t) component, the Aikake Information Criterion (AIC) and the Schwartz Criterion (SC) were calculated from fits with both Eqs. 17 and 18 according to Landaw and DiStefano (20) as
AIC = <IT>N</IT> ln (RSS) + 2<IT>K</IT>
SC = <IT>N</IT> ln (RSS) + <IT>K</IT> ln (<IT>N</IT>) (19)
where N is the number of sampled data points, and K is the number of parameters. Note that the first term on the right-hand side of Eq. 19 is a measure of how well the model fit the data, whereas the second term is a penalty function based on the number of model parameters. Therefore, increasing the number of parameters increases AIC and SC unless there is a more than compensating reduction in RSS. In considering two competing equations such as Eqs. 17 and 18, the better representation is the one with the smallest AIC and SC. To further determine whether significant reduction in the RSS occurred with Eq. 18 compared with Eq. 17, an incremental F-test was used (11).

When fit to &Dgr;P<SUB><IT>f</IT></SUB> (t) of each of the 280 data records obtained in these hearts (14 hearts times 20 records/heart), both Eq. 17, which did not include the &Pdot;r(t) component, and Eq. 18, which did include this component, fit &Dgr;P<SUB><IT>f</IT></SUB> (t) very well with median R2 of 0.980 and 0.981, respectively. The contribution of the &Pdot;r(t) component was quite small as judged by the fact that A2 was always two-orders of magnitude less than A1. However, despite this small contribution, the inclusion of the &Pdot;r(t) component in Eq. 18 consistently reduced the AIC (only one exception in 280 instances; median reduction: -0.78%; range: 0.02 to -10.1%) and SC (14 exceptions in 280 instances; median reduction: -0.63%; range: 0.09 to -9.93%), suggesting that the addition of the two parameters associated with the &Pdot;r(t) component improved the representation of the information in the &Dgr;P<SUB><IT>f</IT></SUB> (t) signals. Furthermore, of the 280-response records analyzed, the incremental F-test generated an F-statistic that was significant at the P < 0.01 level in all 280 instances. Thus, although making only a small contribution in accounting for the total variability in &Dgr;P<SUB><IT>f</IT></SUB> (t), the &Pdot;r(t) component contributed significantly to representing its information content and in reducing the RSS.

To summarize these results, the model predicted that there would be an in-frequency response component with an amplitude rising and falling in proportion to the pressure around which the response occurred [Pr(t)] and another component with an amplitude rising and falling in proportion to the derivative of the pressure around which the response occurred [&Pdot;r(t)]. Analysis of all the response data revealed that the response was dominated by the Pr(t) component, although a small but significant component existed with an amplitude of which was proportional to &Pdot;r(t). Given the very small contribution by the &Pdot;r(t) component, the model was further used to test the importance of including this term in validation procedures described below.

Amplitude Ratio and Phase of Response

Additional model predictions resulted from considering the nature of the response just around the time of peak Pr(t), when &Pdot;r(t) approximated zero. When transients are ignored and it is assumed that steady state had been achieved at this time, an argument can be made that during a short interval around the time of peak pressure, the e0 state does not contribute to the response and the model reduces to
&Dgr;<A><AC>P</AC><AC>ˆ</AC></A><SUB><IT>fa</IT></SUB>(<IT>t</IT>) = <IT>E</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>)&Dgr;<IT>X<SUB>ep</SUB></IT>(<IT>t</IT>)
&Dgr;<IT><A><AC>X</AC><AC>˙</AC></A></IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>) = −<IT>g</IT>&Dgr;<IT>X</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>) + &Dgr;<A><AC>V</AC><AC>˙</AC></A>(<IT>t</IT>) (20)
Input-output relationships between Delta Xep(t) and Delta V(t) for this reduced model may be derived as
<FR><NU>&Dgr;<IT>X</IT><SUB><IT>ep</IT></SUB>(&ohgr;)</NU><DE>&Dgr;V(&ohgr;)</DE></FR> = <FENCE><FR><NU>&Dgr;<IT>X</IT><SUB><IT>ep</IT></SUB>(&ohgr;)</NU><DE>&Dgr;V(&ohgr;)</DE></FR></FENCE><IT>e </IT><SUP><IT>j</IT>&thgr;(&ohgr;)</SUP>
<FENCE><FR><NU>&Dgr;<IT>X</IT><SUB><IT>ep</IT></SUB>(&ohgr;)</NU><DE>&Dgr;V(&ohgr;)</DE></FR></FENCE> = <FR><NU>&ohgr;</NU><DE><RAD><RCD>&ohgr;<SUP>2</SUP> + <IT>g</IT><SUP>2</SUP></RCD></RAD></DE></FR>
&thgr;(&ohgr;) = 90° − tan<SUP>−1</SUP> <FENCE><FR><NU>&ohgr;</NU><DE><IT>g</IT></DE></FR></FENCE> (21)
where omega  is angular frequency in radians and equals 2pi f, |&Dgr;<IT>X</IT><SUB>ep</SUB>(omega )/Delta V(omega )| is the magnitude of an input-output amplitude ratio equivalent to A1/Delta V from the previous sinusoidal fits, and theta (omega ) is the phase difference between output pressure response and input volume sinusoid.

Two predictions result from Eq. 21: 1) the amplitude ratio will increase with frequency up to some plateau, provided the frequencies examined are in the vicinity of the characteristic frequency, g. At frequencies either far below or far above g, the amplitude ratio will change only weakly with frequency. 2) The phase of &Dgr;P<SUB><IT>f</IT></SUB> (t) will lead the phase of Delta V(t) by as much as 90° at low frequencies, but this phase lead will decline and approach zero as frequencies increase above g.

To test these model predictions, the amplitude ratio, A1/Delta V, was evaluated for its dependence on f and Delta V. As noted above, A1/Delta V will change sharply with f over a frequency range around the characteristic frequency of the underlying process. Additionally, A1/Delta V is not expected to be dependent on the amplitude of the Delta V input. Any dependence of A1/Delta V on the amplitude of the input is an indication of nonlinear processes that are not part of the current model. Nonlinearities may also show up as dependence of A1/Delta V on product combinations of f and Delta V. Stepwise regression analysis was used for these determinations. Regression equations were formulated as
<FR><NU>&Dgr;<IT>A</IT><SUB>1</SUB></NU><DE>&Dgr;V</DE></FR> = <IT>c</IT><SUB>0</SUB> + <IT>c</IT><SUB>1</SUB> <IT>f</IT> + <IT>c</IT><SUB>2</SUB>&Dgr;V + <IT>c</IT><SUB>3–7</SUB> ( <IT>f</IT>, &Dgr;V) (22)
where the ci values are regression coefficients and (f,Delta V) represents one or more of four candidate interaction terms: f · (Delta V), root-mean-squared (rms) flow; f  2 · (Delta V), rms acceleration; Delta V2, squared rms volume amplitude; and (f · Delta V)2, squared rms flow amplitude. The regression procedure used dummy variables and effects coding to account for between-subjects differences, and the subject dummy variables were forced into the stepwise regression (11). A candidate predictor variable was considered significant only when the P value for its inclusion was <0.05.

The dependence of A1/Delta V on f at the various commanded Delta Vc for one heart is shown in Fig. 4. At all Delta Vc, A1/Delta V increased with f. Furthermore, at f equal to 25, 50, and 76.9 Hz, A1/Delta V had virtually no dependence on Delta Vc; there was an apparent small dependence on Delta Vc at 100 Hz. Regression analysis of pooled data from all hearts revealed that, of all potential predictor variables, there was a significant dependence on f and Delta V2; the simplest best regression equation (leaving out terms for between subjects variability) was
<FR><NU><IT>A</IT><SUB>1</SUB></NU><DE>&Dgr;V</DE></FR> = 2.42 + 0.0113<IT>f</IT> − 311&Dgr;V<SUP>2</SUP>
<IT>R</IT><SUP>2</SUP>(adjusted) = 0.87 (23)
No interaction terms were found to be significant in this regression. Of the two significant predictor variables, f was the more important variable. For example, at f = 50 Hz and Delta V = 1%, the term in Eq. 23 due to f adds 0.55 units to A1/Delta V, whereas the term due to Delta V2 subtracts only 0.03 units, an 18-fold difference in sensitivity of A1/Delta V on these variables. The three important outcomes of this regression analysis are 1) A1/Delta V increased with f in accordance with Delta V(t) acting as the effective driving function; 2) because of the strong dependence on f, the characteristic frequency of the dynamic processes responsible for the Pr(t) response component lies either within or not far distant from the 25- to 100-Hz frequency range; and 3) the very small dependence of the Pr(t) response component on Delta V indicated that nonlinearities were not a large part of the underlying processes.


View larger version (7K):
[in this window]
[in a new window]
 
Fig. 4.   Amplitude ratio (A1/Delta V) of &Dgr;P<SUB><IT>f</IT></SUB> (t) component proportional to Pr(t) as a function of frequency (f) and commanded amplitude (Delta Vc) of perturbation. open circle , Delta Vc = 0.50% VBL; +, Delta Vc = 0.75% VBL; ×, Delta Vc = 1.00% VBL; * Delta Vc = 1.25% VBL; bullet , Delta Vc = 1.50% VBL. A1/Delta V increases with frequency and has very little dependence on Delta Vc, as predicted by model.

Predictions with regard to the phase lead were not analyzed exhaustively. Rather, single cycles, spanning the time of peak Pr(t) of response to 1% Delta Vc at each frequency, were examined. In these, the phase of &Dgr;P<SUB><IT>f</IT></SUB> (t) led that of Delta V(t) at 25 Hz by ~30-40°, and this phase lead decreased progressively to approach zero at 76.9 and 100 Hz.

To summarize these results, model predictions with regard to frequency dependence of input-output amplitude ratio and phase relations were confirmed. These indirect confirmations of the model suggested a more rigorous model validation test.

    MODEL VALIDATION

Ability to Fit the Data

Model validation was, in part, by evaluating how well the model fit the full time course of the pressure response over an entire cardiac cycle. Model fitting was by the following procedure. Initial values of the four model parameters (g, h, d, and X0) and the two dynamic passive pressure parameters (Ap and phip) were assigned. Derivatives of measured Delta V(t) and Pr(t) were calculated using a five-point Lagrangian polynomial method. Measured Pr(t), the calculated derivatives, and the then-current parameter values were fed into the differential equations Eqs. 11 and 12, allowing these equations to be solved numerically by integrating with a fourth-order Runge-Kutta algorithm (integration step size = 0.0005 s) to obtain predictions of Delta Xe0(t) and Delta Xep(t). These were then used with measured Pr(t) and its derivative to compute a &Dgr;P<SUB><IT>fa</IT></SUB>(t) according to Eq. 13 and a Delta Pp(t) according to Eq. 5, which were then added to obtain a &Dgr;P<SUB><IT>f</IT></SUB> (t). The RSS between predicted and observed &Dgr;P<SUB><IT>f</IT></SUB> (t) was calculated. Values of model parameters (g, h, d, and X0) and dynamic passive pressure parameters (Ap and phip) were then adjusted according to the rules of a Levenberg-Marquardt heuristic search algorithm, and the processes were repeated iteratively until RSS was minimized. Median model parameter values obtained with this procedure in the 14 hearts are reported in Table 1.

                              
View this table:
[in this window]
[in a new window]
 
Table 1.   Median model parameter values

Unlike the sinusoidal approximation of Eq. 18, which could be fit only to individual responses to a single perturbation, model Eqs. 9-13 are more general and could be fit to groups of responses to perturbations of multiple frequencies and amplitudes. By fitting simultaneously to responses to the 20 perturbations imposed in any one heart, we required this single model with a single set of parameters to account for a wide range of behaviors resulting from many different perturbations as in Fig. 2. Successful reproduction of this broad range of dynamic responses was taken as compelling evidence for validity of the basic model.

By all measures, the basic 2-state, 4-parameter model fit the response data very well. An example of this good fit in one heart is shown for a single heartbeat ( f = 50 Hz; Delta Vc = 1% VBL) in Fig. 5. In evaluating Fig. 5, it must be kept in mind that the fit that generated the predicted response was to responses obtained from the complete set of four frequencies and five amplitudes of inputs delivered to each of 20 beats and not just to the single beat shown in Fig. 5. Because rapid cycling at 50 Hz generated a dense pattern on display from which model-predicted and measured waveforms could not be discriminated, individual cycles in the response were identified that 1) spanned the point on the ascending limb of Pr(t) at which Pr(t) = 1/2 its peak value, 2) spanned the peak value of Pr(t), and 3) spanned the point on the descending limb at which Pr(t) = 1/2 its peak value. These three cycles were then expanded in row B of Fig. 5 such that predicted and observed waveforms could be compared. The comparatively small values of the differences between predicted and observed waveforms (residuals) are given in row C of Fig. 5. In this particular example, but not true in all cases, the residuals appeared to be random at all times during the heart period and exhibited no transient systematic character. Systematic patterns in the residuals will exhibit as a periodicity at the frequency of perturbation.


View larger version (28K):
[in this window]
[in a new window]
 
Fig. 5.   A: response to a single perturbation f = 50 Hz; Delta Vc = 1% VBL. Model was fit to complete set of 4 frequencies and 5 amplitudes of inputs delivered to 20 beats, not just to single beat shown. B: model-predicted vs. measured waveforms. These overlain waveforms cannot be discriminated in left panel. Waveforms of single cycles (identified by period between pairs of vertical lines in left panel) are displayed in 3 panels to right. These cycles spanned: 1) the point on ascending limb of Pr(t) at which Pr(t) = 1/2 its peak value, 2) the peak value of Pr(t), and 3) the point on descending limb at which Pr(t) = 1/2 its peak value. C: residuals between model-predicted and observed waveforms. Residuals are small valued and do not have periodicity of perturbation. Thus, in this example, they are apparently random.

To demonstrate the goodness of the fit over the full set of 20 perturbed beats collected in the single heart featured in Fig. 5, predicted and measured values were plotted against one another to generate Fig. 6. In keeping with an R2 value of 0.98 in this heart, all the points cluster tightly around a line that is not clearly distinguishable from the line of identity. It can be seen that deviations of predicted from observed &Dgr;P<SUB><IT>fa</IT></SUB>(t) were never large, and an excellent fit was achieved at all frequencies and amplitudes of perturbation.


View larger version (17K):
[in this window]
[in a new window]
 
Fig. 6.   Model predicted vs. measured &Dgr;P<SUB><IT>fa</IT></SUB>(t) for all 20 perturbed beats in one heart; 40,000 points shown (20 beats, each of 1-s duration, sampled at 2 kHz). Points in densest part of cluster around origin represent data generated at smallest Delta Vc, at lowest f, at lowest values of Pr(t), and during zero crossings at all perturbations. Points in less dense vertices of cluster represent data generated at highest Delta Vc, at highest f, and at peak of Pr(t). All points cluster tightly around line of identity. Deviations of predicted from observed &Dgr;P<SUB><IT>fa</IT></SUB>(t) were never large, and an excellent fit was achieved at all frequencies and amplitudes of perturbation.

An additional demonstration of model goodness of fit comes from comparing the actual responses displayed in Fig. 2 with the corresponding model-predicted responses displayed in the same format in Fig. 7. Visual comparison of these two figures reveals no important differences. Yet another demonstration of correspondence between model prediction and measured responses is seen in Fig. 8, where it is shown that there was good agreement between A1/Delta V for all 20 responses from the sinusoidal analysis and the 20 model-derived equivalents. This agreement was with respect to both absolute magnitudes and to systematic variations with f and Delta V, strong dependence on f, and little or no dependence on Delta V.


View larger version (30K):
[in this window]
[in a new window]
 
Fig. 7.   Alternative view of some of data plotted in Fig. 6. Twelve of 20 model-predicted responses are shown, corresponding to 12 measured responses shown in Fig. 2. A: 0.5%; B: 1.0%; C: 1.5%.


View larger version (7K):
[in this window]
[in a new window]
 
Fig. 8.   bullet , Model-predicted equivalent of A1/Delta V for each of 20 volume perturbations. There is virtually no dependence of model-predicted A1/Delta V on perturbation amplitude but a strong dependence on frequency. open circle , A1/Delta V from sinusoidal analysis, i.e., same data shown in Fig. 4.

That good fits were obtained in all 14 hearts is supported by a median R2 value of 0.978 with a narrow range between 0.968 and 0.987. These R2 values were slightly smaller than the R2 value obtained when the sinusoidal approximation (Eq. 18) was fit to a single response; these R2 values often exceeded 0.98. However, it must be kept in mind that the sinusoidal equation was fit individually to each of the 20 responses in each heart, and the resultant R2 values reflected the ability of the sinusoidal equation to reproduce only the variability seen in a single response with no consideration of the causal relation between Delta V(t) and the response. This contrasts with the R2 obtained with the model in which causal relationships between Delta V(t)and &Dgr;P<SUB><IT>fa</IT></SUB>(t) were stipulated and in which the model was asked to reproduce not just the variability in a single response but rather the variability due to changes in Delta V and f in the input Delta V(t) over the whole family of 20 responses in each heart. Given the resulting range of response variation to these broad ranges in input Delta V and f (Fig. 2), the median R2 value of 0.978 obtained with the basic model is remarkable.

This good fit by the basic model left only 2.2% of the total variation in the data unaccounted. Of this 2.2% it was very difficult to ascertain how much was simply random noise that could not be reduced with model improvement and how much was systematic and subject to reduction with model improvement. Systematic error would appear as systematic variation in the residuals at the frequency of perturbation, whereas random error would appear as residual variation unrelated to the frequency of perturbation. Unfortunately, the residuals did not lend themselves to analysis because systematic error, when it could be observed, did not persist throughout the period of the cardiac cycle. Rather, detectable systematic error would appear briefly during, say, just early systole or just during late systole and then disappear into apparently random error for the rest of the cycle. Such comings and goings of systematic error made it impossible to quantify the degree to which the residuals were random. Qualitatively, our impression was that systematic variation in the residuals was only a small part of the total residual variation. Because there was only a very small systematic error, our impression was that there was very little room for improvement with model modifications. However, results of this analysis do not eliminate the possibility that other features are present or even that a simpler model may be as good.

Basic Model vs. Competing Alternatives

Further model validation was confirmed by comparing the basic model against both simpler and more complex alternatives.

Simpler model, without the prepower stroke state, does not work as well. Because of the very small contribution of the &Pdot;r(t) component, it was important to test whether the 2-state, 4-parameter model, which predicted that there would be a &Pdot;r(t) response component, was actually needed, or whether a simpler model would suffice. A model that discounts the contribution by prepower stroke state e0 and considers only contributions from the postpower stroke state ep is more parsimonious. Such a model predicts a response that is composed of a single component proportional only to Pr(t) and possesses just 1 state and 2 parameters as given by
&Dgr;<IT><A><AC>X</AC><AC>˙</AC></A></IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>) = −<FENCE><FR><NU><IT><A><AC>E</AC><AC>˙</AC></A></IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>)</NU><DE><IT>E<SUB>ep</SUB></IT>(<IT>t</IT>)</DE></FR> + <IT>g</IT></FENCE> &Dgr;<IT>X</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>) + &Dgr;<A><AC>V</AC><AC>˙</AC></A>(<IT>t</IT>) (24)
&Dgr;<A><AC>P</AC><AC>ˆ</AC></A><SUB><IT>fa</IT></SUB>(<IT>t</IT>) = [<IT>E</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>)]&Dgr;<IT>X<SUB>ep</SUB></IT>(<IT>t</IT>) (25)
The simple 1-state, 2-parameter (Eqs. 10 and 24-25) and the basic 2-state, 4-parameter (Eqs. 9-13) models were each fit to the data and compared to determine which one was the best based on reduction in the AIC and SC and on the incremental F-test as described above.

The basic 2-state, 4-parameter model was superior to the simple 1-state, 2-parameter model in fitting the data. The improvement in median R2 for data from the 14 hearts was on the order of 3% (from 0.951 to 0.978). Smaller AIC and SC were associated with the 2-state, 4-parameter model in all 14 hearts with a median reduction of 8.5% (range = -1.22 to -13.6%) and 8.5% (range = -1.22 to -13.6%), respectively. Furthermore, the incremental F-test generated an F-statistic that was significant at P < 0.01 in all 14 hearts. These are strong results for accepting the 2-state, 4-parameter model over the simpler 1-state, 2-parameter model. Another way of stating the outcome of the incremental F-test is that the precision of the estimates of parameters remain within an acceptable range in the 2-state, 4-parameter model. Precision of parameter estimates were quantified by using the standard error of parameter estimate (11) using the RSS and covariance information returned from the optimization algorithm. Over all 14 hearts, the median and range of the standard error of the estimates of X0, g, h, and d were (expressed as percentage of the respective parameter estimate) 0.19 (0.08-0.44), 0.38 (0.00-0.65), 1.88 (0.78-31.32), and 0.60 (0.20-26.93). The heart that generated the upper extremes in the ranges for h and d was an outlier among these 14 hearts because the heart with the next highest standard error of the estimate for these two parameters exhibited values of 8.48 and 8.60%, respectively. Except for estimates of h and d in one heart, the respective parameters were generally estimated with more than adequate precision.

Despite these results, certain caveats in representing the prepower stroke e0 state in the model must be made clear. In contrast to the good correspondence between A1/Delta V and its model-derived equivalent (Fig. 8), which is a response primarily attributable to the postpower stroke state, only a rough correspondence was obtained between A2/Delta V and its model-derived equivalent (Fig. 9), which is a response attributable to the prepower stroke state. The order of magnitude of values of A2/Delta V and its model-derived equivalent were the same, but dependencies of these quantities on f and Delta V were not the same: there was no dependence on f and strong dependence on Delta V for A2/Delta V, whereas there was weak dependence on f and no dependence on Delta V for the model-derived equivalent. Employing the procedures used in the derivation of Eq. 25, we show that an approximate model-derived equivalent of A2/Delta V is
<FR><NU>1</NU><DE><IT>hX</IT><SUB>0</SUB></DE></FR> <FENCE> <FR><NU>&Dgr;<IT>X</IT><SUB><IT>e0</IT></SUB>(<IT>f</IT>)</NU><DE>&Dgr;V(<IT>f</IT>)</DE></FR> </FENCE> (26)
where
<FENCE> <FR><NU>&Dgr;<IT>X</IT><SUB><IT>e0</IT></SUB>(<IT>f</IT>)</NU><DE>&Dgr;V(<IT>f</IT>)</DE></FR> </FENCE> = <FR><NU>2&pgr;<IT>f</IT></NU><DE><RAD><RCD>(2&pgr;<IT>f</IT>)<SUP>2</SUP> + (<IT>h</IT> + <IT>d</IT>)<SUP>2</SUP></RCD></RAD></DE></FR> (27)
Because median h + d was 2,555 s-1 (see Table 1), at 100 Hz
<FR><NU>1</NU><DE><IT>hX</IT><SUB>0</SUB></DE></FR> <FENCE> <FR><NU>&Dgr;<IT>X</IT><SUB><IT>e0</IT></SUB>(<IT>f</IT>)</NU><DE>&Dgr;V(<IT>f</IT>)</DE></FR> </FENCE>< 0.01 (28)
This model approximate value, which is on the order of values found for A2/Delta V, is more than two orders of magnitude less than that for A1/Delta V, implicating virtually no role for the prepower stroke state. However, the contribution of the prepower stroke state is not only determined by the relative values of A1/Delta V and A2/Delta V, but also because the amplitudes of components associated with pre- and postpower stroke states depend on the relative values of Pr(t) and &Pdot;r(t), which change with time over the cycle. A more valid assessment of contributions by the prepower stroke state is in the ratio of A2&Pdot;r(t)/A1Pr(t). This ratio varied with frequency, but representative values ranged between 1 and 5%. Even with these considerations, the contribution by the prepower stroke state is relatively small. Thus despite the statistical evidence suggesting good precision of the estimates in most cases and significant reduction in RSS, AIC, and SC, the prepower stroke state is contributing very little to the response and caution must be exercised in placing unqualified confidence in the estimates of associated parameters h and d.


View larger version (10K):
[in this window]
[in a new window]
 
Fig. 9.   bullet , Model-predicted equivalent of A2/Delta V for each of 20 volume perturbations. There is virtually no dependence of model-predicted A2/Delta V on perturbation amplitude and only a weak dependence on frequency. open circle , A2/Delta V from sinusoidal analysis. There is only a rough correspondence between model-predicted and sinusoidal A2/Delta V.

The reason for the relatively small dynamic contribution of the prepower stroke state in the 25- to 100-Hz range is the very high value of the characteristic frequency of responsible underlying processes h + d = 2,555 s-1. As a consequence, the prepower stroke state is very labile and, once a generator enters this state, it does not last long enough to undergo appreciable distortion when the volume is perturbed at frequencies in the 25- to 100-Hz frequency range. It would take frequencies in the 400-Hz (congruent 2,555/2pi ) range to generate appreciable distortion in the e0 state and make a sizable contribution to the response. Consequently, the e0 state makes only a small contribution to the response over the frequencies used in these studies, and the associated parameters are estimated with less precision than those associated with the postpower stroke state.

There were additional consequences to incorporating the prepower stroke state in the model. Including this state doubled the estimated g from a median of 69.0 s-1 in the 1-state, 2-parameter model to 115.6 s-1 in the 2-state, 4-parameter model. Also, including this state increased the robustness of the parameter estimates, making their values less sensitive to the manner in which the data were prepared for analysis, i.e., values of filter settings, complete removal of trends, etc. (results not shown). We concluded that although it was important to include contributions from the prepower stroke state to account for variability in the data and for achieving robustness in the estimates of the parameters, we could identify only the rough magnitude of the underlying dynamic processes, because it made only small contributions to the overall response at even the highest frequency used in these studies (100 Hz).

Addition of noncontractile series elastance did not improve model. An important assumption in the model was that only cross-bridges were contributing to active chamber properties, i.e., there was no contribution by noncontractile structures in parallel or in series with the cross-bridge force generators. Contributions by noncontractile parallel structures were identified from passive behavior and removed by previously described methods. Our previous findings (4) indicated that there was no important noncontractile series elasticity participating in LV mechanical behavior to rapid ramp changes in volume. However, there was a need to confirm that this was so for sinusoidal volume changes.

To test for the contributions by noncontractile series elastic elements, the basic 2-state, 4-parameter model was compared with an elaborated model containing series elastance and consisting of 3 states and 5 parameters. The third-state variable was required to account for internal shortening during Pr(t), and the fifth parameter was a noncontractile series elastance. Equations for the model with series elastance are given in the APPENDIX. In addition to the extra state and parameter, these equations differed from the basic model Eqs. 9-13 in that the two distortional variables Delta Xe0(t) and Delta Xep(t) were no longer uncoupled, i.e., the differential equations describing the derivatives of each of these variables depended on both variables. The model with series elastance was fit to the data and then compared with the basic model without series elastance using the AIC, SC, and the incremental F-test.

The improvement in fit when series elastance was added to the basic model was scarcely observable, the median R2 value was 0.978 both with and without the series elastance. Furthermore, the AIC was reduced in only 9 of 14 hearts, and the SC was reduced in only 8 of 14 hearts. The incremental F-statistic was significant in 8 of 14 hearts (P < 0.05). In the six hearts where the series elastic element was found not to be significant, the estimated value was >1010 mmHg/ml. In these six hearts, the model was equally well fit with any large value of series elastance, leading to considerable imprecision in parameter estimates and correspondingly insignificant reduction in RSS. In the eight hearts where series elastance was found to be significant, the average estimated value was 17,892 mmHg/ml. This compares with an average peak generator elastance [i.e., the maximal value of Ee0(t) + Eep(t)] of 411 mmHg/ml. Thus when it was found to be significant (in roughly half of all hearts), series elastance was 43 times greater than maximal generator elastance and, consequently, played only a minor role in chamber mechanical properties. Furthermore, the addition of the series elastance element did not appreciably change the estimated value of the other model parameters (Table 1). In accordance with previous findings (4), we concluded that there was no important role for noncontractile series elastance in these hearts and no compelling reason for complicating the basic model by incorporating noncontractile series elastance.

Incorporation of distortion-dependent cross-bridge detachment did not improve model. Another fundamental model assumption was that the kinetic rate constants were truly constants and did not depend on behavior such as distortion of the generators. There are several reasons why this assumption required testing. First, distortion-dependent, cross-bridge rate constants have been a central tenet in much of the theoretical work in muscle mechanics (14) and have been a major feature of several models to explain aspects of mechanical behavior (10, 15, 22, 23, 34). Second, enhanced cross-bridge detachment (i.e., elevated detachment rate constant) secondary to vibration-induced, cross-bridge distortion in cardiac muscle has been hypothesized to be the mechanism by which force was depressed (12, 18, 21, 36, 37) and relaxation period shortened (16) in cardiac muscle, and this depressive response was a prominent part of the response observed in these experiments (as in Fig. 1). Third, our results indicated a small but significant effect of Delta V on A2/delta V through the Delta V2 term in Eq. 23, which suggested amplitude-related nonlinearity such as distortion-dependent detachment.

To test whether distortion-dependent kinetic factors accounted for any of the observed behavior, we incorporated these effects into an expanded form of the basic model and compared the results with those from the basic 2-state, 4-parameter model. To introduce distortional effects, we considered that both distortional variables, Delta Xe0(t) and Delta Xep(t), vary with Delta V(t) according to Eqs. 11 and 12. Because the rate constant governing detachment of the prepower stroke state (h + d) was very large compared with both the frequency of perturbation and the detachment rate constant for the postpower stroke state (g), a Delta V(t) would induce a much greater Delta Xep(t) than Delta Xe0(t) and, because of this, it is reasonable to consider only the effects of Delta V-induced changes in Delta Xep(t) on g, the rate constant of cross-bridge detachment. We used two functional forms of distortion-dependent g, where g varied around its isometric (isovolumic) value (g0) according to either
<IT>g</IT> = <IT>g</IT>(<IT>t</IT>) = <IT>g</IT><SUB>0</SUB> + <IT>g</IT><SUB>1</SUB> <FENCE><FR><NU>&Dgr;<IT>X</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>)</NU><DE><IT>X</IT><SUB>0</SUB></DE></FR></FENCE> (29)
or
<IT>g</IT> = <IT>g</IT>(<IT>t</IT>) = <IT>g</IT><SUB>0</SUB> + <IT>g</IT><SUB>1</SUB> <FENCE><FR><NU>&Dgr;<IT>X<SUB>ep</SUB></IT>(<IT>t</IT>)</NU><DE><IT>X</IT><SUB>0</SUB></DE></FR></FENCE><SUP>2</SUP> (30)
These two formulations differ importantly as follows. In Eq. 29, when Delta Xep(t) is negative, g takes on values less than g0, i.e., reduced generator distortion below baseline would slow rate of detachment, whereas increased distortion would increase rate of detachment. This contrasts with the formulation of Eq. 30, where g is greater than g0 for both positive and negative values of Delta Xep(t), i.e., both reduced and increased distortion increase g and rate of detachment.

A fundamental difference between equations of the model with distortion-dependent g compared with those of the basic model (Eqs. 9-13) was that, with distortion-dependent g, it could no longer be assumed that there were no changes in elastance at the perturbation frequency. Accordingly, it was necessary to add a third-state variable to account for in-frequency variation in Eep(t). Equations for this formulation are given in the APPENDIX.

Comparison of competing model formulations was by fitting the model with each g formulation, then comparing results from Eq. 29 with those from the basic model using AIC, SC, and incremental F-test, and then comparing results from Eq. 30 with those from the basic model using the same criteria.

There was no systematic improvement in accounting for variation in the data (R2) in reducing the AIC and SC or in yielding a significant incremental F-statistic when either form of distortion-dependent detachment (Eqs. 29 and 30) was incorporated into the model (Table 1). Furthermore, the standard error of the estimate of the parameter associated with distortion dependence, in those cases where the parameter was found to be significant by the incremental F-test, was, on the average 5-20 times larger than the standard error of the other parameter estimates. Such lack of precision in the parameter estimate indicated that the parameter played a lesser role in accounting for variability in the data. Finally, the values of the other four parameters were virtually unaffected by whether distortion-dependent detachment was part of the model or not (Table 1). Thus the effect of including distortion dependence in the model was negligible even when statistically significant.

We concluded that the added complexity as a result of incorporating distortion-dependent g into the basic model, in either of the forms tested, was not worth the negligible increase in ability to explain the in-frequency component of the response to high-frequency volume perturbations.

From the ability of the basic 2-state, 4-parameter model to account for a wide variety of observed features and the inability of variants to bring about any significant improvement, we adopted it as the best of several options and as a proper representation of mechanisms responsible for the active part of the in-frequency response. This representation not only satisfies statistical requirements, but it is consistent with all of the major phenomenology seen in this analysis.

    OVERALL SUMMARY AND CONCLUSIONS
Top
Abstract
Introduction
Methods
Conclusions
Appendix
References

Dynamic Processes Between 25 and 100 Hz Are Driven by Derivative of Volume Rather Than Volume Itself

Evidence indicating that LV dynamic processes in the 25- to 100-Hz range are responding to the derivative of volume change comes from the model-independent observations that the phase of &Dgr;P<SUB><IT>fa</IT></SUB>(t) pressure leads Delta V(t) in frequency-dependent manner in this frequency range and that the A1/Delta V amplitude ratio increases with frequency rather than decreases. This is strong evidence that the underlying dynamic processes are obeying behavioral patterns similar to those represented by the distortional Eqs. 11 and 12 of the model. Thus when we employ the model equations to predict the observations, they reproduce these observations with a high degree of fidelity. We argue that this distortion actually resides in the cross-bridges because we can find no evidence in this work and in previous work (4) that there is any effective series elastance in noncontractile structures in these isolated hearts.

Nonlinearities Are Not an Important Part of Response to Small-Amplitude Volume Sinusoids

We found no convincing evidence for important nonlinearity in these responses. First, when we plotted the A1/Delta V amplitude ratio in Fig. 4 and conducted the regression analysis, there was only a slight indication that any nonlinearities were operative. Furthermore, when we fit the data with the model and examined the correspondence between model-fitted and measured pressure as in Fig. 6, there was no indication for unrepresented nonlinearity; unrepresented nonlinearity would have exhibited itself by imposing some crescent shape to the oblong cluster of points in Fig. 6. Absolutely no evidence of crescent shape to this cluster was detectable. Finally, when we introduced the most likely candidate for nonlinearity into the model in the form of distortion-dependent, cross-bridge detachment, we were unable to detect any improvement.

Parameters From Current Sinusoidal Technique Agree With Those From Previous Studies and With Those Obtained From Isolated Muscle

It is instructive to compare results for g and X0 obtained in this sinusoidal analysis with estimates of these parameters from earlier studies done in this laboratory. Using ferret hearts at 30°C, ramp-volume withdrawals at the time of peak isovolumic pressure, and a model equivalent to the 1-state, 2-parameter model of the current study, we (5, 30) found values of g of ~70-80 s-1 and of X0 of 0.31 ml. These compare with values of 115 s-1 for g and 0.29 ml for X0 in the current study. Thus g, as estimated here, was 1.5 times greater than our previous estimates with other techniques, whereas X0 was not different. [Interestingly, the value of g estimated from the 1-state, 2-parameter model in the current sinusoidal study (69 s-1) is the same as the value we previously obtained using essentially this model in the ramp-withdrawal studies.] In results reported elsewhere (unpublished observations) that were obtained using the sinusoidal technique in beating rabbit hearts at a lower temperature of 25°C, we found average g from the 2-state, 4-parameter model to be 49 s-1, which is the same as we previously obtained at 25°C from ferret hearts using ramp-volume withdrawals and the 1-state, 2-parameter model (28). Furthermore, the 49 s-1 value at 25°C in beating rabbit hearts from the current sinusoidal method is the same as a value of 52 s-1 estimated in constantly activated (Ba2+), nonbeating rabbit hearts at 25°C, using responses over a dense spectrum of sinusoidal frequencies between 0.1 and 30 Hz (6). In summary, values of g estimated in the current study at 30°C are higher than values obtained previously with the ramp-withdrawal technique at 30°C, but results with the current sinusoidal technique at 25°C in beating hearts are identical with previous results with the ramp-withdrawal technique in beating hearts and sinusoidal techniques in constantly activated hearts.

In earlier studies with the ramp-withdrawal technique (28), a Q10 for g of 2.15 was obtained. When the 25°C value of g reported elsewhere (unpublished observations) is used with the 30°C value reported here to calculate a Q10 for g from the sinusoidal technique, a value close to 5 is obtained. Although this value is far above our previous estimate, it compares favorably with a Q10 of 4.6 for maximal velocity of shortening (which is proportional to g) in membrane-intact rat trabecular muscle (9). At this time we do not know whether the 2-state, 4-parameter model detects features that respond more strongly to temperature than do features of the 1-state, 2-parameter model. Resolution of these contrasting results and the effects of temperature will have to await the conclusion of on-going studies.

Results obtained here in isolated hearts may also be compared with results from isolated muscle studies by Kawai and co-workers (17, 24, 39). These workers used skinned, maximally activated, papillary and trabecular muscles from ferret (17, 24) and porcine (39) hearts and conducted sinusoidal perturbation experiments at 20°C. Our g and h correspond with their rate constants for the detachment step, 2pi c, and the energy transduction step, 2pi b, respectively. Because 2pi c and 2pi b change dramatically with ATP, ADP, and Pi, which were varied in their experiments, the proper comparison is with values of these constants obtained when fiber-bathing solution concentrations were comparable to a membrane-intact condition. Reported values of comparable 2pi c and 2pi b were approximately 30 and 8.6 s-1 (24) or 37 and 15 s-1 (17), respectively, in ferret heart muscle; and 23 and 15 s-1, respectively, in porcine heart muscle (39). Despite differences in preparations, species, interpretive models, and analysis, there are remarkable similarities between our findings in the beating, isolated heart and those of Kawai's group in skinned, constantly activated, isolated muscle. First, the rate constant of the energy transduction step was on the same order of that of the detachment step in Kawai's studies (2pi b/2pi c = 1/3 to 2/3) as it was in our study (h/g = 2/3). Second, given the potential high Q10 for these processes, it can be argued that the isolated muscle values found at 20°C are in good agreement with those we found in the isolated heart at 30°C. Again, future studies at equivalent temperatures and with the same species are needed to establish how closely the estimates from the isolated heart agree with those from isolated muscle.

Cross-Bridge Mechanisms Can be Observed in Whole Heart Behavior

The most important outcome of this study is that aspects of cross-bridge kinetic behavior may be observed from pressure responses to volume perturbations in the isolated heart. This is based on evidence for the validity of the proposed cross-bridge model. Model validity was demonstrated from concordance between the sinusoidal analysis and model predictions, from how well the model fit the data, and from inability of elaborated models to improve representation of the data.

No model is perfect, including this one. Particularly, it is troubling that one of the model components, the component associated with the prepower stroke state, contributes so very little to the magnitude of the response. Thus a fair amount of caution must be exercised in relying on the accuracy of estimated parameter values associated with this state when estimation is made from responses to sinusoidal frequencies no higher than 100 Hz. We predict that frequencies of 400 Hz or more would be required to precisely estimate parameters associated with the prepower stroke state. More work is needed to secure greater confidence in this parameter estimate. However, given the overall good performance of the model and agreement of our estimated parameters with equivalent parameters found by others in isolated muscle, there is every reason to accept its general form and basic premises and to proceed with additional work for its verification and application.

In summary, pressure responses to small-amplitude, high-frequency volume perturbations in the isolated heart provide dynamic information relative to cross-bridge kinetics. The use of a model to analyze these responses gives estimates of specific kinetic parameters relative to the power stroke and cross-bridge detachment. The application of these methods to isolated heart studies promises to be a powerful experimental tool for elucidating mechanisms involved in changes in LV contractile function in response to a variety of interventions, including inotropic agents, ischemia, and chronic pathophysiological states such as cardiomyopathy.

    APPENDIX
Top
Abstract
Introduction
Methods
Conclusions
Appendix
References

Linear Transformation Between Incremental Wall Stress and Strain and Incremental Chamber Pressure and Volume

To justify the assumption that incremental linear pressure-volume elastance of the LV chamber can coexist with linear stress-strain muscle elasticity in the LV wall, we need to consider the relationship between representative lineal dimensions in the wall and chamber volume and between representative wall stress and chamber pressure. If changes in chamber volume linearly transform into changes in wall lineal dimensions and if changes in chamber pressure linearly transform into changes in wall stress, then incremental volumetric elastance can be taken as linearly related to incremental wall elasticity.

When it is assumed that the LV is a thick-walled sphere with a reference volume of VBL, then it can be shown that a change in the half-mass wall circumference (Delta lm/2) is related to a change in chamber volume (Delta V) according to
&Dgr;<IT>l</IT><SUB><IT>m</IT>/2</SUB> = <RAD><RCD>6&pgr;<SUP>2</SUP></RCD><RDX>3</RDX></RAD> <FENCE><RAD><RCD>&Dgr;V + V<SUB>BL</SUB> + <FR><NU>V<SUB>w</SUB></NU><DE>2</DE></FR></RCD><RDX>3</RDX></RAD> − <RAD><RCD>V<SUB>BL</SUB> + <FR><NU>V<SUB>w</SUB></NU><DE>2</DE></FR></RCD><RDX>3</RDX></RAD> </FENCE> (A1)
For a representative ventricle of VBL = 2 ml and VW = 6 ml, a Delta V of ±0.015 VBL will produce a Delta lm/2 of +0.1996% and -0.2004% of the reference lm/2. With the equivalent positive and negative changes in Delta lm/2 for equal positive and negative changes in Delta V, there is, for all intents and purposes, a linear relationship between Delta lm/2 and Delta V over this range of Delta V. Similarly, for this same ventricle the relationship between average wall stress and chamber pressure is given by
&sfgr; = <FENCE><FR><NU>1</NU><DE><FENCE><FR><NU>V<SUB>w</SUB></NU><DE>&Dgr;V + V<SUB>BL</SUB></DE></FR> + 1</FENCE><SUP>2/3</SUP> − 1</DE></FR></FENCE> P (A2)
In this relationship, pressure transforms approximately linearly into stress if the value of the coefficient on P is not changed appreciably with Delta V. For Delta V equalling ±0.015 VBL in these representative hearts, the coefficient on P varies between ±1.22% its baseline value of 0.658. This small change in the coefficient justifies the approximation of a constant coefficient and a linear transformation between P and sigma  over the range in volume given by Delta V.

Relationship Between Elastance and Measured Pressure

By definition, chamber pressure equals elastance times volumetric distortion. Invoking the parallel generator assumption and the model of Fig. 3
P(<IT>t</IT>) = <IT>E<SUB>e0</SUB></IT>(<IT>t</IT>)<IT>X<SUB>e0</SUB></IT>(<IT>t</IT>) + <IT>E<SUB>ep</SUB></IT>(<IT>t</IT>)<IT>X</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>) (A3)
where Xe0(t) and Xep(t) are the respective average volumetric distortions of Ne0(t) and Nep(t) generators.

Consider that elastances of the perturbed beat consist of isovolumic elastance plus changes in elastance induced by the perturbation. Furthermore, distortion consists of isovolumic distortion plus perturbational-induced in distortion. Thus
P(<IT>t</IT>) = [<IT>E</IT><SUP>iso</SUP><SUB><IT>e0</IT></SUB>(<IT>t</IT>) + &Dgr;<IT>E</IT><SUB><IT>e0</IT></SUB>(<IT>t</IT>)][<IT>X</IT><SUP>iso</SUP><SUB><IT>e0</IT></SUB> + &Dgr;<IT>X</IT><SUB><IT>e0</IT></SUB>(<IT>t</IT>)] 
+ [<IT>E</IT><SUP>iso</SUP><SUB><IT>ep</IT></SUB>(<IT>t</IT>) + &Dgr;<IT>E</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>)][<IT>X</IT><SUP>iso</SUP><SUB><IT>ep</IT></SUB> + &Dgr;<IT>X</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>) (A4)
where isovolumic conditions are indicated with the iso superscript and changes induced by the sinusoidal perturbation are indicated with the Delta . Because Xisoe0 = 0, Eq. A4 can be rewritten as
P(<IT>t</IT>) = {[<IT>E</IT><SUP>iso</SUP><SUB><IT>ep</IT></SUB>(<IT>t</IT>) + &Dgr;<IT>E</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>)]<IT>X</IT><SUP>iso</SUP><SUB><IT>ep</IT></SUB>}
+ {[<IT>E</IT><SUP>iso</SUP><SUB><IT>e0</IT></SUB>(<IT>t</IT>) + &Dgr;<IT>E<SUB>e0</SUB></IT>(<IT>t</IT>)]&Dgr;<IT>X</IT><SUB><IT>e0</IT></SUB>(<IT>t</IT>)}
+ {[<IT>E</IT><SUP>iso</SUP><SUB><IT>ep</IT></SUB>(<IT>t</IT>) + &Dgr;<IT>E</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>)]&Dgr;<IT>X</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>)} (A5)
where the first term in brackets on the right-hand side is the pressure around which the distortional response is occurring
P<SUB>r</SUB>(<IT>t</IT>) = [<IT>E</IT><SUP>iso</SUP><SUB><IT>ep</IT></SUB>(<IT>t</IT>) + &Dgr;<IT>E</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>)]<IT>X</IT><SUP>iso</SUP><SUB><IT>ep</IT></SUB> (A6)
and the second term in brackets is the distortional or in-frequency response
&Dgr;P<SUB><IT>f</IT></SUB>(<IT>t</IT>) = {[<IT>E</IT><SUP>iso</SUP><SUB><IT>e0</IT></SUB>(<IT>t</IT>) + &Dgr;<IT>E</IT><SUB><IT>e0</IT></SUB>(<IT>t</IT>)]&Dgr;<IT>X</IT><SUB><IT>e0</IT></SUB>(<IT>t</IT>)}
+ {[<IT>E</IT><SUP>iso</SUP><SUB><IT>ep</IT></SUB>(<IT>t</IT>) + &Dgr;<IT>E</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>)]&Dgr;<IT>X</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>)} (A7)
Thus
P(<IT>t</IT>) = P<SUB>r</SUB>(<IT>t</IT>) + &Dgr;P<SUB><IT>f</IT></SUB> (<IT>t</IT>)
Comparing Eqs. A5 and A3, we make the following assignments
<IT>E</IT><SUB><IT>e0</IT></SUB>(<IT>t</IT>) = <IT>E</IT><SUP>iso</SUP><SUB><IT>e0</IT></SUB>(<IT>t</IT>) + &Dgr;<IT>E</IT><SUB><IT>e0</IT></SUB>(<IT>t</IT>)
<IT>E</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>) = <IT>E</IT><SUP>iso</SUP><SUB><IT>ep</IT></SUB>(<IT>t</IT>) + &Dgr;<IT>E</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>)
<IT>X</IT><SUB>0</SUB> = <IT>X</IT><SUP>iso</SUP><SUB><IT>ep</IT></SUB> (A8)
Then from Eq. A6, it can be written that
<IT>E</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>) = <FR><NU>P<SUB>r</SUB>(<IT>t</IT>)</NU><DE><IT>X</IT><SUB>0</SUB></DE></FR> (A9)
An expression for Ee0(t) can be developed by considering the differential equation for the ep state. From Fig. 3
<IT><A><AC>N</AC><AC>˙</AC></A></IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>) = <IT>hN<SUB>e0</SUB></IT>(<IT>t</IT>) − <IT>gN<SUB>ep</SUB></IT>(<IT>t</IT>) (A10)
which, after substituting Eqs. 16 and 17 in the text, can be rewritten in terms elastances as
<IT><A><AC>E</AC><AC>˙</AC></A></IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>) = <IT>hE<SUB>e0</SUB></IT>(<IT>t</IT>) − <IT>gE<SUB>ep</SUB></IT>(<IT>t</IT>) (A11)
Solving Eq. A11 for Ee0(t)
<IT>E<SUB>e0</SUB></IT>(<IT>t</IT>) = <FR><NU><IT><A><AC>E</AC><AC>˙</AC></A></IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>) + <IT>gE<SUB>ep</SUB></IT>(<IT>t</IT>)</NU><DE><IT>h</IT></DE></FR> (A12)
Substituting for Eep(t) from Eq. A9 gives the following expression for Ee0(t)
<IT>E<SUB>e0</SUB></IT>(<IT>t</IT>) = <FR><NU><A><AC>P</AC><AC>˙</AC></A><SUB>r</SUB>(<IT>t</IT>) + <IT>gP</IT><SUB>r</SUB>(<IT>t</IT>)</NU><DE><IT>hX</IT><SUB><IT>0</IT></SUB></DE></FR> (A13)

Differential Equations for Distortional Variables

The differential equations for Delta Xep(t) may be derived by considering the summed distortion among all generators within a given state. The total distortion among all ep generators is given by
<IT>Z<SUB>ep</SUB></IT>(<IT>t</IT>) = <IT>N</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>)<IT>X<SUB>ep</SUB></IT>(<IT>t</IT>) (A14)
where Xep(t) is the average distortion among the Nep(t) generators.

At some t + Delta t, Zep(t + Delta t) can be written as Zep(t + Delta t) = Zep(t) + (added distortion due to change in volume over Delta t) + (added distortion due to formation of new units with baseline distortion over Delta t- (lost distortion due to detachment of distorted units over Delta t), where (added distortion due to change in length over Delta t) = (externally imposed Delta V) · (no. of ep generators existent at t) = Delta V Nep(t), and where (added distortion due to formation of new units with X0 distortion over Delta t) = [h · Ne0(t)]Delta t(X0 + Delta V/2), and (lost distortion due to detachment of distorted units over Delta t) = [g · Nep(t)]Delta t[Xep(t)] Therefore
<IT>Z<SUB>ep</SUB></IT>(<IT>t</IT> + &Dgr;<IT>t</IT>) = <IT>Z</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>) + &Dgr;V<IT>N</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>) 
+ [<IT>hN<SUB>e0</SUB></IT>(<IT>t</IT>)]&Dgr;<IT>t</IT> <FENCE><IT>X</IT><SUB>0</SUB> + <FR><NU>&Dgr;V</NU><DE>2</DE></FR></FENCE> − [<IT>gN<SUB>ep</SUB></IT>(<IT>t</IT>)]&Dgr;<IT>tX</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>) (A15)
and
<FR><NU><IT>Z<SUB>ep</SUB></IT>(<IT>t</IT> + &Dgr;<IT>t</IT>) − <IT>Z</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>)</NU><DE>&Dgr;<IT>t</IT></DE></FR> = <IT>N</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>) <FR><NU>&Dgr;V</NU><DE>&Dgr;<IT>t</IT></DE></FR>
+ [<IT>hN</IT><SUB><IT>e0</IT></SUB>(<IT>t</IT>)] <FENCE><IT>X</IT><SUB>0</SUB> + <FR><NU>&Dgr;V</NU><DE>2</DE></FR></FENCE> − [<IT>gN</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>)]<IT>X</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>)  (A16)
Taking the limit as Delta t right-arrow 0 
<IT><A><AC>Z</AC><AC>˙</AC></A></IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>) = <IT>N<SUB>ep</SUB></IT>(<IT>t</IT>)<IT><A><AC>V</AC><AC>˙</AC></A></IT>(<IT>t</IT>) + [<IT>hN</IT><SUB><IT>e0</IT></SUB>(<IT>t</IT>)][<IT>X</IT><SUB>0</SUB>] − <IT>gN<SUB>ep</SUB></IT>(<IT>t</IT>)<IT>X</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>) (A17)
Note that differentiation of Eq. A10 yields
<IT><A><AC>Z</AC><AC>˙</AC></A><SUB>ep</SUB></IT>(<IT>t</IT>) = <IT><A><AC>N</AC><AC>˙</AC></A></IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>)<IT>X</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>) + <IT>N<SUB>ep</SUB></IT>(<IT>t</IT>)<IT><A><AC>X</AC><AC>˙</AC></A></IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>) (A18)
Equating Eqs. A17 and A18, making appropriate substitutions for &Ndot;ep(t) from Eq. A10, and solving for &Xdot;ep(t) gives
<IT><A><AC>X</AC><AC>˙</AC></A><SUB>ep</SUB></IT>(<IT>t</IT>) = −<IT>h</IT> <FR><NU><IT>N<SUB>e0</SUB></IT>(<IT>t</IT>)</NU><DE><IT>N</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>)</DE></FR> [<IT>X</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>) − <IT>X</IT><SUB>0</SUB>] + <A><AC>V</AC><AC>˙</AC></A>(<IT>t</IT>) (A19)
To remove the dependence on Ne0(t), we rearrange Eq. A10 to get
<IT>h</IT> <FR><NU><IT>N<SUB>e0</SUB></IT>(<IT>t</IT>)</NU><DE><IT>N</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>)</DE></FR> = <FR><NU><IT><A><AC>N</AC><AC>˙</AC></A></IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>)</NU><DE><IT>N</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>)</DE></FR> + <IT>g</IT> (A20)
Substituting Eq. A20 into Eq. A19 and writing in terms of the incremental variable, Delta Xep(t) = Xep(t) - X0
&Dgr;<IT><A><AC>X</AC><AC>˙</AC></A><SUB>ep</SUB></IT>(<IT>t</IT>) = <FENCE><FR><NU><IT><A><AC>N</AC><AC>˙</AC></A><SUB>ep</SUB></IT>(<IT>t</IT>)</NU><DE><IT>N</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>)</DE></FR> + <IT>g</IT></FENCE> &Dgr;<IT>X</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>) + &Dgr;<A><AC>V</AC><AC>˙</AC></A>(<IT>t</IT>) (A21)
Finally, substituting for Nep(t) its equivalent in terms of Eep(t), as given by Eq. 8 in the text, gives the final form of the differential equation as given by Eqs. 12 and 11 in the text.

Equations When There is Noncontractile Series Elastance

When there is noncontractile series elastance, there is internal shortening during isovolumic or reference beats. Under assumptions made in the text
<IT>E</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>) = <FR><NU>P<SUB>r</SUB>(<IT>t</IT>)</NU><DE><IT>X</IT><SUB><IT>ep</IT><SUB>r</SUB></SUB>(<IT>t</IT>)</DE></FR> (A22)
where Xepr(t) is no longer X0 as was true when the series elastance was considered to be infinitely stiff. Given a measured Pr(t), it can be shown that when there is internal shortening during a reference beat, generator distortion as a result of internal shortening is given by
<IT><A><AC>X</AC><AC>˙</AC></A></IT><SUB><IT>ep</IT><SUB>r</SUB></SUB>(<IT>t</IT>) = − <FENCE><FENCE><FR><NU><IT><A><AC>P</AC><AC>˙</AC></A></IT><SUB>r</SUB>(<IT>t</IT>)</NU><DE><IT>P</IT><SUB>r</SUB>(<IT>t</IT>)</DE></FR> + <IT>g</IT></FENCE> [<IT>X</IT><SUB><IT>ep</IT><SUB>r</SUB></SUB>(<IT>t</IT>) − <IT>X</IT><SUB>0</SUB>] + <FR><NU><IT><A><AC>P</AC><AC>˙</AC></A></IT><SUB>r</SUB>(<IT>t</IT>)</NU><DE><IT>E</IT><SUB>SE</SUB></DE></FR></FENCE> <FR><NU><IT>X</IT><SUB><IT>ep</IT><SUB>r</SUB></SUB>(<IT>t</IT>)</NU><DE><IT>X</IT><SUB>0</SUB></DE></FR> (A23)
where ESE is the elastance of series-coupled, noncontractile elements. Solving Eq. A23 for Xepr(t) allows calculation of Eep(t) using Eq. A22. Furthermore, the presence of a compliant series elastance changes the state equations as follows
&Dgr;<IT><A><AC>X</AC><AC>˙</AC></A></IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>) = − <FENCE><FR><NU><IT><A><AC>E</AC><AC>˙</AC></A></IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>)</NU><DE><IT>E</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>)</DE></FR> + <IT>g</IT></FENCE> &Dgr; <IT>X</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>) + &Dgr;<A><AC>V</AC><AC>˙</AC></A>(<IT>t</IT>)
− <FR><NU>1</NU><DE><IT>E</IT><SUB><IT>SE</IT></SUB></DE></FR> <FR><NU>d</NU><DE>d<IT>t</IT></DE></FR> [<IT>E</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>)&Dgr; <IT>X</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>) + <IT>E</IT><SUB><IT>e0</IT></SUB>(<IT>t</IT>)&Dgr; <IT>X</IT><SUB><IT>e0</IT></SUB>(<IT>t</IT>)] (A24)
&Dgr;<IT><A><AC>X</AC><AC>˙</AC></A></IT><SUB><IT>e0</IT></SUB>(<IT>t</IT>) = − <FENCE><FR><NU><IT><A><AC>E</AC><AC>˙</AC></A></IT><SUB><IT>e0</IT></SUB>(<IT>t</IT>)</NU><DE><IT>E</IT><SUB><IT>e0</IT></SUB>(<IT>t</IT>)</DE></FR> + <IT>h</IT> + <IT>d</IT></FENCE> &Dgr; <IT>X</IT><SUB><IT>e0</IT></SUB>(<IT>t</IT>) + &Dgr;<A><AC>V</AC><AC>˙</AC></A>(<IT>t</IT>)
− <FR><NU>1</NU><DE><IT>E</IT><SUB><IT>SE</IT></SUB></DE></FR> <FR><NU>d</NU><DE>d<IT>t</IT></DE></FR> [<IT>E</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>)&Dgr; <IT>X</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>) + <IT>E</IT><SUB><IT>e0</IT></SUB>(<IT>t</IT>)&Dgr; <IT>X</IT><SUB><IT>e0</IT></SUB>(<IT>t</IT>)] (A25)
After terms were collected and expressed explicitly in terms of the state variables, it is seen that these are no longer independent but are now coupled equations.

Equations When There is Distortion-Dependent "g"

In addition to incorporating the functional form of g into the differential equations of the basic model, it was necessary to expand these equations because, as g changed its value during the period of a volume sinusoid according to the distortions that were imposed, Nep(t) and, consequently, Eep(t) also changed their values subsequent to these variations in g. These in-frequency variations in Eep(t) were treated as follows. Let
<IT>g</IT> = <IT>g</IT><SUB>0</SUB> + &Dgr;<IT>g</IT>(<IT>t</IT>) and
<IT>E</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>) = <IT>E</IT><SUB><IT>ep</IT><SUB>r</SUB></SUB>(<IT>t</IT>) + &Dgr;<IT>E</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>) (A26)
where g0 and Ėepr(t) are the values around which variations occurred and Delta g(t) and Delta Eep(t) are the variations themselves. From Eq. A11
<IT><A><AC>E</AC><AC>˙</AC></A><SUB>ep</SUB></IT>(<IT>t</IT>) = <IT>hE</IT><SUB><IT>e0</IT></SUB>(<IT>t</IT>) − <IT>gE</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>)
= <IT>hE</IT><SUB><IT>e0</IT></SUB>(<IT>t</IT>) − [<IT>g</IT><SUB>0</SUB> + &Dgr;<IT>g</IT>(<IT>t</IT>)]<IT>E</IT><SUB><IT>ep</IT><SUB><IT>r</IT></SUB></SUB>(<IT>t</IT>) + &Dgr;<IT>E</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>)]
= <IT>hE</IT><SUB><IT>e0</IT></SUB>(<IT>t</IT>) − <IT>g</IT><SUB>0</SUB> <IT>E</IT><SUB><IT>ep</IT><SUB>r</SUB></SUB>(<IT>t</IT>)
− [<IT>g</IT><SUB>0</SUB> + &Dgr;<IT>g</IT>(<IT>t</IT>)]&Dgr;<IT>E<SUB>ep</SUB></IT>(<IT>t</IT>) − &Dgr;<IT>gE</IT><SUB><IT>ep</IT><SUB>r</SUB></SUB>(<IT>t</IT>)
= <IT><A><AC>E</AC><AC>˙</AC></A></IT><SUB><IT>ep</IT><SUB>r</SUB></SUB>(<IT>t</IT>) + &Dgr;<IT><A><AC>E</AC><AC>˙</AC></A></IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>) (A27)
where
<IT><A><AC>E</AC><AC>˙</AC></A></IT><SUB><IT>ep</IT><SUB>r</SUB></SUB>(<IT>t</IT>) = <IT>hE</IT><SUB><IT>e0</IT></SUB>(<IT>t</IT>) − <IT>g</IT><SUB>0</SUB><IT>E</IT><SUB><IT>ep</IT><SUB>r</SUB></SUB>(<IT>t</IT>) (A28)
and
&Dgr;<IT><A><AC>E</AC><AC>˙</AC></A></IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>) = −[<IT>g</IT><SUB>0</SUB> + &Dgr;<IT>g</IT>(<IT>t</IT>)]&Dgr;<IT>E</IT><SUB><IT>ep</IT></SUB>(<IT>t</IT>) − &Dgr;<IT>gE</IT><SUB><IT>ep</IT><SUB>r</SUB></SUB>(<IT>t</IT>) (A29)
We do not need Eq. A28 because Eepr(t) can be calculated from Eepr(t) Pr(t)/X0. Substituting Eepr(t) into Eq. A29 and then solving that equation numerically allows the calculation of &Dgr;<IT>E</IT><SUB><IT>ep</IT></SUB>(t). Thus all elements of Eep(t) are determined and the complete model solution with the distortion-dependent g may be conducted.

    ACKNOWLEDGEMENTS

This work was support by National Heart, Lung, and Blood Institute Grant HL-21462.

    FOOTNOTES

Address for reprint requests: K. Campbell, Depts. of VCAPP & Biological Systems Engineering, Washington State Univ., Pullman, WA 99164.

Received 14 January 1997; accepted in final form 2 July 1997.

    REFERENCES
Top
Abstract
Introduction
Methods
Conclusions
Appendix
References

1.   Boom, H. B., and H. Wijkstra. The step response of the left ventricular pressure to ejection flow: a system oriented approach. Ann. Biomed. Eng. 20: 99-126, 1992[Medline].

2.   Campbell, K. B. Rate constant of muscle force re-development reflects cooperative activation as well as crossbridge kinetics. Biophys. J. 72: 254-262, 1997[Medline].

3.   Campbell, K. B., L. W. Campbell, J. E. Pinto, and T. D. Burton. Contractile-based model interpretation of pressure-volume dynamics in the constantly activated (Ba2+) isolated heart. Ann. Biomed. Eng. 22: 550-567, 1994[Medline].

4.   Campbell, K. B., R. D. Kirkpatrick, A. H. Tobias, H. Taheri, and S. G. Shroff. Series-coupled passive elastance is functionally unimportant in the isolated heart. Cardiovasc. Res. 28: 242-251, 1994[Abstract/Free Full Text].

5.   Campbell, K. B., S. G. Shroff, and R. D. Kirkpatrick. Short time-scale LV systolic dynamics: evidence for a common mechanism in both LV chamber and heart-muscle mechanics. Circ. Res. 68: 1532-1548, 1991[Abstract/Free Full Text].

6.   Campbell, K., B., H. Taheri, R. D. Kirkpatrick, and B. K. Slinker. Single perturbed beat vs. steady-state beats for assessing systolic mechanical function in the isolated heart. Am. J. Physiol. 262 (Heart Circ. Physiol. 31): H1631-H1639, 1992[Abstract/Free Full Text].

7.   Campbell, K. B., H. Taheri, R. D. Kirkpatrick, T. Burton, and W. C. Hunter. Similarities between dynamic elastance of left ventricular chamber and papillary muscle of rabbit heart. Am. J. Physiol. 264 (Heart Circ. Physiol. 33): H1926-H1941, 1993[Abstract/Free Full Text].

9.   De Tombe, P. P., and H. E. D. J. ter Keurs. Force and velocity of sarcomere shortening in trabeculae from rat heart. Effects of temperature. Circ. Res. 66: 1239-1254, 1990[Abstract/Free Full Text].

10.   Eisenberg, E., T. L. Hill, and Y. Chen. Cross-bridge model of muscle contraction. Biophys. J. 29: 195-227, 1980[Medline].

11.   Glantz, S. A., and B. K. Slinker. Primer of Applied Regression and Analysis of Variance. McGraw-Hill, 1990.

12.   Honda, H., Y. Koiwa, and T. Takashima. Mathematical model of the effects of mechanical vibration on crossbridge kinetics in cardiac muscle. Jpn. Circ. J. 58: 416-425, 1994[Medline].

13.   Hunter, W. C., J. S. Janicki, K. T. Weber, and A. Noordergraaf. Systolic mechanical properties of the left ventricle: effects of volume and contractile state. Circ. Res. 52: 319-327, 1989[Abstract/Free Full Text].

14.   Huxley, A. F. Muscle structure and theories of contraction. Prog. Biophys. Mol. Biol. 7: 257-318, 1957.

15.   Huxley, A. F., and R. M. Simmons. Proposed mechanism of force generation in striated muscle. Nature 233: 533-538, 1971[Medline].

16.   Janssen, P. M. L., H. Honda, Y. Koiwa, and K. Shirato. The effect of of diastolic vibration on the relaxation of rat papillary muscle. Cardiovasc. Res. 32: 334-350, 1996.

17.   Kawai, M., Y. Saeki, and Y. Zhao. Crossbridge scheme and the kinetic constants of elementary steps deduced from chemically skinned papillary and trabecular muscles of the ferret. Circ. Res. 73: 35-50, 1993[Abstract].

18.   Koiwa, Y., N. Hoshi, T. Ohyama, T. Takagi, J. Kikuchi, H. Honda, and T. Takashima. Evaluation of left ventricular function by application of external minute vibration. Front. Med. Biol. Eng. 2: 207-210, 1990[Medline].

19.   Kirkpatrick, R. D., K. B. Campbell, D. L. Bell, and H. Taheri. Method for studying arterial wave-transmission effects on left ventricular function. Am. J. Physiol. 260 (Heart Circ. Physiol. 29): H1003-H1012, 1991[Abstract/Free Full Text].

20.   Landaw, E. M., and J. J. DiStefano. Multiexponential, multicompartmental, and noncompartmental modeling. I. Data analysis and statistical considerations. Am. J. Physiol. 246 (Regulatory Integrative Comp. Physiol. 15): R665-R677, 1984.

21.   Nishioka, T., Y. Goto, K. Hata, T. Takasago, A. Saeki, T. W. Taylor, and H. Suga. Mechanoenergetics of negative inotropism of ventricular wall vibration in dog heart. Am. J. Physiol. 270 (Heart Circ. Physiol. 39): H583-H593, 1996[Abstract/Free Full Text].

22.   Pate, E., and R. Cooke. A model of crossbridge action: the effects of ATP, ADP and Pi. J. Muscle Res. Cell Motil. 10: 181-196, 1989[Medline].

23.   Podolski, R. J., and A. C. Nolan. Muscle contraction transients, cross-bridge kinetics, and the Fenn effect. Cold Spring Harb. Symp. Quant. Biol. 37: 661-668, 1974.

24.   Saeki, Y., M. Kawai, and Y. Zhao. Comparison of cross-bridge dynamics between intact and skinned myocardium from ferret right ventricles. Circ. Res. 68: 772-781, 1991[Abstract/Free Full Text].

25.   Schiereck, P., and H. B. K. Boom. Left ventricular active stiffness: dependency on time and inotropic state. Pflügers Arch. 374: 135-143, 1978[Medline].

26.   Schiereck, P., and H. B. K. Boom. Left ventricular force-velocity relations measured from quick volume changes. Pflügers Arch. 379: 251-258, 1979[Medline].

27.   Shroff, S. G., K. B. Campbell, and R. D. Kirkpatrick. Short time-scale LV systolic dynamics: pressure vs. flow clamps and effects of activation. Am. J. Physiol. 264 (Heart Circ. Physiol. 33): H946-H959, 1993[Abstract/Free Full Text].

28.   Shroff, S. G., K. B. Campbell, D. E. Miller, R. D. Kirkpatrick, and H. Taheri. Effect of temperature on short time-scale left ventricular contractile dynamics (Abstract). Circulation 86: I553, 1992.

29.   Shroff, S. G., J. S. Janicki, and K. T. Weber. Left ventricular systolic dynamics in terms of its chamber mechanical properties. Am. J. Physiol. 245 (Heart Circ. Physiol. 14): H110-H124, 1983.

30.   Shroff, S. G., J. S. Janicki, and K. T. Weber. Evidence of quantitation of left ventricle resistance. Am. J. Physiol. 249 (Heart Circ. Physiol. 18): H358-H370, 1985.

31.   Suga, H., L. Demer, and K. Sagawa. Determinants of instantaneous pressure in canine left ventricle. Circ. Res. 46: 314-322, 1980.

32.   Templeton, G. H., and L. R. Nardizzi. Elastic and viscous stiffness of the canine left ventricle. J. Appl. Physiol. 36: 123-127, 1974[Free Full Text].

33.   Templeton, G. H., K. Wildenthal, J. T. Willerson, and J. H. Mitchell. Influence of acute myocardial depression on left ventricular stiffness and its elastic and viscous components. J. Clin. Invest. 56: 278-285, 1975.

34.   Thorson, J., and D. C. S. White. Role of cross-bridge distortion in the small-signal mechanical cynamics of insect and rabbit striated muscle. J. Physiol. (Lond.) 343: 59-84, 1983[Abstract/Free Full Text].

35.   Vaartjes, S. R., and H. B. K. Boom. Left ventricular internal resistance and unloaded ejection flow assessed from pressure-flow relations: a flow-clamp study on isolated rabbit hearts. Circ. Res. 60: 727-737, 1987[Abstract/Free Full Text].

36.   Vukas, M., I. Malek, and A. Hjalmarson. Myocardial depressant effect of vibration in isolated rabbit heart. Scand. J. Clin. Lab. Invest. 38: 421-424, 1978[Medline].

37.   Vukas, M., R. Silvertsson, and B. Ljung. Inhibitory effects of vibration on contractiltiy of isolated rabbit papillary muscle. Scand. J. Clin. Lab. Invest. 38: 415-419, 1978[Medline].

38.   Wijkstra, H., and H. B. Boom. Deactivation in the rabbit left ventricle induced by constant ejection flow. IEEE Trans. Biomed. Eng. 36: 1113-1123, 1989[Medline].

39.   Zhao, Y., and M. Kawai. Inotropic agent EMD-53998 weakens nucleotide and phosphate binding to cross bridges in porcine myocardium. Am. J. Physiol. 271 (Heart Circ. Physiol. 40): H1394-H1406, 1996[Abstract/Free Full Text].


AJP Heart Circ Physiol 273(4):H2044-H2061
0363-6135/97 $5.00 Copyright © 1997 the American Physiological Society



This article has been cited by other articles:


Home page
Am. J. Physiol. Heart Circ. Physiol.Home page
K. B. Campbell, Y. Wu, A. M. Simpson, R. D. Kirkpatrick, S. G. Shroff, H. L. Granzier, and B. K. Slinker
Dynamic myocardial contractile parameters from left ventricular pressure-volume measurements
Am J Physiol Heart Circ Physiol, July 1, 2005; 289(1): H114 - H130.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Regul. Integr. Comp. Physiol.Home page
T. Wronski, P. B. Persson, E. Seeliger, A. Harnath, and B. Flemming
Coupling of left ventricular and aortic volume elasticity in the rabbit
Am J Physiol Regulatory Integrative Comp Physiol, August 1, 2000; 279(2): R539 - R547.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Heart Circ. Physiol.Home page
K. B. Campbell, Y. Wu, R. D. Kirkpatrick, and B. K. Slinker
Myocardial contractile depression from high-frequency vibration is not due to increased cross-bridge breakage
Am J Physiol Heart Circ Physiol, April 1, 1998; 274(4): H1141 - H1151.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Campbell, K. B.
Right arrow Articles by Slinker, B. K.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Campbell, K. B.
Right arrow Articles by Slinker, B. K.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Visit Other APS Journals Online