A numerical model of the renal distal tubule
Hangil
Chang1 and
Toshiro
Fujita2
1 Health Service Center and
2 Fourth Department of Internal Medicine,
University of Tokyo, Tokyo 153-8902, Japan
 |
ABSTRACT |
A numerical model of the rat
distal tubule was developed to simulate water and solute transport in
this nephron segment. This model incorporates the following: 1)
Na-Cl cotransporter, K-Cl cotransporter, Na channel, K channel, and Cl
channel in the luminal membrane; 2) Na-K-ATPase, K channel, and
Cl channel in the basolateral membrane; and 3) conductances for
Na, K, and Cl in the paracellular pathway. Transport rates were
calculated using kinetic equations. Axial heterogeneity was represented
by partitioning the model into two subsegments with different sets of
model parameters. Model equations derived from the principles of mass
conservation and electrical neutrality were solved numerically. Values
of the model parameters were adjusted to minimize a penalty function that was devised to quantify the difference between model
predictions and experimental results. The developed model could
simulate the water and solute transport of the distal tubule in the
normal state, as well as in conditions including thiazide or amiloride application and various levels of sodium load and tubular flow rate.
transporters; ion channels; water and mineral metabolism
 |
INTRODUCTION |
DISTAL TUBULE is defined as the portion of the nephron
between the macula densa region and the first confluence with another tubule. This part of the nephron is the site where most of potassium secretion takes place, and normally as much as 70% of potassium in the
final urine is derived from this segment (44). Furthermore, the distal
tubule reabsorbs ~10% of filtered sodium (45) and is highly
permeable to water in the presence of antidiuretic hormone (11),
concentrating the initially hyposmotic tubular fluid to isosmolality.
Transport properties of the distal tubule have been studied extensively
in the rat using the micropuncture and in vivo microperfusion techniques. Those studies have revealed that 1) the early part (which corresponds to distal convoluted tubule) and the late part (which is comprised of connecting tubule and initial collecting tubule)
have different morphological and physiological characteristics; 2) the early part reabsorbs sodium mostly through the
thiazide-sensitive Na-Cl cotransporter in the luminal cell membrane and
has low permeability to water; 3) the late part is rich in
amiloride-sensitive Na channel, has highly negative luminal electrical
potential that drives potassium secretion, and is highly permeable to
water if antidiuretic hormone is present; 4) total sodium
reabsorption is load dependent; and 5) total potassium
secretion is flow dependent.
Recently, an increasing number of molecules involved in
the transport of water and solutes in the distal tubule have been identified by cDNA cloning. These molecules include amiloride-sensitive Na channel (8, 9, 41), thiazide-sensitive Na-Cl cotransporter (22),
ROMK K channel (32, 70), and water channel (21). Such progress in
molecular biology promotes detailed analysis of transport kinetics by
providing sufficient amounts of these molecules for flux measurements.
Accordingly, information regarding their transport mechanisms is
growing. One of the most effective usages of these data would be to use
them to predict the macroscopic transport in the renal tubule with the
help of numerical models (65). Numerical models suitable for this
purpose had been successfully constructed for proximal tubule (63, 64),
thin descending limb of Henle (46, 57), thick ascending limb of Henle
(20), and cortical collecting tubule (55). However, to our knowledge, models to simulate distal tubule have not been developed. In this study, we report a numerical model that simulates water and solute transport in the distal tubule of the rat. This model was constructed in the theoretical framework similar to the one adopted in the previous
models (55, 64) except that electrodiffusive mass transfer in the
tubular lumen was taken into account to accurately predict the profile
of luminal electrical potential. With appropriate values of the model
parameters, the model could simulate the transport properties of the
distal tubule in the basic state, as well as in experimental conditions
such as thiazide or amiloride application and various levels of sodium
load and tubular flow rate.
Glossary
R |
Gas constant,
J · mmol 1 · K 1 or
mmHg · cm3 · mmol 1 · K 1
|
F |
Faraday constant, C/mmol
|
R1 |
Outer radius of the distal tubule, cm
|
R2 |
Inner radius of the distal tubule, cm
|
x |
Axial coordinate along the tubule, cm
|
TL |
Total length of the distal tubule, cm
|
N |
Number of sections
|
w |
Width of each section, cm (w = TL/N )
|
m |
Superscript used to denote luminal compartment
|
c |
Superscript used to denote cellular compartment
|
s |
Superscript used to denote serosal compartment
|
Imp |
Impermeant solute
|
D |
Diffusion coefficient of solute , cm2/s
|
C pk |
Concentration of solute k in compartment p,
mmol/cm3
|
p |
Electrical potential of compartment p, V
|
Fv |
Tubular flow rate, cm3/s
|
J pqv |
Volume flux from compartment p to compartment q,
cm/s (All the transepithelial fluxes were expressed as transport
rate per tubular outer surface area.)
|
Jv |
Transepithelial volume flux, cm/s
(Jv = J mcv +
J msv)
|
J pqk |
Flux of solute k from compartment p to compartment
q,
mmol · cm 2 · s 1
|
Jk |
Transepithelial flux of solute k,
mmol · cm 2 · s 1
(Jk = J mck +
J msk)
|
zk |
Valence of solute k
|
zpImp |
Valence of impermeant solute in compartment p
|
|
Discretized model variables are coded by supplying an additional
superscript indicating the section number. For example,
Cc,iNa stands for the
concentration of intracellular sodium in the i-th section.
Superscripts 0 and N + 1 indicate the entrance and exit of the
distal tubule, respectively. Thus, Cm,0K and Cm,N + 1K stand for
potassium concentrations of fluid entering and leaving the distal
tubule, respectively.
|
 |
METHODS |
Overview of the Model
Distal tubule was modeled as a rigid cylinder made of a single layer of
epithelial cells (Fig. 1A). Three
compartments were discriminated: luminal, cellular, and serosal
compartments (Fig. 1B). The luminal compartment is the space
in which tubular fluid flows. The cellular compartment is the cytosolic
space of epithelial cells. The serosal compartment is the well-stirred
bath surrounding the tubule. Transports of water, sodium, potassium,
and chloride among these compartments were considered. Transport
between the luminal and cellular compartments is governed by the Na-Cl
cotransporter, K-Cl cotransporter, Na channel, K channel, and Cl
channel (Fig. 2, A and
B). Transport between the
cellular and serosal compartments is governed by Na-K-ATPase, K
channel, and Cl channel. Paracellular pathway represents the (direct)
transport between the luminal and serosal compartments, and it is
comprised of conductive pathways for sodium, potassium, and chloride.
Solute fluxes through ion channels were calculated using the
Goldman-Hodgkin-Katz current equation (26, 31, 33). Solute fluxes
through Na-K-ATPase (69), Na-Cl cotransporter, and K-Cl cotransporter
were calculated using appropriate kinetic models. We assumed water flux
was proportional to the osmotic difference across the membrane.

View larger version (13K):
[in this window]
[in a new window]
|
Fig. 1.
Schematic diagram of model tubule. A: outline of the model.
Distal tubule was modeled as a cylinder made of a sheet of epithelial
cells. In the lumen, tubular fluid flows from the proximal end to the
distal end as indicated by arrows. In the derivation of discretized
model equations, tubule was divided into N sections of equal
width. B: magnified view. There are three distinct
compartments; luminal, cellular, and serosal. Solute concentrations and
electrical potential within each compartment as well as tubular flow
rate in the luminal compartment are taken into account in the model.
|
|


