Skip to main content

Minimally invasive, patient specific, beat-by-beat estimation of left ventricular time varying elastance

Abstract

Background

The aim of this paper was to establish a minimally invasive method for deriving the left ventricular time varying elastance (TVE) curve beat-by-beat, the monitoring of which’s inter-beat evolution could add significant new data and insight to improve diagnosis and treatment. The method developed uses the clinically available inputs of aortic pressure, heart rate and baseline end-systolic volume (via echocardiography) to determine the outputs of left ventricular pressure, volume and dead space volume, and thus the TVE curve. This approach avoids directly assuming the shape of the TVE curve, allowing more effective capture of intra- and inter-patient variability.

Results

The resulting TVE curve was experimentally validated against the TVE curve as derived from experimentally measured left ventricular pressure and volume in animal models, a data set encompassing 46,318 heartbeats across 5 Piétrain pigs. This simulated TVE curve was able to effectively approximate the measured TVE curve, with an overall median absolute error of 11.4% and overall median signed error of −2.5%.

Conclusions

The use of clinically available inputs means there is potential for real-time implementation of the method at the patient bedside. Thus the method could be used to provide additional, patient specific information on intra- and inter-beat variation in heart function.

Background

Cardiovascular disease and dysfunction (CVD) are leading causes of Intensive Care Unit (ICU) admission and mortality worldwide. CVD was responsible for 31% of global deaths in 2013, and the estimated worldwide cost associated with CVD in 2010 was $863 billion USD, accounting for approximately 1.39% of gross world product in the same year [1]. Despite the need with aging populations for optimised cardiovascular care in the ICU, inadequate or incorrect diagnosis of cardiac disturbances resulting in increased length of stay, cost and mortality is an ongoing issue [2, 3]. Cardiac management in the ICU is often informed by measurements taken from catheters placed near the heart. However, despite the rich information available from such catheter waveforms, their use is not necessarily associated with improved clinical outcomes [4, 5]. Hence, improvement in the extraction of cardiac information from these catheter waveforms has the potential to yield value from readily available data that has potentially been under-utilised to date.

Time varying elastance (TVE) is an important means of expressing internal cardiac dynamics and function [6]. The TVE curve represents the active elastance changes in the ventricles that drive heart contraction, thus providing valuable intra-beat information about cardiac behaviour and energetics [6–8]. The TVE curve has a wide range of potential applications. It is often employed as an input in large cardiovascular models that seek to model the entire circulatory system [9–13], and also frequently used in the determining of important cardiac metrics, such as the end-systolic pressure–volume relation (ESPVR) [7].

There has been more limited investigation into the utility of the TVE curve in isolation. The area under the curve of the TVE curve is analogous to work done by the ventricle. It is thus a potentially significant indicator of patient condition. However, the typical normalisation of the TVE curve means only relative changes in work for a given inotropic state can be compared. The TVE curve relies on a combination of ventricular pressure and volume waveforms similar to pressure–volume (P–V) loops. It thus contains information concerning cardiac work [14, 15] and contractility [6, 16], both of which change in response to cardiac dysfunction. As such, it has been suggested that the shape of the ventricular elastance curve itself may have diagnostic use [17–19].

Unfortunately, the TVE curve can only be directly measured by placing catheters into the heart chambers, which is, understandably, not common practice. As such, the traditional approach to non-invasively generate a TVE curve is to fix it to a population based waveform [7, 11, 12]. Most existing work surrounding TVE curves is focused on using the TVE curve to estimate a clinical parameter, such as ejection fraction [20], and, most commonly, end-systolic elastance [7, 21, 22]. However, these studies are validated on correlations with the derived parameter, rather than on the shape and change in shape of the TVE curve itself. Thus, these approaches are not validated for use as part of a larger model of cardiac dynamics, or for direct use as a diagnostic aid.

Previous work specifically focused on experimentally generating a TVE curve and validating it based on its correlation with the analytically derived function showed promise, and noted changes in the TVE curves during pulmonary embolism and septic shock [17]. However, this work was limited by the availability of data for validation and the reliance on an assumed driver shape [17, 18]. Work has also been undertaken in modelling time varying ventricular elastance, split into active and passive components, in humans [19]. Active elastance was shown to compare well with other metrics of contractility, and the properties of these elastance curves were shown to change for different disease states. However, this method required highly invasive ventricular catheterisation and thus is not broadly implementable in a clinical environment. Hence, there is a significant gap created by the current clinical inability to directly measure or estimate the TVE curve every beat.

This paper presents a novel, minimally invasive method for deriving the TVE curve beat-by-beat. The method focuses on combining simple physiological assumptions with clinically available catheter waveforms to individually simulate the pressure and volume components that define the TVE curve, rather than generating the TVE curve itself. Importantly, this approach avoids the need to directly assume a shape for the TVE curve, and is thus better equipped to capture variations in this shape over time and condition, as well as the corresponding alterations in intra-beat cardiac behaviour.

Clinically feasible measurements mean the method has the potential for real-time implementation at the patient bedside, without requiring additional, invasive instrumentation. Such a TVE curve could be used to provide additional, patient specific information on intra-beat behaviour and inter-beat variation in the functioning of the heart. Finally, monitoring the driver’s evolution over time could add significant new data and insight to improve diagnosis and treatment.

Methods

Proposed method

The TVE curve is defined:

$$e(t) = \frac{{P_{lv} \left( t \right)}}{{V_{lv} \left( t \right) - V_{d} }}$$
(1)

where P lv is the pressure in the left ventricle, V lv is the volume in the left ventricle and V d is the ‘dead space’ volume in the ventricle [6, 23]. Thus, the TVE curve is defined by two waveforms (P lv and V lv ) and one constant (V d ). These values can be measured directly, but doing so is not clinically feasible [24].

The proposed method approximates these waveforms (V lv , P lv ) and constant (V d ) using three inputs, as shown in Fig. 1. These inputs include continuously sampled aortic pressure waveforms (P ao ) and heart rate (HR), data which is typically available in a modern ICU. The final input required is baseline end-systolic (V es ) and end-diastolic (V ed ) Volume, obtained from a brief echocardiography reading, which is increasingly available in a clinical setting [25].

