Tutorials › Quantitative Systems Pharmacology (QSP) › Natriuretic Peptides as Cardiorenal Volume Signals

Cardiovascular and Renal Physiology · Tutorial 372 of 799

Natriuretic Peptides as Cardiorenal Volume Signals

Learn how a simplified cardiorenal model links atrial and ventricular stretch to natriuretic-peptide signaling and predicts its renal and vascular consequences.

Advanced 12 min read

What You'll Learn

  • Distinguish atrial natriuretic peptide and B-type natriuretic peptide release signals.
  • Represent stretch-dependent release and receptor activation with a compact model.
  • Calculate predicted changes in sodium excretion from peptide signaling.
  • Estimate a direct vascular contribution to mean arterial pressure.
  • Identify assumptions, limitations, and diagnostic checks for these predictions.

From Cardiac Stretch to a Cardiorenal Signal

A rise in cardiac filling can stretch the atria and ventricles. The heart responds in part by releasing natriuretic peptides: atrial natriuretic peptide (ANP), released predominantly by atrial cardiomyocytes, and B-type natriuretic peptide (BNP), released predominantly by ventricular cardiomyocytes in response to increased wall stress. Their effects oppose sodium and water retention: they can increase renal sodium excretion, inhibit sodium-retaining hormonal signals, and relax vascular smooth muscle.

For a quantitative systems model, the useful question is not whether a peptide is simply “high” or “low,” but how stretch changes its signal, how that signal changes renal and vascular outputs, and under what assumptions those outputs can be interpreted as changes in volume and pressure. The equations below are deliberately reduced models. They represent direction and approximate magnitude, not every physiological pathway.

Definition: A natriuretic-peptide signal is a model quantity representing the effective activation of peptide-responsive pathways, principally through natriuretic peptide receptor A (NPR-A) and its cyclic GMP signal. ANP and BNP can both contribute to this activation. A stretch input is a dimensionless index of cardiac wall stretch relative to a chosen reference condition.

A Compact Release and Signaling Model

Let \(s_A\) and \(s_V\) denote normalized atrial and ventricular stretch indices, respectively, and let \(a\) and \(b\) denote normalized circulating ANP and BNP concentrations. A simple bounded release response is \(f(s)=s/(1+s)\), for \(s\geq 0\). One possible first-order model is

$$ \tau_A \frac{da}{dt}=a_0+a_1 f(s_A)-a,\qquad \tau_B \frac{db}{dt}=b_0+b_1 f(s_V)-b. $$

Here, \(a_0\) and \(b_0\) are baseline normalized concentrations, \(a_1\) and \(b_1\) are release-response scales, and \(\tau_A\) and \(\tau_B\) are positive time constants. At a fixed stretch, the steady-state concentrations are \(a=a_0+a_1f(s_A)\) and \(b=b_0+b_1f(s_V)\). The first-order form represents production approaching a stretch-dependent target; it does not claim that secretion or clearance is literally a single process.

If \(a\) and \(b\) are scaled by their respective affinities for the shared receptor pathway, define \(x\) as their combined dimensionless effective concentration. Under a simple equilibrium binding assumption, the fraction of available receptor activation is

$$ \theta=\frac{x}{1+x}. $$

This expression assumes that the relevant ligand contributions can be summarized by \(x\), receptor binding equilibrates quickly relative to the modeled response, and the maximum receptor activation is normalized to one. It is a useful saturation function, not a detailed receptor-trafficking model. In particular, it does not represent all differences between ANP and BNP processing, clearance, or tissue access.

Worked Example: Translating Two Stretch Inputs Into a Receptor Signal

Suppose a hypothetical model uses \(f(s)=s/(1+s)\), \(a_0=0.20\), \(a_1=0.80\), \(b_0=0.10\), and \(b_1=1.00\). For this illustration, take normalized receptor-affinity-scaled concentrations as \(x=a+b\). Compare a condition with \(s_A=0.50\) and \(s_V=0.25\) with one having \(s_A=1.50\) and \(s_V=1.00\).

At the first condition, \(f(0.50)=0.50/1.50=0.3333\) and \(f(0.25)=0.25/1.25=0.2000\). Thus \(a=0.20+0.80(0.3333)=0.4667\) and \(b=0.10+1.00(0.2000)=0.3000\). The combined signal is \(x=0.7667\), giving \(\theta=0.7667/(1+0.7667)=0.4340\).

At the higher stretch condition, \(f(1.50)=1.50/2.50=0.6000\) and \(f(1.00)=1.00/2.00=0.5000\). Therefore \(a=0.20+0.80(0.6000)=0.6800\), \(b=0.10+1.00(0.5000)=0.6000\), and \(x=1.280\). The receptor activation is \(\theta=1.280/2.280=0.5614\). The model therefore predicts greater effective peptide-pathway activation at higher stretch. These normalized values are hypothetical and are not clinical concentration estimates.

