New insights into diffusion in 3D crowded media by Monte Carlo simulations: effect of size, mobility and spatial distribution of obstacles

Particle diffusion in crowded media was studied through Monte Carlo simulations in 3D obstructed lattices. Three particular aspects affecting the diffusion, not extensively treated in a three-dimensional geometry, were analysed: the relative particle-obstacle size, the relative particle-obstacle mobility and the way of having the obstacles distributed in the simulation space (randomly or uniformly). The results are interpreted in terms of the parameters that characterize the time dependence of the diffusion coefficient: the anomalous diffusion exponent (α), the crossover time from anomalous to normal diffusion regimes (τ) and the long time diffusion coefficient (D*). Simulation results indicate that there are a more anomalous diffusion (smaller α) and a lower long time diffusion coefficient (D*) when obstacle concentration increases, and that, for a given total excluded volume and immobile obstacles, the anomalous diffusion effect is less important for bigger size obstacles. However, for the case of mobile obstacles, this size effect is inverted yielding values that are in qualitatively good agreement with in vitro experiments of protein diffusion in crowded media. These results underline that the pattern of the spatial partitioning of the obstacle excluded volume is a factor to be considered together with the value of the excluded volume itself.     

Anomalous and Normal Diffusion of Tracers in Crowded Environments: E ect of Size Disparity between Tracer and Crowders

The dynamics of tracers in crowded matrix is of interest in various areas of physics, such as the diffusion of proteins in living cells. By using two-dimensional (2D) Langevin dynamics simulations, we investigate the diffusive properties of a tracer of a diameter in crowded environments caused by randomly distributed crowders of a diameter. Results show that the emergence of subdiffusion of a tracer at intermediate time scales depends on the size ratio of the tracer to crowders δ. If δ falls between a lower critical size ratio and a upper one, the anomalous diffusion occurs purely due to the molecular crowding. Further analysis indicates that the physical origin of subdiffusion is the "cage effect". Moreover, the subdiffusion exponent α decreases with the increasing medium viscosity and the degree of crowding, and gets a minimum αmin=0.75 at δ=1. At long time scales, normal diffusion of a tracer is recovered. For δ≤1, the relative mobility of tracers is independent of the degree of crowding. Meanwhile, it is sensitive to the degree of crowding for δ>1. Our results are helpful in deepening the understanding of the diffusive properties of biomacromolecules that lie within crowded intracellular environments, such as proteins, DNA and ribosomes.     

Diffusion of spheres in crowded suspensions of rods

Translational tracer diffusion of spherical macromolecules in crowded suspensions of rodlike colloids is investigated. Experiments are done using several kinds of spherical tracers in fd-virus suspensions. A wide range of size ratios L/2a of the length L of the rods and the diameter 2a of the tracer sphere is covered by combining several experimental methods: fluorescence correlation spectroscopy for small tracer spheres, dynamic light scattering for intermediate sized spheres, and video microscopy for large spheres. Fluorescence correlation spectroscopy is shown to measure long-time diffusion only for relatively small tracer spheres. Scaling of diffusion coefficients with a/xi, predicted for static networks, is not found for our dynamical network of rods (with xi the mesh size of the network). Self-diffusion of tracer spheres in the dynamical network of freely suspended rods is thus fundamentally different as compared to cross-linked networks. A theory is developed for the rod-concentration dependence of the translational diffusion coefficient at low rod concentrations for freely suspended rods. The proposed theory is based on a variational solution of the appropriate Smoluchowski equation without hydrodynamic interactions. The theory can, in principle, be further developed to describe diffusion through dynamical networks at higher rod concentrations with the inclusion of hydrodynamic interactions. Quantitative agreement with the experiments is found for large tracer spheres, and qualitative agreement for smaller spheres. This is probably due to the increasing importance of hydrodynamic interactions as compared to direct interactions as the size of the tracer sphere decreases.     

Correlating anomalous diffusion with lipid bilayer membrane structure using single molecule tracking and atomic force microscopy

