Resolving Cytosolic Diffusive States in Bacteria by Single-Molecule Tracking

The trajectory of a single protein in the cytosol of a living cell contains information about its molecular interactions in its native environment. However, it has remained challenging to accurately resolve and characterize the diffusive states that can manifest in the cytosol using analytical approaches based on simplifying assumptions. Here, we show that multiple intracellular diffusive states can be successfully resolved if sufficient single-molecule trajectory information is available to generate well-sampled distributions of experimental measurements and if experimental biases are taken into account during data analysis. To address the inherent experimental biases in camera-based and MINFLUX-based single-molecule tracking, we use an empirical data analysis framework based on Monte Carlo simulations of confined Brownian motion. This framework is general and adaptable to arbitrary cell geometries and data acquisition parameters employed in two-dimensional or three-dimensional single-molecule tracking. We show that, in addition to determining the diffusion coefficients and populations of prevalent diffusive states, the timescales of diffusive state switching can be determined by stepwise increasing the time window of averaging over subsequent single-molecule displacements. Time-averaged diffusion analysis of single-molecule tracking data may thus provide quantitative insights into binding and unbinding reactions among rapidly diffusing molecules that are integral for cellular functions.     

An Intermittent Model for Intracellular Motions of Gold Nanostars by k-Space Scattering Image Correlation

Anisotropic metallic nanoparticles have been devised as powerful potential tools for in vivo imaging, photothermal therapy, and drug delivery thanks to plasmon-enhanced absorption and scattering cross sections, ease in synthesis and functionalization, and controlled cytotoxicity. The rational design of all these applications requires the characterization of the nanoparticles intracellular trafficking pathways. In this work, we exploit live-cell time-lapse confocal reflectance microscopy and image correlation in both direct and reciprocal space to investigate the intracellular transport of branched gold nanostars (GNSs). Different transport mechanisms, spanning from pure Brownian diffusion to (sub-)ballistic superdiffusion, are revealed by temporal and spatio-temporal image correlation spectroscopy on the tens-of-seconds timescale. According to these findings, combined with numerical simulations and with a Bayesian (hidden Markov model-based) analysis of single particle tracking data, we ascribe the superdiffusive, subballistic behavior characterizing the GNSs dynamics to a two-state switching between Brownian diffusion in the cytoplasm and molecular motor-mediated active transport. For the investigation of intermittent-type transport phenomena, we derive an analytical theoretical framework for Fourier-space image correlation spectroscopy (kICS). At first, we evaluate the influence of all the dynamic and kinetic parameters (the diffusion coefficient, the drift velocity, and the transition rates between the diffusive and the active transport regimes) on simulated kICS correlation functions. Then we outline a protocol for data analysis and employ it to derive whole-cell maps for each parameter underlying the GNSs intracellular dynamics. Capable of identifying even simpler transport phenomena, whether purely diffusive or ballistic, our intermittent kICS approach allows an exhaustive investigation of the dynamics of GNSs and biological macromolecules.     

Local motion analysis reveals impact of the dynamic cytoskeleton on intracellular subdiffusion

Intracellular transport is a complex interplay of ballistic transport along filaments and of diffusive motion, reliably delivering material and allowing for cell differentiation, migration, and proliferation. The diffusive regime, including subdiffusive, Brownian, and superdiffusive motion, is of particular interest for inferring information about the dynamics of the cytoskeleton morphology during intracellular transport. The influence of dynamic cytoskeletal states on intracellular transport are investigated in Dictyostelium discoideum cells by single particle tracking of fluorescent nanoparticles, to relate quantitative motion parameters and intracellular processes before and after cytoskeletal disruption. A local mean-square displacement (MSD) analysis separates ballistic motion phases, which we exclude here, from diffusive nanoparticle motion. In this study, we focus on intracellular subdiffusion and elucidate lag-time dependence, with particular focus on the impact of cytoskeleton compartments like microtubules and actin filaments. This method proves useful for binary motion state distributions. Experimental results are compared to simulations of a data-driven Langevin model with finite velocity correlations that captures essential statistical features of the local MSD algorithm. Specifically, the values of the mean MSD exponent and effective diffusion coefficients can be traced back to negative correlations of the motion's increments. We clearly identify both microtubules and actin filaments as the cause for intracellular subdiffusion and show that actin-microtubule cross talk exerts viscosifying effects at timescales larger than 0.2 s. Our findings might give insights into material transport and information exchange in living cells, which might facilitate gaining control over cell functions.     

A Method for Molecular Dynamics on Curved Surfaces

