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

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.

💡 Research Summary

This paper presents a novel mathematical framework, the Compartmental-Reaction Diffusion (CRD) system, to model the microscale spatiotemporal dynamics of extracellular serotonin in brain tissue. Serotonin, released from varicosities along serotonergic fibers, plays a crucial role in neuroplasticity and is a target for antidepressants like SSRIs. However, its extracellular signaling occurs at sub-micrometer and millisecond scales, making direct experimental observation challenging.

The core model conceptualizes the tissue as a two-dimensional plane containing discrete varicosity sites. Around each varicosity, a well-mixed circular compartment is defined. The serotonin concentration within each compartment (μ_j) is governed by an Ordinary Differential Equation (ODE) incorporating Michaelis-Menten uptake kinetics and pulsed release modeled by a smooth Gaussian function. The extracellular space between compartments is described by a Partial Differential Equation (PDE) for serotonin diffusion and linear degradation. These compartments and the extracellular space are coupled through flux boundary conditions, forming the complete CRD system.

A major technical achievement is the application of strong localized perturbation theory to derive an asymptotically equivalent set of nonlinear Integro-Ordinary Differential Equations (Integro-ODEs). This is achieved by taking the limit of infinite permeability at compartment boundaries. This reformulation preserves the essential diffusive coupling between varicosities while replacing the computationally expensive PDE with a convolution integral involving a modified Green’s function (a metabolic kernel). This transformation drastically improves computational efficiency.

Using this framework, the authors perform several key analyses. They compute period-averaged steady-state concentrations under periodic firing and establish bounds for these averages using Jensen’s inequality. Closed-form solutions for the peak (maximum) and trough (minimum) concentrations following a release pulse are derived, explicitly showing dependence on release amount, uptake parameters, and firing frequency. For numerical implementation, a fast solver based on a sum-of-exponentials approximation of the kernel and a fast marching scheme is developed.

The model is applied to simulated networks of fibers with multiple varicosities. Results demonstrate that varicosities are not independent release sites but form diffusively coupled microdomains. The firing frequency nonlinearly controls the time-averaged serotonin level across the network. Spatial gradients emerge, and nearby varicosities can synchronize their concentration profiles via diffusion, potentially creating spatial “serotonin reservoirs.” The model also clarifies the continuum between purely local synaptic transmission and volume transmission, showing that the effective signaling range is dynamically regulated by the balance of release and reuptake kinetics.

In conclusion, this work provides a rigorous and computationally tractable mathematical foundation for understanding extracellular serotonin dynamics. It offers quantitative predictions for how factors like SSRI administration (which alters uptake parameters V_max and K_m), varicosity density, and firing patterns shape the extracellular serotonin landscape. This framework is poised to become an essential tool for interpreting data from modern high-resolution serotonin sensors and for elucidating the spatial mechanisms of serotonergic modulation in health and disease.