Fig. 1
figure 1

Flowchart of the proposed method. Roman numerals indicate different method regions

The overall goal of this method is to use the clinically available inputs P ao , HR and baseline V es and V ed to determine the outputs P lv , V d and V lv , and thus the TVE curve e(t) as set out in Fig. 1 and Eq. 1. This resulting TVE curve can be experimentally validated against the TVE curves derived from experimentally measured P lv and V lv in animal models. The availability of a nearby continuous pressure measurement (P ao ) forms an effective basis for the continuous approximation of P lv , and the availability of a baseline volume measurement (V es ) forms an effective basis for the approximation of baseline V d . However, there is no continuous volume measurement available for continuous approximation of Vlv, resulting in this process being considerably more involved. Thus, the shaded central Region III in Fig. 1 is considerably more complicated than Regions I and II.

While the overall method involves approximating two output waveforms (V lv , P lv ) from a single input waveform (P ao ), it’s important to note that all three waveforms (P lv , V lv , P ao ) have distinct features. Different regions of behaviour are governed by different physiological phenomenon, and these waveforms have been extensively characterised [26]. As such, all three waveforms are heavily interconnected and information rich, making this task more reasonable than it might first appear.

Determining P lv from P ao (Region I, Fig. 1)

The left ventricle is situated directly upstream from the aorta, separated by the aortic valve. This valve is open during systole and closed during diastole. As such, if aortic valve resistance is neglected, P lv is equivalent to P ao , with a slight phase lag (δ) during the majority of systole (Section P. 1, Fig. 2). While aortic valve resistance is non-negligible in conditions such as aortic stenosis [26], valve dysfunction of any type is present in only 3.61% of CVD mortalities in the US [27]. Further, such conditions are typically chronic in nature, relatively easily diagnosed [28] and evolve slowly, while the method is designed to monitor short term changes in an ICU environment. As such, this assumption should not significantly affect the methods’ applicability to the vast majority of the target cohort.

Fig. 2
figure 2

Simulating left ventricular pressure. Note P ao has been shifted left by δ to account for phase lag

During diastole the aortic valve is closed and little information is available from the aortic waveform about ventricular behaviour. However, the ventricle behaves in a largely passive manner in this region, meaning the TVE curve is typically near zero during diastole [26]. As such, a generic function consisting of two exponentials was used to approximate left ventricular pressure during diastole. In early diastole (Section P. 2, Fig. 2) an exponential decay to a fixed baseline pressure captures ventricular relaxation. In late diastole-early systole (Section P. 3, Fig. 2), an exponential increase captures the beginning of ventricular contraction [26].

While atrial contraction contributes significantly to late diastolic filling, ventricular elastance remains largely passive in late diastole [26]. Thus the TVE curve is typically at its baseline value until the beginning of ventricular contraction in early systole. While atrial function affects the magnitude of the driver function, which is normalised in this work, and may affect its shape, this effect is indirect as the driver function represents the impact of contraction in driving pulsatile blood from the heart to the arterial system. Further, with the increasing unpopularity of pulmonary artery catheters [29] none of the typically available instrumentation in an ICU provides a clear picture of atrial behaviour. As such, while the exponential in section P. 3 is broadly intended to capture ventricular filling, no specific atrial behaviour component is integrated into this model.

Using Fig. 2, the left ventricular pressure for the n th heartbeat is thus defined using P ao :

$$t_{1} = t\left( {\frac{{dP_{ao} }}{dt}_{max} } \right)_{n}$$
(2a)
$$t_{2} = t\left( {\frac{{dP_{ao} }}{dt}_{min } } \right)_{n}$$
(2b)
$$t_{4} = t\left( {\frac{{dP_{ao} }}{dt}_{max } } \right)_{n + 1}$$
(2c)
$$t_{3} = 0.62t_{2} + 0.38t_{4}$$
(2d)
$$P_{lv} (t) = \left\{ {\begin{array}{*{20}l} {P_{ao} \left( {t_{1} + \delta < t < t_{2} + \delta } \right)} \hfill & {t_{1} < t < t_{2} } \hfill \\ {6 + \left( {P_{ao} \left( {t_{2} } \right) - 6} \right)e^{{ - 17.5\left( {t - t_{2} } \right)}} } \hfill & {t_{2} < t < t_{3} } \hfill \\ {P_{lv} \left( {t_{3} } \right) + \left( {P_{ao} \left( {t_{4} } \right) - P_{lv} \left( {t_{3} } \right)} \right)e^{{37.5\left( {t - t_{4} } \right)}} } \hfill & {t_{3} < t < t_{4} } \hfill \\ \end{array} } \right\}$$
(3)

where: \(\delta = 0.008{\text{s}}\)

Determining V d from baseline V es (Region II, Fig. 1)

A recent method has been developed to approximate V d from ventricular volume measurements [30]. This approach relies on linear regression of the Frank-Starling curve (SV–V ed ) and its end-systolic equivalent (SV–V es ) to the point where SV = 0 and ‘the ventricle cannot develop any systolic pressure’, the definition of V d [23]. This work also showed V d for a baseline, healthy pig be an approximately fixed percentage of V es [30]. Defining V d as a percentage of baseline V es allows approximation of baseline V d during the initial echocardiographic reading where measured V es is available. While V d has been shown to change with condition, there is no practical means of capturing this change short of additional echocardiography measurements. Thus, while intermittent measures are feasible, in this study V d is fixed at a baseline value.

$$V_{d} = 0.48 \times V_{es}$$
(4)

Determining V lv from P ao , HR and V d (Region III, Fig. 1)

Unlike pressure, little volume or flow information is readily available from the typical, clinically available instrumentation. Simulating V lv is thus more challenging. The shape of the ventricular volume waveform was approximated by a piecewise sine wave consisting of two sections: systole (Section V. 1, Fig. 3) and diastole (Section V. 2, Fig. 3), with a 90° phase shift at the beginning of systole. The underlying physiological behaviour might be better represented by a series of exponentials [26]. However, using sine waves achieves a similar result with considerably fewer variables involved.

Fig. 3
figure 3

Simulating left ventricular volume

