Mathematical model of an avian urine concentrating mechanism

Mathematical modeling of urea transport in the kidney

Mathematical modeling techniques have been useful in providing insights into biological systems, including the kidney. This article considers some of the mathematical models that concern urea transport in the kidney. Modeling simulations have been conducted to investigate, in the context of urea cycling and urine concentration, the effects of hypothetical active urea secretion into pars recta. Simulation results suggest that active urea secretion induces a "urea-selective" improvement in urine concentrating ability. Mathematical models have also been built to study the implications of the highly structured organization of tubules and vessels in the renal medulla on urea sequestration and cycling. The goal of this article is to show how physiological problems can be formulated and studied mathematically, and how such models may provide insights into renal functions.

Mathematical modeling of kidney transport

In addition to metabolic waste and toxin excretion, the kidney also plays an indispensable role in regulating the balance of water, electrolytes, nitrogen, and acid-base. In this review, we describe representative mathematical models that have been developed to better understand kidney physiology and pathophysiology, including the regulation of glomerular filtration, the regulation of renal blood flow by means of the tubuloglomerular feedback mechanisms and of the myogenic mechanism, the urine concentrating mechanism, epithelial transport, and regulation of renal oxygen transport. We discuss the extent to which these modeling efforts have expanded our understanding of renal function in both health and disease.

Modeling Transport and Flow Regulatory Mechanisms of the Kidney

The kidney plays an indispensable role in the regulation of whole-organism water balance, electrolyte balance, and acid-base balance, and in the excretion of metabolic wastes and toxins. In this paper, we review representative mathematical models that have been developed to better understand kidney physiology and pathophysiology, including the regulation of glomerular filtration, the regulation of renal blood flow by means of the tubuloglomerular feedback mechanisms and of the myogenic mechanism, the urine concentrating mechanism, and regulation of renal oxygen transport. We discuss how such modeling efforts have significantly expanded our understanding of renal function in both health and disease.

Systems biology of kidney diseases

Kidney diseases manifest in progressive loss of renal function, which ultimately leads to complete kidney failure. The mechanisms underlying the origins and progression of kidney diseases are not fully understood. Multiple factors involved in the pathogenesis of kidney diseases have made the traditional candidate gene approach of limited value toward full understanding of the molecular mechanisms of these diseases. A systems biology approach that integrates computational modeling with large-scale data gathering of the molecular changes could be useful in identifying the multiple interacting genes and their products that drive kidney diseases. Advances in biotechnology now make it possible to gather large data sets to characterize the role of the genome, epigenome, transcriptome, proteome, and metabolome in kidney diseases. When combined with computational analyses, these experimental approaches will provide a comprehensive understanding of the underlying biological processes. Multiscale analysis that connects the molecular interactions and cell biology of different kidney cells to renal physiology and pathology can be utilized to identify modules of biological and clinical importance that are perturbed in disease processes. This integration of experimental approaches and computational modeling is expected to generate new knowledge that can help to identify marker sets to guide the diagnosis, monitor disease progression, and identify new therapeutic targets.

Maximum urine concentrating capability in a mathematical model of the inner medulla of the rat kidney

In a mathematical model of the urine concentrating mechanism of the inner medulla of the rat kidney, a nonlinear optimization technique was used to estimate parameter sets that maximize the urine-to-plasma osmolality ratio (U/P) while maintaining the urine flow rate within a plausible physiologic range. The model, which used a central core formulation, represented loops of Henle turning at all levels of the inner medulla and a composite collecting duct (CD). The parameters varied were: water flow and urea concentration in tubular fluid entering the descending thin limbs and the composite CD at the outer-inner medullary boundary; scaling factors for the number of loops of Henle and CDs as a function of medullary depth; location and increase rate of the urea permeability profile along the CD; and a scaling factor for the maximum rate of NaCl transport from the CD. The optimization algorithm sought to maximize a quantity E that equaled U/P minus a penalty function for insufficient urine flow. Maxima of E were sought by changing parameter values in the direction in parameter space in which E increased. The algorithm attained a maximum E that increased urine osmolality and inner medullary concentrating capability by 37.5% and 80.2%, respectively, above base-case values; the corresponding urine flow rate and the concentrations of NaCl and urea were all within or near reported experimental ranges. Our results predict that urine osmolality is particularly sensitive to three parameters: the urea concentration in tubular fluid entering the CD at the outer-inner medullary boundary, the location and increase rate of the urea permeability profile along the CD, and the rate of decrease of the CD population (and thus of CD surface area) along the cortico-medullary axis.

Urine concentration and avian aquaporin water channels

