Structure of fully hydrated fluid phase lipid bilayers with monounsaturated chains

Molecular insight into the effect of lipid bilayer environments on thrombospondin-1 and calreticulin interactions

Thrombospondin-1 (TSP1) binding to cell surface calreticulin (CRT) stimulates the association of CRT with low-density lipoprotein (LDL) receptor-related protein (LRP1) to signal focal adhesion disassembly and engagement of cellular activities. A recent study demonstrated that membrane rafts are necessary for TSP1-mediated focal adhesion disassembly, but the molecular role of membrane rafts in mediating TSP1-CRT-LRP1 signaling is unknown. In this study, we investigated the effect of lipid bilayer environments on TSP1 and CRT interactions via atomically detailed molecular dynamics simulations. Results showed that the microscopic structural properties of lipid molecules and mesoscopic mechanical properties and electrostatic potential of the bilayer were significantly different between a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer and a raftlike lipid bilayer [a POPC/cholesterol (CHOL) raftlike lipid bilayer or a POPC/CHOL/sphingomyelin (SM) raftlike lipid bilayer], and the difference was enhanced by SM lipids in a raftlike lipid bilayer. These bilayer property differences affect the interactions of CRT with the bilayer, further influencing CRT conformation and TSP1-CRT interactions. A raftlike lipid bilayer stabilized CRT conformation as compared to a POPC bilayer environment. TSP1 binding to CRT resulted in a conformation for the CRT N-domain more "open" than that of the CRT P-domain in a raftlike lipid bilayer environment, which could facilitate binding of CRT to LRP1 to engage downstream signaling. The open conformational changes of CRT by binding to TSP1 in a raftlike lipid bilayer were enhanced by SM lipids in a lipid bilayer. The direct interactions of both the N- and P-domains of CRT with the bilayer contribute to the more open conformation of CRT in the TSP1-CRT complex on a raftlike lipid bilayer as compared to that on a POPC bilayer. The interactions of CRT or the TSP1-CRT complex with the lipid bilayer also caused CHOL molecules and/or lipids to be more coordinated and to aggregate into patchlike regions in the raftlike lipid bilayers. The lipid and CHOL molecule coordination and aggregation could in turn affect the interactions of CRT with the membrane raft, thereby altering TSP1-CRT interactions and CRT conformational changes that potentially regulate its interactions with LRP1. This study provides molecular insights into the role of lipid bilayer environments in TSP1-CRT interactions and in the CRT conformational changes that are predicted to facilitate binding of CRT to LRP1 to engage downstream signaling events.

Two homologous neutrophil serine proteases bind to POPC vesicles with different affinities: When aromatic amino acids matter

Neutrophil serine proteases Proteinase 3 (PR3) and human neutrophil elastase (HNE) are homologous antibiotic serine proteases of the polymorphonuclear neutrophils. Despite sharing a 56% sequence identity they have been shown to have different functions and localizations in the neutrophils. In particular, and in contrast to HNE, PR3 has been detected at the outer leaflet of the plasma membrane and its membrane expression is a risk factor in a number of chronic inflammatory diseases. Although a plethora of studies performed in various cell-based assays have been reported, the mechanism by which PR3, and possibly HNE bind to simple membrane models remains unclear. We used surface plasmon resonance (SPR) experiments to measure and compare the affinity of PR3 and HNE for large unilamellar vesicles composed of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC). We also conducted 500-nanosecond long molecular dynamics simulations of each enzyme at the surface of a POPC bilayer to map the interactions between proteins and lipids and rationalize the difference in affinity observed in the SPR experiment. We find that PR3 binds strongly to POPC large unilamellar vesicles (Kd=9.2×10(-7)M) thanks to the insertion of three phenylalanines, one tryptophan and one leucine beyond the phosphate groups of the POPC lipids. HNE binds in a significantly weaker manner (Kd>10(-5)M) making mostly electrostatic interactions via lysines and arginines and inserting only one leucine between the hydrophobic lipid tails. Our results support the early reports that PR3, unlike HNE, is able to directly and strongly anchor directly to the neutrophil membrane.

Zooming in on disordered systems: neutron reflection studies of proteins associated with fluid membranes

Neutron reflectometry (NR) is an emerging experimental technique for the structural characterization of proteins interacting with fluid bilayer membranes under conditions that mimic closely the cellular environment. Thus, cellular processes can be emulated in artificial systems and their molecular basis studied by adding cellular components one at a time in a well-controlled environment while the resulting structures, or structural changes in response to external cues, are monitored with neutron reflection. In recent years, sample environments, data collection strategies and data analysis were continuously refined. The combination of these improvements increases the information which can be obtained from NR to an extent that enables structural characterization of protein-membrane complexes at a length scale that exceeds the resolution of the measurement by far. Ultimately, the combination of NR with molecular dynamics (MD) simulations can be used to cross-validate the results of the two techniques and provide atomic-scale structural models. This review discusses these developments in detail and demonstrates how they provide new windows into relevant biomedical problems. This article is part of a Special Issue entitled: Interfacially Active Peptides and Proteins. Guest Editors: William C. Wimley and Kalina Hristova.

Phase coexistence in single-lipid membranes induced by buffering agents

Recent literature has shown that buffers affect the interaction between lipid bilayers through a mechanism that involves van der Waals forces, electrostatics, hydration forces and membrane bending rigidity. This letter shows an additional peculiar effect of buffers on the mixed chain 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipid bilayers, namely phase coexistence similar to what was reported by Rappolt et al. for alkali chlorides. The data presented suggest that one phase appears to dehydrate below the value in pure water, while the other phase swells as the concentration of buffer is increased. However, since the two phases must be in osmotic equilibrium with one another, this behavior challenges theoretical models of lipid interactions.

On calculation of the electrostatic potential of a phosphatidylinositol phosphate-containing phosphatidylcholine lipid membrane accounting for membrane dynamics

Many signaling events require the binding of cytoplasmic proteins to cell membranes by recognition of specific charged lipids, such as phosphoinositol-phosphates. As a model for a protein-membrane binding site, we consider one charged phosphoinositol phosphate (PtdIns(3)P) embedded in a phosphatidylcholine bilayer. As the protein-membrane binding is driven by electrostatic interactions, continuum solvent models require an accurate representation of the electrostatic potential of the phosphoinositol phosphate-containing membrane. We computed and analyzed the electrostatic potentials of snapshots taken at regular intervals from molecular dynamics simulations of the bilayer. We observe considerable variation in the electrostatic potential of the bilayer both along a single simulation and between simulations performed with the GAFF or CHARMM c36 force fields. However, we find that the choice of GAFF or CHARMM c36 parameters has little effect on the electrostatic potential of a given configuration of the bilayer with a PtdIns(3)P embedded in it. From our results, we propose a remedian averaging method for calculating the electrostatic potential of a membrane system that is suitable for simulations of protein-membrane binding with a continuum solvent model.

