Effect of head group and curvature on binding of the antimicrobial peptide tritrpticin to lipid membranes

In this work we examine the interaction between the 13-residue cationic antimicrobial peptide (AMP) tritrpticin (VRRFPWWWPFLRR, TRP3) and model membranes of variable lipid composition. The effect on peptide conformational properties was investigated by means of CD (circular dichroism) and fluorescence spectroscopies. Based on the hypothesis that the antibiotic acts through a mechanism involving toroidal pore formation, and taking into account that models of toroidal pores imply the formation of positive curvature, we used large unilamellar vesicles (LUV) to mimic the initial step of peptide-lipid interaction, when the peptide binds to the bilayer membrane, and micelles to mimic the topology of the pore itself, since these aggregates display positive curvature. In order to more faithfully assess the role of curvature, micelles were prepared with lysophospholipids containing (qualitatively and quantitatively) head groups identical to those of bilayer phospholipids. CD and fluorescence spectra showed that, while TRP3 binds to bilayers only when they carry negatively charged phospholipids, binding to micelles occurs irrespective of surface charge, indicating that electrostatic interactions play a less predominant role in the latter case. Moreover, the conformations acquired by the peptide were independent of lipid composition in both bilayers and micelles. However, the conformations were different in bilayers and in micelles, suggesting that curvature has an influence on the secondary structure acquired by the peptide. Fluorescence data pointed to an interfacial location of TRP3 in both types of aggregates. Nevertheless, experiments with a water soluble fluorescence quencher suggested that the tryptophan residues are more accessible to the quencher in micelles than in bilayers. Thus, we propose that bilayers and micelles can be used as models for the two steps of toroidal pore formation.     

Conformational energetics of rhodopsin modulated by nonlamellar-forming lipids

Rhodopsin is an important example of a G protein-coupled receptor (GPCR) in which 11-cis-retinal is the ligand and acts as an inverse agonist. Photolysis of rhodopsin leads to formation of the activated meta II state from its precursor meta I. Various mechanisms have been proposed to explain how the membrane composition affects the meta I-meta II conformational equilibrium in the visual process. For rod disk membranes and recombinant membranes containing rhodopsin, the lipid properties have been discussed in terms of elastic deformation of the bilayer. Here we have investigated the relation of nonlamellar-forming lipids, such as dioleoylphosphatidylethanolamine (DOPE), together with dioleoylphosphatidylcholine (DOPC), to the photochemistry of membrane-bound rhodopsin. We conducted flash photolysis experiments for bovine rhodopsin recombined with DOPE/DOPC mixtures (0:100 to 75:25) as a function of pH to explore the dependence of the photochemical activity on the monolayer curvature free energy of the membrane. It is well-known that DOPC forms bilayers, whereas DOPE has a propensity to adopt the nonlamellar, reverse hexagonal (H(II)) phase. In the case of neutral DOPE/DOPC recombinants, calculations of the membrane surface pH confirmed that an increase in DOPE favored the meta II state. Moreover, doubling the PE headgroup content versus the native rod membranes substituted for the polyunsaturated, docosahexaenoic acyl chains (22:6 omega 3), suggesting rhodopsin function is associated with a balance of forces within the bilayer. The data are interpreted by applying a flexible surface model, in which the meta II state is stabilized by lipids tending to form the H(II) phase, with a negative spontaneous curvature. A simple theory, based on principles of surface chemistry, for coupling the energetics of membrane proteins to material properties of the bilayer lipids is described. For rhodopsin, the free energy balance of the receptor and the lipids is altered by photoisomerization of retinal and involves curvature stress/strain of the membrane (frustration). A new biophysical principle is introduced: matching of the spontaneous curvature of the lipid bilayer to the mean curvature of the lipid/water interface adjacent to the protein, which balances the lipid/protein solvation energy. In this manner, the thermodynamic driving force for the meta I-meta II conformational change of rhodopsin is tightly controlled by mixtures of nonlamellar-forming lipids having distinctive material properties.     

