Structure of fully hydrated fluid phase lipid bilayers with monounsaturated chains

Losartan's affinity to fluid bilayers modulates lipid-cholesterol interactions

Losartan is an angiotensin II receptor antagonist mainly used for the regulation of high blood pressure. Since it was anticipated that losartan reaches the receptor site via membrane diffusion, the impact of losartan on model membranes has been investigated by small angle X-ray scattering. For this purpose 2-20 mol% losartan was incorporated into dimyristoyl-phosphatidylcholine (DMPC) and palmitoyl-oleoyl-phosphatidylcholine (POPC) bilayers and into their binary mixtures with cholesterol in the concentration range of 0 to 40 mol%. Effects of losartan on single component bilayers are alike. Partitioning of losartan into the membranes confers a negative charge to the lipid bilayers that causes the formation of unilamellar vesicles and a reduction of the bilayer thickness by 3-4%. Analysis of the structural data resulted in an estimate for the partial area of losartan, A(Los) ≈ 40 Å(2). In the presence of cholesterol, differences between the effects of losartan on POPC and DMPC are striking. Membrane condensation by cholesterol is retarded by losartan in POPC. This contrasts with DMPC, where an increase of the cholesterol content shifts the partitioning equilibrium of losartan towards the aqueous phase, such that losartan gets depleted from the bilayers from 20 mol% cholesterol onwards. This indicates (i) a chain-saturation dependent competition of losartan with lipid-cholesterol interactions, and (ii) the insolubility of losartan in the liquid ordered phase of PCs. Consequently, losartan's action is more likely to take place in fluid plasma membrane patches rather than in domains rich in cholesterol and saturated lipid species such as in membrane rafts.

Comparative molecular dynamics simulations of the antimicrobial peptide CM15 in model lipid bilayers

We report altogether 3-μs molecular dynamics (MD) simulations of the antimicrobial peptide CM15 to systematically investigate its interaction with two model lipid bilayers, pure POPC and mixed POPG:POPC (1:2). Starting with either an α-helical or a random-coil conformation, CM15 is found to insert into both bilayers. Peptide-lipid interaction is stronger with the anionic POPG:POPC than the zwitterionic POPC, which is largely attributed to the electrostatic attraction between CM15 and the negatively charged POPG. Simulations initiated with CM15 as a random coil allowed us to study peptide folding at the lipid-water interface. Interestingly, CM15 folding appears to be faster in POPC than POPG:POPC, which may be explained by a lower activation energy barrier of structural rearrangement in the former system. Our data also suggest that compared with the random-coil conformation, CM15 in a pre-folded α-helix has significantly reduced interactions with the lipids, indicating that peptide initial structures may bias the simulation results considerably on the 100-ns timescale. The implications of this result should be considered when preparing and interpreting future AMP simulations.

Membrane Protein Simulations Using AMBER Force Field and Berger Lipid Parameters

AMBER force fields are among the most commonly used in molecular dynamics (MD) simulations of proteins. Unfortunately, they lack a specific set of lipid parameters, thus limiting its use in membrane protein simulations. In order to overcome this limitation we assessed whether the widely used united-atom lipid parameters described by Berger and co-workers could be used in conjunction with AMBER force fields in simulations of membrane proteins. Thus, free energies of solvation in water and in cyclohexane, and free energies of water to cyclohexane transfer, were computed by thermodynamic integration procedures for neutral amino acid side-chains employing AMBER99, AMBER03, and OPLS-AA amino acid force fields. In addition, MD simulations of three membrane proteins in a POPC lipid bilayer, the β2 adrenergic G protein-coupled receptor, Aquaporin-1, and the outer membrane protein Omp32, were performed with the aim of comparing the AMBER99SB/Berger combination of force fields with the OPLS-AA/Berger combination. We have shown that AMBER99SB and Berger force fields are compatible, they provide reliable free energy estimations relative to experimental values, and their combination properly describes both membrane and protein structural properties. We then suggest that the AMBER99SB/Berger combination is a reliable choice for the simulation of membrane proteins, which links the easiness of ligand parametrization and the ability to reproduce secondary structure of AMBER99SB force field with the largely validated Berger lipid parameters.

Bottom-up design and synthesis of limit size lipid nanoparticle systems with aqueous and triglyceride cores using millisecond microfluidic mixing

Limit size systems are defined as the smallest achievable aggregates compatible with the packing of the molecular constituents in a defined and energetically stable structure. Here we report the use of rapid microfluidic mixing for the controlled synthesis of two types of limit size lipid nanoparticle (LNP) systems, having either polar or nonpolar cores. Specifically, limit size LNP consisting of 1-palmitoyl, 2-oleoyl phosphatidylcholine (POPC), cholesterol and the triglyceride triolein were synthesized by mixing a stream of ethanol containing dissolved lipid with an aqueous stream, employing a staggered herringbone micromixer. Millisecond mixing of aqueous and ethanol streams at high flow rate ratios (FRR) was used to rapidly increase the polarity of the medium, driving bottom-up synthesis of limit size LNP systems by spontaneous assembly. For POPC/triolein systems the limit size structures consisted of a hydrophobic core of triolein surrounded by a monolayer of POPC where the diameter could be rationally engineered over the range 20-80 nm by varying the POPC/triolein ratio. In the case of POPC and POPC/cholesterol (55/45; mol/mol) the limit size systems achieved were bilayer vesicles of approximately 20 and 40 nm diameter, respectively. We further show that doxorubicin, a representative weak base drug, can be efficiently loaded and retained in limit size POPC LNP, establishing potential utility as drug delivery systems. To our knowledge this is the first report of stable triglyceride emulsions in the 20-50 nm size range, and the first time vesicular systems in the 20-50 nm size range have been generated by a scalable manufacturing method. These results establish microfluidic mixing as a powerful and general approach to access novel LNP systems, with both polar or nonpolar core structures, in the sub-100 nm size range.

Modification of Lipid Bilayer Structure by Diacylglycerol: A Comparative Study of Diacylglycerol and Cholesterol

