## Abstract

Right ventricular adaptation to the increased pulmonary arterial load is a key determinant of outcomes in pulmonary hypertension (PH). Pulmonary vascular resistance (PVR) and total arterial compliance (C) quantify resistive and elastic properties of pulmonary arteries that modulate the steady and pulsatile components of pulmonary arterial load, respectively. PVR is commonly calculated as transpulmonary pressure gradient over pulmonary flow and total arterial compliance as stroke volume over pulmonary arterial pulse pressure (SV/PApp). Assuming that there is an inverse, hyperbolic relationship between PVR and C, recent studies have popularised the concept that their product (RC-time of the pulmonary circulation, in seconds) is “constant” in health and diseases. However, emerging evidence suggests that this concept should be challenged, with shortened RC-times documented in post-capillary PH and normotensive subjects. Furthermore, reported RC-times in the literature have consistently demonstrated significant scatter around the mean. In precapillary PH, the true PVR can be overestimated if one uses the standard PVR equation because the zero-flow pressure may be significantly higher than pulmonary arterial wedge pressure. Furthermore, SV/PApp may also overestimate true C. Further studies are needed to clarify some of the inconsistencies of pulmonary RC-time, as this has major implications for our understanding of the arterial load in diseases of the pulmonary circulation.

## Abstract

**Empiric estimates of pulmonary arterial load are prone to errors and may result in overestimation of RC-time** http://ow.ly/QqtAV

## Introduction

In pulmonary hypertension (PH), the increased pulmonary arterial load is a key determinant of right ventricular (RV) dysfunction and failure, which in turn is strongly related to patient survival [1, 2]. Although a complete description of the arterial load can be accomplished by evaluation of the input impedance in the frequency domain [3, 4], this approach is complex and more simple time domain measures of the arterial load are often used, particularly in clinical practice. There are two major haemodynamic components of pulmonary arterial load, a steady component and a pulsatile component, which depend upon different parts and properties of the pulmonary arterial tree [5, 6]. The steady component of load can be represented simplistically by mean pulmonary artery pressure (mPAP), which is influenced by the calibre of small distal resistive arteries as reflected in the pulmonary vascular resistance (PVR). The pulsatile component of load is influenced by the elastic properties of the pulmonary arterial tree, which can be estimated by the total arterial compliance (C) [6, 7]. The product PVR×C characterises the decay of pulmonary arterial pressure in diastole and represents the time-constant of the pulmonary circulation (RC-time, with units of seconds).

Although pulmonary arterial hypertension (PAH) is classically viewed as a disease of the distal resistive vessels, the rise in RV afterload results from an increase in PVR as well as a reduction in C. Thus, accurate assessment of both the steady and pulsatile components of the RV afterload is of critical importance. It has been proposed that a unique feature of the pulmonary circulation is that PVR and total arterial compliance are tightly coupled through an inverse hyperbolic relationship, resulting in a constant RC-time product that prevails in both health and disease [8–21]. A constant RC-time product has important implications for the assessment of RV afterload, namely that PVR and total arterial compliance are redundant measurements and knowledge of one enables the derivation of the other. Exceptions to the constant RC-time product have recently been documented, for example in PH-related to left heart disease [22–25], prompting a re-evaluation of the concept of RC-time constancy [26–29].

The aim of this review is to describe how the steady and pulsatile components of the RV load are quantified, how PVR and total arterial compliance are related in the pulmonary circulation, and the methodological limitations of these measurements in clinical practice.

## Arterial components of right ventricular load

### Overview

The right ventricle ejects blood intermittently into the pulmonary circulation, and this explains why pulmonary arterial pressure is pulsatile. In a hydraulic system where the flow is pulsatile, the relationship between pulsatile flow and pulsatile pressure includes not only the resistance to flow afforded by friction but also the opposition to flow due to both vascular elasticity and the inertia of blood mass [3–7, 30]. Thus, the pulmonary arterial load is influenced by distributed compliant and inertial properties. The pulmonary tree is a system of branching tubes with wave travel and reflection (a transmission line network). The characteristic impedance governs the pressure–flow relationship in the proximal pulmonary arterial until the arrival of the first pressure wave reflection, and is a function of total arterial compliance and blood inertance. Although inertance has a more important role at the pulmonary artery level than at the aortic level, clinical studies most often neglect its contribution to load [7, 31]. It is estimated that the pulsatile component of the pulmonary arterial load accounts for ∼25% of RV power on average, as compared with ∼10% in a systemic arteries–left ventricle coupled system [10, 11].

