Tutorials › Quantitative Systems Pharmacology (QSP) › Integrated Case: Congestion and Diuretic Response in Heart Failure

Cardiovascular and Renal Physiology · Tutorial 388 of 799

Integrated Case: Congestion and Diuretic Response in Heart Failure

Follow a simplified heart-failure model from diuretic-induced sodium loss through volume and pressure changes to renal filtration, while checking what the predicted trajectory can and cannot establish.

Advanced 12 min read

What You'll Learn

  • Define a short-horizon sodium-balance model for diuretic-induced natriuresis in heart failure.
  • Convert cumulative sodium change into an estimated extracellular-volume change.
  • Calculate how arterial and venous pressures jointly determine the renal perfusion gradient.
  • Trace predicted MAP, venous pressure, and GFR over time using a hypothetical case.
  • Compare a stronger natriuretic perturbation with the base case and assess its limitations.
  • Identify sign, unit, and extrapolation errors that can distort cardiorenal predictions.

From Natriuresis to a Cardiorenal Trajectory

In heart failure, the treatment response cannot be inferred from urine sodium loss alone. Removing sodium tends to reduce extracellular volume and congestion, but the resulting change in arterial pressure may reduce renal perfusion. At the same time, lower venous pressure can relieve a back-pressure on the kidney. The net prediction depends on how these effects combine over time.

This tutorial follows one deliberately simplified, hypothetical case from cumulative sodium balance to volume, mean arterial pressure (MAP), renal venous pressure, the renal perfusion gradient, and glomerular filtration rate (GFR). The aim is to make the chain of assumptions and signs explicit. It is not a clinical dosing model or a substitute for the more complete feedback mechanisms represented in cardiorenal physiology.

Definition: Let \(x(t)\) be the cumulative change in body sodium, in mmol, relative to the start of a short simulation. Let \(\Delta V(t)\) be the corresponding change in effective extracellular volume, in L. In this case, \(x<0\) denotes net sodium loss and \(\Delta V<0\) denotes a fall in effective extracellular volume. MAP and renal venous pressure \(P_v\) are in mmHg; the renal perfusion gradient is \(G=\mathrm{MAP}-P_v\), also in mmHg.

Case Model and Assumptions

The earlier material on whole-body sodium and water balance establishes the role of net balance in changing body stores; the earlier material on heart-failure cardiorenal feedback establishes why arterial and venous pressures both matter. Here, a reduced model translates those ideas into a numerical trajectory. Its assumptions are intentionally restrictive: the sodium perturbation is represented by one state, sodium concentration for the volume conversion is held fixed, and pressure and GFR responses are approximated by linear relations around a heart-failure baseline.

For the hypothetical case, the net sodium-balance response after the intervention is represented by a first-order approach to a negative limiting change. The rate constant is \(k=0.0200\ \mathrm{day}^{-1}\), and the limiting cumulative change is \(-2939\ \mathrm{mmol}\). The initial change is zero. Thus the model is

$$ \frac{dx}{dt}=-k\bigl(x+2939\bigr),\qquad x(0)=0. $$

This forcing summarizes the net balance response in the example; it does not represent a detailed drug concentration or nephron model. The remaining relations are

$$ \Delta V(t)=\frac{x(t)}{140\ \mathrm{mmol/L}}, \qquad \mathrm{MAP}=100+2.00\,\Delta V, \qquad P_v=14+2.365\,\Delta V, $$

where the pressure slopes have units of mmHg/L. The baseline has \(\Delta V=0\), MAP \(=100\ \mathrm{mmHg}\), and \(P_v=14\ \mathrm{mmHg}\), so \(G=86\ \mathrm{mmHg}\). GFR is approximated locally as

$$ \mathrm{GFR}=76.0+0.500\,(G-86)\quad \mathrm{mL/min}. $$

The coefficient \(0.500\) has units of \((\mathrm{mL/min})/\mathrm{mmHg}\). These pressure and GFR equations are stipulated local approximations for this case, not universal physiological laws. They assume that the trajectory remains near enough to the baseline for linear relations to be informative, and they omit changes in cardiac output, vascular resistance, autoregulation, hormones, and kidney injury.

Worked Example: Establishing the Baseline and Reading the Signs

At the start of the hypothetical simulation, \(x=0\) and \(\Delta V=0\). Calculate the baseline renal perfusion gradient and GFR, and state what a negative sodium change implies in this model.