View larger version (26K):
[in this window]
[in a new window]
|
Fig. 2.
Transport mechanisms of model tubule. In the luminal cell membrane,
there are Na-Cl cotransporter, K-Cl cotransporter, Na channel, K
channel, and Cl channel. In the basolateral cell membrane, there are
Na-K-ATPase, K channel, and Cl channel. In the paracellular pathway,
which faces luminal and basolateral compartments, there are
conductances for sodium, potassium, and chloride. Axial heterogeneity
of the distal tubule was represented by changing the model parameters
in the early and the late parts of the model tubule. Transport rates
(in units of pmol/min) in the basic state through individual transport
devices in the early (A) or late (B) parts of the
distal tubule are shown by numbers beside arrows that indicate the
direction of transport.
|
|
A system of differential equations was postulated on the basis of
principles of mass conservation and electrical neutrality. To transform
these differential equations to difference equations suitable for
numerical solution, the model tubule was divided into N
sections with equal width (Fig. 1A). Then, the continuous model variables (solute concentrations, electrical potential, and
tubular flow rate) in the differential equations were replaced by
discrete model variables approximating the continuous counterparts at
the center of each section. Two sets of difference equations were
derived depending on the treatment of axial electrodiffusive mass
transfer in the luminal space. In the first set, electrodiffusive terms
were neglected as in the previous models (55, 64), and the derived
equations could be solved efficiently in terms of computational time
and computer memory. In the second set, electrodiffusive terms were
taken into account. When each solution of these sets of equations were
compared, we found that electrodiffusive terms have to be
taken into account to accurately estimate the profile of luminal
electrical potential in the distal tubule, which is comprised of
heterogeneous subsegments.
A computer program that was developed for this project solved
difference equations by the Newton-Raphson method (50) and determined
the transepithelial fluxes of water and solutes as well as profiles of
solute concentrations and electrical potential in the luminal and
cellular compartments. Values of the model parameters were adjusted to
minimize a penalty function that was devised to quantify the difference
between model predictions and experimental measurements. The definition
of the penalty function reflected the results of free-flow
micropuncture experiments of the rat. With the use of Powell's
direction set method in multidimensions (50), we could find a set of
model parameters that minimizes the penalty function.
Model Variables
Model variables are comprised of solute concentrations, electrical
potential and tubular fluid flow rate in the luminal compartment; and
solute concentrations and electrical potential in the cellular compartment. In the continuous form, they are
CmNa, CmK,
CmCl, CmImp,
m, CcNa,
CcK, CcCl,
CcImp,
c, and Fv.
In the discretized form, they are
Cm,iNa (i = 1, ... , N + 1),
Cm,iK
(i = 1, ... , N + 1),
Cm,iCl
(i = 1, ... , N + 1),
Cm,iImp (i = 1, ... , N + 1),
m,i (i = 0, ... ,
N + 1), Cc,iNa
(i = 0, ... , N + 1),
Cc,iK
(i = 0, ... , N + 1),
Cc,iCl (i = 0, ... , N + 1),
Cc,iImp
(i = 0, ... , N + 1),
c,i (i = 0, ... ,
N + 1), and Fiv
(i = 1, ... , N + 1). Total number of model
variables in the discretized form amounts to 11N + 17. Note
that Cm,0Na, Cm,0K, Cm,0Cl, Cm,0Imp,
and F0v (luminal solute concentrations and
flow rate at the entrance of the tubule) are given as boundary
conditions and are not model variables.
Model Equations
Equations describing the steady state of the model are listed in Table
1. Equations 1 to 4
represent conservation of water and solutes in the cellular
compartment. Equation 5 represents electrical neutrality
within the cellular compartment. Equations 6 to
ME-10 represent conservation of water and solutes in the luminal compartment. Equation 11 implies electrical
neutrality within the mucosal compartment.
All these equations could be expressed in terms of the model variables
using the equations listed in Table 2. Flux
of water across the membrane was assumed to be proportional to the
difference in osmolality across the membrane
(Osmq
Osmp) (Eq. 1). We have assumed reflection coefficients of all the solutes
to be 1.0, and we ignored the difference between osmolarity and
osmolality (Eq. 2). Equation 3 is the
Goldman-Hodgkin-Katz current equation (26, 31, 33) which was applied to
all the conductive pathways (ion channels and paracellular pathways) in the present study. Equation 4 is a kinetic equation used for
the electroneutral cation and anion cotransporter with the
stoichiometry of 1:1 (Na-Cl cotransporters and K-Cl cotransporters).
Derivation of Eq. 4 is described in APPENDIX 1. Rate of ATP hydrolysis by Na-K-ATPase is assumed to obey Eq. 5 (69). As one molecule of ATP is hydrolyzed, three molecules of
sodium are transported out of the cell and two molecules of potassium
are transported into the cell.
Across the luminal cell membrane, sodium is transported through Na-Cl
cotransporter and Na channel (Eq. 6); potassium through K-Cl
cotransporter and K channel (Eq. 7); and chloride through Na-Cl cotransporter, K-Cl cotransporter, and Cl channel (Eq. 8). Across the basolateral cell membrane, sodium is transported
through Na-K-ATPase (Eq. 9); potassium through Na-K-ATPase
and K channel (Eq. 10); and chloride through Cl channel
(Eq. 11). Across the paracellular pathway, sodium, potassium
and chloride are transported through selective conductances (Eqs.
12 to 14). Transepithelial transport is comprised of
transcellular transport and paracellular transport (Eqs. 15
to 18). Axial flow of solute
in the luminal compartment is the sum of the convective term
(FvCm
) and the electrodiffusive term (Eq. 19).
Discretization and Solution of the Model Equations Without
Electrodiffusive Terms
Without the electrodiffusive terms, axial flow of solute
(Eq. 19) is reduced to
If Eq. 11 is taken into account, then it follows
that axial electrical current flow is constantly zero (and its
derivative is also zero) along the tubule
Therefore, transepithelial electrical current is also zero,
which can be deduced as follows. From Eqs. 7 to
10
and because
d(
z
F
)/dx = 0 from the above consideration
|
(ME-11`)
|
We substituted Eq. 11` for Eq. 11 in
the following derivation of difference equations. Equation 11` implies that, in the model without electrodiffusion, luminal
electrical potential is determined as the potential invoking no (local)
transepithelial electrical current.
Discretized model equations are listed in Table
3. Applying Eqs. 1 to
5 and 11` to discretized model variables at the
entrance of the tubule (i = 0), Eqs. 1 to
6 are derived. Solving these equations by Newton-Raphson
method, Cc,0Na,
Cc,0K, Cc,0Cl,
Cc,0Imp,
c,0, and
m,0 can be determined. Initial values for the iteration
were as follows (6): CcNa, 12 mM;
CcK, 170 mM; CcCl, 15 mM; CcImp, 108 mM; and
c,
90 mV.
Equations 7 to 17 are derived from Eqs.
1 to 10 and 11`. By sequentially solving
these equations with i = 1 to N + 1, values for all the
discretized model variables can be determined.
Discretization and Solution of the Model Equations With
Electrodiffusive Terms
In the derivation of the difference equations from the model equations
taking into account electrodiffusive terms, additional (discretized)
model variables estimating the derivative of luminal solute
concentrations (dCi
/dx)
and luminal electrical potential (d
i/dx)
have been introduced. Discretized model equations are listed in Table
4. Equations 1 to 5
and Eq. 11 are from Eqs. 1 to 5 and 11. Equations 6 to 10 are difference
equations derived from Eqs. 6 to 10. In these
equations, Fi
can be expressed
in terms of model variables using Eq. 17, which corresponds to Eq. 19 in Table 2. Relationships between
Cm,i
and
dCi
/dx are as in
Eq. DE2-12. Similar relations also apply to
m,i and d
i/dx
(Eq. DE2-13). Equations 14 and 16
represent the boundary conditons. In micropuncture experiments, tubular
fluid leaving the distal tubule is usually collected into a collecting
pipette that can be regarded as electrically insulated. We assumed that
in the steady state, solute concentrations and electrical potential are homogeneous in the collecting pipette including its orifice (Eqs. 14 and 15). We also assumed that the system
providing the tubular fluid to the distal tubule (upstream nephron
segment or perfusion pipette) can be regarded, as a first
approximation, as being electrically insulated and that there is no net
axial electrical current at the beginning of the distal tubule (Eq. 16). This system of equations was solved with the
Newton-Raphson method. In those calculations, solution of the equations
without electrodiffusive terms was used as the initial values for
iteration.
Constants
Constants used in the present study are listed in Table
5. Temperature was assumed to be 37°C.
R1 and R2 were from de
Mello-Aires et al. (15). Total length of the model distal tubule
(TL) was from Good and Wright (28). Distal tubule was divided
into two parts: early and late distal tubules. Lengths of the early and late distal tubules were 0.1 cm (42) and 0.13 cm, respectively. Number
of sections (N ) was chosen to be 23, and width of each section (w) was 0.01 cm. Diffusion coefficients of solutes at infinite dilution were from Hille (31) and Cussler (14). (Diffusion coefficient of impermeant solute was that of urea.) Values for serosal
sodium and potassium concentrations were from Malnic et al. (45), and
those for chloride were from Kunau et al. (39). Concentration of
impermeant solute in the serosal compartment was calculated from the
reported osmolality of the rat plasma (305.3 mosmol/kg) (59)
Because most of the impermeant solutes in the mucosal
compartment can be assumed to be urea, the valence
(zmImp) was considered to be
0. Valence of impermeant solute in the serosal compartment was
calculated from solute concentrations based on the principle of
electrical neutrality.
Penalty Function
Transport characteristics of the model tubule are completely specified
by a set of model parameters. Ideally, values of these parameters
should be determined by experimental measurements. However, since
direct evaluations are not available for most of the parameters, we
attempted to find the values with which the model simulates the real
distal tubule best. Specifically, we devised a penalty function that
gives larger values as the model prediction deviates from
experimental results and systematically searched for the parameter
values which minimize this function.
In this report, we confined our studies to use the results
of free-flow micropuncture experiments in the rat, since experimental conditions were relatively uniform in those studies. We summarized their results as follows. 1) In the basic state, the distal
tubule reabsorbs water, sodium, and chloride; secretes potassium; and has negative luminal electrical potential. 2) The distal tubule has thiazide-sensitive Na-Cl cotransporters. Application of thiazide reduces sodium reabsorption but has a negligible effect on
transepithelial potential difference. 3) The distal tubule has
amiloride-sensitive sodium channels. Application of amiloride
depolarizes luminal electrical potential but has a minimal effect on
total sodium reabsorption. 4) Sodium reabsorption is load
dependent. As the rate of sodium entry increases, sodium reabsorption
increases. 5) Potassium secretion is flow dependent. As the
initial tubular flow rate increases, potassium secretion increases.
Definition of the penalty function incorporating these characteristics
is presented in Table 6. We chose 13 representative experimental measurements (E1 to
E13). E1 to
E5 were from experiments conducted in the basic
state. Total water reabsorption rate (E1) and
total sodium reabsorption rate (E2) were from
Malnic et al. (45). Total potassium secretion rate
(E3) was from Malnic et al. (44). Luminal
electrical potential in the late distal tubule (E5) was from Hayslett et al. (30). Luminal
electrical potential in the early distal tubule
(E4) was based on the experimental observation
that luminal potential in the early distal tubule was depolarized by
23.3 mV relative to that in the late distal tubule (5). Model
predictions corresponding to these values (S1 to
S5) were calculated with the boundary condition
listed in Table 6. Values of this boundary condition were determined from reports of free-flow micropuncture experiments (44, 45). Simulated
values of luminal potential in the early (S4) and
late (S5) distal tubules were the values of
luminal potential at the midpositions of these subsegments
(
m,5 and
m,16).
E6 to E8 are from experiments
conducted in the presence of thiazide. E6
represents total sodium reabsorption rate. E7 and E8 represent luminal potential in the early and
late distal tubules, respectively. These values were based on the
results of Costanzo and Windhager (13), which showed that thiazide
reduced sodium reabsorption by 43% compared with the basic state and
did not affect luminal electrical potential. Model predictions
(S6 to S8) corresponding to
E6 to E8 were calculated
similarly as in the basic state, except that
J mcNa,max (both in the early
and late distal tubules) was reduced by 90%.
E9 to E11 are related to the
distal tubular function when amiloride was applied.
E9 represents total sodium reabsorption rate. This
value was set equal to that of the basic state
(E2) based on the report of Duarte et al. (16) in
which amiloride had essentially no effect on total sodium reabsorption.
E10 and E11 are luminal potential in the early and late distal tubules, respectively, and the values for these were based on the report by
Barratt (4) that amiloride depolarized luminal potential by 1.5 and
20.5 mV in the early and late distal tubules, respectively. Model
predictions (S9 to S11)
corresponding to E9 to E11 were
calculated with the same boundary conditions as in the basic state,
except that P mcNa (both in the
early and late distal tubules) was reduced by 99%.
E12 is from the measurement (37) of total sodium
reabsorption rate when sodium load was increased to 2,160 pmol/min.
Model prediction (S12) corresponding to
E12 was calculated using the boundary condition
that simulates the one observed when the experiment was conducted (37).
E13 is from the measurement of total potassium
secretion rate when tubular flow rate was experimentally increased to
38.1 nl/min (38). S12 is the simulated value
obtained with the boundary condition resembling the one observed in the
experiment (38).
Parameter Search
There were 40 model parameters in total. The minimization algorithm
used in this study (Powell's direction set method) is an iterative
procedure and requires a set of initial (presumptive) values for the
parameters to start with. These initial parameters, listed in Table
7, were deduced from (manual) calculations
that were based on the information relevant to the basic state (Tables 6B and 12). Apical hydraulic conductivity
(Lmcv) in the late distal
tubule had a value (7.48 × 10
8
cm · s
1 · mmHg
1)
similar to those reported for transepithelial hydraulic conductivity of
the nondiuretic rat (6.47 × 10
8
cm · s
1 · mmHg
1,
Ref. 54; and 10.03 × 10
8
cm · s
1 · mmHg
1,
Ref. 58). Intrinsic dissociation constants of Na-Cl cotransporter (KmcNa,NaCl and
KmcCl,NaCl) were selected to
account for the measured kinetic constants of this transporter. With
the listed values, the kinetic equation (Eq. 4 in Table 2)
exhibits effective affinities (53) of 25.0 and 13.6 mM for sodium and
chloride, respectively. These values are identical to the ones reported
with the thiazide-sensitive Na-Cl cotransporter expressed in
Xenopus oocytes (22). Intrinsic dissociation constants for K-Cl
cotransporter (KmcK,KCl and
KmcCl,KCl) were also selected
to account for the kinetic constants reported for the recently cloned
K-Cl cotransporter (25) that is expressed in the kidney. The kinetic equation predicts effective affinities of 25.0 and 50.0 mM for potassium and chloride, respectively. These values are identical to the
experimental results (25). Maximal rate of ATP hydrolysis (Ja,max) was 317 pmol · min
1 · tube
1,
which was similar to the experimental value (34) for dissected rat
distal tubule (257 pmol · min
1 · tube
1).
Value of kinetic parameter of Na-K-ATPase
(KNa,ATPase) was that of Na-K-ATPase
with rat
1-isoform, which is the predominant isoform in the kidney,
reported in a recent study (69). Ion permeabilities through the
paracellular pathway (PmsNa, PmsK, and
PmsCl) were calculated from a
measurement of transepithelial conductance (2.6 mS/cm2)
in the early distal tubule (43). With the listed values, paracellular pathway of the model has conductance of 2.5 mS/cm2.
To calculate initial values for other parameters, we made following
presumptions: cytosolic solute concentrations are close to the values
reported by Beck et al. (6); luminal potassium is nearly at equilibrium
with cytosolic potassium at the end of the distal tubule; the early
distal tubule accounts for 80% of total sodium reabsorption; potassium
transport rate through K channel is equal to chloride transport rate
through Cl channel in the basolateral membrane; potassium secretion
rate in the early distal tubule is small (1 pmol/min); apical K-Cl
cotransporter and K channel contribute equally to potassium secretion
both in the early and late distal tubules; apical Na-Cl cotransporter accounts for 90% of the sodium reabsorption in the early distal tubule; water reabsorption rate in the early distal tubule is 1 nl/min;
hydraulic conductivities in the basolateral membrane (Lscv) are 100-fold larger
than that in the luminal membrane
(Lmcv); and water flux in the
paracellular pathway is negligible.
With the initial parameter values in Table 7, the penalty function was
33.8 (calculation with electrodiffusive terms). Starting with these
values, preliminary attempts to adjust all the model parameters to
minimize the penalty function were unsuccessful. The minimization
procedure required too much computational time (more than 3 days), and
the acquired solutions were mostly unrealistic in terms of predicted
solute concentrations and electrical potential in the cytosolic space.
Hence, we restricted the parameters for adjustment to the ones
representing the maximum transport rates of the luminal transporters
(J mcNaCl,max, PmcNa, and
PmcCl in the early distal
tubule; and J mcNaCl,max, J mcKCl,max,
PmcNa, PmcK, and
PmcCl in the late distal
tubule). These parameters are underlined in Table 7. Exclusion of the
basolateral parameters ensured the acceptable predictions for cytosolic
variables. J mcKCl,max and
PmcK in the early distal tubule
were also excluded from optimization, since it was obvious from
experiments (52) that they should remain very small. Furthermore, one
constraint was imposed in the minimization procedure:
J mcNaCl,max in the early
distal tubule is 10-fold larger than that in the late distal tubule.
Without this constraint, the search procedure converged to a much
larger value of J mcNaCl,max in
the late distal tubule than that of the early distal tubule, which is
contradictory to experimental results (12, 19). Accordingly, among the
eight parameters adjusted, seven of these were independently varied.
Optimized values of these parameters are listed in Table 8, with which the penalty function was
0.327. Thus minimization procedure reduced the value of the penalty
function by a factor of 103. Other parameters are also included in
Table 8 to present the complete set of model parameters used throughout
this study.
Computer Software
A computer program for this project was implemented with C++
programming language (56). Discretized model equations were solved by
Newton-Raphson method in multidimensions (50). Systems of linear
equations were solved by LU decomposition (50).
Minimization of the penalty function was done using direction set
(Powell's) method in multidimensions (50), which do not require
derivatives of the target function. All the calculations were conducted
in 64-bit double precision. The program was compiled by Metrowerks C/C++ compiler (version 1.8) and run on a personal computer with 180 MHz PowerPC 604e CPU running under MacOS.
 |
RESULTS |
Comparison Between Calculations With and Without Electrodiffusive
Terms
To solve the model equations with our computer program, it took ~0.4
and 14 s without and with electrodiffusive terms, respectively. This
difference consequently affected the time required to calculate the
penalty function, and it took 4 and 73 s. Minimization of the penalty
function completed in ~30 min or 12 h depending on whether we neglect
those terms. Necessary heap memory was 0.15 megabytes without
electrodiffusion and 1.7 megabytes with electrodiffusion. Thus,
calculations without electrodiffusive terms could be conducted much
more efficiently.
When there is no electromotive force, we can safely drop the
electrodiffusive terms in modeling the renal tubule (40). To examine
whether similar simplification was acceptable in our case, solutions of
the model equations with or without electrodiffusion were compared.
Profiles of luminal variables of both models are plotted in Fig.
3, and profiles of transepithelal water and
solute fluxes are plotted in Fig. 4. As can
be seen in Fig. 3, B-D, predictions of the two
models were similar with regard to luminal flow rate and solute
concentrations. This was also the case for water and solute fluxes
(Fig. 4) and cytosolic variables (data not shown), although a small
difference in the transepithelial potassium flux was observed in the
late distal tubule (Fig. 4C). However, there was a
significant difference in the luminal potential profile (Fig.
3A) between these models. The model without electrodiffusion (broken line, Fig. 3A) predicted an
unrealistically abrupt change in luminal potential at the junction of
the early and late distal tubules (x = 0.10 cm), in contrast
to the smooth transition discernible in the model with electrodiffusive
terms (solid line, Fig.
3A).1 This was not
unexpected, since in the model without electrodiffusion the luminal
potential is determined as the potential with which there is no local
transepithelial current (note Eq. 11`) and the levels of
luminal potential satisfying this condition are quite different
between the two subsegments. Thus, in modeling the distal tubule,
simplification by dropping the electrodiffusive terms was considered to
be unacceptable. In the following part of this report, only the
results with electrodiffusive terms will be presented.

