Binding of thermalized and active membrane curvature-inducing proteins

Which Coverages of Arc-Shaped Proteins Are Required for Membrane Tubulation?

Arc-shaped BIN/Amphiphysin/Rvs (BAR) domain proteins generate curvature by binding to membranes and induce membrane tubulation at sufficiently large protein coverages. For the amphiphysin N-BAR domain, Le Roux et al., Nat. Commun. 2021, 12, 6550, measured a threshold coverage of 0.44 ± 0.097 for nanotubules emerging from the supported lipid bilayer. In this article, we systematically investigate membrane tubulation induced by arc-shaped protein-like particles with coarse-grained modeling and simulations and determine the threshold coverages at different particle-particle interaction strengths and membrane spontaneous curvatures. In our simulations, the binding of arc-shaped particles induces a membrane shape transition from spherical vesicles to tubules at a particle threshold coverage of about 0.5, which is rather robust to variations of the direct attractive particle interactions or spontaneous membrane curvature in the coarse-grained model. Our study suggests that threshold coverages of around or slightly below 0.5 are a general requirement for membrane tubulation by arc-shaped BAR domain proteins.

Curvature sensing of curvature-inducing proteins with internal structure

Many types of peripheral and transmembrane proteins can sense and generate membrane curvature. Laterally isotropic proteins and crescent proteins with twofold rotational symmetry, such as Bin/Amphiphysin/Rvs superfamily proteins, have been studied theoretically. However, proteins often have an asymmetric structure or a higher rotational symmetry. We studied theoretically the curvature sensing of proteins with asymmetric structures and structural deformations. First, we examined proteins consisting of two rodlike segments. When proteins have mirror symmetry, their sensing ability is similar to that of single-rod proteins; hence, with increasing protein density on a cylindrical membrane tube, a second- or first-order transition occurs at a middle or small tube radius, respectively. As asymmetry is introduced, this transition becomes a continuous change and metastable states appear at high protein densities. Protein with threefold, fivefold, or higher rotational symmetry has laterally isotropic bending energy. However, when a structural deformation is allowed, the protein can have a preferred orientation and stronger curvature sensing.

Membrane domain formation induced by binding/unbinding of curvature-inducing molecules on both membrane surfaces

The domain formation of curvature-inducing molecules, such as peripheral or transmembrane proteins and conical surfactants, is studied in thermal equilibrium and nonequilibrium steady states using meshless membrane simulations. These molecules can bind to both surfaces of a bilayer membrane and also move to the opposite leaflet by a flip-flop. Under symmetric conditions for the two leaflets, the membrane domains form checkerboard patterns in addition to striped and spot patterns. The unbound membrane stabilizes the vertices of the checkerboard. Under asymmetric conditions, the domains form kagome-lattice and thread-like patterns. In the nonequilibrium steady states, a flow of the binding molecules between the upper and lower solutions can occur via flip-flop.

Perspective: New directions in dynamical density functional theory

Classical dynamical density functional theory (DDFT) has become one of the central modeling approaches in nonequilibrium soft matter physics. Recent years have seen the emergence of novel and interesting fields of application for DDFT. In particular, there has been a remarkable growth in the amount of work related to chemistry. Moreover, DDFT has stimulated research on other theories such as phase field crystal models and power functional theory. In this perspective, we summarize the latest developments in the field of DDFT and discuss a variety of possible directions for future research.

Excitable reaction-diffusion waves of curvature-inducing proteins on deformable membrane tubes

Living cells employ excitable reaction-diffusion waves for internal cellular functions, in which curvature-inducing proteins are often involved. However, the role of their mechanochemical coupling is not well understood. Here, we report the membrane deformation induced by the excitable reaction-diffusion waves of curvature-inducing proteins and the alternation in the waves due to the deformation, using a coarse-grained simulation of tubular membranes with a modified FitzHugh-Nagumo model. Protein-propagating waves deform tubular membranes and large deformations induce budding and erase waves. The wave speed and shape are determined by a combination of membrane deformation and spatial distribution of the curvature-inducing protein. Waves are also undulated in the azimuthal direction depending on the condition. Rotationally symmetric waves locally deform the tubes into a symmetric shape but maintain a straight shape on average. Our simulation method can be applied to other chemical reaction models and used to investigate various biomembrane phenomena.

