Lipid Clustering Correlates with Membrane Curvature as Revealed by Molecular Simulations of Complex Lipid Bilayers

Cell membranes are complex multicomponent systems, which are highly heterogeneous in the lipid distribution and composition. To date, most molecular simulations have focussed on relatively simple lipid compositions, helping to inform our understanding of in vitro experimental studies. Here we describe on simulations of complex asymmetric plasma membrane model, which contains seven different lipids species including the glycolipid GM3 in the outer leaflet and the anionic lipid, phosphatidylinositol 4,5-bisphophate (PIP₂), in the inner leaflet. Plasma membrane models consisting of 1500 lipids and resembling the in vivo composition were constructed and simulations were run for 5 µs. In these simulations the most striking feature was the formation of nano-clusters of GM3 within the outer leaflet. In simulations of protein interactions within a plasma membrane model, GM3, PIP₂, and cholesterol all formed favorable interactions with the model α-helical protein. A larger scale simulation of a model plasma membrane containing 6000 lipid molecules revealed correlations between curvature of the bilayer surface and clustering of lipid molecules. In particular, the concave (when viewed from the extracellular side) regions of the bilayer surface were locally enriched in GM3. In summary, these simulations explore the nanoscale dynamics of model bilayers which mimic the in vivo lipid composition of mammalian plasma membranes, revealing emergent nanoscale membrane organization which may be coupled both to fluctuations in local membrane geometry and to interactions with proteins.     

Annular Anionic Lipids Stabilize the Integrin αIIbβ3 Transmembrane Complex

Cationic membrane-proximal amino acids determine the topology of membrane proteins by interacting with anionic lipids that are restricted to the intracellular membrane leaflet. This mechanism implies that anionic lipids interfere with electrostatic interactions of membrane proteins. The integrin αIIbβ3 transmembrane (TM) complex is stabilized by a membrane-proximal αIIb(Arg⁹⁹⁵)-β3(Asp⁷²³) interaction; here, we examine the influence of anionic lipids on this complex. Anionic lipids compete for αIIb(Arg⁹⁹⁵) contacts with β3(Asp⁷²³) but paradoxically do not diminish the contribution of αIIb(Arg⁹⁹⁵)-β3(Asp⁷²³) to TM complex stability. Overall, anionic lipids in annular positions stabilize the αIIbβ3 TM complex by up to 0.50 ± 0.02 kcal/mol relative to zwitterionic lipids in a headgroup structure-dependent manner. Comparatively, integrin receptor activation requires TM complex destabilization of 1.5 ± 0.2 kcal/mol, revealing a sizeable influence of lipid composition on TM complex stability. We implicate changes in lipid headgroup accessibility to small molecules (physical membrane characteristics) and specific but dynamic protein-lipid contacts in this TM helix-helix stabilization. Thus, anionic lipids in ubiquitous annular positions can benefit the stability of membrane proteins while leaving membrane-proximal electrostatic interactions intact.     

Implicit solvent model estimates of the stability of model structures of the alamethicin channel

Alamethicin is a hydrophobic helical peptide of 20 residues, which oligomerizes to form ion-conducting channels in membranes. The behavior of an intact alamethicin channel in POPC bilayers was recently studied, using 2 ns molecular dynamics (MD) simulations of a model hexameric channel. These simulations produced numerous conformations of the channel. In the present study, we used 11 of these channel conformations and carried out continuum-solvent model calculations, similar to those used for the monomers in our previous studies, to investigate the energetics of the channel inside the lipid bilayer. Our results suggest that, out of the 11 channel conformations produced by the MD simulations, only four are stable inside the lipid bilayer, with water-to-membrane free energies of transfer ranging from ~–6 to ~–10 kcal/mol. Analysis of the results suggests two causes for the apparent instability of the remainder of the structures inside the lipid bilayer, both resulting from the desolvation of channel polar groups (i.e. their transfer from the aqueous phase into the bilayer). The first is specific, uncompensated backbone hydrogen bonds, which exist in the region of the channel exposed to the hydrocarbon of the lipid bilayer. The second is exposure of intra-pore water molecules to the surrounding lipid. Thus, the association of these structures with the membrane involves a large electrostatic desolvation free-energy penalty. The apparent conflict between continuum-solvent and MD calculations, and its significance for the interpretation of membrane proteins simulations, are discussed.     