Dynamics simulations of constrained particles can greatly aid in understanding the temporal and spatial evolution of biological processes such as lateral transport along membranes and self-assembly of viruses. Most theoretical efforts in the field of diffusive transport have focused on solving the diffusion equation on curved surfaces, for which it is not tractable to incorporate particle interactions even though these play a crucial role in crowded systems. We show here that it is possible to take such interactions into account by combining standard constraint algorithms with the classical velocity Verlet scheme to perform molecular dynamics simulations of particles constrained to an arbitrarily curved surface. Furthermore, unlike Brownian dynamics schemes in local coordinates, our method is based on Cartesian coordinates, allowing for the reuse of many other standard tools without modifications, including parallelization through domain decomposition. We show that by applying the schemes to the Langevin equation for various surfaces, we obtain confined Brownian motion, which has direct applications to many biological and physical problems. Finally we present two practical examples that highlight the applicability of the method: 1) the influence of crowding and shape on the lateral diffusion of proteins in curved membranes; and 2) the self-assembly of a coarse-grained virus capsid protein model.     

Particle Transport through Hydrogels Is Charge Asymmetric

Transport processes within biological polymer networks, including mucus and the extracellular matrix, play an important role in the human body, where they serve as a filter for the exchange of molecules and nanoparticles. Such polymer networks are complex and heterogeneous hydrogel environments that regulate diffusive processes through finely tuned particle-network interactions. In this work, we present experimental and theoretical studies to examine the role of electrostatics on the basic mechanisms governing the diffusion of charged probe molecules inside model polymer networks. Translational diffusion coefficients are determined by fluorescence correlation spectroscopy measurements for probe molecules in uncharged as well as cationic and anionic polymer solutions. We show that particle transport in the charged hydrogels is highly asymmetric, with diffusion slowed down much more by electrostatic attraction than by repulsion, and that the filtering capability of the gel is sensitive to the solution ionic strength. Brownian dynamics simulations of a simple model are used to examine key parameters, including interaction strength and interaction range within the model networks. Simulations, which are in quantitative agreement with our experiments, reveal the charge asymmetry to be due to the sticking of particles at the vertices of the oppositely charged polymer networks.     

Analysis of Microscopic Parameters of Single-Particle Trajectories in Neurons

We performed a comparative study of the statistical uncertainties that arise when calculating the velocity and diffusion coefficients from single-particle trajectories. We show that a method where particle mean displacement is used to calculate velocity and mean square fluctuation is used to calculate diffusion coefficient offers greater accuracy than analysis of time-dependent mean square displacement. Our assessment of the performance of the two analysis strategies is conducted in two ways. First, we apply each of the methods to simulated trajectories where each parameter term is known. Second, we analyze the motion of previously uncharacterized EphB2 receptors in the membrane of hippocampal neurons. We find that EphB2 receptors display different types of motion mode and transition between these modes. We present our data as a distribution of microscopic diffusion coefficients for each particle trajectory, which we refer to as partial distributions. Partial distributions are summed to form a cumulative distribution of diffusion coefficients for EphB2 receptors in hippocampal neurons. The structure and interpretation of the EphB2 cumulative distribution are discussed.     

Tracer diffusion through F-actin: effect of filament length and cross-linking.

We have determined diffusion coefficients for small (50- to 70-nm diameter) fluorescein-thiocarbamoyl-labeled Ficoll tracers through F-actin as a function of filament length and cross-linking. fx45 was used to regulate filament length and avidin/biotinylated actin or ABP-280 was used to prepare cross-linked actin gels. We found that tracer diffusion was generally independent of filament length in agreement with theoretical predictions for diffusion through solutions of rods. However, in some experiments diffusion was slower through short (< or = 1.0 micron) filaments, although this result was not consistently reproducible. Measured diffusion coefficients through unregulated F-actin and filaments of lengths > 1.0 micron were more rapid than predicted by theory for tracer diffusion through rigid, random networks, which was consistent with some degree of actin bundling. Avidin-induced cross-linking of biotinylated F-actin did not affect diffusion through unregulated F-actin, but in cases where diffusion was slower through short filaments this cross-linking method resulted in enhanced tracer diffusion rates indistinguishable from unregulated F-actin. This finding, in conjunction with increased turbidity of 1.0-micron filaments upon avidin cross-linking, indicated that this cross-linking method induces F-actin bundling. By contrast, ABP-280 cross-linking retarded diffusion through unregulated F-actin and decreased turbidity. Tracer diffusion under these conditions was well approximated by the diffusion theory. Both cross-linking procedures resulted in gel formation as determined by falling ball viscometry. These results demonstrate that network microscopic geometry is dependent on the cross-linking method, although both methods markedly increase F-actin macroscopic viscosity.     