Thus, six points per heartbeat (t 1, t 2, t 3 and (V ed ) n , (V es ) n , (V ed ) n+1 ) are required to define the ventricular volume waveform. The timing associated with systole start (t 1), systole end (t 2), and diastole end (t 3) are readily determined from the aortic pressure waveform (Fig. 3):

$$t_{1} = t\left( {P_{ao_{min}} } \right)_{n}$$
(5a)
$$t_{2} = t\left( {P_{DN} } \right)_{n}$$
(5b)
$$t_{3} = t\left( {P_{ao_{min}} } \right)_{n + 1}$$

Using existing work [31] SV can be approximated beat-to-beat using the aortic waveform. Thus, only one of V es or V ed is required, as SV can be used to convert between the two. The ESPVR allows for determination of V es , and is defined [23]:

$$P_{es} = E_{es} \times \left( {V_{es} - V_{0} } \right)$$
(6)

where E es is the end-systolic elastance and V 0 is the ventricular volume at zero pressure. Equation 6 can be rewritten:

$$P_{DN} = E_{es} \times \left( {V_{es} - V_{d} } \right)$$
(7)

where this change is justified by:

  • The pressures in the ventricle and aorta are roughly equivalent until the aortic valve closes, thus P DN is close to P es

  • V d and V 0 have similar, but distinct, physiological significance and values. The two are often used interchangeably [17, 18, 23]

Finally, it is necessary to account for E es , which changes in response to a number of factors including contractility [6], and loading conditions [32, 33]. Thus, Eq. 7 is modified:

$$P_{DN} = \left( {E_{c} \times HR^{3} } \right) \times \left( {V_{es} - V_{d} } \right)$$
(8)

here E es is defined as a function of HR and a coefficient (E C ), with a cubic selected as it provides the best compromise between simplicity and effective tracking for the data set presented here. In particular, the cardiovascular system responds to most changes in conditions in a number of ways, including changes in heart rate and elastance. As heart rate is easily measured, it provides an easy to obtain, if incomplete, indication of cardiovascular system response, which can be used to inform an approximated elastance [34]. Further supporting evidence is provided in the validation and discussion of results.

During the echocardiography calibration, measurements for P DN , HR and V es are available [35]. Thus, using Eq. 8, a constant value for E C can be defined, allowing approximation of E es and thus determination of V es on a beat-by-beat basis. The beat-to-beat ventricular volume can thus be determined:

$$V_{lv} (t) = \left\{ {\begin{array}{*{20}l} {\left( {V_{ed} } \right)_{n} + \left( {\left( {V_{es} } \right)_{n} - \left( {V_{ed} } \right)_{n} } \right)\sin \left( {\frac{{\pi \left( {t - t_{1} } \right)}}{{2\left( {t_{2} - t_{1} } \right)}}} \right)} \hfill & {t_{1} < t < t_{2} } \hfill \\ {\left( {V_{es} } \right)_{n} - \left( {\left( {V_{ed} } \right)_{n + 1} - \left( {V_{es} } \right)_{n} } \right)\left( {\frac{1}{2}\cos \left( {\frac{{\pi \left( {t - t_{2} } \right)}}{{\left( {t_{3} - t_{2} } \right)}}} \right) - \frac{1}{2}} \right)} \hfill & {t_{2} < t < t_{3} } \hfill \\ \end{array} } \right\}$$
(9)

where:

$$V_{es} = \frac{{P_{DN} }}{{\left( {E_{c} \times HR^{3} } \right)}} + V_{d}$$
(10)
$$V_{ed} = V_{es} + SV$$
(11)

Summary of proposed method

The overall derivation of the TVE curve can be summarised:

Initially or intermittently

  1. 1.

    Calculate V d using Eq. 4 and baseline V es (Region I, Fig. 1)

  2. 2.

    Calculate E C using Eq. 8, P DN , HR and baseline V es (Region III, Fig. 1)

Every heartbeat

  1. 1.

    Simulate P lv using Eq. 2, Eq. 3 and P ao (Region II, Fig. 1)

  2. 2.

    Determine V es using Eq. 10, P DN , HR and E C (Region III, Fig. 1)

  3. 3.

    Determine SV using [31] and P ao (Region III, Fig. 1)

  4. 4.

    Determine V ed using Eq. 11, V es and SV (Region III, Fig. 1)

  5. 5.

    Simulate V lv using Eq. 5, Eq. 9, P ao , V es and V ed (Region III, Fig. 1)

  6. 6.

    Calculate and normalise the TVE curve e(t) using Eq. 1

Analysis and validation

The proposed method was validated on experimentally gathered data. A range of input (P ao ) and output (V lv , P lv ) waveforms were continuously measured via catheter. This data allowed validation of individual model assumptions through comparison with directly measured output waveforms, as well as validation of the overall method through comparisons between the TVE curve as calculated using simulated and directly measured V lv and P lv wave forms.

The data set encompasses 46,318 heartbeats across 5 Piétrain pigs. A diverse clinical protocol provides the ability to assess intra- and inter-subject variability across a large amount of invasive and non-invasive measurements. Together they enable rigorous assessment and validation of the method.

Experimental procedure

The experimental protocol was approved by the Ethics Commission for the Use of Animals at the University of Liège, Belgium. Five male, pure Piétrain pigs weighing between 18.5 and 29 kg were sedated, anaesthetised and mechanically ventilated (GE Engstrom CareStation) with a baseline positive end-expiratory pressure (PEEP) of 5 cmH2O (Fig. 4, Additional file 1: Table S1). Proximal aortic pressure was continually sampled using a pressure catheter (Transonic, NY, USA) with a sampling rate of 250 Hz. To provide direct measurements of P lv and V lv for validation, the heart was accessed via a median sternotomy, and an admittance pressure–volume catheter (Transonic, NY, USA) with a sampling rate of 250 Hz inserted into the left ventricle via an apical stab [36, 37].

Fig. 4
figure 4

The fully instrumented Piétrain Pig. The catheters for measuring P ao , P lv and V lv are positioned to the left of the image