Cryo-EM of the ATP11C flippase reconstituted in Nanodiscs shows a distended phospholipid bilayer inner membrane around transmembrane helix 2

ATP11C is a member of the P4-ATPase flippase family that mediates translocation of phosphatidylserine (PtdSer) across the lipid bilayer. In order to characterize the structure and function of ATP11C in a model natural lipid environment, we revisited and optimized a quick procedure for reconstituting ATP11C into Nanodiscs using methyl-β-cyclodextrin as a reagent for the detergent removal. ATP11C was efficiently reconstituted with the endogenous lipid, or the mixture of endogenous lipid and synthetic dioleoylphosphatidylcholine (DOPC)/dioleoylphosphatidylserine (DOPS), all of which retained the ATPase activity. We obtained 3.4 Å and 3.9 Å structures using single-particle cryo-electron microscopy (cryo-EM) of AlF- and BeF-stabilized ATP11C transport intermediates, respectively, in a bilayer containing DOPS. We show that the latter exhibited a distended inner membrane around ATP11C transmembrane helix 2, possibly reflecting the perturbation needed for phospholipid release to the lipid bilayer. Our structures of ATP11C in the lipid membrane indicate that the membrane boundary varies upon conformational changes of the enzyme and is no longer flat around the protein, a change that likely contributes to phospholipid translocation across the membrane leaflets.     

Finding a Needle in a Haystack: The Role of Electrostatics in Target Lipid Recognition by PH Domains

Interactions between protein domains and lipid molecules play key roles in controlling cell membrane signalling and trafficking. The pleckstrin homology (PH) domain is one of the most widespread, binding specifically to phosphatidylinositol phosphates (PIPs) in cell membranes. PH domains must locate specific PIPs in the presence of a background of approximately 20% anionic lipids within the cytoplasmic leaflet of the plasma membrane. We investigate the mechanism of such recognition via a multiscale procedure combining Brownian dynamics (BD) and molecular dynamics (MD) simulations of the GRP1 PH domain interacting with phosphatidylinositol (3,4,5)-trisphosphate (PI(3,4,5)P(3)). The interaction of GRP1-PH with PI(3,4,5)P(3) in a zwitterionic bilayer is compared with the interaction in bilayers containing different levels of anionic 'decoy' lipids. BD simulations reveal both translational and orientational electrostatic steering of the PH domain towards the PI(3,4,5)P(3)-containing anionic bilayer surface. There is a payoff between non-PIP anionic lipids attracting the PH domain to the bilayer surface in a favourable orientation and their role as 'decoys', disrupting the interaction of GRP1-PH with the PI(3,4,5)P(3) molecule. Significantly, approximately 20% anionic lipid in the cytoplasmic leaflet of the bilayer is nearly optimal to both enhance orientational steering and to localise GRP1-PH proximal to the surface of the membrane without sacrificing its ability to locate PI(3,4,5)P(3) within the bilayer plane. Subsequent MD simulations reveal binding to PI(3,4,5)P(3), forming protein-phosphate contacts comparable to those in X-ray structures. These studies demonstrate a computational framework which addresses lipid recognition within a cell membrane environment, offering a link between structural and cell biological characterisation.     

Characterization of a Membrane-active Peptide from the Bordetella pertussis CyaA Toxin