As summed up by Reuben [32]: “Although the concept of the pulmonary arterial system as a lumped system is an oversimplification of reality, it is believed to be valid to a first approximation (…)”. Indeed, in the vast majority of the clinical studies performed in the field, a simplified model of arterial load based on the two-element Windkessel model is used. In this lumped haemodynamic model, the arterial system is described in terms of two elements: PVR and C. Although the two-element Windkessel model is clearly an oversimplification of reality, it remains invaluable from the clinical viewpoint and is regarded as providing an adequate description of the arterial system, at least during diastole. However, it must be kept in mind that the model has important limitations: mainly that pulse wave velocity, characteristic impedance and wave reflection are not included in this model [7–9, 33]. For these reasons, the two-element Windkessel model does not apply well at high frequency and during the systolic period. In an attempt to improve the description of pressure and flow in systole, the more complex three-element Windkessel model has been used, which takes into account characteristic impedance [7].

As compared with the normal systemic circulation, the normal pulmonary circulation is characterised by its low resistance and high compliance vasculature with great capabilities for distension and recruitment in reserve. Various physiological mechanisms help preserve low pulmonary arterial load and thus protect both RV function and capillary pressure. Increased arterial load has deleterious consequences on RV function and gas exchange. Before discussing the way RC-time is obtained, it is necessary to recall the pathophysiological basis of PVR and total arterial compliance in health and PH states.

### Pulmonary vascular resistance

Physicians are familiar with the concept of PVR, which is a measure of the extent to which the pulmonary circulation “resists” cardiac output (CO). This resistance to mean flow is a function of the blood viscosity and the size and number of vessels [31]. Poiseuille's law predicts that PVR increases in proportion to blood viscosity and inversely with luminal radius to the fourth power. The mean driving pressure of mean pulmonary flow is the transpulmonary pressure gradient, namely the mPAP minus the mean downstream zero-flow pressure (Pzf), which is usually assumed to be equal to pulmonary arterial wedge pressure (PAWP) or left atrial pressure. The driving pressure is related to the frictional resistance of the blood when flow is induced. The standard formula for calculating PVR (based on Ohm's law of electrical circuits) is as follows.1

In precapillary PH, increased PVR occurs as consequence of proliferative remodelling and vasoconstriction of the pulmonary artery in response to chronic pulmonary vascular injury [34]. In pure post-capillary PH, increased mPAP occurs as consequence of increased PAWP while PVR is normal.

The standard PVR equation is integral for the clinical evaluation of pulmonary vascular disease and is a component of the current diagnostic criteria for precapillary PH, which is defined by mPAP ≥25 mmHg, PAWP ≤15 mmHg and PVR >3 Wood units [35]. Whether PAWP represents the true Pzf is fundamental to the physiological validity of the PVR equation. It is not easy to measure Pzf in humans. In health, PAWP appears to be a reasonable estimate of Pzf, which is further supported by the very small gradient between diastolic pulmonary artery pressure (dPAP) and PAWP, since Pzf must be confined between dPAP and PAWP. The same reasoning holds true for venous post-capillary PH where the dPAP–PAWP gradient is small [35, 36]. However, hypoxia and pulmonary vascular diseases leading to precapillary PH may alter the validity of this assumption, with Pzf being significantly above PAWP, as demonstrated when studying either the pressure intercept of the mPAP–CO linear plots when flow is increased or the pulmonary artery pressure (PAP) decay curve after balloon occlusion [37–42]. In cases where Pzf is significantly higher than PAWP, the true PVR will be overestimated if one uses equation 1 [31]. Finally, it must be noted that some studies have focused on total pulmonary resistance, where a zero PAWP is assumed, which also leads to major overestimation of the true PVR.