Diacylglycerols (DAGs) are important second messengers in biomembranes, and they can activate protein kinase C and many other enzymes and receptors. However, their interactions with cholesterol and other lipids have not been previously studied using molecular dynamics (MD) simulation. In this study, nine independent atomistic MD simulations were performed to specifically investigate the interactions between di16:0DAG, 16:0,18:1-phosphatidylcholine (POPC), and cholesterol. Despite of their substantial differences in chemical structure, DAG and cholesterol produce some very similar effects in POPC bilayers: increasing acyl chain order and bilayer thickness, reducing volume-per-lipid, and decreasing lateral diffusion of molecules. More significantly, DAG also produces a strong "condensing effect" in PC bilayers. In comparison, cholesterol is more effective than DAG in producing the above effects. The driving force for the condensing effect is their molecular shape: DAG and cholesterol both have small polar headgroups and large hydrophobic bodies. In a lipid bilayer, in order to avoid the unfavorable exposure of their hydrophobic parts to water, neighboring phospholipid headgroups move toward cholesterol or DAG to provide cover. Thus, seemingly complex interactions between DAG, cholesterol and phospholipid can be clearly explained using the Umbrella Model. Our simulations confirmed the hypothesis that DAG increases the spacing between phospholipid headgroups, which is important for activating protein kinase C and other enzymes. Interestingly, our simulations also show that the conventional wisdom that the spacing created by a DAG is directly above the DAG molecule is incorrect; instead, the largest spacing usually occurs between the first and the second nearest-neighbor PC headgroups from a DAG, due to the umbrella effect.

Nonintercalating nanosubstrates create asymmetry between bilayer leaflets

The physical properties of lipid bilayers can be remodeled by a variety of environmental factors. Here we investigate using molecular dynamics simulations the specific effects of nanoscopic substrates or external contact points on lipid membranes. We expose palmitoyl-oleoyl phosphatidylcholine bilayers unilaterally and separately to various model nanosized substrates differing in surface hydroxyl densities. We find that a surface hydroxyl density as low as 10% is sufficient to keep the bilayer juxtaposed to the substrate. The bilayer interacts with the substrate indirectly through multiple layers of water molecules; however, despite such buffered interaction, the bilayers exhibit certain properties different from unsupported bilayers. The substrates modify transverse lipid fluctuations, charge density profiles, and lipid diffusion rates, although differently in the two leaflets, which creates an asymmetry between bilayer leaflets. Other properties that include lipid cross-sectional areas, component volumes, and order parameters are minimally affected. The extent of asymmetry that we observe between bilayer leaflets is well beyond what has been reported for bilayers adsorbed on infinite solid supports. This is perhaps because the bilayers are much closer to our nanosized finite supports than to infinite solid supports, resulting in a stronger support-bilayer electrostatic coupling. The exposure of membranes to nanoscopic contact points, therefore, cannot be considered as a simple linear interpolation between unsupported membranes and membranes supported on infinite supports. In the biological context, this suggests that the exposure of membranes to nonintercalating proteins, such as those belonging to the cytoskeleton, should not always be considered as passive nonconsequential interactions.

Hydrophobic mismatch of mobile transmembrane helices: Merging theory and experiments

Hydrophobic mismatch still represents a puzzle for transmembrane peptides, despite the apparent simplicity of this concept and its demonstrated validity in natural membranes. Using a wealth of available experimental ((2))H NMR data, we provide here a comprehensive explanation of the orientation and dynamics of model peptides in lipid bilayers, which shows how they can adapt to membranes of different thickness. The orientational adjustment of transmembrane α-helices can be understood as the result of a competition between the thermodynamically unfavorable lipid repacking associated with peptide tilting and the optimization of peptide/membrane hydrophobic coupling. In the positive mismatch regime (long-peptide/thin-membrane) the helices adapt mainly via changing their tilt angle, as expected from simple geometrical predictions. However, the adaptation mechanism varies with the peptide sequence in the flanking regions, suggesting additional effects that modulate hydrophobic coupling. These originate from re-adjustments of the peptide hydrophobic length and they depend on the hydrophobicity of the flanking region, the strength of interfacial anchoring, the structural flexibility of anchoring side-chains and the presence of alternative anchoring residues.

Material properties of matrix lipids determine the conformation and intermolecular reactivity of diacetylenic phosphatidylcholine in the lipid bilayer

Photopolymerizable phospholipid DC(8,9)PC (1,2-bis-(tricosa-10,12-diynoyl)-sn-glycero-3-phosphocholine) exhibits unique assembly characteristics in the lipid bilayer. Because of the presence of the diacetylene groups, DC(8,9)PC undergoes polymerization upon UV (254 nm) exposure and assumes chromogenic properties. DC(8,9)PC photopolymerization in gel-phase matrix lipid 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) monitored by UV-vis absorption spectroscopy occurred within 2 min after UV treatment, whereas no spectral shifts were observed when DC(8,9)PC was incorporated into liquid-phase matrix 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC). Liquid chromatography-tandem mass spectrometry analysis showed a decrease in the amount of DC(8,9)PC monomer in both DPPC and POPC environments without any change in the matrix lipids in UV-treated samples. Molecular dynamics (MD) simulations of DPPC/DC(8,9)PC and POPC/DC(8,9)PC bilayers indicate that the DC(8,9)PC molecules adjust to the thickness of the matrix lipid bilayer. Furthermore, the motions of DC(8,9)PC in the gel-phase bilayer are more restricted than in the fluid bilayer. The restricted motional flexibility of DC(8,9)PC (in the gel phase) enables the reactive diacetylenes in individual molecules to align and undergo polymerization, whereas the unrestricted motions in the fluid bilayer restrict polymerization because of the lack of appropriate alignment of the DC(8,9)PC fatty acyl chains. Fluorescence microscopy data indicates the homogeneous distribution of lipid probe 1,2-dioleoyl-sn-glycero-3-phosphoethanolamine-N-lissamine rhodamine B sulfonyl ammonium salt (N-Rh-PE) in POPC/DC(8,9)PC monolayers but domain formation in DPPC/DC(8,9)PC monolayers. These results show that the DC(8,9)PC molecules cluster and assume the preferred conformation in the gel-phase matrix for the UV-triggered polymerization reaction.

Molecular dynamics simulations of membranes composed of glycolipids and phospholipids

