## Propagation of Enzyme-driven Active Fluctuations in Crowded Milieu

Rik Chakraborty, 1 Arnab Maiti, 1 Diptangshu Paul, 2 Rajnandan Borthakur, 2 K. R. Jayaprakash, 2 Uddipta Ghosh, 2 and Krishna Kanti Dey 1, ∗

2 Department of Mechanical Engineering, Indian Institute of Technology Gandhinagar, Palaj, Gujarat 382055, India

1 Laboratory of Soft and Living Materials, Department of Physics, Indian Institute of Technology Gandhinagar, Palaj, Gujarat 382055, India (Dated: August 2, 2024)

We investigated the energy transfer from active enzymes to their surroundings in crowded environments by measuring the diffusion of passive microscopic tracers in active solutions of ficoll and glycerol. Despite observing lower rates of substrate turnover and relatively smaller enhancement of passive tracer diffusion in artificial crowded media compared to those in aqueous solutions, we found a significantly higher relative diffusion enhancement in crowded environments in the presence of enzymatic activity. Our experimental observations, coupled with supporting analytical estimations, underscored the critical role of the intervening media in facilitating mechanical energy distribution around active enzymes.

## I. INTRODUCTION

The ability of active enzymes to generate mechanical fluctuations during catalytic turnovers has been recently a focus area of investigation in active matter physics. [16] Several experimental [7-10] and theoretical [11-14] studies, have demonstrated that active enzymes, while catalysing substrates, could enhance their diffusion under aqueous environment, the exact mechanism of which remains still to be fully understood. The observations however, opened a new area of investigation in molecular biophysics, where mechanical energy transduction by single molecule was previously believed to be unlikely. [15] While the enhancement of catalytic enzymes' diffusion was subject to both positive and negative arguments [1621], it was theoretically proposed [22-27] and experimentally demonstrated [28-31] that active fluctuations generated by enzymes during catalysis could also influence the dynamics of nearby passive tracers through long-ranged interactions. The observed capability of active enzyme molecules to transduce mechanical energy from catalytic turnover and thereby influence the dynamics of their surroundings, akin to the ATP-driven molecular motors, presents an intriguing avenue for further exploration. The observations were important as molecular motors have been reported to drive diffusive-like, non-thermal motion of cellular components, influencing the overall metabolic state of the cell. [32, 33] The dynamic coupling observed between catalytic enzymes and their surroundings provides compelling evidence that these molecules, beyond their primary role as catalysts in biochemical reactions, have secondary functions and may serve as mediators of dynamic interaction pathways. [34] However, considering that the cytoskeletal transport by motor proteins occurs in highly dissipative and crowded environments, [35-37] it is imperative to investigate whether the mechanical forces generated by the enzymes are perceivable in dense and crowded intracellular environments. A positive answer to this question would also warrant investigation to check the influence of the crowded environment on the generation and propagation of reactioninduced fluctuations, and assess the potential implications of such active processes on intracellular particle dynamics. Successful exploration of such phenomena is expected to uncover new research avenues, significantly advancing our understanding of activity-induced transport, assembly, and other related phenomena. With this broader perspective in mind, in this study, we studied the mechanical energy transfer from active enzymes within a dense and crowded environment by estimating the diffusion of microscopic passive tracers suspended in such media. To mimic intracellular viscosity and crowding, we used molecular crowders such as ficoll 400 and glycerol, as commonly used in other studies. [38, 39] The diffusion of passive tracers was considerably lower in the crowder solutions compared to that in water, primarily due to the higher viscosity of the former. Notably, enzymatic substrate turnover rates in crowded environments were also slower than those in aqueous media, which could further contribute to reducing passive tracer diffusion in these media compared to that in water. Interestingly, upon initiating enzymatic reactions, we observed enhanced diffusion of passive tracers, even in crowded environments, similar to observations in aqueous solutions. [28, 30] However, the relative increase in diffusion in crowded environments, observed with both ficoll and glycerol, was significantly greater than that in aqueous media. We conducted experiments measuring tracer diffusion in three different media (water, ficoll, and glycerol) under different reaction rate conditions. Based on the results obtained in these experiments we hypothesized that artificial crowded environments facilitate a greater transfer of reaction-generated fluctuations around active enzymes compared to dilute aqueous environments.

∗ Author to whom correspondence should be addressed:k.dey@iitgn.ac.in

In order to validate our hypothesis, we developed a model considering that the binding and unbinding of