Anomalous diffusion has been observed abundantly in the plasma membrane of biological cells, but the underlying mechanisms are still unclear. In general, it has not been possible to directly image the obstacles to diffusion in membranes, which are thought to be skeleton bound proteins, protein aggregates, and lipid domains, so the dynamics of diffusing particles is used to deduce the obstacle characteristics. We present a supported lipid bilayer system in which we characterized the anomalous diffusion of lipid molecules using single molecule tracking, while at the same time imaging the obstacles to diffusion with atomic force microscopy. To explain our experimental results, we performed lattice Monte Carlo simulations of tracer diffusion in the presence of the experimentally determined obstacle configurations. We correlate the observed anomalous diffusion with obstacle area fraction, fractal dimension, and correlation length. To accurately measure an anomalous diffusion exponent, we derived an expression to account for the time-averaging inherent to all single molecule tracking experiments. We show that the length of the single molecule trajectories is critical to the determination of the anomalous diffusion exponent. We further discuss our results in the context of confinement models and the generating stochastic process.     

Calculation of diffusion coefficients for carbon dioxide + solute system near the critical conditions by non-equilibrium molecular dynamics simulation

A non-equilibrium molecular dynamics simulation was adopted to calculate the diffusion coefficients for a pseudo-binary system of carbon dioxide and for a carbon dioxide + solute system at 308.2 and 318.2 K. The calculated results were compared with the self- and tracer diffusion coefficients calculated by an equilibrium molecular dynamics simulation. The simulated results for the pseudo-binary system of carbon dioxide by the non-equilibrium molecular dynamics simulation are in good agreement with the results of self diffusion coefficients for pure carbon dioxide by the equilibrium molecular dynamics simulation. The simulated results of mutual diffusion coefficients for the carbon dioxide + solute system by the non-equilibrium molecular dynamics simulation are slightly lower than the results of the tracer diffusion coefficients by the equilibrium molecular dynamics simulation. The anomalous behavior of diffusion coefficients near the critical concentration was represented by the results of the non-equilibrium molecular dynamics simulation.     

Anomalous nanoparticle diffusion in polymer solutions and melts: a mode-coupling theory study

Mode-coupling theory is employed to study diffusion of nanoparticles in polymer melts and solutions. Theoretical results are directly compared with molecular dynamics simulation data for a similar model. The theory correctly reproduces the effects of the nanoparticle size, mass, particle-polymer interaction strength, and polymer chain length on the nanoparticle diffusion coefficient. In accord with earlier experimental, simulation, and theoretical work, it is found that when the polymer radius of gyration exceeds the nanoparticle radius, the Stokes-Einstein relation underestimates the particle diffusion coefficient by as much as an order of magnitude. Within the mode-coupling theory framework, a microscopic interpretation of this phenomenon is given, whereby the total diffusion coefficient is decomposed into microscopic and hydrodynamic contributions, with the former dominant in the small particle limit, and the latter dominant in the large particle limit. This interpretation is in agreement with previous mode-coupling theory studies of anomalous diffusion of solutes in simple dense fluids.     

Theory of asymmetric electrochemical stochastic diffusion

Strict analysis of electrochemical strochastic diffusion due to asymmetric Brownian motion of electric charge in an equilibrium electrochemical ac circuit containing double electric layer capacitance and noisy Faradaic resistance is carried out. Cumulant analogs (for 3rd and 5th order correlations) of the Einstein formula are obtained. It is proved that equilibrium asymmetric (nongauss) stochastic diffusion is in agreement with the central limiting theorem of the probability theory. The Hurst exponent was found in the case of the nongauss components of the process of equilibrium stochastic diffusion. Apart from electrochemistry, the performed stochastic analysis of equilibrium electrochemical nongauss diffusion is also of general theoretical interest, including its application in the stochastic theory of asymmetric anomalous transport and strict theory of fluctuations at large deviations from equilibrium.     

Movement of proteins in an environment crowded by surfactant micelles: anomalous versus normal diffusion

