Diffusive lensing as a mechanism of intracellular transport and compartmentalization

eLife assessment

The authors discuss an effect, "diffusive lensing", by which particles would accumulate in high-viscosity regions – for instance in the intracellular medium. To obtain these results, the authors rely on agent-based simulations using custom rules performed with the Ito stochastic calculus convention. The "lensing effect" discussed is a direct consequence of the choice of the Ito convention without spurious drift which has been discussed before and its adequacy for the intracellular medium is insufficiently discussed and relatively doubtful. Consequently, the relevance of the presented results for biology remain unclear and based on incomplete evidence.

https://doi.org/10.7554/eLife.89794.3.sa0

Strength of evidence:

Incomplete: Main claims are only partially supported

Exceptional
Compelling
Convincing
Solid
Incomplete
Inadequate

During the peer-review process the editor and reviewers write an eLife Assessment that summarises the significance of the findings reported in the article (on a scale ranging from landmark to useful) and the strength of the evidence (on a scale ranging from exceptional to inadequate). Learn more about eLife Assessments

Abstract
Introduction
Results
Discussion
Methods
Appendix 1
Data availability
References
Article and author information
Metrics

Abstract

While inhomogeneous diffusivity has been identified as a ubiquitous feature of the cellular interior, its implications for particle mobility and concentration at different length scales remain largely unexplored. In this work, we use agent-based simulations of diffusion to investigate how heterogeneous diffusivity affects the movement and concentration of diffusing particles. We propose that a nonequilibrium mode of membrane-less compartmentalization arising from the convergence of diffusive trajectories into low-diffusive sinks, which we call ‘diffusive lensing,’ is relevant for living systems. Our work highlights the phenomenon of diffusive lensing as a potentially key driver of mesoscale dynamics in the cytoplasm, with possible far-reaching implications for biochemical processes.

Introduction

Diffusion is a fundamental phenomenon of transport at scales ranging from atoms to galaxies. In cells, diffusion of individual components occurs in a complex, crowded milieu (Ellis, 2001; Luby-Phelps, 1999; van den Berg et al., 2017) that is known to exhibit position-dependent diffusivity (Berret, 2016; Garner et al., 2023; Huang et al., 2022; McLaughlin et al., 2020; Śmigiel et al., 2022; Xiang et al., 2020). Diffusion can occur within or between cellular compartments, where concentrated components carry out chemical reactions. This rich interaction of diffusion and compartmentalization provides the context for cellular biochemistry. Diffusivity varies inversely with viscosity, a key biophysical parameter of the cytoplasm (Bausch et al., 1999; Hu et al., 2017) that dictates translational and rotational mobility of proteins and, by extension, possibly influences their activity (Huang et al., 2022; Lippincott-Schwartz et al., 2001; Pan et al., 2009). While diffusivity has been implicated in modulating or driving a range of cellular processes (Molines et al., 2022; Persson et al., 2020; Xie et al., 2022), the role of inhomogeneous diffusivity in shaping biochemistry by regulating biomolecular concentration and dynamics remains poorly understood. Observation of diverse instances of accumulation across scales motivates our search for uncovering how space-dependent diffusivity affects cell biology. The accumulation of small molecules within the nuclear pore, for instance, has been attributed to diffusion through a viscous region (Ma et al., 2012). At the macroscale, Chladni patterns are an example of particle concentration resulting from inhomogeneous stochastic transport coefficients (Grabec, 2017). The implications of inhomogeneous diffusivity as a nonequilibrium phenomenon occurring at time scales and length scales relevant to biology remain largely unexplored. Theoretically, more information is required to specify the problem than just the diffusion constant: different mathematical interpretations of the stochastic term in diffusion equations with a spatially inhomogeneous diffusion constant result in different physical predictions (see Appendix for more information). Interestingly, diverse mesoscale outcomes are also seen in the case of active biological matter (Bechinger et al., 2016; Needleman and Dogic, 2017; Yeomans, 2017), the density-dependent concentration of active Brownian particles (Cates and Tailleur, 2015) and size-dependent condensation kinetics in the case of C. elegans colony formation (Chen and Ferrell, 2021). While these phenomena focus on motile energy-expending tracers, here we emphasize the underlying space-dependency of a physical property characterizing diffusion. In particular, accumulation arising from inhomogeneous diffusivity may represent a novel mechanism of effective compartmentalization, a key activity for cells in regulating biochemical processes.

In this work, we employ agent-based modeling to explore how position-dependent diffusivity can affect the distribution of tracer particles. We show that under a set of assumptions that relate to the ambiguities intrinsic to modeling inhomogeneous diffusivity (see Appendix), transport due to a diffusivity gradient leads to particle trajectories being biased toward areas of lower diffusivity, leading to effective compartmentalization and the growth of concentration gradients; we call this effect ‘diffusive lensing,’ in non-quantitative analogy to the effects on light rays of media with inhomogeneous refractive index, including refraction and the formation of caustics. Analyzing particle trajectories, we show that diffusive lensing manifests differently from homogeneous diffusion at the emergent scale. We conclude that inhomogeneous diffusivity may have diverse implications for intracellular transport, from sequestering particles to modulating where and when higher-order processes such as clustering happen, in a way that is not predictable from equivalent homogeneous-diffusivity models and could affect biochemical reactions.

Results

Inhomogeneous diffusivity drives particle accumulation

We probed the effect of inhomogeneous diffusivity on particle concentration using agent-based modeling of particle dynamics (Figure 1—figure supplement 1A; see Methods). In our modeling, the expected macroscale behavior is dictated by the Itô interpretation of heterogeneous diffusion (see Appendix) (Volpe and Wehr, 2016). Our model was non-anticipatory in that for a modeled particle traversing an inhomogeneous diffusivity, the step size distribution was defined by the diffusivity at its present position. Other equally consistent interpretations (such as the entirely anticipatory ‘isothermal’ interpretation) produce different macroscale behaviors (Figure 1—figure supplement 1B). The range of physically incompatible possibilities resulting from different interpretations is known as the Itô-Stratonovich dilemma (Lau and Lubensky, 2007; Sokolov, 2010; Tupper and Yang, 2012; Van Kampen, 1988; Volpe and Wehr, 2016). For systems at thermal equilibrium, the isothermal convention best describes transport; however, the non-equilibrium nature of the cellular interior motivates the consideration of non-isothermal conventions; the physically appropriate convention to use depends upon microscopic parameters and timescale hierarchies not captured in a coarse-grained model of diffusion. Note that while the Itô interpretation is deployed here, it is possible to convert from one interpretation to another (Volpe and Wehr, 2016) resulting in different interpretations converging at the same physical outcome (see Appendix). The equation used here is distinguished from the conventional 1D diffusion equation which arises from Fick’s laws and is only unambiguously true for homogeneous diffusion (characterized by constant diffusivity).

Over the course of the simulation, particles accumulated in the low-diffusivity zone (Figure 1A and C), consistent with steady-state closed form Itô-convention solutions (Tupper and Yang, 2012). This accumulation entailed the transient depletion of particles on the high-diffusive side of the interface. A similar accumulation was observed in a smooth diffusivity gradient (Figure 1B and D). In both cases, the results from agent-based modeling were corroborated by predictions of the steady-state analytical forms derived from theory. Thus, agent-based simulations demonstrate that under the Itô convention, areas of decreased diffusivity lead to increases in the concentration of diffusing particles. We term this phenomenon ‘diffusive lensing’.

