From Binding to Circuit Output
A measured or modeled dopamine or serotonin concentration does not by itself specify a neural response. The fraction of receptors occupied depends on affinity and competing ligands; the signal produced by occupied receptors depends on efficacy and cellular coupling; and circuit output can depend nonlinearly on signals from several receptor subtypes. Keeping these stages distinct helps prevent a common modeling error: treating receptor occupancy as if it were already a measure of behavioral or circuit effect.
The framework here is a compact, illustrative mapping. It is useful for connecting transmitter exposure to a hypothesized receptor-mediated response, but its parameters must be interpreted in the context of the receptor, cell type, and circuit being modeled. Dopamine D1-like and D2-like receptors, and the many serotonin receptor subtypes, do not have universal circuit effects independent of their anatomical and cellular context.
Occupancy and Receptor-Specific Signal
These expressions assume a single class of independent binding sites, equilibrium, and known free concentrations at the receptor. They do not use total tissue concentration as a substitute for free concentration. When there is only one ligand, occupancy rises with concentration and approaches 1 as \(C\) becomes large relative to \(K_D\). For example, if \(C=K_D\), then \(\theta=1/2\).
Occupancy is not efficacy. Let \(\epsilon_l\) denote the efficacy of ligand \(l\), scaled so that a reference full agonist has \(\epsilon=1\), a partial agonist has \(0<\epsilon<1\), and a neutral antagonist has \(\epsilon=0\). In a simple model, efficacy-weighted occupancy is \(x=\sum_l \epsilon_l\theta_l\). This is a simplified measure of receptor stimulus, not a universal molecular quantity. An antagonist can occupy receptors while contributing no positive stimulus in this representation; by competing with an endogenous agonist, it can nevertheless reduce that agonist's occupancy and signal.
A saturating transduction function maps receptor stimulus \(x\) to a normalized receptor-specific signal \(r\):
Here, \(r\) ranges from 0 toward 1, and \(\tau\) is a dimensionless transduction-gain parameter that summarizes receptor abundance and downstream coupling in this model. A larger \(\tau\) represents greater effective amplification, often described as receptor reserve: a substantial response can occur without occupying every receptor. The model is an operational approximation, not a claim that receptor number and coupling can always be reduced to one identifiable parameter.
Why the Transduction Function Saturates
For \(r=\tau x/(1+\tau x)\), differentiation with respect to \(x\) gives
So increasing stimulus cannot decrease this model's receptor signal. The slope is strictly positive when \(\tau>0\), and is zero when \(\tau=0\). Also, \(r=0\) at \(x=0\), while \(r\) approaches 1 as \(x\) increases without bound when \(\tau>0\). The diminishing slope at larger \(x\) is the model's representation of saturation. These properties are built into the chosen function; they should not be mistaken for empirical proof that a particular receptor pathway follows this exact curve.
Worked Example: Occupancy and Partial Agonism
Suppose a hypothetical agonist has free concentration \(C=2.0\) nM and \(K_D=3.0\) nM at a receptor. Assume no competing ligands, efficacy \(\epsilon=1\), and transduction gain \(\tau=0.80\). Calculate occupancy and normalized signal.
The single-ligand occupancy is
The efficacy-weighted stimulus is \(x=\epsilon\theta=(1)(0.400)=0.400\). Therefore,
In this model, 40.0% occupancy corresponds to a normalized receptor signal of 0.2424, not 0.400. The signal is unitless. It is not a percentage of neurons firing or a direct measure of a clinical effect.
Now suppose a hypothetical partial agonist instead has occupancy \(0.750\), efficacy \(\epsilon=0.350\), and \(\tau=2.00\). Then \(x=(0.350)(0.750)=0.2625\), and
Although this partial agonist has lower efficacy, its larger transduction gain in this hypothetical setting produces a larger normalized signal than in the first calculation. The comparison illustrates why efficacy and receptor reserve must not be conflated.
Receptor Reserve and Partial Agonism
Receptor reserve means that the system can produce a substantial response without requiring occupancy of all available receptors. In the model, a higher \(\tau\) increases \(r\) at a given positive stimulus \(x\). Consequently, a partial agonist may produce a considerable signal in a system with high reserve, while producing a smaller signal in a system with low reserve. The term “partial agonist” describes efficacy relative to a reference agonist in a specified assay or system; it does not guarantee one fixed response across tissues.
Worked Example: The Same Partial Agonist in Two Systems
Consider a hypothetical partial agonist with efficacy \(\epsilon=0.400\) and occupancy \(\theta=0.800\) in two receptor systems. Assume the same occupancy and efficacy in both systems, but \(\tau=0.500\) in the first and \(\tau=5.00\) in the second.
In both systems, \(x=\epsilon\theta=(0.400)(0.800)=0.320\). In the lower-gain system,
In the higher-gain system,
The same occupancy and efficacy yield different signals because the modeled transduction gains differ. For comparison, a full agonist with \(\epsilon=1\), the same occupancy \(0.800\), and \(\tau=5.00\) has \(x=0.800\), so \(r=(5.00)(0.800)/[1+(5.00)(0.800)]=4.00/5.00=0.800\). The partial agonist produces a smaller signal than that full agonist in the same system, but its response is not simply 40% of the full agonist response.
Combining Opposing Subtype Signals
A circuit model can combine receptor-specific signals with signed weights. A positive weight represents a contribution that increases the chosen circuit input; a negative weight represents a contribution that decreases it. These signs are model-specific. They should be justified for the relevant neuronal populations and circuit endpoint rather than assigned from the transmitter name alone.
One possible nonlinear mapping from a net input \(z\) to a bounded output \(y\) is the logistic function:
Here \(b\) is a baseline input, \(w_j\) is the signed weight for receptor subtype \(j\), and \(r_j\) is its normalized signal. The logistic output lies between 0 and 1. Unless the model explicitly defines \(y\) as a probability, it is better interpreted as a normalized circuit output, not as the probability of an event. A weighted sum is one assumption about integration; real circuits may include additional interactions, delays, and state dependence.
Worked Example: Opposing Dopamine and Serotonin Inputs
Consider a hypothetical circuit with one positively weighted dopamine-receptor signal and one opposing serotonin-receptor signal. For dopamine, let \(\theta_D=1.00\), \(\epsilon_D=1\), and \(\tau_D=1.20\). For serotonin, let \(\theta_S=1.00\), \(\epsilon_S=1\), and \(\tau_S=1.00\). Let their circuit weights be \(w_D=1.40\) and \(w_S=-0.750\), with baseline \(b=0\).
The dopamine signal is
The serotonin signal is \(r_S=(1.00)(1.00)/[1+(1.00)(1.00)]=1/2=0.500\). Its weighted opposing contribution is \((-0.750)(0.500)=-0.375\). Thus the net input, calculated using the unrounded dopamine signal, is
The circuit output is
The normalized output is approximately 0.5960. The positive dopamine contribution exceeds the opposing serotonin contribution in this hypothetical parameterization. The result depends on the chosen weights and transduction parameters; it does not establish that dopamine universally excites, or serotonin universally inhibits, neural circuits.
Assumptions, Diagnostics, and Alternatives
Before applying this mapping, check the assumptions that carry its interpretation. The occupancy calculation assumes receptor-level free concentrations and binding equilibrium. The simple receptor signal assumes that efficacy-weighted occupancy is a useful stimulus summary and that a saturating function with a single gain parameter is adequate. The circuit equation assumes that subtype signals combine through specified weights and that the chosen output transformation is appropriate for the endpoint. These are modeling assumptions to test, not consequences of measuring occupancy.
Useful diagnostics include plotting predicted receptor signal against stimulus and predicted circuit output against net input; checking that the functions behave as intended over the observed range; and comparing predicted responses with independent measures of receptor-proximal signaling or circuit activity. Residual patterns, systematic underprediction at high exposure, or different response curves across cell types may indicate missing mechanisms, a wrong transduction shape, or an inappropriate shared \(\tau\). If the model estimates occupancy and efficacy only through their product, those quantities may be difficult to distinguish from response data alone.
The main alternatives answer different questions. An occupancy-only model stops at binding and is appropriate only when occupancy itself is the endpoint of interest. A linear occupancy-response model is simpler and may be adequate over a narrow range, but it does not represent saturation and can predict unbounded signal if extrapolated. The operational saturating model represents reserve compactly but has limited mechanistic detail. More detailed receptor-state or signaling-cascade models can represent multiple steps, but require more parameters and data. A circuit model with signed subtype weights adds integration assumptions; it should be compared with a model omitting selected pathways or using a different justified integration rule.
A frequent mistake is to read “60% occupancy” as “60% receptor response,” then carry that number directly into a circuit endpoint. The correction is to calculate occupancy, apply efficacy and transduction, and then combine subtype signals using an explicit circuit model. A second pitfall is to label a subtype intrinsically excitatory or inhibitory. The sign of its effect depends on receptor location, cell identity, and the measured circuit output.
Check Your Understanding
Try each question, then compare your reasoning with the solutions below.
- A ligand has free concentration \(C=4.0\) nM and \(K_D=4.0\) nM, with no competitors. What is its occupancy?
- At fixed positive stimulus \(x\), how does increasing \(\tau\) affect \(r=\tau x/(1+\tau x)\)? State the derivative with respect to \(x\), including the condition for a strictly positive derivative.
- Why can a partial agonist produce a substantial signal in one tissue but a smaller signal in another at the same occupancy?
- In a circuit model, one subtype has a negative weight. What does that sign mean, and what does it not establish?
- Why is a logistic output between 0 and 1 not automatically a probability?
Solutions. (1) \(\theta=4.0/(4.0+4.0)=0.500\). (2) For \(x>0\), increasing \(\tau\) increases \(r\), since \(\partial r/\partial\tau=x/(1+\tau x)^2>0\). Also, \(dr/dx=\tau/(1+\tau x)^2\geq0\) for \(\tau\geq0\), and it is strictly positive only when \(\tau>0\). (3) The tissues may have different transduction gains or downstream coupling, represented here by different \(\tau\); efficacy and occupancy alone do not determine the signal. (4) A negative weight means that the subtype signal reduces the specified model input. It does not mean the receptor is always inhibitory in every cell type or circuit. (5) A bounded transformation gives values in that range, but interpreting them as probabilities requires a model and endpoint that define them as probabilities.