Small proteins move in crowded cell compartments by anomalous diffusion. In many of them, e.g., the endoplasmic reticulum, the proteins move between lipid membranes in the aqueous lumen. Molecular crowding in vitro offers a systematic way to study anomalous and normal diffusion in a well controlled environment not accessible in vivo. We prepared a crowded environment in vitro consisting of hexaethylene glycol monododecyl ether (C(12)E(6)) nonionic surfactant and water and observed lysozyme diffusion between elongated micelles. We have fitted the data obtained in fluorescence correlation spectroscopy using an anomalous diffusion model and a two-component normal diffusion model. For a small concentration of surfactant (below 4 wt %) the data can be fitted by single-component normal diffusion. For larger concentrations the normal diffusion fit gave two components: one very slow and one fast. The amplitude of the slow component grows with C(12)E(6) concentration. The ratio of diffusion coefficients (slow to fast) is on the order of 0.1 for all concentrations of surfactant in the solution. The fast diffusion is due to free proteins while the slow one is due to the protein-micelle complexes. The protein-micelle interaction is weak since even in a highly concentrated solution (35% of C(12)E(6)) the amplitude of the slow mode is only 10%, despite the fact that the average distance between the micelles is the same as the size of the protein. The anomalous diffusion model gave the anomaly index (r(2)(t) approximately t(alpha)), alpha monotonically decreasing from alpha = 1 (at 4% surfactant) to alpha = 0.88 (at 37% surfactant). The fits for two-component normal diffusion and anomalous diffusion were of equally good quality, but the physical interpretation was only straightforward for the former.     

Toward observation of single-file diffusion using the tracer zero-length column method

The tracer zero-length column (ZLC) method has been employed to study the diffusion of toluene in one-dimensional ZSM-12 and SAPO-5 zeolites. A significant deviation in the shape of the measured tracer exchange curves from monoexponential behavior was observed for toluene diffusion in both adsorbents in the limit of long-time asymptotes. In contrast, water/ZSM-12 and acetylene/SAPO-5 systems exhibit tracer exchange curves that are close to monoexponential behavior. Monoexponential curves are usually observed for systems obeying normal (Fickian) diffusion. Such diffusion is expected for the latter two systems because the diameters of both sorbates are less than the radii of their corresponding host channels. The differences in the shape of the tracer exchange curves for large and small sorbates can be explained by assuming the occurrence of anomalous, single-file diffusion for large sorbates in narrow, one-dimensional channels.     

Translational Dynamics of Lipidated Ras Proteins in the Presence of Crowding Agents and Compatible Osmolytes

Ras proteins are small GTPases and are involved in transmitting signals that control cell growth, differentiation, and proliferation. Since the cell cytoplasm is crowded with different macromolecules, understanding the translational dynamics of Ras proteins in crowded environments is crucial to yielding deeper insight into their reactivity and function. Herein, the translational dynamics of lipidated N‐Ras and K‐Ras4B is studied in the bulk and in the presence of a macromolecular crowder (Ficoll) and the compatible osmolyte and microcrowder sucrose by fluorescence correlation spectroscopy. The results reveal that N‐Ras forms dimers due to the presence of its lipid moiety in the hypervariable region, whereas K‐Ras4B remains in its monomeric form in the bulk. Addition of a macromolecular crowding agent gradually favors clustering of the Ras proteins. In 20 wt % Ficoll N‐Ras forms trimers and K‐Ras4B dimers. Concentrations of sucrose up to 10 wt % foster formation of N‐Ras trimers and K‐Ras dimers as well. The results can be rationalized in terms of the excluded‐volume effect, which enhances the association of the proteins, and, for the higher concentrations, by limited‐hydration conditions. The results of this study shed new light on the association state of these proteins in a crowded environment. This is of particular interest for the Ras proteins, because their solution state—monomeric or clustered—influences their membrane‐partitioning behavior and their interplay with cytosolic interaction partners.     

Spatio-temporal anomalous diffusion in heterogeneous media by nuclear magnetic resonance

In this paper, we describe nuclear magnetic resonance measurements of water diffusion in highly confined and heterogeneous colloidal systems using an anomalous diffusion model. For the first time, temporal and spatial fractional exponents, α and μ, introduced within the framework of continuous time random walk, are simultaneously measured by pulsed gradient spin-echo NMR technique in samples of micro-beads dispersed in aqueous solution. In order to mimic media with low and high level of disorder, mono-dispersed and poly-dispersed samples are used. We find that the exponent α depends on the disorder degree of the system. Conversely, the exponent μ depends on both bead sizes and magnetic susceptibility differences within samples. The new procedure proposed here may be a useful tool to probe porous materials and microstructural features of biological tissue.     

14CO2 and CO2 Transport in polycarbonate: Measurement of the time lag and permeability