Figure 1 with 1 supplement see all

Download asset Open asset

Low diffusivity leads to accumulation of particles.

(A) Particle distribution at various timesteps of a simulation with a step-like lower-diffusivity region. (B) Particle distribution at various timesteps for a simulation with a diffusivity gradient. (C) Steady-state particle distribution for the simulation in (A). (D) Steady-state particle distribution for the simulation in (B).

Interaction-mediated clustering is affected by heterogenous diffusivity

Diffusive lensing is an interaction-free mode of concentrating particles that stands in contrast to a more typical paradigm of particle accumulation: interaction-driven formation of higher-order structures like protein complexes, gels, crystals, and phase-separated condensates (Banani et al., 2017; Vekilov, 2010; Wu et al., 2023). How might interaction-induced clustering be modulated by inhomogeneous diffusion in a cellular context? To address this question, we heuristically modeled inter-particle interactions via a neighbor-sensing scheme in high and low interaction-strength regimes. The scheme involved using a step size for the modeled particle, which decreases as the number of particles in the vicinity increases (see Methods). At low interaction strength, clustering occurred only at the low-diffusivity end of a gradient (Figure 2A), while the same interaction strength was insufficient to produce clusters in a uniform diffusivity distribution (Figure 2—figure supplement 1A and C). In contrast, a high interaction strength resulted in robust clustering manifesting before particle gradient formation reached the steady-state, leading to clustering towards the high-diffusivity side of the simulation region as well (Figure 2B). At this high interaction strength, the clustering rate remained the same throughout the region in the absence of a gradient (Figure 2—figure supplement 1B and D). Taken together, the results reveal that diffusive lensing can modulate clustering and under certain circumstances cause diffusivity-dependent localized cluster formation, and furthermore that the relative strengths and timescales of each phenomenon quantitatively dictate whether increased clustering will preferentially occur in low-diffusive zones. Similar density-dependent clustering is observed in the case of active Brownian particles during motility-induced phase separation (Cates and Tailleur, 2015). Effects of diffusive lensing on particle concentration may additionally regulate reaction rates and drive stochastic clustering of enzymes (Jilkine et al., 2011).

Figure 2 with 1 supplement see all

Download asset Open asset

Interaction-driven clustering is modulated by heterogenous diffusivity.

(A) Progress of a simulation comprising particles possessing weak interactions ( $k = 0.04$ is the interaction strength; see *Methods*), initialized with a uniform concentration of particles. (B) Progress of a simulation comprising particles possessing strong interactions ( $k = 0.1$ ), initialized with a uniform concentration of particles.

Heterogeneous diffusion alters bulk particle motion as measured by in silico microrheology

The diffusion coefficient is a fundamental biophysical parameter that affects numerous other phenomena, including biochemical reaction rates. To elucidate particle diffusion at the microscale in the context of diffusive lensing, we used an in silico implementation of microrheology to analyze particle trajectories (see Methods; Figure 3—figure supplement 1A). We computed the mean squared displacements (MSDs) for uniform diffusivity simulations (in the case of unencumbered and confined diffusion) and used these to understand how MSD is affected by heterogenous diffusivity in two cases: a continuous diffusivity gradient and a discrete step in diffusivity.

Particle diffusion was unencumbered in the case of large bounds (relative to step size) (Figure 3A) and confined in the case of small bounds (Figure 3B) all in agreement with earlier results (Dix and Verkman, 2008; Saxton, 2007). The MSD at saturation in homogeneously diffusive systems was found to be agnostic to the underlying uniform diffusivity of the system, indicating that it is exclusively determined by the simulation region size. In contrast, particles in a diffusivity gradient exhibited dynamics intermediate to those of homogeneous high and low diffusivity cases, both in the diffusion coefficient and saturation MSD (Figure 3C, inset). The lowering of the saturation MSD reflects particle diffusion occurring within apparent simulation region bounds that confine more than the actual simulation region size. We note that such effective modifications of geometry are also a general feature of optical lensing. Apparent bounds were also found to occur in the two-zone diffusivity case (as in Figure 1A) where, at steady-state, particles populated the simulation region non-uniformly (Figure 3—figure supplement 1B). For most of the diffusivity ratio parameter space, irrespective of whether the smaller zones were more or less diffusive relative to the bulk, a reduction in MSD was seen indicating effectively lower diffusion bounds (Figure 3D). The magnitude of reduction depended on whether most particles resided in the larger or smaller of the two zones. In one observed case ( $\frac{μ_{i}}{μ_{o}} = 0.25$ ), however, the saturation MSD was higher than what was seen in the homogeneous diffusion scenario possibly due to particles robustly populating the bulk milieu followed by directed motion into the low-diffusive zone. The saturation MSD was also found to depend on the location of the low-diffusive zone: a more-centered zone resulted in a lowered saturation value, possibly due to weaker ratchet effects (Figure 3—figure supplement 1C and D). These results point to the insufficiency of using the diffusion coefficient alone to describe diffusion in heterogenous milieu. They also indicate a potentially rich interplay between heterogenous diffusivity and anomalous diffusion that requires further investigation.

Figure 3 with 1 supplement see all

Download asset Open asset

Heterogeneous diffusion alters bulk particle motion as measured by in silico microrheology.

(A) mean squared displacement (MSD) versus time for homogeneous diffusion of 10,000 particles in a 5 mm × 5 mm simulation region. (B) Same as (A) for homogeneous diffusion in a more tightly bounded simulation region (1 μm × 0.45 μm). (C) MSD versus time for inhomogeneous diffusion in a diffusivity gradient versus homogeneous diffusion in the extreme diffusivity cases (simulation region size: 1 μm × 0.45 μm). Inset: zoomed region showing differential saturation of the MSD. (D) MSD versus time for inhomogeneous diffusion due to a stepwise diffusivity distribution with diffusivity ratio $\frac{μ_{i}}{μ_{o}}$ relative to the bulk (simulation region size: 1 μm × 0.45 μm). In all cases, n=10,000 particles for MSD calculation (error bars denote SEM).

In silico FRAP in heterogeneously diffusive environments reveals drivers of mesoscale dynamics