Regulation of sodium channel function by bilayer elasticity: the importance of hydrophobic coupling. Effects of Micelle-forming amphiphiles and cholesterol

Membrane proteins are regulated by the lipid bilayer composition. Specific lipid-protein interactions rarely are involved, which suggests that the regulation is due to changes in some general bilayer property (or properties). The hydrophobic coupling between a membrane-spanning protein and the surrounding bilayer means that protein conformational changes may be associated with a reversible, local bilayer deformation. Lipid bilayers are elastic bodies, and the energetic cost of the bilayer deformation contributes to the total energetic cost of the protein conformational change. The energetics and kinetics of the protein conformational changes therefore will be regulated by the bilayer elasticity, which is determined by the lipid composition. This hydrophobic coupling mechanism has been studied extensively in gramicidin channels, where the channel-bilayer hydrophobic interactions link a "conformational" change (the monomer<-->dimer transition) to an elastic bilayer deformation. Gramicidin channels thus are regulated by the lipid bilayer elastic properties (thickness, monolayer equilibrium curvature, and compression and bending moduli). To investigate whether this hydrophobic coupling mechanism could be a general mechanism regulating membrane protein function, we examined whether voltage-dependent skeletal-muscle sodium channels, expressed in HEK293 cells, are regulated by bilayer elasticity, as monitored using gramicidin A (gA) channels. Nonphysiological amphiphiles (beta-octyl-glucoside, Genapol X-100, Triton X-100, and reduced Triton X-100) that make lipid bilayers less "stiff", as measured using gA channels, shift the voltage dependence of sodium channel inactivation toward more hyperpolarized potentials. At low amphiphile concentration, the magnitude of the shift is linearly correlated to the change in gA channel lifetime. Cholesterol-depletion, which also reduces bilayer stiffness, causes a similar shift in sodium channel inactivation. These results provide strong support for the notion that bilayer-protein hydrophobic coupling allows the bilayer elastic properties to regulate membrane protein function.     

Bilayer edge and curvature effects on partitioning of lipids by tail length: atomistic simulations

The partitioning of lipids among different microenvironments in a bilayer is of considerable relevance to characterization of composition variations in biomembranes. Atomistic simulation has been ill-suited to model equilibrated lipid mixtures because the time required for diffusive exchange of lipids among microenvironments exceeds typical submicrosecond molecular dynamics trajectories. A method to facilitate local composition fluctuations, using Monte Carlo mutations to change lipid structures within the semigrand-canonical ensemble (at a fixed difference in component chemical potentials, Δμ), was recently implemented to address this challenge. This technique was applied here to mixtures of dimyristoylphosphatidylcholine and a shorter-tail lipid, either symmetric (didecanoylphosphatidylcholine (DDPC)) or asymmetric (hexanoyl-myristoylphosphatidylcholine), arranged in two types of structure: bilayer ribbons and buckled bilayers. In ribbons, the shorter-tail component showed a clear enrichment at the highly curved rim, more so for hexanoyl-myristoylphosphatidylcholine than for DDPC. Results on buckled bilayers were variable. Overall, the DDPC content of buckled bilayers tended to exceed by several percent the DDPC content of flat ones simulated at the same Δμ, but only for mixtures with low overall DDPC content. Within the buckled bilayer structure, no correlation could be resolved between the sign or magnitude of the local curvature of a leaflet and the mean local lipid composition. Results are discussed in terms of packing constraints, surface area/volume ratios, and curvature elasticity.     

Cholesterol and phosphatidylethanolamine lipids exert opposite effects on membrane modulations caused by the M2 amphipathic helix