Transport of CO₂ across polycarbonate films has been studied using a diffusion cell technique employing a radioactively labelled tracer (¹⁴CO₂). Because the ¹⁴CO₂ driving force could be established independently of the unlabelled CO₂ driving force, several classes of experiments not possible with conventional techniques were performed. These different classes of experiments showed measurably different time lags. Formally, these experiments all are limiting cases of the more general mixed-gas permeation problem; however, simplifying assumptions in the dual sorption theory are possible because the tracer concentration approaches zero and because the two species in this special mixed-gas problem exhibit the same dual sorption parameters. These simplifications allow derivation of analytical expressions for the time lag for both the unlabelled and labelled gas species. The experimental measurements are in good agreement with the dual sorption model formulated for the mixed-gas tracer diffusion problem.     

Tracer diffusion in fluctuating block copolymer melts

We present a theoretical investigation of the tracer diffusion of diblock copolymers and homopolymers in a thermally fluctuating block copolymer melt above the order-disorder transition (ODT) temperature. Entanglement effects and differences in monomeric friction coefficients are ignored; hence, the theory should be most applicable to short copolymers with rheologically similar blocks. Overall, we find that the diffusion rates of both tracer block copolymers and homopolymers in a block copolymer melt are suppressed when compared with diffusivities in a strictly homogeneous medium with the same average composition. This mobility suppression is due to thermally excited composition fluctuations in block copolymer melts near the ODT; the latter result in transient potential barriers to diffusion. We explore the dependence of the tracer diffusion coefficient on molecular weights and compositions of both matrix and tracer, as well as temperature. A comparison of our theoretical predictions to recent experiments by T. Lodge and coworkers shows qualitative agreement. © 1996 John Wiley & Sons, Inc.     

Atomistic Simulation of Lysozyme in Solutions Crowded by Tetraethylene Glycol: Force Field Dependence

The behavior of biomolecules in crowded environments remains largely unknown due to the accuracy of simulation models and the limited experimental data for comparison. Here we chose a small crowder of tetraethylene glycol (PEG-4) to investigate the self-crowding of PEG-4 solutions and molecular crowding effects on the structure and diffusion of lysozyme at varied concentrations from dilute water to pure PEG-4 liquid. Two Amber-like force fields of Amber14SB and a99SB-disp were examined with TIP3P (fast diffusivity and low viscosity) and a99SB-disp (slow diffusivity and high viscosity) water models, respectively. Compared to the Amber14SB protein simulations, the a99SB-disp model yields more coordinated water and less PEG-4 molecules, less intramolecular hydrogen bonds (HBs), more protein–water HBs, and less protein–PEG HBs as well as stronger interactions and more hydrophilic and less hydrophobic contacts with solvent molecules. The a99SB-disp model offers comparable protein–solvent interactions in concentrated PEG-4 solutions to that in pure water. The PEG-4 crowding leads to a slow-down in the diffusivity of water, PEG-4, and protein, and the decline in the diffusion from atomistic simulations is close to or faster than the hard sphere model that neglects attractive interactions. Despite these differences, the overall structure of lysozyme appears to be maintained well at different PEG-4 concentrations for both force fields, except a slightly large deviation at 370 K at low concentrations with the a99SB-disp model. This is mainly attributed to the strong intramolecular interactions of the protein in the Amber14SB force field and to the large viscosity of the a99SB-disp water model. The results indicate that the protein force fields and the viscosity of crowder solutions affect the simulation of biomolecules under crowding conditions.     

Solute size effects on the solvation structure and diffusion of ions in liquid methanol under normal and cold conditions

We have performed a series of molecular dynamics simulations of alkali metal (Li+, Na+, K+, Rb+, and Cs+) and halide (F-, Cl-, Br-, and I-) ions in liquid methanol at two different temperatures to investigate the effects of ion size on the hydration structure and diffusion of ions in methanol under normal and cold conditions. Simulations are also carried out for some of the larger cations such as I+, (CH3)4N+, and (C2H5)4N+ and also neutral alkali metal atoms in methanol at both temperatures. With the increase of ion size, the diffusion coefficients of both positive and negative ions are found to show anomalous behavior. For cations, it is found that the maximum of the diffusion coefficient versus ion size curve occurs at the rather large cation of (CH3)4N+ unlike in water where the maximum occurs at the relatively smaller ion of Rb+. For halide ions, the anomalous behavior, i.e., the increase of diffusion with ion size, continues up to iodide ion and no maximum is observed. These results are in good agreement with experimental observations. The diffusion coefficients of neutral atoms are found to be greater in methanol than that in water and they decrease monotonically with solute size, whereas the diffusion coefficients of the corresponding ions are found to be smaller in methanol. Accordingly, an ion experiences a smaller Stokes friction and a higher dielectric friction in methanol than in water. These contrasting effects are believed to be responsible for the shift of the maximum of ion diffusion toward a larger ion size when compared with similar anomalous size dependence in liquid water.     