substrates in enzymes largely governed their conformational fluctuations resulting in the generation of stochastic forces. The actual conversion of reactant to products, on the other hand, was likely to contribute to the net thermal energy of the system, which was previously found to be insignificant in enhancing single enzyme diffusion. [40] The conformation change-induced forces propagated through the surrounding medium in the form of high-frequency pressure pulses. We assumed the active fluctuations propagating around enzymes' vicinity as spherical pressure waves and computed the corresponding average energy density at specific distances away from the point of actuation. The analytical estimates found nice qualitative match with the experimental results, suggesting that the intervening media played a crucial role in deciding the energy propagation around catalytically active enzymes, which might have important implications over various biochemical processes in active intracellular milieu.

## II. RESULTS AND DISCUSSIONS

3 µ msized polystyrene latex beads (LB30, 10% solids), catalase from bovine liver (catalogue no: C9322), urease from Canavalia ensiformis (Jack bean) (catalogue no: U1500), urea (catalogue no: U5128), phenol red (catalogue no: P3532), glycerol (catalogue no: G5516), ficoll PM 400 (catalogue no: F4375) were purchased from Sigma-Aldrich. 30% (w/v) hydrogen peroxide (H O ) 2 2 solution (catalogue no: 107209) and 95-97% (w/w) sulphuric acid (H SO ) 2 4 solution (catalogue no: 100731) were purchased from Merck. To prepare the experimental tracer solution, 10 µ L of the stock solution was diluted in a series of steps, with each step having a dilution ratio of 1 : 10, using deionized (DI) water. This resulted in a secondary stock solution with a solid concentration of 0.1%. The secondary stock solution was then further diluted using DI water (pH 7 and resistivity of 18.2 MΩ · cm) or 10% (v/v) glycerol or 3% (w/v) ficoll 400 to prepare the final experimental solutions containing 0.005% polymer particles. For experiments, all measurements were conducted at room temperature of 25 ◦ C. Commercially procured glass slides (Blue Star, PIC-1 microslides, 76 mm L × 26 mm W × 1 35 mm H), were washed . thoroughly using Labolene solution, followed by sonicating them for 10 min. After sonication, the glass slides were treated with piranha solution, which is a mixture of H SO 2 4 and H O 2 2 in the ratio of 7 : 3, for a period of 1 h. Next, the slides were removed carefully from the piranha solution and washed with deionized (DI) water. Finally, they were dried using an air drier. A small circular plastic chamber was then securely attached to the glass slide using suitable adhesive followed by air drying at room temperature.

We used optical video microscopy and single particle tracking technique to characterize the motion of polymer particles over the entire field of view in a quasi-2D con- dition. Experimental videos were recorded using a Nikon Eclipse Ti2 Inverted Microscope, with an objective of 60 × , 10 frames per second. Each video was recorded for 5 min after allowing a 30 min stabilization period to ensure that it had reached a steady state. For active samples, after the stabilization period, an additional 5 min were given to minimize the effects of perturbations arising due to substrate addition. Frames from each recorded video were analyzed and tracked through particle tracking codes developed by Crocker and Grier. [41] Mean squared displacements (MSD) and diffusion coefficients were estimated using a code developed in MATLAB. The results obtained at different experimental conditions were plotted and fitted using OriginPro.

The dynamic viscosities of DI water, 10% (v/v) glycerol and 3% (w/v) ficoll 400 at 25 ◦ C were measured using Anton Paar MCR 702 rheometer. The viscosities were measured at the shear rates ranging from 10 s -1 to 100 s -1 , within which all three experimental solutions were found to be Newtonian. The viscosity values of DI water, 10% (v/v) glycerol and 3% (w/v) ficoll 400, averaged over the shear rates were found to be 0.95 ± 0.01, 1.32 ± 0.02, and 1.48 ± 0.01 mPa-s respectively.

To measure the activity of catalase in different media, the decomposition of its substrate, H 2 O , was monitored 2 using UV-Vis spectroscopy by recording its absorbance at 240 nm. Urease assay was performed by measuring absorbance of the dye phenol red at 560 nm. The enzyme and substrate concentrations were selected to ensure sufficient reaction times during tracer diffusion measurement. Student's t-test was performed to confirm the statistical significance of the experimental results. The measurements with a p-value less than 0.05 were considered to be statistically significant.

