Edyta
Słyk
ab,
Tomasz
Skóra
a and
Svyatoslav
Kondrat
*acde
aInstitute of Physical Chemistry, Polish Academy of Sciences, 01-224 Warsaw, Poland. E-mail: svyatoslav.kondrat@gmail.com; skondrat@ichf.edu.pl
bDepartment of Theoretical Chemistry, Institute of Chemical Sciences, Faculty of Chemistry, Maria Curie-Skłodowska University in Lublin, 20-031 Lublin, Poland
cMax-Planck-Institut für Intelligente Systeme, Heisenbergstraße 3, 70569 Stuttgart, Germany
dIV. Institut für Theoretische Physik, Universität Stuttgart, Pfaffenwaldring 57, 70569 Stuttgart, Germany
eInstitut für Computerphysik, Universität Stuttgart, Allmandring 3, 70569 Stuttgart, Germany
First published on 30th June 2022
Diffusion in a macromolecularly crowded environment is essential for many intracellular processes, from metabolism and catalysis to gene transcription and translation. So far, theoretical and experimental work has focused on anomalous subdiffusion, and the effects of interactions, shapes, and composition, while the compactness or softness of macromolecules has received less attention. Herein, we use Brownian dynamics simulations to study how the softness of crowders affects macromolecular diffusion. We find that in most cases, soft crowders slow down the diffusion less effectively than hard crowders like Ficoll. For instance, at a 30% occupied volume fraction, the diffusion in Ficoll70 is about 20% slower than in soft crowders of the same size. However, our simulations indicate that elongated macromolecules, such as double-stranded DNA pieces, can diffuse comparably or even faster in hard crowders. We relate these effects to the volume excluded by soft and hard crowders to different tracers. Our results show that the softness and shape of macromolecules are crucial factors determining diffusion under crowding, relevant to diverse intracellular environments.
In simulations of macromolecular diffusion and crowding, one often uses coarse-grained or atomistic rigid-molecule models. However, macromolecules generally have flexible, often polymer-branching structures, which determine the interactions between the macromolecules at small separations. Recently, Blanco et al.25 introduced a shoulder-shaped, chain-entanglement softened potential (CESP, the second term in eqn (1); see also Fig. 1) to take into account the macromolecular softness and chain entanglements when modelling Dextran crowders. They found that this model reproduces the experimental data better than the frequently-used hard-sphere-like models. In more recent work, Junker et al.39 used fluorescent correlation spectroscopy (FCS) to investigate the diffusion of variously shaped macromolecules in compact crowders (Ficoll) and in more expanded, softer polyethylene glycol (PEG) and poly(ethylene oxide) (PEO) crowders. The experiments suggested a slower diffusion in PEG/PEO crowders (see, however, below and Fig. S6, ESI†), but we note that the soft and compact crowders had different sizes, which likely affected the macromolecular diffusion too.24,28
![]() | ||
Fig. 1 Model of soft interactions between macromolecules. (a) Schematics of two soft particles (top) and of a soft and a WCA (hard) particle (bottom). The dark spheres represent the hard core and the light grey spheres show the range of the soft shoulder potential given by eqn (1). The light grey spheres can overlap, while the dark cores cannot. (b) Relation between the parameters U0 and (ae − ac)/ac of the interaction potential, eqn (1), providing the same occupied volume vocc ≈ 555.6 nm3 as the hard (WCA) particle with the radius ac = 5.1 nm corresponding to Ficoll70. The occupied volume fractions were calculated based on the Boltzmann weight with the corresponding interaction potentials (Section S1A, ESI†). The symbols indicate the values used in panel (c) with the same color code. (c) Examples of the interaction potentials between two soft particles with the parameters from panel (b) and between two hard (WCA) particles corresponding to Ficoll70. All particles have the same occupied volume vocc ≈ 555.6 nm3. The thin vertical lines indicate the values of the hard core diameter 2ac. |
Herein, we study with Brownian dynamics (BD) simulations how crowders' softness (or compactness) affects the diffusion in crowded environments. We consider soft and hard particles with the same hydrodynamic radius, effectively occupying the same volume, and ask which particle is a better crowder and how the diffusion slowdown depends on the properties of a tracer.
In our simulations, we used the following total interaction potential between macromolecules
![]() | (1) |
To study how diffusion depends on macromolecular softness, we performed BD simulations of mixtures of soft and hard particles, both acting as tracers and crowders. We used a customized version of the package BD_BOX42,43 with the interaction potential given by eqn (1) and took into account hydrodynamic interactions via the generalized Rotne-Prager-Yamakawa (RPY) tensor.44–47 While the RPY approach is ubiquitously used in numerous studies, we note that it neglects many body interactions and near-field lubrication forces, which might be important for diffusion in crowded environments.21,48 We, therefore, restricted our considerations mainly to a relatively low volume fraction of 10% (except for Fig. 2e). We recall that in all simulations, the hydrodynamic radii of soft and hard particles were taken the same (aH = 5.1 nm), implying the same short-time diffusion coefficients. From the BD trajectories, we calculated time-averaged mean-square displacements TAMSD(t) = 〈|r(t) − r(0)|2〉 of the geometrical center r of the macromolecules and extracted the corresponding long-time diffusion coefficients D as detailed in Section S2D of the ESI.† We note that we used TAMSD instead of ensemble averaged (EA) MSD to get better statistics; in the long-time limit, both TAMSD and EAMSD should give the same results (for recent discussions regarding the relations between time and ensemble averaging, see ref. 49 and 50).
![]() | ||
Fig. 2 Effect of macromolecule's softness on diffusion. (a) Snapshots from BD simulations of mixture of soft and hard particles. (b) Long-time diffusion coefficients of hard and soft particles as functions of the softness parameter (eqn (1)). The value ae = ac corresponds to a hard (WCA) particle. The soft particle occupies the same volume as the hard particle (Fig. 1). Occupied volume fraction ϕocc = 10% and molar fraction of hard particles xhard ≈ 27%. (c) Long-time diffusion coefficients of hard and soft particles as functions of molar fraction of hard particles for occupied volume fraction 10%. (d) Snapshots from BD simulations of two systems containing only soft particles and only hard particles. (e) Long-time diffusion coefficients as functions of the occupied volume fraction for systems consisting of hard particles and soft particles. For radial distribution functions, see Fig. S2–S4 (ESI†). See Table S4 (ESI†) for the number of particles in each system. |
We first studied how diffusion depends on the range of entanglement, i.e., ae–ac (ae = ac corresponds to a compact, hard particle). We simulated a mixture of hard and soft macromolecules of the total occupied volume fraction 10% and molar fraction of hard particles xhard ≈ 27%. Fig. 2b shows that in this setup, softer crowders provide a weaker slowdown of the macromolecular diffusion. At low excluded volume fractions, the diffusion coefficient can be roughly estimated by51,52
D(ϕocc) = D0(1 − αϕex(ϕocc)), | (2) |
Our simulations show a relatively weak dependence of the diffusion coefficient on the molar fraction of hard particles in the mixture (Fig. 2c); this is in line with virtually no change in the radial distribution function observed in the same systems (Fig. S3, ESI†). Nevertheless, increasing the fraction of hard particles slightly slowed down the diffusion in all cases. Consistent with Fig. 2b, at all molar fractions, the diffusion of soft particles was faster than the diffusion of hard particles, which is again due to the lower excluded volumes (Table S7, ESI†).
The effect of the total volume fraction in unmixed systems exhibiting self-crowding was more substantial (Fig. 2d). For ϕocc = 30%, the self-diffusion in the system with hard particles was about 20% slower than in the system composed of soft particles. The differences between the two systems are also reflected in the radial distribution functions, which show a monotonic to oscillatory behaviour with increasing ϕocc for hard particles, but only small quantitative changes in the case of soft ones (Fig. S4, ESI†). Somewhat surprisingly, the difference in the excluded volumes in both systems was only a few per cent (Table S8, ESI†), incompatible with the large difference in the diffusion coefficients. This is reflected in the parameter α in eqn (1). Indeed, fitting this equation to the simulation data gives significantly different values α ≈ 0.24 for hard and α ≈ 0.06 for soft particles. In the calculations of the excluded volumes, we measured the interaction energies in the Boltzmann factor with respect to the infinite separation between the particles. However, in such crowded systems, most macromolecules' repulsive shoulders overlap and the particles experience a lower energy difference when approaching each other from the average distance (determined by the occupied volume fraction) compared to infinite dilution. This shift in the energy level likely reduces the excluded volume as seen by a tracer, providing a faster diffusion.
In addition to the diffusion of hard and soft particles due to self-crowding, we also studied how differently shaped macromolecules diffuse in hard and soft crowders. We chose a cylindrically-shaped, 16 nm-long double-stranded (ds) DNA and Y-shaped Immunoglobulin G (IgG) (Fig. 3a). We modelled a dsDNA as a quasi-rigid charged polymer consisting of eight connected spherical beads, each with the hydrodynamic radius of 1.14 nm (see ref. 26). As in our previous work,26 the dsDNA beads interacted electrostatically with each other via the DLVO potential with the screening length 0.95 nm (all other particles were uncharged). For IgG, we took the model that we introduce in the upcoming article53 (Section S1B, ESI†), developed to reproduce the experimentally measured angle distribution and hydrodynamic radius of IgG.
![]() | ||
Fig. 3 Effect of macromolecular shapes and softness on diffusion. (a) Models of dsDNA (≈16 nm long) and Immunoglobulin G (IgG) and snapshots from BD simulations of dsDNA and IgG in hard and soft crowders. Hydrodynamic radii aH have been obtained from BD simulations of single macromolecules (i.e., in infinite dilution). (b) Long-time diffusion coefficients and (c) excluded volume fractions of dsDNA, IgG and Ficoll70 (hard particle) in soft and hard crowders. See Table S5 in the ESI,† for the number of particles in each system. (d) Difference in the excluded volume fractions Δϕex = ϕhardex − ϕsoftex for a dsDNA modelled as a quasi-rigid polymer of length ![]() |
Our simulations show that for both crowders, the diffusion coefficients are reduced more significantly for larger tracers (Fig. 3b), in accord with previous studies.21,30,39 Moreover, both IgG and Ficoll70 diffused considerably faster in soft crowders. This is similar to Fig. 2. In stark contrast, a smaller and more elongated dsDNA was faster (albeit only weakly) in hard crowders. To understand this result, we calculated with MC simulations the volume fractions excluded to these tracers by soft and hard crowders, which showed similar trends (Fig. 3c). We further investigated this effect by calculating the difference in the excluded volumes for hard and soft crowders, Δϕex = ϕhardex − ϕsoftex, as a function of the crowder size and tracer's aspect ratio. Fig. 3d shows the results in the form of a heatmap and demonstrates well-defined regions of positive and negative Δϕex. We observed a similar behaviour in a wide range of occupied volume fractions (Fig. S5, ESI†). The results show that the excluded volume for an elongated dsDNA is larger in soft crowders, provided the crowder is sufficiently large. In contrast to a hard crowder, a soft crowder can interact repulsively with a few dsDNA beads at the same time, and more than one soft crowder can repel a dsDNA bead if the system is sufficiently crowded. Clearly, these effects increase the excluded volume when the length of a polymer increases or when the size of the polymer bead decreases. Fig. 3d suggests larger differences in the diffusion coefficients in the case of bigger crowders or smaller bead sizes of a polymeric tracer (Fig. 3e). Indeed, Fig. 3f shows that soft crowders slow down the diffusion of a polymer with smaller beads more substantially than hard crowders.
It is instructive to compare our results with FCS measurements by Junker et al.39 who studied the diffusion of variously shaped tracers in Ficoll and polyethylene glycol/poly(ethylene oxide) (PEG/PEO), which are softer, more expanded crowders compared to Ficoll. These authors based their calculations of the occupied volumes on specific molecular volumes and found that the diffusion was faster in Ficoll than in PEG. However, recalculating the occupied volume based on the hydrodynamic radii, as in our simulations, gives a roughly comparable diffusion slowdown (Fig. S6, ESI†), unlike our simulation results. Although the precise reason for this discrepancy is unclear, we note that Ficoll and PEG had different hydrodynamic radii, contributing differently to the diffusion's slowdown.24,28 In addition, streptavidin (hydrodynamic radius aH ≈ 3.2 nm) is much smaller than Ficoll70 (aH ≈ 5.1 nm) used as a tracer in our simulations. Finally, Junker et al.39 considered relatively high occupied volume fractions (above 25% after recalculations), at which polyethylene glycol could form polymeric networks, leading to sieving effects, not captured by the CESP potential of ref. 25 that we utilized in our work. In future work, it will be useful to modify the model to take such effects into account on a computationally inexpensive, coarse-grained level and extend experiments to cover a wider range of occupied volume fractions and macromolecular varieties.
In conclusion, we have studied how macromolecular diffusion is influenced by the softness of crowders. We found that in most cases, soft crowders reduced the tracer diffusion to a much lesser extent than hard crowders of the same size (Fig. 2). However, our BD simulations indicated that cylindrically-shaped macromolecules could diffuse similarly or even faster in hard crowders (Fig. 3). We related these effects to the volume fractions excluded by soft and hard crowders to different tracers. These results are relevant to intracellular environments composed of macromolecules of different shapes and softness, and we hope that our work will motivate further more detailed simulation and experimental studies on the relation between macromolecular structure and diffusion in highly crowded systems.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d2sm00357k |
This journal is © The Royal Society of Chemistry 2022 |