HIV-1 Tat membrane interactions probed using X-ray and neutron scattering, CD spectroscopy and MD simulations

We report the effect on lipid bilayers of the Tat peptide Y47GRKKRRQRRR57 from the HIV-1 virus transactivator of translation (Tat) protein. Synergistic use of low-angle X-ray scattering (LAXS) and atomistic molecular dynamic simulations (MD) indicate Tat peptide binding to neutral dioleoylphosphocholine (DOPC) lipid headgroups. This binding induced the local lipid phosphate groups to move 3Å closer to the center of the bilayer. Many of the positively charged guanidinium components of the arginines were as close to the center of the bilayer as the locally thinned lipid phosphate groups. LAXS data for DOPC, DOPC/dioleoylphosphoethanolamine (DOPE), DOPC/dioleoylphosphoserine (DOPS), and a mimic of the nuclear membrane gave similar results. Generally, the Tat peptide decreased the bilayer bending modulus KC and increased the area/lipid. Further indications that Tat softens a membrane, thereby facilitating translocation, were provided by wide-angle X-ray scattering (WAXS) and neutron scattering. CD spectroscopy was also applied to further characterize Tat/membrane interactions. Although a mechanism for translation remains obscure, this study suggests that the peptide/lipid interaction makes the Tat peptide poised to translocate from the headgroup region.

Interfacial partitioning of a loop hinge residue contributes to diacylglycerol affinity of conserved region 1 domains

Conventional and novel isoenzymes of PKC are activated by the membrane-embedded second messenger diacylglycerol (DAG) through its interactions with the C1 regulatory domain. The affinity of C1 domains to DAG varies considerably among PKCs. To gain insight into the origin of differential DAG affinities, we conducted high-resolution NMR studies of C1B domain from PKCδ (C1Bδ) and its W252Y variant. The W252Y mutation was previously shown to render C1Bδ less responsive to DAG (Dries, D. R., Gallegos, L. L., and Newton, A. C. (2007) A single residue in the C1 domain sensitizes novel protein kinase C isoforms to cellular diacylglycerol production. J. Biol. Chem. 282, 826-830) and thereby emulate the behavior of C1B domains from conventional PKCs that have a conserved Tyr at the equivalent position. Our data revealed that W252Y mutation did not perturb the conformation of C1Bδ in solution but significantly reduced its propensity to partition into a membrane-mimicking environment in the absence of DAG. Using detergent micelles doped with a paramagnetic lipid, we determined that both the residue identity at position 252 and complexation with diacylglycerol influence the geometry of C1Bδ-micelle interactions. In addition, we identified the C-terminal helix α1 of C1Bδ as an interaction site with the head groups of phosphatidylserine, a known activator of PKCδ. Taken together, our studies (i) reveal the identities of C1Bδ residues involved in interactions with membrane-mimicking environment, DAG, and phosphatidylserine, as well as the affinities associated with each event and (ii) suggest that the initial ligand-independent membrane recruitment of C1B domains, which is greatly facilitated by the interfacial partitioning of Trp-252, is responsible, at least in part, for the differential DAG affinities.

Highly hydrated deformable polyethylene glycol-tethered lipid bilayers

The realization of a solid-supported lipid bilayer acting as a workbench for the study of membrane processes is a difficult task. For robustness, the bilayer has to be tethered to the substrate. At the same time, diffusion of the lipids and plastic deformations of the membrane should not be obstructed. Furthermore, a highly hydrated surrounding is mandatory. Here, we show that grafting of a polyethylene glycol-lipid construct (PEG2000-DSPE) to a silicon oxide surface via multiple-step silane chemistry and subsequent deposition of lipids by spin-coating result in a cushioned membrane that has the desired properties. Neutron and X-ray reflectometry measurements are combined to access thickness, density, and hydration of the bilayer and the PEG cushion. We observe a spacer of 55 Å thickness between lipid bilayer and silicon-oxide surface with a rather high hydration of up to 90 ± 3% water. While 11.5 ± 3% of the lipids are grafted to the surface, as determined from the neutron data, the diffusion constant of the lipids, as probed by diffusion of 0.5% Texas Red labeled lipids, remains rather large (D = 2.1 ± 0.1 μm(2)/s), which is a reduction of only 12% compared to a supported lipid bilayer reference without immobilized lipids. Finally, AFM indentation confirms the plastic behavior of the membrane against deformation. We show that rupture of the bilayer does not occur before the deformation exceeds 40 Å. Altogether, the presented PEG-tethered lipid bilayer mimics the deformability of natural cell membranes much better than standard solid-supported lipid bilayers.

Interaction of amphotericin B with lipid monolayers

Langmuir isotherm, neutron reflectivity, and Brewster angle microscopy experiments have been performed to study the interaction of amphotericin B (AmB) with monolayers prepared from 1-palmitoyl-2-oleoylphosphatidylcholine (POPC) and mixtures of this lipid with cholesterol or ergosterol to mimic mammalian and fungal cell membranes, respectively. Isotherm data show that AmB causes a more pronounced change in surface pressure in the POPC/ergosterol system than in the POPC and POPC/cholesterol systems, and its interaction with the POPC/ergosterol monolayer is also more rapid than with the POPC and POPC/cholesterol monolayers. Brewster angle microscopy shows that, in interaction with POPC monolayers, AmB causes the formation of small domains which shrink and disappear within a few minutes. The drug also causes domain formation in the POPC/cholesterol and POPC/ergosterol monolayers; in the former case, these are formed more slowly than is seen with the POPC monolayers and are ultimately much smaller; in the latter case, they are formed rather more quickly and are more heterogeneous in size. Neutron reflectivity data show that the changes in monolayer structure following interaction with AmB are the same for all three systems studied: the data are consistent with the drug inserting into the monolayers with its macrocyclic ring intercalated among the lipid acyl chains and sterol ring systems, with its mycosamine moiety colocalizing with the sterol hydroxyl and POPC head groups. On the basis of these studies, it is concluded that AmB inserts in a similar manner into POPC, POPC/cholesterol, and POPC/ergosterol monolayers but does so with differing kinetics and with the formation of quite different in-plane structures. The more rapid time scale for interaction of the drug with the POPC/ergosterol monolayer, its more pronounced effect on monolayer surface pressure, and its more marked changes as regards domain formation are all consistent with the drug's selectivity for fungal vs mammalian cell membranes.

Molecular simulation of ibuprofen passing across POPC membrane

Permeability assessment is an important procedure in the drug development process, and drug partitioning in membrane bilayer is related to permeability. To investigate the pH dependence on drug partitioning, the process of different ionization state of ibuprofen passing across POPC bilayer was studied using molecular dynamics simulation. The results show that both atomic charge scheme and ionization state of the drug affect the value and shape of energy profile when passing across the POPC bilayer. The neutral ibuprofen (ibuprofen under acidic condition) has a much lower energy barrier as compared with the anionic ibuprofen (ibuprofen under basic condition). Meantime, hydrogen bond analysis also certifies that it is easy for neutral ibuprofen to pass from bulk water to bilayer center. Our calculation suggests that the ionization state of ibuprofen may be changed between neutral and anionic state when passing across membrane: it may be ionized outside the membrane and neutralized inside the membrane.