View larger version (18K):
[in this window]
[in a new window]
|
Fig. 3.
Comparison of profiles of luminal variables predicted by models with
(solid lines) and without (broken lines) electrodiffusion. Model
parameters listed in Table 8 were used. A: luminal electrical
potential. B: tubular flow rate. C: luminal sodium
concentration. D: luminal potassium concentration. Note the
abrupt change of electrical potential at the junction of the early and
late distal tubules (x = 0.10 cm) in the solution without
electrodiffusive terms.
|
|

View larger version (10K):
[in this window]
[in a new window]
|
Fig. 4.
Comparison of profiles of water and solute fluxes predicted by models
with (solid lines) and without (broken lines) electrodiffusion.
Calculations were conducted in a manner similar to that in Fig. 3.
A: transepithelial water flux. B: transepithelial
sodium flux. C: transepithelial potassium flux.
|
|
Basic State
The model could simulate total water and solute reabsorptions as well
as luminal potential of the distal tubule in the basic state (Table
9). More detailed description of the
magnitude of transport rate through individual transport devices can be
found in Table 10. Profiles of luminal
variables and transepithelial fluxes are plotted as continuous lines in
Figs. 3 and 4. Most (80%) of water reabsorption took place in the late
distal tubule (Table 10; Fig. 6A) due to the 14-fold higher
water permeability in the late subsegment of the model tubule than that
in the early subsegment. This is consistent with the report by Woodhall
and Tisher (66), which showed that the late subsegment is the primary site of water reabsorption in the distal tubule. Initially hyposmotic (184 mosmol/kg) tubular fluid (Table 6B) was concentrated to near isomolality at the proximal half of the late distal tubule. Total
sodium reabsorption was 389 pmol/min, of which 84% was in the early
distal tubule (Table 10; Fig. 6B). Sodium reabsorption in the
early subsegment occurred almost exclusively through the luminal Na-Cl
cotransporter, whereas both Na-Cl cotransporter and Na channel
contributed in the late subsegment. The low level of luminal sodium
concentration and hyperpolarized luminal potential (Fig.
5) resulted in substantial sodium secretion
through paracellular pathways (46 and 88 pmol/min in the early and late
distal tubules, respectively). Potassium was secreted mainly in the
late distal tubule (Table 10 and Fig.
6C), which is compatible with the
report by Stanton et al. (52) that potassium secretion in the early subsegment was negligible. Potassium secretion was mostly through luminal K channel in the present
model.2 Intratubular potassium
concentration progressively increased along the tubule and reached a
plateau at 14.1 mM (Fig. 5), which is close to the equilibrium
concentration of 13.8 mM. Values of cytosolic variables (Table 11) were
compatible with the experimental results by Beck et al. (6).