Bordetella pertussis, the pathogenic bacteria responsible for whooping cough, secretes several virulence factors, among which is the adenylate cyclase toxin (CyaA) that plays a crucial role in the early stages of human respiratory tract colonization. CyaA invades target cells by translocating its catalytic domain directly across the plasma membrane and overproduces cAMP, leading to cell death. The molecular process leading to the translocation of the catalytic domain remains largely unknown. We have previously shown that the catalytic domain per se, AC384, encompassing residues 1–384 of CyaA, did not interact with lipid bilayer, whereas a longer polypeptide, AC489, spanning residues 1–489, binds to membranes and permeabilizes vesicles. Moreover, deletion of residues 375–485 within CyaA abrogated the translocation of the catalytic domain into target cells. Here, we further identified within this region a peptidic segment that exhibits membrane interaction properties. A synthetic peptide, P454, corresponding to this sequence (residues 454–485 of CyaA) was characterized by various biophysical approaches. We found that P454 (i) binds to membranes containing anionic lipids, (ii) adopts an α-helical structure oriented in plane with respect to the lipid bilayer, and (iii) permeabilizes vesicles. We propose that the region encompassing the helix 454–485 of CyaA may insert into target cell membrane and induce a local destabilization of the lipid bilayer, thus favoring the translocation of the catalytic domain across the plasma membrane.     

Molecular Dynamics Simulation Analysis of Membrane Defects and Pore Propensity of Hemifusion Diaphragms

Membrane fusion often exhibits slow dynamics in electrophysiological experiments, involving prespike foot and fusion pore-flickering, but the structural basis of such phenomena remains unclear. Hemifusion intermediates have been implicated in the early phase of membrane fusion. To elucidate the dynamics of formation of membrane defects and pores within the hemifusion diaphragm (HD), atomistic and coarse-grained models of hemifusion intermediates were constructed using dipalmitoylphosphatidylcholine or dioleoylphosphatidylcholine membranes. The work necessary to displace a lipid molecule to the hydrophobic core of the bilayer was measured. For a lipid within the HD with radius of 4 nm, the work was ∼80 kJ/mol, similar to that in a planar bilayer. The work was much less (∼40 kJ/mol) when the HD was surrounded by a steep stalk, i.e., stalk wings forming a large angle at the junction of three bilayers. In the latter case, the lipid displacement engendered formation of a pore contacting the HD rim. The work was similarly small (40 kJ/mol) for a small HD of 1.5 nm radius, where a pore formed and grew rapidly, quickly generating a toroidal structure (<40 ns). Combining the steep stalk and the small HD decreased the work further, although quantitative analysis was difficult because the latter system was not in a stable equilibrium state. Results suggest that fine tuning of fusion dynamics requires strict control of the HD size and the angle between the expanded stalk and HD. In additional free simulations, the steep stalk facilitated widening of a preformed pore contacting the HD rim.     

Surface Electrostatics of Lipid Bilayers by EPR of a pH-Sensitive Spin-Labeled Lipid

Many biophysical processes such as insertion of proteins into membranes and membrane fusion are governed by bilayer electrostatic potential. At the time of this writing, the arsenal of biophysical methods for such measurements is limited to a few techniques. Here we describe a, to our knowledge, new spin-probe electron paramagnetic resonance (EPR) approach for assessing the electrostatic surface potential of lipid bilayers that is based on a recently synthesized EPR probe (IMTSL-PTE) containing a reversibly ionizable nitroxide tag attached to the lipids’ polar headgroup. EPR spectra of the probe directly report on its ionization state and, therefore, on electrostatic potential through changes in nitroxide magnetic parameters and the degree of rotational averaging. Further, the lipid nature of the probe provides its full integration into lipid bilayers. Tethering the nitroxide moiety directly to the lipid polar headgroup defines the location of the measured potential with respect to the lipid bilayer interface. Electrostatic surface potentials measured by EPR of IMTSL-PTE show a remarkable (within ±2%) agreement with the Gouy-Chapman theory for anionic DMPG bilayers in fluid (48°C) phase at low electrolyte concentration (50 mM) and in gel (17°C) phase at 150-mM electrolyte concentration. This agreement begins to diminish for DMPG vesicles in gel phase (17°C) upon varying electrolyte concentration and fluid phase bilayers formed from DMPG/DMPC and POPG/POPC mixtures. Possible reasons for such deviations, as well as the proper choice of an electrostatically neutral reference interface, have been discussed. Described EPR method is expected to be fully applicable to more-complex models of cellular membranes.     

