Anomalous versus Slowed-Down Brownian Diffusion in the Ligand-Binding Equilibrium

Measurements of protein motion in living cells and membranes consistently report transient anomalous diffusion (subdiffusion) that converges back to a Brownian motion with reduced diffusion coefficient at long times after the anomalous diffusion regime. Therefore, slowed-down Brownian motion could be considered the macroscopic limit of transient anomalous diffusion. On the other hand, membranes are also heterogeneous media in which Brownian motion may be locally slowed down due to variations in lipid composition. Here, we investigate whether both situations lead to a similar behavior for the reversible ligand-binding reaction in two dimensions. We compare the (long-time) equilibrium properties obtained with transient anomalous diffusion due to obstacle hindrance or power-law-distributed residence times (continuous-time random walks) to those obtained with space-dependent slowed-down Brownian motion. Using theoretical arguments and Monte Carlo simulations, we show that these three scenarios have distinctive effects on the apparent affinity of the reaction. Whereas continuous-time random walks decrease the apparent affinity of the reaction, locally slowed-down Brownian motion and local hindrance by obstacles both improve it. However, only in the case of slowed-down Brownian motion is the affinity maximal when the slowdown is restricted to a subregion of the available space. Hence, even at long times (equilibrium), these processes are different and exhibit irreconcilable behaviors when the area fraction of reduced mobility changes.     

Controlling Anomalous Diffusion in Lipid Membranes

Diffusion in cell membranes is not just simple two-dimensional Brownian motion but typically depends on the timescale of the observation. The physical origins of this anomalous subdiffusion are unresolved, and model systems capable of quantitative and reproducible control of membrane diffusion have been recognized as a key experimental bottleneck. Here, we control anomalous diffusion using supported lipid bilayers containing lipids derivatized with polyethylene glycol (PEG) headgroups. Bilayers with specific excluded area fractions are formed by control of PEG lipid mole fraction. These bilayers exhibit a switch in diffusive behavior, becoming anomalous as bilayer continuity is disrupted. Using a combination of single-molecule fluorescence and interferometric imaging, we measure the anomalous behavior in this model over four orders of magnitude in time. Diffusion in these bilayers is well described by a power-law dependence of the mean-square displacement with observation time. Anomaleity in this system can be tailored by simply controlling the mole fraction of PEG lipid, producing bilayers with diffusion parameters similar to those observed for anomalous diffusion in biological membranes.     

Anomalous Extracellular Diffusion in Rat Cerebellum

Extracellular space (ECS) is a major channel transporting biologically active molecules and drugs in the brain. Diffusion-mediated transport of these substances is hindered by the ECS structure but the microscopic basis of this hindrance is not fully understood. One hypothesis proposes that the hindrance originates in large part from the presence of dead-space (DS) microdomains that can transiently retain diffusing molecules. Because previous theoretical and modeling work reported an initial period of anomalous diffusion in similar environments, we expected that brain regions densely populated by DS microdomains would exhibit anomalous extracellular diffusion. Specifically, we targeted granular layers (GL) of rat and turtle cerebella that are populated with large and geometrically complex glomeruli. The integrative optical imaging (IOI) method was employed to evaluate diffusion of fluorophore-labeled dextran (MW 3000) in GL, and the IOI data analysis was adapted to quantify the anomalous diffusion exponent d_w from the IOI records. Diffusion was significantly anomalous in rat GL, where d_w reached 4.8. In the geometrically simpler turtle GL, d_w was elevated but not robustly anomalous (d_w = 2.6). The experimental work was complemented by numerical Monte Carlo simulations of anomalous ECS diffusion in several three-dimensional tissue models containing glomeruli-like structures. It demonstrated that both the duration of transiently anomalous diffusion and the anomalous exponent depend on the size of model glomeruli and the degree of their wrapping. In conclusion, we have found anomalous extracellular diffusion in the GL of rat cerebellum. This finding lends support to the DS microdomain hypothesis. Transiently anomalous diffusion also has a profound effect on the spatiotemporal distribution of molecules released into the ECS, especially at diffusion distances on the order of a few cell diameters, speeding up short-range diffusion-mediated signals in less permeable structures.     

Anomalous Diffusion Induced by Cristae Geometry in the Inner Mitochondrial Membrane