Brownian dynamics simulation of knot diffusion along a stretched DNA molecule

Manipulation of individual DNA molecules by optical tweezers has made it possible to tie these molecules into knots. After stretching the DNA molecules the knots become highly localized. In their recent study, Quake and co-authors investigated diffusion of such knots along stretched DNA molecules. We used these data to test the accuracy of a Brownian dynamics simulation of DNA bending motion. We simulated stretched DNA molecules with knots 3(1), 4(1), and 7(1), and determined their diffusion coefficients. Comparison of the simulated and experimental results shows that Brownian dynamics simulation is capable of predicting the rates of large-scale DNA rearrangements within a factor of 2.     

Characterization of the diffusive properties of biofilms using pulsed field gradient-nuclear magnetic resonance

The mobility of water in intact biofilms was measured with pulsed field gradient nuclear magnetic resonance (PFG-NMR) and used to characterise their diffusive properties. The results obtained with several well-defined systems, viz. pure water, agar, and agar containing inert particles or active bacteria were compared to glucose diffusion coefficients measured with micro-electrodes and those calculated utilising theoretical diffusion models. A good correspondence was observed indicating that PFG-NMR should also enable the measurement of diffusion coefficients in heterogeneous biological systems. Diffusion coefficients of several types of natural biofilms were measured as well and these results were related to the physical biofilm characteristics. The values had a high accuracy and reflected the properties of a sample of ca. 100 biofilms, while non-uniformity or non-geometrical shapes did not negatively influence the results. The monitored PFG-NMR signal contains supplementary information on e.g. cell fraction or spatial organisation but quantitative analysis was not yet possible. Copyright 1998 John Wiley & Sons, Inc.     

Multiple modes in a coral species abundance distribution

Species abundance distributions are an important measure of biodiversity and community structure. These distributions are affected by sampling, and alternative species-abundance models often make similar predictions for small sample sizes. Very large samples reveal the relative abundances of rare species, and thus provide information about species relative abundances that small samples cannot. Here, we present the species-abundance distribution for a sample of > 40,000 coral colonies at a single site, exceeding existing samples of coral local assemblages by over an order of magnitude. This abundance distribution is multimodal when examined on a logarithmic scale. Four different model selection procedures all indicate that the underlying community abundance distribution has at least three modes. We show that the multiple modes are not caused by mixtures of species with different habitat preferences. However, spatial aggregation partially explains our results. We inspect published work on species abundance distributions, and suggest that multimodality may be a common feature of large samples.     

Diffusion on Membrane Domes,Tubes, and Pearling Structures

Diffusion is a fundamental mechanism for protein distribution in cell membranes. These membranes often exhibit complex shapes, which range from shallow domes to elongated tubular or pearl-like structures. Shape complexity of the membrane influences the diffusive spreading of proteins and molecules. Despite the importance membrane geometry plays in these diffusive processes, it is challenging to establish the dependence between diffusion and membrane morphology. We solve the diffusion equation numerically on various static curved shapes representative for experimentally observed membrane shapes. Our results show that membrane necks become diffusion barriers. We determine the diffusive half-time, i.e., the time that is required to reduce the amount of protein in the budded region by one half, and find a quadratic relation between the diffusive half-time and the averaged mean curvature of the membrane shape, which we rationalize by a scaling law. Our findings thus help estimate the characteristic diffusive timescale based on the simple measure of membrane mean curvature.     

Directionality of stripes formed by anisotropic reaction-diffusion models

Turing mechanism explains the formation of striped patterns in a uniform field in which two substances interact locally and diffuse randomly. In a twin paper, to explain the directionality of stripes on fish skin in closely related species, we studied the effect of anisotropic diffusion of the two substances on the direction of stripes, in the cases in which both substances have high diffusivity in the same direction. In this paper, we study the direction of stripes in more general situations in which the diffusive direction may differ between the two substances. We derive a formula for the direction of stripes, based on a heuristic argument of unstable modes of deviation from the uniform steady state. We confirm the accuracy of the formula by computer simulations. When the diffusive direction is different between two substances, the directions of stripes in the spatial pattern change smoothly with the magnitude of anisotropy of two substances. When the diffusive direction of the two substances is the same, the stripes are formed either parallel or perpendicular to the common diffusive direction, depending on the relative magnitude of the anisotropy. The transition between these two phases occurs sharply.     