Membrane curvature remodeling induced by amphipathic helices (AHs) is essential in many biological processes. Here we studied a model amphipathic peptide, M2AH, derived from influenza A M2. We are interested in how M2AH may promote membrane curvature by altering membrane physical properties. We used atomic force microscopy (AFM) to examine changes in membrane topographic and mechanical properties. We used electron paramagnetic resonance (EPR) spectroscopy to explore changes in lipid chain mobility and chain orientational order. We found that M2AH perturbed lipid bilayers by generating nanoscale pits. The structural data are consistent with lateral expansion of lipid chain packing, resulting in a mechanically weaker bilayer. Our EPR spectroscopy showed that M2AH reduced lipid chain mobility and had a minimal effect on lipid chain orientational order. The EPR data are consistent with the surface-bound state of M2AH that acts as a chain mobility inhibitor. By comparing results from different lipid bilayers, we found that cholesterol enhanced the activity of M2AH in inducing bilayer pits and altering lipid chain mobility. The results were explained by considering specific M2AH-cholesterol recognition and/or cholesterol-induced expansion of interlipid distance. Both AFM and EPR experiments revealed a modest effect of anionic lipids. This highlights that membrane interaction of M2AH is mainly driven by hydrophobic forces. Lastly, we found that phosphatidylethanolamine (PE) lipids inhibited the activity of M2AH. We explained our data by considering interlipid hydrogen-bonding that can stabilize bilayer organization. Our results of lipid-dependent membrane modulations are likely relevant to M2AH-induced membrane restructuring.     

Syringomycin E channel: a lipidic pore stabilized by lipopeptide?

Highly reproducible ion channels of the lipopeptide antibiotic syringomycin E demonstrate unprecedented involvement of the host bilayer lipids. We find that in addition to a pronounced influence of lipid species on the open-channel ionic conductance, the membrane lipids play a crucial role in channel gating. The effective gating charge, which characterizes sensitivity of the conformational equilibrium of the syringomycin E channels to the transmembrane voltage, is modified by the lipid charge and lipid dipolar moment. We show that the type of host lipid determines not only the absolute value but also the sign of the gating charge. With negatively charged bilayers, the gating charge sign inverts with increased salt concentration or decreased pH. We also demonstrate that the replacement of lamellar lipid by nonlamellar with the negative spontaneous curvature inhibits channel formation. These observations suggest that the asymmetric channel directly incorporates lipids. The charges and dipoles resulting from the structural inclusion of lipids are important determinants of the overall energetics that underlies channel gating. We conclude that the syringomycin E channel may serve as a biophysical model to link studies of ion channels with those of lipidic pores in membrane fusion.     

Gramicidin channels in phospholipid bilayers with unsaturated acyl chains.

In organic solvents gramicidin A (gA) occurs as a mixture of slowly interconverting double-stranded dimers. Membrane-spanning gA channels, in contrast, are almost exclusively single-stranded beta(6,3)-helical dimers. Based on spectroscopic evidence, it has previously been concluded that the conformational preference of gA in phospholipid bilayers varies as a function of the degree of unsaturation of the acyl chains. Double-stranded pi pi(5,6)-helical dimers predominate (over single-stranded beta(6,3)-helical dimers) in lipid bilayer membranes with polyunsaturated acyl chains. We therefore examined the characteristics of channels formed by gA in 1-palmitoyl-2-oleoylphosphatidylcholine/n-decane, 1,2-dioleoylphosphatidylcholine/n-decane, and 1,2-dilinoleoylphosphatidylcholine/n-decane bilayers. We did not observe long-lived channels that could be conducting double-stranded pi pi(5,6)-helical dimers in any of these different membrane environments. We conclude that the single-stranded beta(6,3)-helical dimer is the only conducting species in these bilayers. Somewhat surprisingly, the average channel duration and channel-forming potency of gA are increased in dilinoleoylphosphatidylcholine/n-decane bilayers compared to 1-palmitoyl-2-oleoylphosphatidylcholine/n-decane and dioleoylphosphatidylcholine/n-decane bilayers. To test for specific interactions between the aromatic side chains of gA and the acyl chains of the bilayer, we examined the properties of channels formed by gramicidin analogues in which the four tryptophan residues were replaced with naphthylalanine (gN), tyrosine (gT), and phenylalanine (gM). The results show that all of these analogue channels experience the same relative stabilization when going from dioleoylphosphatidylcholine to dilinoleoylphosphatidylcholine bilayers.     

Lipid Configurations from Molecular Dynamics Simulations