Lipid membranes composed of 1,2-di-(9Z,12Z,15Z)-octadecatrienoyl-3-O-β-D-galactosyl-sn-glycerol or monogalactosyldiacylglycerol (MGDG) and 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) were studied by means of molecular dynamics (MD) computer simulations. Three lipid compositions were considered: 0%, 20%, and 45% MGDG (by mole) denoted as MG-0, MG-20, and MG-45, respectively. The article is focused on the calculation of NMR dipolar interactions, which were confronted with previously reported experimental couplings. Dynamical processes and orientational distributions relevant for the averaging of dipolar interactions were evaluated. Furthermore, several parameters important for characterization of the bilayer structure, molecular organization, and dynamics were investigated. In general, only a minor change in DMPC properties was observed upon the increased MGDG/DMPC ratio, whereas properties related to MGDG undergo a more pronounced change. This effect was ascribed to the fact that DMPC is a bilayer (L(α)) forming lipid, whereas MGDG prefers a reverse hexagonal (H(II)) arrangement.

Membrane-lipid therapy in operation: the HSP co-inducer BGP-15 activates stress signal transduction pathways by remodeling plasma membrane rafts

Aging and pathophysiological conditions are linked to membrane changes which modulate membrane-controlled molecular switches, causing dysregulated heat shock protein (HSP) expression. HSP co-inducer hydroxylamines such as BGP-15 provide advanced therapeutic candidates for many diseases since they preferentially affect stressed cells and are unlikely have major side effects. In the present study in vitro molecular dynamic simulation, experiments with lipid monolayers and in vivo ultrasensitive fluorescence microscopy showed that BGP-15 alters the organization of cholesterol-rich membrane domains. Imaging of nanoscopic long-lived platforms using the raft marker glycosylphosphatidylinositol-anchored monomeric green fluorescent protein diffusing in the live Chinese hamster ovary (CHO) cell plasma membrane demonstrated that BGP-15 prevents the transient structural disintegration of rafts induced by fever-type heat stress. Moreover, BGP-15 was able to remodel cholesterol-enriched lipid platforms reminiscent of those observed earlier following non-lethal heat priming or membrane stress, and were shown to be obligate for the generation and transmission of stress signals. BGP-15 activation of HSP expression in B16-F10 mouse melanoma cells involves the Rac1 signaling cascade in accordance with the previous observation that cholesterol affects the targeting of Rac1 to membranes. Finally, in a human embryonic kidney cell line we demonstrate that BGP-15 is able to inhibit the rapid heat shock factor 1 (HSF1) acetylation monitored during the early phase of heat stress, thereby promoting a prolonged duration of HSF1 binding to heat shock elements. Taken together, our results indicate that BGP-15 has the potential to become a new class of pharmaceuticals for use in ‘membrane-lipid therapy’ to combat many various protein-misfolding diseases associated with aging.

Interaction of melittin peptides with perfluorocarbon nanoemulsion particles

Melittin, an antimicrobial peptide, forms pores in biological membranes and triggers cell death. Therefore, it has potential as an anticancer therapy. However, until recently, the therapeutic application of melittin has been impractical because a suitable platform for delivery was not available. Recently, we showed that phospholipid-stabilized perfluorooctyl bromide based nanoemulsion particles (PFOB-NEPs) were resistant to destruction by melittin and enabled specific delivery of melittin to tumor cells, killing them and reducing tumor growth. Earlier, prior work also showed that melittin adsorbed onto the stabilizing phospholipid monolayer of PFOB-NEP but did not disrupt the phospholipid monolayer or produce "cracking" of the PFOB-NEPs. The present work identifies the important structural motifs for melittin binding to PFOB-NEPs through a series of atomistic molecular dynamics simulations. The conformational ensemble of melittin bound to PFOB-NEP lipid monolayer was compared to structure from a control simulation of melittin bound to a lipid bilayer to identify several differences in melittin-lipid interactions between the two systems. First, melittin was deeply buried in the hydrophobic tail region of bilayer, while its depth was attenuated in the PFOB-NEP monolayer. Second, a helical conformation was the major secondary structure in the bilayer, but the fraction of helix was reduced in the PFOB-NEP. Finally, the overall pattern for the direct interaction of melittin with surrounding lipids was similar between liposome and PFOB-NEP, but the level of interaction was slightly decreased in the PFOB-NEP. These results suggest that melittin interacts with the monolayer of PFOB-NEP in a way that is similar way to its interaction with bilayers but that deeper penetration into the hydrophobic interior is inhibited.

Revealing the functionality of hypothetical protein KPN00728 from Klebsiella pneumoniae MGH78578: molecular dynamics simulation approaches

Background

Previously, the hypothetical protein, KPN00728 from Klebsiella pneumoniae MGH78578 was the Succinate dehydrogenase (SDH) chain C subunit via structural prediction and molecular docking simulation studies. However, due to limitation in docking simulation, an in-depth understanding of how SDH interaction occurs across the transmembrane of mitochondria could not be provided.

Results

In this present study, molecular dynamics (MD) simulation of KPN00728 and SDH chain D in a membrane was performed in order to gain a deeper insight into its molecular role as SDH. Structural stability was successfully obtained in the calculation for area per lipid, tail order parameter, thickness of lipid and secondary structural properties. Interestingly, water molecules were found to be highly possible in mediating the interaction between Ubiquinone (UQ) and SDH chain C via interaction with Ser27 and Arg31 residues as compared with earlier docking study. Polar residues such as Asp95 and Glu101 (KPN00728), Asp15 and Glu78 (SDH chain D) might have contributed in the creation of a polar environment which is essential for electron transport chain in Krebs cycle.

Conclusions

As a conclusion, a part from the structural stability comparability, the dynamic of the interacting residues and hydrogen bonding analysis had further proved that the interaction of KPN00728 as SDH is preserved and well agreed with our postulation earlier.

The role of membrane thickness in charged protein-lipid interactions