To demonstrate a diverse range of cardiac states, several procedures were performed:

  • A single infusion of endotoxin (lipopolysaccharide from E. Coli, 0.5 mg/kg injected over 30 min) to induce septic shock. Septic shock drives a change in afterload conditions and is associated with a large variety of effects including an inflammatory response and capillary leakage that may lead to hypovolemia, decreased cardiac output, decreased ejection fraction and cardiac failure [38].

  • Several PEEP driven recruitment manoeuvres (RMs), both pre- and post-endotoxin infusion. RMs drive a change in preload conditions and are typically associated with a decrease in mean blood pressure and cardiac output [39].

  • One to four infusions of 500 mL saline solution over 30 min, pre- and post-endotoxin infusion, simulating fluid resuscitation therapy, a key component of hemodynamic resuscitation in patients with severe sepsis, which itself results in a change in circulatory volume [40].

Validation of significant model assumptions

Two major assumptions made in deriving this model are:

  1. 1.

    That V d can be expressed as a function of baseline V es , as in Eq. 4

  2. 2.

    That E es can be expressed as a function of HR, as in Eq. 8

Direct evaluation of the tracking of V es using different forms of ESPVR allows validation of both of these assumptions. In particular, 3 different methods of tracking V es were compared:

  1. 1.

    Fixed E es and neglected V 0 : The standard ESPVR (Eq. 6) with V 0 = 0 (a commonly used assumption [17, 18, 23])

  2. 2.

    Fixed E es and fixed V 0 : The standard ESPVR (Eq. 7) with V 0 = V d (allowing assessment of the validity of Eq. 4)

  3. 3.

    Dynamic E es and fixed V 0 : The ESPVR as used in the proposed method (Eq. 8) with V 0 = V d , and E es as a function of HR (allowing assessment of the validity of Eq. 8)

Validation of overall model

The overall method presented here is designed to simulate the TVE curve beat-by-beat, without requiring invasive instrumentation of the heart or real-time image-based monitoring, neither of which is clinically or ethically feasible in care. As such, validation of the method relies on comparison of the simulated TVE curve to the invasively measured, ‘true’ TVE curve, which is calculated using the catheter measured V lv and P lv waveforms for a single beat. This comparison is achieved by calculating the absolute and signed ‘error area’ between the measured and simulated TVE curve, according to Eqs. 12 and 13:

$$\varepsilon_{abs} = \frac{{\mathop \smallint \nolimits_{t = 0}^{1} \left| {e_{sim} (t) - e_{meas} (t)} \right|}}{{\mathop \smallint \nolimits_{t = 0}^{1} \left( {e_{meas} (t)} \right)}}$$
(12)
$$\varepsilon_{\text{sgn}} = \frac{{\mathop \smallint \nolimits_{t = 0}^{1} \left( {e_{sim} (t) - e_{meas} (t)} \right)}}{{\mathop \smallint \nolimits_{t = 0}^{1} \left( {e_{meas} (t)} \right)}}$$
(13)

where ε sim and ε meas are the simulated and measured TVE curves respectively, t is normalised time set to 1 for every heart beat to enable comparison over different beats, and ε abs and ε sgn and denote the absolute and signed errors respectively. An example TVE curve with an absolute error of 7.8% and bias of −3.4% is shown in Fig. 5, where shading denotes the error area.

Fig. 5
figure 5

Example TVE curve error. This simulated driver has an error of 7.8% and bias of −3.4%

Results

Validation of significant model assumptions

Table 1 shows the percentage error associated with the 3 methods specified to track V es . Moving from neglected V 0 (Method 1) to fixed V 0 (Method 2), via Eq. 4, shows a modest reduction in median and 25th percentile errors (15.9–14.5% and 3.1–1.3%, respectively), but a highly significant reduction in 75th percentile error (43.2–25.5%). Introduction of dynamic E es (Method 3) as opposed to fixed E es (Method 2), via Eq. 8, shows a very significant reduction in median and 75th percentile errors (14.5–4.3% and 25.5–15.0%, respectively). These results support the validity and usefulness of Eqs. 4 and 8, their associated assumptions in the tracking of V es , and thus the use of Method 3 in the overall approximation of the TVE curve.

Table 1 Absolute percentage error associated with different methods of tracking Ves

Note the values quoted from the final row of Table 1 are an average of the 25th percentile, median and 75th percentile error for the pigs rather than an overall 25th percentile, median and 75th percentile error. The decision was made to represent population values in this way as the overall median is less representative than the average median due a differing number of heartbeats recorded for different pigs.

The comparative behaviour of these 3 methods is well illustrated in Fig. 6, which shows the tracking of V es for Pig 1. Method 1 agrees reasonably well with measured V es during normal behaviour, but diverges significantly from measured V es during recruitment manoeuvres and the onset of severe sepsis, where there are large vertical changes in V es . Method 2 accurately captures these large vertical changes in V es , reducing the large 75th percentile error in Method 1. Method 3 retains the tracking of large vertical changes in Method 2, and tracks normal behaviour more effectively by accounting for changes in E es , significantly reducing median error.

Fig. 6
figure 6

Tracking of V es . Comparison between the three methods in 2.2.2 for Pig 1

Validation of overall model

The overall purpose of this model is to track the shape of the TVE curve, and how this shape changes when circulatory behaviour changes. Table 2 shows the area under the curve errors, ϵ abs and ϵ sgn , for all 5 pigs. The median error is relatively small at 11.4%, suggesting the method is effective, and the interquartile range relatively narrow at 9.2–14.7%, suggesting the method is consistent. The bias is also small at −2.5%, with an interquartile range of −6.1 to 0.8%, suggesting only a slight method bias.

Table 2 TVE curve percentage area under the curve errors (ε abs and ε sgn ) associated with proposed method (identifying V es and V ed )

Table 3 shows, as a point of comparison, the area under the curve errors associated with the method if the directly measured values for V es and V ed are used. Using these measured values removes the majority of Region III, Fig. 1 (Eqs. 4 and 6–11) from the overall method, where a lot of relatively significant assumptions are made. The errors in Table 3 are very comparable to those in Table 2, with an overall modest reduction in average median error (11.4–10.2%) and bias (−2.5 to −2.2%), as expected when measuring V es and V ed directly and invasively rather than estimating them. Hence, the assumptions made had very little impact on error when removed.

Table 3 TVE curve percentage area under the curve error (ε abs and ε sgn ) associated with using measured V es and V ed values

