Molecular dynamics simulations of mixed micelles modeling human bile

Extensive molecular dynamics (MD) simulations of binary systems of phospholipids and bile salts, a model for human bile, have been performed. Recent progress in hardware and software development allows simulation of the spontaneous aggregation of the constituents into small mixed micelles, in agreement with experimental observations. The MD simulations reveal the structure of these micelles at atomic detail. The phospholipids are packed radially with their headgroups at the surface and the hydrophobic tails pointing toward the micellar center. The bile salts act as wedges between the phospholipid headgroups, with their hydrophilic sides exposed to the aqueous environment. The structure of the micelles strongly resembles the previously proposed radial shell model. Simulations including small fractions of cholesterol reveal how cholesterol is solubilized inside these mixed micelles without changing their overall structure.     

Partitioning of lipids at domain boundaries in model membranes

Line-active molecules (“linactants”) that bind to the boundary interface between different fluid lipid domains in membranes have a strong potential as regulators of the lateral heterogeneity that is important for many biological processes. Here, we use molecular dynamics simulations in combination with a coarse-grain model that retains near-atomic resolution to identify lipid species that can act as linactants in a model membrane that is segregated into two lipid domains of different fluidity. Our simulations predict that certain hybrid saturated/unsaturated chain lipids can bind to the interface and lower the line tension, whereas cone-shaped lysolipids have a less pronounced effect.     

Mechanosensitive membrane channels in action

The tension-driven gating process of MscL from Mycobacterium tuberculosis, Tb-MscL, has been addressed at near-atomic detail using coarse-grained molecular dynamics simulations. To perform the simulations, a novel coarse-grained peptide model based on a thermodynamic parameterization of the amino-acid side chains has been applied. Both the wild-type Tb-MscL and its gain-of-function mutant V21D embedded in a solvated lipid bilayer have been studied. To mimic hypoosmotic shock conditions, simulations were performed at increasing levels of membrane tension approaching the rupture threshold of the lipid bilayer. Both the wild-type and the mutant channel are found to undergo significant conformational changes in accordance with an irislike expansion mechanism, reaching a conducting state on a microsecond timescale. The most pronounced expansion of the pore has been observed for the V21D mutant, which is consistent with the experimentally shown gain-of-function phenotype of the V21D mutant.     

Areas of Monounsaturated Diacylphosphatidylcholines

We have studied the structural properties of monounsaturated diacylphosphatidylcholine lipid bilayers (i.e., diCn:1PC, where n = 14, 16, 18, 20, 22, and 24 is the number of acyl chain carbons). High-resolution x-ray scattering data were analyzed in conjunction with contrast-varied neutron scattering data using a technique we recently developed. Analyses of the data show that the manner by which bilayer thickness increases with increasing n in monounsaturated diacylphosphatidylcholines is dependent on the double bond's position. For commonly available monounsaturated diacylphosphatidylcholines, this results in the nonlinear behavior of both bilayer thickness and lipid area, whereas for diC18:1PC bilayers, lipid area assumes a maximum value. It is worthwhile to note that compared to previous data, our results indicate that lipid areas are smaller by ∼10%. This observation highlights the need to revisit lipid areas, as they are often used in comparisons with molecular dynamics simulations. Moreover, simulators are encouraged to compare their results not only to x-ray scattering data, but to neutron data as well.     

Structural basis for the enhanced activity of cyclic antimicrobial peptides: the case of BPC194

We report the molecular basis for the differences in activity of cyclic and linear antimicrobial peptides. We iteratively performed atomistic molecular dynamics simulations and biophysical measurements to probe the interaction of a cyclic antimicrobial peptide and its inactive linear analogue with model membranes. We establish that, relative to the linear peptide, the cyclic one binds stronger to negatively charged membranes. We show that only the cyclic peptide folds at the membrane interface and adopts a β-sheet structure characterised by two turns. Subsequently, the cyclic peptide penetrates deeper into the bilayer while the linear peptide remains essentially at the surface. Finally, based on our comparative study, we propose a model characterising the mode of action of cyclic antimicrobial peptides. The results provide a chemical rationale for enhanced activity in certain cyclic antimicrobial peptides and can be used as a guideline for design of novel antimicrobial peptides.     

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_BT 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_BT on the gating energy of MscL. This contribution arises largely from the interfacial tension between hydrophobic and hydrophilic regions in a lipid bilayer.     