Charged amino acids are known to be important in controlling the actions of integral and peripheral membrane proteins and cell disrupting peptides. Atomistic molecular dynamics studies have shed much light on the mechanisms of membrane binding and translocation of charged protein groups, yet the impact of the full diversity of membrane physico-chemical properties and topologies has yet to be explored. Here we have performed a systematic study of an arginine (Arg) side chain analog moving across saturated phosphatidylcholine (PC) bilayers of variable hydrocarbon tail length from 10 to 18 carbons. For all bilayers we observe similar ion-induced defects, where Arg draws water molecules and lipid head groups into the bilayers to avoid large dehydration energy costs. The free energy profiles all exhibit sharp climbs with increasing penetration into the hydrocarbon core, with predictable shifts between bilayers of different thickness, leading to barrier reduction from 26 kcal/mol for 18 carbons to 6 kcal/mol for 10 carbons. For lipids of 10 and 12 carbons we observe narrow transmembrane pores and corresponding plateaus in the free energy profiles. Allowing for movements of the protein and side chain snorkeling, we argue that the energetic cost for burying Arg inside a thin bilayer will be small, consistent with recent experiments, also leading to a dramatic reduction in pK(a) shifts for Arg. We provide evidence that Arg translocation occurs via an ion-induced defect mechanism, except in thick bilayers (of at least 18 carbons) where solubility-diffusion becomes energetically favored. Our findings shed light on the mechanisms of ion movement through membranes of varying composition, with implications for a range of charged protein-lipid interactions and the actions of cell-perturbing peptides. This article is part of a Special Issue entitled: Membrane protein structure and function.

Molecular origin of electron paramagnetic resonance line shapes on β-barrel membrane proteins: the local solvation environment modulates spin-label configuration

In this work, electron paramagnetic resonance (EPR) spectroscopy and X-ray crystallography were used to examine the origins of EPR line shapes from spin-labels at the protein-lipid interface on the β-barrel membrane protein BtuB. Two atomic-resolution structures were obtained for the methanethiosulfonate spin-label derivatized to cysteines on the membrane-facing surface of BtuB. At one of these sites, position 156, the label side chain resides in a pocket formed by neighboring residues; however, it extends from the protein surface and yields a single-component EPR spectrum in the crystal that results primarily from fast rotation about the fourth and fifth bonds linking the spin-label to the protein backbone. In lipid bilayers, site 156 yields a multicomponent spectrum resulting from different rotameric states of the labeled side chain. Moreover, changes in the lipid environment, such as variations in bilayer thickness, modulate the EPR spectrum by modulating label rotamer populations. At a second site, position 371, the labeled side chain interacts with a pocket on the protein surface, leading to a highly immobilized single-component EPR spectrum that is not sensitive to hydrocarbon thickness. This spectrum is similar to that seen at other sites that are deep in the hydrocarbon, such as position 170. This work indicates that the rotameric states of spin-labels on exposed hydrocarbon sites are sensitive to the environment at the protein-hydrocarbon interface, and that this environment may modulate weak interactions between the labeled side chain and the protein surface. In the case of BtuB, lipid acyl chain packing is not symmetric around the β-barrel, and EPR spectra from labeled hydrocarbon-facing sites in BtuB may reflect this asymmetry. In addition to facilitating the interpretation of EPR spectra of membrane proteins, these results have important implications for the use of long-range distance restraints in protein structure refinement that are obtained from spin-labels.

Interaction of a two-transmembrane-helix peptide with lipid bilayers and dodecyl sulfate micelles

To probe structural changes that occur when a membrane protein is transferred from lipid bilayers to SDS micelles, a fragment of bacteriorhodopsin containing transmembrane helical segments A and B was studied by fluorescence spectroscopy, molecular dynamics (MD) simulation, and stopped flow kinetics. In lipid bilayers, Förster resonance energy transfer (FRET) was observed between tyrosine 57 on helix B and tryptophans 10 and 12 on helix A. FRET efficiency decreased substantially when the peptide was transferred to SDS. MD simulation showed no evidence for significant disruption of helix-helix interactions in SDS micelles. However, a cluster of water molecules was observed to form a hydrogen-bonded network with the phenolic hydroxyl group of tyrosine 57, which probably causes the disappearance of tyrosine-to-tryptophan FRET in SDS. The tryptophan quantum yield decreased in SDS, and the change occurred at nearly the same rate as membrane solubilization. The results provide a clear example of the importance of corroborating distance changes inferred from FRET by using complementary methods.

Competing interactions for antimicrobial selectivity based on charge complementarity

Antimicrobial peptides (AMPs) are an evolutionary conserved component of the innate immune system and possible templates for the development of new antibiotics. An important property of antimicrobial peptides is their ability to discriminate bacterial from eucaryotic cells which is attributed to the difference in lipid composition of the outer leaflet of the plasma membrane between the two types of cells. Whereas eucaryotic cells usually expose zwitterionic lipids, procaryotic cells expose also anionic lipids which bind the cationic antimicrobial peptides electrostatically. An example is the antimicrobial peptide NK-2 which is highly cationic and favors binding to anionic membranes. In the present study, the difference in binding affinity of NK-2 for palmitoyl-oleoyl-phosphatidyl-glycerol (POPG) and palmitoyl-oleoyl-phosphatidyl-choline (POPC) is studied using molecular dynamics simulations in conjunction with a coarse grained model and thermodynamic integration, by computing the change in free energy and its components upon the transfer of NK-2 from POPC to POPG. The transfer is indeed found to be highly favorable. Interestingly, the favorable contribution from the electrostatic interaction between the peptide and the anionic lipids is overcompensated by an unfavorable contribution from the change in lipid-cation interactions due to the release of counterions from the lipids. The increase in entropy due to the release of the cations is compensated by other entropic components. The largest favorable contribution arises from the solvation of the counterions. Overall the interaction between NK-2 and POPG is not determined by a single driving force but a subtle balance of competing interactions.

Enhanced Lipid Diffusion and Mixing in Accelerated Molecular Dynamics

Accelerated molecular dynamics (aMD) is an enhanced sampling technique that expedites conformational space sampling by reducing the barriers separating various low-energy states of a system. Here, we present the first application of the aMD method on lipid membranes. Altogether, ∼1.5 μs simulations were performed on three systems: a pure POPC bilayer, a pure DMPC bilayer, and a mixed POPC:DMPC bilayer. Overall, the aMD simulations are found to produce significant speedup in trans–gauche isomerization and lipid lateral diffusion versus those in conventional MD (cMD) simulations. Further comparison of a 70-ns aMD run and a 300-ns cMD run of the mixed POPC:DMPC bilayer shows that the two simulations yield similar lipid mixing behaviors, with aMD generating a 2–3-fold speedup compared to cMD. Our results demonstrate that the aMD method is an efficient approach for the study of bilayer structural and dynamic properties. On the basis of simulations of the three bilayer systems, we also discuss the impact of aMD parameters on various lipid properties, which can be used as a guideline for future aMD simulations of membrane systems.