As mentioned earlier, the mechanical fluctuations generated by active enzymes has previously been reported to be strong enough to influence the diffusion of nearby passive tracers under aqueous environments. A key question is whether such effects are also observable under cytosolic crowded conditions, where reactions catalysed by several enzymes lead to important physiological functions and life-sustaining processes. To check this, we performed experiments with molecules of ficoll 400 and glycerol to mimic cytocellular environments where the diffusion of 3 µ m polystyrene (PS) tracers was studied in the presence of substrate catalysis by different enzymes. Concentrations of the crowders were chosen in such a way that their volume fractions were almost the same (10%) with nearly identical viscosities (1.3-1.5 mPa-s). This is similar to the intracellular viscosity, although the cellular interior is highly inhomogeneous and compartmentalized that are often characterized by spatio-temporal viscosity variations. [42, 43] Importantly, under crowded environment, the rate of substrate turnover by enzymes gets lowered owing to one or more of the following reasons: a) constrained diffusion of substrates to the enzymes' active sites, b) slower rate of conformation change of the enzymes during catalysis and c) accumulation of product

molecules near the enzymes post catalysis [44-46]. In this scenario, it becomes crucial to explore whether the active fluctuations generated by enzymes can overcome dissipative effects and impact tracer diffusion in their vicinity. Our experimental results suggest that fluctuations generated by 1 nM catalase and 10 mM H O 2 2 were sufficient to influence the motion of the polymer microparticles in glycerol and ficoll. Similar observations were recorded with 10 nM urease and 100 mM urea which showed that this effect was likely to be generic and could very well be detectable for other enzyme-substrate pairs.

Initially, the tracer diffusion was measured in DI water where the diffusion coefficient was found to be 7.79 × 10 -14 m 2 / s, which was nearly 54% less than the value predicted by the Stokes-Einstein equation (Fig. 1). This discrepancy may be attributed to the boundary effects on the particles which remained close to the bottom surface of the experimental chamber, at a height, h ∼ 5 µ m. Faucheux and Libchaber characterized this hindered diffusion in terms of a dimensionless length scale parameter γ given by, γ = h -r r [47]. For our experiments with r = 1.5 µ m particles, the value of γ came out to be 2.33, which yielded a diffusion coefficient ratio D D 0 ∼ 0.55. This resulted in the estimated hindered diffusion coefficient to be D = 7.99 × 10 -14 m 2 / s, which matched nicely with our experimental results. The diffusion coefficients of the polymer particles in glycerol and ficoll were measured to be 5.91 × 10 -14 m 2 / s and 5.42 × 10 -14 m 2 / s respectively, which also matched well with the estimated hindered diffusion coefficient values of 6.06 × 10 -14 m 2 / s and 5.40 × 10 -14 m 2 / s respectively (Fig. 1).

FIG. 1. (a) The mean-squared displacement (MSD) profiles, and (b) diffusion coefficients of 3 µ mtracer particles in (A) DI water, (B) 10% (v/v) glycerol, and (C) 3% (w/v) ficoll 400. The symbol ∗∗∗ denotes the significance level of p &lt; 0 001. .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00578v1-figure-f204f498-0dc0-4ecc-a6da-12540d718f29.jpg)

We introduced non-thermal fluctuations in our system using two different biochemical reactions - the disproportionation of H O 2 2 with catalase and hydrolysis of urea with urease. Previous experiments conducted in aqueous environments demonstrated enhanced diffusion of passive tracers in the presence of free enzyme activity, where the increase in diffusion was found to correlate with enzymatic activity. [28] Importantly, consistent with predictions of theoretical studies, [22, 23, 25] our experimental results presented herein confirmed that fluctuations generated by active enzymes are perceptible even in crowded environments and are capable of influencing the dynam- ics of their surroundings. Interestingly, we found that the relative enhancement in tracer diffusion is notably higher in crowded environments compared to that measured under aqueous conditions. In the presence of active catalse, the tracer diffusion was enhanced by 18.6% and 18.3% in 10% glycerol and 3% ficoll respectively (Fig. 2a). Similarly with active urease, diffusion enhancements of 20.6% and 10.7% were recorded in glycerol and ficoll respectively (Fig. 2b). To confirm that the free enzymes in solution did not get adsorbed over the polymer bead surface during experiments and influenced their propulsion, experiments were conducted with microparticles coated with a thin layer of bovine serum albumin (BSA). These coated microparticles also exhibited a similar enhancement in diffusion in the presence of catalysis. Furthermore, both BSA coated and uncoated tracers, when dispersed in fluorescently tagged enzyme solution showed no signs of enzyme attachment on their surfaces.