A systematic molecular dynamics simulation study of temperature dependent bilayer structural properties

Although lipid force fields (FFs) used in molecular dynamics (MD) simulations have proved to be accurate, there has not been a systematic study on their accuracy over a range of temperatures. Motivated by the X-ray and neutron scattering measurements of common phosphatidylcholine (PC) bilayers (Kučerka et al. BBA. 1808: 2761, 2011), the CHARMM36 (C36) FF accuracy is tested in this work with MD simulations of six common PC lipid bilayers over a wide range of temperatures. The calculated scattering form factors and deuterium order parameters from the C36 MD simulations agree well with the X-ray, neutron, and NMR experimental data. There is excellent agreement between MD simulations and experimental estimates for the surface area per lipid, bilayer thickness (DB), hydrophobic thickness (DC), and lipid volume (VL). The only minor discrepancy between simulation and experiment is a measure of (DB-DHH)/2 where DHH is the distance between the maxima in the electron density profile along the bilayer normal. Additional MD simulations with pure water and heptane over a range of temperatures provide explanations of possible reasons causing the minor deviation. Overall, the C36 FF is accurate for use with liquid crystalline PC bilayers of varying chain types and over biologically relevant temperatures.

Peripheral and integral membrane binding of peptides characterized by time-dependent fluorescence shifts: focus on antimicrobial peptide LAH₄

Positioning of peptides with respect to membranes is an important parameter for biological and biophysical studies using model systems. Our experiments using five different membrane peptides suggest that the time-dependent fluorescence shift (TDFS) of Laurdan can help when distinguishing between peripheral and integral membrane binding and can be a useful, novel tool for studying the impact of transmembrane peptides (TMP) on membrane organization under near-physiological conditions. This article focuses on LAH4, a model α-helical peptide with high antimicrobial and nucleic acid transfection efficiencies. The predominantly helical peptide has been shown to orient in supported model membranes parallel to the membrane surface at acidic and, in a transmembrane manner, at basic pH. Here we investigate its interaction with fully hydrated large unilamellar vesicles (LUVs) by TDFS and fluorescence correlation spectroscopy (FCS). TDFS shows that at acidic pH LAH4 does not influence the glycerol region while at basic pH it makes acyl groups at the glycerol level of the membrane less mobile. TDFS experiments with antimicrobial peptides alamethicin and magainin 2, which are known to assume transmembrane and peripheral orientations, respectively, prove that changes in acyl group mobility at the glycerol level correlate with the orientation of membrane-associated peptide molecules. Analogous experiments with the TMPs LW21 and LAT show similar effects on the mobility of those acyl groups as alamethicin and LAH4 at basic pH. FCS, on the same neutral lipid bilayer vesicles, shows that the peripheral binding mode of LAH4 is more efficient in bilayer permeation than the transmembrane mode. In both cases, the addition of LAH4 does not lead to vesicle disintegration. The influence of negatively charged lipids on the bilayer permeation is also addressed.

Recent theoretical advances in elasticity of membranes following Helfrich's spontaneous curvature model

Recent theoretical advances in elasticity of membranes following Helfrich's famous spontaneous curvature model are summarized in this review. The governing equations describing equilibrium configurations of lipid vesicles, lipid membranes with free edges, and chiral lipid membranes are presented. Several analytic solutions to these equations and their corresponding configurations are demonstrated.

Recent developments in the field of bending rigidity measurements on membranes

This review gives a brief overview of experimental approaches used to assess the bending rigidity of membranes. Emphasis is placed on techniques based on the use of giant unilamellar vesicles. We summarize the effect on the bending rigidity of membranes as a function of membrane composition, presence of various inclusions in the bilayer and molecules and ions in the bathing solutions. Examples for the impact of temperature, cholesterol, some peptides and proteins, sugars and salts are provided and the literature data are discussed critically. Future directions, open questions and possible developments in this research field are also included.

Molecular dynamics simulations of homo-oligomeric bundles embedded within a lipid bilayer

Using molecular dynamics simulations, we studied the structure, interhelix interactions, and dynamics of transmembrane proteins. Specifically, we investigated homooligomeric helical bundle systems consisting of synthetic α-helices with either the sequence Ac-(LSLLLSL)3-NH2 (LS2) or Ac-(LSSLLSL)3-NH2 (LS3). The LS2 and LS3 helical peptides are designed to have amphipathic characteristics that form ion channels in membrane. We simulated bundles containing one to six peptides that were embedded in palmitoyl-oleoyl-phosphatidylcholine (POPC) lipid bilayer and placed between two lamellae of water. We aim to provide a fundamental understanding of how amphipathic helical peptides interact with each other and their dynamical behaviors in different homooligomeric states. To understand structural properties, we examined the helix lengths, tilt angles of individual helices and the entire bundle, interhelix distances, interhelix cross-angles, helix hydrophobic-to-hydrophilic vector projections, and the average number of interhelix hydrophilic (serine-serine) contacts lining the pore of the transmembrane channel. To analyze dynamical properties, we calculated the rotational autocorrelation function of each helix and the cross-correlation of the rotational velocity between adjacent helices. The observed structural and dynamical characteristics show that higher order bundles containing four to six peptides are composed of multiple lower order bundles of one to three peptides. For example, the LS2 channel was found to be stable in a tetrameric bundle composed of a "dimer of dimers." In addition, we observed that there is a minimum of two strong hydrophilic contacts between a pair of adjacent helices in the dimer to tetramer systems and only one strong hydrophilic interhelix contact in helix pairs of the pentamer and hexamer systems. We believe these results are general and can be applied to more complex ion channels, providing insight into ion channel stability and assembly.

What are the true values of the bending modulus of simple lipid bilayers?

Values of the bending modulus KC are reviewed, and possible causes for the considerable differences are discussed. One possible cause is the use of glucose and sucrose in the classical micromechanical manipulation and shape analysis methods. New data, using the more recent low angle X-ray method, are presented that do not support an effect of glucose or sucrose on KC. Another possible cause is using an incomplete theory to interpret the data. Adding a tilt term to the theory clearly does not affect the value obtained from the shape analysis method. It is shown that a tilt term, using a value of the modulus Kθ indicated by simulations, theory, and estimated from order parameters obtained from NMR and from the wide angle X-ray method, should also not affect the value obtained using the micromechanical manipulation method, although it does require a small correction when determining the value of the area compressibility modulus KA. It is still being studied whether including a tilt term will significantly affect the values of KC obtained using low angle X-ray data. It remains unclear what causes the differences in the experimental values of KC for simple lipid bilayers.

Improved Coarse-Grained Modeling of Cholesterol-Containing Lipid Bilayers

