1Department of Mathematics, Duke University, Durham, North Carolina 27708; and 2Department of Physiology and Biophysics, State University of New York, Stony Brook, New York 11794
Submitted 18 October 2002 ; accepted in final form 18 June 2003
![]() |
ABSTRACT |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
kidney; renal hemodynamics; autoregulation; mathematical model; nonlinear dynamics
A common and useful method of investigating TGF is to impose sustained perturbations of proximal tubule (PT) fluid flow in a freely flowing nephron where the TGF feedback loop is closed and functional (3, 5, 6, 20, 21, 26). Using this technique, Leyssac and collaborators (5, 17, 18) have demonstrated that LCO can be initiated or extinguished in halothane-anesthetized rats by insertion or removal of PT fluid, findings that indicate that PT fluid flow can have a substantial effect on the stability of the TGF system. Some insight into the basis of this phenomenon is provided by our previous investigations of the TGF-mediated LCO (1216, 22). These modeling studies indicate that LCO will emerge for sufficiently large feedback loop gain magnitude (12), that the parameter regime that supports LCO may overlap the parameter regime of normal TGF operation, that a gain magnitude near that needed for LCO will produce maximal feedback compensation (16), and that LCO may augment NaCl delivery to the distal nephron (16). A finding of fundamental importance is that the emergence of LCO depends on key TGF system parameters and variables that, in turn, depend on the physiological state of the nephron or on the animal as a whole. For example, the parameters may be affected by an increase in systemic blood pressure, a resetting of the TGF response curve, a primary change in tubular transport, or an experimental intervention such as microperfusion of fluid into the PT. We have previously presented numerical simulations that indicate that sustained PT flow perturbations can elicit LCO or extinguish LCO (16). However, such simulations cannot determine how the feedback loop gain or the gain threshold at which LCO emerge is influenced by PT flow perturbations.
The principal objectives of this study were to analyze, in the context of our mathematical model, the effect of sustained PT flow perturbations on feedback loop gain, on the gain threshold at which LCO emerge, and on feedback compensation. By means of this analysis we provide a theoretical framework for understanding how fluctuations in physiological variables impact the stability of the TGF system and its regulatory efficacy.
Glossary
Parameters
Independent Variables
Specified Functions
Dependent Variables
![]() |
METHODS |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Concepts and Terminology
The TGF system is a negative feedback loop in which afferent arteriolar (AA) smooth muscle tension is increased in response to increased chloride concentration in TAL luminal fluid passing by the macula densa (MD) (23). AA vasoconstriction will decrease filtration pressure in glomerular capillaries and thus reduce SNGFR. Reduced SNGFR reduces both fluid flow through the TAL and chloride concentration in luminal fluid passing the MD.
SNGFR and thus TAL flow may be nearly unchanging in time and thus approximate a time-independent steady state. Associated with that steady state is an SNGFR and a TAL chloride concentration at the MD that together are called the "operating point" of the TGF response function. Frequently, that operating point is found to be near the steepest portion of the TGF response function (1); in the standard form of that function (a logistic curve that is equivalent to a scaled hyperbolic tangent function), the point of steepest slope corresponds to the point where the response function changes from concave down to concave up, i.e., to the point of inflection.
Also associated with steady-state flow, and thus with its corresponding operating point, is feedback gain, a measure of feedback signal amplification. Roughly speaking, gain can be measured by breaking the TGF (signal) loop (e.g., breaking the signal loop at the beginning of the thick limb), increasing the feedback loop signal by a small amount (e.g., increasing TAL flow by FAL), and measuring the resulting change in the return signal (e.g., in a short-looped nephron, measuring the change in terminal descending limb flow,
FDL). Gain is then approximated by the resulting change divided by the signal increase (e.g., by
FDL/
FAL).
Two concepts of gain, steady-state gain (GSS) and instantaneous gain (-), have been applied to the TGF loop (13); each can be viewed as an interpretation of
FDL/
FAL (see APPENDIX A). Steady-state gain is the gain that has been measured in laboratory experiments; the concept of instantaneous gain arose from a theoretical analysis of a mathematical model of the TGF system (12, 13). We have previously shown that, under restricted circumstances, GSS is closely approximated by -
(13).
For a given set of model parameters (see model below), there is a particular steady-state feedback loop configuration with an associated operating point and with associated values of -GSS and . In the context of our model's normalized quantities, steady-state gain magnitude -GSS corresponds to the product of the magnitude of the slope of the feedback response function (see Eq. 2) at the operating point and the derivative of chloride concentration in luminal fluid alongside the MD with respect to TAL flow, also at the operating point;
corresponds to the product of the slope of the feedback response function at the operating point, the derivative of chloride concentration at the MD with respect to TAL axial position, and the steady-state TAL fluid transit time.
A model steady state may be either stable or unstable. If a steady state is stable (stable steady state, SSS), model solutions will more and more closely approximate that steady state after a transient perturbation (as, e.g., a transient addition to SNGFR). If a steady state is unstable (unstable steady state, USS), model solutions will more and more closely approximate LCO after a transient perturbation, no matter how small the perturbation.
For a sufficiently long delay in TGF signal transmission at the juxtaglomerular apparatus (JGA), the stability of a steady state is determined by the gain1 and by the critical gain
c (12). If
<
c, then the stable state is a (time-independent) steady state; if
>
c, then the stable state is an LCO. The critical gain
c can be determined by numerical experiments in which one varies
to find the boundary between the stable and unstable steady states or by use of a characteristic equation, which can be derived from our model and which can be used to find a general expression for
c in terms of model parameters.
A steady state (whether stable or not), and thus its associated operating point, may be altered by non-TGF-regulated factors that increase effective SNGFR (and thus TAL flow), e.g., an increase in systemic blood pressure or an experiment in which fluid is inserted into the PT. A steady state may also be altered by a change in the slope of the feedback response curve, if the operating point does not coincide with the point of inflection. A convenient parameter for characterizing the slope of the response curve is the feedback sensitivity , which for the purposes of this study is defined to be the magnitude of the slope of the feedback response curve at its inflection point. We define the critical sensitivity
c to be the sensitivity corresponding to a case in which
=
c. Thus, if
>
c, the model solution will tend toward an LCO; if
<
c, the model solution will tend toward the SSS.
A sustained flow perturbation or other sustained change in model parameters has the potential to change both and
c and thus change the relationship between them (i.e., change the sign of
-
c). Thus a stable state may be affected by a perturbation, so that an SSS may become a USS, resulting in a stable LCO, or a state characterized by LCO may become an SSS. Figure 3 (see below) summarizes the dependence of stable model behavior on gain, operating point, flow perturbations, and TGF sensitivity.
|
Mathematical Model
Our model is given by the following system of coupled equations
![]() | (1) |
![]() | (2) |
![]() | (3) |
Each equation is in nondimensional form (see APPENDIX B). The space variable x is oriented so that it extends from the entrance of the TAL (x = 0), through the outer medulla, and into the cortex to the MD (x = 1). Figure 1 gives a schematic representation of the model.
|
Equation 1 is a partial differential equation for the chloride concentration C in the intratubular fluid of the TAL of a short-looped nephron. We assume that fluid entering the TAL has constant chloride concentration; thus, we assume that C(0, t) always equals 1. At time t = 0, initial concentrations C(x, 0) (for x (0, 1]) and C(1, t) (for t
(-
, 0)) must be specified. For t > 0, the rate of change of that concentration depends on processes represented by the three right-hand terms in Eq. 1. The first term is axial advective chloride transport at intratubular flow speed F. The second is transepithelial efflux of chloride driven by active metabolic pumps situated in the tubular walls; that efflux is assumed to be approximated by Michaelis-Menten kinetics, with maximum transport rate Vmax and Michaelis constant KM. The third term is transtubular chloride backleak, which depends on a specified fixed extratubular chloride concentration profile Ce(x) and on chloride permeability P.
Equation 2 describes fluid speed in the TAL as a function of effective luminal chloride concentration CMD at the MD. This feedback relationship is an empirical equation well established by steady-state experiments (23). The constant CI is the inflection point of the TGF response curve; in our model, it is also the chloride concentration at the MD when F = 1 and the chloride concentration profile in the TAL has assumed a steady state. The positive constants K1 and K2 describe, respectively, the range of the feedback response and its sensitivity to deviations from the steady state. The constant represents a flow perturbation applied into Bowman's space. Owing to the nondimensional model formulation and to our assumption that a fixed fraction of the glomerular filtrate is delivered to the TAL,
can represent a fractional flow perturbation applied into Bowman's space or into the tubular lumen at the entrance of the TAL (in the nondimensional formulation, the TAL flow rate F and the SNGFR Q have the same value, because the dimensional F and Q have both been scaled by their respective base-case values).
Equation 3 represents time delays in the feedback pathway between the luminal fluid chloride concentration at the MD, C(1,t), and an effective MD concentration, CMD(t), which is used to calculate the flow response that is mediated by AA smooth muscle. In a quasi-steady state, Eq. 2 provides a good description of the TGF response. However, dynamic experiments (2) show that a change in MD concentration does not significantly affect AA muscle tension until after a discrete (or pure) delay time p > 0, and then the effect is distributed in time so that a full response requires additional time, with greatest weight in the time interval [t -
p -
, t -
p], where
> 0 is a second delay parameter. To simulate the pure delay followed by a distributed delay, the convolution integral given in Eq. 3 is used to describe the effective signal received by the AA at time t (22). The kernel function
for this integral is given by
![]() | (4) |
With this function, a step change in C results in a sigmoidal increase in CMD over a nondimensional time interval of .
Model parameters. A summary of parameters and variables, with their dimensional units as commonly reported, is given in the Glossary. The base-case parameters, which collectively represent the reference point for our parameter studies, are given in Table 1; the criteria for their selection and supporting references were given in Ref. 12. The extratubular chloride concentration is given in nondimensional form by Ce(x) = Co(A1e-A3x + A2), where A1 = (1 - Ce(1)/Co)/(1 - e-A3), A2 = 1 - A1,A3 = 2, and Ce(1) corresponds to a cortical interstitial concentration of 150 mM.2 A graph of Ce for F = 1 was given in Fig. 1 of Ref. 12. The steady-state operating MD chloride concentration that corresponds to = 0, which is CI, can be calculated numerically using the TAL dimensions and transport parameters, with steady flow F = 1 in Eq. 1 for the time interval equal to the washout time of the TAL at that flow, 1 dimensionless time unit.
|
Model solutions. For the case of no perturbation, i.e., = 0, a steady-state solution to Eqs. 14 may be obtained and sustained by fixing C(1, t) = CI for t
1, fixing F at the base-case value of 1 for the TAL transit time interval (here, t
[0,1]), solving Eq. 1 for that time interval (to obtain C(1, 1) = CI), and then closing the TGF loop so that F is computed via Eqs. 24 for t > 1. Then F(CMD(t)) will equal 1 for all time.
The steady state for
0 is more difficult to obtain because it corresponds to an unknown steady-state value of F, a value that in general we will call Fop (thus
= 0 implies Fop = 1, whereas
0 implies Fop
1). Because we have previously established that, if
p = 0 and
= 0 a model steady state is stable (12), Fop can be found directly by eliminating the delay at the MD (thus, CMD(t) = C(1,t)) and solving Eqs. 1 and 2 until the solutions are sufficiently close to steady state. The long-time values of CMD(t) and F(CMD) will closely approximate Cop and Fop.
Alternatively, one can consider the inverse problem in which one specifies Fop and finds the corresponding . Thus one can fix F = Fop in Eq. 1 for an interval equal to the TAL transit time at flow speed Fop (thus for t
[0,1/Fop]), find Cop = C(1, 1/Fop), assume that C(1, t) = Cop for t
1/Fop, and choose
so that for t > 1/Fop, F(CMD(t)) = Fop in Eq. 2. If the feedback loop is closed at t = 1/Fop, then F(CMD(t)) = Fop for all time.
When Fop is known, Cop may be obtained by solving the time-independent form of Eq. 1
![]() | (5) |
A dynamic solution to Eqs. 14 can be found by taking the corresponding steady-state solution as initial data and perturbing transiently. The long-time solution, depending on model parameters, will tend either to a steady state or to an LCO.
Steady-state gain. The steady-state gain, which has been measured in experiments (1, 3), is expressed (13) in terms of dimensionless model variables as
![]() | (6) |
Characteristic equation and instantaneous gain. Information about the stability of a steady-state solution of the model TGF system (i.e., how it will respond to a transient perturbation) can be obtained from the model's characteristic equation,
![]() | (7) |
![]() | (8) |
This characteristic equation generalizes previous versions (12, 15) by including the effect of Fop 1; previously, we assumed Fop = 1 so that Fop did not appear. The method of derivation of the characteristic equation was explained previously (12, 13, 22).
There are two ways in which the characteristic equation may be used. First, given a particular set of model parameters, one can find the corresponding steady-state solution S(x,Fop), from that solution obtain via Eq. 8, substitute that
in the characteristic equation, and then solve the characteristic equation (by techniques of numerical analysis) for the corresponding real and complex parts of the eigenvalues
n (n = 1, 2, 3,...) that satisfy the characteristic equation when
=
n. If all the real parts of the
n's are negative, then the model solution is an SSS. If the real part of any one of the eigenvalues is positive, then the steady state is a USS.
Conversely, given a set of model parameters, and the corresponding steady-state solution S(x, Fop), one can use the characteristic equation to solve for a particular value of , the critical gain
c, which corresponds to the transition between the SSS and the USS:
c is the smallest
such that the real part of an eigenvalue
n equals 0. Thus, if one is given an external concentration Ce, a steady-state profile S with associated Fop (a steady state based on the feedback response parameters), delay parameters
p and
, and backleak permeability P, then one can determine the critical gain
c, and compare it with the feedback gain
(computed by Eq. 8) to determine whether a steady-state solution is stable; this is how the characteristic equation was used in this study. Because a steady-state solution and its corresponding operating values Cop and Fop are altered by a sustained flow perturbation introduced through
, both
and
c depend on
(see Fig. 3).
Changing the physiological state: sustained flow perturbation and TGF sensitivity. In this study, we assume that the model physiological state is determined by two parameters. The first parameter, the sustained flow perturbation (which appears in Eq. 2), quantifies all non-TGF-related factors affecting SNGFR. Such factors include sustained microperfusion into the PT, changes in mean systemic arterial pressure, and changes in extracellular volume (ECV). The second parameter, the TGF sensitivity
, is defined to be the magnitude of the slope of the TGF response curve at its inflection point, which has been shown to depend on the physiological state (1, 23). By maximizing the slope of the TGF response given by Eq. 2 as a function of CMD, one finds that the dependence of model TGF sensitivity on model parameters is given by
![]() | (9) |
Feedback compensation. The efficacy of model TGF regulation can be quantified by calculating feedback compensation, an index used in experimental investigations (10, 20, 26). Feedback compensation is defined by
![]() | (10) |
![]() | (11) |
In the definition for magnification, Y is a system variable that is regulated by means of the feedback loop (e.g., distal chloride delivery); Y is the change in the system variable Y in response to a change
X (i.e., a perturbation) in another system variable X (e.g., PT flow). The denominator of Eq. 11, (
Y/
X)OL, is the ratio of
Y to
X in the case where the feedback loop is open (open-feedback-loop case, or OL). The numerator (
Y/
X)CL is the corresponding ratio when the loop is closed (closed-feedback-loop case, or CL). If
X tends to zero, M converges to a ratio of derivatives; however, the
-notation is retained because it is consistent with experimental studies, which necessarily entail measurable perturbations.
Because the stable state of the model TGF system may be an LCO or a stable steady state, and the type of stable state (LCO or steady state) may differ for X = 0 and
X
0, "
Y" requires further interpretation to be well defined. We thus adopt two conventions with respect to computing
Y: 1) specific values of Y are time averages that correspond to the stable state of the system under the given perturbation; and 2) changes in Y are normalized with respect to the stable, nonperturbed state of the system. We thus set
Y = (Y
- Y0)/Y0, where Y0 is the time-averaged value of Y obtained from the stable state when
X = 0, and Y
is the time-averaged value of Y obtained from the stable state when
X
0 (if the stable state of the system is an LCO for zero perturbation, Y
is compared to the zero perturbation stable-state average, even if the perturbation
X has resulted in an SSS). We normalize with respect to Y0 because Y0, the time-averaged value of Y obtained from the stable state when
X = 0, may differ for the OL and CL model TGF systems.
If a steady state is the stable state, distal chloride delivery is computed as the product of the steady-state MD chloride concentration and the steady-state TAL fluid flow rate at the MD. If an LCO is the stable state, distal chloride delivery is computed as the time average of the product of the chloride concentration at the MD and the TAL fluid flow rate at the MD.
Numerical methods. Details of the numerical methods used to obtain results are given in APPENDIX C. The base-case parameter values (excepting those for and
) were used in all calculations.
![]() |
RESULTS |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Because we aim to make predictions that are related to experiments, it is incumbent on us to provide evidence that the instantaneous gain that arises in our theoretical analysis can be estimated from experiments that measure steady-state gain.
Previously, we have calculated, for the case of no sustained flow perturbation and for base-case parameters, that the instantaneous gain magnitude exceeds steady-state gain magnitude -GSS by
10% (13). Here we report further calculations to demonstrate good agreement between instantaneous gain and steady-state gain in steady states corresponding to sustained flow perturbations over a range of feedback sensitivities. Some results of those calculations are illustrated in Fig. 2, which shows instantaneous and steady-state gain magnitudes
and -GSS for TGF sensitivities
4 and
8 (sensitivities that yield
= 4 and
= 8 at zero perturbation,
= 0). The value of
used in this figure and subsequently in this section is the dimensional value of a sustained flow perturbation in early PT.
|
The relative difference between the two measures of gain, at zero perturbation, is 10.3%, for both 4 and
8. The maximum relative difference between the two measures, for
4, is 13.4% at
= -12 nl/min; for
8, that difference is 12.5%, also at
= -12 nl/min. The maximum absolute difference between the two measures, for
4, is 0.374, at
= -0.060; for
8, that difference is 0.748 at
= -0.027. We computed the two measures of gain over the range of
= -12 to 12 nl/min (in increments) and for sensitivities ranging from
0.25 to
10 (in increments). As sensitivities increase, we found that the relative differences tend to stabilize near 10%. Also, we found that, although both measures of gain tend to zero as sensitivity tends to zero, the relative difference increases as sensitivity tends to zero and perturbations tend to -12 nl/min. Indeed, for sensitivities greater than
2 and perturbations greater than -8 nl/min, the maximum relative difference was 12.4%. However, in the "corner" region marked off by sensitivities
such that
1
2 and by perturbations
such that -10 nl/min
-8 nl/min, the maximum relative difference was 16.3%. The difference between the two gain values in this corner is not likely to be physiologically significant, however, because in that region both values are less than
1.2 and their difference is never larger than
0.14.
By means of an analysis like that in Ref. 13, we found that exceeds -GSS for all perturbations at a given sensitivity.
Shape of the -Curve
The dependence of on flow perturbations and feedback sensitivity can be understood by considering the two steps involved in calculating
(Fig. 3). Given a flow perturbation
and a feedback sensitivity
, one determines the operating steady-state TAL flow Fop, and then one evaluates the three factors comprising
(given by Eq. 8): F'(Cop), 1/Fop, and
.
The model's steady-state operating point is represented by Qop and Cop, where Qop is the effective SNGFR (which is defined to be the sum Q + ), and Cop is the MD chloride concentration that corresponds to the steady-state TAL flow Fop. The calculation of the steady-state operating point is represented graphically in Fig. 4; the operating point corresponds to the intersection of the MD chloride concentration curve and the TGF response curve. Figure 4, A and B, illustrates the variation of the intersection point with sustained perturbation
for two different feedback sensitivities
. A sustained flow perturbation affects the location of the intersection point by translating the TGF response curve vertically (see Eq. 2). The range of operating points is smaller for the case of higher feedback sensitivity, because the higher sensitivity represents a stronger TGF response and thus a response with greater feedback compensation. Figure 4C explicitly shows the dependence of Qop and Fop on
for the feedback sensitivities used in Fig. 4, A and B. For the higher sensitivity
8, Fop deviates less from the base-case TAL operating flow rate of 6 nl/min, compared with the lower sensitivity
4.
|
|
|
|
Shape of the c-curve. Figure 5A2 illustrates the dependence of critical gain
c on
for sensitivities
4 and
8, a dependence computed by means of the characteristic equation, Eq. 7 (
and
play a role in determining
c through their influence on the steady-state TAL flow Fop; see Fig. 3). The decrease in
c as a function of increasing perturbation largely arises from the accompanying increase in Fop, which increases the ratio of the JGA delay times (
p and
) to TAL fluid transit time, 1/Fop (based on an analysis similar to that given in Ref. 12). Because of the smaller variation of
c, relative to
, as a function of
(note the different scales on the vertical axes of Fig. 5, A1 and A2), variation of
primarily determines the sign of
-
c and thus the stable state of the model system (LCO or SSS).
Dependence of the Sign of -
c on
In Fig. 5, B1 and B2, the gain curves from Fig. 5, A1 and A2, are matched according to the sensitivities 4 or
8. For sustained flow perturbations between a and b or between c and d, the stable state is an LCO (
-
c > 0), and for
outside those regions the model solution tends toward an SSS (
-
c < 0). As TGF sensitivity
increases, the gain
changes more rapidly than critical gain
c, and therefore the range of values for which an LCO is the stable state increases as
increases. These results provide a theoretical explanation for a well-known observation: if the TGF response curve is sufficiently steep at its inflection point, and if the administered PT microperfusion is sufficiently small, then TGF mediates regular, sustained oscillations in tubular flow.
Critical Sensitivity Curve
The curve in Fig. 6A gives critical sensitivity as a function of perturbation (recall that, given a sensitivity
, if there is a sustained perturbation
such that the instantaneous gain
is equal to the critical gain
c, then that sensitivity
is defined to be the critical sensitivity
c for that
; in Fig. 5B1,
4 equals the critical sensitivity for
= a and
= b, and in Fig. 5B2,
8 is the critical sensitivity for
= c and
= d). Because both
and
c depend on
, the curve in Fig. 6A can be calculated directly only by means of exhaustive numerical simulations based on Eqs. 14. Therefore, we used the following implicit procedure to calculate points on the curve in Fig. 6A. We began by selecting a value of Fop; next, we calculated the corresponding
c via the characteristic equation (Eq. 7); finally, Eqs. 8 and 2 were used to determine the sensitivity
and the sustained perturbation
, respectively, needed to achieve the selected value Fop.
The points a, b, c, and d on the perturbation axis in Fig. 6A are the same as those shown in Fig. 5, B1 and B2. The critical sensitivity corresponding to points a and b in Fig. 6A,
4, corresponds to the sensitivity
4 in Fig. 5, and the critical sensitivity
8 corresponding to points c and d in Fig. 6A is the sensitivity
8 used in Fig. 5. Thus, Fig. 6A can be viewed as a summary of the results of many experiments like those shown in Fig. 5, B1 and B2, although the curve given in Fig. 6A was calculated by means of the efficient and accurate technique described in the previous paragraph. Above the curve, the stable solution is an LCO, whereas below and outside the curve the stable solution is a steady state.
Critical Sensitivity Curve: Tests and Examples
Figure 6B gives the results of two numerical experiments conducted to test the critical sensitivity curve in Fig. 6A. Values of the sustained perturbation and the TGF sensitivity
were selected to place the model TGF system just inside (at point S) or just outside (at point T) the critical sensitivity curve. For each case, the model system was initialized at the time-independent steady-state operating point determined by those values (see Fig. 3). Then, a transient flow perturbation was applied at an elapsed time of 6 s. For the state just inside the curve, an LCO developed; however, for the state just outside the curve, the oscillations initiated by the transient flow perturbation rapidly decreased in amplitude and the solution approached an SSS. These tests confirm the explicit analytical results summarized in Figs. 5 and 6A.
Figure 6C1 depicts the effects of two different sustained flow perturbations, = 5.0 and
= 9.5 nl/min, on the TGF response curve and on the corresponding operating points, indicated by the open and closed circles. The critical sensitivity curve in Fig. 6A and the points U and V corresponding to these perturbations predict that the stable state for
= 9.5 is a steady state (marked by the closed circle in Fig. 6C1) and that the stable state for
= 5.0 is an LCO (the associated unstable steady-state operating point is marked by the open circle in Fig. 6C1). The steady-state operating point for the smaller perturbation is closer to the steepest portion of the TGF response curve, compared with the operating point for the larger perturbation. Figure 6C2 illustrates the effect of suddenly changing (at time 4.0 min) the applied perturbation from
= 9.5 to
= 5.0 and then suddenly restoring the original perturbation (at time 8.0 min). Oscillations emerge and are damped out as predicted by the locations of points U and V in Fig. 6A. Note that the flow range of the oscillation, indicated by the thick gray curve in Fig. 6C1, occupies a significant portion of the sloped part of the TGF response function. This model simulation gives results similar to those obtained in experiments by Leyssac and Baumbach (18): in halothane-anesthetized rats, PT pressure oscillations can be initiated or terminated by adding fluid to (or removing fluid from) the late PT via a micropipette.
Figure 6D1, adapted from a study by Holstein-Rathlou and Leyssac (5), illustrates PT pressure from an experiment in a halothane-anesthetized Sprague-Dawley rat. The stable state approximated an LCO, but the removal of PT fluid (a perturbation of -10 nl/min) resulted in an SSS (albeit with a superimposed high-frequency oscillation from respiration). The removal of the perturbation resulted in a return to an LCO, and the response to the perturbation was repeatable. Figure 6D2 contains analogous results from our model. We began in a state, corresponding to X in Fig. 6A, in which the steady-state operating point is at the inflection point of the TGF response curve; moreover, because >
c, the stable state is an LCO. The introduction of a -10 nl/min flow perturbation3 resulted in the relocation of the steady-state operating point to the point corresponding to W in Fig. 6A, and that steady state was stable. The removal of the perturbation resulted in a return to an LCO. The flow pattern in Fig. 6D2 is very similar to the pressure pattern in Fig. 6D1, except that the model oscillation has a higher frequency. However, by a temporal rescaling of our base-case parameters, that aspect of the experimental record could be matched also.
Figure 6E1, also adapted from Ref. 5, illustrates PT pressure from an experiment in which a 7.5 nl/min flow perturbation increased the amplitude of oscillations in a preexisting LCO. Figure 6E2 shows results from our model that predict this pattern. If one initializes the model at a steady-state operating point that corresponds to a sustained negative perturbation and that is inside the LCO regime, but is near the bifurcation boundary (e.g., an operating point corresponding to point Y in Fig. 6A), and then one introduces a positive perturbation that results in a steady-state operating point nearer the point of inflection (e.g., corresponding to point Z in Fig. 6A); then the increased gain results in a stronger feedback response and an increased oscillation amplitude.
Relationship to Previous Work
In our previous model investigations (12, 13, 15, 22), Fop was always set to the (nondimensional) base-case value of 1. Thus Fop did not appear in the earlier versions of the characteristic equation, and the steady-state concentration profile S (which depends on Fop) was a fixed function. In Ref. 12 (with Fop equal to 1), we used the characteristic equation to determine the dependence of c on the JGA delay
.
In the present study, the delays p and
were fixed, and we studied the dependence of
c on Fop (more precisely, we studied how
c depends on
and
, which determine Fop). For each value of Fop, there is a curve, analogous to the n = 1 bifurcation curve in Ref. 12, that gives the dependence of
c on
=
p +
/2, as illustrated in Fig. 7. For
= 0 and any sensitivity, Fop equals 1, and therefore the middle curves in Fig. 7, A and B, are the same (these two middle curves are analogous to the curve labeled n = 1 in Fig. 4 of Ref. 12).
Along the Critical Sensitivity Curve, Chloride Compensation and Distal Chloride Delivery Change Abruptly as One Passes from Stable Steady State to Stable LCO
For the ranges of and
used in Fig. 6A, feedback compensation and delivery of chloride to the MD were computed as described in METHODS and APPENDIX C. Column B in Fig. 8 illustrates results for sensitivity
8. Figure 8B1 is a reproduction of Fig. 5B2, Fig. 8B2 is feedback compensation for chloride delivery to the MD corresponding to sensitivity
8, and Fig. 8B3 is chloride delivery to the MD for sensitivity
8. In Fig. 8, B2 and B3, the solid curve corresponds to the stable-state case (LCO or SSS); the dashed curve corresponds to the steady-state case where LCO have been suppressed in the USS regime by setting
p and
equal to 0.
|
Column A in Fig. 8, which illustrates model results for the stable-state case, is analogous to results in column B but also represents the role of variable , as shown by the axes accompanying Fig. 8A4. Column C illustrates analogous results for the case in which LCO are suppressed.
Figure 8A1 is analogous to Fig. 5, B1 and B2. The (mostly) upper surface is the gain as a function of sustained perturbation
and sensitivity
; the lower surface, largely hidden by the upper surface, is the critical gain
c. The bold line that appears on each surface corresponds to sensitivity
8. The intersection of these surfaces corresponds to the critical sensitivity curve exhibited in Fig 6A; specifically, the projection of that intersection onto the
-
plane, shown in Fig. 8A4, is the critical sensitivity curve. Figure 8C1 is the same figure as Fig. 8A1, but here the intersection of these surfaces is interpreted to mark the boundary between SSS and USS. In a physiological context, the unstable steady state is unrealizable, because frequent transient perturbations, always present in a living system, will result in LCO. However, the unrealizable USS is useful for evaluating the effects of stable oscillations, relative to steady-state operation.
Figure 8A2 illustrates model predictions for TGF compensation, assessed in terms of chloride delivery to the MD, assuming stable-state operation: compensation decreases abruptly when a change in or
moves the model physiological state from just outside to just inside the critical sensitivity curve. In contrast, steady-state operation yields a smooth relationship between compensation and the independent variables
and
, as shown in Fig. 8, B2 and C2. Thus, in the USS parameter regime, the effect of LCO, relative to steady-state operation, is to reduce the system's capacity to keep chloride delivery to the MD within a narrow range around the zero-perturbation chloride delivery.
Figure 8A3 illustrates model predictions for chloride delivery to the site of the MD. Comparison with Fig. 8, B3 and C3, indicates that, in the presence of LCO, perturbations result in an asymmetric delivery change. Positive perturbations result in increased chloride delivery relative to the USS, whereas negative perturbations result in decreased delivery relative to the USS.
Figure 9 shows that water delivery to the MD, as a function of perturbation, varies much less than chloride delivery to the MD, and that the variation of water compensation with respect to perturbation is more symmetrical than that of chloride compensation. For =
8, Fig. 9, A and B, shows water compensation and delivery as a function of perturbation, both for stable-state and for steady state, along with the stable-state result for chloride. Figure 9C shows the maximum difference, over sustained perturbations of -12 to 12 nl/min, between percent stable-state and percent steady-state delivery as a function of sensitivity, for both chloride and water: the model predicts that LCO affect MD chloride delivery much more than they affect water delivery.
|
The results in Figs. 8 and 9 generalize results in Ref. 16 by showing that the pattern reported there holds for a substantial range of physiologically plausible sensitivities and by relating the stability thresholds observed in Ref. 16 to the conceptual framework developed in Refs. 12 and 13.4
![]() |
DISCUSSION |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
In previous studies (1216), we have noted limitations of the model framework used in this study. Many simplifying assumptions are required to obtain a model formulation that is amenable to analytical treatment, e.g., derivation of a characteristic equation (Eq. 7). For example, fractional volume absorption in the PT and descending limb is assumed to be constant and equal to the parameter ; thus our model assumes perfect glomerulotubular balance (GTB). Moreover, by using a simple expression to represent saturable, concentration-dependent epithelial transport, we have avoided the detailed representation of tubular epithelial cells. Nonetheless, our previous studies have suggested that our model is effective in predicting and explaining the results of laboratory experiments. Thus we believe that the model results from this study provide insight into the relationships among key physiological parameters and the state of fluid flow in the nephron.
Because this work is specifically concerned with perturbations in proximal nephron flow, a concern that may be raised is that our assumption of perfect GTB may not be strictly applicable. However, we believe that this simplification is appropriate for our relatively simple model and for the specific goals of our investigation, for two reasons. First, the capability to solve and utilize a characteristic equation is of inestimable value in understanding the fundamental behavior of this system and in verifying numerical methods and simulation strategies. At present, it is unclear whether the analytical treatment afforded by our characteristic equation can be extended to encompass representations of flow-dependent PT absorption. Second, we do not know how PT transepithelial transport may be affected by TGF-mediated flow oscillations, which have short periods relative to the time intervals used in steady-state measurements. In contrast, both experimental and theoretical studies indicate that the TAL has a major influence on the dynamic behavior of the TGF loop in its action as a transducer of flow changes into concentration changes (6, 15).
However, notwithstanding our assumption of perfect GTB, this study offers substantive insight into the effects of changes in PT absorption on the stability and dynamics of TGF. The basis for this assertion is that a sustained perturbation in flow applied into Bowman's capsule will have the same effect on Henle's loop flow and, hence, TGF dynamics, as a change in PT absorption. Thus our graphical results can be used to predict the effects on TGF stability of perturbations that are ascribed in part or in whole to deviations from GTB. For example, a change in absolute PT absorption can be represented in Fig. 8 by means of the perturbation value
= -
/
, which results in a change -
in TAL flow. Nevertheless, an effort to develop a comprehensive, dynamic mathematical model of PT absorption is clearly warranted by recent studies (27, 30, 31), which indicate that the abnormal regulation, chronic adaptation, and dysfunction of PT absorption can play a major role in the pathogenesis of hemodynamic abnormalities in renal disease.
Holstein-Rathlou and Leyssac (5) have also used a mathematical model to simulate the effects of sustained perturbations on the stability of TGF-mediated LCO. However, our work extends their early work on this problem and is distinguishable from theirs in important ways. First, their model study was based on numerical simulations and used informal conceptual reasoning to explain the effects of perturbations. In contrast, in the context of our model and its characteristic equation, we have provided a rigorous analysis and explicit predictions in terms of parameter values for the effects of perturbations; our numerical simulations confirm the accuracy of the analysis. Second, Holstein-Rathlou and Leyssac included in their model a detailed submodel of glomerular filtration and somewhat less detailed submodels of AA dynamics and the PT; the dynamics of the TAL and TGF signaling delays were represented by an implicitly defined third-order lag function. In contrast, we included a detailed model of spatially distributed TAL dynamics, which our previous studies have shown to introduce substantial dynamic complexity (1216). We did not model filtration or AA dynamics but, as in the model by Holstein-Rathlou and Leyssac (5), we incorporated the standard empirical TGF response function, although in terms of tubular chloride concentration at the MD rather than fluid flow into the loop of Henle. The relative simplicity of our model allows us to formulate a rigorous and detailed conceptual framework and make predictions with an economy of assumptions. Finally, none of the stability regions identified by Holstein-Rathlou and Leyssac by means of numerical simulation corresponds to the effects of sustained flow perturbations (i.e., none is analogous to Fig. 6A). However, the experimental pressure records of LCO given by Holstein-Rathlou and Leyssac in Figs. 1014 of Ref. 5 provide invaluable examples of the phenomena predicted in our study, and the model simulations given in those figures are analogous to our simulations; indeed, Figs. 12 and 13 in Ref. 5 correspond to our Fig. 6, D and E, respectively.
|
A significant result of this study is the prediction that the feedback compensation afforded by the TGF system might be significantly reduced, compared with the compensation associated with steady-state5 flow, by LCO. Laboratory experiments appear to support this hypothesis: transfer function analyses relating renal arterial blood pressure to renal blood flow indicate that a decrease in compensation is associated with an 30-mHz oscillation, an oscillation that can reasonably be inferred to arise from TGF-mediated oscillations in many nephrons (8, 9). Our previous simulation studies using a fixed TGF sensitivity also indicated that LCO could lead to reduced TGF regulation of chloride delivery to the distal nephron (16). The present study extends these earlier findings by showing that this phenomenon is robust, in that it occurs over a wide range of physiologically reasonable TGF sensitivities (from
3 to
10, the sensitivities that yield feedback loop gain magnitudes of 3 and 10, respectively, when there is zero PT flow perturbation; see Fig. 8).
The traditional view is that TGF, in conjunction with GTB, acts to match the filtered load with the absorptive capacity of the nephron, thereby stabilizing both SNGFR and the delivery of sodium into the distal nephron. Furthermore, in response to disturbances in extracellular fluid volume, systemic neurohumoral mechanisms are activated that modify TGF sensitivity and operating point, as well as proximal and distal nephron sodium absorption, resulting in an appropriate adaptive change in distal sodium delivery and renal sodium excretion (19, 24, 29). However, the results of the present study suggest that TGF regulation of distal sodium delivery may be more complex, owing to the prediction that TGF regulatory efficacy can change rapidly (within tens of seconds) when the dynamic state of the nephron changes (e.g., from steady flow to LCO).
One example of this additional complexity is that the response to an increase in tubular fluid load will depend on both the initial dynamic state of the individual nephron and the amplitude of the perturbation. If the initial TGF sensitivity is sufficiently high and the TGF system is oscillating, then small perturbations will likely not change the dynamic state of the nephron, and the degree of feedback regulation of distal sodium delivery and SNGFR afforded by TGF will correspond to the reduced TGF system compensation with LCO. In contrast, a large perturbation sufficient to quench LCO will result in an increase in TGF feedback compensation, which should enhance the autoregulatory response, thereby protecting the renal microcirculation from large surges in blood flow and intravascular pressure.
Another, more involved example of the potential physiological importance of the dynamic state of the nephron is TGF adaption to changes in extracellular fluid volume. Figure 10 presents a graphical illustration of the transitions from euvolemia to ECV expansion and contraction. Volume expansion results in renal vasodilation, reduced TGF sensitivity, increased filtration, and reduced proximal and distal nephron sodium absorption. In terms of our model, these effects correspond to an increase in and a decrease in
, whereas opposite changes are associated with volume contraction. As shown in Fig. 10A, when
increases, the TGF response curve translates upward, and when
decreases, the slope of the response curve at its inflection point decreases. During volume expansion, the model physiological state follows a path in the
-
plane similar to that given by the dashed straight line in Fig. 10B. The resulting effect on distal sodium delivery is illustrated in Fig. 10B by the relative thicknesses of the black and gray regions surrounding the dashed straight line. This is a qualitative representation of the difference between the distal NaCl deliveries associated with steady-state and stable oscillatory flow (see Fig. 9B). Given a point (
,
) along the dashed straight line, the width of the black path represents the magnitude of the distal NaCl delivery associated with steady-state flow, and the width of the gray path represents the magnitude of the distal NaCl delivery associated with stable oscillatory flow.
The model predicts that, after volume expansion (or contraction), NaCl delivery changes immediately and significantly, and if the model physiological state does not leave the interior of the critical sensitivity curve after the perturbation, the presence of LCO accelerates the TGF system's response to the initiating challenge. That is, LCO result in distal NaCl delivery changes that are both appropriate for the ECV status and larger than the changes that would occur under steady-state flow conditions, thereby shortening the time required for a return to euvolemia. The open circles in Fig. 10B represent states of volume expansion and contraction that are sufficient to suppress LCO and thus drive the model physiological state outside the critical sensitivity curve (the solid curve in the middle of Fig. 10B). The subsequent response of systemic neurohumoral systems (e.g., renin-angiotensin system, renal sympathetic nerve activity, ANF) will enhance distal sodium delivery. In addition, as illustrated in Fig. 10A, the TGF response curve will be reset horizontally over the following 1060 min, which will act to restore the TGF operating point near the point of maximum TGF sensitivity (26, 28, 29). Concurrent with the horizontal resetting of the TGF response curve, the critical sensitivity curve resets horizontally and changes slightly. Figure 10B illustrates the resetting of the critical sensitivity curve (the dashed curves that represent the reset critical sensitivity are not exact and are for illustrative purposes only). After resetting of the critical sensitivity curve, the open circles are again inside the region of stable oscillations. The reestablishment of LCO will then again facilitate the excretion of the excess sodium during volume expansion and enhance sodium conservation during volume contraction. The model thus suggests that an additional effect of TGF resetting may be to maintain the dynamic state of the nephron when extracellular fluid volume is altered.
As discussed above, volume expansion, as well as acute and chronic hypertension, is associated with reduced proximal absorption. However, in the present version of our model, fractional proximal and descending limb volume absorption (represented by the parameter ) was kept constant. Thus more detailed consideration of the physiological impact of LCO on the TGF system will require
to be variable and flow dependent, a modification that may preclude an analytical treatment of model system stability. We predict that this modification to our model will increase the relative contribution of LCO to the regulation of distal sodium delivery, owing to the synergistic effects of reduced upstream volume absorption and reduced TGF compensation when LCO are present. We also note that the effect of LCO on sodium excretion, and thus on sodium balance, depends on the correlation between distal sodium delivery and sodium excretion. Whether LCO-related changes in distal sodium delivery correspond to changes in sodium excretion depends on factors such as the response of distal tubular epithelial cells to LCO in tubular flow rate (16).
The results presented here suggest that the TGF system may display a wider range of functional behavior than previously thought. The differences in distal sodium delivery associated with steady-state and oscillatory flow may give the TGF system the flexibility needed to participate in the regulation of sodium balance while at the same time homeostatically regulating GFR. Further investigation with more detailed modeling and clever experimentation will be required to evaluate the physiological importance of the nonlinear dynamics of the TGF system. The analytical strategy and the numerical tools developed in this study will greatly facilitate such efforts.
![]() |
APPENDIX A |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Steady-state gain, denoted by GSS, corresponds to a case in which one introduces a sustained step perturbation in the feedback loop signal (here taken to be fluid flow rate) and waits sufficiently long for a new steady state to be established before one quantifies the return signal (13). Then, GSS is approximated by the resulting change divided by the signal increase; that approximation approaches the exact value of GSS as the signal increase tends to zero. Thus, in a short-looped nephron, if the feedback loop is broken at the entrance to the TAL, the feedback signal is the tubular flow rate, and GSS is given by the limit
![]() | (A1) |
Instantaneous gain, denoted by -, corresponds to a case in which one instantaneously increases the feedback response through the JGA by means of a change in the feedback loop signal and obtains the full corresponding change in the return signal, subsequent to delays in the feedback pathway (13). Then, -
is the value of the change in the return signal divided by the signal increase as that increase tends to zero. Thus, if the feedback loop is broken at the entrance to the TAL, one may use the same limit expression as given above in Eq. A1, but 1)
FAL is interpreted to be the fractional magnitude of an instantaneous insertion of fluid in the TAL entrance (i.e., a Dirac delta function perturbation, which results in an instantaneous shift of the fluid column) and
FDL is interpreted to be the fractional response in descending limb flow, assuming that the TAL chloride concentration at the MD arising from the fluid insertion remains fixed; and 2) delays at the JGA are ignored.
![]() |
APPENDIX B |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
The dimensional forms of Eqs. 1 and 2 are given by
![]() | (B1) |
![]() | (B2) |
To express these equations in nondimensional form, let , and
, where L is the TAL length and the quantities subscripted with an "o" are conveniently chosen reference values: Ao =
r2, to = AoL/Fo, Co = C(0, t), Fo = FI =
QI, (Vmax)o = FoCo/(2
rL), Po = Fo/(2
rL), and Qo = QI. With these conventions, to is the filling time (and thus the transit time) of the TAL at flow rate Fo, and (Vmax)o is the rate of solute advection into the inlet of the TAL at flow rate Fo, divided by the area of the walls of the TAL. When Eqs. B1 and B2 are rewritten in dimensionless terms and the tilde symbols are dropped, Eqs. 1 and 2 follow directly. The dimensional forms of Eqs. 3 and 4 are the same as their nondimensional forms.
The instantaneous gain magnitude has dimensional form
![]() | (B3) |
![]() |
APPENDIX C |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Methods are identified by corresponding figure numbers.
Figure 2. Data for these two curves were calculated using Eqs. 6 and 8. The steady-state operating points, which are needed in Eqs. 6 and 8, were calculated from the values of and
using an iterative procedure. Given an initial estimate for Cop, an estimate for Fop was calculated using Eq. 2, and this approximate Fop was then substituted for F in Eq. 5, which was solved using the standard explicit fourth-order Runge-Kutta method (RKM) to obtain S(1, F). The approximate Cop for the next iteration was taken to be a weighted average of the previous estimate and this value of S(1, F), thus guaranteeing convergence of the iterates to the actual steady-state operating values Cop and Fop. The factor
in Eq. 6 was calculated using the five-point centered differentiation formula (CDF), and the factor
in Eq. 8 was calculated using the five-point right-hand differentiation formula (RDF). In both Eqs. 6 and 8, F'(Cop) could be calculated analytically.
Figure 4. FIGURE 4, A and B. The TGF response curves were calculated using Eq. 2, with the nondimensional TAL flow values scaled by Qo = 30 nl/min. The MD chloride concentration curves were calculated by multiplying the effective SNGFR by to yield the TAL flow rate F, substituting this value of F in Eq. 5 and solving for S(1, F), using RKM.
FIGURE 4C. The steady-state operating TAL flow rates Fop were calculated using the iterative method described for Fig. 2. The corresponding values of Qop were given by Qop = Fop/.
Figure 5. FIGURE 5A1. Instantaneous gain values were calculated by the method described for Fig. 2.
FIGURE 5A2. Critical gain values were calculated using Eq. 7. Given values of and
, the steady-state operating TAL flow rate Fop and steady-state TAL chloride concentration profile S(x, Fop) were calculated using the iterative method described for Fig. 2. Equation 7 was then solved for
, after setting
= ±i
, using a method based on that summarized in APPENDIX C of Ref. 13.
FIGURE 5, B1 and B2. These curves are the same as those given in Fig. 5, A1 and A2, but they are matched according to TGF sensitivity .
Figure 6. FIGURE 6A. The critical sensitivities in this figure were calculated using the following algorithm. Given an operating TAL flow rate Fop, calculate the steady-state TAL chloride concentration profile S(x, Fop) (using Eq. 5 and RKM), set Cop = S(1, Fop), and calculate using RDF. Then, after calculating the critical gain
c associated with Fop (by solving Eq. 7 by the method described for Fig. 5A2), calculate
(Eq. 8, with
=
c and m = F'(Cop)); m is the required slope of the TGF response curve at the operating MD chloride concentration Cop. Next, solve the equation m = F'(Cop), using Newton's method, for the necessary sensitivity parameter K2 and set
c = K1K2. Finally, calculate the sustained perturbation
associated with Fop and
c using Eq. 2; i.e., let
= Fop - (1 + K1 tanh(K2(CI - Cop))).
FIGURE 6B. To obtain data for Fig. 6B, Eq. 1 was solved for values of and
that correspond to points S and T in Fig. 6A. After initialization at the time-independent steady state associated with points S and T, an additional transient Gaussian flow perturbation was applied at t
6 s. Equation 1 was solved by means of a second-order, essentially nonoscillatory (ENO) scheme, coupled with Heun's method for the time advance. This algorithm yields solutions that exhibit second-order convergence in both space and time (11). The numerical time and space steps in normalized units were
x = 1/640 and
t = (320 x t0)-1, where t0 is the base-case steady-state TAL transit time in seconds. These mesh widths, which correspond to dimensional values of
x = 7.8125 x 10-4 cm and
t = 3.125 x 10-3 s, were used for all dynamic calculations. The high degree of numerical grid refinement was required both to faithfully represent the nonlinearities that are embodied in the model equations (15) and to compute with sufficient accuracy the time-averaged chloride and water delivery rates used in the calculation of compensation values.
FIGURE 6 C1. The TGF response curves and the MD chloride concentration curve were calculated using the method described for Fig. 4,A and B.
FIGURE 6, C2, D2, and E2. Data for these figures were calculated by solving Eq. 1 for values of and
that correspond to the points U, V, W, X, Y, and Z in Fig. 6A, using the method described for Fig. 6B. For each of these figures, the model solution was stable at time t = 0.
Figure 7. For each fixed value of , data for these curves were calculated from the given values of
and
using the method described for Fig. 5A2. The value of
=
p +
/2 was changed by making proportional changes in the values of
p and
; in particular, the ratio
p/
= 2/3 was maintained, so that
p = 4
/7 and
= 6
/7.
Figure 8. FIGURE 8, A1, B1, and C1. Data for these figures were calculated using the methods described for Figs. 2 and 5A2.
FIGURE 8, A2 and A3. For a given sustained perturbation and TGF sensitivity
, average stable-state distal chloride deliveries were used to calculate chloride compensation (according to Eq. 10) and chloride percent distal delivery (100% x (Y
- Y0)/Y0). Closed-loop stable states were found by solving Eq. 1 for sufficiently large simulated time; i.e., for enough time either to allow transient behavior to die out or let stable oscillatory behavior fully develop (Eq. 1 was solved using the method described for Fig. 6B). If an LCO was the stable state, the waveform was recorded only after the oscillation had reached full amplitude and was indistinguishable from the limit-cycle trajectory (the model equations were solved for a minimum of 10 simulated minutes, and then two-cycle averages were calculated at 2-min intervals until the relative change in those averages was, in most cases, <0.0001). Open-loop stable states were found by solving Eq. 5, using RKM, with F = 1 +
. For stable states that were LCOs, the average chloride delivery was calculated, using the trapezoidal rule, over two complete cycles of the oscillation. For stable steady states, the average chloride delivery was equal to the stable constant value. The zero-perturbation value of chloride compensation was interpreted to be the average of the compensation values for
= ±0.25 nl/min.
FIGURE 8, C2 and C3. For a given sustained perturbation and TGF sensitivity
, steady-state distal chloride deliveries were used to calculate steady-state chloride compensation (according to Eq. 10) and steady-state chloride percent delivery (100% x (Y
- Y0)/Y0). Closed-loop steady states were found using the method described for Fig. 2, and open-loop steady states were found by solving Eq. 5, using RKM, with F = 1 +
. The zero-perturbation value of steady-state chloride compensation was found by interpreting Eq. 11 as the ratio of two derivatives and using CDF to calculate the magnification M.
FIGURE 8, B2 and B3. These curves were taken from Fig. 8,A2, A3, C2, and C3 ( =
8).
FIGURE 8, A4 and C4. This curve is the critical sensitivity curve from Fig. 6A.
Figure 9. FIGURE 9, A and B. Compensation and percent delivery values were calculated as described for Fig. 8.
FIGURE 9C. Maximum difference was calculated between percent stable-state delivery and percent steady-state delivery, for both chloride and water.
![]() |
DISCLOSURES |
---|
![]() |
ACKNOWLEDGMENTS |
---|
Portions of this work were presented in poster format at Experiment Biology 2001 (Abstract No. 404.14, FASEB J 15: A449, 2001).
![]() |
FOOTNOTES |
---|
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
1 In a negative feedback system, gain is negative by convention. However, it is usually more convenient to think in terms of gain magnitude , and for simplicity we will frequently refer to
as the "gain."
2 This value, which is somewhat higher than typical cortical interstitial Cl- concentration, and which approximates interstitial Na+ concentration, is used because we model transport of only a single solute, taken to be NaCl, and we assume that transport fluxes of Cl- and Na+ are identical.
3 Because our flow perturbation is applied into Bowman's space, our -10 nl/min perturbation would correspond to a somewhat smaller negative flow perturbation applied in the late PT, owing to proximal absorption. Thus, in our simulations, the TGF system is perturbed somewhat less than in Ref. 5; this difference, however, does not affect our results qualitatively.
4 A small difference in the calculation algorithm for stable-state compensation and percent delivery, relative to the calculations in Ref. 16, should be acknowledged. For stable-state compensation and percent delivery in this study, in calculating the difference between the deliveries associated with perturbed and nonperturbed states, the zero-perturbation reference was always taken to be the delivery associated with the zero-perturbation stable state. In contrast, in Ref. 16, the state of the perturbed solution (steady state or LCO) was noted and then compared with the corresponding state (steady state or LCO) at zero perturbation. For steady-state compensation and percent delivery, in both this study and in Ref. 16, the time-independent steady states (stable or not) were compared with steady-state zero-perturbation reference values. As a result of this small methodological difference, the stable-state compensation and delivery values do not agree exactly with the steady-state values for perturbations that push the system outside the critical sensitivity curve, but the disagreement, which is visible at the left and right of Fig. 8, B2 and B3, (for perturbations of magnitude greater than 8 nl/min) is slight. Although it is not clear which standard is more appropriate for comparison, the results indicate that the differing standards yield similar results and essentially the same predictions.
5 The steady-state compensation curves shown in Figs. 8, B2 (dashed curve) and C2, and 9A (dashed curve) are roughly comparable to the "fractional compensation" curves published by Thomson and Blantz (25). Both our compensation curves and Thomson and Blantz's fractional compensation curves imply that, assuming steady-state flow, when the operating point of the TGF system is near the inflection point of the TGF response curve, the TGF system exhibits its strongest feedback and is most able to compensate for flow perturbations.
![]() |
REFERENCES |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|