If one is interested in quantifying the severity of pulmonary vascular disease *per se*, as reflected in the geometric reduction in cross-sectional luminal area (*i.e.* hindrance), adjustment for haematocrit (which is the main determinant of blood viscosity) might be desirable particularly when extremes of haematocrit levels are encountered [37]. Other conceptual limitations of the standard PVR equation include assumptions that blood behaves as a Newtonian fluid and that the arteries behave as rigid tubes while they are both collapsible and distensible. Because normal pulmonary vessels are highly distensible, a more complex curvilinear fit of the mPAP *versus* CO relationship may be used [43], but it must be noted that the multipoint pressure–flow relationship remains linear over physiological ranges of flow under most physiological, pharmacological and pathological conditions [31, 37].

### Compliance

During systole, a fraction of the ejected stroke volume (SV) is stored in the compliant pulmonary arteries, which is equivalent to charging a capacitor in electronics. During diastole, the blood stored in the compliant pulmonary arteries is released. The systolic storage of this reservoir volume dampens pulsatile RV ejection and decreases both RV stroke work (and thus energy expenditure) and mean RV ejection pressure [7, 44, 45]. Total arterial compliance (C) can be viewed as the storage capacity of the entire pulmonary arterial tree and is defined as the increase in blood volume (ΔV) in the arterial system that produces a unit increase in arterial transmural pressure (ΔP):2It can be appreciated that total arterial compliance is difficult to measure *in vivo* based on the above formula. Although ΔP can be readily approximated using the pulmonary arterial pulse pressure (PApp), direct measurement of ΔV is not possible since a fraction of the stroke volume (SV) is discharged out of the arterial system during RV ejection. Despite this obvious limitation, total arterial compliance is often estimated empirically in clinical practice as:3Given that SV exceeds the reservoir volume of the compliant pulmonary arteries, previous studies have documented that the SV/PApp ratio may significantly overestimate total arterial compliance [46, 47]. Indeed, Segers *et al*. [46] first showed that SV/PApp overestimates total arterial compliance by 81% in pigs and 60% in dogs when total arterial compliance was calculated using the reference method, namely the pulse pressure method. Similarly, Muthurangu *et al*. [47] also reported that SV/PApp overestimates total arterial compliance by 61% in patients with suspected PH or congenital heart disease. As a result, authors initially warned against the use of the SV/PApp ratio as an estimate of total arterial compliance [7, 33]. Finally, total arterial compliance must be distinguished from local arterial compliance, which can be estimated by modalities such as magnetic resonance imaging and intravascular ultrasound [48].

Arterial compliance is determined by the prevailing mPAP and by the elasticity of the pulmonary arterial wall, which is mainly determined by the content of elastin (elastic fibres) and collagen (stiffer fibres) in the wall. Arterial compliance decreases when mPAP increases, and this pressure-dependence of pulmonary arterial wall stiffness (1/C) is due to the nonlinear stress–strain relationship [48, 49]. In other words, a decrease in total arterial compliance can result from an increase in mPAP without a true change in the elastic properties of the pulmonary arterial wall. Arterial compliance can also decrease when the pulmonary arterial wall stiffens due to a decreased elastin/collagen ratio as a result of vascular remodelling [48, 49].

## Resistance–compliance relationship of the pulmonary circulation

A common finding among recent studies is an inverse relationship between PVR and total arterial compliance in both health and disease states including PAH, chronic thromboembolic PH and left heart dysfunction [8–29]. This confirms previous findings from the landmark studies in the field [3, 4, 32, 50–55]. The inverse relationship between PVR and total arterial compliance is usually described by a hyperbolic fit, meaning that a rise in PVR is followed by an inverse reduction in total arterial compliance, and *vice versa*. This relationship has now been reported by numerous groups, and is preserved during the course of therapy with PAH-specific agents. Because of this hyperbolic relationship between PVR and total arterial compliance, authors have advocated that the product of PVR and total arterial compliance, which has units of time (referred to as RC-time), remains constant in the pulmonary circulation across the entire range of PAP [8–11]. The RC-time represents the time-constant of the pulmonary circulation, which characterises the decay of pulmonary arterial pressure in diastole. A constant RC-time has important implications for quantification of RV afterload. If RC-time is indeed constant, it suggests that PVR and total arterial compliance are essentially redundant measures and knowledge of one is sufficient to derive the value of the other (see later).

