Stability of Asymmetric Lipid Bilayers Assessed by Molecular Dynamics Simulations

Designing nanoparticle translocation through membranes by computer simulations

Nanoparticle penetration into cells is an important process in drug/gene delivery. Here, we successfully design one type of novel nanoparticles with ligands decorating its surface by dynamic bonds and find that the nanoparticle can spontaneously penetrate through membranes by using dissipative particle dynamics simulations. Moreover, the physical parameters of both ligands (for example, ligand type and density) and nanoparticles (such as size and shape) have significant effects on penetration efficiency and translocation time. Especially for nanoparticles with anisotropic shapes or asymmetric surface decoration, the penetration efficiency may reach about 80%. We also provide insights into the interaction between nanoparticles and asymmetric membranes and find that the membrane asymmetry can even increase the penetration efficiency to above 90%. The present study suggests a potential way to translocate novel nanoparticles through membranes, which may provide new ideas for future experimental nanoparticle design and drug delivery.

Water Permeation Through DMPC Lipid Bilayers using Polarizable Charge Equilibration Force Fields

We investigate permeation energetics of water entering a model dimyristoylphosphatidylcholine (DMPC) bilayer via molecular dynamics simulations using polarizable Charge Equilibration (CHEQ) models. Potentials of mean force show 4.5-5.5 kcal/mol barriers for water permeation into bilayers. Barriers are highest when water coordination within the bilayer is prevented, and also when using force fields that accurately reproduce experimental alkane hydration free energies. The magnitude of the average water dipole moment decreases from 2.6 Debye (in bulk) to 1.88 Debye (in membrane interior). This variation correlates with the change in a water molecule's coordination number.

Protein shape change has a major effect on the gating energy of a mechanosensitive channel

Increasing experimental evidence has shown that membrane protein functionality depends on molecular composition of cell membranes. However, the origin of this dependence is not fully understood. It is reasonable to assume that specific lipid-protein interactions are important, yet more generic effects due to mechanical properties of lipid bilayers likely play a significant role too. Previously it has been demonstrated using models for elastic properties of membranes and lateral pressure profiles of lipid bilayers that the mechanical properties of a lipid bilayer can contribute as much as ∼10 k(B)T to the free energy difference associated with a change in protein conformational state. Here, we extend those previous approaches to a more realistic model for a large mechanosensitive channel (MscL). We use molecular dynamics together with the MARTINI model to simulate the open and closed states of MscL embedded in a DOPC bilayer. We introduce a procedure to calculate the mechanical energy change in the channel gating using a three-dimensional pressure distribution inside a membrane, computed from the molecular dynamics simulations. We decompose the mechanical energy to terms associated with area dilation and shape contribution. Our results highlight that the lateral pressure profile of a lipid bilayer together with the shape change in gating can induce a contribution of ∼30 k(B)T on the gating energy of MscL. This contribution arises largely from the interfacial tension between hydrophobic and hydrophilic regions in a lipid bilayer.

A lipocentric view of peptide-induced pores

Although lipid membranes serve as effective sealing barriers for the passage of most polar solutes, nonmediated leakage is not completely improbable. A high activation energy normally keeps unassisted bilayer permeation at a very low frequency, but lipids are able to self-organize as pores even in peptide-free and protein-free membranes. The probability of leakage phenomena increases under conditions such as phase coexistence, external stress or perturbation associated to binding of nonlipidic molecules. Here, we argue that pore formation can be viewed as an intrinsic property of lipid bilayers, with strong similarities in the structure and mechanism between pores formed with participation of peptides, lipidic pores induced by different types of stress, and spontaneous transient bilayer defects driven by thermal fluctuations. Within such a lipocentric framework, amphipathic peptides are best described as pore-inducing rather than pore-forming elements. Active peptides bound to membranes can be understood as a source of internal surface tension which facilitates pore formation by diminishing the high activation energy barrier. This first or immediate action of the peptide has some resemblance to catalysis. However, the presence of membrane-active peptides has the additional effect of displacing the equilibrium towards the pore-open state, which is then maintained over long times, and reducing the size of initial individual pores. Thus, pore-inducing peptides, regardless of their sequence and oligomeric organization, can be assigned a double role of increasing the probability of pore formation in membranes to high levels as well as stabilizing these pores after they appear.