The extent to which current force fields faithfully reproduce conformational properties of lipids in bilayer membranes, and whether these reflect the structural principles established for phospholipids in bilayer crystals, are central to biomembrane simulations. We determine the distribution of dihedral angles in palmitoyl-oleoyl phosphatidylcholine from molecular dynamics simulations of hydrated fluid bilayer membranes. We compare results from the widely used lipid force field of Berger et al. with those from the most recent C36 release of the CHARMM force field for lipids. Only the CHARMM force field produces the chain inequivalence with sn-1 as leading chain that is characteristic of glycerolipid packing in fluid bilayers. The exposure and high partial charge of the backbone carbonyls in Berger lipids leads to artifactual binding of Na⁺ ions reported in the literature. Both force fields predict coupled, near-symmetrical distributions of headgroup dihedral angles, which is compatible with models of interconverting mirror-image conformations used originally to interpret NMR order parameters. The Berger force field produces rotamer populations that correspond to the headgroup conformation found in a phosphatidylcholine lipid bilayer crystal, whereas CHARMM36 rotamer populations are closer to the more relaxed crystal conformations of phosphatidylethanolamine and glycerophosphocholine. CHARMM36 alone predicts the correct relative signs of the time-average headgroup order parameters, and reasonably reproduces the full range of NMR data from the phosphate diester to the choline methyls. There is strong motivation to seek further experimental criteria for verifying predicted conformational distributions in the choline headgroup, including the ³¹P chemical shift anisotropy and ¹⁴N and CD₃ NMR quadrupole splittings.     

Atomic view of calcium-induced clustering of phosphatidylserine in mixed lipid bilayers

Membranes play key regulatory roles in biological processes, with bilayer composition exerting marked effects on binding affinities and catalytic activities of a number of membrane-associated proteins. In particular, proteins involved in diverse processes such as vesicle fusion, intracellular signaling cascades, and blood coagulation interact specifically with anionic lipids such as phosphatidylserine (PS) in the presence of Ca(2+) ions. While Ca(2+) is suspected to induce PS clustering in mixed phospholipid bilayers, the detailed structural effects of this ion on anionic lipids are not established. In this study, combining magic angle spinning (MAS) solid-state NMR (SSNMR) measurements of isotopically labeled serine headgroups in mixed lipid bilayers with molecular dynamics (MD) simulations of PS lipid bilayers in the presence of different counterions, we provide site-resolved insights into the effects of Ca(2+) on the structure and dynamics of lipid bilayers. Ca(2+)-induced conformational changes of PS in mixed bilayers are observed in both liposomes and Nanodiscs, a nanoscale membrane mimetic of bilayer patches. Site-resolved multidimensional correlation SSNMR spectra of bilayers containing (13)C,(15)N-labeled PS demonstrate that Ca(2+) ions promote two major PS headgroup conformations, which are well resolved in two-dimensional (13)C-(13)C, (15)N-(13)C, and (31)P-(13)C spectra. The results of MD simulations performed on PS lipid bilayers in the presence or absence of Ca(2+) provide an atomic view of the conformational effects underlying the observed spectra.     

How cholesterol is distributed between monolayers in asymmetric lipid membranes

The distribution of cholesterol in asymmetric lipid bilayers was studied by extensive coarse-grained molecular dynamics simulations. The effects of the lipid head group charge, acyl chain saturation, spontaneous membrane curvature and surface tension of the membrane were investigated. Four asymmetric bilayers containing DOPC, DOPS, DSPC or DSPS lipids were simulated on a time scale extended to tens of microseconds. We show that cholesterol strongly prefers anionic lipids to neutral and saturated lipid tails to unsaturated with a distribution ratio of ~0.7 in neutral/anionic bilayers and of ~0.4 in unsaturated/saturated bilayers. Multiple flip-flop transitions of cholesterol were observed directly, and their mean times ranged from 80 to 250?ns. It was shown that the distribution of cholesterol in the asymmetric membrane depends not only on the type of lipid, but also on the local membrane curvature and the surface tension. The membrane curvature enhances the influence of the lipid head groups on cholesterol distribution, while non-optimal surface tension caused by different areas per lipid in different monolayers increases the effect of the lipid tail saturation. It was clearly seen that the monolayers of asymmetric bilayers are interdependent. Mean distances from the bilayer center to cholesterol molecules depend not only on the type of the lipid in the considered monolayer but also on the composition of the opposite monolayer.     