A number of explanations have been put forward for the inverse relationship between PVR and total arterial compliance. It is possible that the prevailing mPAP is itself the main determinant of total arterial compliance. In other words, an increase in PVR will lead to an increase in pressure within the vasculature, which in turn leads to a decrease in compliance due to the nonlinear pressure–diameter relationship of the pulmonary arteries [56]. Alternatively, other groups have hypothesised that the inverse relationship may be explained by the longitudinal distribution of resistance and compliance vessels of the pulmonary circulation [10, 11]. In the pulmonary circulation, the proximal pulmonary arteries account for proportionally less of the total arterial compliance compared with the aorta in the systemic circulation. Thus, it has been suggested that distal pulmonary arteries have not only a resistive function but also a significant capacitive function and contribute to overall compliance.

An often neglected issue is that the accurate interpretation of the precise functional relationship between two calculated variables such as PVR and SV/PApp requires caution. Vascular resistance is a haemodynamic variable that is frequently mathematically coupled, probably because it is a derived variable [57–60]. Furthermore, there may be a strong mathematical coupling related to the fact that PVR and SV/PApp share SV in the numerator of one and the denominator of the other. In other words, the SV *versus* 1/SV relationship is by definition hyperbolic and the relationship between resistance and total arterial compliance is prone to mathematical coupling because SV appears in both the dependent and independent variable. An additional source of mathematical coupling is that mPAP and PApp have an almost 1:1 ratio.

## Prognostic value of PVR and SV/PApp in PH

Both PVR and SV/PApp have been demonstrated to be of prognostic value in PH, although the data supporting PVR is conflicting. The original National Institutes of Health registry, conducted between 1985 and 1988, and the modern era French registry did not find PVR to be an independent predictor of prognosis in PAH [60, 61]. However, the large American REVEAL registry showed that extremely high PVR >32 Wood units was independently associated with worse survival in PAH [62]. PVR may not be a consistent marker of prognosis in PAH since the main determinant of outcomes is RV function, or the ability of the RV to adapt to the increase in afterload. Accordingly, haemodynamic markers of RV function such as CO and right atrial pressure have been consistently demonstrated to be prognostic in PAH [1, 2].

SV/PApp ratio is a strong independent predictor of mortality in patients with PAH [63, 64] and in scleroderma-related PAH [65]. In congestive heart failure patients, SV/PApp provides stronger prediction than PVR (univariate analysis), and SV/PApp remains an independent predictor with multivariate analysis while PVR does not [23]. It has been suggested that SV/PApp is an integrated measure of RV afterload that may be superior to PVR in explaining RV failure and in predicting adverse outcomes in the setting of left-sided heart failure [23, 66]. SV/Papp can also be viewed as lumped parameter that takes into account the contributions of wave reflections and pulsatile load because of the inclusion of PApp in its denominator. When wave reflections return from sites of impedance mismatch to the RV, the RV ejection systolic pressure must increase to exceed this added load resulting in an increase in PApp [6, 26].

## RC-time: definition, calculation and reference values

### Physiological meaning of RC-time and the Windkessel model

RC-time is a time constant that describes the diastolic decay of the pulmonary arterial pressure curve. The concept of RC-time comes from the Windkessel model where the decrease in diastolic pressure is governed by PVR and total arterial compliance. In 1827, Weber first baptised the Windkessel elastic function of aorta. A Windkessel was an air (top) and water (bottom) filled reservoir, similar to that behind the pump of an old-fashioned fire hose, the purpose of which was to transform the rhythmic water output of the plunger strokes into an almost continuous waterjet [7, 44]. It was Frank who popularised the model [45] and viewed the arterial system as a lumped, two-element Windkessel, with resistance and compliance elements disposed in parallel.