The in silico microrheology analysis we performed provided insights into dynamics at the single-particle level (i.e. the microscale). To explore collective, emergent behaviors at the mesoscale while continuing to make contact with feasible experiments, we employed an in silico version of fluorescence recovery after photobleaching (in silico FRAP) (Figure 4—figure supplement 1A and B), in more cell-like inhomogeneous environments. In particular, we modeled low-diffusive patches/granules in a cell using a three-parameter disc-packing setup comprising granule radius ( $r$ ), packing density ( $ϕ$ ), and the ratio of granule diffusivity to bulk diffusivity ( $\frac{μ_{i}}{μ_{o}}$ ) (see Methods). We investigated the effect on dynamics of varying these parameters individually, with the goal of gaining an understanding of the effects of varying the amount, nature, and distribution of viscogens in cells. In all cases, the in silico ‘photobleaching’ event was conducted after the steady-state was attained (Figure 4—figure supplement 1C, D and E). To explain observed changes in the recovery time that would be measured in a FRAP-type experiment, we probed how the mean dwell time of particles in low-diffusive granules varies as a function of these parameters. A decrease in the diffusivity ratio ( $\frac{μ_{i}}{μ_{o}}$ ) at fixed $ϕ$ and $r$ resulted in a decline in measured particle mobility, as characterized by an increase in the simulated FRAP t_1/2 values (Figure 4A). Decreasing $\frac{μ_{i}}{μ_{o}}$ from 1 to 0.1 caused an approximate doubling of t_1/2 (or halving of diffusivity). Similar reduction in mobility was observed upon variation of $ϕ$ or $r$ separately, keeping the diffusivity ratio constant (Figure 4B and C). The decrease in average mobility in all three cases arose from changes in flux between the low-diffusive and bulk zones, as reflected by an increase in mean dwell times of particles within low-diffusive granules (Figure 4—figure supplement 1F, G and H). Furthermore, such reductions in mobility were emergent in that they arose from the interplay between granular diffusivity and bulk-granule fluxes, as the regions of interest in the simulated photobleaching events comprised granules and the surrounding bulk environment. To investigate whether particle dynamics is affected by the underlying topography realizing the system’s diffusivity, we averaged the granular and bulk diffusivity values to produce weighted-average diffusivity values, and compared in silico recovery in these simulations to that of the equivalent granule-comprising simulations. Such an averaging of the diffusivity to cause an effective uniform mobility for all resident particles resulted in slower dynamics than that of the equivalent granule-comprising simulations (Figure 4D). We conclude that inhomogeneity in diffusivity drives rapid effective dynamics via fluxes between the granular (interior) and bulk (exterior) environments, creating ‘diffusive highways’ for particles to move rapidly between low-diffusive regions. The diffusive lensing of particles into low-diffusive zones, and their consequent dwelling in these regions, can be tuned by modulating the underlying diffusivity distribution in myriad ways.

Figure 4 with 1 supplement see all

Download asset Open asset

A decrease in granule diffusivity, an increase in granule radius, or packing density slows down mesoscale dynamics.

(A) Simulated fluorescence recovery after photobleaching (FRAP) t_1/2 as a function of granule:bulk diffusivity ratio ( $r = 0.1 μ m$ , $ϕ = o .6$ ). (B) Simulated FRAP as a function of granule radius ( $\frac{μ_{i}}{μ_{o}} = 0.05$ , $ϕ = 0.6$ ). (C) Simulated FRAP t_1/2 as a function of granule packing density ( $\frac{μ_{i}}{μ_{o}} = 0.05$ , $ϕ = 0.01 μ m$ ). (D) Simulated FRAP for homogeneous and inhomogeneous diffusivity setups realizing the same effective diffusivities ( $\frac{μ_{i}}{μ_{o}} = 0.05$ , $ϕ = 0.01 μ m$ ). In all cases, n=3 ROIs were chosen for the simulated photobleaching (error bars denote SEM).

Discussion

The complex milieu of the cellular interior has been recently shown to feature heterogeneous diffusivity (Garner et al., 2023; Huang et al., 2022; McLaughlin et al., 2020; Parry et al., 2014; Śmigiel et al., 2022; Xiang et al., 2020), yet the consequences of such inhomogeneity on compartmentalization and mesoscale molecular dynamics have remained unclear. Through agent-based modeling of diffusion using the Itô integration convention, we show that heterogenous diffusivity can lead to simulated particle trajectories converging into low-diffusive hotspots, causing the accumulation of diffusing particles into membrane-less compartments defined by the lower-diffusivity zones. We term this mode of transport ‘diffusive lensing.’ The underlying conclusions from our 2D simulations extend to 3D directly (see Methods). Diffusive lensing has wide-ranging effects on particle distribution and dynamics and, furthermore, it can occur across a wide parameter space. We, therefore, speculate that diffusive lensing is a ubiquitous phenomenon in living systems.

We found that inhomogeneous diffusivity allows for particle mobility at the microscale and mesoscale to be different from that expected in the presence of homogeneous diffusion. Such an expectation is in line with predicted and observed deviations from normal diffusion in cells (Bancaud et al., 2012; Baum et al., 2014). The relative strengths of diffusive lensing and inter-particle interactions (if any) determined the extent to which clustering was modulated by diffusive lensing: this interplay may be important for determining the effects of inhomogeneous diffusivity on biochemical reaction rates. In these simulations of clustering, particle concentration did not affect diffusivity. In the case that particle concentration decreases diffusivity (for example in the case of branched polysaccharides like glycogen), diffusive lensing may create a positive feedback loop that drives particles into areas where low diffusivity has been nucleated. The effect of diffusive lensing on runaway pathological processes like protein aggregation is a potential direction for future work.

Spatially-averaged effective diffusion timescales were found to depend on the microscopic diffusivity distribution: the same average diffusivity can give rise to slower or faster dynamics depending on whether it is realized via homogeneous or heterogenous diffusivity distributions. In the latter case, the bulk region interspersed between the low-diffusive hotspots provides ‘diffusive highways’ that contribute to large fluxes at the diffusivity interface, thereby accounting for the faster dynamics. Such expressways and their associated fluxes may impact reaction kinetics by altering substrate turnover rates, congruent with the model of unusual transport processes potentially modifying reaction kinetics (Bénichou et al., 2010). In the context of subcellular low-diffusive regions (Garner et al., 2023), cells may compensate for geometry-imposed constraints on packing density and size of these regions by altering the diffusivity ratio (against the bulk milieu) instead. To map the detailed effects of inhomogeneous diffusivity on reaction rates, however, our work suggests that a key prerequisite is to chart a suitable set of meta parameters that provide an adequate description of inhomogeneous diffusion (Jin and Verkman, 2007), as a one-parameter description relying exclusively on the average diffusion coefficient is insufficient to fully specify the dynamics.

Changes in viscosity have been shown to occur in the context of cellular processes including cell death (Kuimova et al., 2008), stress adaptation (Persson et al., 2020) and protein aggregation (Thompson et al., 2015). At any given time point, intracellular transport dynamics arise emergently from contributions across length scales ranging from crowding in the bulk milieu due to proteins (Wang et al., 2010), and large biomolecules (Delarue et al., 2018) to cytoskeleton (Carlini et al., 2020; Chaubet et al., 2020) and active flows in the cytoplasm (Arcizet et al., 2008), all leading to unusual anomalous diffusive behaviors at the mesoscale (Banks and Fradin, 2005; Bressloff, 2014; Dix and Verkman, 2008; Höfling and Franosch, 2013; Kuznetsova et al., 2015; Swaminathan et al., 1997; Zhou et al., 2008). These diffusive behaviors cannot be decoupled from the intrinsic heterogeneity in biomolecular properties themselves (Heald and Cohen-Fix, 2014; Milo and Phillips, 2015). The effects of all of these subcellular determinants and energy-dependent processes on how position-dependent diffusivity is maintained in a cell remains unclear.

