Scale-Dependent Emergence of Hindered Diffusion in the Brain Extracellular Space

Abstract

Diffusion in living tissues governs essential physiological processes and is well studied within cells. Yet how extracellular molecular motion emerges from the structural complexity of tissues remains unresolved. In the brain, molecules move extensively through the extracellular space (ECS) enabling key functions, with effective diffusivities reduced by factors of 2 to 5 relative to free solution. This slowing has traditionally been captured by the phenomenological concept of tortuosity, but tortuosity does not specify the microscopic mechanisms responsible for diffusion hindrance. Here we directly visualize three dimensional extracellular diffusion in brain tissue using ultrashort single walled carbon nanotubes as nearinfrared tracers, achieving nanometric spatial precision and video rate temporal resolution. We find that motion is locally Brownian and that transport does not require scale free stochastic dynamics. Instead, hindered diffusion emerges from a geometry controlled crossover: free diffusion at short length scales gives way to constrained transport beyond a characteristic structural scale of the ECS. Thus, tortuosity arises as an emergent, scale dependent property rather than an intrinsic material constant. Beyond its biological implications, this behavior places extracellular transport within the broader physics of diffusion in disordered media. Brain tissue acts as a natural realization of geometry constrained transport phenomena observed in porous materials and random obstacle systems, linking molecular motion in living matter to the general case of structurally heterogeneous environments.

Figures (8)

• Figure 1: Three-dimensional single-particle tracking in the hippocampal extracellular space. (a) Optical scheme for three-dimensional single-particle tracking based on a double-helix point-spread function (DH-PSF), including the theoretical phase mask and the experimentally retrieved phase profile. (b) Experimental (Exp) and phase-retrieved (PR) DH-PSFs of an individual ultrashort carbon nanotube at different axial positions. (c) Schematic of the hippocampal CA3 layers and pyramidal neuron morphology within organotypic hippocampal slices. (d) Representative three-dimensional trajectories of ultrashort carbon nanotubes diffusing in the extracellular space of the CA3 pyramidal layer. Trajectories are color-coded by time.
• Figure 2: Trajectory-level heterogeneity of extracellular transport in the CA3 pyramidal layer. (a) Time-averaged mean-squared displacements (tMSDs) of individual ultrashort carbon nanotube trajectories recorded in the CA3 pyramidal layer, shown on log--log axes. (b) Joint distribution of the anomalous diffusion exponent $\alpha$ and generalized diffusion coefficient $K_\alpha$, extracted from tMSDs at lag time $\tau = 10$. Marginal distributions at different timelags are shown along the axes.
• Figure 3: Scale-dependent transport dynamics of the fast subpopulation in the CA3 pyramidal layer. (a) Time--ensemble averaged mean-squared displacement (teMSD, black circles) together with individual time-averaged MSDs (gray lines). (b) Probability density function of rescaled displacements at the shortest lag ($\tau=1$), with Gaussian reference (solid line). Inset: non-Gaussian parameter $g$ as a function of lag time. (c) Displacement autocorrelation function $C(\xi)$ as a function of normalized lag $\xi$. Inset: scale-dependent effective exponent $\alpha_C$. (d) Distribution of turning angles $\theta$ between successive displacement for increasing lag times $\tau$. The black curve indicates the expected distribution of Brownian motion.
• Figure 4: Layer-dependent, scale-dependent extracellular transport in CA3. (a) Schematic of the laminar organization of the CA3 hippocampus, highlighting the pyramidal cell layer and the stratum radiatum. (b) Comparison of ensemble mean-squared displacements (MSD) for uCCNT trajectories in the pyramidal layer (black) and stratum radiatum (red). Thin lines show individual time-averaged MSDs and open circles show the time--ensemble averaged MSD. The top subpanel shows the difference $\Delta=\log_{10}(\mathrm{MSD}_{rad})-\log_{10}(\mathrm{MSD}_{pyr})$, highlighting the small but systematic divergence between layers at longer lag times. (c) Scale-dependent displacement-correlation exponent $\alpha_C$ extracted from velocity autocorrelation function fits as a function of lag time (points indicate mean values; error bars denote s.e.m.). (d) Asphericity $\mathcal{A}$ of particle trajectories as a function of lag time for the two layers; the dashed line indicates the isotropic Brownian reference.
• Figure 5: Dynamical signatures of the slow subpopulation in the CA3 pyramidal layer. (a) Time--ensemble averaged mean-squared displacement (teMSD, black circles) together with individual time-averaged MSDs (gray lines). (b)Probability density function of rescaled displacements at the shortest lag ($\tau=1$), with Gaussian reference (solid line). Inset: non-Gaussian parameter $g$ as a function of lag time. (c) Time-dependent effective diffusion coefficient $D_{te}$ as a function of lagtime.