pith. sign in

arxiv: 2512.10983 · v2 · pith:H5JXZ7TCnew · submitted 2025-12-04 · 🧬 q-bio.TO · q-bio.CB· q-bio.NC

Compartmental-reaction diffusion framework for microscale dynamics of extracellular serotonin in brain tissue

Merlin Pelz , Skirmantas Janusonis , Gregory Handy This is my paper

Pith reviewed 2026-05-21 18:34 UTC · model grok-4.3

classification 🧬 q-bio.TO q-bio.CBq-bio.NC
keywords serotonin dynamicsreaction-diffusion modelsvaricositiesmicrodomainsextracellular signalingperturbation theorybrain tissue
0
0 comments X

The pith

Varicosities form diffusively coupled microdomains that generate spatial serotonin reservoirs

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a mathematical framework for the microscale dynamics of extracellular serotonin by formulating a two-dimensional compartmental-reaction diffusion system. It applies strong localized perturbation theory to derive an asymptotically equivalent set of nonlinear integro-ODEs that retain diffusive coupling for efficient computation. Analysis of steady states, bounds via Jensen's inequality, and closed-form expressions for spike maxima and minima show how firing frequency, varicosity geometry, and uptake kinetics control serotonin levels. The results indicate that varicosities create local reservoirs, clarifying distinctions between local and volume transmission while generating testable predictions for imaging and reuptake inhibitor effects.

Core claim

The framework formulates serotonin signaling as a two-dimensional compartmental-reaction diffusion system and uses strong localized perturbation theory to obtain an asymptotically equivalent set of nonlinear integro-ODEs. Analysis of period-averaged steady states, application of Jensen's inequality for bounds, and derivation of closed-form spike maxima and minima reveal that varicosities form diffusively coupled microdomains capable of generating spatial serotonin reservoirs. This clarifies aspects of local versus volume transmission and provides predictions for high-resolution serotonin imaging and the actions of selective serotonin-reuptake inhibitors.

What carries the argument

Two-dimensional compartmental-reaction diffusion system reduced via strong localized perturbation theory to nonlinear integro-ODEs that preserve diffusive coupling

If this is right

  • Firing frequency, varicosity geometry, and uptake kinetics quantitatively shape extracellular serotonin concentrations.
  • Varicosities act as diffusively coupled microdomains that generate spatial serotonin reservoirs.
  • The framework distinguishes local from volume transmission of serotonin.
  • Predictions aid interpretation of high-resolution serotonin imaging data.
  • The model yields insights into microscale actions of selective serotonin-reuptake inhibitors.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The reduced integro-ODE system could enable efficient simulation of serotonin dynamics across larger tissue volumes to examine network integration effects.
  • Comparable compartmental reaction-diffusion reductions may apply to modeling other small-molecule neurotransmitters with similar release site distributions.
  • Experimental tests using optogenetic control of serotonergic firing combined with fast imaging could directly probe the predicted reservoir stability.

Load-bearing premise

The central derivations assume that strong localized perturbation theory applied to the two-dimensional compartmental-reaction diffusion system produces an asymptotically equivalent set of nonlinear integro-ODEs that preserve diffusive coupling.

What would settle it

High-resolution imaging measurements of extracellular serotonin concentration gradients around individual varicosities, compared against the model's predicted spike maxima, minima, and reservoir formation, would test whether the microdomains exist as described.

Figures

Figures reproduced from arXiv: 2512.10983 by Gregory Handy, Merlin Pelz, Skirmantas Janusonis.

