Small scale membrane mechanics

Large scale changes to lipid bilayer shapes are well represented by the Helfrich model. However, there are membrane processes that take place at smaller length scales that this model cannot address. In this work, we present a one-dimensional continuum model that captures the mechanics of the lipid bilayer membrane at the length scale of the lipids themselves. The model is developed using the Cosserat theory of surfaces with lipid orientation, or ‘tilt’, as the fundamental degree of freedom. The Helfrich model can be recovered as a special case when the curvatures are small and the lipid tilt is everywhere zero. We use the tilt model to study local membrane deformations in response to a protein inclusion. Parameter estimates and boundary conditions are obtained from a coarse-grained molecular model using dissipative particle dynamics (DPD) to capture the same phenomenon. The continuum model is able to reproduce the membrane bending, stretch and lipid tilt as seen in the DPD model. The lipid tilt angle relaxes to the bulk tilt angle within 5–6 nm from the protein inclusion. Importantly, for large tilt gradients induced by the proteins, the tilt energy contribution is larger than the bending energy contribution. Thus, the continuum model of tilt accurately captures behaviors at length scales shorter than the membrane thickness.     

Mechanics of curved plasma membrane vesicles: resting shapes, membrane curvature, and in-plane shear elasticity

Highly curved cell membrane structures, such as plasmalemmal vesicles (caveolae) and clathrin-coated pits, facilitate many cell functions, including the clustering of membrane receptors and transport of specific extracellular macromolecules by endothelial cells. These structures are subject to large mechanical deformations when the plasma membrane is stretched and subject to a change of its curvature. To enhance our understanding of plasmalemmal vesicles we need to improve the understanding of the mechanics in regions of high membrane curvatures. We examine here, theoretically, the shapes of plasmalemmal vesicles assuming that they consist of three membrane domains: an inner domain with high curvature, an outer domain with moderate curvature, and an outermost flat domain, all in the unstressed state. We assume the membrane properties are the same in these domains with membrane bending elasticity as well as in-plane shear elasticity. Special emphasis is placed on the effects of membrane curvature and in-plane shear elasticity on the mechanics of vesicle during unfolding by application of membrane tension. The vesicle shapes were computed by minimization of bending and in-plane shear strain energy. Mechanically stable vesicles were identified with characteristic membrane necks. Upon stretch of the membrane, the vesicle necks disappeared relatively abruptly leading to membrane shapes that consist of curved indentations. While the resting shape of vesicles is predominantly affected by the membrane spontaneous curvatures, the membrane shear elasticity (for a range of values recorded in the red cell membrane) makes a significant contribution as the vesicle is subject to stretch and unfolding. The membrane tension required to unfold the vesicle is sensitive with respect to its shape, especially as the vesicle becomes fully unfolded and approaches a relative flat shape.     

Gating mechanism of mechanosensitive channel of large conductance: a coupled continuum mechanical-continuum solvation approach

Gating transition of the mechanosensitive channel of large conductance (MscL) represents a good example of important biological processes that are difficult to describe using atomistic simulations due to the large (submicron) length scale and long (millisecond) time scale. Here we develop a novel computational framework that tightly couples continuum mechanics with continuum solvation models to study the detailed gating behavior of E. coli-MscL. The components of protein molecules are modeled by continuum elements that properly describe their shape, material properties and physicochemical features (e.g., charge distribution). The lipid membrane is modeled as a three-layer material in which the lipid head group and tail regions are treated separately, taking into account the fact that fluidic lipid bilayers do not bear shear stress. Coupling between mechanical and chemical responses of the channel is realized by an iterative integration of continuum mechanics (CM) modeling and continuum solvation (CS) computation. Compared to previous continuum mechanics studies, the present model is capable of capturing the most essential features of the gating process in a much more realistic fashion: due mainly to the apolar solvation contribution, the membrane tension for full opening of MscL is reduced substantially to the experimental measured range. Moreover, the pore size stabilizes constantly during gating because of the intricate interactions of the multiple components of the system, implying the mechanism for sub-conducting states of MscL gating. A significant fraction (\(\sim \)2/3) of the gating membrane strain is required to reach the first sub-conducting state of our model, which is featured with a relative conductance of 0.115 to the fully opened state. These trends agree well with experimental observations. We anticipate that the coupled CM/CS modeling framework is uniquely suited for the analysis of many biomolecules and their assemblies under external mechanical stimuli.     

