Biomedical subjects

Harold E Layton

Publications and source records attributed to Harold E Layton.

12 recordsLinked to original sources

PubMed2006

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.

Algorithms↗

PubMed2005

Multistability in tubuloglomerular feedback and spectral complexity in spontaneously hypertensive rats.

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.

Animals↗

PubMed2005

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.

Animals↗

PubMed2005

A region-based mathematical model of the urine concentrating mechanism in the rat outer medulla. II. Parameter sensitivity and tubular inhomogeneity.

In a companion study (Layton AT and Layton HE. Am J Physiol Renal Physiol 289: F1346-F1366, 2005), a region-based mathematical model was formulated for the urine concentrating mechanism (UCM) in the outer medulla (OM) of the rat kidney. In the present study, we quantified the sensitivity of that model to several structural assumptions, including the degree of regionalization and the degree of inclusion of short descending limbs (SDLs) in the vascular bundles of the inner stripe (IS). Also, we quantified model sensitivity to several parameters that have not been well characterized in the experimental literature, including boundary conditions, short vasa recta distribution, and ascending vasa recta (AVR) solute permeabilities. These studies indicate that regionalization elevates the osmolality of the fluid delivered into the inner medulla via the collecting ducts; that model predictions are not significantly sensitive to boundary conditions; and that short vasa recta distribution and AVR permeabilities significantly impact concentrating capability. Moreover, we investigated, in the context of the UCM, the functional significance of several aspects of tubular segmentation and heterogeneity: SDL segments in the IS that are likely to be impermeable to water but highly permeable to urea; a prebend segment of SDLs that may be functionally like thick ascending limb (TAL); differing IS and outer stripe Na(+) active transport rates in TAL; and potential active urea secretion into the proximal straight tubules. Model calculations predict that these aspects of tubular of segmentation and heterogeneity generally enhance solute cycling or promote effective UCM function.

Animals↗

PubMed2004

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.

Animals↗

PubMed2004

Feedback-mediated dynamics in two coupled nephrons.

Previously, we developed a dynamic model for the tubuloglomerular feedback (TGF) system in a single, short-looped nephron of the mammalian kidney. In that model, a semi-linear hyperbolic partial differential equation was used to represent two fundamental processes of solute transport in the nephron's thick ascending limb (TAL): chloride advection by fluid flow along the TAL lumen and transepithelial chloride transport from the lumen to the interstitium. An empirical function and a time delay were used to relate glomerular filtration rate to the chloride concentration at the macula densa of the TAL. Analysis of the model equations indicated that stable limit-cycle oscillations (LCO) in nephron fluid flow and chloride concentration can emerge for sufficiently large feedback gain magnitude and time delay. In this study, the single-nephron model was extended to two nephrons, which were coupled through their filtration rates. Explicit analytical conditions were obtained for bifurcation loci corresponding to two special cases: (1) identical time delays but differing feedback gains, and (2) identical gains but differing delays. Similar to the case of a single nephron, our analysis indicates that stable LCO can emerge in coupled nephrons for sufficiently large gains and delays. However, these LCO may emerge at lower values of the feedback gain, relative to a single (i.e., uncoupled) nephron, or at shorter delays, provided the delays are sufficiently close. These results suggest that, in vivo, if two nephrons are sufficiently similar, then coupling will tend to increase the likelihood of LCO.

Animals↗

PubMed2003

Effect of sustained flow perturbations on stability and compensation of tubuloglomerular feedback.

A mathematical model previously formulated by us predicts that limit-cycle oscillations (LCO) in nephron flow are mediated by tubuloglomerular feedback (TGF) and that the LCO arise from a bifurcation that depends heavily on the feedback gain magnitude, gamma, and on its relationship to a theoretically determined critical value of gain, gammac. In this study, we used that model to show how sustained perturbations in proximal tubule flow, a common experimental maneuver, can initiate or terminate LCO by changing the values of gamma and gammac, thus changing the sign of gamma - gammac. This result may help explain experiments in which intratubular pressure oscillations were initiated by the sustained introduction or removal of fluid from the proximal tubule (Leyssac PP and Baumbach L. Acta Physiol Scand 117: 415-419, 1983). In addition, our model predicts that, for a range of TGF sensitivities, sustained perturbations that initiate or terminate LCO can yield substantial and abrupt changes in both distal NaCl delivery and NaCl delivery compensation, changes that may play an important role in the response to physiological challenge.

