Determining the local shear viscosity of a lipid bilayer system by reverse non-equilibrium molecular dynamics simulations

Self-Assembly of Lipid Molecules under Shear Flows: A Dissipative Particle Dynamics Simulation Study

The self-assembly of lipid molecules in aqueous solution under shear flows was investigated using the dissipative particle dynamics simulation method. Three cases were considered: zero shear flow, weak shear flow and strong shear flow. Various self-assembled structures, such as double layers, perforated double layers, hierarchical discs, micelles, and vesicles, were observed. The self-assembly behavior was investigated in equilibrium by constructing phase diagrams based on chain lengths. Results showed the remarkable influence of chain length, shear flow and solution concentration on the self-assembly process. Furthermore, the self-assembly behavior of lipid molecules was analyzed using the system energy, particle number and shape factor during the dynamic processes, where the self-assembly pathways were observed and analyzed for the typical structures. The results enhance our understanding of biomacromolecule self-assembly in a solution and hold the potential for applications in biomedicine.

Rheology of sliding leaflets in coarse-grained DSPC lipid bilayers

Amphiphilic lipid bilayers modify the friction properties of the surfaces on top of which they are deposited. In particular, the measured sliding friction coefficient can be significantly reduced compared with the native surface. We investigate in this work the friction properties of a numerical coarse-grained model of DSPC (1,2-distearoyl-sn-glycero-3-phosphocholine) lipid bilayer subject to longitudinal shear. The interleaflet friction coefficient is obtained from out-of-equilibrium pulling or from relaxation simulations. In particular, we gain access to the transient viscoelastic response of a sheared bilayer. The bilayer mechanical response is found to depend significantly on the membrane physical state, with evidence in favor of a linear response regime in the fluid but not in the gel region. The linear response validity domain is established, and the timescales appearing in the membrane response discussed.

Unexpected hydrophobicity on self-assembled monolayers terminated with two hydrophilic hydroxyl groups

Current major approaches to access surface hydrophobicity include directly introducing hydrophobic nonpolar groups/molecules onto the surface or elaborately fabricating surface roughness. Here, for the first time, molecular dynamics simulations show an unexpected hydrophobicity with a contact angle of 82° on a flexible self-assembled monolayer terminated only with two hydrophilic OH groups ((OH)₂-SAM). This hydrophobicity, verified by a water slip phenomenon characterizing the friction on the (OH)₂-SAM surface, is attributed to the formation of a hexagonal-ice-like H-bonding structure in the OH matrix of (OH)₂-SAM, which sharply reduces the hydrogen bonds between the surface and the water molecules above. The unique simple interface presented here offers a significant molecular-level platform for examining the bio-interfacial interactions ranging from biomolecule binding to cell adhesion.

Asymmetric Lipid Membranes under Shear Flows: A Dissipative Particle Dynamics Study

We investigate the phase behavior of the asymmetric lipid membranes under shear flows, using the dissipative particle dynamics simulation. Two cases, the weak and strong shear flows, are considered for the asymmetric lipid microstructures. Three typical asymmetric structures, the membranes, tubes, and vesicle, are included in the phase diagrams, where the effect of two different types of lipid chain length on the formation of asymmetric membranes is evaluated. The dynamic processes are demonstrated for the asymmetric membranes by calculating the average radius of gyration and shape factor. The result indicates that different shear flows will affect the shape of the second type of lipid molecules; the shape of the first type of lipid molecules is more stable than that of the second type of lipid molecules. The mechanical properties are investigated for the asymmetric membranes by analyzing the interface tension. The results reveal an absolute pressure at the junctions of different types of particles under the weak shear flow; the other positions are almost in a state of no pressure; there is almost no pressure inside the asymmetric lipid membrane structure under the strong shear flow. The findings will help us to understand the potential applications of asymmetric lipid microstructures in the biological and medical fields.

Regulating Lipid Composition Rationalizes Acyl Tail Saturation Homeostasis in Ectotherms