Probability of alamethicin conductance states varies with nonlamellar tendency of bilayer phospholipids.

With few exceptions, membrane lipids are usually regarded as a kind of filler or passive solvent for membrane proteins. Yet, cells exquisitely control membrane composition. Many phospholipids found in plasma membrane bilayers favor packing into inverted hexagonal bulk phases. It was suggested that the strain of forcing such lipids into a bilayer may affect membrane protein function, such as the operation of transmembrane channels. To investigate this, we have inserted the peptide alamethicin into bilayer membranes composed of lipids of empirically determined inverted hexagonal phase "spontaneous radii" Ro, which will have expectably different degrees of strain when forced into bilayer form. We observe a correlation between measured Ro and the relative probabilities of different conductance states. States of higher conductance are more probable in dioleoylphosphatidylethanolamine, the lipid of highest curvature, 1/Ro, than in dioleoylphosphatidylcholine, the lipid of lowest curvature.     

Phosphatidylethanolamine enhances rhodopsin photoactivation and transducin binding in a solid supported lipid bilayer as determined using plasmon-waveguide resonance spectroscopy

Flash photolysis studies have shown that the membrane lipid environment strongly influences the ability of rhodopsin to form the key metarhodopsin II intermediate. Here we have used plasmon-waveguide resonance (PWR) spectroscopy, an optical method sensitive to both mass and conformation, to probe the effects of lipid composition on conformational changes of rhodopsin induced by light and due to binding and activation of transducin (G(t)). Octylglucoside-solubilized rhodopsin was incorporated by detergent dilution into solid-supported bilayers composed either of egg phosphatidylcholine or various mixtures of a nonlamellar-forming lipid (dioleoylphosphatidylethanolamine; DOPE) together with a lamellar-forming lipid (dioleoylphosphatidylcholine; DOPC). Light-induced proteolipid conformational changes as a function of pH correlated well with previous flash photolysis studies, indicating that the PWR spectral shifts monitored metarhodopsin II formation. The magnitude of these effects, and hence the extent of the conformational transition, was found to be proportional to the DOPE content. Our data are consistent with previous suggestions that lipids having a negative spontaneous curvature favor elongation of rhodopsin during the activation process. In addition, measurements of the G(t)/rhodopsin interaction in a DOPC/DOPE (25:75) bilayer at pH 5 demonstrated that light activation increased the affinity for G(t) from 64 nM to 0.7 nM, whereas G(t) affinity for dark-adapted rhodopsin was unchanged. By contrast, in DOPC bilayers the affinity of G(t) for light-activated rhodopsin was only 18 nM at pH 5. Moreover exchange of GDP for GTP gamma S was also monitored by PWR spectroscopy. Only the light-activated receptor was able to induce this exchange which was unaffected by DOPE incorporation. These findings demonstrate that nonbilayer-forming lipids can alter functionally linked conformational changes of G-protein-coupled receptors in membranes, as well as their interactions with downstream effector proteins.     

Effect of lipid peroxidation on the properties of lipid bilayers: a molecular dynamics study