Not all cases of heterogeneous diffusivity will lead to diffusive lensing. This ambiguity is captured by the so-called Itô-Stratonovich dilemma (Lau and Lubensky, 2007; Sokolov, 2010; Tupper and Yang, 2012; Van Kampen, 1988; Volpe and Wehr, 2016). Any mathematical conceptualization of diffusion in the presence of position-dependent diffusivity must confront this dilemma, according to which the steady-state concentration distribution of a diffusing tracer depends not only on the localized diffusivity distribution but also on conventions based on microscopic parameters not captured in a coarse-grained model of diffusion; these parameters might, for example, include correlation lengths and times of viscogens or physical characteristics of polymers (Bo et al., 2021; Kupferman et al., 2004; Lau and Lubensky, 2007; Sokolov, 2010; Tupper and Yang, 2012; Van Kampen, 1988; Vishen et al., 2019). We speculate that any source of heterogeneity in diffusivity (including, but not limited to: mesh size experienced by the diffusing tracer, temperature changes, viscogen identity, and concentration) can, in turn, modulate diffusive lensing by means of altering either the particle or the environment-induced noise relaxation time. While the Itô convention is deployed here to model the nonequilibrium cellular interior (Gnesotto et al., 2018; Phillips et al., 2012), in some cases the isothermal convention may be better suited for modeling transport. The choice of the convention (and the effect of the dilemma, by extension) may also be subverted by recasting the dynamics into an alternate convention by taking suitable drift terms into consideration (see Appendix). Indeed, diverse conventions have been used to model experimentally observed accumulation arising from varied sources of such position-dependent noise (Bringuier, 2011; Pesce et al., 2013; Volpe and Wehr, 2016).

Our work underscores the need to not only examine diffusivity distributions in vivo as a function of local composition and the environment, but also to study their time evolution in response to external stimuli. More speculatively, we suggest that diffusive lensing serves as a potential candidate for a rudimentary mode of pre-biotic compartmentalization. Lensing-driven accumulation of diverse biomolecules may have served to produce chemically enriched spaces, acting as an antecedent of more sophisticated, membrane-bound, and membraneless organizational modalities; such a protocell organization is orthogonal to currently studied models (Monnard and Walde, 2015). This work demonstrates that diffusive lensing can have strong effects on transport and may be common in cellular contexts, modulating both passive and active flows. Future experimental and theoretical work will elucidate the extent of lensing inside and outside of cells and its effects on the biochemical reactions that sustain life.

Methods

Agent-based modeling (random walk simulations)

Agent-based modeling of diffusion was conducted via 2D random walk simulations. Irrespective of how the Itô-Stratonovich dilemma is interrogated, the underlying diffusion equations contain additive, separable contributions from each dimension, and this extends to 3D as well. Calculations were, therefore, carried out in 2D for simplicity and visualizability. Non-interacting point particles were initialized uniformly in a 2D simulation region with an aspect ratio matching that of an E. coli bacterium (Phillips et al., 2024). During each time step (also termed epoch or frame), every particle was moved along each coordinate by step sizes sampled from a uniform distribution, $U (- S, S)$ , where denotes the step size limit. Across a large number of steps, the distribution of displacements converges to the normal distribution by virtue of the central limit theorem. While sampling was not performed via the normal distribution directly by using the diffusion coefficient (D) as a parameter, the diffusion coefficient was instead arrived at as an emergent property of trajectories comprising a simulation, in a ground-up fashion. Reflection off the wall was modeled using a mirror-image rule. To model a zone of differential diffusivity relative to bulk diffusivity (either a fluid or a diffusivity zone), particle step sizes were sampled from zones characterized by different diffusivities, noting that the diffusion coefficient and diffusivity are inversely related Phillips et al., 2024 and $S \infty \sqrt{D}$ . At all times, step sizes were sampled from distributions defined by the diffusivity around the present position in accordance with the Itô interpretation of multiplicative noise (Volpe and Wehr, 2016) (for theoretical predictions of the steady-state behaviors, see Numerical methods for the diffusion equations). In all simulations, a set seed of 1 was used for the random number generator. Simulations were run on MATLAB R2020a on Sherlock (a high-performance computing cluster at Stanford).

In the simulations which included inter-particle interactions, these interactions were modeled via a neighbor-sensing approach. The step size limit was modified as per the relation, $S_{e f f} = S e^{-} k n$ , wheredenotes the sensing strength and denotes the number of neighbors (defined as those particles lying within a cutoff span around the particle in question). Such a rule-based approach modeled an effective attractive potential for the inter-particle interactions. Local density calculation used the same cutoff and the data were normalized to the mean local density of particles during initialization. Considering the computational work due to neighbor-sensing, a smaller number of particles (10³) were deployed, for a longer period of 2×10⁴ epochs.

In the low-diffusive granule simulations, the granules were modeled as disks with randomly initialized centers and fixed radii (r), covering the simulation region up to a desired packing density, $ϕ$ . The algorithm saturated for $ϕ \geq 0.6$ , in which case, the disks were generated as per cubic close packing, and their positions were incrementally deviated over steps to reduce local ordering as much as possible. The ratio of diffusivity inside the granules to diffusivity outside the granules ( $\frac{μ_{i}}{μ_{0}}$ ) was the third parameter under consideration. No two disks were allowed to overlap and all disks were kept confined within the boundaries of the simulation region. The default setup is as follows:r = 0.01 μm (uniform), $ϕ$ =0.6 (that is, 60% of the simulation region is covered by the granules) and $\frac{μ_{i}}{μ_{0}}$ = 0.05. Titration of one of these three parameters involved keeping the other two at the specified levels.

Numerical methods for the diffusion equations

The Fokker-Planck equations corresponding to the Ito, Stratonovich, and isothermal interpretations of inhomogeneous diffusion are as follows Gardiner, 2004; Tupper and Yang, 2012 (here $c (x, t)$ denotes the concentration distribution and $D (x)$ denotes the position-dependent diffusivity):

Itô interpretation: $\frac{\partial c}{\partial x} = \frac{\partial^{2} c D}{\partial x^{2}}$
Stratonovich interpretation: $\frac{\partial c}{\partial x} = \frac{\partial}{\partial x} (\sqrt{D} \frac{\partial c \sqrt{D}}{\partial x})$
Isothermal interpretation: $\frac{\partial}{\partial x} = \frac{\partial}{\partial x} (D \frac{\partial c}{\partial x})$

These equations were numerically evaluated via forward time-centered space (FTCS) schemes, with length and time increments set as 10^-3 and 5×10^-7 arbitrary units, respectively, and the number of time steps was set to 10⁵. A Gaussian well profile was used for the diffusion coefficient and the initial condition for the concentration distribution was a uniform distribution (Figure 1—figure supplement 1B). For the theoretical prediction in each case, the following relation is used: $c (x, t) D (x)^{1 - α} = c o n s t a n t$ in steady-state, where $α$ denotes the integration convention used (Tupper and Yang, 2012). Analysis and data visualization were performed on MATLAB R2019a.

In silico microrheology

Analysis of particle trajectories was carried out via quantifying the mean squared displacements (MSD). These were calculated from 10⁴ trajectories (each 10⁵ timesteps in duration) per simulation. The timestep was set as 50 μs so that the diffusion coefficient was $\approx 5 μ m^{2} / s$ order of magnitude for a small protein’s mobility in the E. coli cytoplasm (Milo and Phillips, 2015).

