Antimicrobial peptides (AMPs) act as host defenses against microbial pathogens. Here we investigate the interactions of SVS-1 (KVKVKVKVDPLPTKVKVKVK), an engineered AMP and anti-cancer β-hairpin peptide, with lipid bilayers using spectroscopic studies and atomistic molecular dynamics simulations. In agreement with literature reports, simulation and experiment show preferential binding of SVS-1 peptides to anionic over neutral bilayers. Fluorescence and circular dichroism studies of a Trp-substituted SVS-1 analogue indicate, however, that it will bind to a zwitterionic DPPC bilayer under high-curvature conditions and folds into a hairpin. In bilayers formed from a 1:1 mixture of DPPC and anionic DPPG lipids, curvature and lipid fluidity are also observed to promote deeper insertion of the fluorescent peptide. Simulations using the CHARMM C36m force field offer complementary insight into timescales and mechanisms of folding and insertion. SVS-1 simulated at an anionic mixed POPC/POPG bilayer folded into a hairpin over a microsecond, the final stage in folding coinciding with the establishment of contact between the peptide's valine sidechains and the lipid tails through a “flip and dip” mechanism. Partial, transient folding and superficial bilayer contact are seen in simulation of the peptide at a zwitterionic POPC bilayer. Only when external surface tension is applied does the peptide establish lasting contact with the POPC bilayer. Our findings reveal the influence of disruption to lipid headgroup packing (via curvature or surface tension) on the pathway of binding and insertion, highlighting the collaborative effort of electrostatic and hydrophobic interactions on interaction of SVS-1 with lipid bilayers.
The desire to create cell-like models for fundamental science and applications has spurred extensive effort toward creating giant unilamellar vesicles (GUVs). However, a route to selectively self-assemble GUVs in bulk has remained elusive. In bulk solution, membrane-forming molecules such as phospholipids, single-tailed surfactants, and block copolymers typically self-assemble into multilamellar, onion-like structures. So although self-assembly processes can form nanoscale unilamellar vesicles, scaffolding by droplets or surfaces is required to create GUVs. Here we show that it is possible to bulk self-assemble cell-sized GUVs with almost complete selectivity over other vesicle topologies. The seemingly paradoxical pair of features that enables this appears to be having very dynamic molecules at the nanoscale that create unusually rigid membranes. The resultant self-assembly pathway enables encapsulation of molecules and colloids and can also generate model primitive cells that can grow and divide.
Abstract
Mixtures of long- and short-tail phosphatidylcholine lipids are known to self-assemble into a variety of aggregates combining flat bilayerlike and curved micellelike features, commonly called bicelles. Atomistic simulations of bilayer ribbons and perforated bilayers containing dimyristoylphosphatidylcholine (DMPC, di-C14 tails) and dihexanoylphosphatidylcholine (DHPC, di-C6 tails) have been carried out to investigate the partitioning of these components between flat and curved microenvironments and the stabilization of the bilayer edge by DHPC. To approach equilibrium partitioning of lipids on an achievable simulation timescale, configuration-bias Monte Carlo mutation moves were used to allow individual lipids to change tail length within a semigrand-canonical ensemble. Since acceptance probabilities for direct transitions between DMPC and DHPC were negligible, a third component with intermediate tail length (didecanoylphosphatidylcholine, di-C10 tails) was included at a low concentration to serve as an intermediate for transitions between DMPC and DHPC. Strong enrichment of DHPC is seen at ribbon and pore edges, with an excess linear density of ∼3 nm−1. The simulation model yields estimates for the onset of edge stability with increasing bilayer DHPC content between 5% and 15% DHPC at 300 K and between 7% and 17% DHPC at 323 K, higher than experimental estimates. Local structure and composition at points of close contact between pores suggest a possible mechanism for effective attractions between pores, providing a rationalization for the tendency of bicelle mixtures to aggregate into perforated vesicles and perforated sheets.
To understand the effects of lipid composition on membrane protein function in a mixture as complex as a biomembrane, one must know whether the lipid composition local to the protein differs from the mean lipid composition. In this study, we simulated the transmembrane domain of a β-barrel protein, OmpA, in mixtures of lipids of different tail lengths under conditions of negative hydrophobic mismatch, i.e., local bilayer thinning. We modeled the influence of OmpA on the local lipid composition both at a coarse-grained (CG) resolution using conventional molecular dynamics, and at an atomistic resolution within the semi-grand canonical ensemble using mutation moves to rapidly approach an equilibrium lateral distribution of lipids. Moderate enrichment of the shorter tail component (either DDPC in DDPC/DMPC mixtures or DMPC in DMPC/DSPC mixtures) extending 2–3 nm away from the protein surface was observed with both the atomistic and CG models. The similarity in trends suggests that the more computationally economical CG models capture the essential features of lipid sorting induced by hydrophobic mismatch.
The dynamics of the gel to fluid phase transformation in 100 nm large unilamellar vesicles (LUV) of 1,2-dipalmitoyl(d62)-sn-glycero-3-phosphocholine (d62-DPPC), has been studied by laser-induced temperature-jump initiation coupled with time-resolved infrared spectroscopy and by MD simulations. The infrared transients that characterize the temperature dependent phase transformation are complex, extending from the nanosecond to the millisecond time scales. An initial fast (submicrosecond) component can be modeled by partial melting of the gel domains, initiated at pre-existing defects at the edges of the faceted structure of the gel phase. Molecular dynamics simulations support the model of fast melting from edge defects. The extent of melting during the fast phase is limited by the area expansion on melting, which leads to a surface pressure that raises the effective melting temperature. Subsequent melting is observed to follow highly stretched exponential kinetics, consistent with collective relaxation of the surface pressure through a hierarchy of surface undulations with different relaxation times. The slowest step is water diffusion through the bilayer to allow the vesicle volume to grow along with its expanded surface area. The results demonstrate that the dominant relaxation in the gel to fluid phase transformation in response to a large T-jump perturbation (compared to the transition width) is fast (submicrosecond), which has important practical and fundamental consequences.