Fluid phase lipid areas and bilayer thicknesses of commonly used phosphatidylcholines as a function of temperature

The structural parameters of fluid phase bilayers composed of phosphatidylcholines with fully saturated, mixed, and branched fatty acid chains, at several temperatures, have been determined by simultaneously analyzing small-angle neutron and X-ray scattering data. Bilayer parameters, such as area per lipid and overall bilayer thickness have been obtained in conjunction with intrabilayer structural parameters (e.g. hydrocarbon region thickness). The results have allowed us to assess the effect of temperature and hydrocarbon chain composition on bilayer structure. For example, we found that for all lipids there is, not surprisingly, an increase in fatty acid chain trans-gauche isomerization with increasing temperature. Moreover, this increase in trans-gauche isomerization scales with fatty acid chain length in mixed chain lipids. However, in the case of lipids with saturated fatty acid chains, trans-gauche isomerization is increasingly tempered by attractive chain-chain van der Waals interactions with increasing chain length. Finally, our results confirm a strong dependence of lipid chain dynamics as a function of double bond position along fatty acid chains.

Identification of the GPR55 agonist binding site using a novel set of high-potency GPR55 selective ligands

Marijuana is the most widely abused illegal drug, and its spectrum of effects suggests that several receptors are responsible for the activity. Two cannabinoid receptor subtypes, CB1 and CB2, have been identified, but the complex pharmacological properties of exogenous cannabinoids and endocannabinoids are not fully explained by their signaling. The orphan receptor GPR55 binds a subset of CB1 and CB2 ligands and has been proposed as a cannabinoid receptor. This designation, however, is controversial as a result of recent studies in which lysophosphatidylinositol (LPI) was identified as a GPR55 agonist. Defining a biological role for GPR55 requires GPR55 selective ligands that have been unavailable. From a β-arrestin, high-throughput, high-content screen of 300000 compounds run in collaboration with the Molecular Libraries Probe Production Centers Network initiative (PubChem AID1965), we identified potent GPR55 selective agonists. By modeling of the GPR55 activated state, we compared the GPR55 binding conformations of three of the novel agonists obtained from the screen, CID1792197, CID1172084, and CID2440433 (PubChem Compound IDs), with that of LPI. Our modeling indicates the molecular shapes and electrostatic potential distributions of these agonists mimic those of LPI; the GPR55 binding site accommodates ligands that have inverted-L or T shapes with long, thin profiles that can fit vertically deep in the receptor binding pocket while their broad head regions occupy a horizontal binding pocket near the GPR55 extracellular loops. Our results will allow the optimization and design of second-generation GPR55 ligands and provide a means for distinguishing GPR55 selective ligands from those interacting with cannabinoid receptors.

Lipid-controlled peptide topology and interactions in bilayers: structural insights into the synergistic enhancement of the antimicrobial activities of PGLa and magainin 2

To gain further insight into the antimicrobial activities of cationic linear peptides, we investigated the topology of each of two peptides, PGLa and magainin 2, in oriented phospholipid bilayers in the presence and absence of the other peptide and as a function of the membrane lipid composition. Whereas proton-decoupled (15)N solid-state NMR spectroscopy indicates that magainin 2 exhibits stable in-plane alignments under all conditions investigated, PGLa adopts a number of different membrane topologies with considerable variations in tilt angle. Hydrophobic thickness is an important parameter that modulates the alignment of PGLa. In equimolar mixtures of PGLa and magainin 2, the former adopts transmembrane orientations in dimyristoyl-, but not 1-palmitoyl-2-oleoyl-, phospholipid bilayers, whereas magainin 2 remains associated with the surface in all cases. These results have important consequences for the mechanistic models explaining synergistic activities of the peptide mixtures and will be discussed. The ensemble of data suggests that the thinning of the dimyristoyl membranes caused by magainin 2 tips the topological equilibrium of PGLa toward a membrane-inserted configuration. Therefore, lipid-mediated interactions play a fundamental role in determining the topology of membrane peptides and proteins and thereby, possibly, in regulating their activities as well.

Comparative model studies of gastric toxicity of nonsteroidal anti-inflammatory drugs

A high percentage of people treated with a long-term nonsteroidal anti-inflammatory drug (NSAID) therapy suffer NSAID-induced gastrointestinal-tract-related side effects. A current hypothesis states that the side effects are related to the topical action of NSAID molecules on gastric mucus that lowers its resistance to luminal acid. The main lipids in human mucus are palmitoyloleoylphosphatidylcholine (POPC) and cholesterol (Chol). In this study, both X-ray diffraction and molecular dynamics (MD) simulation methods were employed to investigate the effects of selected NSAIDs in protonated and deprotonated states on the structural parameters of a POPC-Chol bilayer. The drugs were three commonly used NSAIDs with apparently different gastric toxicity: ketoprofen (KET), aspirin (ASP), and piroxicam (PXM). Both methods revealed that the effects of the NSAIDs on the POPC-Chol bilayer parameters were moderate and only slightly differentiated among the drugs. Much larger differences among the drugs were noticed in their interactions with interfacial water and Na(+) as well as with the polar groups of POPC and Chol, mainly via H-bonds. Of the three NSAIDs, KET interacted with POPC and water the most extensively, whereas ASP interacted with Chol and Na(+) more than did the other two. Interactions of PXM with POPC and Chol polar groups as well as with water and Na(+) were limited.

Membrane thickness varies around the circumference of the transmembrane protein BtuB

