s ndor j kov cs phd md washington university st louis
Download
Skip this Video
Download Presentation
Sándor J Kovács PhD MD Washington University, St. Louis

Loading in 2 Seconds...

play fullscreen
1 / 76

Sándor J Kovács PhD MD Washington University, St. Louis - PowerPoint PPT Presentation


  • 64 Views
  • Uploaded on

Discovering (predicting) new cardiac physiology/function from cardiac imaging, mathematical modeling and first principles. Sándor J Kovács PhD MD Washington University, St. Louis. UCLA/IPAM 2/6/06. Imaging and modeling allows us to go beyond. correlation to…. causality!.

loader
I am the owner, or an agent authorized to act on behalf of the owner, of the copyrighted work described.
capcha
Download Presentation

PowerPoint Slideshow about 'Sándor J Kovács PhD MD Washington University, St. Louis' - phong


An Image/Link below is provided (as is) to download presentation

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.


- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -
Presentation Transcript
s ndor j kov cs phd md washington university st louis
Discovering (predicting) new cardiac physiology/function from cardiac imaging, mathematical modeling and first principles

Sándor J Kovács PhD MD

Washington University, St. Louis

UCLA/IPAM 2/6/06

slide2
Imaging and modeling

allows us to go beyond

correlation to…

causality!

UCLA/IPAM 2/6/06

slide3
Focus: How the Heart Works When it Fills

The physiologic process by which the heart fills has confused cardiologists, physiologists, biomedical engineers, medical students and graduate students for generations.

UCLA/IPAM 2/6/06

slide4
How the Heart Works When it Fills

Why does it matter?

The recent recognition that up to 50% of patients

admitted to hospitals with congestive heart failure

have ‘normal systolic function’ as reflected by

ejection fraction, has further emphasized the

need to more fully understand the physiology

of diastole.

UCLA/IPAM 2/6/06

slide5
How the Heart Works When it Fills

In an effort to quantitate diastolic function using

a number or an index, the filling process has been

characterized via correlations of selected features of

either fluid (blood) flow or tissue displacement or

motion to LV ejection fraction, end-diastolic pressure

and other observables or clinical correlates such as

exercise tolerance or mortality.

UCLA/IPAM 2/6/06

slide6
How the Heart Works When it Fills

What do we know?

Anatomy

UCLA/IPAM 2/6/06

how the heart works anatomy
How the Heart Works:anatomy

Pericardial anatomy

UCLA/IPAM 2/6/06

how the heart works anatomy1
How the Heart Works: anatomy

Pericardial anatomy

UCLA/IPAM 2/6/06

slide9
How the Heart Works When it Fills

Anatomy and terminology

UCLA/IPAM 2/6/06

how the heart works anatomy2
How the Heart Works: anatomy

Pericardial anatomy

UCLA/IPAM 2/6/06

slide11
How the Heart Works When it Fills

What else do we know?

Physiology

UCLA/IPAM 2/6/06

slide12
Doppler echocardiography reveals physiology:

Method by which transmitral

Doppler flow velocity data

is acquired

UCLA/IPAM 2/6/06

slide13
Echocardiographically observed patterns of filling:

S2 = second heart sound,

IR = isovolumic relaxation,

AT = acceleration time,

DT= deceleration time.

(Note: velocity scales differ slightly among images)

Waveform features (Epeak, E/A, DT, …) are correlated

with clinical aspects.

UCLA/IPAM 2/6/06

slide14
Cardiac catheterization reveals physiology:

Simultaneous, high fidelity LAP, LVP and transmitral Doppler

in closed chest canine. Note reversal of sign of A-V pressure gradient

As flow accelerates

(LAP > LVP) and decelerates (LAP < LVP).

Isovolumic

Relaxation

Rapid

Filling

Atrial

Systole

Diastasis

UCLA/IPAM 2/6/06

slide15
Cardiac catheterization reveals physiology:

Simultaneous aortic root,

LV pressure and LV volume as a function of time for one cardiac cycle

as measured in the cardiac catheterization laboratory.

dP/dV<0 at MVO

UCLA/IPAM 2/6/06

slide16
Cardiac catheterization reveals physiology:

AO

AVC

IVR

AVO

MVO

MVC

LA

LV

diastasis

atrial systole

Doppler A-wave

rapid filling

