Computational modelling identifies key determinants of subregion-specific dopamine dynamics in the striatum

Aske Ejdrup, Jakob Kisbye Dreyer, Matthew D Lycas, Søren H Jørgensen, Trevor W Robbins, Jeffrey Dalley, Freja Herborg, Ulrik Gether · eLife
L3 · Re-executed scripts Ran deposited figure-generation scripts with matplotlib patch. Vmax parameter sweep timed out (600s) — 3 simulation claims unverified.

Hypotheses The paper tests whether regional differences in striatal dopamine dynamics arise from release-related phenomena or from differential dopamine-transporter activity, and whether the classical phasic/tonic dichotomy corresponds to distinct receptor-level readouts by D1 and D2 receptors. It further asks whether DAT nanoclustering is sufficient to generate the observed dorsoventral gradient.

Claims A large-scale three-dimensional simulation of extracellular dopamine reproduces isolated hotspots without tonic baseline in dorsal striatum and a pervasive tonic-like signal in ventral striatum at 4 Hz pacemaker activity; parameter sweeps identify Vmax as the sole parameter generating the regional difference, consistent with immunostaining showing elevated DAT in DS and dSTORM showing greater DAT nanoclustering in VS. D1 receptor occupancy tracks phasic dopamine with ~50 ms delay, while D2 receptor occupancy integrates over seconds and is insensitive to brief firing pauses.

Inferences The prevailing view that phasic and tonic dopamine represent distinct release modes is refuted: both emerge from a single release process filtered through region-specific DAT activity, with nanoclustering proposed as the upstream regulator. D2 receptors cannot resolve individual bursts or pauses on behavioral timescales, constraining theories that assign D2 signaling a role in rapid event-level computation.

25 claims 13 verified

Logic of the Claims

Summary
TESTS —

Hypotheses tested

H1 HYPOTHESIS

D1 and D2 receptors operate on distinct temporal scales — D1 tracks burst DA with millisecond delay; D2 integrates over seconds and cannot resolve brief pauses.

Tested by
H1.1 EMPIRICAL fig1H, fig1I ✓ verified

D1R occupancy closely tracks extracellular DA with approximately 50 ms delay during burst firing; D1R occupancy is negligible during pacemaker activity because the EC50 of 1000 nM far exceeds tonic [DA] of ~10 nM in DS, making burst events the effective threshold for D1R engagement.

H1.3 EMPIRICAL fig1H ✓ verified

D2R occupancy takes at least 5 s to return to baseline after a burst due to slow off-kinetics (k_off = 0.2 s⁻¹), making D2R incapable of temporally separating closely linked bursts; the directional conclusion is robust but absolute occupancy values are sensitive to the initialization assumption.

H1.2 EMPIRICAL fig1J ✓ verified

A complete 1 s pause in firing reduces D2R occupancy from approximately 0.55 to 0.45 only; this finding is robust across an order of magnitude of D2R affinity (2–20 nM).

H2 HYPOTHESIS

DAT nanoclustering is a possible regulator of effective transporter activity, with denser clusters lowering effective Vmax via diffusion-limited substrate access — proposed as one contributor to the regional Vmax difference.

Tested by
H2.2 EMPIRICAL fig4C, fig4D, fig4E, fig4F, fig4G

Dense DAT nanoclusters (20 nm diameter) take approximately 400 ms to clear a 100 nM DA bolus compared to approximately 200 ms for unclustered DAT, because local [DA] at the cluster surface drops near zero creating a diffusion-limited bottleneck.

H2.1 EMPIRICAL fig4M, fig4N