BtuB is a large outer-membrane β-barrel protein that belongs to a class of active transport proteins that are TonB-dependent. These TonB-dependent transporters are based upon a 22-stranded antiparallel β-barrel, which is notably asymmetric in its length. Here, site-directed spin labeling and simulated annealing were used to locate the membrane lipid interface surrounding BtuB when reconstituted into phosphatidylcholine bilayers. Positions on the outer facing surface of the β-barrel and the periplasmic turns were spin-labeled and distances from the label to the membrane interface estimated by progressive power saturation of the electron paramagnetic resonance spectra. These distances were then used as atom-to-plane distance restraints in a simulated annealing routine, to dock the protein to two independent planes and produce a model representing the average position of the lipid phosphorus atoms at each interface. The model is in good agreement with the experimental data; however, BtuB is mismatched to the bilayer thickness and the resulting planes representing the bilayer interface are not parallel. In the model, the membrane thickness varies by 11 Å around the circumference of the protein, indicating that BtuB distorts the bilayer interface so that it is thinnest on the short side of the protein β-barrel.

Conjugated double bonds in lipid bilayers: a molecular dynamics simulation study

Conjugated linoleic acids (CLA) are found naturally in dairy products. Two isomers of CLA, that differ only in the location of cis and trans double bonds, are found to have distinct and different biological effects. The cis 9 trans 11 (C9T11) isomer is attributed to have the anti-carcinogenic effects, while the trans 10 cis 12 (T10C12) isomer is believed to be responsible for the anti-obesity effects. Since dietary CLA are incorporated into membrane phospholipids, we have used Molecular Dynamics (MD) simulations to investigate the comparative effects of the two isomers on lipid bilayer structure. Specifically, simulations of phosphatidylcholine lipid bilayers in which the sn-2 chains contained one of the two isomers of CLA were performed. Force field parameters for the torsional potential of double bonds were obtained from ab initio calculations. From the MD trajectories we calculated and compared structural properties of the two lipid bilayers, including areas per molecule, density profiles, thickness of bilayers, tilt angle of tail chains, order parameters profiles, radial distribution function (RDF) and lateral pressure profiles. The main differences found between bilayers of the two CLA isomers, are (1) the order parameter profile for C9T11 has a dip in the middle of sn-2 chain while the profile for T10C12 has a deeper dip close to terminal of sn-2 chain, and (2) the lateral pressure profiles show differences between the two isomers. Our simulation results reveal localized physical structural differences between bilayers of the two CLA isomers that may contribute to different biological effects through differential interactions with membrane proteins or cholesterol.

Structural models of TREK channels and their gating mechanism

Mechanosensitive TREK channels belong to the family of K2P channels, a family of widely distributed, well modulated channels that uniquely have two similar or identical subunits, each with two TM1-P-TM2 motifs. Our goal is to build viable structural models of TREK channels, as representatives of K2P channels family. The structures available to be used as templates belong to the 2TM channels superfamily. These have low sequence similarity and different structural features: four symmetrically arranged subunits, each having one TM1-P-TM2 motif. Our model building strategy used two subunits of the template (KcsA) to build one subunit of the target (TREK-1). Our models of the Closed channel were adjusted to differ substantially from those of the template, e.g., TM2 of the 2nd repeat is near the axis of the pore whereas TM2 of the 1st repeat is far from the axis. Segments linking the two repeats and immediately following the last TM segment were modeled ab initio as α-helices based on helical periodicities of hydrophobic and hydrophilic residues, highly conserved and poorly conserved residues, and statistically related positions from multiple sequence alignments. The models were further refined by two-fold symmetry-constrained MD simulations using a protocol we developed previously. We also built models of the Open state and suggest a possible tension-activated gating mechanism characterized by helical motion with two-fold symmetry. Our models are consistent with deletion/truncation mutagenesis and thermodynamic analysis of gating described in the accompanying paper.

Raman Spectroscopic Imaging of Cholesterol and Docosahexaenoic Acid Distribution in the Retinal Rod Outer Segment

Raman vibrational spectroscopic imaging was performed on retinal rod cells isolated from bullfrogs (Rana catesbeiana). The Raman spectra enable determination of the lipid and protein rich rod outer segment (ROS) from the nucleus and inner segment of the cell. Peak fitting analysis of spectra obtained from individual rod photoreceptor cells show characteristic vibrational modes that can be associated with cholesterol and docosahexaenoic acid containing lipids. These results provide direct observations of biomolecular gradients in the rod photoreceptor cells, which, thus far, have been based on indirect detergent extracts and histochemical analysis with indicators such as filipin. The detected biomolecules are associated with regulation of the integral membrane protein rhodopsin, and methods capable direct observation of these biomolecules offer new routes to exploring their role in the regulation of cellular processes.

Combination of the CHARMM27 force field with united-atom lipid force fields

Computer simulations offer a valuable way to study membrane systems, from simple lipid bilayers to large transmembrane protein complexes and lipid-nucleic acid complexes for drug delivery. Their accuracy depends on the quality of the force field parameters used to describe the components of a particular system. We have implemented the widely used CHARMM22 and CHARMM27 force fields in the GROMACS simulation package to (i) combine the CHARMM22 protein force field with two sets of united-atom lipids parameters; (ii) allow comparisons of the lipid CHARMM27 force field with other lipid force fields or lipid-protein force field combinations. Our tests do not show any particular issue with the combination of the all-atom CHARMM22 force field with united-atoms lipid parameters, although pertinent experimental data are lacking to assess the quality of the lipid-protein interactions. The conversion utilities allow automatic generation of GROMACS simulation files with CHARMM nucleic acids and protein parameters and topologies, starting from pdb files using the standard GROMACS pdb2gmx method. CMAP is currently not implemented.

Lipid bilayer molecular dynamics study of lipid-derived agonists of the putative cannabinoid receptor, GPR55

Both L-α-lysophosphatidylinositol (LPI) and 2-arachidonoyl-sn-glycero-3-phosphoinositol (2-AGPI) have been reported to activate the putative cannabinoid receptor, GPR55. Recent microsecond time-scale molecular dynamics (MD) simulations and isothiocyanate covalent labeling studies have suggested that a transmembrane helix 6/7 (TMH6/7) lipid pathway for ligand entry may be necessary for interaction with cannabinoid receptors. Because LPI and 2-AGPI are lipid-derived ligands, conformations that each assumes in the lipid bilayer are therefore likely important for their interaction with GPR55. We report here the results of 70 ns NAMD molecular dynamics (MD) simulations of LPI and of 2-AGPI in a fully hydrated bilayer of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC). These simulations are compared with a 70 ns simulation of the cannabinoid CB1 receptor endogenous ligand, N-arachidonoylethanolamine (anandamide, AEA) in a POPC bilayer. These simulations revealed that (1) LPI and 2-AGPI sit much higher in the bilayer than AEA, with inositol headgroups that can at times be solvated completely by water; (2) the behavior of the acyl chains of AEA and 2-AGPI are similar in their flexibilities in the bilayer, while the acyl chain of LPI has reduced flexibility; and (3) both 2-AGPI and LPI can adopt a tilted headgroup orientation by hydrogen bonding to the phospholipid phosphate/glycerol groups or via intramolecular hydrogen bonding. This tilted head group conformation (which represents over 40% of the conformer population of LPI (42.2 ± 3.3%) and 2-AGPI (43.7 ± 1.4%)) may provide a low enough profile in the lipid bilayer for LPI and 2-AGPI to enter GPR55 via the putative TMH6/7 entry port.