After closure of the pulmonary valve, there is no inflow into the pulmonary arterial system. When pressure and flow are in steady-state, the model assumes that the diastolic pulmonary arterial pressure decay is well fitted by a single exponential function of time, P(t), with the time constant (RC, expressed in units of time), according to the following equation [59]:4P_{0} is the pressure at the onset of analysis (t=0, corresponding to the time at or immediately after the pulmonary arterial dicrotic notch pressure) and P_{asymptote} is the theoretical pressure at the infinite time when the exponential decline is completed [59]. Thus, the RC-time represents the time necessary for pressure to decrease in value to 36.8% (1/e) of the P_{0}–P_{asymptote} difference. It will take theoretically three RC-times to complete 95% of this pressure change (fig. 1).

### Calculation of RC-time

The RC-time product is most commonly calculated using an empiric formula where equation 1 is used to calculate PVR and where total arterial compliance is estimated from SV/PApp (equation 3). Doing this, the basic principles underlying the Windkessel function, namely the exponential decay of pulmonary arterial pressure in diastole (equation 4) is neglected, and PVR and total arterial compliance values are fixed *a priori*. It is assumed that both PVR (equation 1) and SV/PApp are reliable estimates of true PVR and total arterial compliance, respectively, and that their product accurately reflects the RC-time constant of the Windkessel model. From this empiric method, several groups have documented RC-time values [8–29]. The main studies, corresponding study populations and empirical estimates of RC-time are detailed in table 1.

Alternative methods of calculating RC-time have been described. To the best of our knowledge, the pulmonary arterial RC-time was first calculated in 1959 by Engelberg and DuBois [50] using the semi-logarithmic method in rabbit lungs. In 1971, Reuben [32] used the same method to calculate RC-time in humans with various forms of cardiac and pulmonary diseases. The semi-logarithmic method consists of a log-transformation of equation 4 and the RC-time is calculated as follows:5

The PAWP or the left ventricular (LV) filling pressure (when available) is used as P_{asymptote}, and this also allows calculation of PVR. Then total arterial compliance is calculated as the RC-time/PVR ratio. The main studies, corresponding study populations and semi-logarithmic estimates of RC-time are detailed in table 2.

### Conflicting results on RC-time value in health and diseases

Recent studies have proposed the notion of constant RC-time (empiric estimate) of ∼0.5–0.7 s in both health and disease [8–11]. Subsequently, it has emerged that various conditions are associated with significant alterations of RC-time. Recently, Tedford *et al*. [22] documented shortened empiric RC-times in post-capillary PH, thus suggesting an increased pulsatile afterload when compared to precapillary PH at a given level of PVR. This has been subsequently confirmed by other authors [23–25]. These studies still documented an inverse hyperbolic relationship between PVR and SV/Papp, and supported the concept that the mean PVR×C product was in the 0.25–0.53 s range in post-capillary PH [22–25]. Very recently, Tedford [26] has further challenged the paradigm of constant RC-times in normotensive subjects *versus* PAH patients. Indeed, a shortened pulmonary RC-time was documented in 415 subjects with mPAP ≤25 mmHg (mean RC-time=0.36 s) as compared with 593 subjects with PAH (RC-time=0.53 s) [26].

Using the semi-logarithmic method (table 2), the recovered RC-time averaged 0.16–0.39 s in normal subjects but was higher in those with PH, averaging 0.21–0.83 s (table 2) [32, 51–55]. The shorter RC-times derived from the semi-logarithmic method would be more consistent with more complete expression of the Windkessel function compared with RC-times calculated from empiric estimates. In some studies performed in patients with PH, the RC-times calculated by the semi-logarithmic method are actually similarly as long as RC-times obtained with the empirical method.

## RC-time: a critical evaluation and new proposal

### Is the RC-time truly constant?

While it is clear that RC-times can no longer be considered constant when precapillary PH, post-capillary PH and controls are compared, constancy of RC-time has been suggested within a given disease subgroup. Unfortunately, a critical analysis of the data may also challenge this notion. Constancy implies a narrow range of RC-time values. This is not supported by most of the data, which show a large sd and thus a large coefficient of variation (100×sd/mean). In numerous studies (table 1), the RC-time coefficient of variation is >30%, which reflects great variation of RC-time within the sample. One may simply consider the meaning of a 30% coefficient of variation of say resting heart rate in a given population with a mean of 80 beats·min^{−1}, assuming a normal distribution. The 95% confidence interval of 33–127 beats·min^{−1} certainly points to major variation of resting heart rate in the sample rather than “constancy”, with the upper limit of the 95% confidence interval being four-times greater than the lower limit. Our group has recently documented major scattering of RC-time in the precapillary PH subgroup (fig. 2) [27]. In post-capillary PH, RC-time decreases linearly with increasing left atrial pressure [22–24]. Finally, in precapillary PH, it has been suggested that at high values of PVR, compliance reaches a minimum and that further increases in PVR would result in an increase in RC-time [17]. The above observations have led to concerns and cast doubt regarding the RC-time constancy paradigm even among specific disease subgroups.