View larger version (20K):
[in this window]
[in a new window]
|
Fig. 5.
Luminal variables in presence of thiazide (solid lines) or in the basic
state (broken lines). A: luminal electrical potential.
B: tubular flow rate. C: luminal sodium concentration.
D: luminal potassium concentration.
|
|

View larger version (11K):
[in this window]
[in a new window]
|
Fig. 6.
Transepithelial water and solute fluxes in presence of thiazide (solid
lines) or in the basic state (broken lines). A: transepithelial
water flux. B: transepithelial sodium flux. C:
transepithelial potassium flux. Positive values indicate net
reabsorption.
|
|
Thiazide
The model could simulate the effects of thiazide on total reabsorption
of sodium and luminal potential of the distal tubule (Table 9). In the
simulation, thiazide reduced total sodium reabsorption by 49% and had
only small effects on luminal potential, which is consistent with
experimental results. To present the results of the simulation in
further detail, fluxes through individual transport mechanisms are
listed in Table 10, and profiles of luminal variables and
transepithelial fluxes are plotted in Figs. 5 and 6, respectively. As
can be seen in Table 10, the reduction in sodium reabsorption occurred
largely in the early distal tubule. This reflects the property of the
model distal tubule in which ~90% of thiazide-sensitive Na-Cl
cotransporters are in the early distal tubule. In the late distal
tubule, in addition to the relative sparsity of the luminal Na-Cl
cotransporter, compensatory increase in sodium reabsorption through the
luminal Na channel contributed to the smaller reduction in sodium
reabsorption. As a result of the inhibition of luminal Na-Cl
cotransporters, luminal concentrations of both sodium and chloride
tended to increase (Fig. 5), which should have contributed to the
thiazide-induced 140% increase in potassium secretion (Table 10) by
two mechanisms. In the first place, the consequent increase in
osmolality of the tubular fluid hindered water reabsorption and led to
higher flow rates throughout distal tubule. Indeed, tubular flow rate
at the end of the distal tubule increased by a factor of 1.5 in
response to thiazide. Since potassium secretion is a flow-dependent
process as shown below, the rise in tubular flow rate should bring
about an increase in potassium secretion. In the second place,
increased luminal sodium chloride concentration depolarized cytosolic
potential and hyperpolarized luminal potential in the late distal
tubule, increasing the electromotive force of potassium secretion. This
is not an intuitively evident result, since increases in luminal sodium
and chloride concentrations would have opposite effects. Increase in
luminal sodium would depolarize the apical cell membrane and
hyperpolarize luminal potential, and an increase in luminal chloride
would induce the opposite effects. Consequently, direction of changes
depends on other factors including relative magnitude of the
conductances of these ions. Simulation of thiazide application
predicted cytosolic depolarization (Table 11) and luminal
hyperpolarization. In summary, the main effects of thiazide in the
simulation were 49% reduction in total sodium reabsorption due to the
inhibition of Na-Cl cotransporter and 140% increase in total potassium
secretion due to increased tubular flow rate, depolarized cytosolic
potential in the late distal tubule, and hyperpolarized luminal potential.
Amiloride
The model could simulate the effects of amiloride on total sodium
reabsorption and luminal potential of the distal tubule (Table 9). In
the simulation, amiloride reduced total sodium reabsorption by only
5.4%, but markedly depolarized luminal potential. The latter effect
was more prominent in the late distal tubule, and at its end, luminal
potential depolarized by 16.1 mV in response to amiloride. To present
the results of the simulation in further detail, fluxes through
individual transport mechanisms are listed in Table 10, and profiles of
luminal variables and transepithelial fluxes are plotted in Figs.
7 and 8,
respectively. Although amiloride abolished sodium reabsorption through
luminal Na channels (Table 10), compensatory increase in sodium
reabsorption through luminal Na-Cl cotransporters and decrease in
sodium secretion through paracellular pathways prevented a large change
in total sodium reabsorption. Inhibition of (electrogenic) luminal Na
channel led to depolarization of luminal potential and
hyperpolarization of cytosolic potential, which caused a marked
inhibition of potassium secretion by reducing the electromotive force
driving potassium secretion. Specifically, luminal potential at the end
of distal tubule was
8.9 mV, whereas it was
25.1 mV in the basic
state, and average cytosolic potential in the late distal tubule
hyperpolarized by 1.6 mV in response to amiloride (Table 11). These
changes lowered equilibrium concentration of luminal potassium to 7.7 mM, which is close to the luminal potassium concentration at the end of the distal tubule (7.2 mM). In summary, the main effect of amiloride in
the simulation was depolarization of luminal potential, which led to a
marked reduction of tubular potassium secretion.