Gradient. Subtract renal venous pressure from MAP:

$$ G=100-14=86\ \mathrm{mmHg}. $$

GFR. Substituting this gradient into the local relation gives

$$ \mathrm{GFR}=76.0+0.500(86-86)=76.0\ \mathrm{mL/min}. $$

Sign interpretation. If net sodium balance becomes negative, then \(x<0\), so \(\Delta V=x/140<0\). The stipulated equations then predict both MAP and \(P_v\) will fall. Whether \(G\) rises or falls depends on which pressure falls more; it cannot be decided from the sodium-loss sign alone.

Deriving and Evaluating the Sodium Trajectory

The sodium equation can be solved by shifting the state by its limiting value. Define \(y(t)=x(t)+2939\). Then \(dy/dt=dx/dt=-ky\), with \(y(0)=2939\). The solution to this first-order decay equation is \(y(t)=2939e^{-kt}\). Subtracting 2939 gives the cumulative change:

$$ x(t)=2939e^{-kt}-2939 =-2939\bigl(1-e^{-kt}\bigr)\ \mathrm{mmol}. $$

This derivation shows why the trajectory starts at zero and becomes increasingly negative for \(t>0\), while approaching the model's negative limiting change. Dividing by \(140\ \mathrm{mmol/L}\) converts a sodium amount into an estimated volume change in litres. The conversion assumes a fixed effective sodium concentration; it is a simplifying balance approximation, not a claim that all body water is in one compartment.

Worked Example: Three-Day Congestion and Pressure Trajectory

Calculate the sodium change, volume change, and hemodynamic predictions at day 3. Use the unrounded volume change for the downstream calculations, then report results to three significant figures.

First evaluate the exponential:

$$ e^{-0.0200(3)}=e^{-0.0600}\approx0.941765. $$

The sodium change is therefore

$$ x(3)=-2939(1-0.941765) \approx-171.15\ \mathrm{mmol}. $$

The volume change is

$$ \Delta V(3)=\frac{-171.15\ \mathrm{mmol}}{140\ \mathrm{mmol/L}} \approx-1.2225\ \mathrm{L}\approx-1.223\ \mathrm{L}. $$

Now use the unrounded value, approximately \(-1.222529\ \mathrm{L}\), in each pressure equation:

$$ \mathrm{MAP}(3)=100+2.00(-1.222529) \approx97.55\ \mathrm{mmHg}, $$
$$ P_v(3)=14+2.365(-1.222529) \approx11.11\ \mathrm{mmHg}. $$

The renal perfusion gradient is \(97.5549-11.1087\approx86.4462\ \mathrm{mmHg}\), or \(86.45\ \mathrm{mmHg}\). Thus

$$ \mathrm{GFR}(3)=76.0+0.500(86.4462-86) \approx76.22\ \mathrm{mL/min}. $$

In this specified case, venous pressure falls more than MAP, so the pressure gradient increases slightly even though MAP falls. The model consequently predicts a small GFR increase relative to its baseline value of \(76.0\ \mathrm{mL/min}\). This is a model-based local prediction, not evidence that diuresis generally improves GFR.

Comparing a Stronger Natriuretic Perturbation

A useful model comparison changes one assumption at a time. Keep the limiting sodium change, pressure relations, and GFR relation fixed, but let the rate constant increase from \(0.0200\) to \(0.0300\ \mathrm{day}^{-1}\). In this parameterization, that change represents a faster approach toward the same limiting sodium deficit. It does not assert that a particular dose produces a particular rate constant.

Worked Example: Faster Natriuresis at Day 3

For \(k=0.0300\ \mathrm{day}^{-1}\), calculate the day-3 predictions and compare them with the base case.

The sodium change is

$$ x(3)=-2939\bigl(1-e^{-0.0300(3)}\bigr) =-2939(1-e^{-0.0900}) \approx-253.0\ \mathrm{mmol}. $$

Therefore \(\Delta V(3)=-2939(1-e^{-0.0900})/140\approx-1.806830\ \mathrm{L}\approx-1.807\ \mathrm{L}\). Substitution gives MAP \(\approx100+2(-1.806830)=96.39\ \mathrm{mmHg}\), and \(P_v\approx14+2.365(-1.806830)=9.73\ \mathrm{mmHg}\). The gradient is approximately \(96.39-9.73=86.66\ \mathrm{mmHg}\), giving