Doppler E-wave

UCLA/IPAM 2/6/06

slide17
Mechanics of filling:

Ventricle fills in 2 phases:

1) Early, rapid-filling (dP/dV< 0)

2) Atrial filling (dP/dV > 0)

(Actually, diastole has 4 phases: isovolumic relaxation, early rapid filling, diastasis, atrial contraction)

UCLA/IPAM 2/6/06

slide18
Catheterization and echo -combined

Abnormal

Pseudo-

Restriction

Restriction

relaxation

normalization

(reversible)

(irreversible)

Normal

40

0

N-

Mean LAP

TAU

NYHA

I-II

II-III

III-IV

IV

Grade

I

II

III

IV

UCLA/IPAM 2/6/06

slide19
How the Heart Works When it Fills

Recall key physiologic fact:

At -(and for a while after) - MVO, the LV simultaneously decreases its pressure while increasing its volume!

UCLA/IPAM 2/6/06

slide20
How the Heart Works When it Fills

We must therefore conclude that:

The heart is a suction pump in early diastole!

UCLA/IPAM 2/6/06

slide21
To go from correlation to causality devise a

kinematic model of suction initiated filling:

Newton’s Law: md2x/dt2 + c dx/dt + kx = 0

Initial conditions: x(0) = xo  stored elastic strain to power suction

v(0) = 0  no flow prior to valve opening

Recall SHO has 3 regimes of motion,underdampedc2-4mk<0, critically dampedc2=4mk, overdampedc2 - 4mk>0.

VALIDATION: Compare model-predicted velocity of oscillator

to velocity of blood entering the ventricle through mitral valve.

UCLA/IPAM 2/6/06

slide22
Model of suction initiated filling:

Block-diagram of operational steps

Result: 1) re-express all E-and A-waves in terms of parameters AND

2) compute physiologic indexes

UCLA/IPAM 2/6/06

slide23
Model of suction initiated filling: does it fit the data?

Examples of model’s ability to fit in-vivo Doppler data

UCLA/IPAM 2/6/06

slide24
Model prediction compared to actual data:

Observed patterns of mitral valve inflow and superimposed model fits

S2 = second heart sound,

IR = isovolumic relaxation,

AT = acceleration time,

DT= deceleration time.

(Note: velocity scales differ slightly among images)

UCLA/IPAM 2/6/06

slide25
Kinematic model of suction initiated filling compared

to non-linear, coupled PDE models of filling:

Comparison of the PDF (red),

Meisner (blue) and Thomas (green) models for a clinical Doppler image. Note that all three models reproduce the contour of the image with comparable accuracy, and that the three models’ predictions are

essentially indistinguishable

graphically from one another.

UCLA/IPAM 2/6/06

slide26
Kinematic model of suction initiated filling:

Indexes from model parameters:

Mechanical Physiologic

kxo Force in spring Maximum A-V pressure

k Spring constant Chamber stiffness

1/2kxo2Stored energy Stored elastic strain

xoSpring displacement Velocity-time integral of E-wave

c2-4mk Regime of motion Stiff vs. delayed relaxation

UCLA/IPAM 2/6/06

slide27
Kinematic model of suction initiated filling:

Predictions from kinematic modeling:

1) The spring is linear and it is bi-directional

2) Underdamped, critically damped, overdamped regimes

3) Existence of ‘load independentindex’ of filling

4) Equilibrium volume of LV is diastasis

5) Tissue oscillations

6) Resonance

UCLA/IPAM 2/6/06

slide28
Kinematic model of suction initiated filling:

Physiologic analog and prediction of model:

Q: What is the spring?

UCLA/IPAM 2/6/06

what is the spring
What is the ‘spring’?

How the experiment that shows that cells can push was done!

Titin Develops Restoring Force in

Rat Cardiac Myocytes

Michiel Helmes, Károly Trombitás, Henk Granzier

Circulation Research. 1996;79:619-626.

UCLA/IPAM 2/6/06

what is the spring1
What is the ‘spring’?

Experimental data proving that titin acts as a linear, bi-directional spring

It is hinged between thick and thin filaments.

UCLA/IPAM 2/6/06

slide31
Model of suction initiated filling:

Model can be used to

fit and (?) explain

heretofore unexplained

mechanism of biphasic

E-waves.

Early portion is governed by