View larger version (19K):
[in this window]
[in a new window]
|
Fig. 7.
Luminal variables in presence of amiloride (solid lines) or in the
basic state (broken lines). A: luminal electrical potential.
B: tubular flow rate. C: luminal sodium concentration.
D: luminal potassium concentration.
|
|

View larger version (10K):
[in this window]
[in a new window]
|
Fig. 8.
Transepithelial water and solute fluxes in presence of amiloride (solid
lines) or in the basic state (broken lines). A: transepithelial
water flux. B: transepithelial sodium flux. C:
transepithelial potassium flux. Positive values indicate net
reabsorption.
|
|
Sodium Load
The model prediction of total sodium reabsorption under the condition
of high sodium load was comparable with experimental results (Table 9).
Details of the simulation are presented in Table 10 and Figs.
9 and 10.
Water reabsorption increased to 11.51 nl/min (Table 10). The late
distal tubule appears to contribute to most of the increase (Fig.
10A). To simulate high sodium load, we raised initial luminal
sodium concentration by 1.2-fold and tubular flow rate by 4.4-fold
(Table 6). At this level of flow rate, osmolality at the end of the
tubule no longer achieved isosmolality and was significantly hyposmotic
(256 mosmol/kg). Sodium reabsorption increased to 687 pmol/min (Table
10). Both luminal Na channel and Na-Cl cotransporter contributed to
this increase in sodium reabsorption, which was driven by higher
luminal sodium concentration in increased sodium load. Changes in other
factors including intracellular sodium concentration, luminal
potential, and cytosolic potential were in the opposite directions to
explain the increased sodium reabsorption. Potassium secretion was
markedly increased to 106.3 pmol/min. Hyperpolarized luminal potential,
as well as lower luminal potassium concentration, contributed to this
increase in potassium secretion. Figure
11 presents total sodium reabsorptions in
various levels of sodium load. The solid line in Fig. 11 indicates the model prediction, and the solid circles indicate experimental results
from Kunau et al. (37) and Costanzo and Windhager (13). The model
predicted the positive correlation between sodium load and sodium
reabsorption observed in experimental measurements.

View larger version (19K):
[in this window]
[in a new window]
|
Fig. 9.
Luminal variables with high sodium load (solid lines) or in the basic
state (broken lines). A: luminal electrical potential.
B: tubular flow rate. C: luminal sodium concentration.
D: luminal potassium concentration.
|
|

View larger version (10K):
[in this window]
[in a new window]
|
Fig. 10.
Transepithelial water and solute fluxes with high sodium load (solid
lines) or in the basic state (broken lines). A: transepithelial
water flux. B: transepithelial sodium flux. C:
transepithelial potassium flux. Positive values indicate net
reabsorption.
|
|

View larger version (12K):
[in this window]
[in a new window]
|
Fig. 11.
Relationship between sodium reabsorption and sodium load. Effect of
sodium load on total sodium reabsorption was examined by solving the
model equations with various amounts of sodium load. Sodium load was
manipulated by altering initial flow rate and tubular sodium
concentration by means of linear interpolation of the basic state and
the sodium-loaded state used in the definition of the penalty function
(Table 6). Specifically, Fv and
Cm,0Na were changed using the following
equations with x varying from 0.15 to 1.0: Fv = 7.99 + (35.02 7.99) x (in nl/min);
Cm,0Na = 51.06 + (61.74 51.06) x (in mM); where x is a dummy
parameter with x = 0 and x = 1 corresponding to the
basic state and the sodium-loaded state (as used in the penalty
function), respectively. Solid circles indicate experimental results
from Refs. 13 and 37.
|
|
High Flow Rate
The boundary condition of high flow rate was similar to that of high
sodium load (Table 6). Accordingly, simulation results in high flow
rate were similar to those in high sodium load, thus their details are
not presented here. Figure 12 presents
the relationship between potassium secretion and initial flow rate; the
solid line indicates the model prediction, and the solid circles
indicate the experimental results from Kunau et al. (38) and Reineck et
al. (51). The model prediction is in good agreement with the
experimental observation.

View larger version (11K):
[in this window]
[in a new window]
|
Fig. 12.
Relationship between potassium secretion and tubular flow rate. Effect
of tubular flow rate on total potassium secretion was examined by
solving the model equations with various amounts of initial flow rate.
Initial flow rate as well as initial tubular potassium concentration
were altered by means of interpolation of the basic state and the high
flow state used in the definition of the penalty function (Table 6).
Specifically, Fv and Cm,0K were
changed using the following equations with x varying from
0.15 to 1.0: Fv = 7.99 + (38.06 7.99) x
(in nl/min); Cm,0K = 2.30 + (1.86 2.30) x (in mM); where x is a dummy
parameter with x = 0 and x = 1 corresponding to the
basic state and the high flow state (as used in the penalty function),
respectively. Solid circles indicate experimental results from Refs. 38
and 51.
|
|
 |
DISCUSSION |
In summary, we developed a numerical model of the rat distal tubule
that could simulate the experimental results in various conditions
including the basic state, in the presence of thiazide or amiloride,
and with increased sodium load and flow rate. Validity of the present
model can be discussed from three aspects: validity of the experimental
data to simulate; validity of the assumptions made to formulate the
model; and validity of the model
calculation.
Validity of the Experimental Data to Simulate
As an estimate of the validity of experimental data on which we based
our numerical model, a summary of collected experimental measurements
is presented in Table 12. Values adopted
in the definition of the penalty function are underlined in Table 12.
Several points deserve mention with regard to the construction of the
present model. 1) Length of distal tubule (Table 12C)
increases as the rat grows. The value used in this simulation was
typical of rats weighing 200 to 300 g. 2) There is a
controversy about the value of luminal electrical potential in the
distal tubule (Table 12, H and I ). Some researchers
(1, 2, 4, 5) reported values significantly more depolarized than
the values reported by others (43, 67). The origin of
such a difference is assumed to be due to the difference in types of
electrode (Ling-Gerard vs. large-tipped) used to impale the tubule (2).
We adopted a value from a study that conducted simultaneous electrical
recording with both Ling-Gerard and large-tipped electrodes (30).
3) Intracellular potential of distal tubular cells had been
reported to be
67.4 to
92 mV (Table 12S). We speculate
that these values tend to underestimate the magnitude of intracellular
hyperpolarization. Some measurements (35, 36) were conducted with
double-barreled electrodes, and their primary aim was measurement of
intracellular ionic activities. Usage of such large-pore electrodes
with mammalian cells likely caused depolarization of intracellular
potential due to large leak current and cell damage. Even in the case
of measurements conducted with Ling-Gerard electrode,
"they [potential recordings] declined rapidly, and were
rarely stable for longer than a few seconds" (24), which
suggests that distal tubular cells were probably damaged by
microelectrode impalement. Thus we have determined basolateral
parameters to make the intracellular potential of model tubular cells
close to the most hyperpolarized value reported (
92 mV). This value
is also close to the intracellular potential (
90 to
95
mV) of a mathematical model of the cortical collecting
tubule (55), which is a nephron segment functionally similar to the
late distal tubule.
Validity of the Assumptions Made to Formulate the Model
We have considered only fluxes of water, sodium, potassium, and
chloride in the present model. Certainly, this is a substantial simplification in view of the fact that distal tubule also transports urea, hydrogen, bicarbonate, and calcium. However, we concentrated on
transport of sodium, potassium, and chloride in this study, because in
magnitude the solute flux ascribable to them (more than 700 pmol/min)
overwhelms fluxes of other solutes [urea, 70 pmol/min (3);
hydrogen/bicarbonate, 50 pmol/min (10); and calcium, 4 pmol/min (13)].
Disregarding the transport of hydrogen and bicarbonate, we could
exclude, from the model, intercalated cells, which are involved
primarily in hydrogen/bicarbonate transport. Consequently, cells in the
present model can be regarded as distal convoluted tubular cells in the
early subsegment and connecting tubular cells or principal cells in the
late subsegment.
We have assumed that the intercellular space has the same hydraulic
pressure and solute concentrations as the basolateral space, because
almost no information is currently available concerning the transport
involved in the intercellular space. This assumption may not be valid
under extreme experimental conditions, since massive water
reabsorption, induced in diabetes insipidus rats with antidiuretic
hormone, is known to be accompanied by expansion of the intercellular
space (66), which suggests the presence of a pressure gradient between
this space and the basolateral space. However, such an expansion has
not been observed in Wistar and Sprague-Dawley rats in the normal state
(66).
Existence of the thiazide-sensitive Na-Cl cotransporter in the luminal
membrane of the distal tubule is supported by physiological (12, 19) as
well as histochemical (47, 49, 68) studies. For example, in in vivo
microperfusion experiments Ellison et al. (19) have shown that thiazide
diuretics abolished sodium reabsorption in the early distal tubule.
However, effects of thiazide on the late distal tubule were less
remarkable and were statistically insignificant. This is consistent
with the histochemical finding that expression of the
thiazide-sensitive Na-Cl cotransporter was considerably weaker in the
late distal tubule than in the early distal tubule (47, 49). Such a
heterogeneous distribution of Na-Cl cotransporter was represented in
our model by enforcing J mcNaCl,max in the early
distal tubule to be 10-fold larger than that in the late distal tubule.
There is ample evidence for the existence of the amiloride-sensitive Na
channel in the apical membrane of the distal tubule. In physiological
experiments (4, 16), luminal amiloride application depolarized the
luminal electrical potential, which is explainable by suppression of
the (inherently electrogenic) amiloride-sensitive Na channel. Although
effects of amiloride on total sodium reabsorption in the distal tubule
were not readily discernible in micropuncture experiments (16), in vivo
microperfusion experiments revealed that there was a modest but
statistically significant reduction of sodium reabsorption when
amiloride was included in the perfusate (62). On the other hand,
histochemical studies also demonstrated that all three subunits (8, 9,
41), which are required to reproduce the fully functional channel, are
detectable in the apical membrane throughout the distal tubule (17). In
the present model, we have assumed separate values of maximal transport
capacities of the apical Na channel
(PmcNa) in the early and late
subsegments. We varied their values independently while optimizing the
model parameters, and the result showed 11-fold larger luminal Na
conductance in the late than in the early subsegments (Table 8). Future
experiments may discover whether such an axial heterogeneity is
actually present. If it were not the case, then we should have to add
another condition (for example, enforcing equal value of
PmcNa in the two subsegments) while searching the optimal model parameters.
Potassium secretion in the distal tubule is considered to be mostly
through K channel in the luminal cell membrane (23). In cortical
collecting tubule (48), the luminal K channel was identified to be ROMK
inward rectifier (32, 70). Also in the distal tubule, ROMK channel was
demonstrated to be present (7), justifying the luminal K channel in the
model. Besides the luminal K channel, it is also suggested that there
exists K-Cl cotransport in the apical membrane of the distal tubule by
in vivo microperfusion experiments (18, 60). Consequently, we
incorporated both K channel and K-Cl cotransporter in the luminal
membrane. However, optimization procedure (i.e., minimizing the penalty
function) resulted in a small value of maximal transport capacity of
apical K-Cl cotransporter in the late distal tubule (Table 8). (Maximal transport capacity of the K-Cl cotransporter in the early distal tubule
was set to a small fixed value in the present model.)
Validity of the Model Calculation
Because of the nonlinearity of the model equations, it is difficult to
demonstrate the appropriateness of the differencing scheme adopted in
this study. As an alternative, we have tested whether the numerical
solutions converge to a limit as we shorten the section width. We made
four models with different section widths: 1) 0.033 cm in the
early and 0.043 cm in the late distal tubule; 2) 0.020 cm in
the early and 0.022 cm in the late distal tubule; 3) 0.010 cm
in the early and the late distal tubule (reported model); and
4) 0.005 cm in the early and the late distal tubule. Using the
same parameters as listed in Table 8, we solved model equations taking
into account electrodiffusive terms. Results are plotted in Fig.
13. Solutions can be seen to converge to
a limit as we take digitization intervals smaller. The solution with
0.010-cm section width and the one with 0.005-cm section width
coincided with less than a 0.5-mV difference, which validates that the
sampling interval adopted in the present model formulation was
practically small enough.