Size-Dependent Diffusion of Dextrans in Excised Porcine Corneal Stroma

Delivery of therapeutic agents to the eye requires efficient transport through cellular and extracellular barriers. We evaluated the rate of diffusive transport in excised porcine corneal stroma using fluorescently labeled dextran molecules with hydrodynamic radii ranging from 1.3 to 34 nm. Fluorescence correlation spectroscopy (FCS) was used to measure diffusion coefficients of dextran molecules in the excised porcine corneal stroma. The preferential sensitivity of FCS to diffusion along two dimensions was used to differentially probe diffusion along the directions parallel to and perpendicular to the collagen lamellae of the corneal stroma. In order to develop an understanding of how size affects diffusion in cornea, diffusion coefficients in cornea were compared to diffusion coefficients measured in a simple buffer solution. Dextran molecules diffuse more slowly in cornea as compared to buffer solution. The reduction in diffusion coefficient is modest however (67% smaller), and is uniform over the range of sizes that we measured. This indicates that, for dextrans in the 1.3 to 34 nm range, the diffusion landscape of corneal stroma can be represented as a simple liquid with a viscosity approximately 1.5 times that of water. Diffusion coefficients measured parallel vs. perpendicular to the collagen lamellae were indistinguishable. This indicates that diffusion in the corneal stroma is not highly anisotropic. Our results support the notion that the corneal stroma is highly permeable and isotropic to transport of hydrophilic molecules and particles with hydrodynamic radii up to at least 34 nm.     

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.     

Raster image correlation spectroscopy in live cells

Raster image correlation spectroscopy (RICS) is a noninvasive technique to detect and quantify events in a live cell, including concentration of molecules and diffusion coefficients of molecules; in addition, by measuring changes in diffusion coefficients, RICS can indirectly detect binding. Any specimen containing fluorophores that can be imaged with a laser scanning microscope can be analyzed using RICS. There are other techniques to measure diffusion coefficients and binding; however, RICS fills a unique niche. It provides spatial information and can be performed in live cells using a conventional confocal microscope. It can measure a range of diffusion coefficients that is not accessible with any other single optical correlation-based technique. In this article we describe a protocol to obtain raster scanned images with an Olympus FluoView FV1000 confocal laser scanning microscope using Olympus FluoView software to acquire data and SimFCS software to perform RICS analysis. Each RICS measurement takes several minutes. The entire procedure can be completed in ～2 h. This procedure includes focal volume calibration using a solution of fluorophores with a known diffusion coefficient and measurement of the diffusion coefficients of cytosolic enhanced green fluorescent protein (EGFP) and EGFP-paxillin.     

Monte Carlo analysis of obstructed diffusion in three dimensions: application to molecular diffusion in organelles.

Molecular transport in the aqueous lumen of organelles involves diffusion in a confined compartment with complex geometry. Monte Carlo simulations of particle diffusion in three dimensions were carried out to evaluate the influence of organelle structure on diffusive transport and to relate experimental photobleaching data to intrinsic diffusion coefficients. Two organelle structures were modeled: a mitochondria-like long closed cylinder containing fixed luminal obstructions of variable number and size, and an endoplasmic reticulum-like network of interconnected cylinders of variable diameter and density. Trajectories were computed in each simulation for >10(5) particles, generally for >10(5) time steps. Computed time-dependent concentration profiles agreed quantitatively with analytical solutions of the diffusion equation for simple geometries. For mitochondria-like cylinders, significant slowing of diffusion required large or wide single obstacles, or multiple obstacles. In simulated spot photobleaching experiments, a approximately 25% decrease in apparent diffusive transport rate (defined by the time to 75% fluorescence recovery) was found for a single thin transverse obstacle occluding 93% of lumen area, a single 53%-occluding obstacle of width 16 lattice points (8% of cylinder length), 10 equally spaced 53% obstacles alternately occluding opposite halves of the cylinder lumen, or particle binding to walls (with mean residence time = 10 time steps). Recovery curve shape with obstacles showed long tails indicating anomalous diffusion. Simulations also demonstrated the utility of measurement of fluorescence depletion at a spot distant from the bleach zone. For a reticulum-like network, particle diffusive transport was mildly reduced from that in unobstructed three-dimensional space. In simulated photobleaching experiments, apparent diffusive transport was decreased by 39-60% in reticular structures in which 90-97% of space was occluded. These computations provide an approach to analyzing photobleaching data in terms of microscopic diffusive properties and support the paradigm that organellar barriers must be quite severe to seriously impede solute diffusion.