Despite the irrefutable inverse relationship between PVR and SV/PApp, we suggest that there remains significant scatter around the hyperbolic curve fits and RC-time [26]. The accurate measurement of PApp using fluid-filled catheters is especially challenging and a lot of this scatter may also be due to measurement error [26].

### Are estimates of RC-time too long?

As previously discussed, RC-time values in the range of 0.50–0.70 s have been widely documented in numerous recent studies performed in precapillary PH (table 1). Some studies using the semi-logarithmic method have also documented a similarly long RC-time in PH (table 2). Assuming a heart period of 0.80 s (RV ejection time ∼0.30 s, diastolic time ∼0.50 s) and a RC-time of 0.60 s this implies: 1) that <63.2% (*i.e.* 100–36.8%) of the Windkessel function has been manifested at end-diastole; and 2) that completion of the Windkessel function may theoretically only be observed after 1.8 s. This supports incomplete expression of the Windkessel function in PH, with RV ejection occurring at a time when pulmonary arterial emptying is not completed. Because the accurate measurement of PApp using fluid-filled catheters is especially challenging, one possibility is that RC-time may have been overestimated due to measurement error [26]. This is unlikely to be the main explanation as short RC-time has been documented in non-PH subjects and in post-capillary PH. The non-normal distribution of PVR and SV/PApp could also contribute to an abnormally high and unreliable value of the product of their mean values [27]. Initial studies performed in precapillary PH have used total pulmonary resistance (mPAP/CO) instead of PVR and this may have also contributed to overestimation of true PVR and thus, to abnormally long RC-time reference values.

In keeping with the hypothesis of a useful Windkessel function in PH, and in keeping with the widely admitted notion that pulmonary flow is near-zero at end-diastolic pressure [67, 68], here we suggest that the RC-time may have been overestimated, mainly due to the marked overestimation of resistance by using PVR, and, to a lesser extent, due to overestimation of total arterial compliance using SV/PApp. Indeed, if Pzf pressure is higher than PAWP in the setting of pulmonary vascular disease, especially precapillary PH, the true PVR will be overestimated if one uses equation 1.

In the setting of PAH, increases in resting dPAP are responsible for increased mPAP [67, 68], and are related to pulmonary rarefaction, vasoconstriction, vascular remodelling and intravascular plexiform or thrombotic lesions at the level of small pulmonary arteries [34]. The loss of pulmonary vessels, the major narrowing of open vessels, and the heightened sensitivity to hypoxia and other vasoconstrictive agents [34] may well underlie an increase in the critical closing pressure at which the flow dramatically decreases and stops (*i.e.* Pzf). Due to the fact that pulmonary vessels are collapsible [69], the high Pzf exceeds PAWP and Pzf becomes the effective outflow pressure.

Similarly, SV/PApp may also overestimate total arterial compliance when SV exceeds the Windkessel reservoir volume [44, 46, 47]. Finally, some studies have used total pulmonary resistance as an estimation of PVR, implying that Pzf is 0 mmHg. Clearly, Pzf cannot be lower than PAWP (or left atrial pressure) and we concur with other authors [26, 28, 70] who have advised against the use of a 0 mmHg Pzf, which leads to unrealistically high RC-time values.

### Improving the estimation of RC-time