Dual Action of BPC194: A Membrane Active Peptide Killing Bacterial Cells

Membrane active peptides can perturb the lipid bilayer in several ways, such as poration and fusion of the target cell membrane, and thereby efficiently kill bacterial cells. We probe here the mechanistic basis of membrane poration and fusion caused by membrane-active, antimicrobial peptides. We show that the cyclic antimicrobial peptide, BPC194, inhibits growth of Gram-negative bacteria and ruptures the outer and inner membrane at the onset of killing, suggesting that not just poration is taking place at the cell envelope. To simplify the system and to better understand the mechanism of action, we performed Förster resonance energy transfer and cryogenic transmission electron microscopy studies in model membranes and show that the BPC194 causes fusion of vesicles. The fusogenic action is accompanied by leakage as probed by dual-color fluorescence burst analysis at a single liposome level. Atomistic molecular dynamics simulations reveal how the peptides are able to simultaneously perturb the membrane towards porated and fused states. We show that the cyclic antimicrobial peptides trigger both fusion and pore formation and that such large membrane perturbations have a similar mechanistic basis.     

Proteolytic activity at neutral pH in bovine spleen

Proteolytic activity at neutral pH can be demonstrated in extracts from beef spleen. This activity is completely due to an aminopeptidase - or a number of aminopeptides - which is able to hydrolyse proteins completely to amino acids. No evidence was found for the presence of endopeptidases active at neutral pH. The enzyme resembles to some extent swine kidney aminopeptides.     

Alternative Mechanisms for the Interaction of the Cell-Penetrating Peptides Penetratin and the TAT Peptide with Lipid Bilayers

Cell-penetrating peptides (CPPs) have recently attracted much interest due to their apparent ability to penetrate cell membranes in an energy-independent manner. Here molecular-dynamics simulation techniques were used to study the interaction of two CPPs: penetratin and the TAT peptide with 1,2-Dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) and 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) phospolipid bilayers shed light on alternative mechanisms by which these peptides might cross biological membranes. In contrast to previous simulation studies of charged peptides interacting with lipid bilayers, no spontaneous formation of transmembrane pores was observed. Instead, the simulations suggest that the peptides may enter the cell by micropinocytosis, whereby the peptides induce curvature in the membrane, ultimately leading to the formation of small vesicles within the cell that encapsulate the peptides. Specifically, multiple peptides were observed to induce large deformations in the lipid bilayer that persisted throughout the timescale of the simulations (hundreds of nanoseconds). Pore formation could be induced in simulations in which an external potential was used to pull a single penetratin or TAT peptide into the membrane. With the use of umbrella-sampling techniques, the free energy of inserting a single penetratin peptide into a DPPC bilayer was estimated to be ∼75 kJmol⁻¹, which suggests that the spontaneous penetration of single peptides would require a timescale of at least seconds to minutes. This work also illustrates the extent to which the results of such simulations can depend on the initial conditions, the extent of equilibration, the size of the system, and the conditions under which the simulations are performed. The implications of this with respect to the current systems and to simulations of membrane-peptide interactions in general are discussed.     

Lipids on the move: Simulations of membrane pores, domains, stalks and curves

In this review we describe the state-of-the-art of computer simulation studies of lipid membranes. We focus on collective lipid-lipid and lipid-protein interactions that trigger deformations of the natural lamellar membrane state, showing that many important biological processes including self-aggregation of membrane components into domains, the formation of non-lamellar phases, and membrane poration and curving, are now amenable to detailed simulation studies.