Cholesterol trafficking, which is an essential function in mammalian cells, is intimately connected to molecular-scale interactions through cholesterol modulation of membrane structure and dynamics and interaction with membrane receptors. Since these effects of cholesterol occur on micro- to millisecond timescales, it is essential to develop accurate coarse-grained simulation models that can reach these timescales. Cholesterol has been shown experimentally to thicken the membrane and increase phospholipid tail order between 0-40% cholesterol, above which these effects plateau or slightly decrease. Here, we showed that the published MARTINI coarse-grained force-field for phospholipid (POPC) and cholesterol fails to capture these effects. Using reference atomistic simulations, we systematically modified POPC and cholesterol bonded parameters in MARTINI to improve its performance. We showed that the corrections to pseudo-bond angles between glycerol and the lipid tails and around the oleoyl double bond particle (the "angle-corrected model") slightly improves the agreement of MARTINI with experimentally measured thermal, elastic, and dynamic properties of POPC membranes. The angle-corrected model improves prediction of the thickening and ordering effects up to 40% cholesterol but overestimates these effects at higher cholesterol concentration. In accordance with prior work that showed the cholesterol rough face methyl groups are important for limiting cholesterol self-association, we revised the coarse-grained representation of these methyl groups to better match cholesterol-cholesterol radial distribution functions from atomistic simulations. In addition, by using a finer-grained representation of the branched cholesterol tail than MARTINI, we improved predictions of lipid tail order and bilayer thickness across a wide range of concentrations. Finally, transferability testing shows that a model incorporating our revised parameters into DOPC outperforms other CG models in a DOPC/cholesterol simulation series, which further argues for its efficacy and generalizability. These results argue for the importance of systematic optimization for coarse-graining biologically important molecules like cholesterol with complicated molecular structure.

Membrane structure correlates to function of LLP2 on the cytoplasmic tail of HIV-1 gp41 protein

Mutation studies previously showed that the lentivirus lytic peptide (LLP2) sequence of the cytoplasmic C-terminal tail of the HIV-1 gp41 envelope protein inhibited viral-initiated T-cell death and T-cell syncytium formation, at which time in the HIV life cycle the gp41 protein is embedded in the T-cell membrane. In striking contrast, the mutants did not affect virion infectivity, during which time the gp41 protein is embedded in the HIV envelope membrane. To examine the role of LLP2/membrane interactions, we applied synchrotron x-radiation to determine structure of hydrated membranes. We focused on WT LLP2 peptide (+3 charge) and MX2 mutant (-1 charge) with membrane mimics for the T-cell and the HIV-1 membranes. To investigate the influence of electrostatics, cholesterol content, and peptide palmitoylation, we also studied three other LLP2 variants and HIV-1 mimics without negatively charged lipids or cholesterol as well as extracted HIV-1 lipids. All LLP2 peptides bound strongly to T-cell membrane mimics, as indicated by changes in membrane structure and bending. In contrast, none of the weakly bound LLP2 variants changed the HIV-1 membrane mimic structure or properties. This correlates well with, and provides a biophysical basis for, previously published results that reported lack of a mutant effect in HIV virion infectivity in contrast to an inhibitory effect in T-cell syncytium formation. It shows that interaction of LLP2 with the T-cell membrane modulates biological function.

Solid-state NMR studies of full-length BamA in lipid bilayers suggest limited overall POTRA mobility

The outer membrane protein BamA is the key player in β-barrel assembly in Gram-negative bacteria. Despite the availability of high-resolution crystal structures, the dynamic behavior of the transmembrane domain and the large periplasmic extension consisting of five POTRA (POlypeptide-TRansport-Associated) domains remains unclear. We demonstrate reconstitution of full-length BamA in proteoliposomes at low lipid-to-protein ratio, leading to high sensitivity and resolution in solid-state NMR (ssNMR) experiments. We detect POTRA domains in ssNMR experiments probing rigid protein segments in our preparations. These results suggest that the periplasmic region of BamA is firmly attached to the β-barrel and does not experience fast global motion around the angle between POTRA 2 and 3. We show that this behavior holds at lower protein concentrations and elevated temperatures. Chemical shift variations observed after reconstitution in lipids with different chain lengths and saturation levels are compatible with conformational plasticity of BamA's transmembrane domain. Electron microscopy of the ssNMR samples shows that BamA can cause local disruptions of the lipid bilayer in proteoliposomes. The observed interplay between protein-protein and protein-lipid interactions may be critical for BamA-mediated insertion of substrates into the outer membrane.

X-ray structure, thermodynamics, elastic properties and MD simulations of cardiolipin/dimyristoylphosphatidylcholine mixed membranes

Cardiolipins (CLs) are important biologically for their unique role in biomembranes that couple phosphorylation and electron transport like bacterial plasma membranes, chromatophores, chloroplasts and mitochondria. CLs are often tightly coupled to proteins involved in oxidative phosphorylation. The first step in understanding the interaction of CL with proteins is to obtain the pure CL structure, and the structure of mixtures of CL with other lipids. In this work we use a variety of techniques to characterize the fluid phase structure, material properties and thermodynamics of mixtures of dimyristoylphosphatidylcholine (DMPC) with tetramyristoylcardiolipin (TMCL), both with 14-carbon chains, at several mole percentages. X-ray diffuse scattering was used to determine structure, including bilayer thickness and area/lipid, the bending modulus, KC, and SXray, a measure of chain orientational order. Our results reveal that TMCL thickens DMPC bilayers at all mole percentages, with a total increase of ∼6 Å in pure TMCL, and increases AL from 64 Å(2) (DMPC at 35 °C) to 109 Å(2) (TMCL at 50 °C). KC increases by ∼50%, indicating that TMCL stiffens DMPC membranes. TMCL also orders DMPC chains by a factor of ∼2 for pure TMCL. Coarse grain molecular dynamics simulations confirm the experimental thickening of 2 Å for 20mol% TMCL and locate the TMCL headgroups near the glycerol-carbonyl region of DMPC; i.e., they are sequestered below the DMPC phosphocholine headgroup. Our results suggest that TMCL plays a role similar to cholesterol in that it thickens and stiffens DMPC membranes, orders chains, and is positioned under the umbrella of the PC headgroup. CL may be necessary for hydrophobic matching to inner mitochondrial membrane proteins. Differential scanning calorimetry, SXray and CGMD simulations all suggest that TMCL does not form domains within the DMPC bilayers. We also determined the gel phase structure of TMCL, which surprisingly displays diffuse X-ray scattering, like a fluid phase lipid. AL=40.8 Å(2) for the ½TMCL gel phase, smaller than the DMPC gel phase with AL=47.2 Å(2), but similar to AL of DLPE=41 Å(2), consistent with untilted chains in gel phase TMCL.

Conformation and lipid interaction of the fusion peptide of the paramyxovirus PIV5 in anionic and negative-curvature membranes from solid-state NMR