Ca(2+) bridging of apposed phospholipid bilayers

In an effort to provide insight into the mechanism of Ca(2+)-induced fusion of lipid vesicles, molecular dynamics simulations in the isobaric-isothermal ensemble are used to investigate interactions of Ca(2+) with apposed lipid bilayers in close proximity. Simulations reveal the formation of a Ca(2+)-phospholipid "anhydrous complex" between apposed bilayers, whereas similar calculations performed with Na(+) display only complexation between neighboring lipids within the same bilayer. The binding of Ca(2+) to apposed phospholipids brings large regions of the bilayers into close contact (<4 Å), displacing water from phospholipid head groups in the process and creating regions of local dehydration. Dehydration of the apposed bilayers leads to ordering of the phospholipid tails, which is partially disrupted by the presence of Ca(2+)-phospholipid bridges.

Pores formed by Baxα5 relax to a smaller size and keep at equilibrium

Pores made by amphipathic cationic peptides (e.g., antimicrobials and fragments of pore-forming proteins) are typically studied by examining the kinetics of vesicle leakage after peptide addition or obtaining structural measurements in reconstituted peptide-lipid systems. In the first case, the pores have been considered transient phenomena that allow the relaxation of the peptide-membrane system. In the second, they correspond to equilibrium structures at minimum free energy. Here we reconcile both approaches by investigating the pore activity of the α5 fragment from the proapoptotic protein Bax (Baxα5) before and after equilibrium of peptide/vesicle complexes. Quenching assays on suspensions of large unilamellar vesicles suggest that in the presence of Baxα5, the vesicles maintain a leaky state for hours under equilibrium conditions. We proved and analyzed stable pores on single giant unilamellar vesicles (GUVs) in detail by monitoring the entrance of dyes added at different times after incubation with the peptide. When the GUVs came in contact with Baxα5, leakage started stochastically, was delayed for various periods of time, and in the majority of cases proceeded rapidly to completion. After hours in the presence of the peptide, the same individual GUVs that refilled completely at first instance maintained a porated state, which could be observed in subsequent leak-in events for serially added dyes. However, these long-term pores were smaller in size than the initial equilibration pores. Stable pores were also detected in GUVs made in the presence of Baxα5. The latter pores can be considered equilibrium states and may correspond to structures measured previously in bilayer stacks. Although pore formation may occur as a kinetic process, equilibrium pores may also be functionally relevant structures, especially in highly regulated systems such as the apoptotic mitochondrial pores induced by Bax.

Reactivity of a phospholipid monolayer model under periodic boundary conditions: a density functional theory study of the Schiff base formation between phosphatidylethanolamine and acetaldehyde

A mechanism for the formation of the Schiff base between an acetaldehyde and an amine-phospholipid monolayer model based on Dmol3/density functional theory calculations under periodic boundary conditions was constructed. This is the first time such a system has been modeled to examine its chemical reactivity at this computation level. Each unit cell contains two phospholipid molecules, one acetaldehyde molecule, and nine water molecules. One of the amine-phospholipid molecules in the cell possesses a neutral amino group that is used to model the nucleophilic attack on the carboxyl group of acetaldehyde, whereas the other has a charged amino group acting as a proton donor. The nine water molecules form a hydrogen bond network along the polar heads of the phospholipids that facilitates very fast proton conduction at the interface. Using periodic boundary conditions afforded proton transfer between different cells. The reaction takes place in two steps, namely, (1) formation of a carbinolamine and (2) its dehydration to the Schiff base. The carbinolamine is the primary reaction intermediate, and dehydration is the rate-determining step of the process, consistent with available experimental evidence for similar reactions. On the basis of the results, the cell membrane surface environment may boost phospholipid glycation via a neighboring catalyst effect.