Extracting Curvature Preferences of Lipids Assembled in Flat Bilayers Shows Possible Kinetic Windows for Genesis of Bilayer Asymmetry and Domain Formation in Biological Membranes

Studies on the assembly of pure lipid components allow mechanistic insights toward understanding the structural and functional aspects of biological membranes. Molecular dynamic (MD) simulations on membrane systems provide molecular details on membrane dynamics that are difficult to obtain experimentally. A large number of MD studies have remained somewhat disconnected from a key concept of amphipathic assembly resulting in membrane structures—shape parameters of lipid molecules in those structures in aqueous environments. This is because most of the MD studies have been done on flat lipid membranes. With the above in view, we analyzed MD simulations of 26 pure lipid patches as a function of (1) lipid type(s) and (2) time of MD simulations along with 35–40 ns trajectories of five pure lipids. We report, for the first time, extraction of curvature preferences of lipids from MD simulations done on flat bilayers. Our results may lead to mechanistic insights into the possible origins of bilayer asymmetries and domain formation in biological membranes.     

Interactions of the Anticancer Drug Tamoxifen with Lipid Membranes

Interactions of the hydrophobic anticancer drug tamoxifen (TAM) with lipid model membranes were studied using calcein-encapsulated vesicle leakage, attenuated total reflection Fourier transform infrared (FTIR) spectroscopy, small-angle neutron scattering (SANS), atomic force microscopy (AFM) based force spectroscopy, and all-atom molecular dynamics (MD) simulations. The addition of TAM enhances membrane permeability, inducing calcein to translocate from the interior to the exterior of lipid vesicles. A large decrease in the FTIR absorption band’s magnitude was observed in the hydrocarbon chain region, suggesting suppressed bond vibrational dynamics. Bilayer thickening was determined from SANS data. Force spectroscopy measurements indicate that the lipid bilayer area compressibility modulus K_A is increased by a large amount after the incorporation of TAM. MD simulations show that TAM decreases the lipid area and increases chain order parameters. Moreover, orientational and positional analyses show that TAM exhibits a highly dynamic conformation within the lipid bilayer. Our detailed experimental and computational studies of TAM interacting with model lipid membranes shed new light on membrane modulation by TAM.     

Insertion of short amino-functionalized single-walled carbon nanotubes into phospholipid bilayer occurs by passive diffusion

Carbon nanotubes have been proposed to be efficient nanovectors able to deliver genetic or therapeutic cargo into living cells. However, a direct evidence of the molecular mechanism of their translocation across cell membranes is still needed. Here, we report on an extensive computational study of short (5 nm length) pristine and functionalized single-walled carbon nanotubes uptake by phospholipid bilayer models using all-atom molecular dynamics simulations. Our data support the hypothesis of a direct translocation of the nanotubes through the phospholipid membrane. We find that insertion of neat nanotubes within the bilayer is a "nanoneedle" like process, which can often be divided in three consecutive steps: landing and floating, penetration of the lipid headgroup area and finally sliding into the membrane core. The presence of functional groups at moderate concentrations does not modify the overall scheme of diffusion mechanism, provided that their deprotonated state favors translocation through the lipid bilayer.     

E.?coli Outer Membrane and Interactions with OmpLA

