skip to main content


Title: The combined hydrodynamic and thermodynamic effects of immobilized proteins on the diffusion of mobile transmembrane proteins
The plasma membranes of cells are thin viscous sheets in which some transmembrane proteins have two-dimensional mobility and some are immobilized. Previous studies have shown that immobile proteins retard the short-time diffusivity of mobile particles through hydrodynamic interactions and that steric effects of immobile proteins reduce the long-time diffusivity in a model that neglects hydrodynamic interactions. We present a rigorous derivation of the long-time diffusivity of a single mobile protein interacting hydrodynamically and thermodynamically with an array of immobile proteins subject to periodic boundary conditions. This method is based on a finite element method (FEM) solution of the probability density of the mobile protein diffusing with a position-dependent mobility determined through a multipole solution of Stokes equations. The simulated long-time diffusivity in square arrays decreases as the spacing in the array approaches the particle size in a manner consistent with a lubrication analysis. In random arrays, steric effects lead to a percolation threshold volume fraction above which long-time diffusion is arrested. The FEM/multipole approach is used to compute the long-time diffusivity far away from this threshold. An approximate analysis of mobile protein diffusion through a network of pores connected by bonds with resistances determined by the FEM/multipole calculations is then used to explore higher immobile area fractions and to evaluate the finite simulation cell size scaling behaviour of diffusion near the percolation threshold. Surprisingly, the ratio of the long-time diffusivity to the spatially averaged short-time diffusivity in these two-dimensional fixed arrays is higher in the presence of hydrodynamic interactions than in their absence. Finally, the implications of this work are discussed, including the possibility of using the methods developed here to investigate more complex diffusive phenomena observed in cell membranes.  more » « less
Award ID(s):
1803156
NSF-PAR ID:
10162407
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
877
ISSN:
0022-1120
Page Range / eLocation ID:
648 to 681
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The zeta potential (ζ) is the effective charge energy of a solvated protein, describing the magnitude of electrostatic interactions in solution. It is commonly used in the assessment of adsorption processes and dispersion stability. Predicting ζ from molecular structure would be useful to the structure‐based molecular design of drugs, proteins, and other molecules that hold charge‐dependent function while remaining suspended in solution. One challenge in predicting ζ is identifying the location of the slip plane (XSP), a distance from the protein surface where ζ is theoretically defined. This study tests the hypothesis that theXSPcan be estimated by the Stokes–Einstein hydrodynamic radius (Rh), using globular hen egg white lysozyme as a model system. Although theXSPandRhdiffer in their theoretical definitions, with theXSPbeing the position of the ζ during electrokinetic phenomena (e.g., electrophoresis) and theRhbeing a radius pertaining to the edge of solvation during diffusion, they both represent the point where water and ions no longer adhere to a molecule. This work identifies the limited range of ionic strengths in which theXSPcan be determined using diffusivity measurements and the Stokes–Einstein equation. In addition, a computational protocol is developed for determining the ζ from a protein crystal structure. At low ionic strengths, a hyperdiffusivity regime exists, requiring direct measurement of electrophoretic mobility to determine ζ. This work, therefore, supports a basic tenant of EDL theory that the electric double layer during diffusion and electrophoresis are equivalent in the Stokes–Einstein regime.

     
    more » « less
  2. Abstract

    We investigate the effect of steric and hydrodynamic interactions (HI) on quiescent diffusion and flow‐driven transport of finite‐sized nanoparticles through periodic 2D (two‐dimensional) and 3D (three‐dimensional) nanopost arrays using Stokesian dynamics simulations. We find that steric and HI hinder particle diffusivity under quiescent conditions and enhance longitudinal dispersion under flow. Moreover, the presence of HI leads to a power‐law increase in the longitudinal dispersion coefficient with Pe due to spatial variations in the fluid velocity. Lastly, our simulations reveal that longitudinal particle dispersion coefficients behave similarly in 2D and 3D when HI are included.

     
    more » « less
  3. Recent single-molecule measurements [Schoch et al., Proc. Natl. Acad. Sci. U. S. A. 118, e2113202118 (2021)] have observed dynamic lipid–lipid correlations in membranes with submicrometer spatial resolution and submillisecond temporal resolution. While short from an instrumentation standpoint, these length and time scales remain long compared to microscopic molecular motions. Theoretical expressions are derived to infer experimentally measurable correlations from the two-body diffusion matrix appropriate for membrane-bound bodies coupled by hydrodynamic interactions. The temporal (and associated spatial) averaging resulting from finite acquisition times has the effect of washing out correlations as compared to naive predictions (i.e., the bare elements of the diffusion matrix), which would be expected to hold for instantaneous measurements. The theoretical predictions are shown to be in excellent agreement with Brownian dynamics simulations of experimental measurements. Numerical results suggest that the experimental measurement of membrane protein diffusion, in complement to lipid diffusion measurements, might help to resolve the experimental ambiguities encountered for certain black lipid membranes. 
    more » « less
  4. null (Ed.)
    Autonomous active, elastic filaments that interact with each other to achieve cooperation and synchrony underlie many critical functions in biology. The mechanisms underlying this collective response and the essential ingredients for stable synchronization remain a mystery. Inspired by how these biological entities integrate elasticity with molecular motor activity to generate sustained oscillations, a number of synthetic active filament systems have been developed that mimic oscillations of these biological active filaments. Here, we describe the collective dynamics and stable spatiotemporal patterns that emerge in such biomimetic multi-filament arrays, under conditions where steric interactions may impact or dominate the collective dynamics. To focus on the role of steric interactions, we study the system using Brownian dynamics, without considering long-ranged hydrodynamic interactions. The simulations treat each filament as a connected chain of self-propelling colloids. We demonstrate that short-range steric inter-filament interactions and filament roughness are sufficient – even in the absence of inter-filament hydrodynamic interactions – to generate a rich variety of collective spatiotemporal oscillatory, traveling and static patterns. We first analyze the collective dynamics of two- and three-filament clusters and identify parameter ranges in which steric interactions lead to synchronized oscillations and strongly occluded states. Generalizing these results to large one-dimensional arrays, we find rich emergent behaviors, including traveling metachronal waves, and modulated wavetrains that are controlled by the interplay between the array geometry, filament activity, and filament elasticity. Interestingly, the existence of metachronal waves is non-monotonic with respect to the inter-filament spacing. We also find that the degree of filament roughness significantly affects the dynamics – specifically, filament roughness generates a locking-mechanism that transforms traveling wave patterns into statically stuck and jammed configurations. Taken together, simulations suggest that short-ranged steric inter-filament interactions could combine with complementary hydrodynamic interactions to control the development and regulation of oscillatory collective patterns. Furthermore, roughness and steric interactions may be critical to the development of jammed spatially periodic states; a spatiotemporal feature not observed in purely hydrodynamically interacting systems. 
    more » « less
  5. Atomistic molecular dynamics simulations of concentrated protein solutions in the presence of a phospholipid bilayer are presented to gain insights into the dynamics and interactions at the cytosol–membrane interface. The main finding is that proteins that are not known to specifically interact with membranes are preferentially excluded from the membrane, leaving a depletion zone near the membrane surface. As a consequence, effective protein concentrations increase, leading to increased protein contacts and clustering, whereas protein diffusion becomes faster near the membrane for proteins that do occasionally enter the depletion zone. Since protein–membrane contacts are infrequent and short-lived in this study, the structure of the lipid bilayer remains largely unaffected by the crowded protein solution, but when proteins do contact lipid head groups, small but statistically significant local membrane curvature is induced, on average. 
    more » « less