Accurate, causality-based quantification of global diastolic function has been achieved by kinematic modeling-based analysis of transmitral flow via the Parametrized Diastolic Filling (PDF) formalism. PDF generates unique stiffness, relaxation, and load parameters and elucidates 'new' physiology while providing sensitive and specific indexes of dysfunction.
Quantitative cardiac function assessment remains a challenge for physiologists and clinicians. Although historically invasive methods have comprised the only means available, the development of noninvasive imaging modalities (echocardiography, MRI, CT) having high temporal and spatial resolution provide a new window for quantitative diastolic function assessment. Echocardiography is the agreed upon standard for diastolic function assessment, but indexes in current clinical use merely utilize selected features of chamber dimension (M-mode) or blood/tissue motion (Doppler) waveforms without incorporating the physiologic causal determinants of the motion itself. The recognition that all left ventricles (LV) initiate filling by serving as mechanical suction pumps allows global diastolic function to be assessed based on laws of motion that apply to all chambers. What differentiates one heart from another are the parameters of the equation of motion that governs filling. Accordingly, development of the Parametrized Diastolic Filling (PDF) formalism has shown that the entire range of clinically observed early transmitral flow (Doppler E-wave) patterns are extremely well fit by the laws of damped oscillatory motion. This permits analysis of individual E-waves in accordance with a causal mechanism (recoil-initiated suction) that yields three (numerically) unique lumped parameters whose physiologic analogues are chamber stiffness (k), viscoelasticity/relaxation (c), and load (xo). The recording of transmitral flow (Doppler E-waves) is standard practice in clinical cardiology and, therefore, the echocardiographic recording method is only briefly reviewed. Our focus is on determination of the PDF parameters from routinely recorded E-wave data. As the highlighted results indicate, once the PDF parameters have been obtained from a suitable number of load varying E-waves, the investigator is free to use the parameters or construct indexes from the parameters (such as stored energy 1/2kxo2, maximum A-V pressure gradient kxo, load independent index of diastolic function, etc.) and select the aspect of physiology or pathophysiology to be quantified.
Pioneering studies by Katz1 in 1930 revealed that the mammalian left ventricle initiates filling by being a mechanical suction pump, and much effort since then has been devoted to unraveling the workings of diastole. For many years, invasive methods were the only options available for clinical or research assessment of diastolic function (DF)2-16. In the 1970s, however, technical advancements and developments in echocardiography finally gave cardiologists and physiologists practical tools for noninvasive characterization of DF.
Without a unifying causal theory or paradigm for diastole regarding how the heart works when it fills, researchers proposed numerous phenomenologic indexes based on correlation with clinical features. The curvilinear, rapidly rising and falling shape of the transmitral blood flow velocity contour during early, rapid filling, for example, was approximated as a triangle and diastolic function indexes were defined from geometric features (height, width, area, etc.) of that triangle. Technical advancements in echocardiography have allowed tissue motion, strain, and strain rate during filling to be measured, for example, and each technical advancement brought with it a new crop of phenomenological indexes to be correlated with clinical features. However, the indexes remain correlative and not causal and many indexes are different measures of the same underlying physiology. It’s not surprising, therefore, that currently employed clinical indexes of DF have limited specificity and sensitivity.
To overcome these limitations the Parametrized Diastolic Filling (PDF) formalism, a causal kinematic, lumped parameter model of left ventricular filling that is motivated by and incorporates the suction-pump physiology of diastole was developed and validated17. It models diastolic function (as manifested by the curvilinear shapes of transmitral flow contours) in accordance with the rules of damped harmonic oscillatory motion. The equation for damped harmonic oscillatory motion is based on Newton’s Second Law and can be written, per unit mass, as:
Equation 1
This linear 2nd order differential equation has three parameters: k– chamber stiffness, c– viscoelasticity/relaxation, and xo– the oscillator’s initial displacement/preload. The model predicts that the different clinically observed diastolic filling patterns are the result of variation in the numerical value of these three model parameters. Based on the PDF formalism and classical mechanics, E-waves can be classified as being determined by under-damped or over-damped regimes of motion. Numerous studies17-21 have validated that clinically recorded E-wave contours and PDF model predicted contours show superb agreement and have elucidated the hemodynamic/physiologic analogues of the three PDF parameters21. The process for extracting model parameters from clinically recorded E-wave data is detailed in the methods below.
Unlike typical indexes of DF in current clinical use, the PDF model’s three parameters are causality based. As discussed in the methods below, additional indexes of diastolic physiology can be derived from these fundamental parameters and from application of the PDF formalism to aspects of diastole other than transmitral flow. In this work, methods of PDF-based analysis of transmitral flow and the physiologic relations that can be drawn from the PDF approach, its parameters and the derived indexes are described. Additionally, it is shown that the PDF parameters or indexes derived from them can tease apart intrinsic chamber properties from the external effects of load can provide correlates to traditional invasively defined parameters and can differentiate between normal and pathologic groups.
The procedure for acquiring echocardiographic images and analyzing them to obtain the PDF parameters is detailed below. Although cardiac catheterization is mentioned in the subject selection portion below, the methodology described applies only to the echocardiographic portion. The description of the catheterization portion was included for independent validation of model based predictions and is unrelated to the analysis of E-waves via the PDF formalism. Prior to data acquisition, all subjects provide signed, informed consent for participation in the study in accordance with the Institutional Review Board (Human Research Protection Office) at Washington University School of Medicine.
NOTE: All software programs (along with tutorials on how to use them) described in this section can be downloaded from http://cbl1.wustl.edu/SoftwareAgreement.htm
1. Subject Selection
NOTE: All subjects in the Cardiovascular Biophysics Laboratory Database had simultaneous echocardiography and cardiac catheterization performed and were referred by their physicians for diagnostic cardiac catheterization. The database inclusion criteria are: 1) absence of any significant valvular abnormalities, 2) absence of wall motion abnormalities or bundle branch block on ECG, 3) presence of a satisfactory echocardiographic window with clearly identifiable E- and A-waves.
2. Echocardiographic Data Acquisition
3. Doppler Image Processing and Conventional Analysis
NOTE: This section describes two custom MATLAB programs. The first program is described in step 3.1 and the second program is described in steps 3.2-3.5. All software programs (along with tutorials on how to use them) can be downloaded from http://cbl1.wustl.edu/SoftwareAgreement.htm
4. Automated Fitting of Transmitral Flow Using the PDF Formalism
Save the data when the appropriate PDF fit has been generated. NOTE: The program is written to automatically save the data in image and text files containing the PDF parameters and the contour information.
The PDF parameters obtained from the procedure described above can be used to elucidate new physiology and distinguish between normal and pathological physiology as detailed in the Representative Results section below.
Doppler waveforms representative of the four different types of filling patterns (normal, pseudonormal, delayed relaxation, constrictive-restrictive) using the method detailed above are shown in Figure 2. Figure 2A shows the normal pattern, which, by itself is indistinguishable from the pseudonormal pattern. Figure 2B shows a delayed relaxation and Figure 2C shows a constrictive-restrictive pattern associated with severe diastolic dysfunction. For clarity, the PDF model-predicted fits are overlaid on the images. The conventional echo parameters (Epeak, Apeak, E-wave AT, and E-wave DT) and the PDF parameters (k, c, xo) are listed below each image. As the figures indicate, the PDF formalism fits (predicts) all three of these filling patterns very well. The PDF parameters also provide information on chamber properties. The delayed relaxation pattern (Figure 2B) typically has higher viscoelasticity/relaxation PDF parameter c than the normal pattern (Figure 2A). Constrictive-restrictive pattern (Figure 2C) typically has a higher stiffness (PDF parameter k) than the normal pattern.
Analysis of Doppler E-waves using the PDF formalism has been used to differentiate between normal and pathological groups and to discover new physiology. Listed below are some selected published results of PDF formalism based DF analysis intended to differentiate between pathologic and normal physiology and selected applications of the PDF formalism to elucidate new physiology.
DIABETES
The method has been shown to quantify the differences in DF between diabetic and age-matched control subjects. While the conventional indexes such as E-wave deceleration time- DT, Epeak or the time constant of isovolumic relaxation were unable to differentiate between the filling patterns of the groups, the PDF viscoelasticity/relaxation parameter c was significantly different between groups22. Additionally, the peak atrio-ventricular pressure gradient, which can be calculated from the PDF parameters as kxo23 was significantly higher in the diabetic group. Also, see kinematic filling efficiency, applied to diabetics below.
HYPERTENSION
The method has been used to analyze transmitral filling patterns in hypertensive subjects compared to controls24. Conventional Doppler derived indexes were unable to differentiate between groups but the PDF parameter c was significantly higher in the hypertensive subjects group compared to nonhypertensive controls.
CALORIC RESTRICTION SLOWS CARDIAC AGING
The method assessed the effect of caloric restriction on DF in humans25. DF was assessed in subjects practicing caloric restriction by measuring transmitral flow and comparing to age- matched controls. DF was significantly better in the caloric restriction group as quantified by higher value of E/A and higher early filling (E-wave) fraction. Additionally, the PDF parameter k, representing LV chamber stiffness, and c, representing viscoelasticity, was significantly lower in caloric restriction subjects. Since the Epeak was not significantly different between the two groups, the control group expends more energy to achieve the same peak filling velocity. This revealed that caloric restriction is associated with more efficient DF. Moreover the filling in elderly caloric restricted subjects was comparable to a younger normal cohort, suggesting that caloric restriction slows cardiac aging26.
PRESENCE VS. ABSENCE OF MITRAL ANNULAR OSCILLATIONS
The PDF formalism has also been used to analyze mitral annular oscillations (MAO) after the E’-wave (the E”-wave, E’’’-wave, etc.). This ‘ringing’ of the mitral annulus has been observed in humans20 but characterization of the presence and absence of the subsequent oscillations was lacking. The method allowed the hypothesis to be tested that the absence of MAO is explained by increased viscoelastic effects due to less or slower effective relaxation. By comparing 35 subjects with MAO to 20 subjects without MAO, it was found that the longitudinal stiffness (k’) and the longitudinal viscoelasticity/relaxation (c’) were higher in the group without MAO. The initial recoil force and the stored recoil energy both were higher in the group with MAO. In addition, it was shown that the absence of MAO was concordant with relaxation-related diastolic dysfunction27. Hence the PDF analysis of tissue Doppler E’- waves reveals that the absence of MAO indicates relaxation related diastolic dysfunction.
DIASTATIC STIFFNESS FROM E-WAVE ANALYSIS
While the slope of the end-diastolic pressure-volume relationship (EDPVR) provides the familiar stiffness-based index, the slope (ΔP/ΔV) of the diastatic pressure-volume (P-V) relationship (D-PVR) provides the in-vivo stiffness of the relaxed LV. Echocardiographic, (i.e. Doppler E-wave), analysis can provide only relative, rather than absolute pressure information. Accordingly, it has been shown that the relaxed (diastatic) stiffness of the LV can be computed directly from E-wave analysis alone28. Using the PDF formalism and Bernoulli’s equation pressure and volume at diastasis (end of E-wave) is derived. The derived P, V points when fit via linear regression generate the D-PVR from E-wave analysis (D-PVRE-wave) whose slope, diastatic stiffness KE-wave was computed. The results yielded excellent correlation (R2=0.92) between diastatic stiffness from PDF based E-wave analysis (KE-wave) and the simultaneous gold standard measurement of diastatic stiffness from simultaneous P-V data (KCATH) in 30 subjects (444 total cardiac cycles) with normal LVEF (LVEF>55%).
KINEMATIC FILLING EFFICIENCY INDEX
From a kinematic modeling perspective, an increased relaxation/viscosity constant c generates increased resistance to filling. Hence a natural choice for idealized ventricular filling is a scenario due to recoil only and complete relaxation, i.e. no damping (c = 0). The kinematic filling efficiency index (KFEI) was defined and derived29 as the dimensionless ratio of actual volume entering the left ventricle (LV) (velocity time integral [VTI] of real E-wave with PDF parameters c, k, xo) to the ideal volume (VTI of ideal E-wave having same k and xo but with no resistance filling [c = 0]). In 36 patients with normal ventricular function (17 diabetic and 19 well-matched nondiabetic controls) it was shown that30 KFEI of E-waves in diabetic patients (49.1±3.3%) was significantly lower than in normal patients (55.8±3.3%). This means that even when LVEF is normal, filling efficiency is impaired in diabetics compared to nondiabetics.
FILLING EFFICIENCY DETERIORATES WITH AGE
In light of the ability of kinematic filling efficiency index (KFEI)29 to assess filling in diabetes vs. nondiabetic controls, the age dependence of KFEI was determined. It was shown that KFEI, decreases in magnitude with age and correlates very strongly with age (R2=0.80) by analyzing 72 control subjects with normal LVEF (LVEF>55%) and without cardiovascular pathology30. The age dependence of other conventional parameters of DF was also evaluated. In concordance with other noninvasive DF measures known to decrease with age, KFEI decreases and correlates very strongly with age (R2=0.80). Multivariate analysis showed that age is the single most important contributor to KFEI (p=0.003).
LOAD INDEPENDENT INDEX OF DIASTOLIC FUNCTION
E-wave contours demonstrate beat-by-beat changes in response to respiration and hence demonstrate strong load dependence. Indeed all indexes of DF are load-dependent. This is problematic because it calls into question whether observed differences in DF indexes are the result of load variation or the result of intrinsic chamber property variation. Theoretical prediction and experimental validation of a load independent index of diastolic function (LIIDF) has been a long sought unsolved problem in physiology/cardiology. To address the question of load dependence, the PDF formalism was applied to E-waves measured at variable loads. Through kinematic modeling and mathematical derivation, a load independent index was derived, which is conserved between E-waves measured at different loads. For each measured E-wave, the PDF parameters k and xo are multiplied to yield kxo, the model predicted peak force value analogous to the peak instantaneous pressure gradient driving flow, and the PDF parameter c is multiplied by the peak velocity Epeak to yield a value for the peak force resisting filling. Plotting kxo vs. cEpeak as an ordered pair for each E-wave generates a highly linear relationship whose (dimensionless) slope M is the sought after load independent index and remains conserved despite load generated changes in E-waves.
For validation E-waves recorded while load was varied via tilt table (head up, horizontal, and head down) in 16 healthy volunteers was analyzed. The results33 yielded the very high correlation (R2 = 0.98) between kxo and cEpeak as predicted. The ability of M to differentiate between normal and diastolic dysfunction subjects was also assessed by analysis of simultaneous cath-echo data in diastolic dysfunction subjects vs. controls. Average M for the diastolic dysfunction group (M = 0.98±0.07) was significantly lower than controls (M = 1.17±0.05, P<0.001)33.
Figure 1. Sequence of operational steps for fitting (A) an E-wave and (B) an E’-wave via the PDF formalism. A) From Left to Right- Transmitral flow image is cropped to obtain Doppler velocity profile. E-wave maximum velocity envelope (MVE) to be fit is selected (shown in green with time limits in blue). Error minimizing PDF fit is obtained via Levenberg- Marquardt algorithm resulting in PDF parameters and a measure of goodness of fit. B) Similar procedure for the tissue Doppler image. Image is inverted after cropping. See text for details. Please click here to view a larger version of this figure.
Figure 2. Three E-wave patterns with PDF fits. A) Normal/Pseudonormal filling pattern. B) Delayed relaxation pattern. C) Constrictive-restrictive pattern. See text for details. Please click here to view a larger version of this figure.
In keeping with our methodologic focus, the key aspects of the methods that facilitate obtaining accurate and meaningful results are highlighted.
ECHOCARDIOGRAPHY
The American Society of Echocardiography (ASE) has guidelines for the performance of transthoracic studies16. During an echo exam, there are a multitude of factors that affect image quality. Factors that are beyond the control of the sonographer include: technical capabilities of the imager being used, heart rate, patient body habitus, individual variation in the location, orientation of anatomical structures, and quality of ‘echo window’, referring to the characteristics of ultrasound transmission in a given subject’s tissue. Factors that are directly controllable by the sonographer include machine settings, including choice of transducer. Since the fidelity of the PDF analysis is dependent on the echo image quality, care should be taken during the image acquisition process to obtain the best possible images.
For optimal E-wave image quality for PDF analysis, maximizing E-wave size relative to the display and setting the sweep speed to 100 mm/sec are desirable. High sweep speed and use of the full display size in determining maximum velocity scale provides increased temporal resolution (i.e. more points to be fit) along both time and velocity axes. Baseline filter settings can also be better determined with higher sweep speed settings. The number of cardiac cycles recorded is highly variable between echo labs. For meaningful PDF analysis continuous recording through several (3 or 4) respiratory cycles is most desirable. At a typical resting heart rate of 75 beats/min, and 12 respirations/min, 4 respiratory cycles amount to 20 sec of continuous recording that should provide 25 cardiac cycles. Recording this number of cycles is justified because of the load varying consequence of quiet respiration, so that the LIIDF can be computed if desired. Note, that computing values for xo, c, and k based on the 25 beat average is a legitimate way to characterize diastole. Load variation can also be generated during clinical recording by the Valsalva or Mueller maneuvers, or by passive leg elevation using a 30° foam wedge.
PDF PARAMETER DETERMINATION
ALGORITHMIC DETAILS
The equation of motion for a damped harmonic oscillator and its mathematical solution is standard course content in engineering mathematics, in physics and mechanics34. The choice of computer language (C++, Fortran, LabView, MATLAB, etc.) by which it is implemented is also at the discretion of the user/investigator. Standard numerical methods exist and are well known35. Other groups have implemented the PDF formalism by writing their own numerical algorithm and have independently replicated our results, including numerical values for the PDF parameters36 in a large study involving well over 1,000 patients. While ongoing work includes developing web-based PDF analysis tools, the optimal, broad reaching benefit of the method could be best achieved by incorporation of the PDF formalism into the proprietary analysis package of commercial echocardiographic imagers.
OPERATOR DEPENDENT ASPECTS
Once the E-wave image has been imported and cropped (see Figure 1) determination of the maximum velocity envelope, i.e. the actual set of points to which the solution of damped harmonic oscillatory velocity is to be fit by the method, is determined. As shown by the sequence of panels and operational steps in Figure 1 and discussed above, baseline noise as well as extraneous noise that affects the contour is often part of the image. The operator can determine the continuous set of points to be fit, as shown in Figure 1, by adjusting the position of the vertical blue lines which define the start and end of the points to be fit. The method displays the fit directly over the imported image and the operator can easily assess if it is meaningful or not.
Heart rate has an effect on the duration of diastole and the features of the E-wave37, and care must be taken to interpret the results of the fitting algorithm in the context of patient heart rate. At typical heart rates below 80 beats/min, in sinus rhythm E-and A-waves are separated by a brief period of diastasis. This facilitates inclusion of the deceleration portion of the E-wave. As heart rates increase, diastasis diminishes and disappears, since A-wave onset occurs before E-wave termination. At fast heart rates, above 90 beats/min, the A-wave overlies the deceleration portion of the E-wave and PDF analysis of the E-wave becomes unreliable because of the limited number of MVE points available to be fit. For meaningful analysis at least 1/2 to 2/3 of the total deceleration E-wave waveform should be available for fitting.
The authors have nothing to disclose.
This work was supported in part by the Alan A. and Edith L. Wolff Charitable Trust, St. Louis, and the Barnes-Jewish Hospital Foundation. L. Shmuylovich and E. Ghosh were partially supported by predoctoral fellowship awards from the Heartland Affiliate of the American Heart Association. S. Zhu received partial support from the Washington University Compton Scholars Program and the College of Arts and Sciences’ Summer Undergraduate Research Award. S. Mossahebi received partial support from the Department of Physics.
Name of Equipment/ Software | Company | Catalog Number | Comments/Description |
Philips iE33 | Philips (Andover, MA.) | ||
LabVIEW 6.0 | National Instruments | Version 6.0.2 | |
MATLAB | MathWorks | Version R2010b |