Figure 1
Figure 1. Figure 1: A: A typical meshwork of serotonergic fibers in the mouse forebrain (the hippocampus region). Individual [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: A system in R 2 with five serotonergic fibers (axons) in magenta and varicosity neighborhoods in cyan (left) gets straightened (middle), viewing deviations from a straight line as perturbations. We then focus only on a 1D section (black dashed in middle; black on right) orthogonal to the fibers with equidistantly placed varicosities, in order to solely analyze the contribution of neighboring firing varicos… view at source ↗
Figure 3
Figure 3. Figure 3: Top row: three snapshots of the dynamic serotonin distribution at the first firing event with [PITH_FULL_IMAGE:figures/full_fig_p016_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: As for Figure [PITH_FULL_IMAGE:figures/full_fig_p017_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: As for Figure [PITH_FULL_IMAGE:figures/full_fig_p017_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Left: The dynamics of total concentrations in the five varicosity neighborhoods along the 1-D section of [PITH_FULL_IMAGE:figures/full_fig_p018_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: The difference of the true, though numerically computed, [PITH_FULL_IMAGE:figures/full_fig_p018_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: The running maxima of the serotonin concentrations at all five varicosities along the section in Figure [PITH_FULL_IMAGE:figures/full_fig_p019_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: The dynamics of total concentrations in the five varicosity neighborhoods along the 1-D section in 2-D [PITH_FULL_IMAGE:figures/full_fig_p020_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Top: The dimensional serotonin concentration [PITH_FULL_IMAGE:figures/full_fig_p021_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Same as for Figure [PITH_FULL_IMAGE:figures/full_fig_p022_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Top: The dimensional serotonin concentration [PITH_FULL_IMAGE:figures/full_fig_p023_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: As in [48], the deformed Bromwich contour s = χ(1 − sin(α + z)) with discretization points placed at sl = χ(1 − sin(α + ilh)) and equidistant step h in the complex s-plane, zoomed out (left) and zoomed closer to s = 0 (right). In this way, only a few sl are needed for a quadrature for the exponentially decaying summands as Re(s) → −∞. Parameters: χ = 0.00089, α = 0.8, and h = 0.24706. Hence, the approxima… view at source ↗
read the original abstract

Serotonin (5-hydroxytryptamine) is a major neurotransmitter whose release from densely distributed serotonergic varicosities shapes plasticity and network integration throughout the brain, yet its extracellular dynamics remain poorly understood due to the sub-micrometer and millisecond scales involved. We develop a mathematical framework that captures the coupled reaction-diffusion processes governing serotonin signaling in realistic tissue microenvironments. Formulating a two-dimensional compartmental-reaction diffusion system, we use strong localized perturbation theory to derive an asymptotically equivalent set of nonlinear integro-ODEs that preserve diffusive coupling while enabling efficient computation. We analyze period-averaged steady states, establish bounds using Jensen's inequality, obtain closed-form spike maxima and minima, and implement a fast marching-scheme solver based on sum-of-exponentials kernels. These mathematical results provide quantitative insight into how firing frequency, varicosity geometry, and uptake kinetics shape extracellular serotonin. The model reveals that varicosities form diffusively coupled microdomains capable of generating spatial "serotonin reservoirs," clarifies aspects of local versus volume transmission, and yields predictions relevant to interpreting high-resolution serotonin imaging and the actions of selective serotonin-reuptake inhibitors.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 2 minor

Summary. The paper develops a two-dimensional compartmental-reaction diffusion model for extracellular serotonin dynamics around serotonergic varicosities. Using strong localized perturbation theory, it derives an asymptotically equivalent system of nonlinear integro-ODEs that retain diffusive coupling. The work then analyzes period-averaged steady states, derives Jensen bounds and closed-form expressions for spike maxima and minima, and implements a fast marching-scheme numerical solver based on sum-of-exponentials kernels. These results are used to examine how firing frequency, varicosity geometry, and uptake kinetics shape extracellular serotonin profiles, leading to the claim that varicosities form diffusively coupled microdomains that generate spatial serotonin reservoirs with implications for local versus volume transmission, high-resolution imaging, and SSRI pharmacology.

Significance. If the strong localized perturbation reduction is shown to be accurate for biologically plausible parameter regimes, the framework supplies an analytically tractable yet spatially resolved description of microscale serotonin dynamics. The closed-form extrema and Jensen bounds constitute a genuine strength, offering quantitative predictions that go beyond purely numerical studies and could directly inform interpretation of fast serotonin imaging data and the spatial effects of uptake inhibitors.

major comments (1)
  1. [Abstract (derivation of integro-ODEs)] The central claim that varicosities form diffusively coupled microdomains generating spatial serotonin reservoirs rests on the strong localized perturbation reduction from the 2D compartmental-reaction diffusion system to the nonlinear integro-ODEs (described in the abstract). This equivalence must hold for the subsequent period-averaged steady states, Jensen bounds, and closed-form spike extrema to be reliable. Explicit conditions on the localization strength, scale separation relative to varicosity radius and diffusion length, and uptake rates are required; without them the reduction may distort reservoir formation and the local-versus-volume transmission conclusions.
minor comments (2)
  1. [Abstract] The abstract states that uptake kinetics and varicosity geometry shape the results but does not indicate how these parameters are selected or calibrated against experimental data; a short statement on this choice would help readers evaluate the quantitative predictions.
  2. A direct numerical comparison between the reduced integro-ODE system and the original 2D PDE (e.g., in a supplementary figure) would strengthen in the approximation accuracy for the reported parameter ranges.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their detailed and constructive report. The major comment correctly identifies that the validity of the strong localized perturbation reduction underpins the subsequent analysis and biological conclusions. We address this point directly below and will incorporate the requested clarifications.

read point-by-point responses

  1. Referee: [Abstract (derivation of integro-ODEs)] The central claim that varicosities form diffusively coupled microdomains generating spatial serotonin reservoirs rests on the strong localized perturbation reduction from the 2D compartmental-reaction diffusion system to the nonlinear integro-ODEs (described in the abstract). This equivalence must hold for the subsequent period-averaged steady states, Jensen bounds, and closed-form spike extrema to be reliable. Explicit conditions on the localization strength, scale separation relative to varicosity radius and diffusion length, and uptake rates are required; without them the reduction may distort reservoir formation and the local-versus-volume transmission conclusions.

    Authors: We agree that explicit conditions on the applicability of the reduction are necessary to support the claims about microdomains and reservoirs. The derivation in Section 2 already invokes the strong-localization regime (varicosity radius much smaller than the characteristic diffusion length) and assumes the perturbation parameter ε is small, but these are stated implicitly through the asymptotic analysis rather than as a set of explicit inequalities. In the revised manuscript we will add a new subsection (2.4) that states the precise conditions: (i) localization strength ε = a / L_D ≪ 1 where a is varicosity radius and L_D = √(D / k_u) is the uptake-diffusion length; (ii) scale separation between inter-varicosity distance and L_D; and (iii) bounds on the uptake rate k_u relative to the firing frequency so that the quasi-steady approximation remains valid. We will also include a short numerical comparison (new Figure S1) between the full 2D reaction-diffusion PDE and the reduced integro-ODE system for a representative biologically plausible parameter set (a = 0.5 μm, D = 0.8 μm²/ms, k_u = 0.1–1 ms⁻¹). These additions will make the domain of validity transparent and directly address the concern that the reduction could distort the reservoir and local-versus-volume transmission conclusions. revision: yes

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The paper formulates a standard two-dimensional compartmental-reaction diffusion system from first-principles diffusion and reaction kinetics, then applies strong localized perturbation theory to derive an asymptotically equivalent nonlinear integro-ODE system that retains diffusive coupling. This reduction is a conventional asymptotic technique rather than a self-definitional or fitted-input mapping; subsequent steps such as period-averaged steady states, Jensen bounds, and closed-form spike extrema are obtained directly from the derived equations. No load-bearing self-citations, uniqueness theorems imported from prior author work, or ansatzes smuggled via citation appear in the provided derivation outline. Uptake kinetics and geometry parameters are treated as external inputs for calibration, not as quantities that define the predictions by construction. The overall chain remains self-contained against external mathematical benchmarks.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The model rests on standard physical assumptions of diffusion and reaction kinetics together with geometric idealizations of brain tissue; no new physical entities are postulated, but several kinetic and geometric quantities must be supplied from external measurements or chosen by hand.

free parameters (2)
  • uptake rate constants
    Kinetic parameters governing serotonin reuptake are required to close the reaction terms and are expected to be taken from or fitted to experimental values.
  • varicosity spacing and size
    Geometric parameters defining compartment dimensions and varicosity placement are introduced to represent realistic tissue microenvironments.
axioms (2)
  • domain assumption Extracellular serotonin dynamics can be captured by a two-dimensional compartmental reaction-diffusion system.
    Invoked at the outset to formulate the governing equations.
  • domain assumption Strong localized perturbation theory yields an asymptotically equivalent set of nonlinear integro-ODEs that preserve diffusive coupling.
    Used to reduce the PDE system to a computationally tractable form.

pith-pipeline@v0.9.0 · 5742 in / 1476 out tokens · 63090 ms · 2026-05-21T18:34:20.934553+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

7 extracted references · 7 canonical work pages

  1. [1]

    Adell, P

    [1]A. Adell, P. Celada, M. Abellán, and F. Artigas,Origin and functional role of the extracellular serotonin in the midbrain raphe nuclei, Brain Research Reviews, 39 (2002), pp. 154–180,https://doi.org/10.1016/S0165-0173(02) 00182-0. [2]K. Allers and T. Sharp,Neurochemical and anatomical identification of fast- and slow-firing neurones in the rat dorsal r...

    work page doi:10.1016/s0165-0173(02 2002

  2. [2]

    [21]S. V. F araone, C. L. W ard, M. Boucher, R. Elbekai, and E. Brunner,Role of serotonin in psychiatric and somatic comorbidities of attention-deficit/hyperactivity disorder: A systematic literature review, Neuroscience & Biobehavioral Reviews, 176 (2025), p. 106275,https://doi.org/10.1016/j.neubiorev.2025.106275. [22]L. Franke, M. Schmidtmann, A. Riedl,...

    work page doi:10.1016/j.neubiorev.2025.106275 2025

  3. [3]

    Handy, S

    [31]G. Handy, S. D. Lawley, and A. Borisyuk,Receptor recharge time drastically reduces the number of captured particles, PLOS Computational Biology, 14 (2018), p. e1006015,https://doi.org/10.1371/journal.pcbi.1006015. 32 [32]K. Hashimoto, Y. Yamawaki, K. Yamaoka, T. Yoshida, K. Okada, W. Tan, M. Yamasaki, Y. Matsumoto- Makidono, R. Kubo, H. Nakayama, T. K...

    work page doi:10.1371/journal.pcbi.1006015 2018

  4. [4]

    [42]B. L. Jacobs and E. C. Azmitia,Structure and function of the brain serotonin system, Physiological reviews, 72 (1992), pp. 165–229. [43]B. L. Jacobs and E. C. Azmitia,Structure and function of the brain serotonin system, Physiological Reviews, 72 (1992), pp. 165–229,https://doi.org/10.1152/physrev.1992.72.1.165. [44]S. Janušonis,Serotonin dynamics in ...

    work page doi:10.1152/physrev.1992.72.1.165 1992

  5. [5]

    [63]K. C. Mays, J. H. Haiman, and S. Janušonis,An experimental platform for stochastic analyses of single serotonergic fibers in the mouse brain, Frontiers in Neuroscience, 17 (2023), p. 1241919. [64]K. Miyazaki, K. W. Miyazaki, and K. Doya,Activation of Dorsal Raphe Serotonin Neurons Underlies Waiting for Delayed Rewards, The Journal of Neuroscience, 31 ...

    work page doi:10.1523/jneurosci.3714-10 2023

  6. [6]

    Mlinar, A

    [65]B. Mlinar, A. Montalbano, L. Piszczek, C. Gross, and R. Corradetti,Firing Properties of Genetically Identified Dorsal Raphe Serotonergic Neurons in Brain Slices, Frontiers in Cellular Neuroscience, 10 (2016),https://doi.org/10. 3389/fncel.2016.00195. [66]L. L. Moroz, D. Y. Romanova, and A. B. Kohn,Neural versus alternative integrative systems: molecul...

    work page doi:10.1098/rstb.2019.0762 2016

  7. [7]

    Pelz and M

    34 [77]M. Pelz and M. J. W ard,The emergence of spatial patterns for compartmental reaction kinetics coupled by two bulk diffusing species with comparable diffusivities, Philosophical Transactions of the Royal Society A, 381 (2023), p. 20220089. [78]M. Pelz and M. J. W ard,Symmetry-breaking bifurcations for compartmental reaction kinetics coupled by two b...

    work page doi:10.1007/s10571-021-01064-9 2023