Viral fusion proteins catalyze the merger of the virus envelope and the target cell membrane through multiple steps of protein conformational changes. The fusion peptide domain of these proteins is important for membrane fusion, but how it causes membrane curvature and dehydration is still poorly understood. We now use solid-state NMR spectroscopy to investigate the conformation, topology, and lipid and water interactions of the fusion peptide of the PIV5 virus F protein in three lipid membranes, POPC/POPG, DOPC/DOPG, and DOPE. These membranes allow us to investigate the effects of lipid chain disorder, membrane surface charge, and intrinsic negative curvature on the fusion peptide structure. Chemical shifts and spin diffusion data indicate that the PIV5 fusion peptide is inserted into all three membranes but adopts distinct conformations: it is fully α-helical in the POPC/POPG membrane, adopts a mixed strand/helix conformation in the DOPC/DOPG membrane, and is primarily a β-strand in the DOPE membrane. (31)P NMR spectra show that the peptide retains the lamellar structure and hydration of the two anionic membranes. However, it dehydrates the DOPE membrane, destabilizes its inverted hexagonal phase, and creates an isotropic phase that is most likely a cubic phase. The ability of the β-strand conformation of the fusion peptide to generate negative Gaussian curvature and to dehydrate the membrane may be important for the formation of hemifusion intermediates in the membrane fusion pathway.

Lipid14: The Amber Lipid Force Field

The AMBER lipid force field has been updated to create Lipid14, allowing tensionless simulation of a number of lipid types with the AMBER MD package. The modular nature of this force field allows numerous combinations of head and tail groups to create different lipid types, enabling the easy insertion of new lipid species. The Lennard-Jones and torsion parameters of both the head and tail groups have been revised and updated partial charges calculated. The force field has been validated by simulating bilayers of six different lipid types for a total of 0.5 μs each without applying a surface tension; with favorable comparison to experiment for properties such as area per lipid, volume per lipid, bilayer thickness, NMR order parameters, scattering data, and lipid lateral diffusion. As the derivation of this force field is consistent with the AMBER development philosophy, Lipid14 is compatible with the AMBER protein, nucleic acid, carbohydrate, and small molecule force fields.

The effect of hydroxylated fatty acid-containing phospholipids in the remodeling of lipid membranes

The synthetic fatty acid 2-hydroxyoleic acid (2OHOA) is an antitumor drug that regulates membrane lipid composition and structure. An important effect of this drug is the restoration of sphingomyelin (SM) levels in cancer cell membranes, where the SM concentration is lower than in non-tumor cells. It is well known that free fatty acid concentration in cell membranes is lower than 5%, and that fatty acid excess is rapidly incorporated into phospholipids. In a recent work, we have considered the effect of free 2OHOA in model membranes in liquid ordered (Lo) and liquid disordered (Ld) phases, by using all-atom molecular dynamics. This study concerns membranes that are modified upon incorporation of 2OHOA into different phospholipids. 2OHOA-containing phospholipids have a permanent effect on lipid membranes, making a Ld membrane surface more compact and less hydrated, whereas the opposite effect is observed in Lo domains. Moreover, the hydroxyl group of fatty acid chains increases the propensity of Ld model membranes to form hexagonal or other non-lamellar structures. This article is part of a Special Issue entitled: Membrane Structure and Function: Relevance in the Cell's Physiology, Pathology and Therapy.

Structural properties of POPC monolayers under lateral compression: computer simulations analysis

1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), a lipid comprising a saturated and an unsaturated acyl chain, belongs to the class of glycerophosphatidylcholines, major lipids in eukaryotic cell membranes. To get insight into the structural properties of this lipid within monolayers as membrane models, we performed molecular dynamics (MD) simulations of POPC monolayers under compression at the air/water interface. MD simulations were carried out at 300 K and at different surface pressures using the all-atom general Amber force field (GAFF). A good agreement was found between the simulated data and experimental isotherms. At surface pressures greater than 15 mN/m, two orientations of the head groups clearly appear: one nearly parallel to the monolayer interface and another one pointing toward the water. On the basis of the analysis of headgroup dihedral angles, we propose that the conformational variations around the bonds connecting the phosphorus atom to the adjacent oxygens are involved in these two orientations of the headgroup. The glycerol group orientation is characterized by a large distribution centered around 50° with respect to the monolayer normal. The acyl chains are predominantly in trans configuration from 7.5 to 43 mN/m surface pressures. Moreover, the calculated order parameter profiles of both chains suggest an independent behavior of the saturated and unsaturated chains that could be correlated with the formation of chain-type clusters observed along the simulated trajectories.

Molecular conformation and bilayer pores in a nonionic surfactant lamellar phase studied with 1H-13C solid-state NMR and molecular dynamics simulations

The structure of the lamellar phase of aqueous pentaethylene glycol mono-n-dodecyl ether (C12E5) surfactant at various temperatures and molar fractions is studied by using united atom molecular dynamics simulations and nuclear magnetic resonance measurements. Namely, the simulation model is used to interpret the magnitude and temperature dependence of experimental C-H order parameter profiles in terms of the molecular conformation and orientation. Our simulations suggest that the low order parameters that are generally measured in poly(ethylene oxide) surfactant bilayers are due to the presence of bilayer pores throughout the entire lamellar phase region.

Buffers affect the bending rigidity of model lipid membranes

In biophysical and biochemical studies of lipid bilayers the influence of the used buffer is often ignored or assumed to be negligible on membrane structure, elasticity, or physical properties. However, we here present experimental evidence, through bending rigidity measurements performed on giant vesicles, of a more complex behavior, where the buffering molecules may considerably affect the bending rigidity of phosphatidylcholine bilayers. Furthermore, a synergistic effect on the bending modulus is observed in the presence of both salt and buffer molecules, which serves as a warning to experimentalists in the data interpretation of their studies, since typical lipid bilayer studies contain buffer and ion molecules.

Differential effect of 2-hydroxyoleic acid enantiomers on protein (sphingomyelin synthase) and lipid (membrane) targets

The complex dual mechanism of action of 2-hydroxyoleic acid (2OHOA), a potent anti-tumor compound used in membrane lipid therapy (MLT), has yet to be fully elucidated. It has been demonstrated that 2OHOA increases the sphingomyelin (SM) cell content via SM synthase (SGMS) activation. Its presence in membranes provokes changes in the membrane lipid structure that induce the translocation of PKC to the membrane and the subsequent overexpression of CDK inhibitor proteins (e.g., p21(Cip1)). In addition, 2OHOA also induces the translocation of Ras to the cytoplasm, provoking the silencing of MAPK and its related pathways. These two differential modes of action are triggered by the interactions of 2OHOA with either lipids or proteins. To investigate the molecular basis of the different interactions of 2OHOA with membrane lipids and proteins, we synthesized the R and S enantiomers of this compound. A molecular dynamics study indicated that both enantiomers interact similarly with lipid bilayers, which was further confirmed by X-ray diffraction studies. By contrast, only the S enantiomer was able to activate SMS in human glioma U118 cells. Moreover, the anti-tumor efficacy of the S enantiomer was greater than that of the R enantiomer, as the former can act through both MLT mechanisms. The present study provides additional information on this novel therapeutic approach and on the magnitude of the therapeutic effects of type-1 and type-2 MLT approaches. This article is part of a Special Issue entitled: Membrane Structure and Function: Relevance in the Cell's Physiology, Pathology and Therapy.