A multiscale model for red blood cell mechanics

The objective of this article is the derivation of a continuum model for mechanics of red blood cells via multiscale analysis. On the microscopic level, we consider realistic discrete models in terms of energy functionals defined on networks/lattices. Using concepts of Γ-convergence, convergence results as well as explicit homogenisation formulae are derived. Based on a characterisation via energy functionals, appropriate macroscopic stress–strain relationships (constitutive equations) can be determined. Further, mechanical moduli of the derived macroscopic continuum model are directly related to microscopic moduli. As a test case we consider optical tweezers experiments, one of the most common experiments to study mechanical properties of cells. Our simulations of the derived continuum model are based on finite element methods and account explicitly for membrane mechanics and its coupling with bulk mechanics. Since the discretisation of the continuum model can be chosen freely, rather than it is given by the topology of the microscopic cytoskeletal network, the approach allows a significant reduction of computational efforts. Our approach is highly flexible and can be generalised to many other cell models, also including biochemical control.     

The receptor-mediated endocytosis of nonspherical particles

Enveloped viruses and nanosized biomimetic particles for drug and gene delivery enter target cells mainly through receptor-mediated endocytosis. A few models have been presented to elucidate the mechanics of particle engulfment by the cell membrane, showing how size and surface chemico-physical properties favor or oppose internalization. In this work, the effect of particle nonsphericity is addressed considering elliptical cylindrical particles with aspect ratio Γ. Using a continuum energetic approach, three different conditions have been identified: for sufficiently small Γ, the particle is not even wrapped by the cell membrane; for sufficiently large Γ, the particle is partially wrapped (“frustrated endocytosis”); and for intermediate values of Γ, the particle is fully wrapped and eventually internalized. Given the pleomorphism of viruses and the broad spectrum of shapes for nanosized biomimetic particles, the results presented may be of interest to virologists, pharmacologists, toxicologists, and nanotechnologists.     

Mechanics of stabilized intercellular bridges

Numerous engineered and natural systems form through reinforcement and stabilization of a deformed configuration that was generated by a transient force. An important class of such structures arises during gametogenesis, when a dividing cell undergoes incomplete cytokinesis, giving rise to daughter cells that remain connected through a stabilized intercellular bridge (ICB). ICBs can form through arrest of the contractile cytokinetic furrow and its subsequent stabilization. Despite knowledge of the molecular components, the mechanics underlying robust ICB assembly and the interplay between ring contractility and stiffening are poorly understood. Here, we report joint experimental and theoretical work that explores the physics underlying robust ICB assembly. We develop a continuum mechanics model that reveals the minimal requirements for the formation of stable ICBs, and validate the model’s equilibrium predictions through a tabletop experimental analog. With insight into the equilibrium states, we turn to the dynamics: we demonstrate that contractility and stiffening are in dynamic competition and that the time intervals of their action must overlap to ensure assembly of ICBs of biologically observed proportions. Our results highlight a mechanism in which deformation and remodeling are tightly coordinated—one that is applicable to several mechanics-based applications and is a common theme in biological systems spanning several length scales.     

Formation and maintenance of tubular membrane projections require mechanical force, but their elongation and shortening do not require additional force