In silico FRAP

In silico fluorescence recovery after photobleaching (FRAP) studies were performed on the diffusion simulations to quantify emergent dynamics at the mesoscale. 10⁵ particles were deployed for a total duration of 0.5 s (10⁴ epochs). Circular regions (radius of 0.2 μm) were chosen as the regions of interest (ROIs). In silico photobleaching was instantaneously performed and involved assigning the particles in the ROI the photobleach status. The background was chosen from a uniform diffusivity setup to ensure that the normalization is standardized. The outward turnover of these particles and the simultaneous inward flux of unbleached particles were captured via t_1/2, the time taken for recovery up to 50% of the steady-state level of unbleached particles in the ROI (Sprague and McNally, 2005). In these simulations, t_1/2 connotes the time taken for the number of ‘unbleached’ particles in the ROI to reach 50% of the steady-state value. To dissect particles’ behavior during the simulation (in terms of bias towards inhabiting the low-diffusive granules), we calculated the mean dwell time across all particles, per simulation. This involved averaging the periods (of any duration) spent by particles inside low-diffusive granules. For normalization, the total simulation duration was used (0.5 s).

Appendix 1

Itô-Stratonovich dilemma and derivation of the generalized flux for heterogeneous diffusion

The Itô-Stratonovich dilemma is an ambiguity which arises directly from any analytical attempt to solve a stochastic differential equation comprising multiplicative (or position-dependent) noise (Amir, 2020; Gardiner, 2004; Van Kampen, 1988; Volpe and Wehr, 2016). The underlying mathematical reason for the ambiguity is the non-differentiable nature of the stochastic term, which causes different Riemann sum conventions to give quantitatively different results (Gardiner, 2004). From a more physical point of view, the different conventions correspond to different hierarchies of ‘small’ scales (correlation times and lengths of driving terms or viscogens, for example; see Pesce et al., 2013 for a discussion of competing timescales, Sokolov, 2010 for arriving at a context-dependent choice of interpretation). While the dilemma exists at the level of choosing an interpretation, it does not negate the appearance of diffusive lensing; indeed, only the magnitude of lensing is altered by the choice of the interpretation. This, too, may be subverted by adding an appropriate drift term to recast a stochastic differential equation (SDE) into one abiding by an alternate convention, as detailed later in this section. In summary, even the effect of the convention on the magnitude of diffusive lensing may be altered by recasting the SDE.

Formally, the dilemma can be boiled down to a choice of the parameter, $α (0 \leq α \leq 1)$ , in a diffusion equation like

\frac{\partial c}{\partial t} = \frac{- \partial J}{\partial x} = \frac{\partial}{\partial x} [(1 - α) c \frac{d D}{d x} + D \frac{\partial c}{\partial x}],

where $c (x, t)$ denotes the concentration distribution and $D (x)$ denotes the position-dependent diffusivity (see below for the derivation of the generalized flux). Clearly, for a position-dependent diffusivity, different values of $α$ will give rise to different physical predictions. In thermal equilibrium, $α$ must be equal to 1, but away from equilibrium it can take on different values depending on microscopic details of the physical system (Gardiner, 2004; Kupferman et al., 2004; Pesce et al., 2013; Smythe et al., 1983; Turelli, 1977; Vishen et al., 2019; Wang et al., 2020). Living cells are inherently removed from equilibrium by energy-driven processes breaking detailed balance (Gnesotto et al., 2018; Phillips et al., 2012). This motivates the calculations in this work, which use the Itô ( $α = 0$ ) integration convention, but in fact qualitatively similar results would be observed for any value of $α$ other than 1. Taken together with the ability to recast an SDE into one confining to an alternate convention, the dilemma does not take away from our principal thesis on how low diffusivity zones can accrete particles.

Consider a diffusing tracer whose trajectory is specified by the stochastic process, $X (t)$ . The stochastic process is continuous but is nowhere differentiable: the simplistic case of such a mathematical object is the Wiener process $W (t)$ , also known as Brownian motion (Gardiner, 2004). Note that in general, however, the stochastic process is characterized by a deterministic drift $a (x, t)$ and diffusion $b (x, t)$ . Indeed, $X (t)$ is in fact the solution of the stochastic differential equation (SDE):

d X (t) = a (x, t) d t + b (x, t) d W (t),

where $d W (t) = η (t) d t$ denotes the Wiener increment and $η (t)$ is the Gaussian white noise.

The stochastic differential equation here comprises multiplicative noise: that is, a position-dependent function, $b (x, t)$ , factors into the noise term of the SDE. This necessitates confronting the Itô-Stratonovich dilemma; we consider a general convention $α \in [0,1]$ , where $α = 0,0.5,1$ correspond to the Itô, Stratonovich, and isothermal conventions, respectively.

The SDE detailed above corresponds to a Fokker-Planck equation describing the time-evolution of the probability distribution of the tracer, $p (x, t)$ :

\frac{\partial p (x, t)}{\partial t} = \frac{- \partial}{\partial x} (a (x, t) p (x, t) + α b (x, t) \frac{\partial b (x, t)}{\partial x} p (x, t)) + \frac{1}{2} \frac{\partial^{2}}{\partial x^{2}} [b^{2} (x, t) p (x, t)]

(Substituting for $α$ yields specific cases; see Eqn. 4.3.20 in Gardiner, 2004, Eqns. 2.2.25, 2.2.26 in Bressloff, 2014, Table 2 and Discussion in Tupper and Yang, 2012.) Setting $a (x, t)$ (drift) to zero, and $b^{2} (x, t)$ (diffusion) to $2 D$ , a constant, yields Brownian motion, i.e., the canonical 1D diffusion equation which may also be derived from Fick’s laws. Keeping drift as zero, but setting diffusion to $2 D (x)$ (space-dependent) yields:

\frac{\partial p (x, t)}{\partial t} = \frac{- \partial}{\partial x} (α \frac{d D (x)}{d x} p (x, t)) + \frac{\partial^{2}}{\partial x^{2}} [D (x) p (x, t)]

which describes our scenarios in question. This equation can be simplified as follows:

\frac{\partial p (x, t)}{\partial t} = \frac{- \partial}{\partial x} [α \frac{d D (x)}{d x} p (x, t) - \frac{\partial}{\partial x} (D (x) p (x, t))]

\Rightarrow \frac{\partial p (x, t)}{\partial t} = \frac{- \partial}{\partial x} [- (1 - α) p (x, t) \frac{d D (x)}{d x} - D (x) \frac{\partial p (x, t)}{\partial x}] = \frac{- \partial}{\partial x} J (x, t),

where $J = - (1 - α) p (x, t) \frac{d D (x)}{d x} - D (x) \frac{\partial p (x, t)}{\partial x}$ is the generalized flux. Re-writing in terms of concentration yields the final expression, $J = - (1 - α) c \frac{d D}{d x} - D \frac{\partial c}{\partial x}$ . Substituting for $α$ and calculating $\frac{- \partial}{\partial x} J (x, t)$ yields the equations deployed in Methods: Numerical methods for the diffusion equations.