Super-resolution dSTORM imaging shows DAT is significantly more nanoclustered in VS than DS (p=0.012, Welch's two-sample t-test, n=12 DS, n=13 VS), consistent across cluster sizes 20–200 nm.

H3 HYPOTHESIS

Regional differences in striatal DA dynamics (DS hotspots vs VS pervasive tonic coverage) are driven principally by DAT Vmax differences, not by release-side parameters.

Tested by
H3.1 EMPIRICAL fig1D, fig1E, fig1F ✓ verified

During 4 Hz pacemaker activity, the dorsal striatum produces partially segregated DA hotspots with large fractions of the simulated volume devoid of DA — there is no pervasive tonic baseline.

H3.6 EMPIRICAL fig2A, fig2B, fig2C ✓ verified

With DAT Vmax reduced to 33% of DS and terminal density at 90%, VS produces a diffuse tonic-like DA level throughout the simulated volume rather than segregated hotspots during 4 Hz pacemaker activity.

H3.4 EMPIRICAL fig3B, fig3C, fig3E, fig3F, fig3G, fig3K, fig3L ✓ verified

Across parameter sweeps of active terminal fraction, quantal size, release probability, and firing rate, only DAT Vmax generates differential responses between DS and VS; all other parameters shift both regions proportionally without altering the regional contrast.

H3.3 EMPIRICAL fig3K ✓ verified

A ±50% change in DAT Vmax shifts tonic DA by 38 nM in VS but only 11 nM in DS, indicating VS operates closer to the Km saturation regime and is more sensitive to DAT modulation.

H3.5 EMPIRICAL fig3B

VS at 5% active terminals produces a spatial DA distribution resembling DS at 100% active terminals, demonstrating VS operates in a low-focality high-coverage regime while DS requires dense terminal participation for equivalent spatial reach.

H3.2 EMPIRICAL fig2E

Simulated FSCV responses to 10, 30, and 60 Hz stimulation closely replicate May & Wightman (1989): VS reaches considerably higher peak DA than DS at all three stimulation frequencies.

Scope

SCOPE fig2A (implied throughout) ✓ verified

The 3:1 DS:VS DAT Vmax ratio (DS = 6 µM·s⁻¹, VS = 2 µM·s⁻¹) is assumed from published literature rather than directly measured in this study; the immunostaining gradient (Figure 2—supplement 1) corroborates this assumption at the protein level but does not directly establish the functional Vmax ratio.

Dissociations

H1.1 EMPIRICAL fig1H, fig1I ✓ verified

D1R occupancy closely tracks extracellular DA with approximately 50 ms delay during burst firing; D1R occupancy is negligible during pacemaker activity because the EC50 of 1000 nM far exceeds tonic [DA] of ~10 nM in DS, making burst events the effective threshold for D1R engagement.

H1.3 EMPIRICAL fig1H ✓ verified

D2R occupancy takes at least 5 s to return to baseline after a burst due to slow off-kinetics (k_off = 0.2 s⁻¹), making D2R incapable of temporally separating closely linked bursts; the directional conclusion is robust but absolute occupancy values are sensitive to the initialization assumption.

H3.1 EMPIRICAL fig1D, fig1E, fig1F ✓ verified

During 4 Hz pacemaker activity, the dorsal striatum produces partially segregated DA hotspots with large fractions of the simulated volume devoid of DA — there is no pervasive tonic baseline.

H3.6 EMPIRICAL fig2A, fig2B, fig2C ✓ verified

With DAT Vmax reduced to 33% of DS and terminal density at 90%, VS produces a diffuse tonic-like DA level throughout the simulated volume rather than segregated hotspots during 4 Hz pacemaker activity.

Eliminations & validating controls

C1 CONTROL fig2—supplement 1A, fig2—supplement 1B, fig2—supplement 1C

VMAT2 immunostaining shows no significant dorsoventral gradient in striatum (p=0.0086, one-sided t-test, n=4 mice), ruling out differential release capacity as the primary explanation for regional DA differences.

Standalone empirical findings

E1 EMPIRICAL fig2G

D2R occupancy during pacemaker activity is approximately 0.8 in VS versus approximately 0.55 in DS, consistent with higher prevailing tonic DA in VS.

E2 EMPIRICAL fig2—supplement 1B, fig2—supplement 1C

DAT expression is significantly higher in dorsal than ventral striatum (p=0.0021, one-sided t-test, n=4 mice), corroborating the 3:1 Vmax ratio assumed in the model.

E3 EMPIRICAL fig1G ✓ verified

3 APs at 10 Hz generates no significant DA spillover outside the burst zone; 6 APs at 20 Hz and 12 APs at 40 Hz cause frequency-dependent spillover exposing 10× and 30× the burst volume to concentrations above 100 nM respectively.

E4 EMPIRICAL fig2D ✓ verified

Even the lowest DA concentration percentiles in VS exceed 10 nM during 4 Hz pacemaker activity.

Methodological warrants

M1 METHODOLOGICAL fig1H ✓ verified

D2 receptor occupancy is initialized at 0.4 in all receptor dynamics simulations without derivation from steady state; at the modeled EC50 of 7 nM and simulated tonic [DA] of ~10 nM in DS, equilibrium occupancy would be approximately 0.59. No sensitivity analysis over this initialization is reported.

Scope qualifiers

Sc1 SCOPE fig4C, fig4D, fig4E, fig4F ✓ verified

The nanoclustering simulations hold total DAT Vmax constant across clustered and unclustered conditions — the per-voxel rate is multiplied by a normalization factor so total integrated uptake capacity is identical; if DAT nanoclustering co-occurs with increased total DAT expression in biology, the clearance-slowing result would not hold.

Sc2 SCOPE fig4C, fig4D, fig4E, fig4F ✓ verified

The nanoclustering simulations (Figure 4C–F) operate in a standalone varicosity-scale model (1.8 × 1.8 µm domain, 0.02 µm voxels) architecturally separate from the tissue-level model used in Figures 1–3 (100 µm domain, 1 µm voxels); no formal coupling exists between the two models and no effective-Vmax output from the nanoclustering simulation feeds into the tissue simulation.

All claims (alphabetical)

Abstract mapped to claims

The paper's abstract is shown with each sentence linked to the claim(s) it represents in the dependency graph. Hover or click a sentence to highlight the corresponding claim cards. Below: what the graph contains that the abstract leaves out, and vice versa.

Abstract

1Striatal dopamine (DA) release regulates reward-related learning and motivation and is believed to consist of a short-lived phasic and continuous tonic component. 2Here, we build a large-scale three-dimensional model of extracellular DA dynamics in dorsal (DS) and ventral striatum (VS). 3The model predicts rapid dynamics in DS with little to no basal DA and slower dynamics in the VS enabling build-up of tonic DA levels. 4These regional differences do not reflect release-related phenomena but rather differential dopamine transporter (DAT) activity. 5Interestingly, our simulations posit DAT nanoclustering as a possible regulator of this activity. 6Receptor binding simulations show that D1 receptor occupancy follows extracellular DA concentration with milliseconds delay, while D2 receptors do not respond to brief pauses in firing but rather integrate DA signal over seconds. 7Summarised, our model distills recent experimental observations into a computational framework that challenges prevailing paradigms of striatal DA signalling.

[1]
background / framing — not a paper-specific claim
[2]
background / framing — not a paper-specific claim
[3]
combines multiple claims → H3.1 · During 4 Hz pacemaker activity, the dorsal striatum produces partially segregate, H3.6 · With DAT Vmax reduced to 33% of DS and terminal density at 90%, VS produces a di
[4]
synthesis across claims → H3 · Regional differences in striatal DA dynamics (DS hotspots vs VS pervasive tonic , H3.4 · Across parameter sweeps of active terminal fraction, quantal size, release proba, C1 · VMAT2 immunostaining shows no significant dorsoventral gradient in striatum (p=0
[5]
synthesis across claims → H2 · DAT nanoclustering is a possible regulator of effective transporter activity, wi, H2.2 · Dense DAT nanoclusters (20 nm diameter) take approximately 400 ms to clear a 100, H2.1 · Super-resolution dSTORM imaging shows DAT is significantly more nanoclustered in
[6]
combines multiple claims → H1 · D1 and D2 receptors operate on distinct temporal scales — D1 tracks burst DA wit, H1.1 · D1R occupancy closely tracks extracellular DA with approximately 50 ms delay dur, H1.3 · D2R occupancy takes at least 5 s to return to baseline after a burst due to slow, H1.2 · A complete 1 s pause in firing reduces D2R occupancy from approximately 0.55 to
[7]
background / framing — not a paper-specific claim
Claims in the graph not surfaced in the abstract
  • ds-vs-vmax-ratio-assumed fig2A (implied throughout)
    The 3:1 DS:VS DAT Vmax ratio (DS = 6 µM·s⁻¹, VS = 2 µM·s⁻¹) is assumed from published literature rather than directly measured in this study; the immunostaining gradient (Figure 2—supplement 1) corroborates this assumption at the protein level but does not directly establish the functional Vmax ratio.
  • M1 d2r-initialization-unjustified fig1H
    D2 receptor occupancy is initialized at 0.4 in all receptor dynamics simulations without derivation from steady state; at the modeled EC50 of 7 nM and simulated tonic [DA] of ~10 nM in DS, equilibrium occupancy would be approximately 0.59. No sensitivity analysis over this initialization is reported.
  • Sc2 nanoclustering-model-varicosity-scale fig4C, fig4D, fig4E, fig4F
    The nanoclustering simulations (Figure 4C–F) operate in a standalone varicosity-scale model (1.8 × 1.8 µm domain, 0.02 µm voxels) architecturally separate from the tissue-level model used in Figures 1–3 (100 µm domain, 1 µm voxels); no formal coupling exists between the two models and no effective-Vmax output from the nanoclustering simulation feeds into the tissue simulation.
  • Sc1 nanoclustering-constant-vmax-constraint fig4C, fig4D, fig4E, fig4F
    The nanoclustering simulations hold total DAT Vmax constant across clustered and unclustered conditions — the per-voxel rate is multiplied by a normalization factor so total integrated uptake capacity is identical; if DAT nanoclustering co-occurs with increased total DAT expression in biology, the clearance-slowing result would not hold.
  • E3 low-burst-no-spillover-high-burst-does fig1G
    3 APs at 10 Hz generates no significant DA spillover outside the burst zone; 6 APs at 20 Hz and 12 APs at 40 Hz cause frequency-dependent spillover exposing 10× and 30× the burst volume to concentrations above 100 nM respectively.
  • E1 d2r-occupancy-higher-in-vs fig2G
    D2R occupancy during pacemaker activity is approximately 0.8 in VS versus approximately 0.55 in DS, consistent with higher prevailing tonic DA in VS.
  • E4 vs-lowest-percentiles-above-10nm fig2D
    Even the lowest DA concentration percentiles in VS exceed 10 nM during 4 Hz pacemaker activity.
  • H3.2 fscv-matches-may-wightman-1989 fig2E
    Simulated FSCV responses to 10, 30, and 60 Hz stimulation closely replicate May & Wightman (1989): VS reaches considerably higher peak DA than DS at all three stimulation frequencies.
  • H3.3 vmax-modulation-larger-impact-in-vs fig3K
    A ±50% change in DAT Vmax shifts tonic DA by 38 nM in VS but only 11 nM in DS, indicating VS operates closer to the Km saturation regime and is more sensitive to DAT modulation.
  • H3.5 vs-low-active-fraction-resembles-ds-distribution fig3B
    VS at 5% active terminals produces a spatial DA distribution resembling DS at 100% active terminals, demonstrating VS operates in a low-focality high-coverage regime while DS requires dense terminal participation for equivalent spatial reach.
  • E2 dat-immunostaining-dorsoventral-gradient fig2—supplement 1B, fig2—supplement 1C
    DAT expression is significantly higher in dorsal than ventral striatum (p=0.0021, one-sided t-test, n=4 mice), corroborating the 3:1 Vmax ratio assumed in the model.

Argument from the graphv3

An LLM was given only this paper's enriched claim graph — claims, panel references, roles, and the relations between them — with no access to the abstract, paper prose, or any external context. It was asked to reconstruct the argument. Below: the reconstructed argument, then a comparison to the published abstract above. The two prompts that produced these texts are at the bottom, available for inspection.

Reconstructed argument

Using a tissue-scale stochastic simulation of striatal dopamine, this work argues that the contrast between dorsal and ventral striatal dopamine dynamics arises principally from differences in dopamine transporter (DAT) Vmax rather than from release-side parameters, and that this single uptake difference is sufficient to explain a cascade of downstream consequences for receptor signaling. The central argument proceeds by dissociation: at 4 Hz pacemaker firing, the dorsal striatum produces segregated DA hotspots with large unoccupied volumes, while the ventral striatum — modeled with DAT Vmax reduced to one third of DS — produces a pervasive tonic field throughout the same volume. A parameter sweep over active terminal fraction, quantal size, release probability, and firing rate establishes that only Vmax generates differential DS/VS responses, ruling out each release-side alternative in turn; this elimination is corroborated by an immunohistochemistry control showing a DAT dorsoventral gradient with no matched VMAT2 gradient, which both validates the assumed 3:1 Vmax ratio and disconfirms the differential-release explanation. Modeling FSCV responses and replicating May & Wightman (1989) provides predictive cross-validation. A second dissociation argues that D1 and D2 receptors operate on distinct temporal scales: D1R occupancy tracks burst DA with ~50 ms delay and is silent during pacemaker firing because tonic [DA] sits well below its EC50, whereas D2R integrates over seconds and cannot resolve a complete 1 s firing pause — though the authors flag that the D2R initialization is set without steady-state derivation, scoping the absolute occupancy claims (the directional conclusion is reported as robust). A third hypothesis proposes DAT nanoclustering as a candidate regulator of effective Vmax, supported by varicosity-scale simulations showing diffusion-limited clearance and dSTORM evidence of greater DAT clustering in VS; however, this strand is bounded by two explicit scope conditions — total Vmax is held constant in the cluster simulations, and the cluster model is architecturally uncoupled from the tissue model — and the empirical clustering test is flagged as failed:mismatch, leaving the link between molecular architecture and tissue-scale Vmax as a structurally proposed but not demonstrated bridge. The whole argument is qualified by the literature-derived assumption of the 3:1 DS:VS Vmax ratio.

Compared to the published abstract

The v3 synthesis and the published abstract converge on the headline finding — DS hotspots vs VS pervasive tonic, DAT Vmax as the regional driver rather than release-side parameters, D1 ms-tracking vs D2 second-integration, and DAT nanoclustering as a candidate Vmax regulator — but they differ sharply in inferential transparency. The synthesis leads with the argumentative form ("the central argument proceeds by dissociation") and explicitly names the moves: parameter-sweep elimination of release-side alternatives, VMAT2 control as disconfirmation of differential release, May & Wightman cross-replication as predictive validation, and the literature-derived 3:1 Vmax ratio as a global scope qualifier. The abstract collapses this scaffolding into bald conjunctions ("do not reflect release-related phenomena but rather differential DAT activity"), absorbing three distinct claims into one mechanistic clause. The synthesis also surfaces the architectural disconnect between the varicosity-scale nanoclustering model and the tissue-scale model and flags the dSTORM cross-region test as failed:mismatch — both invisible in the abstract, which presents nanoclustering with a single "possible" hedge. Conversely, the abstract adds programmatic framing absent from the graph ("challenges prevailing paradigms"), and asserts dSTORM-style empirical clustering only as a model posit. Figure 3 (the parameter sweep that does the elimination work) is essentially invisible in the abstract; it is the structural backbone of the synthesis.

Diagnostic divergence

The abstract's clause "These regional differences do not reflect release-related phenomena but rather differential DAT activity" is one sentence; in the claim graph it is an entire argument-by-elimination, comprising a parameter sweep that rules out four release-side alternatives (active-fraction, quantal size, release probability, firing rate) plus a VMAT2 immunostaining control that rules out differential vesicular release capacity, the latter doubling as validation of the assumed 3:1 Vmax ratio. The abstract preserves the conclusion but erases the eliminative reasoning and the control structure that warrant it — the reader cannot tell from the abstract that this is a disjunctive argument with explicit alternatives ruled out. This is the diagnostic case for what enriched edges recover that prose compresses: rules-out and validates relations are precisely the joints that make a finding defensible, and they are exactly what an abstract drops first.

On inferential structure recovered

The v3 synthesis recovers the paper's inferential architecture in a way the abstract cannot: three nested hypotheses with explicit evidential, control, and scope edges. Argument by dissociation (DS hotspots vs VS tonic; D1 ms-tracking vs D2 seconds-integration) is named and uses the dissociates-with structure rather than a flat conjunction. Argument by elimination is preserved with its alternatives enumerated. Validation edges (VMAT2 control, FSCV cross-replication with May & Wightman 1989) are surfaced as warrant, not as findings. Scope qualifiers — the inherited 3:1 ratio, the varicosity-vs-tissue architectural gap, the constant-Vmax constraint in the cluster simulations, and the unjustified D2R initialization — are explicitly placed as boundary conditions rather than absorbed into prose hedges. The result is that the synthesis tells the reader why each conclusion is licensed (and where it is brittle), whereas the abstract tells the reader only what the conclusions are. The abstract's single "possible regulator" hedge for nanoclustering happens to track the architectural caveat the graph carries, but does so by collapse rather than explication; the failed:mismatch status of the dSTORM cross-region test is invisible in the abstract entirely. Figure 3's role as eliminative backbone is the most striking loss.

▸ Show traceback (8 synthesis sentences)
  1. Central hypothesis: DAT Vmax explains DS/VS regional difference
  2. Hypothesis: D1/D2 temporal dissociation
  3. Hypothesis: nanoclustering regulates effective Vmax (proposed bridge)
  4. Argument by elimination — release-side alternatives
  5. Argument by elimination — vesicular release capacity
  6. Global scope qualifier
  7. Standalone empirical observations (not in a hypothesis loop)
  8. Methodological caveat
▸ Show structured comparison table (14 rows)
KindAbstract sentenceSynthesis sentenceCategoryNote
Generic field framing; no corresponding claim in the graph.
Architectural scaffolding; the graph captures only the inherited Vmax ratio and the varicosity-vs-tissue scale separation, not a tissue-scale model identity claim.
Maps to ds-lacks-pervasive-tonic-da and vs-maintains-pervasive-tonic-da — the central dissociation. Synthesis names it as a dissociation and as the empirical legs of hypothesis-vmax-explains-regional-difference; abstract presents it as a flat empirical pair.
Three claims compressed to one: hypothesis-vmax-explains-regional-difference + vmax-only-parameter-driving-regional-difference (parameter-sweep elimination over four release-side alternatives) + vmat2-gradient-absent (control disconfirming differential vesicular release capacity). Synthesis preserves the eliminative and control structure; abstract collapses to a bald 'not X but Y'.
The abstract calls nanoclustering a 'simulation posit' — a model claim. The graph and synthesis carry both the simulation leg (dat-nanoclustering-slows-clearance) and the dSTORM observation leg (dat-clustering-greater-in-vs, status failed:mismatch). The abstract under-states empirical content and silently elides the failed test status.
Maps to hypothesis-d1-d2-temporal-distinction with its three empirical legs (d1r-tracks-da-50ms-delay, d2r-integrates-over-seconds, d2r-insensitive-to-brief-pauses). Synthesis frames as dissociation; abstract as conjunction. Abstract drops the d2r-initialization-unjustified scope flag.
Self-positioning rhetoric beyond any claim in the graph.
Same headline as abstract sentences 3–4 but states the thesis at hypothesis level, not as flat empirical pair.
Same claims as abstract sentence 3, but the synthesis uses the dissociates-with edge to name the rhetorical move; the abstract collapses to conjunction.
The rules-out and validates edges are absent from the abstract. Figure 3 is essentially invisible in the abstract; here it carries the argument.
fscv-matches-may-wightman-1989 has no abstract counterpart.
Aligned with abstract sentence 6 on substance; the scope flag (d2r-initialization-unjustified) is synthesis-only.
Where the abstract says 'simulations posit ... possible regulator' the synthesis says 'structurally proposed but not demonstrated', surfaces dSTORM as observation, and exposes the failed cross-region test and the architectural gap. The abstract's single hedge tracks but does not articulate this.
ds-vs-vmax-ratio-assumed scopes ALL claims; the abstract carries no scope qualifier of any kind.
▸ Show synthesizer prompt
You are reconstructing the argument of a scientific paper from its decomposed claim structure.

You have only the claims and the relations between them. You do not have the paper's title, abstract, prose, authors, or interpretive framing. You see the claim sentences, the panels they're tied to, their argumentative role, and the structural relations between them.

The claim graph carries multiple kinds of relation, each representing a different argumentative move:

- **`requires`** — A depends on B being true. Mechanistic / hierarchical chain.
- **`entails` / `derived-from`** — Hypothesis → prediction. Deductive entailment.
- **`tests`** — Empirical claim → prediction it tests.
- **`supports` / `refutes`** — Empirical claim → hypothesis it supports or refutes. Abductive inference.
- **`rules-out`** — A's evidence eliminates an alternative. Argument by elimination.
- **`dissociates-with`** — A and B jointly establish a dissociation. Argument by contrast.
- **`validates`** — A is a control or sign-flip that strengthens B. Argument by disconfirmation.
- **`predicts` / `confirms`** — predictive validation across model and experiment.
- **`scopes`** — A is a boundary condition on B (or on all claims). Argument by qualified scope.
- **`interprets`** — A reframes empirical B through theoretical / literature lens. Argument by reframing — not derivation, but an act of mapping.
- **`enables-method`** — A is the methodological capability that warrants B's interpretability.

Each claim has a role: `hypothesis`, `prediction`, `empirical`, `synthesis`, `interpretation`, `methodological`, `control`, or `scope`.

Scientific argument typically combines three reasoning forms:
- **Deduction** — `entails`/`derived-from` edges.
- **Induction** — `requires`/`supports` edges.
- **Abduction** — `supports`/`refutes` from empirical back to hypothesis.

Your task: write a paragraph (200–400 words) articulating what this paper is arguing, derived from the structure alone, in the style of a scientific abstract.

Use the right rhetorical move for the right structural relation. Honor epistemic markers and roles. Don't add background framing or literature you don't have. Don't speculate beyond claims. The structure of the argument should be visible in the prose.

Output:
1. Synthesis paragraph
2. Traceback

Claim graph follows.

---

# Ejdrup et al. — Claim Graph (reconstructed)

_Total claims: 23_


## Hypotheses


### Hypothesis: `hypothesis-d1-d2-temporal-distinction` (panel hypothesis; epistemic hypothesis; status unknown)
    > D1 and D2 receptors operate on distinct temporal scales — D1 tracks burst DA with millisecond delay; D2 integrates over seconds and cannot resolve brief pauses.

  - **entails (predictions/observations)**:
    - `d1r-tracks-da-50ms-delay` (panel fig1H, fig1I; epistemic moderate; status verified)
    > D1R occupancy closely tracks extracellular DA with approximately 50 ms delay during burst firing; D1R occupancy is negligible during pacemaker activity because the EC50 of 1000 nM far exceeds tonic [DA] of ~10 nM in DS, making burst events the effective threshold for D1R engagement.
    - `d2r-integrates-over-seconds` (panel fig1H; epistemic weak; status verified)
    > D2R occupancy takes at least 5 s to return to baseline after a burst due to slow off-kinetics (k_off = 0.2 s⁻¹), making D2R incapable of temporally separating closely linked bursts; the directional conclusion is robust but absolute occupancy values are sensitive to the initialization assumption.
    - `d2r-insensitive-to-brief-pauses` (panel fig1J; epistemic weak; status verified)
    > A complete 1 s pause in firing reduces D2R occupancy from approximately 0.55 to 0.45 only; this finding is robust across an order of magnitude of D2R affinity (2–20 nM).
  - **tested by**:
    - `d1r-tracks-da-50ms-delay` (panel fig1H, fig1I; epistemic moderate; status verified)
    > D1R occupancy closely tracks extracellular DA with approximately 50 ms delay during burst firing; D1R occupancy is negligible during pacemaker activity because the EC50 of 1000 nM far exceeds tonic [DA] of ~10 nM in DS, making burst events the effective threshold for D1R engagement.
    - `d2r-insensitive-to-brief-pauses` (panel fig1J; epistemic weak; status verified)
    > A complete 1 s pause in firing reduces D2R occupancy from approximately 0.55 to 0.45 only; this finding is robust across an order of magnitude of D2R affinity (2–20 nM).
    - `d2r-integrates-over-seconds` (panel fig1H; epistemic weak; status verified)
    > D2R occupancy takes at least 5 s to return to baseline after a burst due to slow off-kinetics (k_off = 0.2 s⁻¹), making D2R incapable of temporally separating closely linked bursts; the directional conclusion is robust but absolute occupancy values are sensitive to the initialization assumption.

### Hypothesis: `hypothesis-nanoclustering-regulates-vmax` (panel hypothesis; epistemic hypothesis; status unknown)
    > DAT nanoclustering is a possible regulator of effective transporter activity, with denser clusters lowering effective Vmax via diffusion-limited substrate access — proposed as one contributor to the regional Vmax difference.

  - **entails (predictions/observations)**:
    - `dat-nanoclustering-slows-clearance` (panel fig4C, fig4D, fig4E, fig4F, fig4G; epistemic moderate; status partial:zenodo-data-downloaded)
    > Dense DAT nanoclusters (20 nm diameter) take approximately 400 ms to clear a 100 nM DA bolus compared to approximately 200 ms for unclustered DAT, because local [DA] at the cluster surface drops near zero creating a diffusion-limited bottleneck.
    - `dat-clustering-greater-in-vs` (panel fig4M, fig4N; epistemic moderate; status failed:mismatch)
    > Super-resolution dSTORM imaging shows DAT is significantly more nanoclustered in VS than DS (p=0.012, Welch's two-sample t-test, n=12 DS, n=13 VS), consistent across cluster sizes 20–200 nm.
  - **tested by**:
    - `dat-clustering-greater-in-vs` (panel fig4M, fig4N; epistemic moderate; status failed:mismatch)
    > Super-resolution dSTORM imaging shows DAT is significantly more nanoclustered in VS than DS (p=0.012, Welch's two-sample t-test, n=12 DS, n=13 VS), consistent across cluster sizes 20–200 nm.
    - `dat-nanoclustering-slows-clearance` (panel fig4C, fig4D, fig4E, fig4F, fig4G; epistemic moderate; status partial:zenodo-data-downloaded)
    > Dense DAT nanoclusters (20 nm diameter) take approximately 400 ms to clear a 100 nM DA bolus compared to approximately 200 ms for unclustered DAT, because local [DA] at the cluster surface drops near zero creating a diffusion-limited bottleneck.
  - **scoped by (boundary conditions)**:
    - `nanoclustering-constant-vmax-constraint` (panel fig4C, fig4D, fig4E, fig4F; epistemic moderate; status verified)
    > The nanoclustering simulations hold total DAT Vmax constant across clustered and unclustered conditions — the per-voxel rate is multiplied by a normalization factor so total integrated uptake capacity is identical; if DAT nanoclustering co-occurs with increased total DAT expression in biology, the clearance-slowing result would not hold.
    - `nanoclustering-model-varicosity-scale` (panel fig4C, fig4D, fig4E, fig4F; epistemic moderate; status verified)
    > The nanoclustering simulations (Figure 4C–F) operate in a standalone varicosity-scale model (1.8 × 1.8 µm domain, 0.02 µm voxels) architecturally separate from the tissue-level model used in Figures 1–3 (100 µm domain, 1 µm voxels); no formal coupling exists between the two models and no effective-Vmax output from the nanoclustering simulation feeds into the tissue simulation.

### Hypothesis: `hypothesis-vmax-explains-regional-difference` (panel hypothesis; epistemic hypothesis; status unknown)
    > Regional differences in striatal DA dynamics (DS hotspots vs VS pervasive tonic coverage) are driven principally by DAT Vmax differences, not by release-side parameters.

  - **entails (predictions/observations)**:
    - `ds-lacks-pervasive-tonic-da` (panel fig1D, fig1E, fig1F; epistemic moderate; status verified)
    > During 4 Hz pacemaker activity, the dorsal striatum produces partially segregated DA hotspots with large fractions of the simulated volume devoid of DA — there is no pervasive tonic baseline.
    - `vs-maintains-pervasive-tonic-da` (panel fig2A, fig2B, fig2C; epistemic moderate; status verified)
    > With DAT Vmax reduced to 33% of DS and terminal density at 90%, VS produces a diffuse tonic-like DA level throughout the simulated volume rather than segregated hotspots during 4 Hz pacemaker activity.
    - `vmax-only-parameter-driving-regional-difference` (panel fig3B, fig3C, fig3E, fig3F, fig3G, fig3K, fig3L; epistemic strong; status verified)
    > Across parameter sweeps of active terminal fraction, quantal size, release probability, and firing rate, only DAT Vmax generates differential responses between DS and VS; all other parameters shift both regions proportionally without altering the regional contrast.
    - `vmax-modulation-larger-impact-in-vs` (panel fig3K; epistemic moderate; status verified)
    > A ±50% change in DAT Vmax shifts tonic DA by 38 nM in VS but only 11 nM in DS, indicating VS operates closer to the Km saturation regime and is more sensitive to DAT modulation.
    - `vs-low-active-fraction-resembles-ds-distribution` (panel fig3B; epistemic moderate; status unverified:compute-infeasible)
    > VS at 5% active terminals produces a spatial DA distribution resembling DS at 100% active terminals, demonstrating VS operates in a low-focality high-coverage regime while DS requires dense terminal participation for equivalent spatial reach.
    - `fscv-matches-may-wightman-1989` (panel fig2E; epistemic moderate; status unverified:compute-infeasible)
    > Simulated FSCV responses to 10, 30, and 60 Hz stimulation closely replicate May & Wightman (1989): VS reaches considerably higher peak DA than DS at all three stimulation frequencies.
  - **tested by**:
    - `ds-lacks-pervasive-tonic-da` (panel fig1D, fig1E, fig1F; epistemic moderate; status verified)
    > During 4 Hz pacemaker activity, the dorsal striatum produces partially segregated DA hotspots with large fractions of the simulated volume devoid of DA — there is no pervasive tonic baseline.
    - `fscv-matches-may-wightman-1989` (panel fig2E; epistemic moderate; status unverified:compute-infeasible)
    > Simulated FSCV responses to 10, 30, and 60 Hz stimulation closely replicate May & Wightman (1989): VS reaches considerably higher peak DA than DS at all three stimulation frequencies.
    - `vmax-modulation-larger-impact-in-vs` (panel fig3K; epistemic moderate; status verified)
    > A ±50% change in DAT Vmax shifts tonic DA by 38 nM in VS but only 11 nM in DS, indicating VS operates closer to the Km saturation regime and is more sensitive to DAT modulation.
    - `vmax-only-parameter-driving-regional-difference` (panel fig3B, fig3C, fig3E, fig3F, fig3G, fig3K, fig3L; epistemic strong; status verified)
    > Across parameter sweeps of active terminal fraction, quantal size, release probability, and firing rate, only DAT Vmax generates differential responses between DS and VS; all other parameters shift both regions proportionally without altering the regional contrast.
    - `vs-low-active-fraction-resembles-ds-distribution` (panel fig3B; epistemic moderate; status unverified:compute-infeasible)
    > VS at 5% active terminals produces a spatial DA distribution resembling DS at 100% active terminals, demonstrating VS operates in a low-focality high-coverage regime while DS requires dense terminal participation for equivalent spatial reach.
    - `vs-maintains-pervasive-tonic-da` (panel fig2A, fig2B, fig2C; epistemic moderate; status verified)
    > With DAT Vmax reduced to 33% of DS and terminal density at 90%, VS produces a diffuse tonic-like DA level throughout the simulated volume rather than segregated hotspots during 4 Hz pacemaker activity.
  - **abductively supported by (additional)**:
    - `dat-immunostaining-dorsoventral-gradient` (panel fig2—supplement 1B, fig2—supplement 1C; epistemic moderate; status unverified:no-data)
    > DAT expression is significantly higher in dorsal than ventral striatum (p=0.0021, one-sided t-test, n=4 mice), corroborating the 3:1 Vmax ratio assumed in the model.
    - `vmat2-gradient-absent` (panel fig2—supplement 1A, fig2—supplement 1B, fig2—supplement 1C; epistemic moderate; status unverified:no-data)
    > VMAT2 immunostaining shows no significant dorsoventral gradient in striatum (p=0.0086, one-sided t-test, n=4 mice), ruling out differential release capacity as the primary explanation for regional DA differences.
  - **validated by (control/sign-flip)**:
    - `dat-immunostaining-dorsoventral-gradient` (panel fig2—supplement 1B, fig2—supplement 1C; epistemic moderate; status unverified:no-data)
    > DAT expression is significantly higher in dorsal than ventral striatum (p=0.0021, one-sided t-test, n=4 mice), corroborating the 3:1 Vmax ratio assumed in the model.
    - `vmat2-gradient-absent` (panel fig2—supplement 1A, fig2—supplement 1B, fig2—supplement 1C; epistemic moderate; status unverified:no-data)
    > VMAT2 immunostaining shows no significant dorsoventral gradient in striatum (p=0.0086, one-sided t-test, n=4 mice), ruling out differential release capacity as the primary explanation for regional DA differences.

## Eliminations / Controls


### Control: `vmat2-gradient-absent` (panel fig2—supplement 1A, fig2—supplement 1B, fig2—supplement 1C; epistemic moderate; status unverified:no-data)
    > VMAT2 immunostaining shows no significant dorsoventral gradient in striatum (p=0.0086, one-sided t-test, n=4 mice), ruling out differential release capacity as the primary explanation for regional DA differences.

  - **rules-out**:
    - differential VMAT2 expression / vesicular release capacity as the explanation for the DS/VS DA difference
  - **supports**:
    - `hypothesis-vmax-explains-regional-difference` (panel hypothesis; epistemic hypothesis; status unknown)
    > Regional differences in striatal DA dynamics (DS hotspots vs VS pervasive tonic coverage) are driven principally by DAT Vmax differences, not by release-side parameters.
  - **validates**:
    - `hypothesis-vmax-explains-regional-difference` (panel hypothesis; epistemic hypothesis; status unknown)
    > Regional differences in striatal DA dynamics (DS hotspots vs VS pervasive tonic coverage) are driven principally by DAT Vmax differences, not by release-side parameters.

## Methodological caveats


### Method note: `d2r-initialization-unjustified` (panel fig1H; epistemic weak; status verified)
    > D2 receptor occupancy is initialized at 0.4 in all receptor dynamics simulations without derivation from steady state; at the modeled EC50 of 7 nM and simulated tonic [DA] of ~10 nM in DS, equilibrium occupancy would be approximately 0.59. No sensitivity analysis over this initialization is reported.

  - **scopes / qualifies**:
    - `d2r-integrates-over-seconds` (panel fig1H; epistemic weak; status verified)
    > D2R occupancy takes at least 5 s to return to baseline after a burst due to slow off-kinetics (k_off = 0.2 s⁻¹), making D2R incapable of temporally separating closely linked bursts; the directional conclusion is robust but absolute occupancy values are sensitive to the initialization assumption.
    - `d2r-insensitive-to-brief-pauses` (panel fig1J; epistemic weak; status verified)
    > A complete 1 s pause in firing reduces D2R occupancy from approximately 0.55 to 0.45 only; this finding is robust across an order of magnitude of D2R affinity (2–20 nM).
    - `d2r-occupancy-higher-in-vs` (panel fig2G; epistemic weak; status unverified:compute-infeasible)
    > D2R occupancy during pacemaker activity is approximately 0.8 in VS versus approximately 0.55 in DS, consistent with higher prevailing tonic DA in VS.

## Scope qualifiers


### Scope: `ds-vs-vmax-ratio-assumed` (panel fig2A (implied throughout); epistemic moderate; status verified)
    > The 3:1 DS:VS DAT Vmax ratio (DS = 6 µM·s⁻¹, VS = 2 µM·s⁻¹) is assumed from published literature rather than directly measured in this study; the immunostaining gradient (Figure 2—supplement 1) corroborates this assumption at the protein level but does not directly establish the functional Vmax ratio.

  - **scopes**: ALL claims in graph

### Scope: `nanoclustering-constant-vmax-constraint` (panel fig4C, fig4D, fig4E, fig4F; epistemic moderate; status verified)
    > The nanoclustering simulations hold total DAT Vmax constant across clustered and unclustered conditions — the per-voxel rate is multiplied by a normalization factor so total integrated uptake capacity is identical; if DAT nanoclustering co-occurs with increased total DAT expression in biology, the clearance-slowing result would not hold.

  - **scopes**:
    - `dat-nanoclustering-slows-clearance` (panel fig4C, fig4D, fig4E, fig4F, fig4G; epistemic moderate; status partial:zenodo-data-downloaded)
    > Dense DAT nanoclusters (20 nm diameter) take approximately 400 ms to clear a 100 nM DA bolus compared to approximately 200 ms for unclustered DAT, because local [DA] at the cluster surface drops near zero creating a diffusion-limited bottleneck.
    - `hypothesis-nanoclustering-regulates-vmax` (panel hypothesis; epistemic hypothesis; status unknown)
    > DAT nanoclustering is a possible regulator of effective transporter activity, with denser clusters lowering effective Vmax via diffusion-limited substrate access — proposed as one contributor to the regional Vmax difference.

### Scope: `nanoclustering-model-varicosity-scale` (panel fig4C, fig4D, fig4E, fig4F; epistemic moderate; status verified)
    > The nanoclustering simulations (Figure 4C–F) operate in a standalone varicosity-scale model (1.8 × 1.8 µm domain, 0.02 µm voxels) architecturally separate from the tissue-level model used in Figures 1–3 (100 µm domain, 1 µm voxels); no formal coupling exists between the two models and no effective-Vmax output from the nanoclustering simulation feeds into the tissue simulation.

  - **scopes**:
    - `dat-nanoclustering-slows-clearance` (panel fig4C, fig4D, fig4E, fig4F, fig4G; epistemic moderate; status partial:zenodo-data-downloaded)
    > Dense DAT nanoclusters (20 nm diameter) take approximately 400 ms to clear a 100 nM DA bolus compared to approximately 200 ms for unclustered DAT, because local [DA] at the cluster surface drops near zero creating a diffusion-limited bottleneck.
    - `hypothesis-nanoclustering-regulates-vmax` (panel hypothesis; epistemic hypothesis; status unknown)
    > DAT nanoclustering is a possible regulator of effective transporter activity, with denser clusters lowering effective Vmax via diffusion-limited substrate access — proposed as one contributor to the regional Vmax difference.

## Standalone empirical claims (not in a hypothesis loop)


- `d2r-occupancy-higher-in-vs` (panel fig2G; epistemic weak; status unverified:compute-infeasible)
    > D2R occupancy during pacemaker activity is approximately 0.8 in VS versus approximately 0.55 in DS, consistent with higher prevailing tonic DA in VS.
  - requires: ['vs-maintains-pervasive-tonic-da', 'd2r-initialization-unjustified']

- `low-burst-no-spillover-high-burst-does` (panel fig1G; epistemic moderate; status verified)
    > 3 APs at 10 Hz generates no significant DA spillover outside the burst zone; 6 APs at 20 Hz and 12 APs at 40 Hz cause frequency-dependent spillover exposing 10× and 30× the burst volume to concentrations above 100 nM respectively.
  - requires: ['ds-lacks-pervasive-tonic-da']

- `vs-lowest-percentiles-above-10nm` (panel fig2D; epistemic moderate; status verified)
    > Even the lowest DA concentration percentiles in VS exceed 10 nM during 4 Hz pacemaker activity.
  - requires: ['vs-maintains-pervasive-tonic-da']

## Dissociations


- **d1r-tracks-da-50ms-delay ⊥ d2r-integrates-over-seconds**
    - A: D1R occupancy closely tracks extracellular DA with approximately 50 ms delay during burst firing; D1R occupancy is negligible during pacemaker activity because the EC50 of 1000 nM far exceeds tonic [DA] of ~10 nM in DS, making burst events the effective threshold for D1R engagement.
    - B: D2R occupancy takes at least 5 s to return to baseline after a burst due to slow off-kinetics (k_off = 0.2 s⁻¹), making D2R incapable of temporally separating closely linked bursts; the directional conclusion is robust but absolute occupancy values are sensitive to the initialization assumption.

- **ds-lacks-pervasive-tonic-da ⊥ vs-maintains-pervasive-tonic-da**
    - A: During 4 Hz pacemaker activity, the dorsal striatum produces partially segregated DA hotspots with large fractions of the simulated volume devoid of DA — there is no pervasive tonic baseline.
    - B: With DAT Vmax reduced to 33% of DS and terminal density at 90%, VS produces a diffuse tonic-like DA level throughout the simulated volume rather than segregated hotspots during 4 Hz pacemaker activity.

## Key mechanistic 'requires' dependencies

- `d1r-tracks-da-50ms-delay` requires: ['ds-lacks-pervasive-tonic-da']
- `d2r-insensitive-to-brief-pauses` requires: ['d2r-integrates-over-seconds', 'd2r-initialization-unjustified']
- `d2r-integrates-over-seconds` requires: ['ds-lacks-pervasive-tonic-da', 'd2r-initialization-unjustified']
- `d2r-occupancy-higher-in-vs` requires: ['vs-maintains-pervasive-tonic-da', 'd2r-initialization-unjustified']
- `dat-nanoclustering-slows-clearance` requires: ['nanoclustering-model-varicosity-scale', 'nanoclustering-constant-vmax-constraint']
- `ds-lacks-pervasive-tonic-da` requires: ['ds-vs-vmax-ratio-assumed']
- `fscv-matches-may-wightman-1989` requires: ['vs-maintains-pervasive-tonic-da', 'ds-vs-vmax-ratio-assumed']
- `low-burst-no-spillover-high-burst-does` requires: ['ds-lacks-pervasive-tonic-da']
- `vmax-modulation-larger-impact-in-vs` requires: ['vmax-only-parameter-driving-regional-difference', 'ds-vs-vmax-ratio-assumed']
- `vmax-only-parameter-driving-regional-difference` requires: ['vs-maintains-pervasive-tonic-da']
- `vs-low-active-fraction-resembles-ds-distribution` requires: ['vs-maintains-pervasive-tonic-da']
- `vs-lowest-percentiles-above-10nm` requires: ['vs-maintains-pervasive-tonic-da']
- `vs-maintains-pervasive-tonic-da` requires: ['ds-lacks-pervasive-tonic-da', 'ds-vs-vmax-ratio-assumed']
▸ Show comparator prompt
You are comparing two articulations of the same paper's argument: (a) the published abstract, and (b) a synthesis reconstructed from the paper's enriched claim graph by an agent with no access to the abstract. The claim graph carries explicit edges representing argumentative forms — hypothesis (`entails`), prediction (`derived-from`), test (`tests`), abductive support (`supports`/`refutes`), elimination (`rules-out`), dissociation (`dissociates-with`), control validation (`validates`), interpretation (`interprets`), methodological warrant (`enables-method`), and scope qualification (`scopes`).

Your job: identify what each surfaces, what each hides, and what the divergences tell us — particularly about the inferential structure of the argument and how it's compressed in the abstract.

Look for:
- Inferential structure surfaced in synthesis but flattened in abstract
- Claims surfaced in synthesis but absent from abstract — categorize
- Assertions in abstract but not in synthesis — categorize
- Claims framed differently
- Ordering and emphasis
- Treatment of scope and caveats

Output:
1. 200-word comparison paragraph
2. Structured table classifying each abstract sentence and each synthesized sentence
3. Single most diagnostic divergence
4. Note on inferential structure recovery

ABSTRACT:
Striatal dopamine (DA) release regulates reward-related learning and motivation and is believed to consist of a short-lived phasic and continuous tonic component. Here, we build a large-scale three-dimensional model of extracellular DA dynamics in dorsal (DS) and ventral striatum (VS). The model predicts rapid dynamics in DS with little to no basal DA and slower dynamics in the VS enabling build-up of tonic DA levels. These regional differences do not reflect release-related phenomena but rather differential dopamine transporter (DAT) activity. Interestingly, our simulations posit DAT nanoclustering as a possible regulator of this activity. Receptor binding simulations show that D1 receptor occupancy follows extracellular DA concentration with milliseconds delay, while D2 receptors do not respond to brief pauses in firing but rather integrate DA signal over seconds. Summarised, our model distills recent experimental observations into a computational framework that challenges prevailing paradigms of striatal DA signalling.

SYNTHESIS V3:
Using a tissue-scale stochastic simulation of striatal dopamine, this work argues that the contrast between dorsal and ventral striatal dopamine dynamics arises principally from differences in dopamine transporter (DAT) Vmax rather than from release-side parameters, and that this single uptake difference is sufficient to explain a cascade of downstream consequences for receptor signaling. The central argument proceeds by dissociation: at 4 Hz pacemaker firing, the dorsal striatum produces segregated DA hotspots with large unoccupied volumes, while the ventral striatum — modeled with DAT Vmax reduced to one third of DS — produces a pervasive tonic field throughout the same volume. A parameter sweep over active terminal fraction, quantal size, release probability, and firing rate establishes that only Vmax generates differential DS/VS responses, ruling out each release-side alternative in turn; this elimination is corroborated by an immunohistochemistry control showing a DAT dorsoventral gradient with no matched VMAT2 gradient, which both validates the assumed 3:1 Vmax ratio and disconfirms the differential-release explanation. Modeling FSCV responses and replicating May & Wightman (1989) provides predictive cross-validation. A second dissociation argues that D1 and D2 receptors operate on distinct temporal scales: D1R occupancy tracks burst DA with ~50 ms delay and is silent during pacemaker firing because tonic [DA] sits well below its EC50, whereas D2R integrates over seconds and cannot resolve a complete 1 s firing pause — though the authors flag that the D2R initialization is set without steady-state derivation, scoping the absolute occupancy claims (the directional conclusion is reported as robust). A third hypothesis proposes DAT nanoclustering as a candidate regulator of effective Vmax, supported by varicosity-scale simulations showing diffusion-limited clearance and dSTORM evidence of greater DAT clustering in VS; however, this strand is bounded by two explicit scope conditions — total Vmax is held constant in the cluster simulations, and the cluster model is architecturally uncoupled from the tissue model — and the empirical clustering test is flagged as failed:mismatch, leaving the link between molecular architecture and tissue-scale Vmax as a structurally proposed but not demonstrated bridge. The whole argument is qualified by the literature-derived assumption of the 3:1 DS:VS Vmax ratio.