Lipid peroxidation plays an important role in cell membrane damage. We investigated the effect of lipid peroxidation on the properties of 1-palmitoyl-2-linoleoyl-sn-glycero-3-phosphatidylcholine (PLPC) lipid bilayers using molecular dynamics simulations. We focused on four main oxidation products of linoleic acid with either a hydroperoxide or an aldehyde group: 9-trans, cis-hydroperoxide linoleic acid, 13-trans, cis-hydroperoxide linoleic acid, 9-oxo-nonanoic acid, and 12-oxo-9-dodecenoic acid. These oxidized chains replaced the sn-2 linoleate chain. The properties of PLPC lipid bilayers were characterized as a function of the concentration of oxidized lipids, with concentrations from 2.8% to 50% for each oxidation product. The introduction of oxidized functional groups in the lipid tail leads to an important conformational change in the lipids: the oxidized tails bend toward the water phase and the oxygen atoms form hydrogen bonds with water and the polar lipid headgroup. This conformational change leads to an increase in the average area per lipid and, correspondingly, to a decrease of the bilayer thickness and the deuterium order parameters for the lipid tails, especially evident at high concentrations of oxidized lipid. Water defects are observed in the bilayers more frequently as the concentration of the oxidized lipids is increased. The changes in the structural properties of the bilayer and the water permeability are associated with the tendency of the oxidized lipid tails to bend toward the water interface. Our results suggest that one mechanism of cell membrane damage is the increase in membrane permeability due to the presence of oxidized lipids.     

Lipid lateral diffusion and membrane heterogeneity

The pulsed field gradient (pfg)-NMR method for measurements of translational diffusion of molecules in macroscopically aligned lipid bilayers is described. This technique is proposed to have an appreciable potential for investigations in the field of lipid and membrane biology. Transport of molecules in the plane of the bilayer can be successfully studied, as well as lateral phase separation of lipids and their dynamics within the bilayer organizations. Lateral diffusion coefficients depend on lipid packing and acyl chain ordering and investigations of order parameters of perdeuterated acyl chains, using (2)H NMR quadrupole splittings, are useful complements. In this review we summarize some of our recent achievements obtained on lipid membranes. In particular, bilayers exhibiting two-phase coexistence of liquid disordered (l(d)) and liquid ordered (l(o)) phases are considered in detail. Methods for obtaining good oriented lipid bilayers, necessary for the pfg-NMR method to be efficiently used, are also briefly described. Among our major results, besides determinations of l(d) and l(o) phases, belongs the finding that the lateral diffusion is the same for all components, independent of the molecular structure (including cholesterol (CHOL)), if they reside in the same domain or phase in the membrane. Furthermore, quite unexpectedly CHOL seems to partition into the l(d)and l(o) phases to roughly the same extent, indicating that CHOL has no strong preference for any of these phases, i.e. CHOL seems to have similar interactions with all of the lipids. We propose that the lateral phase separation in bilayers containing one high-T(m) and one low-T(m) lipid together with CHOL is driven by the increasing difficulty of incorporating an unsaturated or prenyl lipid into the highly ordered bilayer formed by a saturated lipid and CHOL, i.e. the phase transition is entropy driven to keep the disorder of the hydrocarbon chains of the unsaturated lipid.     

Behavioral differences between phosphatidic acid and phosphatidylcholine in the presence of the nicotinic acetylcholine receptor

It has been found experimentally that negatively charged phosphatidic acid (PA) lipids and cholesterol molecules stabilize the nicotinic acetylcholine receptor (nAChR) in a functional resting state that can participate in an agonist-induced conformational change. In this study, we compare phosphatidylcholine (PC) and PA lipid behavior in the presence of the nAChR to determine why PC lipids do not support a functional nAChR. For lipids that are located within 1.0 nm of the protein, both PC and PA lipids have very similar order-parameter and bilayer-thickness values, which indicate that the annular lipid properties are protein-dependent. The most significant difference between the PC and PA bilayers is the formation of a lipid domain around the protein, which is visible in the PA bilayer but not the PC bilayer. This suggests that the PA domain may help stabilize the nAChR resting state. The PA lipids in the microdomain have a decreased order compared to a homogeneous PA bilayer and the lipids near the protein attempt to increase the free space in their vicinity by residing in multiple lateral planes.     

Specific interactions of tryptophan with phosphatidylcholine and digalactosyldiacylglycerol in pure and mixed bilayers in the dry and hydrated state