Figure 6 shows a variety of measured and modelled TVE curve shapes for each of the 5 pigs, capturing both intra- and inter-subject variability. For each figure, Panel I shows a baseline waveform, Panel II a waveform during a RM and Panel III a waveform after endotoxin infusion has occurred and sepsis develops. The range of shapes and inter-pig variability indicated that assuming a generic population or cohort TVE curve is not particularly valid, and is further supported by the range of intra-pig variability evidenced. There are a variety of shoulder heights, relative gradient and maximum/minimum gradient timings that result in TVE curves with distinct shapes, variations which are well captured by the method. These TVE curves were also selected to demonstrate a range of error values similar to the interquartile error range for that pig (Table 2). Thus, these drivers are reasonably representative of the overall ability of the method to capture TVE curve shapes, despite different overall shapes for each pig and changes in shape as condition changes.

Discussion

Validation of significant model assumptions (Table 1)

A direct evaluation of the effectiveness of three different ESPVR equations set out in 2.2.2 in tracking V es allowed validation of the assumptions that:

  • That V d can be derived from V es as in Eq. 4 (implemented in Method 2)

  • That E es can be expressed as a function of HR as in Eq. 8 (implemented in Method 3)

Method 1 serves as a control and uses a simplified ESPVR, assuming V 0 = 0 and E es is constant. This method tracks V es reasonably well, yielding an overall median error of 15.9% across all 5 pigs compared to V es as directly measured. The assumption that V 0 can be neglected is often used due to a combination of V 0 being difficult to measure directly, as it requires a significant artificial reduction in ventricular pressure [41], and V 0 typically being relatively small [42]. This assumption is largely supported by these results, as a 15.9% median error seems acceptable when weighed against the type of highly invasive and involved protocol traditionally required to determine V 0.

However, as shown in Fig. 6, this 15.5% median error fails to capture the extremely large inaccuracies associated with using Method 1 to track sudden changes in V es due, for example, to recruitment manoeuvres or the onset of severe sepsis. This failure to track sudden changes leads to a significantly larger 75th percentile error of 43.2% for Method 1. This high error when sudden changes occur is of concern in an ICU or cardiac surgery clinical scenario, where sudden changes and accurate, rapid determination of patient responses to these sudden changes is extremely important [26].

Method 2 introduces the assumption that V d can be derived from baseline V es , as in Eq. 4, and that V d can be used as a surrogate for V 0. There is a minor, but important physiological distinction between the two values: V 0 is the ventricular volume at 0 pressure, while V d is the volume at which the ventricle cannot develop any systolic pressure [23]. However, the purpose of both terms is to account for the subject-specific inactive volume within the ventricle, and the two values have been shown to be similar for a given subject [23].

Method 2 results in a notable decrease in 75th percentile error compared to Method 1 (43.2–25.5%), and is able to track sudden changes in V es significantly more effectively, as shown in Fig. 6. However, Method 2 yields only a modest reduction in median error compared to Method 1 (15.9–14.5%), as both methods fail to capture the dynamic nature of E es . Thus, large errors are reduced, but overall accuracy is not greatly improved. Regardless, these results provide support for the validity of Eq. 4, and the use of V d as a surrogate for V 0 here. The fact that these reductions in error are sustained over significant changes in cardiac output and ejection fraction as sepsis develops suggests that the absolute value of V d does not change significantly enough under such conditions to detract from method accuracy.

Method 3 further assumes that E es can be expressed as a function of HR. This assumption is unusual, but is supported by the results. The full method sees a further, significant reduction in 75th percentile error compared to Methods 1 and 2 (43.2% and 25.5–15.0%, Table 1), and, most importantly, a very large reduction in median error compared to both Methods 1 and 2 (15.9% and 14.5–4.3%, Table 1). This result suggests general tracking of trends in V es is being significantly improved, also supporting the validity of expressing E es as a function of HR. This behaviour can also be observed in Fig. 6. This result, combined with the minimal addition in method complexity required to include HR, which is very easy to measure, provides a strong case for the use of Method 3.

However, it is still important to note that the relationship between E es and HR expressed in Eq. 8 is a significant simplification of actual cardiac behaviour. The cardiac system uses a large variety of responses to maintain cardiac output. The cubic used to approximate E es changes as a function of HR attempts to mathematically approximate the sympathetic nature of some of these responses, but, inevitably, the relationship between HR and E es varies between subjects (accounted for by calibration), and as time and condition changes. While the cubic approximation of this relationship has been demonstrated to remain effective across the full progression of septic shock in the data presented here, further validation of this relationship, and especially the use of a cubic, across other conditions is desired.

Validation of overall model

Simulating the TVE curve using the proposed method resulted in a relatively low median absolute error area, ranging from 9.4 to 13.4% (Table 2) across all pigs. This narrow range of median absolute errors implies the method is able to consistently and effectively capture inter-subject variations in TVE curve behaviour, suggesting it is generalizable to other subjects. This inter-subject variability is demonstrated well by various baseline drivers shown in Fig. 7. All 5 pigs demonstrate relatively distinct baseline TVE curve shapes, showing a considerable level of inter-subject variability impossible to capture with a generic cohort or population based TVE curve that relies on a basic assumed shape [17, 18].

The simulated TVE curve also had a relatively narrow average interquartile error range of 9.2–14.7% (Table 2), suggesting the method is able to relatively consistently and effectively capture inter-subject variations in TVE curve behaviour. This intra-subject variability, and inter-quartile error range, is illustrated in Fig. 7. Accurate capture of relative changes and trends in patient behaviour is extremely important, as these changes are critical in assessing whether a patient is recovering, responding to treatment, or is in need of a change in treatment. Most pigs demonstrated notable changes and intra-pig variability in driver shape during recruitment manoeuvres and as sepsis developed. A couple of the more unusual driver shapes were not fully captured, for example Panel 2 of Pig 4 in Fig. 7, where the driver measured driver displayed two peaks and the simulated driver only one.

Fig. 7
figure 7

Example TVE curves for each pig. A range of error values and cardiac states are shown

