Why Sampling Location Belongs in the Model
A neurotransmitter released into a synaptic cleft does not instantly have the same concentration throughout surrounding tissue. It can diffuse out of the cleft, be taken up from extracellular fluid into cells, and be metabolized in more than one location. A sample taken from the cleft, extracellular tissue, or a tissue homogenate therefore need not represent the same concentration, even when all measurements concern the same transmitter.
A useful model keeps these spaces distinct and accounts for every transfer and loss. This extends the mass-balance approach used in transmitter turnover models: instead of treating concentration as if it occupied one uniform space, we assign an amount to each compartment and specify how material moves between them.
Coupled Mass-Balance Model
Consider a constant release rate \(R\) into the cleft. Let \(G\) be the cleft-to-extracellular diffusion conductance, with units of volume per time. The signed diffusion flux from cleft to extracellular tissue is \(G(C_s-C_e)\); if \(C_e>C_s\), this flux is negative and transport is in the opposite direction. Let \(CL_u\) be an uptake clearance from extracellular fluid into cells. Let \(k_s\), \(k_e\), and \(k_i\) be first-order metabolic rate constants in the cleft, extracellular tissue, and intracellular tissue. The balances are:
Every term on the right-hand side is an amount per time. For example, \(CL_uC_e\) has units of volume per time multiplied by amount per volume. It is a loss from extracellular tissue and, in the model, the identical gain to the intracellular compartment. Uptake is therefore a transfer, not a whole-system loss; intracellular metabolism is the loss that follows it.
The model assumes fixed compartment volumes, uniform concentration within each compartment, linear diffusion and uptake, and first-order metabolism. It also assumes that release enters the cleft and that no other source, transfer, or elimination pathway matters over the time interval considered. These are modeling assumptions, not universal properties of neural tissue. The linear formulation is most appropriate when transport and metabolism can be approximated as proportional to concentration across the range of interest.
Steady State and Its Mass-Balance Check
At steady state, all three concentration derivatives are zero. A compact derivation shows how the steady concentrations depend on the model parameters. Define the first-order loss capacities \(a=k_sV_s\) and \(b=k_eV_e\), both in volume per time. From the intracellular balance, \(CL_uC_e=k_iV_iC_i\). The extracellular balance gives \(G(C_s-C_e)=(b+CL_u)C_e\), so:
The cleft balance requires release to equal cleft metabolism plus net diffusion out of the cleft. Substituting the extracellular balance into that requirement gives the extracellular steady concentration:
For positive \(G\), nonnegative loss capacities, and positive \(R\), the denominator is positive and the steady concentrations are nonnegative. Then \(C_s\) follows from the concentration ratio above, and \(C_i=CL_uC_e/(k_iV_i)\) when \(k_iV_i>0\). If intracellular metabolism is absent while uptake continues, intracellular amount accumulates rather than reaching this finite steady state; the formula for \(C_i\) does not apply.
Worked Example: Solving a Three-Compartment Steady State
Consider a hypothetical model with \(V_s=1\) µL, \(V_e=1\) µL, and \(V_i=0.50\) µL. Set \(R=12\) fmol/min, \(G=4\) µL/min, \(k_s=k_e=1\) min\(^{-1}\), \(CL_u=2\) µL/min, and \(k_i=0.40\) min\(^{-1}\). Because \(1\) nM in \(1\) µL corresponds to \(1\) fmol, the concentrations can be reported in nM and fluxes in fmol/min.
At steady state, the extracellular balance is \(4(C_s-C_e)=2C_e+1C_e\), or \(4C_s=7C_e\). The cleft balance is \(12=1C_s+4(C_s-C_e)\). Using the extracellular balance, this becomes \(12=C_s+3C_e\). Substituting \(C_s=7C_e/4\) gives \(12=(19/4)C_e\), so \(C_e=48/19=2.526\) nM and \(C_s=84/19=4.421\) nM.
The intracellular balance gives \(2C_e=0.40(0.50)C_i\), so \(C_i=10C_e=480/19=25.26\) nM. Check the transfers and losses using the unrounded concentrations: diffusion is \(4(84/19-48/19)=144/19=7.579\) fmol/min. Cleft metabolism is \(84/19=4.421\) fmol/min, and these sum to \(228/19=12\) fmol/min of release. In extracellular tissue, uptake is \(2(48/19)=96/19=5.053\) fmol/min and metabolism is \(48/19=2.526\) fmol/min; together they equal the incoming diffusion flux, \(144/19\) fmol/min. Intracellular metabolism is \(0.40(0.50)(480/19)=96/19=5.053\) fmol/min, matching uptake. These equalities verify both compartment balances and whole-system conservation.
What a Sample Measures
The model predicts different local concentrations, not one interchangeable “tissue concentration.” A sample that isolates extracellular fluid is associated with \(C_e\), whereas a sample of intracellular material is associated with \(C_i\). A homogenate that combines extracellular and intracellular material instead measures a volume-weighted concentration, provided the assay recovers the analyte equally from both spaces. If the cleft is excluded, that concentration is:
With the steady-state values in the example, the extracellular amount is \(1(48/19)=2.526\) fmol and the intracellular amount is \(0.50(480/19)=12.63\) fmol. Thus the predicted homogenate concentration is \((2.526+12.63)/(1+0.50)=10.11\) nM. It is neither the extracellular concentration of 2.526 nM nor the intracellular concentration of 25.26 nM. This calculation excludes cleft contents by assumption; if a sampling procedure includes the cleft, its amount and volume must also be included.
Worked Example: Uptake Changes the Compartments Differently
Use the same hypothetical parameters, but reduce \(CL_u\) from \(2\) to \(0.50\) µL/min. This represents a change in the assumed uptake process, with all other model inputs and parameters held fixed. At steady state, the extracellular balance is \(4(C_s-C_e)=(1+0.50)C_e\), so \(C_s=(11/8)C_e\). The cleft balance then gives \(12=C_s+1.50C_e=(23/8)C_e\).
Therefore \(C_e=96/23=4.174\) nM and \(C_s=132/23=5.739\) nM. The intracellular balance gives \(C_i=0.50C_e/[0.40(0.50)]=2.5C_e=240/23=10.43\) nM. Relative to the original steady state, extracellular concentration rises, while intracellular concentration falls. It would be incorrect to infer that all compartments must move in the same direction when uptake changes: uptake affects both the extracellular loss balance and the intracellular input.
A flux check confirms the new steady state. Diffusion is \(4(132/23-96/23)=144/23=6.261\) fmol/min. Extracellular metabolism is \(96/23=4.174\) fmol/min and uptake is \(0.50(96/23)=48/23=2.087\) fmol/min; these sum to \(144/23\) fmol/min. Cleft metabolism is \(132/23=5.739\) fmol/min, so cleft metabolism plus diffusion is \((132+144)/23=12\) fmol/min, equal to release. Intracellular metabolism, \(0.40(0.50)(240/23)=48/23\) fmol/min, equals uptake.
Worked Example: A Change in Release Has Location-Specific Timing
Suppose the first example is at its steady state, and release is then switched from \(12\) fmol/min to zero. Concentrations cannot jump in this model because the compartments have finite volumes. Immediately after the switch, \(C_s=84/19\) nM and \(C_e=48/19\) nM still hold. The cleft balance now gives \(V_s(dC_s/dt)=0-84/19-144/19=-12\) fmol/min, so \(dC_s/dt=-12\) nM/min because \(V_s=1\) µL.
At that same instant the extracellular balance is still zero: diffusion \(144/19\) fmol/min is balanced by uptake \(96/19\) fmol/min and metabolism \(48/19\) fmol/min. The intracellular balance is also initially zero because uptake still equals metabolism. Thus the cleft concentration starts declining immediately, while the other two concentrations initially have zero slope. Their slopes can change subsequently as the cleft-to-tissue gradient and the uptake input change. A measurement time and sampling location are both needed to interpret a concentration after a transient perturbation.
Diagnostics, Alternatives, and Limitations
A practical diagnostic is to check each compartment separately before interpreting a predicted concentration. At steady state, every compartment’s total inflow must equal its total outflow. In a transient simulation, the difference between inflow and outflow must equal \(V_j\,dC_j/dt\) for that compartment. Summing the balances cancels internal transfers: diffusion and uptake must not appear as net whole-system sources or losses. If an internal flux is counted as a loss but not as a gain elsewhere, the model violates mass conservation.
Other warning signs include negative predicted concentrations from a model with nonnegative inputs and losses, incompatible units in a flux, or a claimed finite intracellular steady state despite uptake and no intracellular elimination. A mismatch between a predicted concentration and an assay result does not by itself identify which parameter is wrong. The sample may combine compartments, include a different amount of cleft material, or reflect spatial variation that the well-mixed assumptions omit.
A single well-mixed compartment is simpler and may be adequate when internal concentration gradients equilibrate rapidly relative to the process being studied, or when the intended observation is an average that can be mapped to that compartment. It cannot, by itself, predict different cleft and tissue concentrations. The coupled compartments add this distinction while remaining relatively interpretable. When concentration varies continuously with distance, a spatial diffusion model may be more appropriate than a small set of well-mixed spaces, but it requires spatial information and additional parameters. The three-compartment model also cannot represent saturable transport, binding, explicit cellular subtypes, or changing volumes unless those features are added and supported by data.
Check Your Understanding
Try each problem, then compare your reasoning with the solutions below.
- If extracellular concentration exceeds cleft concentration, what is the sign of \(G(C_s-C_e)\), and what direction does it represent?
- In the first worked example, what is the total modeled amount at steady state, excluding any material outside the three compartments?
- Why does an extracellular concentration not necessarily equal a tissue-homogenate concentration?
- For the first example, immediately after release switches off, which compartment has a nonzero concentration slope, and what is that slope?
Solutions. (1) The flux is negative because \(C_s-C_e<0\); the signed convention means net diffusion is from extracellular tissue into the cleft. (2) The amounts are \(V_sC_s=84/19\) fmol, \(V_eC_e=48/19\) fmol, and \(V_iC_i=240/19\) fmol. Their sum is \(372/19=19.58\) fmol. (3) A homogenate combines amounts from its included spaces and divides by their combined volume; with the example’s extracellular and intracellular spaces it is \(192/19=10.11\) nM, not \(C_e=48/19=2.526\) nM. (4) Immediately after the switch, the cleft slope is \(-12\) nM/min. The extracellular and intracellular slopes are initially zero because their fluxes are still balanced at the pre-switch concentrations.