Note that it is also possible to convert between different interpretations of the multiplicative noise term by adding suitable drift terms as detailed in Volpe and Wehr, 2016; Tupper and Yang, 2012. In cases where thermal equilibrium is relevant, the Itô- and Stratonovich-conforming equations may be modified by adding drift terms to conform to that of the isothermal interpretation. Conversely, it is also possible to convert an equation based upon the isothermal convention to an Itô- or Stratonovich-conforming equation by modifying the drift term in the latter. Potential use cases as well as the consequences of following such conversions are detailed in these three example prescriptions below:

Appendix 1—table 1

Converting between stochastic integration conventions.

Data availability

The current manuscript is a computational study, so no data have been generated for this manuscript. Modelling code is uploaded as supplementary material and at https://zenodo.org/records/7957931.

The following data sets were generated

1. Achuthan RV
2. Kathy L
3. David W
4. Onn B
(2023) Zenodo
Diffusive Lensing Code.

https://doi.org/10.5281/zenodo.7957931

References

Book
1. Amir A
(2020) Thinking Probabilistically: Stochastic Processes, Disordered Systems, and Their Applications
Cambridge University Press.

https://doi.org/10.1017/9781108855259
- Google Scholar
(2008) Temporal analysis of active and passive transport in living cells
Physical Review Letters 101:248103.

https://doi.org/10.1103/PhysRevLett.101.248103
- PubMed
- Google Scholar
1. Banani SF
2. Lee HO
3. Hyman AA
4. Rosen MK
(2017) Biomolecular condensates: organizers of cellular biochemistry
Nature Reviews. Molecular Cell Biology 18:285–298.

https://doi.org/10.1038/nrm.2017.7
- PubMed
- Google Scholar
(2012) A fractal model for nuclear organization: current evidence and biological implications
Nucleic Acids Research 40:8783–8792.

https://doi.org/10.1093/nar/gks586
- PubMed
- Google Scholar
1. Banks DS
2. Fradin C
(2005) Anomalous diffusion of proteins due to molecular crowding
Biophysical Journal 89:2960–2971.

https://doi.org/10.1529/biophysj.104.051078
- PubMed
- Google Scholar
1. Baum M
2. Erdel F
3. Wachsmuth M
4. Rippe K
(2014) Retrieving the intracellular topology from multi-scale protein mobility mapping in living cells
Nature Communications 5:4494.

https://doi.org/10.1038/ncomms5494
- PubMed
- Google Scholar
(1999) Measurement of local viscoelasticity and forces in living cells by magnetic tweezers
Biophysical Journal 76:573–579.

https://doi.org/10.1016/S0006-3495(99)77225-5
- PubMed
- Google Scholar
(2016) Active particles in complex and crowded environments
Reviews of Modern Physics 88:045006.

https://doi.org/10.1103/RevModPhys.88.045006
- Google Scholar
(2010) Geometry-controlled kinetics
Nature Chemistry 2:472–477.

https://doi.org/10.1038/nchem.622
- PubMed
- Google Scholar
1. Berret JF
(2016) Local viscoelasticity of living cells measured by rotational magnetic spectroscopy
Nature Communications 7:10134.

https://doi.org/10.1038/ncomms10134
- PubMed
- Google Scholar
1. Bo S
2. Hubatsch L
3. Bauermann J
4. Weber CA
5. Jülicher F
(2021) Stochastic dynamics of single molecules across phase boundaries
Physical Review Research 3:043150.

https://doi.org/10.1103/PhysRevResearch.3.043150
- Google Scholar
Book
1. Bressloff PC
(2014) Stochastic Processes in Cell Biology
Springer International Publishing.

https://doi.org/10.1007/978-3-319-08488-6
- Google Scholar
1. Bringuier E
(2011) Particle diffusion in an inhomogeneous medium
European Journal of Physics 32:975–992.

https://doi.org/10.1088/0143-0807/32/4/012
- Google Scholar
(2020) Microtubules enhance mesoscale effective diffusivity in the crowded metaphase cytoplasm
Developmental Cell 54:574–582.

https://doi.org/10.1016/j.devcel.2020.07.020
- PubMed
- Google Scholar
1. Cates ME
2. Tailleur J
(2015) Motility-induced phase separation
Annual Review of Condensed Matter Physics 6:219–244.

https://doi.org/10.1146/annurev-conmatphys-031214-014710
- Google Scholar
(2020) Dynamic actin cross-linking governs the cytoplasm’s transition to fluid-like behavior
Molecular Biology of the Cell 31:1744–1752.

https://doi.org/10.1091/mbc.E19-09-0504
- PubMed
- Google Scholar
1. Chen Y
2. Ferrell JE
(2021) C. elegans colony formation as a condensation phenomenon
Nature Communications 12:4947.

https://doi.org/10.1038/s41467-021-25244-9
- Google Scholar
1. Delarue M
2. Brittingham GP
3. Pfeffer S
4. Surovtsev IV
5. Pinglay S
6. Kennedy KJ
7. Schaffer M
8. Gutierrez JI
9. Sang D
10. Poterewicz G
11. Chung JK
12. Plitzko JM
13. Groves JT
14. Jacobs-Wagner C
15. Engel BD
16. Holt LJ
(2018) mTORC1 controls phase separation and the biophysical properties of the cytoplasm by tuning Crowding
Cell 174:338–349.

https://doi.org/10.1016/j.cell.2018.05.042
- Google Scholar
1. Dix JA
2. Verkman AS
(2008) Crowding effects on diffusion in solutions and cells
Annual Review of Biophysics 37:247–263.

https://doi.org/10.1146/annurev.biophys.37.032807.125824
- PubMed
- Google Scholar
1. Ellis RJ
(2001) Macromolecular crowding: obvious but underappreciated
Trends in Biochemical Sciences 26:597–604.

https://doi.org/10.1016/s0968-0004(01)01938-7
- PubMed
- Google Scholar
Book
1. Gardiner CW
(2004) Handbook of Stochastic Methods for Physics, Chemistry and the Natural Sciences
Springer.

https://doi.org/10.1007/978-3-662-05389-8
- Google Scholar
(2023) Vast heterogeneity in cytoplasmic diffusion rates revealed by nanorheology and Doppelgänger simulations
Biophysical Journal 122:767–783.

https://doi.org/10.1016/j.bpj.2023.01.040
- PubMed
- Google Scholar
(2018) Broken detailed balance and non-equilibrium dynamics in living systems: a review
Reports on Progress in Physics. Physical Society 81:066601.

https://doi.org/10.1088/1361-6633/aab3ed
- PubMed
- Google Scholar
1. Grabec I
(2017) Vibration driven random walk in a Chladni experiment
Physics Letters A 381:59–64.

https://doi.org/10.1016/j.physleta.2016.10.059
- Google Scholar
1. Heald R
2. Cohen-Fix O
(2014) Morphology and function of membrane-bound organelles
Current Opinion in Cell Biology 26:79–86.

https://doi.org/10.1016/j.ceb.2013.10.006
- PubMed
- Google Scholar
1. Höfling F
2. Franosch T
(2013) Anomalous transport in the crowded world of biological cells
Reports on Progress in Physics. Physical Society 76:046602.

