Open Access Research

Modeling left ventricular diastolic dysfunction: classification and key indicators

Chuan Luo1, Deepa Ramachandran1, David L Ware2, Tony S Ma34 and John W Clark1*

Author Affiliations

1 Dept. Electrical and Computer Engineering, Rice University, Houston, TX 77005, USA

2 Div. Cardiology, University of Texas Medical Branch, Galveston, TX 77555, USA

3 Div. Cardiology, VA Medical Center, Houston, Texas 77030, USA

4 Baylor College of Medicine, One Baylor Plaza, Houston, Texas 77030, USA

For all author emails, please log on.

Theoretical Biology and Medical Modelling 2011, 8:14  doi:10.1186/1742-4682-8-14

Published: 9 May 2011



Mathematical modeling can be employed to overcome the practical difficulty of isolating the mechanisms responsible for clinical heart failure in the setting of normal left ventricular ejection fraction (HFNEF). In a human cardiovascular respiratory system (H-CRS) model we introduce three cases of left ventricular diastolic dysfunction (LVDD): (1) impaired left ventricular active relaxation (IR-type); (2) increased passive stiffness (restrictive or R-type); and (3) the combination of both (pseudo-normal or PN-type), to produce HFNEF. The effects of increasing systolic contractility are also considered. Model results showing ensuing heart failure and mechanisms involved are reported.


We employ our previously described H-CRS model with modified pulmonary compliances to better mimic normal pulmonary blood distribution. IR-type is modeled by changing the activation function of the left ventricle (LV), and R-type by increasing diastolic stiffness of the LV wall and septum. A 5th-order Cash-Karp Runge-Kutta numerical integration method solves the model differential equations.


IR-type and R-type decrease LV stroke volume, cardiac output, ejection fraction (EF), and mean systemic arterial pressure. Heart rate, pulmonary pressures, pulmonary volumes, and pulmonary and systemic arterial-venous O2 and CO2 differences increase. IR-type decreases, but R-type increases the mitral E/A ratio. PN-type produces the well-described, pseudo-normal mitral inflow pattern. All three types of LVDD reduce right ventricular (RV) and LV EF, but the latter remains normal or near normal. Simulations show reduced EF is partly restored by an accompanying increase in systolic stiffness, a compensatory mechanism that may lead clinicians to miss the presence of HF if they only consider LVEF and other indices of LV function. Simulations using the H-CRS model indicate that changes in RV function might well be diagnostic. This study also highlights the importance of septal mechanics in LVDD.


The model demonstrates that abnormal LV diastolic performance alone can result in decreased LV and RV systolic performance, not previously appreciated, and contribute to the clinical syndrome of HF. Furthermore, alterations of RV diastolic performance are present and may be a hallmark of LV diastolic parameter changes that can be used for better clinical recognition of LV diastolic heart disease.