$$ \mathrm{GFR}\approx76.0+0.500(86.6595-86) \approx76.33\ \mathrm{mL/min}. $$

Compared with the base case at day 3, the faster response predicts a larger fall in volume, MAP, and venous pressure. Because the stipulated venous-pressure response is steeper than the MAP response, the gradient and modelled GFR are slightly higher. This comparison illustrates competing influences, not a general recommendation for faster or more intensive diuresis.

Checking the Trajectory and Choosing an Approach

A trajectory should be checked at each link in the calculation, rather than judged from its final GFR alone. In this case, \(x(t)\) must start at zero and become negative; division by a positive sodium concentration must give a negative volume change; the pressure equations must then be evaluated with that signed volume change; and \(G\) must equal MAP minus \(P_v\). A useful diagnostic is to compare all predicted quantities with their baseline values and check that units are consistent. For example, sodium in mmol divided by mmol/L gives L, not mL.

The model's algebraic pressure and GFR relations are a reduced alternative to a fully coupled dynamic model. They are convenient when the purpose is to trace the direction and approximate short-term size of effects under fixed assumptions. A more mechanistic approach would represent cardiac output, vascular resistance, renal blood flow, filtration, and compensatory feedback as interacting states or functions. Such a model can express delayed responses and nonlinearities, but requires more data and parameters and can be harder to identify. Neither level of detail automatically makes a prediction reliable: assumptions and calibration still matter.

Several failure modes deserve attention:

  • Reversing the venous-pressure sign. The gradient is MAP minus \(P_v\), not MAP plus \(P_v\). A fall in \(P_v\) raises the gradient if MAP is held fixed.
  • Using rounded intermediate values. Rounding \(\Delta V\) before calculating pressure can shift downstream results. Retain the unrounded trajectory for calculations, then round reported outputs.
  • Extrapolating a local linear model. The volume equation approaches a deficit of \(2939/140\approx21.0\ \mathrm{L}\), an extreme value that is not a credible long-term clinical prediction. It signals that the forcing and linear approximations are intended only for a short illustrative horizon.
  • Equating modelled GFR with observed clinical response. The equation omits autoregulation and other responses. A predicted pressure-gradient change does not by itself establish the measured change in filtration or patient benefit.

The earlier material on pressure-natriuresis and compensatory sodium retention emphasizes that renal output can respond to pressure and that sodium balance is subject to feedback. Here, those processes are compressed into a prescribed net sodium trajectory. If compensation changes substantially during the simulated interval, the assumed \(x(t)\) may be wrong, and the pressure predictions will inherit that error. This is an important limitation, not a minor numerical detail.

Key takeaway: In this case, diuretic-induced sodium loss lowers effective volume, MAP, and venous pressure. The renal perfusion gradient depends on the difference between MAP and venous pressure, so its direction—and the corresponding local GFR prediction—depends on their relative changes. Preserve signs, units, and unrounded intermediate values, and treat short-horizon linear predictions as conditional on their assumptions.

Check Your Understanding

Answer each question using the definitions and equations in this case. Solutions are included below the questions.

  1. If \(x=-140\ \mathrm{mmol}\), what is \(\Delta V\) under the stated conversion?
  2. If MAP is \(98\ \mathrm{mmHg}\) and \(P_v\) is \(12\ \mathrm{mmHg}\), what is the renal perfusion gradient?
  3. Why can a fall in MAP occur at the same time as a rise in the renal perfusion gradient?
  4. What does increasing \(k\) change in this model when the limiting sodium change is held fixed?
  5. Give one reason the day-3 GFR prediction should not be treated as a guaranteed clinical outcome.

Solutions. (1) \(\Delta V=-140/140=-1.00\ \mathrm{L}\). (2) \(G=98-12=86\ \mathrm{mmHg}\). (3) The gradient is a difference: if \(P_v\) falls by more than MAP, then \(\mathrm{MAP}-P_v\) can rise even as MAP itself falls. (4) A larger \(k\) makes \(x(t)\) approach its limiting negative change faster; it does not change the specified limiting value. (5) For example, the GFR equation is a local linear approximation that omits autoregulation and other physiological feedback, so the model prediction is conditional rather than guaranteed.

Clinical Biostats