The outer membrane of Gram-negative bacteria is a unique asymmetric lipid bilayer composed of phospholipids (PLs) in the inner leaflet and lipopolysaccharides (LPSs) in the outer leaflet. Its function as a selective barrier is crucial for the survival of bacteria in many distinct environments, and it also renders Gram-negative bacteria more resistant to antibiotics than their Gram-positive counterparts. Here, we report the structural properties of a model of the Escherichia coli outer membrane and its interaction with outer membrane phospholipase A (OmpLA) utilizing molecular dynamics simulations. Our results reveal that given the lipid composition used here, the hydrophobic thickness of the outer membrane is ∼3 Å thinner than the corresponding PL bilayer, mainly because of the thinner LPS leaflet. Further thinning in the vicinity of OmpLA is observed due to hydrophobic matching. The particular shape of the OmpLA barrel induces various interactions between LPS and PL leaflets, resulting in asymmetric thinning around the protein. The interaction between OmpLA extracellular loops and LPS (headgroups and core oligosaccharides) stabilizes the loop conformation with reduced dynamics, which leads to secondary structure variation and loop displacement compared to that in a DLPC bilayer. In addition, we demonstrate that the LPS/PL ratios in asymmetric bilayers can be reliably estimated by the per-lipid surface area of each lipid type, and there is no statistical difference in the overall membrane structure for the outer membranes with one more or less LPS in the outer leaflet, although individual lipid properties vary slightly.     

Passive transport of Ca2+ ions through lipid bilayers imaged by widefield second harmonic microscopy

In biology, release of Ca²⁺ ions in the cytosol is essential to trigger or control many cell functions. Calcium signaling acutely depends on lipid membrane permeability to Ca²⁺. For proper understanding of membrane permeability to Ca²⁺, both membrane hydration and the structure of the hydrophobic core must be taken into account. Here, we vary the hydrophobic core of bilayer membranes and observe different types of behavior in high-throughput wide-field second harmonic imaging. Ca²⁺ translocation is observed through mono-unsaturated (DOPC:DOPA) membranes, reduced upon the addition of cholesterol, and completely inhibited for branched (DPhPC:DPhPA) and poly-unsaturated (SLPC:SLPA) lipid membranes. We propose, using molecular dynamics simulations, that ion transport occurs through ion-induced transient pores, which requires nonequilibrium membrane restructuring. This results in different rates at different locations and suggests that the hydrophobic structure of lipids plays a much more sophisticated regulating role than previously thought.     

Membrane perturbation by the antimicrobial peptide PMAP-23: A fluorescence and molecular dynamics study

Several bioactive peptides exert their biological function by interacting with cellular membranes. Structural data on their location inside lipid bilayers are thus essential for a detailed understanding of their mechanism of action. We propose here a combined approach in which fluorescence spectroscopy and molecular dynamics (MD) simulations were applied to investigate the mechanism of membrane perturbation by the antimicrobial peptide PMAP-23. Fluorescence spectra, depth-dependent quenching experiments, and peptide-translocation assays were employed to determine the location of the peptide inside the membrane. MD simulations were performed starting from a random mixture of water, lipids and peptide, and following the spontaneous self-assembly of the bilayer. Both experimental and theoretical data indicated a peptide location just below the polar headgroups of the membrane, with an orientation essentially parallel to the bilayer plane. These findings, together with experimental results on peptide-induced leakage from large and giant vesicles, lipid flip-flop and peptide exchange between vesicles, support a mechanism of action consistent with the “carpet” model. Furthermore, the atomic detail provided by the simulations suggested the occurrence of an additional, more specific and novel mechanism of bilayer destabilization by PMAP-23, involving the unusual insertion of charged side chains into the hydrophobic core of the membrane.     

Partitioning of 2,6-Bis(1H-Benzimidazol-2-yl)pyridine fluorophore into a phospholipid bilayer: complementary use of fluorescence quenching studies and molecular dynamics simulations

