No abstract

Layton, Anita T., and Harold E. Layton. A region-based mathematical model of the urine concentrating mechanism in the rat outer medulla. I. Formulation and base-case results. Am J Physiol Renal Physiol 289: F1346 -F1366, 2005. First published May 24, 2005 doi:10.1152/ajprenal.00346.2003.-We have developed a highly detailed mathematical model for the urine concentrating mechanism (UCM) of the rat kidney outer medulla (OM). The model simulates preferential interactions among tubules and vessels by representing four concentric regions that are centered on a vascular bundle; tubules and vessels, or fractions thereof, are assigned to anatomically appropriate regions. Model parameters, which are based on the experimental literature, include transepithelial transport properties of short descending limbs inferred from immunohistochemical localization studies. The model equations, which are based on conservation of solutes and water and on standard expressions for transmural transport, were solved to steady state. Model simulations predict significantly differing interstitial NaCl and urea concentrations in adjoining regions. Active NaCl transport from thick ascending limbs (TALs), at rates inferred from the physiological literature, resulted in model osmolality profiles along the OM that are consistent with tissue slice experiments. TAL luminal NaCl concentrations at the corticomedullary boundary are consistent with tubuloglomerular feedback function. The model exhibited solute exchange, cycling, and sequestration patterns (in tubules, vessels, and regions) that are generally consistent with predictions in the physiological literature, including significant urea addition from long ascending vasa recta to inner-stripe short descending limbs. In a companion study (Layton AT and Layton HE. Am J Physiol Renal Physiol 289: F1367-F1381, 2005, the impact of model assumptions, medullary anatomy, and tubular segmentation on the UCM was investigated by means of extensive parameter studies. kidney; countercurrent multiplication; countercurrent exchange; NaCl transport; urea transport CONCENTRATED URINE, i.e., urine having an osmolality exceeding that of blood plasma, is produced by the absorption of water, in excess of solute, from the medullary collecting ducts (CDs). In the outer medulla (OM), water absorption from CDs is driven by vigorous active transport of NaCl from the thick ascending limbs (TALs) into the interstitium; at each medullary level, this transport results in an osmolality difference between the TALs and other medullary structures. This difference, frequently called the "single effect," is augmented (or "multiplied") by the countercurrent flow configuration of renal tubules and vasa recta, to generate a substantial osmolality gradient along the corticomedullary axis, a gradient that is believed to be common to all OM structures. In the inner medulla (IM), however, the means by which water is absorbed from CDs remains undetermined (20, 67).Substantial effort has been directed to constructing mathematical models...

Recent micropuncture studies in rats have demonstrated the existence of oscillatory states in nephron filtration mediated by tubuloglomerular feedback (TGF). We develop a minimal mathematical model of the TGF system, consisting of a first-order hyperbolic partial differential equation describing thick ascending limb (TAL) NaCl reabsorption and an empirical feedback relation. An analytic bifurcation analysis of this model provides fundamental insight into how oscillatory states depend on the physiological parameters of the model. In the special case of no solute backleak in the TAL, the emergence of oscillations explicitly depends on two nondimensional parameters. The first corresponds to the delay time of the TGF response across the juxtaglomerular apparatus, and the second corresponds to the product of the slope of the TGF response curve at the steady-state operating point and the space derivative of the steady-state NaCl concentration profile in the TAL at the macula densa. Numerical calculations for the case without TAL backleak are consistent with this result. Numerical simulation of the more general case with TAL backleak shows that the bifurcation analysis still provides useful predictions concerning nephron dynamics. With typical parameter values, the analysis predicts that the TGF system will be in oscillatory state. However, the system is near enough to the boundary of the nonoscillatory region so that small changes in parameter values could result in nonoscillatory behavior.

We used a mathematical model of the urine concentrating mechanism of rat inner medulla (IM) to investigate the implications of experimental studies in which immunohistochemical methods were combined with three-dimensional computerized reconstruction of renal tubules. The mathematical model represents a distribution of loops of Henle with loop bends at all levels of the IM, and the vasculature is represented by means of the central core assumption. Based on immunohistochemical evidence, descending limb portions that reach into the papilla are assumed to be only moderately water permeable or to be water impermeable, and only prebend segments and ascending thin limbs are assumed to be NaCl permeable. Model studies indicate that this configuration favors the targeted delivery of NaCl to loop bends, where a favorable gradient, sustained by urea absorption from collecting ducts, promotes NaCl absorption. We identified two model modes that produce a significant axial osmolality gradient. One mode, suggested by preliminary immunohistochemical findings, assumes that aquaporin-1-null portions of loops of Henle that reach into the papilla have very low urea permeability. The other mode, suggested by perfused tubule experiments from the literature, assumes that these same portions of loops of Henle have very high urea permeabilities. Model studies were conducted to determine the sensitivity of these modes to parameter choices. Model results are compared with extant tissue-slice and micropuncture studies.

Single-nephron proximal tubule pressure in spontaneously hypertensive rats (SHR) can exhibit highly irregular oscillations similar to deterministic chaos. We used a mathematical model of tubuloglomerular feedback (TGF) to investigate potential sources of the irregular oscillations and the corresponding complex power spectra in SHR. A bifurcation analysis of the TGF model equations, for nonzero thick ascending limb (TAL) NaCl permeability, was performed by finding roots of the characteristic equation, and numerical simulations of model solutions were conducted to assist in the interpretation of the analysis. These techniques revealed four parameter regions, consistent with TGF gain and delays in SHR, where multiple stable model solutions are possible: 1) a region having one stable, time-independent steady-state solution; 2) a region having one stable oscillatory solution only, of frequency f1; 3) a region having one stable oscillatory solution only, of frequency f2, which is approximately equal to 2f1; and 4) a region having two possible stable oscillatory solutions, of frequencies f1 and f2. In addition, we conducted simulations in which TAL volume was assumed to vary as a function of time and simulations in which two or three nephrons were assumed to have coupled TGF systems. Four potential sources of spectral complexity in SHR were identified: 1) bifurcations that permit switching between different stable oscillatory modes, leading to multiple spectral peaks and their respective harmonic peaks; 2) sustained lability in delay parameters, leading to broadening of peaks and of their harmonics; 3) episodic, but abrupt, lability in delay parameters, leading to multiple peaks and their harmonics; and 4) coupling of small numbers of nephrons, leading to multiple peaks and their harmonics. We conclude that the TGF system in SHR may exhibit multistability and that the complex power spectra of the irregular TGF fluctuations in this strain may be explained by switching between multiple dynamic modes, temporal variation in TGF parameters, and nephron coupling.