Amphiphilic solutes play an important role in the desiccation tolerance of plant cells, because they can reversibly partition into cellular membranes during dehydration. Their effects on membrane stability depend on their chemical structure, but also on the lipid composition of the host membrane. We have shown recently that tryptophan destabilizes liposomes during freezing. The degree of destabilization depends on the presence of glycolipids in the membranes, but not on the phase preference (bilayer or non-bilayer) of the lipids in mixtures with the bilayer lipid phosphatidylcholine. Here, we have investigated the influence of tryptophan on the phase behavior and intermolecular interactions in dry and hydrated bilayers made from the phospholipid egg phosphatidylcholine and the plant chloroplast glycolipid digalactosyldiacylglycerol, or from a mixture (1:1) of these lipids, using Fourier-transform infrared spectroscopy. To distinguish effects of the hydrophobic ring structure of tryptophan from those of the amino acid moiety, we also performed experiments with the hydrophilic amino acid glycine. Our data show that there are specific interactions between tryptophan and either phospholipid or glycolipid in the dry state, as well as H-bonding interactions between the lipids and both solutes. In the rehydrated state, the H-bonding interactions between amino acids and lipids are mostly replaced by interactions between water and lipids, while the hydrophobic interactions between lipids and tryptophan mostly persist.     

Silica xerogel/aerogel-supported lipid bilayers: Consequences of surface corrugation

The objective of this paper was to review our recent investigations of silica xerogel and aerogel-supported lipid bilayers. These systems provide a format to observe relationships between substrate curvature and supported lipid bilayer formation, lipid dynamics, and lipid mixtures phase behavior and partitioning. Sensitive surface techniques such as quartz crystal microbalance and atomic force microscopy are readily applied to these systems. To inform current and future investigations, we review the experimental literature involving the impact of curvature on lipid dynamics, lipid and phase-separated lipid domain localization, and membrane-substrate conformations and we review our molecular dynamics simulations of supported lipid bilayers with the atomistic and molecular information they provide.     

Theoretical conformational analysis of phospholipids bilayers

We present a computational approach describing the conformation of lipid molecules (1-2-dipalmitoyl-sn-glycero-3 phosphocholine (DPPC)) organized in bilayers. The classical semi-empirical method used in peptide conformational analysis has been extended successfully to lipids. The excellent agreement between our theoretical predictions and recent experimental data on the molecular organization of lipid bilayers suggests that the method could be a valuable tool in the lipid conformational analysis but also in the prediction of orientation and mode of insertion of amphiphilic molecules into the lipid bilayer.     

On the role of anionic lipids in charged protein interactions with membranes

We investigate the role of anionic lipids in the binding to, and subsequent movement of charged protein groups in lipid membranes, to help understand the role of membrane composition in all membrane-active protein sequences. We demonstrate a small effect of phosphatidylglycerol (PG) lipids on the ability of an arginine (Arg) side chain to bind to, and cross a lipid membrane, despite possessing a neutralizing charge. We observe similar membrane deformations in lipid bilayers composed of phosphatidylcholine (PC) and PC/PG mixtures, with comparable numbers of water and lipid head groups pulled into the bilayer hydrocarbon core, and prohibitively large ~20 kcal/mol barriers for Arg transfer across each bilayer, dropping by just 2-3 kcal/mol due to the binding of PG lipids. We explore the causes of this small effect of introducing PG lipids and offer an explanation in terms of the limited membrane interaction for the choline groups of PC lipids bound to the translocating ion. Our calculations reveal a surprising lack of preference for Arg binding to PG lipids themselves, but a small increase in interfacial binding affinity for lipid bilayers containing PG lipids. These results help to explain the nature of competitive lipid binding to charged protein sequences, with implications for a wide range of membrane binding domains and cell perturbing peptides.     

Lipid polymorphism and biomembrane function

In recent years, several major developments have taken place in the biology, physical chemistry and technology of polymorphism of membrane lipids. These include the identification of polymorphic regulation of membrane lipid composition in Escherichia coli, the importance of nonbilayer lipids for protein functioning, the special packing properties of bilayers containing these lipids, and the crystalization of a membrane protein out of three dimensional bilayer networks (lipid cubic phases). These exciting developments bring us closer to understanding the paradox of the lipid bilayer structure of biomembranes and the molecular basis of membrane protein structure and function.