Successful use of fluorescence sensing in elucidating the biophysical properties of lipid membranes requires knowledge of the distribution and location of an emitting molecule in the bilayer. We report here that 2,6-bis(1H-benzimidazol-2-yl)pyridine (BBP), which is almost non-fluorescent in aqueous solutions, reveals a strong emission enhancement in a hydrophobic environment of a phospholipid bilayer, making it interesting for fluorescence probing of water content in a lipid membrane. Comparing the fluorescence behavior of BBP in a wide variety of solvents with those in phospholipid vesicles, we suggest that the hydrogen bonding interactions between a BBP fluorophore and water molecules play a crucial role in the observed “light switch effect”. Therefore, the loss of water-induced fluorescence quenching inside a membrane are thought to be due to deep penetration of BBP into the hydrophobic, water-free region of a bilayer. Characterized by strong quenching by transition metal ions in solution, BBP also demonstrated significant shielding from the action of the quencher in the presence of phospholipid vesicles. We used the increase in fluorescence intensity, measured upon titration of probe molecules with lipid vesicles, to estimate the partition constant and the Gibbs free energy (ΔG) of transfer of BBP from aqueous buffer into a membrane. Partitioning BBP revealed strongly favorable ΔG, which depends only slightly on the lipid composition of a bilayer, varying in a range from − 6.5 to − 7.0 kcal/mol. To elucidate the binding interactions of the probe with a membrane on the molecular level, a distribution and favorable location of BBP in a POPC bilayer were modeled via atomistic molecular dynamics (MD) simulations using two different approaches: (i) free, diffusion-driven partitioning of the probe molecules into a bilayer and (ii) constrained umbrella sampling of a penetration profile of the dye molecule across a bilayer. Both of these MD approaches agreed with regard to the preferred location of a BBP fluorophore within the interfacial region of a bilayer, located between the hydrocarbon acyl tails and the initial portion of the lipid headgroups. MD simulations also revealed restricted permeability of water molecules into this region of a POPC bilayer, determining the strong fluorescence enhancement observed experimentally for the membrane-partitioned form of BBP.     

Multiscale Simulations Suggest a Mechanism for the Association of the Dok7 PH Domain with PIP-Containing Membranes

Dok7 is a peripheral membrane protein that is associated with the MuSK receptor tyrosine kinase. Formation of the Dok7/MuSK/membrane complex is required for the activation of MuSK. This is a key step in the complex exchange of signals between neuron and muscle, which lead to neuromuscular junction formation, dysfunction of which is associated with congenital myasthenic syndromes. The Dok7 structure consists of a Pleckstrin Homology (PH) domain and a Phosphotyrosine Binding (PTB) domain. The mechanism of the Dok7 association with the membrane remains largely unknown. Using multi-scale molecular dynamics simulations we have explored the formation of the Dok7 PH/membrane complex. Our simulations indicate that the PH domain of Dok7 associates with membranes containing phosphatidylinositol phosphates (PIPs) via interactions of the β1/β2, β3/β4, and β5/β6 loops, which together form a positively charged surface on the PH domain and interact with the negatively charged headgroups of PIP molecules. The initial encounter of the Dok7 PH domain is followed by formation of additional interactions with the lipid bilayer, and especially with PIP molecules, which stabilizes the Dok7 PH/membrane complex. We have quantified the binding of the PH domain to the model bilayers by calculating a density landscape for protein/membrane interactions. Detailed analysis of the PH/PIP interactions reveal both a canonical and an atypical site to be occupied by the anionic lipid. PH domain binding leads to local clustering of PIP molecules in the bilayer. Association of the Dok7 PH domain with PIP lipids is therefore seen as a key step in localization of Dok7 to the membrane and formation of a complex with MuSK.     

Nanopore-facilitated, voltage-driven phosphatidylserine translocation in lipid bilayers--in cells and in silico

Nanosecond, megavolt-per-meter pulses--higher power but lower total energy than the electroporative pulses used to introduce normally excluded material into biological cells--produce large intracellular electric fields without destructively charging the plasma membrane. Nanoelectropulse perturbation of mammalian cells causes translocation of phosphatidylserine (PS) to the outer face of the cell, intracellular calcium release, and in some cell types a subsequent progression to apoptosis. Experimental observations and molecular dynamics (MD) simulations of membranes in pulsed electric fields presented here support the hypothesis that nanoelectropulse-induced PS externalization is driven by the electric potential that appears across the lipid bilayer during a pulse and is facilitated by the poration of the membrane that occurs even during pulses as brief as 3 ns. MD simulations of phospholipid bilayers in supraphysiological electric fields show a tight association between PS externalization and membrane pore formation on a nanosecond time scale that is consistent with experimental evidence for electropermeabilization and anode-directed PS translocation after nanosecond electric pulse exposure, suggesting a molecular mechanism for nanoelectroporation and nanosecond PS externalization: electrophoretic migration of the negatively charged PS head group along the surface of nanometer-diameter electropores initiated by field-driven alignment of water dipoles at the membrane interface.     

