It is now well established that cell interiors are significantly crowded by macromolecules, which impede diffusion and enhance binding rates. However, it is not fully appreciated that levels of crowding are heterogeneous, and can vary substantially between subcellular regions. In this article, starting from a microscopic model, we derive coupled nonlinear partial differential equations for the concentrations of two populations of large and small spherical particles with steric volume exclusion. By performing an expansion in the ratio of the particle sizes, we find that the diffusion of a small particle in the presence of large particles obeys an advection–diffusion equation, with a reduced diffusion coefficient and a velocity directed towards less crowded regions. The interplay between advection and diffusion leads to behaviour that differs significantly from Brownian diffusion. We show that biologically plausible distributions of macromolecules can lead to highly non-Gaussian probability densities for the small particle position, including asymmetrical and multimodal densities. We confirm all our results using hard-sphere Brownian dynamics simulations.

The reaction-diffusion master equation (RDME) is a standard modelling approach for understanding stochastic and spatial chemical kinetics. An inherent assumption is that molecules are point-like. Here, we introduce the excluded volume reaction-diffusion master equation (vRDME) which takes into account volume exclusion effects on stochastic kinetics due to a finite molecular radius. We obtain an exact closed form solution of the RDME and of the vRDME for a general chemical system in equilibrium conditions. The difference between the two solutions increases with the ratio of molecular diameter to the compartment length scale. We show that an increase in the fraction of excluded space can (i) lead to deviations from the classical inverse square root law for the noise-strength, (ii) flip the skewness of the probability distribution from right to left-skewed, (iii) shift the equilibrium of bimolecular reactions so that more product molecules are formed, and (iv) strongly modulate the Fano factors and coefficients of variation. These volume exclusion effects are found to be particularly pronounced for chemical species not involved in chemical conservation laws. Finally, we show that statistics obtained using the vRDME are in good agreement with those obtained from Brownian dynamics with excluded volume interactions.

The Turing instability paradigm is revisited in the context of a multispecies diffusion scheme derived from a self-consistent microscopic formulation. The analysis is developed with reference to the case of two species. These latter share the same spatial reservoir and experience a degree of mutual interference due to the competition for the available resources. Turing instability can set in for all ratios of the main diffusivities, also when the (isolated) activator diffuses faster then the (isolated) inhibitor. This conclusion, at odd with the conventional vision, is here exemplified for the Brusselator model and ultimately stems from having assumed a generalized model of multispecies diffusion, fully anchored to first principles, which also holds under crowded conditions.

Gene expression occurs in an environment in which both stochastic and diffusive effects are significant. Spatial stochastic simulations are computationally expensive compared with their deterministic counterparts, and hence little is currently known of the significance of intrinsic noise in a spatial setting. Starting from the reaction–diffusion master equation (RDME) describing stochastic reaction–diffusion processes, we here derive expressions for the approximate steady-state mean concentrations which are explicit functions of the dimensionality of space, rate constants and diffusion coefficients. The expressions have a simple closed form when the system consists of one effective species. These formulae show that, even for spatially homogeneous systems, mean concentrations can depend on diffusion coefficients: this contradicts the predictions of deterministic reaction–diffusion processes, thus highlighting the importance of intrinsic noise. We confirm our theory by comparison with stochastic simulations, using the RDME and Brownian dynamics, of two models of stochastic and spatial gene expression in single cells and tissues.