https://doi.org/10.1088/0034-4885/76/4/046602
- PubMed
- Google Scholar
1. Hu J
2. Jafari S
3. Han Y
4. Grodzinsky AJ
5. Cai S
6. Guo M
(2017) Size- and speed-dependent mechanical behavior in living mammalian cytoplasm
PNAS 114:9529–9534.

https://doi.org/10.1073/pnas.1702488114
- PubMed
- Google Scholar
(2022) Cytoplasmic organization promotes protein diffusion in Xenopus extracts
Nature Communications 13:5599.

https://doi.org/10.1038/s41467-022-33339-0
- PubMed
- Google Scholar
(2011) A density-dependent switch drives stochastic clustering and polarization of signaling molecules
PLOS Computational Biology 7:e1002271.

https://doi.org/10.1371/journal.pcbi.1002271
- PubMed
- Google Scholar
1. Jin S
2. Verkman AS
(2007) Single particle tracking of complex diffusion in membranes: simulation and detection of barrier, raft, and interaction phenomena
The Journal of Physical Chemistry. B 111:3625–3632.

https://doi.org/10.1021/jp067187m
- PubMed
- Google Scholar
(2008) Molecular rotor measures viscosity of live cells via fluorescence lifetime imaging
Journal of the American Chemical Society 130:6672–6673.

https://doi.org/10.1021/ja800570d
- PubMed
- Google Scholar
(2004) Itô versus Stratonovich white-noise limits for systems with inertia and colored multiplicative noise
Physical Review E 70:036120.

https://doi.org/10.1103/PhysRevE.70.036120
- Google Scholar
(2015) Beyond the excluded volume effects: mechanistic complexity of the crowded milieu
Molecules 20:1377–1409.

https://doi.org/10.3390/molecules20011377
- PubMed
- Google Scholar
1. Lau AWC
2. Lubensky TC
(2007) State-dependent diffusion: thermodynamic consistency and its path integral formulation
Physical Review E 76:011123.

https://doi.org/10.1103/PhysRevE.76.011123
- Google Scholar
(2001) Studying protein dynamics in living cells
Nature Reviews. Molecular Cell Biology 2:444–456.

https://doi.org/10.1038/35073068
- PubMed
- Google Scholar
1. Luby-Phelps K
(1999) Cytoarchitecture and physical properties of cytoplasm: volume, viscosity, diffusion, intracellular surface area
International Review of Cytology 192:189–221.

https://doi.org/10.1016/s0074-7696(08)60527-6
- PubMed
- Google Scholar
1. Ma J
2. Goryaynov A
3. Sarma A
4. Yang W
(2012) Self-regulated viscous channel in the nuclear pore complex
PNAS 109:7326–7331.

https://doi.org/10.1073/pnas.1201724109
- PubMed
- Google Scholar
(2020) Spatial heterogeneity of the cytosol revealed by machine learning-based 3D particle tracking
Molecular Biology of the Cell 31:1498–1511.

https://doi.org/10.1091/mbc.E20-03-0210
- PubMed
- Google Scholar
Book
1. Milo R
2. Phillips R
(2015) Cell Biology by the Numbers
Garland Science.

https://doi.org/10.1201/9780429258770
- Google Scholar
1. Molines AT
2. Lemière J
3. Gazzola M
4. Steinmark IE
5. Edrington CH
6. Hsu CT
7. Real-Calderon P
8. Suhling K
9. Goshima G
10. Holt LJ
11. Thery M
12. Brouhard GJ
13. Chang F
(2022) Physical properties of the cytoplasm modulate the rates of microtubule polymerization and depolymerization
Developmental Cell 57:466–479.

https://doi.org/10.1016/j.devcel.2022.02.001
- PubMed
- Google Scholar
1. Monnard PA
2. Walde P
(2015) Current ideas about prebiological compartmentalization
Life 5:1239–1263.

https://doi.org/10.3390/life5021239
- PubMed
- Google Scholar
1. Needleman D
2. Dogic Z
(2017) Active matter at the interface between materials science and cell biology
Nature Reviews Materials 2:1–14.

https://doi.org/10.1038/natrevmats.2017.48
- Google Scholar
1. Pan W
2. Filobelo L
3. Pham NDQ
4. Galkin O
5. Uzunova VV
6. Vekilov PG
(2009) Viscoelasticity in homogeneous protein solutions
Physical Review Letters 102:058101.

https://doi.org/10.1103/PhysRevLett.102.058101
- Google Scholar
(2014) The bacterial cytoplasm has glass-like properties and is fluidized by metabolic activity
Cell 156:183–194.

https://doi.org/10.1016/j.cell.2013.11.028
- Google Scholar
(2020) Cellular control of viscosity counters changes in temperature and energy availability
Cell 183:1572–1585.

https://doi.org/10.1016/j.cell.2020.10.017
- PubMed
- Google Scholar
1. Pesce G
2. McDaniel A
3. Hottovy S
4. Wehr J
5. Volpe G
(2013) Stratonovich-to-Itô transition in noisy systems with multiplicative feedback
Nature Communications 4:2733.

https://doi.org/10.1038/ncomms3733
- PubMed
- Google Scholar
Book
1. Phillips R
2. Kondev J
3. Theriot J
4. Garcia HG
5. Orme N
(2012) Physical Biology of the Cell
Garland Science.

https://doi.org/10.1201/9781134111589
- Google Scholar
Book
1. Phillips R
2. Kondev J
3. Theriot J
4. Garcia HG
5. Orme N
(2024)
Physical Biology of the Cell

Routledge.
- Google Scholar
1. Saxton MJ
(2007) A biological interpretation of transient anomalous subdiffusion. I. Qualitative model
Biophysical Journal 92:1178–1191.

https://doi.org/10.1529/biophysj.106.092619
- PubMed
- Google Scholar
1. Śmigiel WM
2. Mantovanelli L
3. Linnik DS
4. Punter M
5. Silberberg J
6. Xiang L
7. Xu K
8. Poolman B
(2022) Protein diffusion in Escherichia coli cytoplasm scales with the mass of the complexes and is location dependent
Science Advances 8:eabo5387.

https://doi.org/10.1126/sciadv.abo5387
- PubMed
- Google Scholar
(1983) Observation of a noise-induced phase transition with an analog simulator
Physical Review Letters 51:1062–1065.

https://doi.org/10.1103/PhysRevLett.51.1062
- Google Scholar
1. Sokolov IM
(2010) Ito, Stratonovich, Hänggi and all the rest: The thermodynamics of interpretation
Chemical Physics 375:359–363.

https://doi.org/10.1016/j.chemphys.2010.07.024
- Google Scholar
1. Sprague BL
2. McNally JG
(2005) FRAP analysis of binding: proper and fitting
Trends in Cell Biology 15:84–91.

https://doi.org/10.1016/j.tcb.2004.12.001
- PubMed
- Google Scholar
(1997) Photobleaching recovery and anisotropy decay of green fluorescent protein GFP-S65T in solution and cells: cytoplasmic viscosity probed by green fluorescent protein translational and rotational diffusion
Biophysical Journal 72:1900–1907.

https://doi.org/10.1016/S0006-3495(97)78835-0
- PubMed
- Google Scholar
(2015) Molecular rotors provide insights into microscopic structural changes during protein aggregation
The Journal of Physical Chemistry. B 119:10170–10179.