k dominance, (underdamped)

later portion is governed by

c dominance (overdamped).

UCLA/IPAM 2/6/06

slide32
Kinematic modeling of filling:

“When you solve one difficulty, other new difficulties arise. You then try to solve them. You can never solve all difficulties at once.” P.A.M. Dirac

UCLA/IPAM 2/6/06

slide33
Modeling how the heart works:

Recall a physiologic fact -

Although the heart is an oscillator:

It is possible to remain (essentially) motionless!

UCLA/IPAM 2/6/06

slide34
Modeling how the heart works:

Hence:

The four-chambered heart is a

constant- volume pump!

UCLA/IPAM 2/6/06

how the heart works constant volume
How the Heart Works :(constant volume)
  • Constant-volume attribute of the four-chambered heart -
    • Hamilton and Rompf -1932

Hamilton W, Rompf H. Movements of the Base of the Ventricle and the Relative Constancy of the Cardiac Volume. Am J Physiol. 1932;102:559-65.

    • Hoffman and Ritman -1985

Hoffman EA, Ritman E. Invariant Total Heart Volume in the Intact Thorax. Am J Physiol. 1985;18:H883-H890. Also showed that Left heart and Right heart are very nearly constant volume!

    • Bowman and Kovács - 2003

Bowman AW, Kovács SJ. Assessment and consequences of the constant-volume attribute of the four-chambered heart. American Journal of Physiology, Heart and Circulatory Physiology 285:H2027-H2033, 2003.

UCLA/IPAM 2/6/06

how the heart works when it fills constant volume
How the Heart Works When it Fills : (constant volume)

Cardiac MRI Cine Loop

‘four-chamber view”

Note relative absence

of ‘radial’ or ‘longitudinal’

pericardial surface

displacement or motion

UCLA/IPAM 2/6/06

how the heart works when it fills constant volume1
How the Heart Works When it Fills : (constant volume)

Cardiac MRI Cine Loop

‘LV outflow track view”

Note relative absence

of ‘radial’ or ‘longitudinal’

pericardial surface

displacement or motion

UCLA/IPAM 2/6/06

how the heart works when it fills constant volume2
How the Heart Works When it Fills : (constant volume)

Cardiac MRI Cine Loop

‘short-axis view”

Note slight ‘radial’ motion

of pericardial surface

UCLA/IPAM 2/6/06

how the heart works constant volume1
How the Heart Works:(constant volume)

Cardiac MRI Cine Loop

‘four-chamber view”

Normal, human

UCLA/IPAM 2/6/06

how the heart works constant volume2
How the Heart Works:(constant volume)

Cardiac MRI Cine Loop

‘short-axis view”

UCLA/IPAM 2/6/06

how the heart works constant volume3
How the Heart Works:(constant volume)

Plot of # of pixels vs. frame number for 4-chamber slice

Diastole

Systole

Area (in Pixels)

Frame #

UCLA/IPAM 2/6/06

how the heart works constant volume4
How the Heart Works:(constant volume)

Plot of # of voxels vs. fraction R-R interval for 3-D data set

Voxels

Fraction of R-R Interval

UCLA/IPAM 2/6/06

how the heart works constant volume5
How the Heart Works:(constant volume)

Constant-Volume Attribute of the Four-Chambered Heart

Via MRI - how are images analyzed? (with Bowman, Caruthers, Watkins)

Conclusion: In normal, healthy subjects, the total volume enclosed within the pericardial sack remains constant to within a few percent. The pericardial surface exhibits only slight radial displacement throughout the cardiac cycle most notably along its diaphragmatic aspect.

UCLA/IPAM 2/6/06

slide45
How the Heart Works:(constant volume)

Cine MRI loop of pericardium for one cardiac cycle

UCLA/IPAM 2/6/06

how the heart works constant volume6
How the Heart Works:(constant volume)

Right heart vs. left heart (n=20)

UCLA/IPAM 2/6/06

how the heart works constant volume7
How the Heart Works:(constant volume)

What are predictable consequences of a constant volume, four-chambered heart

as they pertain to diastole?

(In light of the previous slide showing that the volumes of left and right heart are also independently constant.)

UCLA/IPAM 2/6/06

how the heart works constant volume8
How the Heart Works:(constant volume)

Consider the motion of the mitral valve plane relative to the fixed apex and base.