Accurate tracer diffusion coefficients of Na+ and Cl− ions in dilute aqueous sodium chloride solutions measured with the closed capillary method

A simple and absolute closed capillary method is introduced for measuring tracer diffusion coefficients in liquids with both - and -active tracers. In this method a narrow capillary is partially filled with a labelled solution and an unlabelled solution is used to fill the rest of the capillary. The diffusion coefficient is obtained from the time dependence of the monitored activity when the length of the capillary is known. The method has been tested by remeasuring the tracer diffusion coefficients of²²NaCl and Na³⁶Cl for solutions covering a wide range of total NaCl concentration. The average precision of the measurements was about 0.3%. The results obtained for Na³⁶Cl tracer diffusion are in good agreement with the data found in literature. The²²NaCl tracer diffusion coefficients that were measured in dilute solutions agree well with those obtained using the continuous open-ended capillary method but differ from results for solutions between 0.1 and 1M obtained with the diaphragm cell method.     

Tracer diffusion in colloidal gels

Computer simulations were done of the mean square displacement (MSD) of tracer particles in colloidal gels formed by diffusion or reaction limited aggregation of hard spheres. The diffusion coefficient was found to be determined by the volume fraction accessible to the spherical tracers (phi a) independent of the gel structure or the tracer size. In all cases, critical slowing down was observed at phi a approximately 0.03 and was characterized by the same scaling laws reported earlier for tracer diffusion in a Lorentz gas. Strong heterogeneity of the MSD was observed at small phi a and was related to the size distribution of pores.     

Heat Effects in ZLC Experiments

The problem of nonisothermal desorption in a zero length column (ZLC) experiment is considered theoretically. Simple analytical expressions for the ZLC desorption curve are derived for certain limiting situations in which the governing equations reduce to a linear form. More general numerical solutions are calculated for a wide range of experimental conditions assuming both negligible mass transfer resistance and finite mass transfer resistance controlled by intraparticle diffusion. A simple criterion for negligible thermal effects is developed. It is shown that when the ZLC technique is applied to the measurement of diffusion in unaggregated zeolite crystals, as originally intended, heat effects are generally insignificant. However, when applied to the measurement of macropore diffusion in relatively large adsorbent particles heat effects can become important and may cause major modification of both the desorption rate and the shape of the desorption curve. A recent experimental ZLC study carried out with commercial adsorbent particles, under conditions of macropore diffusion control, showed an anomalous dependence of the desorption rate on both temperature and particle size. These effects can be qualitatively explained by the nonisothermal model. A more precise quantitative representation of these experiments will require a more refined model incorporating a nonlinear equilibrium isotherm as well as intraparticle diffusional resistance.     

Effect of Molecular Crowding on the Temperature–Pressure Stability Diagram of Ribonuclease A

FT‐IR spectroscopic and thermodynamic measurements were designed to explore the effect of a macromolecular crowder, dextran, on the temperature and pressure‐dependent phase diagram of the protein Ribonuclease A (RNase A), and we compare the experimental data with approximate theoretical predictions based on configuration entropy. Exploring the crowding effect on the pressure‐induced unfolding of proteins provides insight in protein stability and folding under cell‐like dense conditions, since pressure is a fundamental thermodynamic variable linked to molecular volume. Moreover, these studies are of relevance for understanding protein stability in deep‐sea organisms, which have to cope with pressures in the kbar range. We found that not only temperature‐induced equilibrium unfolding of RNase A, but also unfolding induced by pressure is markedly prohibited in the crowded dextran solutions, suggesting that crowded environments such as those found intracellularly, will also oppress high‐pressure protein unfolding. The FT‐IR spectroscopic measurements revealed a marked increase in unfolding pressure of 2 kbar in the presence of 30 wt % dextran. Whereas the structural changes upon thermal unfolding of the protein are not significantly influenced in the presence of the crowding agent, through stabilization by dextran the pressure‐unfolded state of the protein retains more ordered secondary structure elements, which seems to be a manifestation of the entropic destabilization of the unfolded state by crowding.