Renal Actions and Sodium Excretion

ANP and BNP can promote natriuresis through several linked processes. Their receptor-mediated cyclic GMP signaling can alter renal vascular tone and filtration, reduce tubular sodium reabsorption, and oppose renin and aldosterone signaling. The relative contribution of each route depends on physiological state and model structure. A reduced output model can represent their combined effect without claiming that all of the effect passes through one nephron segment.

For example, let \(E_{\mathrm{Na}}\) be sodium excretion in mmol/day and use the phenomenological relationship

$$ E_{\mathrm{Na}}=E_0+E_{\max}\theta, $$

where \(E_0\) is excretion at zero modeled peptide activation under the specified background conditions, and \(E_{\max}\) is the additional excretion at full activation. This relationship assumes that other determinants—such as sodium intake, renal function, pressure-natriuresis, and other hormones—are held fixed or absorbed into the fitted parameters. It is a local, illustrative response function, not a universal excretion law.

Derivation: For \(x\geq 0\), the activation function is increasing because
$$ \frac{d\theta}{dx}=\frac{1}{(1+x)^2}>0. $$

If \(E_{\max}>0\), then \(dE_{\mathrm{Na}}/dx=E_{\max}/(1+x)^2>0\). Thus, within this model and with the background held fixed, increasing effective peptide signal increases predicted sodium excretion. The saturation is also explicit: as \(x\) becomes large, \(\theta\) approaches one and the modeled excretion approaches \(E_0+E_{\max}\), rather than increasing without bound.

Worked Example: Predicting Natriuresis From a Higher Peptide Signal

Consider a hypothetical fixed-background model with \(E_0=80\) mmol/day and \(E_{\max}=40\) mmol/day. Compare \(x=0.300\) with \(x=1.625\).

At \(x=0.300\), \(\theta=0.300/(1+0.300)=0.230769\). Therefore,

$$ E_{\mathrm{Na}}=80+40\left(\frac{0.300}{1.300}\right) =80+9.23077=89.23077\ \text{mmol/day}. $$

At \(x=1.625\), \(\theta=1.625/(1+1.625)=0.619048\), so

$$ E_{\mathrm{Na}}=80+40\left(\frac{1.625}{2.625}\right) =80+24.76190=104.76190\ \text{mmol/day}. $$

Using the unrounded predictions, the difference is \(40[(1.625/2.625)-(0.300/1.300)]=15.5311\) mmol/day, or 15.53 mmol/day to four significant figures. This is the predicted increase in sodium excretion under the stated fixed-background assumption. It is not by itself a prediction of the eventual volume change: that also depends on sodium intake, water balance, and the responses of other systems.

Vascular Effects and Pressure Predictions

Natriuretic-peptide signaling can relax vascular smooth muscle and lower systemic vascular resistance. To isolate this direct route, suppose cardiac output (CO) and central venous pressure (CVP) are fixed over the comparison. In a lumped circulation, mean arterial pressure (MAP) is approximated by \(MAP=CVP+CO\,R_{\mathrm{sys}}\), where \(R_{\mathrm{sys}}\) is systemic vascular resistance. A simple signal-dependent resistance is

$$ R_{\mathrm{sys}}=R_0(1-\gamma\theta), $$

where \(R_0\) is resistance without the modeled peptide effect and \(0\leq\gamma\leq 1\) is the maximum fractional reduction in this simplified formulation. The model requires \(1-\gamma\theta\geq 0\). It omits changes in cardiac output, reflex responses, and the renal contribution to pressure.

Worked Example: Isolating the Direct Vascular Pressure Effect

Use a hypothetical \(R_0=18\) mmHg·min/L, \(\gamma=0.25\), \(CO=5.0\) L/min, and \(CVP=5.0\) mmHg. Compare the receptor activation values displayed to three decimal places, \(\theta=0.231\) and \(\theta=0.619\), and carry those values through the calculation.

At the lower activation, \(R_{\mathrm{sys}}=18[1-0.25(0.231)]=16.9605\) mmHg·min/L. Hence \(MAP=5.0+(5.0)(16.9605)=89.8025\) mmHg, or 89.8 mmHg.

At the higher activation, \(R_{\mathrm{sys}}=18[1-0.25(0.619)]=15.2145\) mmHg. Hence \(MAP=5.0+(5.0)(15.2145)=81.0725\) mmHg, or 81.1 mmHg. The difference calculated from these displayed activation values is \(89.8025-81.0725=8.7300\) mmHg: a predicted direct-route decrease of 8.73 mmHg. CO and CVP were held constant, so this is not a prediction of the total pressure change after all cardiovascular and renal feedbacks.