Living cells develop their own characteristic shapes depending on their physiological functions, and their morphologies are based on the mechanical characteristics of the cytoskeleton and of membranes. To investigate the role of lipid membranes in morphogenesis, we constructed a simple system that can manipulate liposomes and measure the forces required to transform their shapes. Two polystyrene beads (1 microm in diameter) were encapsulated in giant liposomes and were manipulated using double-beam laser tweezers. Without any specific interaction between the lipid membrane and beads, mechanical forces could be applied to the liposome membrane from the inside. Spherical liposomes transformed into a lemon shape with increasing tension, and tubular membrane projections were subsequently generated in the tips at either end. This process is similar to the liposomal transformation caused by elongation of encapsulated cytoskeletons. In the elongation stage of lemon-shaped liposomes, the force required for the transformation became larger as the end-to-end length increased. Just before the tubular membrane was generated, the force reached the maximum strength (approximately 11 pN). However, immediately after the tubular membrane developed, the force suddenly decreased and was maintained at a constant strength (approximately 4 pN) that was independent of further tube elongation or shortening, even though there was no excess membrane reservoir as occurs in living cells. When the tube length was shortened to approximately 2 microm, the liposome reversed to a lemon shape and the force temporarily increased (to approximately 7 pN). These results indicate that the simple application of mechanical force is sufficient to form a protrusion in a membrane, that a critical force and length is needed to form and to maintain the protrusion, and suggest that the lipid bilayer itself has the ability to buffer the membrane tension.     

Modulation of cellular mechanics during osteogenic differentiation of human mesenchymal stem cells

Recognition of the growing role of human mesenchymal stem cells (hMSC) in tissue engineering and regenerative medicine requires a thorough understanding of intracellular biochemical and biophysical processes that may direct the cell's commitment to a particular lineage. In this study, we characterized the distinct biomechanical properties of hMSCs, including the average Young's modulus determined by atomic force microscopy (3.2 +/- 1.4 kPa for hMSC vs. 1.7 +/- 1.0 kPa for fully differentiated osteoblasts), and the average membrane tether length measured with laser optical tweezers (10.6 +/- 1.1 microm for stem cells, and 4.0 +/- 1.1 microm for osteoblasts). These differences in cell elasticity and membrane mechanics result primarily from differential actin cytoskeleton organization in these two cell types, whereas microtubules did not appear to affect the cellular mechanics. The membrane-cytoskeleton linker proteins may contribute to a stronger interaction of the plasma membrane with F-actins and shorter membrane tether length in osteoblasts than in stem cells. Actin depolymerization or ATP depletion caused a two- to threefold increase in the membrane tether length in osteoblasts, but had essentially no effect on the stem-cell membrane tethers. Actin remodeling in the course of a 10-day osteogenic differentiation of hMSC mediates the temporally correlated dynamical changes in cell elasticity and membrane mechanics. For example, after a 10-day culture in osteogenic medium, hMSC mechanical characteristics were comparable to those of mature bone cells. Based on quantitative characterization of the actin cytoskeleton remodeling during osteodifferentiation, we postulate that the actin cytoskeleton plays a pivotal role in determining the hMSC mechanical properties and modulation of cellular mechanics at the early stage of stem-cell osteodifferentiation.     

Coupling field theory with continuum mechanics: a simulation of domain formation in giant unilamellar vesicles

Domain formation is modeled on the surface of giant unilamellar vesicles using a Landau field theory model for phase coexistence coupled to elastic deformation mechanics (e.g., membrane curvature). Smooth particle applied mechanics, a form of smoothed particle continuum mechanics, is used to solve either the time-dependent Landau-Ginzburg or Cahn-Hilliard free-energy models for the composition dynamics. At the same time, the underlying elastic membrane is modeled using smooth particle applied mechanics, resulting in a unified computational scheme capable of treating the response of the composition fields to arbitrary deformations of the vesicle and vice versa. The results indicate that curvature coupling, along with the field theory model for composition free energy, gives domain formations that are correlated with surface defects on the vesicle. In the case that external deformations are included, the domain structures are seen to respond to such deformations. The present simulation capability provides a significant step forward toward the simulation of realistic cellular membrane processes.     

Actin cross-linkers and the shape of stereocilia

Stereocilia are actin-based cellular protrusions essential for hearing. We propose that they are shaped by the detachment dynamics of actin cross-linkers, in particular espin. We account for experimentally observed stereocilium shapes, treadmilling velocity to length relationship, espin 1 localization profile, and microvillus length to espin level relationship. If the cross-linkers are allowed to reattach, our model yields a dynamical phase transition toward unbounded growth. Considering the simplified case of a noninteracting, one-filament system, we calculate the length probability distribution in the growing phase and its stationary form in a continuum approximation of the finite-length phase. Numerical simulations of interacting filaments suggest an anomalous power-law divergence of the protrusion length at the growth transition, which could be a universal feature of cross-linked depolymerizing systems.     