The simulated TVE curve demonstrated consistently low signed area error, at −2.5% (IQR −6.1 to 0.8%), shown in Table 2. This result shows the method only slightly underestimating the TVE curve, supported further by 2 pigs having a slight positive bias and the other 3 a slight negative bias. This outcome further supports the ability of the TVE curve to accurately capture both intra- and inter-subject variability over time and condition.

In assessing the impact of assumptions on the TVE curve, a comparison of Table 2, using approximated V es and V ed , and Table 3, using measured V es and V ed , show very similar error values. For example, overall median error fell from 11.4 to 10.2% and overall bias from −2.5 to −2.2%, only a modest reduction in error. This implies that the body of assumptions and equations in Region III, Fig. 1 concerned with the approximating the V es and V ed for simulation of V lv do not result in a large increase in error compared to using measured V es and V ed . The assumptions made in areas of the method not involving simulating V es and V ed are relatively minimal, mostly involving using P ao to determine waveform timing. As such, it would seem that much of the error associated with this method is the result of the necessity of assuming equations for parts the two waveforms being reconstructed (P lv and V lv ), as shown in Figs. 2 and 3. Thus, the assumptions employed in Region III, Fig. 1 appear to function as intended, and further reduction in error would probably require increased method complexity or an increase in the clinically available data.

An important point to consider is the fact that the TVE curve is consistently normalised to a duration and amplitude of 1.0, as it is designed specifically as an indicator of how the heart is behaving relatively over the course of a beat, to be coupled with a lumped metric (E es ) indicating the overall strength of that heartbeat [6, 17, 18]. This normalisation does mean that some of the errors associated with the various assumptions and approximations made throughout the method are negated, and that indicators of absolute cardiac work and its changes, for example, are not able to be directly extracted from the TVE curve created. Fortunately, the driver does not exist alone, the intent is that the shape of the TVE curve, indicative of transient, relative cardiac behaviour, be used alongside other existing metrics, such as Cardiac Output or Stroke Volume [31, 43], indicative of lumped, absolute cardiac behaviour, to provide further diagnostic information.

Limitations

There are study limitations that should be considered. First, all data presented is derived from a single protocol involving a single, but complex and varied [38], condition (sepsis). This data set encompasses several pigs, a full progression from healthy, baseline behaviour to cardiac failure and clinically standard ventilation and fluid interventions. Nevertheless, there is a much larger range of possible cardiac conditions, and further validation over several of these would be beneficial. For example, the method would benefit from validation on contractility altering drugs such as dobutamine [44], which may alter the behaviour of the elastance term in Eq. 8. However, the method already detects changes in haemodynamics, including those due to circulatory or cardiac muscle changes during sepsis in this study. Thus, the ability to detect changes due to inotropes should be similar to what is presented for the range of behaviours already observed in sepsis. Further, given inotrope infusions are determined and performed by a clinician, it would be possible to recalibrate the method directly after such an infusion to adjust to the new inotropic state if it proved necessary. Overall, the underlying physiology and data supporting the development of this method has been discussed in detail, and would be expected to generalise well to a wider range of conditions, as there are no intervention or condition specific assumptions made.

The method also requires validation on human subjects to ensure the methodology as presented here remains physiologically accurate, though the strong similarities between porcine and human physiology and the effectiveness of porcine models are well established [45, 46]. Equally, only an animal model, as used here, allows the direct validation against cardiac measured PV loops, which would not be possible in humans. Thus, only an animal trial allows this important first validation.

The method does require an initial calibration via echocardiography or similar means. Echocardiography equipment is increasingly available in modern ICUs [25]. Further, echocardiography is non-invasive and the calibration period required is relatively short, requiring approximately 10 heartbeats. However, the requirement of such a calibration still prevents the method from being fully implementable without modest additional clinical workload using normal ICU instrumentation.

Conclusion

The TVE curve is an important, but difficult to clinically measure, expression of internal cardiac dynamics that captures the heart’s ability as a pump and can evolve over time, condition and patients. A novel, minimally invasive method for deriving the TVE curve beat-by-beat, by combining simple physiological assumptions with readily available catheter waveforms to individually simulate the components of the TVE curve, is proposed. This method was assessed across a cohort of 5 Piétrain pigs undergoing a progression from healthy behaviour to cardiac failure due to sepsis. The TVE curve generated by the method was shown to effectively track a directly measured function throughout the experiments, with low overall median absolute (11.4%) and signed (−2.5%) area under the curve errors. There is the potential for this method to provide real time, patient specific information on intra-beat behaviour and inter-beat variation in the heart, at the patient bedside, without requiring additional, invasive instrumentation.

Abbreviations

CVD:

cardiovascular disease and dysfunction

ICU:

Intensive Care Unit

ESPVR:

end-systolic pressure–volume relation

P–V:

pressure–volume

PEEP:

positive end-expiratory pressure

RM:

recruitment manoeuvre

