Quantification of Global Diastolic Function by Kinematic Modeling-based Analysis of Transmitral Flow via the Parametrized Diastolic Filling Formalism

Podsumowanie

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.

Streszczenie

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 (x_o). 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/2kx_o², maximum A-V pressure gradient kx_o, load independent index of diastolic function, etc.) and select the aspect of physiology or pathophysiology to be quantified.

Wprowadzenie

Pioneering studies by Katz¹ 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 validated¹⁷. 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 2^nd order differential equation has three parameters: k- chamber stiffness, c- viscoelasticity/relaxation, and x_o- 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 studies^17-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 parameters²¹. 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.

Access restricted. Please log in or start a trial to view this content.

Protokół

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

1. Record a complete 2D/echo-Doppler study for all subjects in accordance with American Society of Echocardiography criteria¹⁶. NOTE: The screening echocardiograms were recorded on a standard clinical imager by a sonographer. If desired, additional transthoracic echocardiographic recording can be performed for verification purposes after a suitable, high fidelity catheter is advanced into the LV to measure LV hemodynamics simultaneously.
2. Image subjects in the supine position. In a nonresearch setting, standard left lateral positioning can be used without loss of generality of the method. Obtain apical four-chamber views using a 2.5 MHz transducer, with the sample volume gated at 1.5-5 mm directed between the tips of the mitral valve leaflets and orthogonal to the MV plane (to minimize alignment effects as seen on color M-mode Doppler), the wall filter set at 1 (125 Hz) or 2 (250 Hz), the baseline adjusted to take advantage of the full height of the display and the velocity scale adjusted to exploit the dynamic range of the output without aliasing.
3. Perform Doppler tissue imaging with the sample volume gated at 2.5 mm and positioned at the lateral and septal potions of the mitral annulus.
4. Save Doppler examinations in DICOM format in the echo machine and record on DVD with simultaneously recorded electrocardiogram (ECG).

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

1. Convert images from the DICOM format and video to bitmap (.bmp) files (using a custom MATLAB program). NOTE: The procedure described below to fit Doppler E-waves and tissue Doppler E’-waves is shown in Figure 1.
2. Load the bitmap image files on another custom MATLAB program to measure conventional transmitral flow parameters such as E_peak, A_peak, E_dur, E’_peak, A’_peak, etc. and crop the images for PDF analysis. Select images with discernible transmitral flow contour and complete cardiac cycle as indicated by ECG for analysis.
3. Mark the time sampling rate (measured in pixels/s on the horizontal axis) and velocity sampling rate (measured in pixels/(m/sec) along the vertical axis) in the images. Identify the complete cardiac cycle by noting and marking consecutive R peaks (or any distinct feature of the ECG) on the image.
4. Mark the transmitral Doppler E- and A-wave or tissue Doppler E’- and A’- wave in the selected cardiac cycle.
   1. Select the Doppler E-wave peak point i.e. E_peak, (or E’_peak) and mark the start of the wave using the line connecting the peak to the start as a guide to match the acceleration slope of the E-wave (or E’-wave). The start of the wave is used to calculate the interval from start to peak flow denoted as the E-wave (or E’-wave) acceleration time (AT).
   2. Mark the end of the E-wave (or E’-wave) using the line connecting the peak to the end as a guide to match the deceleration slope. This is used to calculate the interval from the peak to the baseline denoted as the deceleration time (DT). The interval from start to end of the wave is the duration of the E-wave (E_dur = AT+DT). The program guides the user through the entire process with appropriate instructions.
5. Mark the A-wave using a similar procedure as the E-wave. With both the E- and A-waves marked the program calculates the E_peak/A_peak ratio.
   NOTE: The program saves the marked waves as cropped images containing the E- and A-waves only. The program also creates a data file with the cropping and measured parameters for each beat.

4. Automated Fitting of Transmitral Flow Using the PDF Formalism

1. The automated fitting of Doppler E- and A-wave and tissue Doppler E’- and A’- wave contours is done using a custom LabView program^18,19.
   1. Load the cropped image, and the program automatically calculates the maximum velocity envelope (MVE). Select the MVE by setting the threshold such that MVE approximates transmitral flow as shown in Figure 1. The onset and termination of the points that define the MVE can be selected along the time-axis by the operator such that only MVE points that provide good correspondence to the actual selected portion of the wave are used as input for the subsequent fitting.
2. NOTE: The user-selected MVE points are the input to the computer program that automatically fits the PDF model solution for velocity as a function of time using a Levenberg- Marquardt (iterative) algorithm. The fitting is accomplished with the requirement that the mean square error between the clinical (input) data (MVE) and the PDF model predicted contour be minimized. Since the model is linear, a unique set of parameters is obtained for each Doppler E-wave derived MVE used as input. Thus numerically unique k, c, and x_o values are generated for each E-wave and k’, c’, and x_o’ for each E’-wave.
3. In the event the fit is obviously suboptimal when the fit is superimposed on the E-wave (or E’-wave) image (i.e. the algorithm attempted to fit noise included in the MVE for example) modify the MVE by using more/less points, thereby modifying the model predicted contour with consequent modification of PDF parameters to achieve a better fit.

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.

Access restricted. Please log in or start a trial to view this content.

Wyniki

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 c...

Access restricted. Please log in or start a trial to view this content.

Dyskusje

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 studies¹⁶. 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 bein...

Access restricted. Please log in or start a trial to view this content.

Materiały

Name	Company	Catalog Number	Comments
Philips iE33	Philips (Andover, MA)
LabView 6.0	National Instruments		Version 6.0.2
MATLAB	MathWorks		Version R2010b