The Rate Constants of Valinomycin-Mediated Ion Transport through Thin Lipid Membranes

Electrical relaxation experiments have been performed with phosphatidylinositol bilayer membranes in the presence of the ion carrier valinomycin. After a sudden change of the voltage a relaxation of the membrane current with a time constant of about 20 μsec is observed. Together with previous stationary conductance data, the relaxation amplitude and the relaxation time are used to evaluate the rate constants of valinomycin-mediated potassium transport across the lipid membrane. It is found that the rate constants of translocation of the free carrier S and the carrier-ion complex MS⁺ are nearly equal (2·10⁴ sec^-1) and are of the same order as the dissociation rate constant of MS⁺ in the membrane-solution interface (5·10⁴ sec^-1). The equilibrium constant of the heterogeneous association reaction M⁺ (solution) + S (membrane) → MS⁺ (membrane) is found to be ~ 1 M^-1, about 10⁶ times smaller than the association constant in ethanolic solution.     

A Coarse Grained Model for a Lipid Membrane with Physiological Composition and Leaflet Asymmetry

The resemblance of lipid membrane models to physiological membranes determines how well molecular dynamics (MD) simulations imitate the dynamic behavior of cell membranes and membrane proteins. Physiological lipid membranes are composed of multiple types of phospholipids, and the leaflet compositions are generally asymmetric. Here we describe an approach for self-assembly of a Coarse-Grained (CG) membrane model with physiological composition and leaflet asymmetry using the MARTINI force field. An initial set-up of two boxes with different types of lipids according to the leaflet asymmetry of mammalian cell membranes stacked with 0.5 nm overlap, reliably resulted in the self-assembly of bilayer membranes with leaflet asymmetry resembling that of physiological mammalian cell membranes. Self-assembly in the presence of a fragment of the plasma membrane protein syntaxin 1A led to spontaneous specific positioning of phosphatidylionositol(4,5)bisphosphate at a positively charged stretch of syntaxin consistent with experimental data. An analogous approach choosing an initial set-up with two concentric shells filled with different lipid types results in successful assembly of a spherical vesicle with asymmetric leaflet composition. Self-assembly of the vesicle in the presence of the synaptic vesicle protein synaptobrevin 2 revealed the correct position of the synaptobrevin transmembrane domain. This is the first CG MD method to form a membrane with physiological lipid composition as well as leaflet asymmetry by self-assembly and will enable unbiased studies of the incorporation and dynamics of membrane proteins in more realistic CG membrane models.     

On the role of lipid in colicin pore formation

Insights into the protein-membrane interactions by which the C-terminal pore-forming domain of colicins inserts into membranes and forms voltage-gated channels, and the nature of the colicin channel, are provided by data on: (i) the flexible helix-elongated state of the colicin pore-forming domain in the fluid anionic membrane interfacial layer, the optimum anionic surface charge for channel formation, and voltage-gated translocation of charged regions of the colicin domain across the membrane; (ii) structure-function data on the voltage-gated K(+) channel showing translocation of an arginine-rich helical segment through the membrane; (iii) toroidal channels formed by small peptides that involve local participation of anionic lipids in an inverted phase. It is proposed that translocation of the colicin across the membrane occurs through minimization of the Born charging energy for translocation of positively charged basic residues across the lipid bilayer by neutralization with anionic lipid head groups. The resulting pore structure may consist of somewhat short, ca. 16 residues, trans-membrane helices, in a locally thinned membrane, together with surface elements of inverted phase lipid micelles.