References

  1. Mozaffarian D, Benjamin EJ, Go AS, Arnett DK, Blaha MJ, Cushman M, Das SR, de Ferranti S, Després J-P, Fullerton HJ. Heart Disease and Stroke Statistics—2016 update: a report from the American Heart Association. Circulation. 2015;. doi:10.1161/CIR.0000000000000152.

    Google Scholar 

  2. Angus DC, Linde-Zwirble WT, Lidicker J, Clermont G, Carcillo J, Pinsky MR. Epidemiology of severe sepsis in the United States: analysis of incidence, outcome, and associated costs of care. Crit Care Med. 2001;29(7):1303–10.

    Article  Google Scholar 

  3. Pineda LA, Hathwar VS, Grant BJ. Clinical suspicion of fatal pulmonary embolism. CHEST J. 2001;120(3):791–5.

    Article  Google Scholar 

  4. Frazier SK, Skinner GJ. Pulmonary artery catheters: state of the controversy. J Cardiovasc Nurs. 2008;23(2):113–21.

    Article  Google Scholar 

  5. Chatterjee K. The swan-ganz catheters: past, present, and future a viewpoint. Circulation. 2009;119(1):147–52.

    Article  Google Scholar 

  6. Suga H, Sagawa K, Shoukas AA. Load independence of the instantaneous pressure–volume ratio of the canine left ventricle and effects of epinephrine and heart rate on the ratio. Circ Res. 1973;32(3):314–22.

    Article  Google Scholar 

  7. Senzaki H, Chen C-H, Kass DA. Single-beat estimation of end-systolic pressure–volume relation in humans a new method with the potential for noninvasive application. Circulation. 1996;94(10):2497–506.

    Article  Google Scholar 

  8. Sunagawa K, Sagawa K, Maughan WL. Ventricular interaction with the loading system. Ann Biomed Eng. 1984;12(2):163–89.

    Article  Google Scholar 

  9. Revie JA, Stevenson D, Chase JG, Pretty CJ, Lambermont BC, Ghuysen A, Kolh P, Shaw GM, Desaive T. Evaluation of a model-based hemodynamic monitoring method in a porcine study of septic shock. Comput Math Methods Med. 2013;2013:505417.

    Article  Google Scholar 

  10. Revie JA, Stevenson DJ, Chase JG, Hann CE, Lambermont BC, Ghuysen A, Kolh P, Shaw GM, Heldmann S, Desaive T. Validation of subject-specific cardiovascular system models from porcine measurements. Comput Methods Programs Biomed. 2013;109(2):197–210.

    Article  Google Scholar 

  11. Chung D, Niranjan S, Clark J, Bidani A, Johnston W, Zwischenberger J, Traber D. A dynamic model of ventricular interaction and pericardial influence. Am J Physiol Heart Circ Physiol. 1997;272(6):H2942–62.

    Google Scholar 

  12. Starfinger C, Chase J, Hann C, Shaw G, Lambermont B, Ghuysen A, Kolh P, Dauby P, Desaive T. Model-based identification and diagnosis of a porcine model of induced endotoxic shock with hemofiltration. Math Biosci. 2008;216(2):132–9.

    Article  MathSciNet  MATH  Google Scholar 

  13. Smith BW, Chase JG, Nokes RI, Shaw GM, Wake G. Minimal haemodynamic system model including ventricular interaction and valve dynamics. Med Eng Phys. 2004;26(2):131–9.

    Article  Google Scholar 

  14. Suga H. Ventricular energetics. Physiol Rev. 1990;70(2):247–77.

    Google Scholar 

  15. Burkhoff D, Sagawa K. Ventricular efficiency predicted by an analytical model. Am J Physiol Regul Integr Comp Physiol. 1986;250(6):R1021–7.

    Google Scholar 

  16. Broscheit J-A, Weidemann F, Strotmann J, Steendijk P, Karle H, Roewer N, Greim C-A. Time-varying elastance concept applied to the relation of carotid arterial flow velocity and ventricular area. J Cardiothorac Vasc Anesth. 2006;20(3):340–6.

    Article  Google Scholar 

  17. Stevenson D, Revie J, Chase JG, Hann CE, Shaw GM, Lambermont B, Ghuysen A, Kolh P, Desaive T. Beat-to-beat estimation of the continuous left and right cardiac elastance from metrics commonly available in clinical settings. Biomed Eng Online. 2012;11:73.

    Article  Google Scholar 

  18. Stevenson D, Revie J, Chase JG, Hann CE, Shaw GM, Lambermont B, Ghuysen A, Kolh P, Desaive T. Algorithmic processing of pressure waveforms to facilitate estimation of cardiac elastance. Biomed Eng Online. 2012;11(1):1–16.

    Article  Google Scholar 

  19. Zhong L, Ghista DN, Ng EY, Lim ST. Passive and active ventricular elastances of the left ventricle. Biomed Eng Online. 2005;4(1):10.

    Article  Google Scholar 

  20. Swamy G, Kuiper J, Gudur MS, Olivier NB, Mukkamala R. Continuous left ventricular ejection fraction monitoring by aortic pressure waveform analysis. Ann Biomed Eng. 2009;37(6):1055–68.

    Article  Google Scholar 

  21. Ten Brinke E, Klautz R, Verwey H, Van Der Wall E, Dion R, Steendijk P. Single-beat estimation of the left ventricular end-systolic pressure–volume relationship in patients with heart failure. Acta Physiol. 2010;198(1):37–46.

    Article  Google Scholar 

  22. Shishido T, Hayashi K, Shigemi K, Sato T, Sugimachi M, Sunagawa K. Single-beat estimation of end-systolic elastance using bilinearly approximated time-varying elastance curve. Circulation. 2000;102(16):1983–9.

    Article  Google Scholar 

  23. Sagawa K. Editorial: the end-systolic pressure–volume relation of the ventricle: definition, modifications and clinical use. Circulation. 1981;63(6):1223–7.

    Article  Google Scholar 

  24. Kastrup M, Markewitz A, Spies C, Carl M, Erb J, Grosse J, Schirmer U. Current practice of hemodynamic monitoring and vasopressor and inotropic therapy in post-operative cardiac surgery patients in Germany: results from a postal survey. Acta Anaesthesiol Scand. 2007;51(3):347–58.

    Article  Google Scholar 

  25. Vieillard-Baron A, Slama M, Cholley B, Janvier G, Vignon P. Echocardiography in the intensive care unit: from evolution to revolution? Intensive Care Med. 2008;34(2):243–9.

    Article  Google Scholar 

  26. Hall JE. Guyton and Hall textbook of medical physiology. Philadelphia: Elsevier Health Sciences; 2010.

    Google Scholar 

  27. Mozaffarian D, Benjamin EJ, Go AS, Arnett DK, Blaha MJ, Cushman M, Das SR, de Ferranti S, Després J-P, Fullerton HJ. Heart disease and stroke statistics—2016 update. Circulation. 2016;133(4):e38–360.

    Article  Google Scholar 

  28. Bouvier E, Logeart D, Sablayrolles J-L, Feignoux J, Scheublé C, Touche T, Thabut G, Cohen-Solal A. Diagnosis of aortic valvular stenosis by multislice cardiac computed tomography. Eur Heart J. 2006;27(24):3033–8.

    Article  Google Scholar 

  29. Harvey S, Harrison DA, Singer M, Ashcroft J, Jones CM, Elbourne D, Brampton W, Williams D, Young D, Rowan K. Assessment of the clinical effectiveness of pulmonary artery catheters in management of patients in intensive care (PAC-Man): a randomised controlled trial. Lancet. 2005;366(9484):472–7.

    Article  Google Scholar 

  30. Davidson S, Pretty C, Pironet A, Desaive T, Jannsen N, Lambermont B, Morimont P, Chase JG: Estimation of ventricular dead space volume through use of Frank-Starling curves. PLoS ONE. 2017 (In Review).

  31. Kamoi S, Pretty C, Docherty P, Squire D, Revie J, Chiew YS, Desaive T, Shaw GM, Chase JG. Continuous stroke volume estimation from aortic pressure using zero dimensional cardiovascular model: proof of concept study from porcine experiments. PLoS ONE. 2014;9(7):e102476.

    Article  Google Scholar 

  32. Burkhoff D, De Tombe PP, Hunter WC. Impact of ejection on magnitude and time course of ventricular pressure-generating capacity. Am J Physiol Heart Circ Physiol. 1993;265(3):H899–909.

    Google Scholar 

  33. Baan J, Van Der Velde ET. Sensitivity of left ventricular end-systolic pressure–volume relation to type of loading intervention in dogs. Circ Res. 1988;62(6):1247–58.

    Article  Google Scholar 

  34. Davidson S, Pretty C, Kamoi S, Balmer J, Desaive T, Chase JG. Real-time, minimally invasive, beat-to-beat estimation of end-systolic volume using a modified end-systolic pressure-volume relation. In: The 20th world congress of the international federation of automatic control, 9–14 July 2017; 2017.

  35. Lang RM, Badano LP, Mor-Avi V, Afilalo J, Armstrong A, Ernande L, Flachskampf FA, Foster E, Goldstein SA, Kuznetsova T. Recommendations for cardiac chamber quantification by echocardiography in adults: an update from the American Society of Echocardiography and the European Association of Cardiovascular Imaging. J Am Soc Echocardiogr. 2015;28(1):1–39.e14.

    Article  Google Scholar 

  36. Porterfield JE, Kottam AT, Raghavan K, Escobedo D, Jenkins JT, Larson ER, Treviño RJ, Valvano JW, Pearce JA, Feldman MD. Dynamic correction for parallel conductance, GP, and gain factor, α, in invasive murine left ventricular volume measurements. J Appl Physiol. 2009;107(6):1693–703.

    Article  Google Scholar 

  37. Fujiwara H, Ashraf M, Sato S, Millard RW. Transmural cellular damage and blood flow distribution in early ischemia in pig hearts. Circ Res. 1982;51(6):683–93.

    Article  Google Scholar 

  38. Nguyen HB, Rivers EP, Abrahamian FM, Moran GJ, Abraham E, Trzeciak S, Huang DT, Osborn T, Stevens D, Talan DA. Severe sepsis and septic shock: review of the literature and emergency department management guidelines. Ann Emerg Med. 2006;48(1):54. e51.

    Article  Google Scholar 

  39. Jardin F, Farcot J-C, Boisante L, Curien N, Margairaz A, Bourdarias J-P. Influence of positive end-expiratory pressure on left ventricular performance. N Engl J Med. 1981;304(7):387–92.

    Article  Google Scholar 

  40. Vincent J-L, Gerlach H. Fluid resuscitation in severe sepsis and septic shock: an evidence-based review. Crit Care Med. 2004;32(11):S451–4.

    Article  Google Scholar 

  41. Sunagawa K, Maughan W, Friesinger G, Chang M, Sagawa K. Coronary perfusion-pressure and left-ventricular endsystolic pressure–volume relation. In: Circulation. American Heart Association; 1980. p. 203–203.

  42. Grossman W, Braunwald E, Mann T, McLaurin L, Green L. Contractile state of the left ventricle in man as evaluated from end-systolic pressure–volume relations. Circulation. 1977;56(5):845–52.

    Article  Google Scholar 

  43. Tuchschmidt J, Fried J, Astiz M, Rackow E. Elevation of cardiac output and oxygen delivery improves outcome in septic shock. Chest. 1992;102(1):216–20.

    Article  Google Scholar 

  44. Tuttle RR, Mills J. Dobutamine: development of a new catecholamine to selectively increase cardiac contractility. Circ Res. 1975;36(1):185–96.

    Article  Google Scholar 

  45. Weaver ME, Pantely GA, Bristow JD. D LADLEY H: a quantitative study of the anatomy and distribution of coronary arteries in swine in comparison with other animals and man. Cardiovasc Res. 1986;20(12):907–17.

    Article  Google Scholar 

  46. Hasenfuss G. Animal models of human cardiovascular disease, heart failure and hypertrophy. Cardiovasc Res. 1998;39(1):60–76.

    Article  MathSciNet  Google Scholar 