FIG. 2. Diffusion coefficients of tracer particles in the absence and presence of (a) catalase activity (1 nM catalase and 10 mM H O ), and (b) urease activity (10 nM urease 2 2 and 100 mM urea) in (A) DI water, (B) 10% (v/v) glycerol, and (C) 3% (w/v) ficoll 400 respectively. The symbols ∗∗ and ∗∗∗ denote the significance levels of p &lt; 0 01 and . p &lt; 0 001 . respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00578v1-figure-d83a80e4-fa6a-471a-8736-7093ae45025e.jpg)

Considering the correlation observed between enzyme reaction rate and tracer diffusion in aqueous solutions, one possible explanation for greater enhancement in tracer diffusion within crowded environments could be the higher enzyme activity in crowded media, potentially boosting tracer diffusion therein. To assess the catalytic activity of catalase and urease in glycerol and ficoll, we conducted enzyme activity assays to determine the catalytic reaction rates r. We found that the substrate catalysis rates for catalase decreased by 19.9% and 34.7% in 10% glycerol and 3% ficoll respectively, compared to their rates in DI water (Fig. 3a). Similarly, for urease, the reaction rate decreased by 40.8% and 96.9% in glycerol and ficoll respectively (Fig. 3b). We further investigated the relative tracer diffusion, denoted as, ∆D rel = D D -0 D 0 × 100% corresponding to different catalase reaction rates r, in different media (Fig. 3c). Assuming linearity between ∆D rel and r, we calculated the rate of change of ∆D rel with the reaction rate r from the slopes of (Fig. 3c). These slopes were subsequently normalized to obtain the relative enhancement in tracer diffusion in different crowded media compared to water, denoted as Φ norm , for identical reaction rate conditions. The varia-

tion of Φ norm for different media is depicted in (Fig. 3d), which clearly indicates that the catalytic reaction rate was not primarily responsible for the higher energy transfer from enzymes to the tracers under artificial crowded environments.

FIG. 3. Normalized reaction rates for (a) active catalase (1 nM catalase and 10 mM H O ), and (b) active urease 2 2 (10 nM urease and 100 mM urea) in (A) DI water, (B) 10% (v/v) glycerol, and (C) 3% (w/v) ficoll 400. All the reaction rates were measured during 300-600 s from the beginning of the reactions. (c) Relative diffusion enhancement ∆D rel of the tracer particles for different reaction rates r for active catalase in (A) DI water, (B) 10% (v/v) glycerol, and (C) 3% (w/v) ficoll 400. Experimental data were linearly fitted using OriginPro software. (d) Normalized relative tracer diffusion enhancement per unit reaction rate (Φ norm ), estimated for different media. The symbols ∗∗ , and ∗∗∗ denote the significance levels of p &lt; 0 05 and . p &lt; 0 001 respectively. . NS signifies not significant.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00578v1-figure-638a9ea9-a612-40ca-8825-9817836a9597.jpg)

The reaction generated heat can also be shown to be inefficient in boosting tracer diffusion in crowded solutions. In line with the estimations reported earlier, [40] simple calculations showed that upon each catalytic cycle, the increase in temperature across a single catalase molecule could amount to ∆T ≈ 10 -5 K. Considering that the entire amount of heat generated upon substrate turnover was responsible for influencing nearby tracer diffusion, the corresponding thermophoretic tracer velocity could be estimated as V ST ≃ D S 0 T ∆T R , where D , 0 S T and R are the normal diffusion coefficient, Soret coefficient and radius of the tracer particle. Considering S T ∼ 150 K -1 , [48] the change in the effective tracer diffusion coefficient due to thermophoresis is given by ∆D ≃ V 2 ST D r , where D r is the rotational diffusion coefficient of the tracer. This yields, ∆D ≈ 10 -20 m 2 / s giving ∆D D 0 ≈ 10 -4 %. This is too small to account for the experimentally observed enhancement in tracer diffusion upon catalysis. Moreover, considering that the crowder solutions used were sufficiently dilute, their thermal conductivities and resultant increase in tracer diffusion are expected to be similar to those of DI water.

Since neither the enzyme catalysis nor the reaction generated heat was found to be primarily responsible for greater enhancement in tracer diffusion in crowded environment, we next investigated the role of intervening media in facilitating energy transfer from the enzyme's active sites to the tracers. We developed a simple analytical model, where the binding and unbinding of substrates, resulting in conformational fluctuations in enzymes, were assumed to generate stochastic forces in the media. These forces were considered to be propagating around active enzymes in the form of high-frequency pressure pulses. For different actuation frequencies, we calculated the corresponding average energy densities at specific distances away from the point of actuation. The analysis finds its motivation from the fact that the intracellular environement is inherently viscoelastic, inhomegeneous and crowded, which is expected to play a key role in the propagation of enzyme-generated active fluctuations in its vicinity.