Diffusion of inner membrane proteins is a prerequisite for correct functionality of mitochondria. The complicated structure of tubular, vesicular or flat cristae and their small connections to the inner boundary membrane impose constraints on the mobility of proteins making their diffusion a very complicated process. Therefore we investigate the molecular transport along the main mitochondrial axis using highly accurate computational methods. Diffusion is modeled on a curvilinear surface reproducing the shape of mitochondrial inner membrane (IM). Monte Carlo simulations are carried out for topologies resembling both tubular and lamellar cristae, for a range of physiologically viable crista sizes and densities. Geometrical confinement induces up to several-fold reduction in apparent mobility. IM surface curvature per se generates transient anomalous diffusion (TAD), while finite and stable values of projected diffusion coefficients are recovered in a quasi-normal regime for short- and long-time limits. In both these cases, a simple area-scaling law is found sufficient to explain limiting diffusion coefficients for permeable cristae junctions, while asymmetric reduction of the junction permeability leads to strong but predictable variations in molecular motion rate. A geometry-based model is given as an illustration for the time-dependence of diffusivity when IM has tubular topology. Implications for experimental observations of diffusion along mitochondria using methods of optical microscopy are drawn out: a non-homogenous power law is proposed as a suitable approach to TAD. The data demonstrate that if not taken into account appropriately, geometrical effects lead to significant misinterpretation of molecular mobility measurements in cellular curvilinear membranes.     

Anomalous Surface Diffusion of Protons on Lipid Membranes

The cellular energy machinery depends on the presence and properties of protons at or in the vicinity of lipid membranes. To asses the energetics and mobility of a proton near a membrane, we simulated an excess proton near a solvated DMPC bilayer at 323 K, using a recently developed method to include the Grotthuss proton shuttling mechanism in classical molecular dynamics simulations. We obtained a proton surface affinity of −13.0 ± 0.5 kJ mol⁻¹. The proton interacted strongly with both lipid headgroup and linker carbonyl oxygens. Furthermore, the surface diffusion of the proton was anomalous, with a subdiffusive regime over the first few nanoseconds, followed by a superdiffusive regime. The time- and distance dependence of the proton surface diffusion coefficient within these regimes may also resolve discrepancies between previously reported diffusion coefficients. Our simulations show that the proton anomalous surface diffusion originates from restricted diffusion in two different surface-bound states, interrupted by the occasional bulk-mediated long-range surface diffusion. Although only a DMPC membrane was considered in this work, we speculate that the restrictive character of the on-surface diffusion is highly sensitive to the specific membrane conditions, which can alter the relative contributions of the surface and bulk pathways to the overall diffusion process. Finally, we discuss the implications of our findings for the energy machinery.     

Anomalous Surface Diffusion of Protons on Lipid Membranes

The cellular energy machinery depends on the presence and properties of protons at or in the vicinity of lipid membranes. To asses the energetics and mobility of a proton near a membrane, we simulated an excess proton near a solvated DMPC bilayer at 323 K, using a recently developed method to include the Grotthuss proton shuttling mechanism in classical molecular dynamics simulations. We obtained a proton surface affinity of −13.0 ± 0.5 kJ mol⁻¹. The proton interacted strongly with both lipid headgroup and linker carbonyl oxygens. Furthermore, the surface diffusion of the proton was anomalous, with a subdiffusive regime over the first few nanoseconds, followed by a superdiffusive regime. The time- and distance dependence of the proton surface diffusion coefficient within these regimes may also resolve discrepancies between previously reported diffusion coefficients. Our simulations show that the proton anomalous surface diffusion originates from restricted diffusion in two different surface-bound states, interrupted by the occasional bulk-mediated long-range surface diffusion. Although only a DMPC membrane was considered in this work, we speculate that the restrictive character of the on-surface diffusion is highly sensitive to the specific membrane conditions, which can alter the relative contributions of the surface and bulk pathways to the overall diffusion process. Finally, we discuss the implications of our findings for the energy machinery.     

Spontaneous Detachment of Colloids from Primary Energy Minima by Brownian Diffusion

The Derjaguin-Landau-Verwey-Overbeek (DLVO) interaction energy profile has been frequently used to interpret the mechanisms controlling colloid attachment/detachment and aggregation/disaggregation behavior. This study highlighted a type of energy profile that is characterized by a shallow primary energy well (i.e., comparable to the average kinetic energy of a colloid) at a small separation distance and a monotonic decrease of interaction energy with separation distance beyond the primary energy well. This energy profile is present due to variations of height, curvature, and density of discrete physical heterogeneities on collector surfaces. The energy profile indicates that colloids can be spontaneously detached from the shallow primary energy well by Brownian diffusion. The spontaneous detachment from primary minima was unambiguously confirmed by conducting laboratory column transport experiments involving flow interruptions for two model colloids (polystyrene latex microspheres) and engineered nanoparticles (fullerene C₆₀ aggregates). Whereas the spontaneous detachment has been frequently attributed to attachment in secondary minima in the literature, our study indicates that the detached colloids could be initially attached at primary minima. Our study further suggests that the spontaneous disaggregation from primary minima is more significant than spontaneous detachment because the primary minimum depth between colloid themselves is lower than that between a colloid and a collector surface.     

