Many of the cell membrane's vital functions are achieved by the self-organization of the proteins and biopolymers embedded in it. The protein dynamics is in part determined by its drag. A large number of these proteins can polymerize to form filaments.
Hydrodynamics of a single filament moving in a fluid spherical membrane
Dynamic organization of the cytoskeletal filaments and rod-like proteins in the cell membrane and other biological interfaces occurs in many cellular processes. Previous modeling studies have considered the dynamics of a single rod on fluid planar membranes. We extend these studies to the more physiologically relevant case of a single filament moving in a spherical membrane. Specifically, we use a slender-body formulation to compute the translational and rotational resistance of a single filament of length L moving in a membrane of radius R and 2D viscosity ηm, and surrounded on its interior and exterior with Newtonian fluids of viscosities η− and η+. We first discuss the case where the filament's curvature is at its minimum κ=1/R. We show that the boundedness of spherical geometry gives rise to flow confinement effects that increase in strength with increasing the ratio of filament's length to membrane radius L/R. These confinement flows only result in a mild increase in filament's resistance along its axis, ξ∥, and its rotational resistance, ξΩ. As a result, our predictions of ξ∥ and ξΩ can be quantitatively mapped to the results on a planar membrane. In contrast, we find that the drag in perpendicular direction, ξ⊥, increases superlinearly with the filament's length, when L/R>1 and ultimately ξ⊥→∞ as L/R→π. Next, we consider the effect of the filament's curvature, κ, on its parallel motion, while fixing the membrane's radius. We show that the flow around the filament becomes increasingly more asymmetric with increasing its curvature. These flow asymmetries induce a net torque on the filament, coupling its parallel and rotational dynamics. This coupling becomes stronger with increasing L/R and κ.
more »
« less
- Award ID(s):
- 1944156
- NSF-PAR ID:
- 10318632
- Date Published:
- Journal Name:
- ArXivorg
- ISSN:
- 2331-8422
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
In vitro studies of protein–membrane interactions often involve using rigid beads coated with lipid bilayers, as a model for the cell membrane. Motivated by this, we use slender-body theory to compute the translational and rotational resistance of a single filamentous protein embedded in the outer layer of a supported bilayer membrane and surrounded on the exterior by a Newtonian fluid. We first consider the regime where the two layers are strongly coupled through their inter-leaflet friction. We find that the drag along the parallel direction grows linearly with the filament's length and quadratically with the length for the perpendicular and rotational drag coefficients. These findings are explained using scaling arguments and by analysing the velocity fields around the moving filament. We then present and discuss the qualitative differences between the drag of a filament moving in a freely suspended bilayer and a supported membrane as a function of the membrane's inter-leaflet friction. Finally, we briefly discuss how these findings can be used in experiments to determine membrane rheology. In summary, we present a formulation that allows computation of the effects of membrane properties (its curvature, viscosity and inter-leaflet friction), and the exterior and interior three-dimensional fluids’ depth and viscosity on the drag of a rod-like/filamentous protein, all in a unified theoretical framework. -
more » « less
The shape assumed by a slender elastic structure is a function both of the geometry of the space in which it exists and the forces it experiences. We explore, by experiments and theoretical analysis, the morphological phase space of a filament confined to the surface of a spherical bubble. The morphology is controlled by varying bending stiffness and weight of the filament, and its length relative to the bubble radius. When the dominant considerations are the geometry of confinement and elastic energy, the filament lies along a geodesic and when gravitational energy becomes significant, a bifurcation occurs, with a part of the filament occupying a longitude and the rest along a curve approximated by a latitude. Far beyond the transition, when the filament is much longer than the diameter, it coils around the selected latitudinal region. A simple model with filament shape as a composite of two arcs captures the transition well. For better quantitative agreement with the subcritical nature of bifurcation, we study the morphology by numerical energy minimization. Our analysis of the filament’s morphological space spanned by one geometric parameter, and one parameter that compares elastic energy with body forces, may provide guidance for packing slender structures on complex surfaces.
-
Using a combination of theory and experiments we study the interface between two immiscible domains in a colloidal membrane composed of rigid rods of different lengths. Geometric considerations of rigid rod packing imply that a domain of sufficiently short rods in a background membrane of long rods is more susceptible to twist than the inverse structure, a long-rod domain in a short-rod membrane background. The tilt at the inter-domain edge forces splay, which in turn manifests as a spontaneous edge curvature whose energetics are controlled by the length asymmetry of constituent rods. A thermodynamic model of such tilt-curvature coupling at inter-domain edges explains a number of experimental observations, including a non-monotonic dependence of the edge twist on the domain radius, and annularly shaped domains of long rods. Our work shows how coupling between orientational and compositional degrees of freedom in two-dimensional fluids give rise to complex shapes and thermodynamics of domains, analogous to shape transitions in 3D fluid vesicles.more » « less
-
We perform direct numerical simulations of a gas bubble dissolving in a surrounding liquid. The bubble volume is reduced due to dissolution of the gas, with the numerical implementation of an immersed boundary method, coupling the gas diffusion and the Navier–Stokes equations. The methods are validated against planar and spherical geometries’ analytical moving boundary problems, including the classic Epstein–Plesset problem. Considering a bubble rising in a quiescent liquid, we show that the mass transfer coefficient $k_L$ can be described by the classic Levich formula $k_L = (2/\sqrt {{\rm \pi} })\sqrt {\mathscr {D}_l\,U(t)/d(t)}$ , with $d(t)$ and $U(t)$ the time-varying bubble size and rise velocity, and $\mathscr {D}_l$ the gas diffusivity in the liquid. Next, we investigate the dissolution and gas transfer of a bubble in homogeneous and isotropic turbulence flow, extending Farsoiya et al. ( J. Fluid Mech. , vol. 920, 2021, A34). We show that with a bubble size initially within the turbulent inertial subrange, the mass transfer coefficient in turbulence $k_L$ is controlled by the smallest scales of the flow, the Kolmogorov $\eta$ and Batchelor $\eta _B$ microscales, and is independent of the bubble size. This leads to the non-dimensional transfer rate ${Sh}=k_L L^\star /\mathscr {D}_l$ scaling as ${Sh}/{Sc}^{1/2} \propto {Re}^{3/4}$ , where ${Re}$ is the macroscale Reynolds number ${Re} = u_{rms}L^\star /\nu _l$ , with $u_{rms}$ the velocity fluctuations, $L^*$ the integral length scale, $\nu _l$ the liquid viscosity, and ${Sc}=\nu _l/\mathscr {D}_l$ the Schmidt number. This scaling can be expressed in terms of the turbulence dissipation rate $\epsilon$ as ${k_L}\propto {Sc}^{-1/2} (\epsilon \nu _l)^{1/4}$ , in agreement with the model proposed by Lamont & Scott ( AIChE J. , vol. 16, issue 4, 1970, pp. 513–519) and corresponding to the high $Re$ regime from Theofanous et al. ( Intl J. Heat Mass Transfer , vol. 19, issue 6, 1976, pp. 613–624).more » « less
-
Slender fibres are ubiquitous in biology, physics and engineering, with prominent examples including bacterial flagella and cytoskeletal fibres. In this setting, slender body theories (SBTs), which give the resistance on the fibre asymptotically in its slenderness $\epsilon$ , are useful tools for both analysis and computations. However, a difficulty arises when accounting for twist and cross-sectional rotation: because the angular velocity of a filament can vary depending on the order of magnitude of the applied torque, asymptotic theories must give accurate results for rotational dynamics over a range of angular velocities. In this paper, we first survey the challenges in applying existing SBTs, which are based on either singularity or full boundary integral representations, to rotating filaments, showing in particular that they fail to consistently treat rotation–translation coupling in curved filaments. We then provide an alternative approach which approximates the three-dimensional dynamics via a one-dimensional line integral of Rotne–Prager–Yamakawa regularized singularities. While unable to accurately resolve the flow field near the filament, this approach gives a grand mobility with symmetric rotation–translation and translation–rotation coupling, making it applicable to a broad range of angular velocities. To restore fidelity to the three-dimensional filament geometry, we use our regularized singularity model to inform a simple empirical equation which relates the mean force and torque along the filament centreline to the translational and rotational velocity of the cross-section. The single unknown coefficient in the model is estimated numerically from three-dimensional boundary integral calculations on a rotating, curved filament.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- The DOI is not currently available.