## III. ANALYTICAL STUDIES

We developed a qualitative understanding of the enhanced energy transfer in 10% glycerol and 3% ficoll, compared to that in water based on fundamental principles of fluid dynamics. [49] We considered the disturbances produced by catalytic reactions occurring at an enzyme site propagating through the media and reaching the tracer particles as spherical pressure waves. The energy density of these waves around an enzyme has been estimated. This is calculated by considering a simple case where an enzyme molecule was assumed to interact with single nearby tracer through reaction generated force fields, as shown in Fig. 4a. Starting from the equation of continuity, and using the Navier-Stokes equation for an isentropic compressible fluid with irrotational flows and negligible temperature gradient, we used the spherical polar coordinates to obtain the energy density as function of both space and time.

Although the timescale of conformational fluctuations corresponding to different domains of a single enzyme has been reported to vary over an extended range, [50] for our estimation, we assumed that the stress generated in the surrounding media following substrate binding and product release was related to one primary conformational change (which may or may not correspond to the catalytic step). We assumed the angular frequency of this primary conformational change to be ω in water, which generated an acoustic wave with a maximum pressure amplitude of p 0 at the surface of the enzyme. We also considered the pressure amplitude of this wave to be consistent in crowded conditions, although it propagated with a lower angular frequency in these media, primarily determined by the viscosity of the medium. This assumption found its rationale in previous experimental studies where the rate of conformation changes in enzymes was

found to decrease with the increase in viscosity of the media η . [51] For our case, we took the angular frequency corresponding to the primary conformational fluctuation in crowded media to decrease as 1 η . To estimate the wave number k corresponding to the acoustic wave propagation in different media, we also considered their density ρ 0 , velocity of sound c 0 therein, and viscosity µ 0 values. The values considered and the corresponding angular frequencies (Ω) for wave propagation in these media are given in Table. I.

TABLE I. Properties of different experimental solutions and their assumed effect on primary conformational fluctuation frequency of enzymes. The density and velocity of sound values for 10% glycerol and 3% ficoll are taken from [52, 53] and [54], respectively. The viscosity values of these solutions were measured experimentally.

| Medium       |   ρ 0 [kg/m 3 |   c 0 [m/s] |   µ 0 [mPa-s] | Ω [s - 1 ]   |
|--------------|---------------|-------------|---------------|--------------|
| Water        |          1000 |        1481 |          1    | ω            |
| 10% glycerol |          1026 |        1532 |          1.32 | ω /1.32      |
| 3% ficoll    |          1007 |        1504 |          1.48 | ω /1.48      |

FIG. 4. (a) Schematic diagram of the enzyme catalase and a tracer particle. Spatial variation of the (b) pressure, and (c) the energy density at fixed instance of time in different medium. (d) Normalized root mean square value of the average energy density ( ⟨ e a ( ) t ⟩ norm ) in (A) water, (B) 10% glycerol, and (C) 3% ficoll. Catalase structure images is taken from rcsb.org (PDB IDs: 1TGU).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00578v1-figure-e9b985b0-feff-453d-9d51-8359b117d69d.jpg)

The estimated wave number k led to the calculation of enzyme fluctuation generated pressure p ′ ( r, t ) and the velocity u ′ ( r, t ) fields, along with the expression for the spatiotemporal variation of energy density e ′ ( r, t ).

$$\sum _ { \substack { p ^ { \prime } ( r, t ) = \text{Re} \left [ \frac { p _ { 0 } r _ { c } } { r } e ^ { i \{ k ( r - r _ { c } ) - \Omega t \} } \right ] } } ( 1 )$$