View larger version (10K):
[in this window]
[in a new window]
|
Fig. 13.
Numerical solutions with different digitization intervals. Profiles of
luminal electrical potential obtained in numerical models with
digitization intervals of 0.033, 0.020, 0.010, and 0.005 cm (in the
early distal tubule) were plotted.
|
|
Future Extension of the Model
This study is a first attempt to develop a model of the distal tubule,
and we have concentrated largely on mathematical formulation and
implementation of an efficient program. Accordingly, the model has
several restrictions that should be improved in future revisions. 1) The model was optimized to simulate only results of
micropuncture experiments. It would be interesting to test whether this
model also simulates the results from in vivo microperfusion
experiments (13, 18, 19, 27-29, 60, 61) in which solute
concentrations and flow rate had been manipulated more extensively.
2) The model handles transport of only water, sodium,
potassium, and chloride. Other solutes such as hydrogen, bicarbonate,
magnesium, and calcium should be incorporated, and this is possible
within the mathematical framework and numerical methods used in the
present model. 3) The program lacks flexible user interface and
is not ideally easy to use. This problem can be solved by providing an
intuitive graphical user interface to the program, which would allow
for researchers to test the appropriateness of the model by themselves
as well as to use the model to predict water and solute transport in
physiological and pathophysiological situations.
 |