Download references

Authors’ contributions

SD was responsible for the conception and design of the study, analysis and interpretation of data, and drafting the manuscript; CP and JCG were responsible for the conception and design of the study and revising the manuscript; SK and JB were responsible for analysis and interpretation of data and revising the manuscript; AP and TD were responsible for the acquisition of data and revising the manuscript. All authors read and approved the final manuscript.

Acknowledgements

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

Ethics approval and consent to participate

All experimental procedures and protocols used in this investigation were reviewed and approved by the Institutional Animal Care and Use Ethics Committee of the University of Liège, Belgium (Reference Number 14-1726). Their guidelines conform completely with the Guide for the Care and Use of Laboratory Animals published by the US National Institutes of Health (NIH Publication No. 85-23, revised 1996), as well as EU DIRECTIVE 2010/63/EU on the protection of animals used for scientific purposes.

Funding

This work was supported by the Engineering Technology-based Innovation in Medicine (eTIME) consortium grant [eTIME 318943]; the EU FP7 International Research Staff Exchange Scheme (IRSES) grant [#PIRSES-GA-2012-318943]; the Fonds de la Recherche Scientifique (F.R.S) grant; and the University of Canterbury ‘Canterbury Scholarship’ grant. These funding bodies had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Shaun Davidson.

Additional file

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Davidson, S., Pretty, C., Pironet, A. et al. Minimally invasive, patient specific, beat-by-beat estimation of left ventricular time varying elastance. BioMed Eng OnLine 16, 42 (2017). https://doi.org/10.1186/s12938-017-0338-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12938-017-0338-7

Keywords