Bridging molecular docking to membrane molecular dynamics to investigate GPCR-ligand recognition: the human A₂A adenosine receptor as a key study

G protein-coupled receptors (GPCRs) represent the largest family of cell-surface receptors and about one-third of the actual targets of clinically used drugs. Following the progress made in the field of GPCRs structural determination, docking-based screening for novel potent and selective ligands is becoming an increasingly adopted strategy in the drug discovery process. However, this methodology is not yet able to anticipate the "bioactive" binding mode and discern it among other conformations. In the present work, we present a novel approach consisting in the integration of molecular docking and membrane MD simulations with the aim to merge the rapid sampling of ligand poses into in the binding site, typical of docking algorithms, with the thermodynamic accuracy of MD simulations in describing, at the molecular level, the stability a GPCR-ligand complex embedded into explicit lipid-water environment. To validate our approach, we have chosen as a key study the human A(2A) adenosine receptor (hA(2A) AR) and selected four receptor-antagonist complexes and one receptor-agonist complex that have been recently crystallized. In light of the obtained results, we believe that our novel strategy can be extended to other GPCRs and might represent a valuable tool to anticipate the "bioactive" conformation of high-affinity ligands.

CHARMM36 united atom chain model for lipids and surfactants

Molecular simulations of lipids and surfactants require accurate parameters to reproduce and predict experimental properties. Previously, a united atom (UA) chain model was developed for the CHARMM27/27r lipids (Hénin, J., et al. J. Phys. Chem. B. 2008, 112, 7008-7015) but suffers from the flaw that bilayer simulations using the model require an imposed surface area ensemble, which limits its use to pure bilayer systems. A UA-chain model has been developed based on the CHARMM36 (C36) all-atom lipid parameters, termed C36-UA, and agreed well with bulk, lipid membrane, and micelle formation of a surfactant. Molecular dynamics (MD) simulations of alkanes (heptane and pentadecane) were used to test the validity of C36-UA on density, heat of vaporization, and liquid self-diffusion constants. Then, simulations using C36-UA resulted in accurate properties (surface area per lipid, X-ray and neutron form factors, and chain order parameters) of various saturated- and unsaturated-chain bilayers. When mixed with the all-atom cholesterol model and tested with a series of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC)/cholesterol mixtures, the C36-UA model performed well. Simulations of self-assembly of a surfactant (dodecylphosphocholine, DPC) using C36-UA suggest an aggregation number of 53 ± 11 DPC molecules at 0.45 M of DPC, which agrees well with experimental estimates. Therefore, the C36-UA force field offers a useful alternative to the all-atom C36 lipid force field by requiring less computational cost while still maintaining the same level of accuracy, which may prove useful for large systems with proteins.

Voltage-dependent structural changes of the membrane-bound anion channel hVDAC1 probed by SEIRA and electrochemical impedance spectroscopy

SEIRA spectroscopy provides insight into the potential-induced structural changes of the anion channel hVDAC1 embedded in the tethered bilayer lipid membrane.

A microscopic multiphase diffusion model of viable epidermis permeability

A microscopic model of passive transverse mass transport of small solutes in the viable epidermal layer of human skin is formulated on the basis of a hexagonal array of cells (i.e., keratinocytes) bounded by 4-nm-thick, anisotropic lipid bilayers and separated by 1-μm layers of extracellular fluid. Gap junctions and tight junctions with adjustable permeabilities are included to modulate the transport of solutes with low membrane permeabilities. Two keratinocyte aspect ratios are considered to represent basal and spinous cells (longer) and granular cells (more flattened). The diffusion problem is solved in a unit cell using a coordinate system conforming to the hexagonal cross section, and an efficient two-dimensional treatment is applied to describe transport in both the cell membranes and intercellular spaces, given their thinness. Results are presented in terms of an effective diffusion coefficient, D¯(epi), and partition coefficient, K¯(epi/w), for a homogenized representation of the microtransport problem. Representative calculations are carried out for three small solutes-water, L-glucose, and hydrocortisone-covering a wide range of membrane permeability. The effective transport parameters and their microscopic interpretation can be employed within the context of existing three-layer models of skin transport to provide more realistic estimates of the epidermal concentrations of topically applied solutes.

A microscopic view of phospholipid insertion into biological membranes

Understanding the process of membrane insertion is an essential step in developing a detailed mechanism, for example, for peripheral membrane protein association and membrane fusion. The highly mobile membrane mimetic (HMMM) has been used to accelerate the membrane association and binding of peripheral membrane proteins in simulations by increasing the lateral diffusion of phospholipid headgroups while retaining an atomistic description of the interface. Through a comparative study, we assess the difference in insertion rate of a free phospholipid into an HMMM as well as into a conventional phospholipid bilayer and develop a detailed mechanistic model of free phospholipid insertion into biological membranes. The mechanistic insertion model shows that successful irreversible association of the free phospholipid to the membrane interface, which results in its insertion, is the rate-limiting step. Association is followed by independent, sequential insertion of the acyl tails of the free phospholipid into the membrane, with splayed acyl tail intermediates. Use of the HMMM is found to replicate the same intermediate insertion states as in the full phospholipid bilayer; however, it accelerates overall insertion by approximately a factor of 3, with the probability of successful association of phospholipid to the membrane being significantly enhanced.

Monolayer spontaneous curvature of raft-forming membrane lipids

Monolayer spontaneous curvatures for cholesterol, DOPE, POPE, DOPC, DPPC, DSPC, POPC, SOPC, and egg sphingomyelin were obtained using small-angle X-ray scattering (SAXS) on inverted hexagonal phases (HII). Spontaneous curvatures of bilayer forming lipids were estimated by adding controlled amounts to a HII forming template following previously established protocols. Spontaneous curvatures of both phosphatidylethanolamines and cholesterol were found to be at least a factor of two more negative than those of phosphatidylcholines, whose J0 values are closer to zero. Interestingly, a significant positive J0 value was retrieved for DPPC. We further determined the temperature dependence of the spontaneous curvatures J0(T) in the range from 15 to 55 °C, resulting in a quite narrow distribution of -1 to -3 × 10-3 (nm °C)-1 for most investigated lipids. The data allowed us to estimate the monolayer spontaneous curvatures of ternary lipid mixtures showing liquid ordered/liquid disordered phase coexistence. We report spontaneous curvature phase diagrams for DSPC/DOPC/Chol, DPPC/DOPC/Chol and SM/POPC/Chol and discuss effects on protein insertion and line tension.

Self-assembly of squalene-based nucleolipids: relating the chemical structure of the bioconjugates to the architecture of the nanoparticles