The many faces of membrane tension: Challenges across systems and scales

Our understanding of the role of membrane tension in the field of membrane biophysics is rapidly evolving from a passive construct to an active player in a variety of cellular phenomena. Membrane tension has been shown to be a key regulator of many cellular processes ranging including trafficking, ion channel activation, and the invasion of red blood cells by malaria parasites. Recent experimental advances in cells, including the development of a fluorescent tension reporter, have shown that membrane tension is heterogeneous. In this mini-review, I summarize the recent advances in membrane tension measurements and discuss the contributions from different cellular constituents such as the cortical cytoskeleton. Then, I will explore how these different complexities can be considered in biophysical models of different scales. Finally, I will elaborate on the need for iterations between models and experiments as technologies in both fields advance to enable us to obtain critical insights into the physiological role of membrane tension as a critical component of mechanotransduction.

Curvature dependence of BAR protein membrane association and dissociation kinetics

BAR (Bin/Amphiphysin/Rvs) domain containing proteins function as lipid bilayer benders and curvature sensors, and they contribute to membrane shaping involved in cell signaling and metabolism. The mechanism for their membrane shape sensing has been investigated by both equilibrium binding and kinetic studies. In prior research, stopped-flow spectroscopy has been used to deduce a positive dependence on membrane curvature for the binding rate constant, k_on, of a BAR protein called endophilin. However, the impact of bulk diffusion of endophilin, on the kinetic binding parameters has not been thoroughly considered. Employing similar methods, and using lipid vesicles of multiple sizes, we obtained a linear dependence of k_on on vesicle curvature. However, we found that the observed relation can be explained without considering the local curvature sensing ability of endophilin in the membrane association process. In contrast, the diffusion-independent unbinding rate constant (k_off) obtained from stopped-flow measurements shows a negative dependence on membrane curvature, which is controlled/mediated by endophilin-membrane interactions. This latter dependency, in addition to protein-protein interactions on the membrane, explains the selective binding of BAR proteins to highly curved membranes in equilibrium binding experiments.

Binding of anisotropic curvature-inducing proteins onto membrane tubes

Bin/Amphiphysin/Rvs superfamily proteins and other curvature-inducing proteins have anisotropic shapes and anisotropically bend biomembranes. Here, we report how the anisotropic proteins bind the membrane tube and are orientationally ordered using mean-field theory including an orientation-dependent excluded volume. The proteins exhibit a second-order or first-order nematic transition with increasing protein density depending on the radius of the membrane tube. The tube curvatures for the maximum protein binding and orientational order are different and varied by the protein density and rigidity. As the external force along the tube axis increases, a first-order transition from a large tube radius with low protein density to a small radius with high density occurs once, and subsequently, the protein orientation tilts to the tube-axis direction. When an isotropic bending energy is used for the proteins with an elliptic shape, the force-dependence curves become symmetric and the first-order transition occurs twice. This theory quantitatively reproduces the results of meshless membrane simulation for short proteins, whereas deviations are seen for long proteins owing to the formation of protein clusters.

Field-mediated interactions of passive and conformation-active particles: multibody and retardation effects

Particles in soft matter interact through the deformation field they create, as in the "cheerios" effect or the curvature-mediated interactions of membrane proteins. Using a simple model for field-mediated interactions between passive particles, or active particles that switch conformation randomly or synchronously, we derive generic results concerning multibody interactions, activity driven patterns, and retardation effects.

Universal properties of active membranes

We put forward a general field theory for nearly flat fluid membranes with embedded activators and analyze their critical properties using renormalization group techniques. Depending on the membrane-activator coupling, we find a crossover between acoustic and diffusive scaling regimes, with mean-field dynamical critical exponents z=1 and 2, respectively. We argue that the acoustic scaling, which is exact in all spatial dimensions, leads to an early-time behavior, which is representative of the spatiotemporal patterns observed at the leading edge of motile cells, such as oscillations superposed on the growth of the membrane width. In the case of mean-field diffusive scaling, one-loop corrections to the mean-field exponents reveal universal behavior distinct from the Kardar-Parisi-Zhang scaling of passive interfaces and signs of strong-coupling behavior.