Caltech 3/10/05

how the heart works constant volume9
How the Heart Works:(constant volume)

One dimensional analog of mitral valve plane motion

atriumventricle

UCLA/IPAM 2/6/06

slide50
How the Heart Works:(constant volume)

Normalized MVP displacement vs. cardiac cycle

Normalized displacement

Percentage of cardiac cycle

UCLA/IPAM 2/6/06

modeling how the heart works constant volume
atrioventricular

cross-section =A cm

2

Mitral valve area MVA cm

2

atrium

Mitral valve plane

-in diastole

myocardium

Mitral valve plane

velocity - Vmvp

Mitral valve plane

-in systole

ventricle

Modeling how the heart works:(constant volume)

Concept: Consider a

simplified 2-chamber

constant-volume geometry

Application:Derive

the mitral annular velocity

(E’) to Doppler E-wave

(filling velocity) relation

UCLA/IPAM 2/6/06

modeling how the heart works constant volume1
Modeling how the heart works:(constant volume)

Conservation of volume for the upper and lower portions of the cylinder also imply tissue volume is conserved. How does the idealized LV chamber appear as it fills?

UCLA/IPAM 2/6/06

modeling how the heart works constant volume2
Modeling how the heart works:(constant volume)

Conservation of volume for the upper and lower portions of the cylinder imply:

Amvp Vmvp = Amv VE

At every instant during early rapid filling (Doppler E-wave)!

Note: Amvp and Amv are constant!!

UCLA/IPAM 2/6/06

modeling how the heart works constant volume3
Modeling how the heart works:(constant volume)

Conservation of volume means:

Amvp Vmvp = Amv VE

At every instant during early rapid filling (Doppler E-wave)!

Note: time varying quantity = time varying quantity

Rewrite as:

Amvp /Amv =VE /Vmvp

constant = constant

UCLA/IPAM 2/6/06

modeling how the heart works validation
Modeling how the heart works: validation

Transmitral Doppler

Is constancy of

Amvp /Amv =VE /Vmvp really true?

UCLA/IPAM 2/6/06

modeling how the heart works validation1
Modeling how the heart works: validation

Mitral valve annular velocity via DTI

Is constancy of

Amvp /Amv =VE /Vmvp really true?

UCLA/IPAM 2/6/06

modeling how the heart works validation2
Modeling how the heart works: validation

Re-express the E and MVP velocity

Contours in terms of equivalent

contours using

PDF model and MBIP

(Note:time scale for lower picture

is expanded, velocity is plotted inverted)

UCLA/IPAM 2/6/06

modeling how the heart works validation3
Modeling how the heart works: validation

Overlay of VE and Vmvp on

the same velocity vs..

time coordinate axes.

Q1: What is peculiar about this?

Q2: What does it mean?

UCLA/IPAM 2/6/06

modeling how the heart works validation4
Modeling how the heart works: validation

mitral valve plane velocity (Vmvp) to Doppler E-wave (VE) relation - normal hearts

UCLA/IPAM 2/6/06

modeling how the heart works validation5
Modeling how the heart works: validation

Mitral valve plane velocity (Vmvp) to Doppler E-wave (VE) relation - data for

enlarged hearts!

UCLA/IPAM 2/6/06

modeling how the heart works validation prediction
Modeling how the heart works: validation + prediction

VE /VmvpRatio vs. LVEDP

Q1: it appears reasonably linear-

(r = 0.9196)WHY??

ANSWER: (- Hooke’s Law )

A =  LVEDP

VE/ VMVP = (/MVA) LVEDP

UCLA/IPAM 2/6/06

how the heart works constant volume10
How the Heart Works :(constant volume)

Relationship of [VE]max/ [Vmvp]max = E/E’ to

left ventricular end-diastolic pressure during simultaneous catheterization and echocardiography.

The ‘constant volume pump’ model predicted linear relationship is well fit by the data.

Best linear fit is provided by

E/E’ = 0.1753LVEDP + 1.8949 with r = 0.9196.

UCLA/IPAM 2/6/06

how the heart works constant volume11
How the Heart Works :(constant volume)

E/E’ tabulated for 24 normal

subjects and 3 subjects with clinical CHF.