Cell membranes mainly consist of lipid bilayers with an actively regulated composition. The underlying processes are still poorly understood, in particular, how the hundreds of components are controlled. Cholesterol has been found to correlate with phospholipid saturation for reasons that remain unclear. To better understand the link between cell membrane regulation and chemical composition, we establish a computational framework based on chemical reaction networks, resulting in multiple semigrand canonical ensembles. By running computer simulations, we show that regulating the chemical potential of lipid species is sufficient to reproduce the experimentally observed increase in acyl tail saturation with added cholesterol. Our model proposes a different picture of lipid regulation in which components can be regulated passively instead of actively. In this picture, phospholipid acyl tail composition naturally adapts to added molecules such as cholesterol or proteins. A comparison between regulated membranes with commonly studied ternary model membranes shows stark differences: for instance, correlation lengths and viscosities observed are independent of lipid chemical affinity.

Shear-induced microstructures and dynamics processes of phospholipid cylinders in solutions

Shear-induced microstructures and their corresponding dynamic processes are investigated for phospholipid cylinders in aqueous solution by dissipative particle dynamic simulation. Various phospholipid cylinders with cross-sections, which are formed under shear-free flow, are selected to examine the effects of shear flow on their structures and dynamic processes. Shear flow induces the transition from cylinders into vesicles at weak rate and the transition into vesicle–lamella mixtures with increased shear rate and lamella structures at the strong shear rate. Then, the average radius of gyration and shape factors of the polymer chains in the dynamic processes are discussed in detail. Results show that shear flow causes the structure of the polymer chains to be elongated along the shear direction, and the configuration of the polymer chain can be rapidly transformed into an ellipsoid structure under strong shear.

Self-assembly of phospholipid molecules in solutions under shear flows: Microstructures and phase diagrams

Shear-induced microstructures and their phase diagrams were investigated for phospholipid molecules in aqueous solution by dissipative particle dynamic simulation. Self-assembled microstructures, including spherical and cylindrical micelles, spherical vesicles, lamellae, undulated lamellae, perforated lamellae, and continuous networks, were observed under various shear flows and phospholipid concentrations, where the spatial inhomogeneity and symmetry were analysed. A series of phase diagrams were constructed based on the chain lengths under various phospholipid concentrations. The phase distributions showed that the structures with spherical symmetry could be shear-induced to structures with cylindrical symmetry in the dilute solutions. In the semi-concentrated solutions, the lamellae were located in most spaces under zero shear flows, which could be shear-induced into undulated lamellae and then into cylindrical micelles. For the concentrated solutions, the strong shear flows oriented the directions of multilayer lamellae and phase transitions appeared between several cylindrical network structures. These observations on shear-induced microstructures and their distributions revealed a promising approach that could be used to design bio-microstructures based on phospholipid molecules under shear flows.

Adhesion dynamics of confined membranes

We report on the modeling of the dynamics of confined lipid membranes. We derive a thin film model in the lubrication limit which describes an inextensible liquid membrane with bending rigidity confined between two adhesive walls. The resulting equations share similarities with the Swift-Hohenberg model. However, inextensibility is enforced by a time-dependent nonlocal tension. Depending on the excess membrane area available in the system, three different dynamical regimes, denoted as A, B and C, are found from the numerical solution of the model. In regime A, membranes with small excess area form flat adhesion domains and freeze. Such freezing is interpreted by means of an effective model for curvature-driven domain wall motion. The nonlocal membrane tension tends to a negative value corresponding to the linear stability threshold of flat domain walls in the Swift-Hohenberg equation. In regime B, membranes with intermediate excess areas exhibit endless coarsening with coexistence of flat adhesion domains and wrinkle domains. The tension tends to the nonlinear stability threshold of flat domain walls in the Swift-Hohenberg equation. The fraction of the system covered by the wrinkle phase increases linearly with the excess area in regime B. In regime C, membranes with large excess area are completely covered by a frozen labyrinthine pattern of wrinkles. As the excess area is increased, the tension increases and the wavelength of the wrinkles decreases. For large membrane area, there is a crossover to a regime where the extrema of the wrinkles are in contact with the walls. In all regimes after an initial transient, robust localised structures form, leading to an exact conservation of the number of adhesion domains.

Thixotropy and shear thinning of lubricated contacts with confined membranes

We have modeled the nonlinear dynamics and the rheological behavior of a system under shear containing a membrane confined between two attractive walls. The presence of the membrane induces additional tangential forces on the walls that always increase the global friction. At low shear rates, the membrane exhibits chaotic dynamics with slow coarsening leading to thixotropy, i.e. to a slow decrease of the membrane-induced tangential forces on the walls. At intermediate shear rates, the membrane profile presents stationary periodic patterns. At higher shear rates, membrane dynamics are governed by a nonlinear evolution equation which is similar to the Kuramoto-Sivashinski equation, but with a sixth-order stabilizing term. The membrane experiences chaotic dynamics without coarsening. As a consequence of the nonlinear dynamics of the membrane at intermediate and large shear rates, the system exhibits shear thinning.