Distinct membrane mechanical properties of human mesenchymal stem cells determined using laser optical tweezers

The therapeutic efficacy of mesenchymal stem cells (MSCs) in tissue engineering and regenerative medicine is determined by their unique biological, mechanical, and physicochemical characteristics, which are yet to be fully explored. Cell membrane mechanics, for example, has been shown to critically influence MSC differentiation. In this study, we used laser optical tweezers to measure the membrane mechanics of human MSCs and terminally differentiated fibroblasts by extracting tethers from the outer cell membrane. The average tether lengths were 10.6+/-1.1 microm (hMSC) and 3.0+/-0.5 microm (fibroblasts). The tether extraction force did not increase during tether formation, which suggests existence of a membrane reservoir intended to buffer membrane tension fluctuations. Cytoskeleton disruption resulted in a fourfold tether length increase in fibroblasts but had no effect in hMSCs, indicating weak association between the cell membrane and hMSC actin cytoskeleton. Cholesterol depletion, known to decrease lipid bilayer stiffness, caused an increase in the tether length both in fibroblasts and hMSCs, as does the treatment of cells with DMSO. We postulate that whereas fibroblasts use both the membrane rigidity and membrane-cytoskeleton association to regulate their membrane reservoir, hMSC cytoskeleton has only a minor impact on stem cell membrane mechanics.     

Curvature Generation and Pressure Profile Modulation in Membrane by Lysolipids: Insights from Coarse-Grained Simulations

Although many membrane additives are known to modulate the activities of membrane proteins via perturbing the properties of lipid membrane, the underlying mechanism is often not precisely understood. In this study, we investigate the impact of asymmetrically incorporating single-tailed lysophosphatidylcholine (LPC) into a membrane bilayer using coarse-grained molecular dynamics simulations. Using a simple computational protocol designed to approximately mimic a micropipette setting, we show that asymmetric incorporation of LPC can lead to significant curvature in a bilayer. Detailed analysis of geometrical and mechanical properties (pressure profile) of the resulting mound structure indicates that the degree of pressure profile perturbation is determined not by the local curvature per se but by the packing of lipid headgroups (i.e., area-per-lipid). The findings help provide a concrete basis for understanding the activation mechanism of mechanosensitive channels by asymmetric incorporation of LPC into membrane patches in patch-clamp experiments. The calculated local pressure profiles are valuable to the construction of realistic membrane models for the analysis of mechanosensation in a continuum mechanics framework.     

How can nematodes mate without spicules? Function of the male gonoduct glands in the roundworm Myolaimus

Males of roundworms (Nematoda) usually possess cuticular copulatory organs (spicules) that are inserted in the female's vulva to attach the male to the female and to widen the vulva against the inner body pressure for sperm transfer. Among free-living nematodes, the only exception of this rule is Myolaimus where the males lack spicules. Until now there exist no reports on how mating is achieved in Myolaimus. Here we show that sperm transfer in Myolaimus apparently involves at least six different secretions of the male gonoduct that are pumped into a sack-like cuticular protrusion of the female's vulva to form a spermatophore-like capsule. The role of gonoduct glands in male nematodes (even in the model organism Caenorhabditis elegans) is poorly understood. Here we present the first study explaining the role of different vas deferens gland products in nematodes and argue that Myolaimus males lost their spicules as a result of sperm competition.     

Flickering Analysis of Erythrocyte Mechanical Properties: Dependence on Oxygenation Level, Cell Shape, and Hydration Level