Adaptation, Physiological↗

PubMed2003

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.

Animals↗

PubMed2003

An inverse algorithm for a mathematical model of an avian urine concentrating mechanism.

A nonlinear optimization technique, in conjunction with a single-nephron, single-solute mathematical model of the quail urine concentrating mechanism, was used to estimate parameter sets that optimize a measure of concentrating mechanism efficiency, viz., the ratio of the free-water absorption rate to the total NaCl active transport rate. The optimization algorithm, which is independent of the numerical method used to solve the model equations, runs in a few minutes on a 1000 MHz desktop computer. The parameters varied were: tubular permeabilities to water and solute; maximum active solute transport rates of the ascending limb of Henle and the collecting duct (CD); length of the prebend enlargement (PBE) of the descending limb; fractional solute delivery to the CD; solute concentration of tubular fluid entering the CD at the cortico-medullary boundary; and rate of exponential CD population decrease along the medullary cone. Using a base-case parameter set and parameter bounds suggested by physiologic experiments, the optimization algorithm identified a maximum-efficiency set of parameter values that increased efficiency by 40% above base-case efficiency; a minimum-efficiency set reduced efficiency by about 41%. When maximum-efficiency parameter values were computed as medullary length varied over the physiologic range, the PBE was found to make up 88% of a short medullary cone but only 8% of a long medullary cone.

Algorithms↗

PubMed2003

A region-based model framework for the rat urine concentrating mechanism.

The highly structured organization of tubules and blood vessels in the outer medulla of the mammalian kidney is believed to result in preferential interactions among tubules and vessels; such interactions may promote solute cycling and enhance urine concentrating capability. In this study, we formulate a new model framework for the urine concentrating mechanism in the outer medulla of the rat kidney. The model simulates preferential interactions among tubules and vessels by representing two concentric regions and by specifying the fractions of tubules and vessels assigned to each of the regions. The model equations are based on standard expressions for transmural transport and on solute and water conservation. Model equations, which are derived in dynamic form, are solved to obtain steady-state solutions by means of a stable and efficient numerical method, based on the semi-Lagrangian semi-implicit method and on Newton's method. In this application, the computational cost scales as O(N2), where N is the number of spatial subintervals along the medulla. We present representative solutions and show that the method generates approximations that are second-order accurate in space and that exhibit mass conservation.

Animals↗

PubMed2003

Advective transport of nitric oxide in a mathematical model of the afferent arteriole.

Endothelium-derived nitric oxide (NO) is thought to be short-lived in blood because of rapid removal from plasma, mainly by binding to Hb. The extent to which removal limits NO advection is unclear, especially for blood flow in the renal afferent arteriole (AA), which has a transit time of 3-30 ms. A mathematical model of AA fluid dynamics and myogenic response that includes NO diffusion, advection, degradation, and vasorelaxant action was used to estimate NO advective transport. Model simulations indicate that advective transport of locally produced NO is sufficient to yield physiologically significant NO concentrations along much of the AA. Advective transport is insensitive to NO scavenging by Hb because the NO-Hb binding rate is slow relative to AA transit time. Hence, plasma NO concentration near the vessel wall is influenced by both diffusion from endothelial cells and advection from upstream sites. Simulations also suggest that NO advection may constitute a mechanism to stabilize arteriolar flow in response to a localized vasoconstriction accompanied by enhanced NO release.

Animals↗

PubMed2002

A numerical method for renal models that represent tubules with abrupt changes in membrane properties.

The urine concentrating mechanism of mammals and birds depends on a counterflow configuration of thousands of nearly parallel tubules in the medulla of the kidney. Along the course of a renal tubule, cell type may change abruptly, resulting in abrupt changes in the physical characteristics and transmural transport properties of the tubule. A mathematical model that faithfully represents these abrupt changes will have jump discontinuities in model parameters. Without proper treatment, such discontinuities may cause unrealistic transmural fluxes and introduce suboptimal spatial convergence in the numerical solution to the model equations. In this study, we show how to treat discontinuous parameters in the context of a previously developed numerical method that is based on the semi-Lagrangian semi-implicit method and Newton's method. The numerical solutions have physically plausible fluxes at the discontinuities and the solutions converge at second order, as is appropriate for the method.

Animals↗