Lipid and Peptide Diffusion in Bilayers: The Saffman-Delbrück Model and Periodic Boundary Conditions

The periodic Saffman-Delbrück (PSD) model, an extension of the Saffman-Delbrück model developed to describe the effects of periodic boundary conditions on the diffusion constants of lipids and proteins obtained from simulation, is tested using the coarse-grained Martini and all-atom CHARMM36 (C36) force fields. Simulations of pure Martini dipalmitoylphosphatidylcholine (DPPC) bilayers and those with one embedded gramicidin A (gA) dimer or one gA monomer with sizes ranging from 512 to 2048 lipids support the PSD model. Underestimates of D^∞ (the value of the diffusion constant for an infinite system) from the 512-lipid system are 35% for DPPC, 45% for the gA monomer, and 70% for the gA dimer. Simulations of all-atom DPPC and dioleoylphosphatidylcholine (DOPC) bilayers yield diffusion constants not far from experiment. However, the PSD model predicts that diffusion constants at the sizes of the simulation should underestimate experiment by approximately a factor of 3 for DPPC and 2 for DOPC. This likely implies a deficiency in the C36 force field. A Bayesian method for extrapolating diffusion constants of lipids and proteins in membranes obtained from simulation to infinite system size is provided.

Strong influence of periodic boundary conditions on lateral diffusion in lipid bilayer membranes

The Saffman-Delbrück hydrodynamic model for lipid-bilayer membranes is modified to account for the periodic boundary conditions commonly imposed in molecular simulations. Predicted lateral diffusion coefficients for membrane-embedded solid bodies are sensitive to box shape and converge slowly to the limit of infinite box size, raising serious doubts for the prospects of using detailed simulations to accurately predict membrane-protein diffusivities and related transport properties. Estimates for the relative error associated with periodic boundary artifacts are 50% and higher for fully atomistic models in currently feasible simulation boxes. MARTINI simulations of LacY membrane protein diffusion and LacY dimer diffusion in DPPC membranes and lipid diffusion in pure DPPC bilayers support the underlying hydrodynamic model.

Mixed Mechanism of Lubrication by Lipid Bilayer Stacks

Although the key role of lipid bilayer stacks in biological lubrication is generally accepted, the mechanisms underlying their extreme efficiency remain elusive. In this article, we report molecular dynamics simulations of lipid bilayer stacks undergoing load and shear. When the hydration level is reduced, the velocity accommodation mechanism changes from viscous shear in hydration water to interlayer sliding in the bilayers. This enables stacks of hydrated lipid bilayers to act as efficient boundary lubricants for various hydration conditions, structures, and mechanical loads. We also propose an estimation for the friction coefficient; thanks to the strong hydration forces between lipid bilayers, the high local viscosity is not in contradiction with low friction coefficients.

Frozen states and order-disorder transition in the dynamics of confined membranes

The adhesion dynamics of a membrane confined between two permeable walls is studied using a two-dimensional hydrodynamic model. The membrane morphology decomposes into adhesion patches on the upper and the lower walls and obeys a nonlinear evolution equation that resembles that of phase-separation dynamics, which is known to lead to coarsening, i.e., to the endless growth of the adhesion patches. However, due to the membrane bending rigidity, the system evolves toward a frozen state without coarsening. This frozen state exhibits an order-disorder transition when increasing the permeability of the walls.

Interleaflet sliding in lipidic bilayers under shear flow: comparison of the gel and fluid phases using reversed non-equilibrium molecular dynamics simulations

The friction between two rubbing surfaces lubricated by water can be diminished if they are coated with phospholipidic bilayers or brushes of polyelectrolytes. In the case of a coating by lipid membranes, the friction is lower when the lipids are in the gel phase rather than in the liquid phase. We investigated the response of fluid or gel bilayers to a mechanical load or under shear using non-equilibrium molecular dynamics simulations (NEMD) to understand whether this difference could come from intermonolayer sliding. The system is composed of a single fully hydrated bilayer of coarse grained phospholipids under a parallel shear with vorticity parallel to the bilayer. In both the liquid and the gel phases, an intermonolayer slip was measured in the velocity profile. In the liquid phase this slip is proportional to the shear stress. In the tilted gel phase of our model the stress is not systematically linear and relaxes differently when the shear is in the direction of the tilt or perpendicular to it. The impact of surface tension (or load) on the friction is different for the liquid and gel phases, but grossly the slip remains of the same order of magnitude.