Erythrocytes (red blood cells) play an essential role in the respiratory functions of vertebrates, carrying oxygen from lungs to tissues and CO₂ from tissues to lungs. They are mechanically very soft, enabling circulation through small capillaries. The small thermally induced displacements of the membrane provide an important tool in the investigation of the mechanics of the cell membrane. However, despite numerous studies, uncertainties in the interpretation of the data, and in the values derived for the main parameters of cell mechanics, have rendered past conclusions from the fluctuation approach somewhat controversial. Here we revisit the experimental method and theoretical analysis of fluctuations, to adapt them to the case of cell contour fluctuations, which are readily observable experimentally. This enables direct measurements of membrane tension, of bending modulus, and of the viscosity of the cell cytoplasm. Of the various factors that influence the mechanical properties of the cell, we focus here on: 1), the level of oxygenation, as monitored by Raman spectrometry; 2), cell shape; and 3), the concentration of hemoglobin. The results show that, contrary to previous reports, there is no significant difference in cell tension and bending modulus between oxygenated and deoxygenated states, in line with the softness requirement for optimal circulatory flow in both states. On the other hand, tension and bending moduli of discocyte- and spherocyte-shaped cells differ markedly, in both the oxygenated and deoxygenated states. The tension in spherocytes is much higher, consistent with recent theoretical models that describe the transitions between red blood cell shapes as a function of membrane tension. Cell cytoplasmic viscosity is strongly influenced by the hydration state. The implications of these results to circulatory flow dynamics in physiological and pathological conditions are discussed.     

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.     

New Pudicinae (Trichostrongylina, Heligmosomoidea) coparasites of Proechimys longicaudatus (Caviomorpha) from Bolivia. I--Description of Pudica ginsburgi n. sp. and Heligmostrongylus chiarae n. sp

Two new Pudicinae (Trichostrongylina, Heligmosomoidea, Heligmonellidae) coparasites of Proechimys longicaudatus (Caviomorph rodent) from Bolivia are described: Pudica ginsburgi n. sp. is differentiated from all the other species of the genus by the great length of the spicules and of the vestibule and by a different ratio of the length of the spicules on the length of the body. Heligmostrongylus chiarae n. sp. has very developed rays 4, strongly divergent from rays 5 as occurs in three other species, H. almeidai (Durette-Desset & Tcheprakoff, 1969), a parasite of Trichomys apereoides (= Cercomys cunicularius) from Brazil, H. squamastrongylus (Travassos, 1937), a parasite of Proechimys oris from Brazil and H. proechimysi Durette-Desset, 1970, a parasite of Proechimys semispinosus from Columbia. The specimens from Bolivia are differentiated from the three species by the division of the dorsal ray (anterior half versus posterior half) and by a different ratio of the length of the spicules on the length of the body.     

Membrane bending is critical for the stability of voltage sensor segments in the membrane

The interaction between membrane proteins and the surrounding membrane is becoming increasingly appreciated for its role in regulating protein function, protein localization, and membrane morphology. In particular, recent studies have suggested that membrane deformation is needed to stably accommodate proteins harboring charged amino acids in their transmembrane (TM) region, as it is energetically prohibitive to bury charge in the hydrophobic core of the bilayer. Unfortunately, current computational methods are poorly equipped for describing such deformations, as atomistic simulations are often too short to observe large-scale membrane reorganization and most continuum approaches assume a flat membrane. Previously, we developed a method that overcomes these shortcomings by using elasticity theory to characterize equilibrium membrane distortions in the presence of a TM protein, while using traditional continuum electrostatic and nonpolar energy models to determine the energy of the protein in the membrane. Here, we linked the elastostatics, electrostatics, and nonpolar numeric solvers to permit the calculation of energies for nontrivial membrane deformations. We then coupled this procedure to a robust search algorithm that identifies optimal membrane shapes for a TM protein of arbitrary chemical composition. This advance now permits us to explore a host of biological phenomena that were beyond the scope of our original method. We show that the energy required to embed charged residues in the membrane can be highly nonadditive, and our model provides a simple mechanical explanation for this nonadditivity. Our results also predict that isolated voltage sensor segments do not insert into rigid membranes, but membrane bending dramatically stabilizes these proteins in the bilayer despite their high charge content. Additionally, we use the model to explore hydrophobic mismatch with regard to nonpolar peptides and mechanosensitive channels. Our method is in quantitative agreement with molecular dynamics simulations at a tiny fraction of the computational cost.