Because pulmonary flow is near-zero at end-diastole [67, 68], near complete expression of the Windkessel function is thought to occur during diastole, and it is anticipated that the order of magnitude of RC-time should be much shorter than diastolic time (perhaps approximately one third of diastolic time since three RC-times correspond to 95% of the pressure fall to asymptotic pressure) [71]. In this setting, RC-time will be constrained by diastolic time if near complete expression of the Windkessel function occurs at varying heart rates. At the systemic circulation level, a meta-analysis involving 36 major studies on Windkessel parameters has documented that the diastolic time/systemic RC-time ratio was independent of body mass in mammals, and this may be viewed as a contributory factor in the left ventricular–arterial coupling [72]. In the right-sided counterpart, it seems logical that diastolic time (and thus heart rate) and the mechanical properties of the pulmonary arterial tree (PVR and total arterial compliance) may also be mutually constrained. The right ventricle and pulmonary arterial tree form a strongly coupled unit. Recently, an unusually strong linear relationship between heart rate and mPAP has been described in PH patients performing exercise [73]. In this respect, a fixed diastolic time/RC-time ratio may be especially important on exercise, as delayed Windkessel function may compromise RV ejection if pulmonary arteries emptying is not completed.

As recalled by Basnet *et al*. [53]: “The time constant can be calculated from the pressure curve between any two points during diastole”. Indeed, the fundamental property of a monoexponential function is that, at each time, the rate of change of a given variable (here pressure) is proportional to the value of the variable itself at this time. In other words, in cases where decay of diastolic pressure is truly a monoexponential function of time, there must be a significant negative linear relationship between the instantaneous rate of pressure decline and pressure itself. The slope of the relationship (−1/Tau) can be calculated using the so-called derivative method popularised by Raff and Glantz [74] to help calculate the LV relaxation time constant. Overall, numerous, mathematical tools that take the pressure asymptote into account are available to reliably calculate the time constant (Tau) of exponential decay processes [59, 74, 75]. Cardiac physiologists are well aware of such tools allowing the accurate calculation the LV relaxation time constant, as has been documented in many previously published studies. To our knowledge, these mathematical tools remain to be used to study pulmonary arterial RC-time. These mathematical tools will also allow the validity of the exponential model applied to the PAP decay in diastole to be checked [76], and thus the validity of the RC-time model in both health and disease. Because pulmonary capillary pressure/flow remain pulsatile in the normal pulmonary circulation, the classic Windkessel concept of fully dampening pulsatility and providing steady pressure/flow at the capillary level may not be fully applicable to healthy subjects.

## Conclusion

In PH, the increased pulmonary arterial load is a key determinant of RV dysfunction and thus prognosis. PVR and total arterial compliance modulate the steady and pulsatile component of pulmonary arterial load, respectively. In the pulmonary circulation, PVR is inversely related to total arterial compliance, such that an increase or decrease in PVR is followed by a respective decrease or increase in total arterial compliance. Assuming that there is a hyperbolic relationship between PVR and total arterial compliance, recent studies have popularised the concept that the RC-time of the pulmonary circulation is “constant” in health and diseases. However, emerging evidence suggests that this concept should be challenged, and different pathophysiological states are clearly associated with significant alteration in RC-times, such as the shortened RC-times documented in post-capillary PH and in normotensive subjects. Furthermore, reported RC-times in the literature have consistently demonstrated significant scatter around the mean.

Another potential problem with pulmonary RC-time is the methodology that is employed to calculate it. Empiric estimates of PVR and total arterial compliance may be prone to errors. In precapillary PH, the true PVR can be overestimated if one uses the standard PVR equation because Pzf may be significantly higher than PAWP. Furthermore, SV/PApp may also be associated with a potential overestimation of true total arterial compliance. With these considerations, the recently published empiric RC-times might be overestimates and are too long to be consistent with a useful Windkessel function. The inaccuracies associated with empiric PVR and total arterial compliance estimates may well change opinions on the distribution of compliance and resistance in the lung, and this point deserves further study.

A re-evaluation of the current concept of RC-time is needed, using alternative methodologies that have been employed to study the systemic circulation, and high-fidelity pressure instruments to improve measurement precision. Further studies are needed to clarify some of the inconsistencies of pulmonary RC-time, as this has major implications for our understanding of the arterial load in diseases of the pulmonary circulation.

## Footnotes

Conflict of interest: None declared.

- Received May 12, 2015.
- Accepted July 8, 2015.

- Copyright ©ERS 2015