Intramolecular hydrogen bonding in articaine can be related to superior bone tissue penetration: a molecular dynamics study

Local anesthetics (LAs) are drugs that cause reversible loss of nociception during surgical procedures. Articaine is a commonly used LA in dentistry that has proven to be exceptionally effective in penetrating bone tissue and induce anesthesia on posterior teeth in maxilla and mandibula. In the present study, our aim was to gain a deeper understanding of the penetration of articaine through biological membranes by studying the interactions of articaine with a phospholipid membrane. Our approach involves Langmuir monolayer experiments combined with molecular dynamics simulations. Membrane permeability of LAs can be modulated by pH due to a titratable amine group with a pKa value close to physiological pH. A change in protonation state is thus known to act as a lipophilicity switch in LAs. Our study shows that articaine has an additional unique lipophilicity switch in its ability to form an intramolecular hydrogen bond. We suggest this intramolecular hydrogen bond as a novel and additional solvent-dependent mechanism for modulation of lipophilicity of articaine which may enhance its diffusion through membranes and connective tissue.

Characterization of perfluorooctylbromide-based nanoemulsion particles using atomistic molecular dynamics simulations

Perfluorocarbon-based nanoemulsion particles have arisen as promising platforms for the cellular delivery of imaging and therapeutic agents to specific targets. However, current knowledge of the agent delivery mechanism is limited to qualitative and phenomenological models. Lack of detail at the molecular level has hence delayed optimizing or customizing nanoemulsion particles for therapeutic and imaging applications. Here we report the first atomistic structural details of a perfluorooctylbromide-based (PFOB-based) nanoemulsion particle (NEP) with a 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) lipid emulsifier. Newly developed PFOB force-field parameters were used in molecular dynamics simulations to model the PFOB-NEP interface in a planar configuration. These PFOB force field parameters were developed and tested to reproduce the characteristics of bulk PFOB as well as PFOB at interfaces with water and emulsifying phospholipids. The modeled PFOB-NEP interface demonstrated significant intercalation of PFOB into the emulsifying lipid monolayer and consequent changes in the structural, electrostatic, and mechanical properties of the POPC monolayer and PFOB. This intercalation provides an explanation for experimental data demonstrating melittin tryptophan fluorescence quenching upon binding to the nanoemulsion particles through the observation of direct contact between the melittin tryptophan and the PFOB bromine. Additionally, the atomistic details of the PFOB-NEP interface structure provided by our simulations are used to suggest the influence of each component on PFOB-NEP delivery function which will be tested in future coarse-grained simulations.

Elliptical structure of phospholipid bilayer nanodiscs encapsulated by scaffold proteins: casting the roles of the lipids and the protein

Phospholipid bilayers host and support the function of membrane proteins and may be stabilized in disc-like nanostructures, allowing for unprecedented solution studies of the assembly, structure, and function of membrane proteins (Bayburt et al. Nano Lett. 2002, 2, 853-856). Based on small-angle neutron scattering in combination with variable-temperature studies of synchrotron small-angle X-ray scattering on nanodiscs in solution, we show that the fundamental nanodisc unit, consisting of a lipid bilayer surrounded by amphiphilic scaffold proteins, possesses intrinsically an elliptical shape. The temperature dependence of the curvature of the nanodiscs prepared with two different phospholipid types (DLPC and POPC) shows that it is the scaffold protein that determines the overall elliptical shape and that the nanodiscs become more circular with increasing temperature. Our data also show that the hydrophobic bilayer thickness is, to a large extent, dictated by the scaffolding protein and adjusted to minimize the hydrophobic mismatch between protein and phospholipid. Our conclusions result from a new comprehensive and molecular-based model of the nanodisc structure and the use of this to analyze the experimental scattering profile from nanodiscs. The model paves the way for future detailed structural studies of functional membrane proteins encapsulated in nanodiscs.

Lateral diffusion of membrane proteins: consequences of hydrophobic mismatch and lipid composition

Biological membranes are composed of a large number lipid species differing in hydrophobic length, degree of saturation, and charge and size of the headgroup. We now present data on the effect of hydrocarbon chain length of the lipids and headgroup composition on the lateral mobility of the proteins in model membranes. The trimeric glutamate transporter (GltT) and the monomeric lactose transporter (LacY) were reconstituted in giant unilamellar vesicles composed of unsaturated phosphocholine lipids of varying acyl chain length (14-22 carbon atoms) and various ratios of DOPE/DOPG/DOPC lipids. The lateral mobility of the proteins and of a fluorescent lipid analog was determined as a function of the hydrophobic thickness of the bilayer (h) and lipid composition, using fluorescence correlation spectroscopy. The diffusion coefficient of LacY decreased with increasing thickness of the bilayer, in accordance with the continuum hydrodynamic model of Saffman-Delbrück. For GltT, the mobility had its maximum at diC18:1 PC, which is close to the hydrophobic thickness of the bilayer in vivo. The lateral mobility decreased linearly with the concentration of DOPE but was not affected by the fraction of anionic lipids from DOPG. The addition of DOPG and DOPE did not affect the activity of GltT. We conclude that the hydrophobic thickness of the bilayer is a major determinant of molecule diffusion in membranes, but protein-specific properties may lead to deviations from the Saffman-Delbrück model.

A systematically coarse-grained solvent-free model for quantitative phospholipid bilayer simulations