APPENDIX |
A Kinetic Model of Cation-Anion Cotransporters
The model transporter is a cation-anion cotransporter with a
stoichiometry of 1:1, which was applied to both Na-Cl cotransporter and
K-Cl cotransporter. The state diagram of the model is illustrated in
Fig. 14, in which following assumptions
were made: 1) the transporter has one cation binding site and
one anion binding site; 2) binding to each site occurs in any
order and independently of each other; 3) affinities of binding
sites are symmetrical; 4) transition of the transporter across
the membrane is a much slower process compared with other binding or
dissociation processes, and bound and unbound forms of the transporters
are essentially in equilibrium on each side; and 5) rate
constants of transition of the transporter across the membrane are
symmetrical. The whole transport cycle can be outlined as follows:
binding of cation (C) and anion (A) to an unloaded transporter (T)
converts the transporter to half loaded states (TC or TA) and then to
the fully loaded state (TAC); fully loaded transporter can cross the
membrane to the other side (T'A'C'); and after dissociation of both
ions, the empty transporter (T') crosses the membrane back to return
to the original state (T). It should be noted that transitions between
T and T' or between TAC and T'A'C' do not necessarily imply actual
translocation of the transporter molecule but likely correspond to some
conformational change that forces the binding sites to face the other
side of the membrane.

View larger version (19K):
[in this window]
[in a new window]
|
Fig. 14.
State diagram of cation-anion cotransporter. T, transporter molecule; C
and A, cation and anion molecule, respectively. Prime (') indicates
molecules in side 2. Transition rates are shown beside
arrows.
|
|
From the assumption 4 (see above), we can assume that states T,
TA, TC, and TAC are near equilibrium, and if we use the same symbol for
each kinetic state in Figure 16 and for its concentration, then
|
(AK-1)
|
|
(AK-2)
|
|
(AK-3)
|
where
1 = k1A/k2, and
2 = k3C/k4.
Similarly
|
(AK-4)
|
|
(AK-5)
|
|
(AK-6)
|
where
'1 = k1A'/k2, and
'2 = k3C'/k4.
Sum of the concentrations of all the states are constant
(S)
|
(AK-7)
|
Finally, the total number of transporters in one side of the membrane
should be constant in steady state
|
(AK-8)
|
These eight Eqs. 1 to 8 can be solved
for T, TA, TC, TAC, T', T'A', T'C', and T'A'C', and flux of cation
(FC) or anion (FA) can be expressed as
where KA = k2/k1 is the intrinsic
dissociation constant of anion binding site, and KC = k4/k3 is the intrinsic
dissociation constant of cation binding site. If we further assume that
k5 = k6, the above equation
reduces to Eq. 4 in Table 2, in which J pq
,max = k6S. We note that maximal transport rate is
1/2 × J pq
,max in this notation.
 |
ACKNOWLEDGEMENTS |
We thank Drs. N. Yamashita and K. Cho for critical reading of the
manuscript and helpful discussions.
 |
FOOTNOTES |
1
A practical way to lessen the difference between
the two models may be to change parameter values gradually at the
junction of the two subsegments rather than abruptly as in the present models. Such an arrangement of parameters is justified, since transition between these subsegments is a gradual process in the case
of the rat (42).
2
However, ad hoc exploration revealed that models
with significant (up to 8 pmol/min) potassium secretion through luminal K-Cl cotransporter could also be constructed by readjusting the parameters, although the penalty function increased to the level of
0.5.
Address for reprint requests and other correspondence: H. Chang, Health
Service Center, Univ. of Tokyo, 3-8-1 Komaba, Meguro-ku, Tokyo
153-8902, Japan.
Received 17 February 1998; accepted in final form 28 January 1999.
 |
REFERENCES |
1.
Allen, G. G.,
and
L. J. Barratt.
Effect of adolesterone on the transepithelial potential difference of the rat distal tubule.
Kidney Int.
19:
678-686,
1981[Medline].
2.
Allen, G. G.,
and
L. J. Barratt.
Origin of positive transepithelial potential difference in early distal segments of rat kidney.
Kidney Int.
27:
622-629,
1985[Medline].
3.
Armsen, T.,
and
H. W. Reinhardt.
Transtubular movement of urea at different degrees of water diuresis.
Pflügers Arch.
326:
270-280,
1971[Medline].
4.
Barratt, L. J.
The effect of amiloride on the transepithelial potential difference of the distal tubule of the rat kidney.
Pflügers Arch.
361:
251-254,
1976[Medline].
5.
Barratt, L. J.,
F. C. Rector, Jr.,
J. Kokko,
C. C. Tisher,
and
D. W. Seldin.
Transepithelial potential difference profile of the distal tubule of the rat kidney.
Kidney Int.
8:
368-375,
1975[Medline].
6.
Beck, F.-X.,
A. Dörge,
R. Rick,
M. Schramm,
and
K. Thurau.
Effect of potassium adaptation on the distribution of potassium, sodium, and chloride across the apical membrane of renal tubular cells.
Pflügers Arch.
409:
477-485,
1987[Medline].
7.
Boim, M. A.,
K. Ho,
M. E. Shuck,
M. J. Bienkowski,
J. Block,
J. L. Slightom,
Y. Yang,
B. M. Brenner,
and
S. C. Hebert.
ROMK inwardly rectifying ATP-sensitive K+ channel. II. Cloning and distribution of alternative forms.
Am. J. Physiol.
268 (Renal Fluid Electrolyte Physiol. 37):
F1132-F1140,
1995[Abstract/Free Full Text].
8.
Canessa, C. M.,
J. D. Horisberger,
and
B. C. Rossier.
Epithelial sodium channel related to proteins involved in neurodegeneration.
Nature
361:
467-470,
1993[Medline].
9.
Canessa, C. M.,
L. Schild,
G. Buell,
B. Thorens,
I. Gautschi,
J.-D. Horisberger,
and
B. C. Rossier.
Amiloride-sensitive epithelial Na+ channel is made of three homologous subunits.
Nature
367:
463-467,
1994[Medline].
10.
Capasso, G.,
G. Malnic,
T. Wang,
and
G. Giebisch.
Acidification in mammalian cortical distal tubule.
Kidney Int.
45:
1543-1554,
1994[Medline].
11.
Clapp, J. R.,
and
R. R. Robinson.
Osmolality of distal tubular fluid in the dog.
J. Clin. Invest.
45:
1847-1853,
1966[Medline].
12.
Costanzo, L. S.
Localization of diuretic action in microperfused rat distal tubules: Ca and Na transport.
Am. J. Physiol.
248 (Renal Fluid Electrolyte Physiol. 17):
F527-F535,
1985[Medline].
13.
Costanzo, L. S.,
and
E. E. Windhager.
Calcium and sodium transport by the distal convoluted tubule of the rat.
Am. J. Physiol.
235 (Renal Fluid Electrolyte Physiol. 4):
F492-F506,
1978[Abstract/Free Full Text].
14.
Cussler, E. L.
Diffusion: Mass Transfer in Fluid Systems (2nd ed.). Cambridge, UK: Cambridge University Press, 1997.
15.
De Mello-Aires, M.,
G. Giebisch,
G. Malnic,
and
P. F. Curran.
Kinetics of potassium transport across single distal tubules of rat kidney.
J. Physiol. Lond.
232:
47-70,
1973[Medline].
16.
Duarte, C. G.,
F. Chomety,
and
G. Giebisch.
Effect of amiloride, ouabain, and furosemide on distal tubule function in the rat.
Am. J. Physiol.
221:
632-639,
1971[Medline].
17.
Duc, C.,
N. Farman,
C. M. Canessa,
J.-P. Bonvalet,
and
B. C. Rossier.
Cell-specific expression of epithelial sodium channel
,
, and
subunits in adolsterone-responsive epithelia from the rat: localization by in situ hybridization and immunocytochemistry.
J. Cell Biol.
127:
1907-1921,
1994[Abstract].
18.
Ellison, D. H.,
H. Velázquez,
and
F. S. Wright.
Stimulation of distal potassium secretion by low lumen chloride in the presence of barium.
Am. J. Physiol.
248 (Renal Fluid Electrolyte Physiol. 17):
F638-F649,
1985[Abstract/Free Full Text].
19.
Ellison, D. H.,
H. Velázquez,
and
F. S. Wright.
Thiazide-sensitive sodium chloride cotransport in early distal tubule.
Am. J. Physiol.
253 (Renal Fluid Electrolyte Physiol. 22):
F546-F554,
1987[Abstract/Free Full Text].
20.
Fernandes, P. L.,
and
H. G. Ferreira.
A mathematical model of rabbit cortical thick ascending limb of the Henle's loop.
Biochim. Biophys. Acta
1064:
111-123,
1991[Medline].
21.
Fushimi, K.,
S. Uchida,
Y. Hara,
Y. Hirata,
F. Marumo,
and
S. Sasaki.
Cloning and expression of apical membrane water channel of rat kidney collecting tubule.
Nature
361:
549-552,
1993[Medline].
22.
Gamba, G.,
S. N. Saltzberg,
M. Lombardi,
A. Miyanoshita,
J. Lytton,
M. A. Hediger,
B. M. Brenner,
and
S. C. Hebert.
Primary structure and functional expression of a cDNA encoding the thiazide-sensitive, electroneutral sodium-chloride cotransporter.
Proc. Natl. Acad. Sci. USA
90:
2749-2753,
1993[Abstract].
23.
Giebisch, G.
Renal potassium channles: an overview.
Kidney Int.
48:
1004-1009,
1995[Medline].
24.
Giebisch, G.,
G. Malnic,
R. M. Klose,
and
E. E. Windhager.
Effect of ionic substitutions on distal potential differences in rat kidney.
Am. J. Physiol.
211:
560-568,
1966[Medline].
25.
Gillen, C. M.,
S. Brill,
J. A. Payne,
and
B. Forbush III.
Molecular cloning and functional expression of the K-Cl cotransporter from rabbit, rat, and human.
J. Biol. Chem.
271:
16237-16244,
1996[Abstract/Free Full Text].
26.
Goldman, D. E.
Potential, impedance, and rectification in membranes.
J. Gen. Physiol.
27:
37-60,
1943[Free Full Text].
27.
Good, D. W.,
H. Velázquez,
and
F. S. Wright.
Luminal influences on potassium secretion: low sodium concentration.
Am. J. Physiol.
246 (Renal Fluid Electrolyte Physiol. 15):
F609-F619,
1984[Abstract/Free Full Text].
28.
Good, D. W.,
and
F. S. Wright.
Luminal influences on potassium secretion: sodium concentration and fluid flow rate.
Am. J. Physiol.
236 (Renal Fluid Electrolyte Physiol. 5):
F192-F205,
1979[Abstract/Free Full Text].
29.
Good, D. W.,
and
F. S. Wright.
Luminal influences on potassium secretion: transepithelial voltage.
Am. J. Physiol.
239 (Renal Fluid Electrolyte Physiol. 8):
F289-F298,
1980[Medline].
30.
Hayslett, J. P.,
E. L. Boulpaep,
and
G. H. Giebisch.
Factors influencing transepithelial potential difference in mammalian distal tubule.
Am. J. Physiol.
234 (Renal Fluid Electrolyte Physiol. 3):
F182-F191,
1978[Medline].
31.
Hille, B.
Ionic Channels of Excitable Membranes (2nd ed.). Sunderland, MA: Sinauer, 1992.
32.
Ho, K.,
C. G. Nichols,
W. J. Lederer,
J. Lytton,
P. M. Vassilev,
M. V. Kanazirska,
and
S. C. Hebert.
Cloning and expression of an inwardly rectifying ATP-regulated potassium channel.
Nature
362:
31-38,
1993[Medline].
33.
Hodgkin, A. L.,
and
B. Katz.
The effect of sodium ions on the electrical activity of the giant axon of the squid.
J. Physiol. (Lond.)
108:
37-77,
1949.
34.
Katz, A. I.,
A. Doucet,
and
F. Morel.
Na-K-ATPase activity along the rabbit, rat, and mouse nephron.
Am. J. Physiol.
237 (Renal Fluid Electrolyte Physiol. 6):
F114-F120,
1979[Abstract/Free Full Text].
35.
Khuri, R. N.,
S. K. Agulian,
and
K. Bogharian.
Electrochemical potentials of chloride in distal renal tubule of the rat.
Am. J. Physiol.
227:
1352-1355,
1974[Medline].
36.
Khuri, R. N.,
S. K. Agulian,
and
A. Kalloghlian.
Intracellular potassium in cells of the distal tubule.
Pflügers Arch.
335:
297-308,
1972[Medline].
37.
Kunau, R. T., Jr.,
H. L. Webb,
and
S. C. Borman.
Characteristics of sodium reabsorption in the loop of Henle and distal tubule.
Am. J. Physiol.
227:
1181-1191,
1974[Medline].
38.
Kunau, R. T., Jr.,
H. L. Webb,
and
S. C. Borman.
Characteristics of the relationship between the flow rate of tubular fluid and potassium transport in the distal tubule of the rat.
J. Clin. Invest.
54:
1488-1495,
1974[Medline].
39.
Kunau, R. T., Jr.,
D. R. Weller,
and
H. L. Webb.
Clarification of the site of action of chlorothiazide in the rat nephron.
J. Clin. Invest.
56:
401-407,
1975[Medline].
40.
Layton, H. E.,
and
E. B. Pitman.
A dynamic numerical method for model of renal tubules.
Bull. Math. Biol.
56:
547-565,
1994[Medline].
41.
Lingueglia, E.,
N. Voilley,
R. Waldmann,
M. Lazdunski,
and
P. Barbry.
Expression cloning of an epithelial amiloride-sensitive Na+ channel.
FEBS Lett.
318:
95-99,
1993[Medline].
42.
Madsen, K. M.,
and
C. C. Tisher.
Structural-functional relationships along the distal nephron.
Am. J. Physiol.
250 (Renal Fluid Electrolyte Physiol. 19):
F1-F15,
1986[Abstract/Free Full Text].
43.
Malnic, G.,
and
G. Giebisch.
Some electrical properties of distal tubular epithelium in the rat.
Am. J. Physiol.
223:
797-808,
1972[Medline].
44.
Malnic, G.,
R. M. Klose,
and
G. Giebisch.
Micropuncture study of renal potassium excretion in the rat.
Am. J. Physiol.
206:
674-686,
1964.
45.
Malnic, G.,
R. M. Klose,
and
G. Giebisch.
Micropuncture study of distal tubular potassium and sodium transport in rat nephron.
Am. J. Physiol.
211:
529-547,
1966[Medline].
46.
Mejia, R.,
M. F. Flessner,
and
M. A. Knepper.
Model of ammonium and bicarbonate transport along LDL: implications of alkalinization of luminal fluid.
Am. J. Physiol.
264 (Renal Fluid Electrolyte Physiol. 33):
F397-F403,
1993[Abstract/Free Full Text].
47.
Obermüller, N.,
P. Bernstein,
H. Velázquez,
R. Reilly,
D. Moser,
D. H. Ellison,
and
S. Bachmann.
Expression of the thiazide-sensitive Na-Cl cotransporter in rat and human kidney.
Am. J. Physiol.
269 (Renal Fluid Electrolyte Physiol. 38):
F900-F910,
1995[Abstract/Free Full Text].
48.
Palmer, L. G.,
H. Choe,
and
G. Frindt.
Is the secretory K channel in the rat CCT ROMK?
Am. J. Physiol.
273 (Renal Physiol. 42):
F404-F410,
1997[Abstract/Free Full Text].
49.
Plotkin, M. D.,
M. R. Kaplan,
J. W. Verlander,
W.-S. Lee,
D. Brown,
E. Poch,
S. R. Gullans,
and
S. C. Hebert.
Localization of the thiazide sensitive Na-Cl cotransporter, rTSC1, in the rat kidney.
Kidney Int.
50:
174-183,
1996[Medline].
50.
Press, W. H.,
B. P. Flannery,
S. A. Teukolsky,
and
W. T. Vetterling.
Numerical Recipes in C. Cambridge, NY: Cambridge Univ. Press, 1988.
51.
Reineck, H. J.,
R. W. Osgood,
T. F. Ferris,
and
J. H. Stein.
Potassium transport in the distal tubule and collecting duct of the rat.
Am. J. Physiol.
229:
1403-1409,
1975[Medline].
52.
Stanton, B. A.,
D. Biemesderfer,
J. B. Wade,
and
G. Giebisch.
Structural and functional study of the rat distal nephron: Effects of potassium adaptation and depletion.
Kidney Int.
19:
36-48,
1981[Medline].
53.
Stein, W. D.
Intrinsic, apparent, and effective affinities of co- and countertransport systems.
Am. J. Physiol.
250 (Cell Physiol. 19):
C523-C533,
1986[Abstract/Free Full Text].
54.
Stolte, H.,
J. P. Brecht,
M. Wiederholt,
and
K. Hierholzer.
Einfluss von Adrenalektomie und Glukocorticoiden auf die Wasserpermeabilität cortikaler Nephronabschnitte der Rattenniere.
Pflügers Arch.
299:
99-127,
1968.
55.
Strieter, J.,
J. L. Stephenson,
G. Giebisch,
and
A. M. Weinstein.
A mathematical model of the rabbit cortical collecting tubule.
Am. J. Physiol.
263 (Renal Fluid Electrolyte Physiol. 32):
F1063-F1075,
1992[Abstract/Free Full Text].
56.
Stroustrup, B.
The C++ Programming Language (3rd ed.). Reading, MA: Addison-Wesley, 1997.
57.
Taniguchi, J.,
K. Tabei,
and
M. Imai.
Profiles of water and solute transport along long-loop descending limb: analysis by mathematical model.
Am. J. Physiol.
252 (Renal Fluid Electrolyte Physiol. 21):
F393-F402,
1987[Abstract/Free Full Text].
58.
Ullrich, K. J.,
G. Rumrich,
and
G. Fuchs.
Wasserpermeabilität und transtubulärer Wasserfluss cortikaler Nephronabschnitte bei verschiedenen Diurezeständen.
Pflügers Arch.
280:
99-119,
1964.
59.
Ullrich, K. J.,
B. Schmidt-Nielsen,
R. O'Dell,
G. Pehling,
C. W. Gottschalk,
W. E. Lassiter,
and
M. Mylle.
Micropuncture study of composition of proximal and distal tubular fluid in rat kidney.
Am. J. Physiol.
204:
527-531,
1963.
60.
Velázquez, H.,
D. H. Ellison,
and
F. S. Wright.
Chloride-dependent potassium secretion in early and late renal distal tubules.
Am. J. Physiol.
253 (Renal Fluid Electrolyte Physiol. 22):
F555-F562,
1987[Abstract/Free Full Text].
61.
Velázquez, H.,
D. W. Good,
and
F. S. Wright.
Mutual dependence of sodium and chloride absorption by renal distal tubule.
Am. J. Physiol.
247 (Renal Fluid Electrolyte Physiol. 16):
F904-F911,
1984[Medline].
62.
Velázquez, H.,
and
F. S. Wright.
Effects of diuretic drugs on Na, Cl, and K transport by rat distal tubule.
Am. J. Physiol.
250 (Renal Fluid Electrolyte Physiol. 19):
F1013-F1023,
1986[Medline].
63.
Verkman, A. S.,
and
R. J. Alpern.
Kinetic transport model for cellular regulation of pH and solute concentration in the renal proximal tubule.
Biophys. J.
51:
533-546,
1987[Abstract].
64.
Weinstein, A. M.
A mathematical model of the rat proximal tubule.
Am. J. Physiol.
250 (Renal Fluid Electrolyte Physiol. 19):
F860-F873,
1986[Medline].
65.
Weinstein, A. M.
Mathematical models of tubular transport.
Annu. Rev. Physiol.
56:
691-709,
1994[Medline].
66.
Woodhall, P. B.,
and
C. C. Tisher.
Response of the distal tubule and cortical collecting duct to vasopressin in the rat.
J. Clin. Invest.
52:
3095-3108,
1973[Medline].
67.
Wright, F. G.
Increasing magnitude of electrical potential along the renal distal tubule.
Am. J. Physiol.
220:
624-638,
1971[Medline].
68.
Yang, T.,
Y. G. Huang,
I. Singh,
J. Schnermann,
and
J. P. Briggs.
Localization of bumetanide- and thiazide-sensitive Na-K-Cl cotransporters along the rat nephron.
Am. J. Physiol.
271 (Renal Fluid Electrolyte Physiol. 40):
F931-F939,
1996[Abstract/Free Full Text].
69.
Zahler, R.,
Z.-T. Zhang,
M. Manor,
and
W. F. Boron.
Sodium kinetics of Na,K-ATPase
isoforms in intact transfected HeLa cells.
J. Gen. Physiol.
110:
201-213,
1997[Abstract/Free Full Text].
70.
Zhou, H.,
S. S. Tate,
and
L. G. Palmer.
Primary structure and functional properties of an epithelial K channel.
Am. J. Physiol.
266 (Cell Physiol. 35):
C809-C824,
1994[Abstract/Free Full Text].
Am J Physiol Renal Physiol 276(6):F931-F951
0002-9513/99 $5.00
Copyright © 1999 the American Physiological Society