Squalene-based nucleolipids, including anticancer or antiviral prodrugs, gave rise to nanoparticles displaying a diversity of structures upon nanoprecipitation in water. Synchrotron small-angle X-ray scattering and cryo-TEM imaging revealed that both the nature of the nucleoside and the position of the squalene moiety relative to the nucleobase determined the self-assembly of the corresponding bioconjugates. It was found that small chemical differences resulted in major differences in the self-organization of nucleolipids when squalene was grafted onto the nucleobase whereas only lamellar phases were observed when squalene was linked to the sugar moiety. The key role of hydrogen bonds between nucleobases in the formation of the lamellar phases was suggested, in agreement with molecular simulations. These findings provide a way to fine tune the supramolecular organization of squalene-based prodrugs, with the aim of improving their pharmacological activity.

Effects of PEGylation on the binding interaction of magainin 2 and tachyplesin I with lipid bilayer surface

Poly(ethylene glycol) (PEG)-grafted magainin 2 and tachyplesin I were simulated with lipid bilayers. In the simulations of PEGylated magainin 2 and tachyplesin I in water, both peptides are wrapped by PEG chains. The α-helical structure of PEGylated magainin 2 is broken in water, while the β-sheet of PEGylated tachyplesin I keeps stable, similar to the structural behavior of unPEGylated peptides, in agreement with experiments. Simulations of PEGylated peptides with lipid bilayers show that PEG chains block the electrostatic interaction between cationic residues of peptides and anionic phosphates of lipids, leading to the less binding of the peptide to the bilayer surface, which is observed more significantly for magainin 2 than for tachyplesin I. Since the random-coiled magainin 2 can be more completely covered by PEGs than does the β-sheet tachyplesin I, the PEGylation effect on the decreased binding is larger for magainin 2, showing the dependence of PEGylation on the peptide structure. These simulation findings qualitatively support the experimental observation of the different extents of the reduced membrane-permeabilizing activity for PEGylated magainin 2 and tachyplesin I.

Influence of ester-modified lipids on bilayer structure

Lipid membranes function as barriers for cells to prevent unwanted chemicals from entering the cell and wanted chemicals from leaving. Because of their hydrophobic interior, membranes do not allow water to penetrate beyond the headgroup region. We performed molecular simulations to examine the effects of ester-modified lipids, which contain ester groups along their hydrocarbon chains, on bilayer structure. We chose two lipids from those presented in Menger et al. [J. Am. Chem. Soc. 2006, 128, 14034] with ester groups in (1) the upper half of the lipid chain (MEPC) and (2) the middle and end of the lipid chain (MGPC). MGPC (30%)/POPC bilayers formed stable water pores of diameter 5-7 Å, but MGPC (22%)/POPC and MEPC (30%)/POPC bilayers did not form these defects. These pores were similar to those formed during electroporation; i.e., the head groups lined the pore and allowed water and ions to transport across the bilayer. However, we found that lateral organization of the MGPC lipids into clusters, instead of an electric field or charge disparity as in electroporation, was essential for pore formation. On the basis of this, we propose an overall mechanism for pore formation. The similarities between the ester-modified lipids and byproducts of lipid peroxidation with multiple hydrophilic groups in the middle of the chain suggest that free radical reactions with unsaturated lipids and sterols result in fundamental changes that may be similar to what is seen in bilayers with ester-modified lipids.

APL@Voro: a Voronoi-based membrane analysis tool for GROMACS trajectories

APL@Voro is a new program developed to aid in the analysis of GROMACS trajectories of lipid bilayer simulations. It can read a GROMACS trajectory file, a PDB coordinate file, and a GROMACS index file to create a two-dimensional geometric representation of a bilayer. Voronoi diagrams and Delaunay triangulations--generated for different selection models of lipids--support the analysis of the bilayer. The values calculated on the geometric structures can be visualized in a user-friendly interactive environment and, then, plotted and exported to different file types. APL@Voro supports complex bilayers with a mix of various lipids and proteins. For the calculation of the projected area per lipid, a modification of the well-known Voronoi approach is presented as well as the presentation of a new approach for including atoms into an existing triangulation. The application of the developed software is discussed for three example systems simulated with GROMACS. The program is written in C++, is open source, and is available free of charge.

Mechanism of membrane permeation induced by synthetic β-hairpin peptides

We have investigated the membrane destabilizing properties of synthetic amphiphilic cationic peptides, MAX1 and MAX35, which have the propensity to form β-hairpin structures under certain conditions, and a control non-β-hairpin-forming peptide MAX8V16E. All three peptides bind to liposomes containing a mixture of zwitterionic POPC and negatively charged POPS lipids as determined by Zeta potential measurements. Circular dichroism measurements indicated folding of MAX1 and MAX35 in the presence of the POPC/POPS liposomes, whereas no such folding was observed with MAX8V16E. There was no binding or folding of these peptides to liposomes containing only POPC. MAX1 and MAX35 induced release of contents from negatively charged liposomes, whereas MAX8V16E failed to promote solute release under identical conditions. Thus, MAX1 and MAX35 bind to, and fold at the surface of negatively charged liposomes adopting a lytic conformation. We ruled out leaky fusion as a mechanism of release by including 2 mol % PEG-PE in the liposomes, which inhibits aggregation/fusion but not folding of MAX or MAX-induced leakage. Using a concentration-dependent quenching probe (calcein), we determined that MAX-induced leakage of liposome contents was an all-or-none process. At MAX1 concentrations, which cause release of ~50% of the liposomes that contain small (R(h) <1.5 nm) markers, only ~15% of those liposomes release a fluorescent dextran of 40 kDa. A multimeric model of the pore is presented based on these results. Atomistic molecular dynamics simulations show that barrels consisting of 10 β-hairpin MAX1 and MAX35 peptides are relatively more stable than MAX8V16E barrels in the bilayer, suggesting that barrels of this size are responsible for the peptides lytic action.

Adsorption of α-synuclein to supported lipid bilayers: positioning and role of electrostatics

An amyloid form of the protein α-synuclein is the major component of the intraneuronal inclusions called Lewy bodies, which are the neuropathological hallmark of Parkinson's disease (PD). α-Synuclein is known to associate with anionic lipid membranes, and interactions between aggregating α-synuclein and cellular membranes are thought to be important for PD pathology. We have studied the molecular determinants for adsorption of monomeric α-synuclein to planar model lipid membranes composed of zwitterionic phosphatidylcholine alone or in a mixture with anionic phosphatidylserine (relevant for plasma membranes) or anionic cardiolipin (relevant for mitochondrial membranes). We studied the adsorption of the protein to supported bilayers, the position of the protein within and outside the bilayer, and structural changes in the model membranes using two complementary techniques-quartz crystal microbalance with dissipation monitoring, and neutron reflectometry. We found that the interaction and adsorbed conformation depend on membrane charge, protein charge, and electrostatic screening. The results imply that α-synuclein adsorbs in the headgroup region of anionic lipid bilayers with extensions into the bulk but does not penetrate deeply into or across the hydrophobic acyl chain region. The adsorption to anionic bilayers leads to a small perturbation of the acyl chain packing that is independent of anionic headgroup identity. We also explored the effect of changing the area per headgroup in the lipid bilayer by comparing model systems with different degrees of acyl chain saturation. An increase in area per lipid headgroup leads to an increase in the level of α-synuclein adsorption with a reduced water content in the acyl chain layer. In conclusion, the association of α-synuclein to membranes and its adsorbed conformation are of electrostatic origin, combined with van der Waals interactions, but with a very weak correlation to the molecular structure of the anionic lipid headgroup. The perturbation of the acyl chain packing upon monomeric protein adsorption favors association with unsaturated phospholipids preferentially found in the neuronal membrane.