Binding of curvature-inducing proteins onto tethered vesicles

A tethered vesicle, which consists of a cylindrical membrane tube and a spherical vesicle, is produced by a mechanical force that is experimentally imposed by optical tweezers and a micropipette. This tethered vesicle is employed for examining the curvature sensing of curvature-inducing proteins. In this study, we clarify how the binding of proteins with a laterally isotropic spontaneous curvature senses and generates the membrane curvatures of the tethered vesicle using mean-field theory and meshless membrane simulation. The force-dependence curves of the protein density in the membrane tube and the tube curvature are reflection symmetric and point symmetric, respectively, from the force point, in which the tube has a sensing curvature. The bending rigidity and spontaneous curvature of the bound proteins can be estimated from these force-dependence curves. First-order transitions can occur between low and high protein densities in the tube at both low and high force amplitudes. The simulation results of the homogeneous phases agree very well with the theoretical predictions. In addition, beaded-necklace-like tubes with microphase separation are found in the simulation.

Thermodynamics and Kinetics of Aggregation of Flexible Peripheral Membrane Proteins

Biomembrane remodeling is essential for cellular trafficking, with membrane-binding peripheral proteins playing a key role in it. Significant membrane remodeling as in endo- and exocytosis is often due to aggregates of many proteins with direct or membrane-mediated interactions. Understanding this process via computer simulations is extremely challenging: protein-membrane systems involve time and length scales that make atomistic simulations impractical, while most coarse-grained models fall short in resolving dynamics and physical effects of protein and membrane flexibility. Here, we develop a coarse-grained model of the bilayer membrane bestrewed with rotationally symmetric flexible proteins, parametrized to reflect local curvatures and lateral dynamics of proteins. We investigate the kinetics, equilibrium distributions, and the free energy landscape governing the formation and breakup of protein clusters on the surface of the membrane. We demonstrate how the flexibility of the proteins as well as their surface concentration play deciding roles in highly selective macroscopic aggregation behavior.

Vesicle budding induced by binding of curvature-inducing proteins

Vesicle budding induced by protein binding that generates an isotropic spontaneous curvature is studied using a mean-field theory. Many spherical buds are formed via protein binding. As the binding chemical potential increases, the proteins first bind to the buds and then to the remainder of the vesicle. For a high spontaneous curvature and/or high bending rigidity of the bound membrane, it is found that a first-order transition occurs between a small number of large buds and a large number of small buds. These two states coexist around the transition point. The proposed scheme is simple and easily applicable to many interaction types, so we investigate the effects of interprotein interactions, the protein-insertion-induced changes in area, the variation of the saddle-splay modulus, and the area-difference-elasticity energy. The differences in the preferred curvatures for curvature sensing and generation are also clarified.

Reaction-diffusion waves coupled with membrane curvature

The reaction-diffusion waves of proteins are known to be involved in fundamental cellular functions, such as cell migration, cell division, and vesicular transportation. In some of these phenomena, pattern formation on the membranes is induced by the coupling between membrane deformation and the reaction-diffusion system through curvature-inducing proteins that bend the biological membranes. Although the membrane shape and the dynamics of the curvature-inducing proteins affect each other in these systems, the effect of such mechanochemical feedback loops on the waves has not been studied in detail. In this study, reaction-diffusion waves coupled with membrane deformation are investigated using simulations combining a dynamically triangulated membrane model with the Brusselator model extended to include the effect of membrane curvature. It is found that the propagating wave patterns change into nonpropageting patterns and spiral wave patterns due to the mechanochemical effects. Moreover, the wave speed is positively or negatively correlated with the local membrane curvature depending on the spontaneous curvature and bending rigidity. In addition, self-oscillation of the vesicle shape occurs, associated with the reaction-diffusion waves of curvature-inducing proteins. This agrees with the experimental observation of GUVs with a reconstituted Min system, which plays a key role in the cell division of Escherichia coli. The findings of this study demonstrate the importance of mechanochemical coupling in biological phenomena.