Transport efficiency and workload distribution in a mathematical model of the thick ascending limb.

The thick ascending limb (TAL) is a major NaCl reabsorbing site in the nephron. Efficient reabsorption along that segment is thought to be a consequence of the establishment of a strong transepithelial potential that drives paracellular Na(+) uptake. We used a multicell mathematical model of the TAL to estimate the efficiency of Na(+) transport along the TAL and to examine factors that determine transport efficiency, given the condition that TAL outflow must be adequately dilute. The TAL model consists of a series of epithelial cell models that represent all major solutes and transport pathways. Model equations describe luminal flows, based on mass conservation and electroneutrality constraints. Empirical descriptions of cell volume regulation (CVR) and pH control were implemented, together with the tubuloglomerular feedback (TGF) system. Transport efficiency was calculated as the ratio of total net Na(+) transport (i.e., paracellular and transcellular transport) to transcellular Na(+) transport. Model predictions suggest that 1) the transepithelial Na(+) concentration gradient is a major determinant of transport efficiency; 2) CVR in individual cells influences the distribution of net Na(+) transport along the TAL; 3) CVR responses in conjunction with TGF maintain luminal Na(+) concentration well above static head levels in the cortical TAL, thereby preventing large decreases in transport efficiency; and 4) under the condition that the distribution of Na(+) transport along the TAL is quasi-uniform, the tubular fluid axial Cl(-) concentration gradient near the macula densa is sufficiently steep to yield a TGF gain consistent with experimental data.

Fluid dilution and efficiency of Na(+) transport in a mathematical model of a thick ascending limb cell.

Thick ascending limb (TAL) cells are capable of reducing tubular fluid Na(+) concentration to as low as ~25 mM, and yet they are thought to transport Na(+) efficiently owing to passive paracellular Na(+) absorption. Transport efficiency in the TAL is of particular importance in the outer medulla where O(2) availability is limited by low blood flow. We used a mathematical model of a TAL cell to estimate the efficiency of Na(+) transport and to examine how tubular dilution and cell volume regulation influence transport efficiency. The TAL cell model represents 13 major solutes and the associated transporters and channels; model equations are based on mass conservation and electroneutrality constraints. We analyzed TAL transport in cells with conditions relevant to the inner stripe of the outer medulla, the cortico-medullary junction, and the distal cortical TAL. At each location Na(+) transport efficiency was computed as functions of changes in luminal NaCl concentration ([NaCl]), [K(+)], [NH(4)(+)], junctional Na(+) permeability, and apical K(+) permeability. Na(+) transport efficiency was calculated as the ratio of total net Na(+) transport to transcellular Na(+) transport. Transport efficiency is predicted to be highest at the cortico-medullary boundary where the transepithelial Na(+) gradient is the smallest. Transport efficiency is lowest in the cortex where luminal [NaCl] approaches static head.

Parameter estimation for mathematical models of NKCC2 cotransporter isoforms.

An optimization problem, formulated using a nonlinear least-squares approach, was used to estimate parameters for kinetic models of the three isoforms of the kidney-specific Na-K-2Cl (NKCC2) cotransporter. Specifically, the optimization problem estimates the magnitude of model parameters (i.e., off-binding and translocation rate constants) by minimizing the distance between model unidirectional fluxes and published unidirectional (86)Rb(+) uptake curves for the A, B, and F isoforms of the NKCC2 cotransporter obtained in transfected Xenopus oocytes. By using different symmetry assumptions, NKCC2 models with five, six, seven, or eight parameters were evaluated. The optimization method identified parameter sets that yielded computed unidirectional fluxes consistent with the uptake data. However, the parameter values were not unique, in that systematic exploration of the parameter space revealed alternative parameter sets that fit the data with similar accuracy. Finally, we demonstrate that the optimization method can identify parameter sets for the three transporter isoforms that differ only in ion binding affinities, a result that is consistent with a published mutagenesis analysis of the molecular and structural bases for the differences in (86)Rb(+) uptake among the A, B, and F isoforms. These NKCC2 cotransporter models will facilitate the development of larger scale models of ion transport by thick ascending limb cells.