Proapoptotic Bax and Bak proteins form stable protein-permeable pores of tunable size

The Bcl-2 proapoptotic proteins Bax and Bak mediate the permeabilization of the mitochondrial outer membrane during apoptosis. Current models consider that Bax and Bak form pores at the mitochondrial outer membrane that are responsible for the release of cytochrome c and other larger mitochondrial apoptotic factors (i.e. Smac/DIABLO, AIF, and endoglycosidase G). However, the properties and nature of Bax/Bak apoptotic pores remain enigmatic. Here, we performed a detailed analysis of the membrane permeabilizing activity of Bax and Bak at the single vesicle level. We directly visualized that cBid-activated Bax and BakΔC21 can form membrane pores large enough to release not only cytochrome c, but also allophycocyanine, a protein of 104 kDa. Interestingly, the size of Bax and BakΔC21 pores is not constant, as typically observed in purely proteinaceous channels, but evolves with time and depends on protein concentration. We found that Bax and BakΔC21 formed long-lived pores, whose areas changed with the amount of Bax/BakΔC21 but not with cardiolipin concentration. Altogether, our results demonstrate that Bax and BakΔC21 follow similar mechanisms of membrane permeabilization characterized by the formation of protein-permeable pores of dynamic size, in agreement with the proteolipidic nature of these apoptotic pores.

Familial Alzheimer's disease Osaka mutant (ΔE22) β-barrels suggest an explanation for the different Aβ1-40/42 preferred conformational states observed by experiment

An unusual ΔE693 mutation in the amyloid precursor protein (APP) producing a β-amyloid (Aβ) peptide lacking glutamic acid at position 22 (Glu22) was recently discovered, and dabbed the Osaka mutant (ΔE22). Previously, several point mutations in the Aβ peptide involving Glu22 substitutions were identified and implicated in the early onset of familial Alzheimer's disease (FAD). Despite the absence of Glu22, the Osaka mutant is also associated with FAD, showing a recessive inheritance in families affected by the disease. To see whether this aggregation-prone Aβ mutant could directly relate to the Aβ ion channel-mediated mechanism as observed for the wild type (WT) Aβ peptide in AD pathology, we modeled Osaka mutant β-barrels in a lipid bilayer. Using molecular dynamics (MD) simulations, two conformer ΔE22 barrels with the U-shaped monomer conformation derived from NMR-based WT Aβ fibrils were simulated in explicit lipid environment. Here, we show that the ΔE22 barrels obtain the lipid-relaxed β-sheet channel topology, indistinguishable from the WT Aβ1-42 barrels, as do the outer and pore dimensions of octadecameric (18-mer) ΔE22 barrels. Although the ΔE22 barrels lose the cationic binding site in the pore which is normally provided by the negatively charged Glu22 side chains, the mutant pores gain a new cationic binding site by Glu11 at the lower bilayer leaflet, and exhibit ion fluctuations similar to the WT barrels. Of particular interest, this deletion mutant suggests that toxic WT Aβ1-42 would preferentially adopt a less C-terminal turn similar to that observed for Aβ17-42, and explains why the solid state NMR data for Aβ1-40 point to a more C-terminal turn conformation. The observed ΔE22 barrels conformational preferences also suggest an explanation for the lower neurotoxicity in rat primary neurons as compared to WT Aβ1-42.

Molecular basis of cannabinoid CB1 receptor coupling to the G protein heterotrimer Gαiβγ: identification of key CB1 contacts with the C-terminal helix α5 of Gαi

The cannabinoid (CB1) receptor is a member of the rhodopsin-like G protein-coupled receptor superfamily. The human CB1 receptor, which is among the most expressed receptors in the brain, has been implicated in several disease states, including drug addiction, anxiety, depression, obesity, and chronic pain. Different classes of CB1 agonists evoke signaling pathways through the activation of specific subtypes of G proteins. The molecular basis of CB1 receptor coupling to its cognate G protein is unknown. As a first step toward understanding CB1 receptor-mediated G protein signaling, we have constructed a ternary complex structural model of the CB1 receptor and Gi heterotrimer (CB1-Gi), guided by the x-ray structure of β2-adrenergic receptor (β2AR) in complex with Gs (β2AR-Gs), through 824-ns duration molecular dynamics simulations in a fully hydrated 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine bilayer environment. We identified a group of residues at the juxtamembrane regions of the intracellular loops 2 and 3 (IC2 and IC3) of the CB1 receptor, including Ile-218(3.54), Tyr-224(IC2), Asp-338(6.30), Arg-340(6.32), Leu-341(6.33), and Thr-344(6.36), as potential key contacts with the extreme C-terminal helix α5 of Gαi. Ala mutations of these residues at the receptor-Gi interface resulted in little G protein coupling activity, consistent with the present model of the CB1-Gi complex, which suggests tight interactions between CB1 and the extreme C-terminal helix α5 of Gαi. The model also suggests that unique conformational changes in the extreme C-terminal helix α5 of Gα play a crucial role in the receptor-mediated G protein activation.

Partitioning of oleic acid into phosphatidylcholine membranes is amplified by strain

Partitioning of fatty acids into phospholipid membranes is studied on giant unilamellar vesicles (GUVs) utilizing phase-contrast microscopy. With use of a micropipet, an individual GUV is transferred from a vesicle suspension in a mixed glucose/sucrose solution into an isomolar glycerol solution with a small amount of oleic acid added. Oleic acid molecules intercalate into the phospholipid membrane and thus increase the membrane area, while glycerol permeates into the vesicle interior and thus via osmotic inflation causes an increase of the vesicle volume. The conditions are chosen at which a vesicle swells as a sphere. At sufficiently low oleic acid concentrations, when the critical membrane strain is reached, the membrane bursts and part of vesicle content is ejected, upon which the membrane reseals and the swelling commences again. The radius of the vesicle before and after the burst is determined at different concentrations of oleic acid in suspension. The results of our experiments show that the oleic acid partitioning increases when the membrane strain is increased. The observed behavior is interpreted on the basis of a tension-dependent intercalation of oleic acid into the membrane.