Selection Shapes the Robustness of Ligand-Binding Amino Acids

The phenotypes of biological systems are to some extent robust to genotypic changes. Such robustness exists on multiple levels of biological organization. We analyzed this robustness for two categories of amino acids in proteins. Specifically, we studied the codons of amino acids that bind or do not bind small molecular ligands. We asked to what extent codon changes caused by mutation or mistranslation may affect physicochemical amino acid properties or protein folding. We found that the codons of ligand-binding amino acids are on average more robust than those of non-binding amino acids. Because mistranslation is usually more frequent than mutation, we speculate that selection for error mitigation at the translational level stands behind this phenomenon. Our observations suggest that natural selection can affect the robustness of very small units of biological organization.     

The Anomalous Diffusion of a Tumor Invading with Different Surrounding Tissues

We simulated the invasion of a proliferating, diffusing tumor within different surrounding tissue conditions using a hybrid mathematical model. The in silico invasion of a tumor was addressed systematically for the first time within the framework of a generalized diffusion theory. Our results reveal that a tumor not only migrates using typical Fickian diffusion, but also migrates more generally using subdiffusion, superdiffusion, and even ballistic diffusion, with increasing mobility of the tumor cell when haptotaxis and chemotaxis toward the host tissue surrounding the proliferative tumor are involved. Five functional terms were included in the hybrid model and their effects on a tumor''s invasion were investigated quantitatively: haptotaxis toward the extracellular matrix tissue that is degraded by matrix metalloproteinases; chemotaxis toward nutrients; cell-cell adhesion; the proliferation of the tumor; and the immune response toward the tumor. Haptotaxis and chemotaxis, which are initiated by extracellular matrix and nutrient supply (i.e., glucose) respectively, as well as cell-cell adhesions all drastically affect a tumor''s diffusion mode when a tumor invades its surrounding host tissue and proliferates. We verified the in silico invasive behavior of a tumor by analyzing experimental data gathered from the in vitro culturing of different tumor cells and clinical imaging observations that used the same approach as was used to process the simulation data. The different migration modes of a tumor suggested by the simulations generally conform to the results observed in cell cultures and in clinical imaging. Our study not only discloses some migration modes of a tumor that proliferates and invades under different host tissues conditions, but also provides a heuristic method to characterize the invasion of a tumor in clinical medical imaging analysis.     

High Conformational Variability in the GluK2 Kainate Receptor Ligand-Binding Domain

1. Download high-res image (277KB)
2. Download full-size image

     

Defining a Canonical Ligand-Binding Pocket in the Orphan Nuclear Receptor Nurr1

1. Download high-res image (248KB)
2. Download full-size image

     

Anomalous Diffusion Reports on the Interaction of Misfolded Proteins with the Quality Control Machinery in the Endoplasmic Reticulum

A multitude of transmembrane proteins enters the endoplasmic reticulum (ER) as unfolded polypeptide chains. During their folding process, they interact repetitively with the ER's quality control machinery. Here, we have used fluorescence correlation spectroscopy to probe these interactions for a prototypical transmembrane protein, VSVG ts045, in vivo. While both folded and unfolded VSVG ts045 showed anomalous diffusion, the unfolded protein had a significantly stronger anomaly. This difference subsided when unfolded VSVG ts045 was in a complex with its chaperone calnexin, or when a mutant form of VSVG ts045 with only one glycan was used. Our experimental data and accompanying simulations suggest that the folding sensor of the quality control (UGT1) oligomerizes unfolded VSVG ts045, leading to a more anomalous/obstructed diffusion. In contrast, calnexin dissolves the oligomers, rendering unfolded VSVG ts045 more mobile, and hence prevents poisoning of the ER.     

A New Coarse-Grained Model for E. coli Cytoplasm: Accurate Calculation of the Diffusion Coefficient of Proteins and Observation of Anomalous Diffusion

A new coarse-grained model of the E. coli cytoplasm is developed by describing the proteins of the cytoplasm as flexible units consisting of one or more spheres that follow Brownian dynamics (BD), with hydrodynamic interactions (HI) accounted for by a mean-field approach. Extensive BD simulations were performed to calculate the diffusion coefficients of three different proteins in the cellular environment. The results are in close agreement with experimental or previously simulated values, where available. Control simulations without HI showed that use of HI is essential to obtain accurate diffusion coefficients. Anomalous diffusion inside the crowded cellular medium was investigated with Fractional Brownian motion analysis, and found to be present in this model. By running a series of control simulations in which various forces were removed systematically, it was found that repulsive interactions (volume exclusion) are the main cause for anomalous diffusion, with a secondary contribution from HI.