Interpreting Changes in Volume and Pressure

An increase in peptide signaling can affect pressure by at least two routes in this reduced picture. First, vascular relaxation can lower resistance directly. Second, increased renal sodium excretion can contribute to a decline in body sodium and extracellular volume over time if excretion exceeds intake. A lower effective circulating volume may then influence venous return, cardiac output, and pressure. These routes operate on different time scales and can be modified by compensatory responses.

Accordingly, the sign of a short-term direct vascular effect is easier to infer from the resistance equation than the size or timing of a long-term MAP change. Likewise, greater sodium excretion during a time interval does not alone establish a final volume or pressure. The interpretation requires comparing excretion with intake and considering how filtration, tubular transport, hormonal signals, and hemodynamics adapt.

1
Specify stretch.
Define atrial and ventricular stretch relative to a reference, and state whether they are measured inputs or predicted states.
2
Calculate peptide signal.
Use the stated release and clearance representation to obtain normalized concentrations, then combine them into the effective signal.
3
Predict outputs conditionally.
Calculate renal and vascular changes with the other modeled inputs held fixed; label those predictions as pathway-specific, not whole-system outcomes.
4
Check the time scale and feedbacks.
Distinguish an immediate vascular response from accumulated sodium and volume changes, and assess whether omitted feedbacks could alter the prediction.

Assumptions, Diagnostics, and Alternatives

The model is most defensible when stretch lies within the range used to define its response, the normalized peptide signal is nonnegative, the receptor activation function is appropriate for the intended time scale, and the renal and vascular parameters are calibrated for the background condition being studied. Fixed-background calculations are useful for isolating a pathway, but they should not be mistaken for a closed-loop prediction. Changes in kidney function, receptor responsiveness, or other hormonal inputs can change the response.

Several checks help identify a faulty or overextended model. Confirm that higher stretch produces a nondecreasing release target, that \(\theta\) remains between zero and one, that increasing \(\theta\) increases sodium excretion when \(E_{\max}>0\), and that the resistance expression remains nonnegative. Simulated concentrations or outputs that diverge despite fixed inputs may indicate a sign error, an incorrectly specified clearance term, or an unstable implementation. Compare predictions with plausible ranges for the intended application; a mathematically bounded response is not automatically physiologically calibrated.

The receptor-occupancy model is a compact alternative to separate ANP and BNP receptor-binding, processing, and clearance equations. The direct renal response is a compact alternative to explicitly representing filtration, segmental tubular transport, and hormonal suppression. These reduced forms require fewer parameters and make pathway direction transparent, but they cannot identify which mechanism generated an observed change. More detailed models are preferable when data distinguish the mechanisms or when the question concerns timing, tissue-specific effects, or altered renal function.

Common mistake: “The model predicts sodium excretion is higher, so extracellular volume and MAP must immediately fall by a calculable amount.” This confuses an excretion rate with a net balance and a pathway-specific response with a closed-loop outcome. The correction is to state what was held fixed, compare excretion with intake over time before inferring a volume trend, and account for vascular, cardiac, renal, and hormonal feedbacks before claiming a total pressure response.

Key takeaway: Stretch-dependent ANP and BNP signaling provides a plausible link from cardiac loading to renal natriuresis and vascular relaxation. Quantitative predictions are conditional on the release, response, and background assumptions; the sign of an isolated pathway does not determine the final integrated cardiorenal state.

Check Your Understanding

Work each problem from the stated model and assumptions. Solutions follow the questions.

  1. If \(x=0.5\), calculate \(\theta\). Does doubling \(x\) double \(\theta\)?
  2. With \(E_0=70\) mmol/day, \(E_{\max}=30\) mmol/day, and \(x=1\), calculate predicted sodium excretion.
  3. If peptide activation lowers systemic resistance while CVP and CO are held fixed, what happens to MAP in the lumped pressure-flow relation? Name one reason this may not equal the total long-term pressure change.
  4. A simulation predicts \(\theta=1.08\) from \(x=2\) using \(\theta=x/(1+x)\). Identify the inconsistency and give the correct activation.

Solutions. (1) \(\theta=0.5/1.5=0.3333\). Doubling \(x\) gives \(1/2=0.5000\), not twice the original activation; the response saturates. (2) \(\theta=1/(1+1)=0.5\), so \(E_{\mathrm{Na}}=70+30(0.5)=85\) mmol/day. (3) With \(MAP=CVP+CO\,R_{\mathrm{sys}}\) and fixed CVP and CO, lower resistance lowers MAP. Long-term pressure can differ because renal sodium handling, cardiac output, and compensatory feedbacks may change. (4) For \(x=2\), \(\theta=2/(1+2)=0.6667\). The reported value exceeds the required range zero to one and signals an equation or implementation error.

Clinical Biostats