$$u ^ { \prime } ( r, t ) = \text{Re} \left [ \frac { p _ { 0 } \Omega r _ { c } } { \rho _ { 0 } c _ { 0 } ^ { 2 } } \left ( \frac { 1 } { k r } + \frac { \mu \dot { i } } { k ^ { 2 } r ^ { 2 } } \right ) e ^ { i \{ k ( r - r _ { c } ) - \Omega t \} } \dot { \Big | } \ \ ( 2 ) \quad \text{$O_{u}$}$$

$$e ^ { \prime } ( r, t ) = \frac { \rho _ { 0 } } { 2 } u ^ { \prime } ( r, t ) ^ { 2 } + \frac { 1 } { 2 \rho _ { 0 } c _ { 0 } ^ { 2 } } p ^ { \prime } ( r, t ) ^ { 2 } \ ( 3 )$$

Here, r c is the radius of the enzyme. We also considered r t to be the distance between the center of the enzyme and the tracer surface (Fig. 4a). The average energy density over a volume bounded by [ r , r c t ] is then given by

$$e _ { a } ( t ) = \frac { 4 \pi } { V _ { 0 } } \int _ { r _ { c } } ^ { r _ { t } } e ^ { \prime } ( r, t ) r ^ { 2 } d r$$

where, V 0 = 4 π r ( 3 t -r 3 c ) / 3 is the volume of the spherical shell with its center coinciding with the center of the enzyme (Fig. 4a).

The integral in Eq. (4) was evaluated numerically using the Gauss-Legendre quadrature method. For quantitative estimates, we considered ω = 10 8 Hz and r t = 1000 µ m. Considering the force generated per catalytic turnover ∼ 10 pN, as reported earlier [8, 55] and enzyme radius r c ∼ 5 nm, the instantaneous pressure generated per catalytic turnover came out to be ∼ 0.1 MPa. We took the value of maximum pressure amplitude p 0 to be of this order, which resulted in the spatial variation of the pressure and energy density at a fixed instance of time, as shown in Fig. 4b and 4c, respectively. We then obtained the root mean square value of the average energy density, ⟨ e a ( ) t ⟩ , over a time period corresponding to the actuation frequency. The ratio of ⟨ e a ( ) t ⟩ for glycerol and ficoll with respect to water, represented as ⟨ e a ( ) t ⟩ norm , has been shown in Fig. 4d. The estimated trend of ⟨ e a ( ) t ⟩ norm variation qualitatively matched well with the experimentally measured energy transfer, characterized by Φ norm , as shown in Fig. 3d. We performed the analysis for two other values of Ω (10 7 and 10 9 Hz), which also resulted in similar trend of ⟨ e a ( ) t ⟩ norm in different experimental media. We finally noted that the estimated energy transfer values (measured by ⟨ e a ( ) t ⟩ norm ) were relatively lesser in magnitude compared to the experimentally measured values of Φ norm . The difference in the magnitudes may be attributed to the assumptions made in the model - such as the consideration of (1) interaction between a single enzyme and tracer and disregard for potential influences from multiple particles, (2) only a single frequency of conformational fluctuations responsible for energy transfer from the enzyme to the tracer. It may be noted that the viscosity of the medium played a key role only in reducing the actuation frequency of the enzyme as mentioned earlier. Owing to the small length scales of the system, the attenuation of the pressure wave had negligible dependence on the viscosity of the medium. The sharp decrement in the spatial components of p ′ ( r, t ), u ′ ( r, t ), and e ′ ( r, t ) with distance were solely due to the spherical nature of the reaction generated pressure waves.

## IV. CONCLUSIONS

Our experimental findings suggested that the change in tracer diffusion (∆D) per unit enzyme reaction rate

for 10% (v/v) glycerol and 3% (w/v) ficoll 400 are 4.38 and 4.71 times higher than that measured in aqueous environment (Fig. 3d). Analytical estimates based on the theory of wave propagation in fluids also suggested that the energy transfer from active enzymes to their immediate neighbourhood enhanced in glycerol and ficoll compared to that in aqueous media (Fig. 4d). Thus our observation led to the insights into the propagation of enzymatic forces in active crowded environments. The augmented transport of molecules and particles through enzyme generated non-thermal mechanical fluctuations carries significant scientific and technological implications. The generation of forces induced by catalysis may provide an adequate explanation for the observed stochastic motion of the cytoplasm, cytoplasmic glass transitions, and improved mixing observed during metabolic transformations in bacterial cells. [33, 56] The movement of proteins

- [1] H. Kitahata and A. S. Mikhailov, J. Phys. Soc. Jpn. 92 , 054801 (2023).
- [2] S. Ghosh, A. Somasundar, and A. Sen, Annu. Rev. Condens. Matter Phys. 12 , 177 (2021).
- [3] M. Feng and M. K. Gilson, Annu. Rev. Biophys. 49 , 87 (2020).
- [4] Y. Zhang and H. Hess, ACS Cent. Sci. 5 , 939 (2019).
- [5] A.-Y. Jee, Y.-K. Cho, S. Granick, and T. Tlusty, Proc. Natl. Acad. Sci. U.S.A. 115 , E10812 (2018).
- [6] R. D. Astumian, ACS Nano 8 , 11917 (2014).
- [7] H. S. Muddana, S. Sengupta, T. E. Mallouk, A. Sen, and P. J. Butler, J. Am. Chem. Soc. 132 , 2110 (2010).
- [8] C. Riedel, R. Gabizon, C. A. Wilson, K. Hamadani, K. Tsekouras, S. Marqusee, S. Press´ e, and C. Bustamante, Nature 517 , 227 (2015).
- [9] A.-Y. Jee, T. Tlusty, and S. Granick, Proc. Natl. Acad. Sci. U.S.A. 117 , 29435 (2020).
- [10] M. Xu, J. L. Ross, L. Valdez, and A. Sen, Phys. Rev. Lett. 123 , 128101 (2019).
- [11] S. Kondrat and M. N. Popescu, Phys. Chem. Chem. Phys. 21 , 18811 (2019).
- [12] A. Ryabov and M. Tasinkevych, Soft Matter 18 , 3234 (2022).
- [13] D. R. Slochower and M. K. Gilson, Biophys. J. 114 , 2174 (2018).
- [14] P. Illien, T. Adeleke-Larodo, and R. Golestanian, Europhys. Lett. 119 , 40002 (2017).
- [15] E. M. Purcell, Am. J. Phys. 45 , 3 (1977).
- [16] X. Bai and P. G. Wolynes, J. Chem. Phys. 143 , 165101 (2015).
- [17] J.-P. G¨nther, M. B¨rsch, u o and P. Fischer, Acc. Chem. Res. 51 , 1911 (2018).
- [18] Y. Zhang, M. J. Armstrong, N. M. B. Kazeruni, and H. Hess, Nano Lett. 18 , 8025 (2018).
- [19] A.-Y. Jee, K. Chen, T. Tlusty, J. Zhao, and S. Granick, J. Am. Chem. Soc. 141 , 20062 (2019).
- [20] Z. Chen, A. Shaw, H. Wilson, M. Woringer, X. Darzacq, S. Marqusee, Q. Wang, and C. Bustamante, Proc. Natl. Acad. Sci. U.S.A. 117 , 21328 (2020).

within the cytoplasmic matrix is vital for their effective function across various organelles, and enzyme-generated forces could potentially serve as key orchestrators of this intricate dynamic interplay.

## ACKNOWLEDGMENTS

KKD thanks Science and Engineering Research Board (SERB), India (ECR/2017/002649, CRG/2023/007588), Department of Science and Technology (DST), India (DST/ICD/BRICS/PilotCall3/BioTheraBubble/2019), Ministry of Education, Government of India (MoESTARS/STARS-2/2023-0620) and IIT Gandhinagar for financial support. We also thank Sanchari Swarupa, Aditi Saha, Ayushi Bhatt, Sirsha Ganguly, and Nividha for insightful discussions. RC thanks IIT Gandhinagar for a research fellowship.

- [21] A. A. Choi, H. H. Park, K. Chen, R. Yan, W. Li, and K. Xu, J. Am. Chem. Soc. 144 , 4839 (2022).
- [22] A. S. Mikhailov and R. Kapral, Proc. Natl. Acad. Sci. U.S.A. 112 , E3639 (2015).
- [23] R. Kapral and A. S. Mikhailov, Physica D 318 , 100 (2016).
- [24] M. Dennison, R. Kapral, and H. Stark, Soft Matter 13 , 3741 (2017).
- [25] Y. Koyano, H. Kitahata, and A. S. Mikhailov, Europhys. Lett. 128 , 40003 (2020).
- [26] A. K. Tripathi, T. Das, G. Paneru, H. K. Pak, and T. Tlusty, Commun. Phys. 5 , 101 (2022).
- [27] Y. Hosaka, S. Komura, and A. S. Mikhailov, Soft Matter 16 , 10734 (2020).
- [28] X. Zhao, K. K. Dey, S. Jeganathan, P. J. Butler, U. M. C´ ordova-Figueroa, and A. Sen, Nano Lett. 17 , 4807 (2017).
- [29] K. K. Dey, Angew. Chem. Int. Ed. 58 , 2208 (2019).
- [30] X. Lin and Y. He, Anal. Chem. 94 , 7158 (2022).
- [31] A. Maiti, Y. Koyano, H. Kitahata, and K. K. Dey, Phys. Rev. E 109 , 054607 (2024).
- [32] C. P. Brangwynne, G. H. Koenderink, F. C. MacKintosh, and D. A. Weitz, J. Cell Biol. 183 , 583 (2008).
- [33] M. Guo, A. J. Ehrlicher, M. H. Jensen, M. Renz, J. R. Moore, R. D. Goldman, J. Lippincott-Schwartz, F. C. Mackintosh, and D. A. Weitz, Cell 158 , 822 (2014).
- [34] N. Tokuriki and D. S. Tawfik, Science 324 , 203 (2009).
- [35] I. Goychuk, V. O. Kharchenko, and R. Metzler, PLoS ONE 9 , e91700 (2014).
- [36] A. Raguin, N. Kern, and A. Parmeggiani, J. Theor. Biol. 505 , 110370 (2020).
- [37] C. Leduc, K. Padberg-Gehle, V. Varga, D. Helbing, S. Diez, and J. Howard, Proc. Natl. Acad. Sci. U.S.A. 109 , 6100 (2012).
- [38] X. Ge, D. Luo, and J. Xu, PLoS ONE 6 , e28707 (2011).
- [39] N. O. Junker, F. Vaghefikia, A. Albarghash, H. H¨fig, o D. Kempe, J. Walter, J. Otten, M. Pohl, A. Katranidis, S. Wiegand, and J. Fitter, J. Phys. Chem. B 123 , 4477 (2019).
- [40] R. Golestanian, Phys. Rev. Lett. 115 , 108102 (2015).

- [41] J. C. Crocker and D. G. Grier, J. Colloid Interface Sci. 179 , 298 (1996).
- [42] A. Srivastava and G. Krishnamoorthy, Arch. Biochem. Biophys. 340 , 159 (1997).
- [43] O. O. Adeniba, E. A. Corbin, A. Ganguli, Y. Kim, and R. Bashir, Sci. Rep. 10 , 12803 (2020).
- [44] H. Rastogi and P. K. Chowdhury, Biochim. Biophys. Acta Proteins Proteom. 1869 , 140699 (2021).
- [45] S. S. Paul, P. Sil, R. Chakraborty, S. Haldar, and K. Chattopadhyay, Biochemistry 55 , 2332 (2016).
- [46] S. Shahid, M. I. Hassan, A. Islam, and F. Ahmad, Biochim. Biophys. Acta Gen. Subj. 1861 , 178 (2017).
- [47] L. P. Faucheux and A. J. Libchaber, Phys. Rev. E 49 , 5158 (1994).
- [48] L. Helden, R. Eichhorn, and C. Bechinger, Soft Matter 11 , 2379 (2015).
- [49] S. Temkin, Elements of acoustics (John Wiley &amp; Sons, 1981).
- [50] M. A. Maria-Solano, E. Serrano-Herv´s, a A. RomeroRivera, J. Iglesias-Fern´ndez, a and S. Osuna, Chem. Commun. 54 , 6622 (2018).
- [51] A. Sekhar, M. P. Latham, P. Vallurupalli, and L. E. Kay, J. Phys. Chem. B 118 , 4546 (2014).
- [52] A. Volk and C. J. K¨hler, Exp. Fluids a 59 , 75 (2018).
- [53] Nividha, A. K. Shukla, A. Maiti, B. Bharat, K. Jayaprakash, S. Khatua, P. Thareja, and K. K. Dey, ChemistrySelect 8 , e202301570 (2023).
- [54] W. Qiu, J. T. Karlsen, H. Bruus, and P. Augustsson, Phys. Rev. Appl. 11 , 024018 (2019).
- [55] P. J. Butler, K. K. Dey, and A. Sen, Cell. Mol. Bioeng. 8 , 106 (2015).
- [56] B. R. Parry, I. V. Surovtsev, M. T. Cabeen, C. S. O'Hern, E. R. Dufresne, and C. Jacobs-Wagner, Cell 156 , 183 (2014).
Rik Chakraborty Arnab Maiti Diptangshu Paul Rajnandan Borthakur K. R. Jayaprakash Uddipta Ghosh Krishna Kanti Dey
Propagation of Enzyme-driven Active Fluctuations in Crowded Milieu
We investigated the energy transfer from active enzymes to their surroundings
in crowded environments by measuring the diffusion of passive microscopic
tracers in active solutions of ficoll and glycerol. Despite observing lower
rates of substrate turnover and relatively smaller enhancement of passive
tracer diffusion in artificial crowded media compared to those in aqueous
solutions, we found a significantly higher relative diffusion enhancement in
crowded environments in the presence of enzymatic activity. Our experimental
observations, coupled with supporting analytical estimations, underscored the
critical role of the intervening media in facilitating mechanical energy
distribution around active enzymes.