Although birds and mammals have evolved from primitive tetrapods and advanced divergently, both can conserve water by producing hyperosmotic urine. Unique aspects in the avian system include the presence of loopless and looped nephrons, lack of the thin ascending limb of Henle's loop, a corticomedullary osmotic gradient primarily consisting of NaCl without contribution of urea, and significant postrenal modification of final urine. The countercurrent multiplier mechanism operates between the descending and ascending limbs of Henle via recycling of a single solute (NaCl) with no water accompaniment, forming an osmotic gradient along the medullary cone. Bird kidneys and developing rat kidneys share morphological and functional characteristics. Avian kidneys express aquaporin (AQP) 1, 2, and 4 homologues that share considerable homology with mammalian counterparts, but their distribution and function may not be the same. AQP2 expression in Japanese quail (q) evolves in the collecting duct of early metanephric kidneys and continues to increase in intensity and distribution during nephrogenesis and maturation. qAQP2 mRNA and protein are increased by arginine vasotocin (avian ADH), but vasotocin-induced enhancement of cAMP production and water permeability are less marked than in mammalian kidneys. Nephrogenesis is delayed by insufficient nutrition in avian embryos and newborns and results in fewer nephrons and an impaired water balance in adults. Diabetes insipidus quail with homozygous autosomal recessive mutation and an unaffected vasotocin system have low AQP2 expression, underdeveloped medullary cones. Comparative studies will provide important insight into integrative, cellular, and molecular mechanisms of epithelial water transport and its control by humoral, neural, and hemodynamic mechanisms.

An optimization algorithm for a distributed-loop model of an avian urine concentrating mechanism

To better understand how the avian kidney's morphological and transepithelial transport properties affect the urine concentrating mechanism (UCM), an inverse problem was solved for a mathematical model of the quail UCM. In this model, a continuous, monotonically decreasing population distribution of tubes, as a function of medullary length, was used to represent the loops of Henle, which reach to varying levels along the avian medullary cones. A measure of concentrating mechanism efficiency - the ratio of the free-water absorption rate (FWA) to the total NaCl active transport rate (TAT) - was optimized by varying a set of parameters within bounds suggested by physiological experiments. Those parameters include transepithelial transport properties of renal tubules, length of the prebend enlargement of the descending limb (DL), DL and collecting duct (CD) inflows, plasma Na(+) concentration, length of the cortical thick ascending limbs, central core solute diffusivity, and population distribution of loops of Henle and of CDs along the medullary cone. By selecting parameter values that increase urine flow rate (while maintaining a sufficiently high urine-to-plasma osmolality ratio (U/P)) and that reduce TAT, the optimization algorithm identified a set of parameter values that increased efficiency by approximately 60% above base-case efficiency. Thus, higher efficiency can be achieved by increasing urine flow rather than increasing U/P. The algorithm also identified a set of parameters that reduced efficiency by approximately 70% via the production of a urine having near-plasma osmolality at near-base-case TAT. In separate studies, maximum efficiency was evaluated as selected parameters were varied over large ranges. Shorter cones were found to be more efficient than longer ones, and an optimal loop of Henle distribution was found that is consistent with experimental findings.

Role of structural organization in the urine concentrating mechanism of an avian kidney

The organization of tubules and blood vessels in the quail medullary cone is highly structured. This structural organization may result in preferential interactions among tubules and vessels, interactions that may enhance urine concentrating capability. In this study, we formulate a model framework for the urine concentrating mechanism of the quail kidney. The model simulates preferential interactions among renal tubules by representing two concentric cores and by specifying the fractions of tubules assigned to each of the concentric cores. The model equations are based on standard expressions for transmural transport and on solute and water conservation. Model results suggest that the preferential interactions among tubules enhance the urine concentration capacity of short medullary cones by reducing the diluting effect of the descending limbs on the region of the interstitium where the collecting ducts are located; however, the effects on longer cones are unclear.

A region-based mathematical model of the urine concentrating mechanism in the rat outer medulla. I. Formulation and base-case results

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.

Two modes for concentrating urine in rat inner medulla

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.

An efficient numerical method for distributed-loop models of the urine concentrating mechanism

In this study we describe an efficient numerical method, based on the semi-Lagrangian (SL) semi-implicit (SI) method and Newton's method, for obtaining steady-state (SS) solutions of equations arising in distributed-loop models of the urine concentrating mechanism. Dynamic formulations of these models contain large systems of coupled hyperbolic partial differential equations (PDEs). The SL method advances the solutions of these PDEs in time by integrating backward along flow trajectories, thus allowing large time steps while maintaining stability. The SI approach controls stiffness arising from transtubular transport terms by averaging these terms in time along flow trajectories. An approximate SS solution of a dynamic formulation obtained via the SLSI method can be used as an initial guess for a Newton-type solver, which rapidly converges to a highly accurate numerical approximation to the solution of the ordinary differential equations that arise in the corresponding SS model formulation. In general, it is difficult to specify a priori for a Newton-type solver an initial guess that falls within the radius of convergence; however, the initial guess generated by solving the dynamic formulation via the SLSI method can be made sufficiently close to the SS solution to avoid numerical instability. The combination of the SLSI method and the Newton-type solver generates stable and accurate solutions with substantially reduced computation times, when compared to previously applied dynamic methods.