The model predicted value of @ 4 for the E/E’ relationship for the normal group is well fit by the data showing E/E’ = 4.4 ± 1.15. Three subjects with known CHF have greater than normal E/E’ @ 10, in accordance with model prediction.

UCLA/IPAM 2/6/06

modeling how the heart works prediction of load independent index of filling
Modeling how the heart works: prediction of load independent index of filling

Initial conditions

  • m = inertia
  • c = damping
  • k = spring constant
  • xo = initial displacement of spring

UCLA/IPAM 2/6/06

predicted load independent index
[1]

[2]

kxo=M(cEpeak)+b

Peak AV gradient (kxo)

Peak viscous force (cEpeak)

PredictedLoad Independent Index

Changes in preload change the shape of the E-wave, and thus must cause changes in k, c, and xo

The equation of motion, however, is obeyed regardless of changes in preload:

Consider the equation of motion at time of the E wave peak, t = tpeak

While 2 is true of any SHO, we invoke physiology:

[3]

Which implies:

[4]

Thus the maximum initial driving force (kxo) to peak attained viscous force (cEpeak) relation is predicted to be linear and load independent.

UCLA/IPAM 2/6/06

modeling how the heart works validation6
Modeling how the heart works: validation
  • 15 healthy subjects (ages 20-30) with no history of heart disease and on no prescribed medication
  • Subjects were positioned at three predetermined angles on a tilt-table. Data was acquired after transient heart rate changes resolved.
  • E- and A-waves were recorded from subjects in supine, 90° head-up and 90° head-down positions.

UCLA/IPAM 2/6/06

modeling how the heart works validation7
Modeling how the heart works: validation

Load independent index of filling from kinematic modeling

UCLA/IPAM 2/6/06

modeling how the heart works validation8
Modeling how the heart works: validation

Load independent index of filling

UCLA/IPAM 2/6/06

modeling how the heart works prediction validation
Modeling how the heart works:prediction+validation

Load independent index of filling

Constant slope means

that the response to a

change in peak A-V

gradient is linear.

(data from all 15,

healthy volunteers)

UCLA/IPAM 2/6/06

physiologic interpretation of slope
High load filling regime

Supine load regime

Low load filling regime

Physiologic Interpretation of slope

High load

filling regime

Low load

filling regime

Low slope implies relatively larger increase in viscous loss for the same increase in peak driving force

Higher slope indicates greater efficiency in conversion of initial pressure gradient to attained filling volume.

slide71
Conclusions regarding load independent index
  • Filling patterns change as load is altered, but changed filling patterns obey the same equation of motion. (F=ma)
  • Proposed load independent index M obtainable from non-invasive Doppler Echo.
  • Load independent index is defined by ratio of maximum driving force (kxo peak AV-gradient) to peak viscous force attained (cEpeak).
  • The effect of pathology on M is unknown (so far) - but is predicted to be Mpathologic < M normal
  • M is not expected to be uniquely associated with specific pathology, but will be different from normal.
  • Greatest utility will be in comparing subjects to themselves in response to therapy
summary conclusions
Summary conclusions:
  • Unexplained correlations can be causallyexplained, and
  • new cardiac physiology can be predicted from mathematical
  • modeling and cardiac imaging.
    • E-wave shapes predicted by SHO motion
    • Bi-rectional, linear spring drives filling (TITIN)
    • Constant-volume explains E’/E to LVEDP relation
    • Load Independent index of filling, …

UCLA/IPAM 2/6/06

modeling imaging and function
Modeling, Imaging and Function
  • Unsolved problems remain: (very incomplete listing)
  • Relation between global and segmental indexes of filling
  • What are the eigenvalues of diastolic function
  • Can ‘optimal’ fillling function be defined
    • Relation between model-parameters and biology
    • Relation between model-parameters and pathology
    • Relation between model-parameters and therapy
    • Can you predict ‘stability’ vs ‘instability’ of oscillator?

UCLA/IPAM 2/6/06

slide75
ACKNOWLEDGEMENTS:

NIH

AHA

VETERANS ADMINISTRATION

WHITAKER FOUNDATION

BARNES-JEWISH HOSPITAL FOUNDATION

ALAN A. AND EDITH L.WOLFF CHARITABLE TRUST

UCLA/IPAM 2/6/06

modeling imaging and function1
Modeling, Imaging and Function

THE END

UCLA/IPAM 2/6/06

ad