https://doi.org/10.1021/acs.jpcb.5b05099
- PubMed
- Google Scholar
1. Tupper PF
2. Yang X
(2012) A paradox of state-dependent diffusion and how to resolve it
Proceedings of the Royal Society A 468:3864–3881.

https://doi.org/10.1098/rspa.2012.0259
- Google Scholar
1. Turelli M
(1977) Random environments and stochastic calculus
Theoretical Population Biology 12:140–178.

https://doi.org/10.1016/0040-5809(77)90040-5
- PubMed
- Google Scholar
(2017) Microorganisms maintain crowding homeostasis
Nature Reviews. Microbiology 15:309–318.

https://doi.org/10.1038/nrmicro.2017.17
- PubMed
- Google Scholar
1. Van Kampen NG
(1988) Diffusion in inhomogeneous media
Journal of Physics and Chemistry of Solids 49:673–677.

https://doi.org/10.1016/0022-3697(88)90199-0
- Google Scholar
1. Vekilov PG
(2010) Phase transitions of folded proteins
Soft Matter 6:5254.

https://doi.org/10.1039/c0sm00215a
- Google Scholar
1. Vishen AS
2. Prost J
3. Rao M
(2019) Breakdown of effective temperature, power law interactions, and self-propulsion in a momentum-conserving active fluid
Physical Review. E 100:062602.

https://doi.org/10.1103/PhysRevE.100.062602
- PubMed
- Google Scholar
1. Volpe G
2. Wehr J
(2016) Effective drifts in dynamical systems with multiplicative noise: a review of recent progress
Reports on Progress in Physics 79:053901.

https://doi.org/10.1088/0034-4885/79/5/053901
- PubMed
- Google Scholar
1. Wang Y
2. Li C
3. Pielak GJ
(2010) Effects of proteins on protein diffusion
Journal of the American Chemical Society 132:9392–9397.

https://doi.org/10.1021/ja102296k
- Google Scholar
1. Wang W
2. Cherstvy AG
3. Liu X
4. Metzler R
(2020) Anomalous diffusion and nonergodicity for heterogeneous diffusion processes with fractional Gaussian noise
Physical Review E 102:012146.

https://doi.org/10.1103/PhysRevE.102.012146
- Google Scholar
Preprint
1. Wu T
2. King MR
3. Farag M
4. Pappu RV
5. Lew MD
(2023) Single fluorogen imaging reveals spatial inhomogeneities within biomolecular condensates
bioRxiv.

https://doi.org/10.1101/2023.01.26.525727
- Google Scholar
1. Xiang L
2. Chen K
3. Yan R
4. Li W
5. Xu K
(2020) Single-molecule displacement mapping unveils nanoscale heterogeneities in intracellular diffusivity
Nature Methods 17:524–530.

https://doi.org/10.1038/s41592-020-0793-0
- PubMed
- Google Scholar
1. Xie J
2. Najafi J
3. Le Borgne R
4. Verbavatz JM
5. Durieu C
6. Sallé J
7. Minc N
(2022) Contribution of cytoplasm viscoelastic properties to mitotic spindle positioning
PNAS 119:2115593119.

https://doi.org/10.1073/pnas.2115593119
- Google Scholar
1. Yeomans JM
(2017) Nature’s engines: active matter
Europhysics News 48:21–25.

https://doi.org/10.1051/epn/2017204
- Google Scholar
(2008) Macromolecular crowding and confinement: biochemical, biophysical, and potential physiological consequences
Annual Review of Biophysics 37:375–397.

https://doi.org/10.1146/annurev.biophys.37.032807.125817
- PubMed
- Google Scholar

Article and author information

Author details

Achuthan Raja Venkatesh
1. Department of Biochemistry, Stanford University, Stanford, United States
2. Department of Biological Sciences, Indian Institute of Science Education and Research (IISER) Mohali, Mohali, India
Contribution
Conceptualization, Software, Formal analysis, Investigation, Visualization, Methodology, Writing - original draft, Writing - review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0009-0003-1868-4068
Kathy H Le

Department of Biochemistry, Stanford University, Stanford, United States

Contribution
Conceptualization, Investigation, Methodology, Writing - original draft, Writing - review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0009-0001-5575-7153
David M Weld

Department of Physics, University of California, Santa Barbara, Santa Barbara, United States

Contribution
Conceptualization, Resources, Formal analysis, Supervision, Investigation, Writing - original draft, Project administration, Writing - review and editing

For correspondence
weld@physics.ucsb.edu

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-4574-9491
Onn Brandman

Department of Biochemistry, Stanford University, Stanford, United States

Contribution
Conceptualization, Resources, Formal analysis, Supervision, Investigation, Writing - original draft, Project administration, Writing - review and editing

For correspondence
onn@stanford.edu

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-2084-154X

Funding

National Institutes of Health (R35GM153301)

Onn Brandman

National Science Foundation (DGE-1656518)

Kathy H Le

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We thank the Brandman Lab, Grant M Rotskoff, J E Ferrell Jr., Z Dogic, A Chaudhuri, and G Chu for helpful discussions. We thank P Guptasarma for helpful discussions and for facilitating the arrangement between IISER Mohali, UCSB, and Stanford. Simulations conducted in this study were run on the Sherlock high-performance computing cluster maintained by the Stanford Research Computing Center. ARV. is supported by the KVPY fellowship. K.H.L is supported by the NSF Graduate Research Fellowship Program. O B is funded by the National Institutes of Health grant R35GM153301.

Version history

Preprint posted: May 31, 2023
Sent for peer review: June 8, 2023
Reviewed Preprint version 1: September 28, 2023
Reviewed Preprint version 2: May 20, 2024
Version of Record published: June 18, 2024

Cite all versions

You can cite all versions using the DOI https://doi.org/10.7554/eLife.89794. This DOI represents all versions, and will always resolve to the latest one.

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

792

views
53

downloads
2

citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Mendeley

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Achuthan Raja Venkatesh
Kathy H Le
David M Weld
Onn Brandman

(2024)

Diffusive lensing as a mechanism of intracellular transport and compartmentalization

eLife 12:RP89794.

https://doi.org/10.7554/eLife.89794.3


Active noise/memory effects/other sources of nonequilibrium behaviors	Nonequilibrium in case of Itô and Stratonovich, equilibrium in case of isothermal	Consistent with equilibrium; thermal noise dominates
Diffusive lensing occurs	Diffusive lensing is seen in the Itô and Stratonovich cases (Figure 1—figure supplement 1B)	Diffusive lensing does not occur

Share this article

Cite this article

Low diffusivity leads to accumulation of particles.

Interaction-driven clustering is modulated by heterogenous diffusivity.

Heterogeneous diffusion alters bulk particle motion as measured by in silico microrheology.

A decrease in granule diffusivity, an increase in granule radius, or packing density slows down mesoscale dynamics.

Converting between stochastic integration conventions.

Author details

Achuthan Raja Venkatesh

Contribution

Competing interests

Kathy H Le

Contribution

Competing interests

David M Weld

Contribution

For correspondence

Competing interests

Onn Brandman

Contribution

For correspondence

Competing interests

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Categories and tags

Further reading