We present an implicit solvent coarse-grained (CG) model for quantitative simulations of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayers. The absence of explicit solvent enables membrane simulations on large length and time scales at moderate computational expense. Despite improved computational efficiency, the model preserves chemical specificity and quantitative accuracy. The bonded and nonbonded interactions together with the effective cohesion mimicking the hydrophobic effect were systematically tuned by matching structural and mechanical properties from experiments and all-atom bilayer simulations, such as saturated area per lipid, radial distribution functions, density and pressure profiles across the bilayer, P(2) order, etc. The CG lipid model is shown to self-assemble into a bilayer starting from a random dispersion. Its line tension and elastic properties, such as bending and stretching modulus, are semiquantitatively consistent with experiments. The effects of (i) reduced molecular friction and (ii) more efficient integration combine to an overall speed-up of 3-4 orders of magnitude compared to all-atom bilayer simulations. Our CG lipid model is especially useful for studies of large-scale phenomena in membranes that nevertheless require a fair description of chemical specificity, e.g., membrane patches interacting with movable and transformable membrane proteins and peptides.

Systematic implicit solvent coarse-graining of bilayer membranes: lipid and phase transferability of the force field

We study the lipid and phase transferability of our recently developed systematically coarse-grained solvent-free membrane model. The force field was explicitly parameterized to describe a fluid 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) bilayer at 310 K with correct structure and area per lipid, while gaining at least three orders of magnitude in computational efficiency (see Wang and Deserno 2010 J. Phys. Chem. B 114 11207-20). Here, we show that exchanging CG tails, without any subsequent re-parameterization, creates reliable models of 1,2-dioleoylphosphatidylcholine (DOPC) and 1,2-dipalmitoylphosphatidylcholine (DPPC) lipids in terms of structure and area per lipid. Furthermore, all CG lipids undergo a liquid-gel transition upon cooling, with characteristics like those observed in experiments and all-atom simulations during phase transformation. These studies suggest a promising transferability of our force field parameters to different lipid species and thermodynamic state points, properties that are a prerequisite for even more complex systems, such as mixtures.

HIV fusion peptide penetrates, disorders, and softens T-cell membrane mimics

This work investigates the interaction of N-terminal gp41 fusion peptide (FP) of human immunodeficiency virus type 1 (HIV-1) with model membranes in order to elucidate how FP leads to fusion of HIV and T-cell membranes. FP constructs were (i) wild-type FP23 (23 N-terminal amino acids of gp41), (ii) water-soluble monomeric FP that adds six lysines on the C-terminus of FP23 (FPwsm), and (iii) the C-terminus covalently linked trimeric version (FPtri) of FPwsm. Model membranes were (i) LM3 (a T-cell mimic), (ii) 1,2-dioleoyl-sn-glycero-3-phosphocholine, (iii) 1,2-dioleoyl-sn-glycero-3-phosphocholine/30 mol% cholesterol, (iv) 1,2-dierucoyl-sn-glycero-3-phosphocholine, and (v) 1,2-dierucoyl-sn-glycero-3-phosphocholine/30 mol% cholesterol. Diffuse synchrotron low-angle x-ray scattering from fully hydrated samples, supplemented by volumetric data, showed that FP23 and FPtri penetrate into the hydrocarbon region and cause membranes to thin. Depth of penetration appears to depend upon a complex combination of factors including bilayer thickness, presence of cholesterol, and electrostatics. X-ray data showed an increase in curvature in hexagonal phase 1,2-dioleoyl-sn-glycero-3-phosphoethanolamine, which further indicates that FP23 penetrates into the hydrocarbon region rather than residing in the interfacial headgroup region. Low-angle x-ray scattering data also yielded the bending modulus K(C), a measure of membrane stiffness, and wide-angle x-ray scattering yielded the S(xray) orientational order parameter. Both FP23 and FPtri decreased K(C) and S(xray) considerably, while the weak effect of FPwsm suggests that it did not partition strongly into LM3 model membranes. Our results are consistent with the HIV FP disordering and softening the T-cell membrane, thereby lowering the activation energy for viral membrane fusion.

Electrically addressable, biologically relevant surface-supported bilayers

The assembly of electrically addressable, planar supported bilayers composed of biologically relevant lipids, such as those used in vesicular systems, will greatly enhance the experimental capabilities in membrane and membrane protein research. Here we assess the electrical properties of bilayers composed of a wide range of physiologically relevant lipids and lipid combinations. We demonstrate that robust, biologically relevant, planar supported bilayers with high resistance composed of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and 25 mol % cholesterol can be constructed with high reproducibility. Furthermore, to enable studies of pore-forming peptides, which are commonly cationic, we demonstrate the construction of bilayers with biologically relevant outer leaflets incorporating up to 10 mol % negatively charged lipids. Unique features of the platform are that (1) the substrate is commercially available, atomically smooth, single-crystal silicon, (2) the polymer cushion allows for the natural incorporation of membrane proteins, and (3) the platform is highly reproducible.

Water replacement hypothesis in atomic details: effect of trehalose on the structure of single dehydrated POPC bilayers

We present molecular dynamics (MD) simulations to study the plausibility of the water replacement hypothesis (WRH) from the viewpoint of structural chemistry. A total of 256 2-oleoyl-1-palmitoyl-sn-glycero-3-phosphocholine (POPC) lipids were modeled for 400 ns at 11.7 or 5.4 waters/lipid. To obtain a single dehydrated bilayer relevant to the WRH, simulations were performed in the NP(xy)h(z)T ensemble with h(z) > 8 nm, allowing interactions between lipids in the membrane plane and preventing interactions between neighboring membranes via periodic boundary conditions. This setup resulted in a stable single bilayer in (or near) the gel state. Trehalose caused a concentration-dependent increase of the area per lipid (APL) accompanied by fluidizing the bilayer core. This mechanism has been suggested by the WRH. However, dehydrated bilayers in the presence of trehalose were not structurally identical to fully hydrated bilayers. The headgroup vector was in a more parallel orientation in dehydrated bilayers with respect to the bilayer plane and maintained this orientation in the presence of trehalose in spite of APL increase. The total dipole potential changed sign in dehydrated bilayers and remained slightly positive in the presence of trehalose. The model of a dehydrated bilayer presented here allows the study of the mechanisms of membrane protection against desiccation by different compounds.