Fluctuating hydrodynamics simulations of coarse-grained lipid membranes under steady-state conditions and in shear flow

The stochastic Eulerian-Lagrangian method (SELM) is used to simulate coarse-grained lipid membrane models under steady-state conditions and in shear flow. SELM is an immersed boundary method which combines the efficiency of particle-based simulations with the realistic solvent dynamics provided by fluctuating hydrodynamics. Membrane simulations in SELM are shown to give structural properties in accordance with equilibrium statistical mechanics and dynamic properties in agreement with previous simulations of highly detailed membrane models in explicit solvent. Simulations of sheared membranes are used to calculate surface shear viscosities and inter-monolayer friction coefficients. The membrane models are shown to be shear thinning under a wide range of applied shear rates.

Molecular dynamics simulations of depth distribution of spin-labeled phospholipids within lipid bilayer

Spin-labeled lipids are commonly used as fluorescence quenchers in studies of membrane penetration of dye-labeled proteins and peptides using depth-dependent quenching. Accurate calculations of depth of the fluorophore rely on the use of several spin labels placed in the membrane at various positions. The depth of the quenchers (spin probes) has to be determined independently; however, experimental determination of transverse distributions of spin probe depths is difficult. In this Article, we use molecular dynamics (MD) simulations to study the membrane behavior and depth distributions of spin-labeled phospholipids in a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer. To probe different depths within the bilayer, a series containing five Doxyl-labeled lipids (n-Doxyl PC) has been studied, in which a spin moiety was covalently attached to nth carbon atoms (where n = 5, 7, 10, 12, and 14) of the sn-2 stearoyl chain of the host phospholipid. Our results demonstrate that the chain-attached spin labels are broadly distributed across the model membrane and their environment is characterized by a high degree of mobility and structural heterogeneity. Despite the high thermal disorder, the depth distributions of the Doxyl labels were found to correlate well with their attachment positions, indicating that the distribution of the spin label within the model membrane is dictated by the depth of the nth lipid carbon atom and not by intrinsic properties of the label. In contrast, a much broader and heterogeneous distribution was observed for a headgroup-attached Tempo spin label of Tempo-PC lipids. MD simulations reveal that, due to the hydrophobic nature, a Tempo moiety favors partitioning from the headgroup region deeper into the membrane. Depending on the concentration of Tempo-PC lipids, the probable depth of the Tempo moiety could span a range from 14.4 to 18.2 Å from the membrane center. Comparison of the MD-estimated immersion depths of Tempo and n-Doxyl labels with their suggested experimental depth positions allows us to review critically the possible sources of error in depth-dependent fluorescence quenching studies.

Membrane-bound myo1c powers asymmetric motility of actin filaments

Class I myosins are molecular motors that link cellular membranes to the actin cytoskeleton and play roles in membrane tension generation, membrane dynamics, and mechanosignal transduction. The widely expressed myosin-Ic (myo1c) isoform binds tightly to phosphatidylinositol 4,5-bisphosphate [PtdIns(4,5)P(2)] via a pleckstrin homology domain located in the myo1c tail, which is important for its proper cellular localization. In this study, we found that myo1c can power actin motility on fluid membranes composed of physiological concentrations of PtdIns(4,5)P(2) and that this motility is inhibited by high concentrations of anionic phospholipids. Strikingly, this motility occurs along curved paths in a counterclockwise direction (i.e., the actin filaments turn in leftward circles). A biotinylated myo1c construct containing only the motor domain and the lever arm anchored via streptavidin on a membrane containing biotinylated lipid can also generate asymmetric motility, suggesting that the tail domain is not required for the counterclockwise turning. We found that the ability to produce counterclockwise motility is not a universal characteristic of myosin-I motors, as membrane-bound myosin-Ia (myo1a) and myosin-Ib (myo1b) are able to power actin gliding, but the actin gliding has no substantial turning bias. This work reveals a possible mechanism for establishing asymmetry in relationship to the plasma membrane.