A Lipid Pathway for Ligand Binding Is Necessary for a Cannabinoid G Protein-coupled Receptor

Recent isothiocyanate covalent labeling studies have suggested that a classical cannabinoid, (−)-7′-isothiocyanato-11-hydroxy-1′,1′dimethylheptyl-hexahydrocannabinol (AM841), enters the cannabinoid CB2 receptor via the lipid bilayer (Pei, Y., Mercier, R. W., Anday, J. K., Thakur, G. A., Zvonok, A. M., Hurst, D., Reggio, P. H., Janero, D. R., and Makriyannis, A. (2008) Chem. Biol. 15, 1207–1219). However, the sequence of steps involved in such a lipid pathway entry has not yet been elucidated. Here, we test the hypothesis that the endogenous cannabinoid sn-2-arachidonoylglycerol (2-AG) attains access to the CB2 receptor via the lipid bilayer. To this end, we have employed microsecond time scale all-atom molecular dynamics (MD) simulations of the interaction of 2-AG with CB2 via a palmitoyl-oleoyl-phosphatidylcholine lipid bilayer. Results suggest the following: 1) 2-AG first partitions out of bulk lipid at the transmembrane α-helix (TMH) 6/7 interface; 2) 2-AG then enters the CB2 receptor binding pocket by passing between TMH6 and TMH7; 3) the entrance of the 2-AG headgroup into the CB2 binding pocket is sufficient to trigger breaking of the intracellular TMH3/6 ionic lock and the movement of the TMH6 intracellular end away from TMH3; and 4) subsequent to protonation at D3.49/D6.30, further 2-AG entry into the ligand binding pocket results in both a W6.48 toggle switch change and a large influx of water. To our knowledge, this is the first demonstration via unbiased molecular dynamics that a ligand can access the binding pocket of a class A G protein-coupled receptor via the lipid bilayer and the first demonstration via molecular dynamics of G protein-coupled receptor activation triggered by a ligand binding event.     

What is the natural ligand of GPR55?

GPR55 is a seven transmembrane G protein-coupled receptor and was originally identified as a putative third cannabinoid receptor. Recently, lysophosphatidylinositol (LPI) was reported to be a GPR55 ligand. Stimulation of GPR55 by LPI activates G(12/13) and G(q/11) proteins, induces phosphorylation of the extracellular signal-regulated kinase and increases intracellular calcium concentration. Lysophospholipids are molecularly quite diverse across species and tissues. A recent report showed that the predominant fatty acyl moiety of LPI in rat brain is stearic acid followed by arachidonic acid. The biological activity of arachidonic acid-containing LPI species towards GPR55 was shown to be markedly higher than that of LPI species containing other fatty acyl groups, suggesting that 2-arachidonolyl LPI is the most likely natural ligand of GPR55.     

Saturation with cholesterol increases vertical order and smoothes the surface of the phosphatidylcholine bilayer: a molecular simulation study

Molecular dynamics (MD) simulations of a mono-cis-unsaturated 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) bilayer and a POPC bilayer containing 50mol% cholesterol (POPC-Chol50) were carried out for 200ns to compare the spatial organizations of the pure POPC bilayer and the POPC bilayer saturated with Chol. The results presented here indicate that saturation with Chol significantly narrows the distribution of vertical positions of the center-of-mass of POPC molecules and POPC atoms in the bilayer. In the POPC-Chol50 bilayer, the same moieties of the lipid molecules are better aligned at a given bilayer depth, forming the following clearly separated membrane regions: the polar headgroup, the rigid core consisting of steroid rings and upper fragments of the acyl chains, and the fluid hydrocarbon core consisting of Chol chains and the lower fragments of POPC chains. The membrane surface of the POPC-Chol50 bilayer is smooth. The results have biological significance because the POPC-Chol50 bilayer models the bulk phospholipid portion of the fiber-cell membrane in the eye lens. It is hypothesized that in the eye lens cholesterol-induced smoothing of the membrane surface decreases light-scattering and helps to maintain lens transparency.     

Coupling Molecular Dynamics Simulations with Experiments for the Rational Design of Indolicidin-Analogous Antimicrobial Peptides

Antimicrobial peptides (AMPs) have attracted much interest in recent years because of their potential use as new-generation antibiotics. Indolicidin (IL) is a 13-residue cationic AMP that is effective against a broad spectrum of bacteria, fungi, and even viruses. Unfortunately, its high hemolytic activity retards its clinical applications. In this study, we adopted molecular dynamics (MD) simulations as an aid toward the rational design of IL analogues exhibiting high antimicrobial activity but low hemolysis. We employed long-timescale, multi-trajectory all-atom MD simulations to investigate the interactions of the peptide IL with model membranes. The lipid bilayer formed by the zwitterionic 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) was chosen as the model erythrocyte membrane; lipid bilayers formed from a mixture of POPC and the negatively charged 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol were chosen to model bacterial membranes. MD simulations with a total simulation time of up to 4 μs revealed the mechanisms of the processes of IL adsorption onto and insertion into the membranes. The packing order of these lipid bilayers presumably correlated to the membrane stability upon IL adsorption and insertion. We used the degree of local membrane thinning and the reduction in the order parameter of the acyl chains of the lipids to characterize the membrane stability. The order of the mixed 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol/POPC lipid bilayer reduced significantly upon the adsorption of IL. On the other hand, although the order of the pure-POPC lipid bilayer was perturbed slightly during the adsorption stage, the value was reduced more dramatically upon the insertion of IL into the membrane's hydrophobic region. The results imply that enhancing IL adsorption on the microbial membrane may amplify its antimicrobial activity, while the degree of hemolysis may be reduced through inhibition of IL insertion into the hydrophobic region of the erythrocyte membrane. In addition, through simulations, we identified the amino acids that are most responsible for the adsorption onto or insertion into the two model membranes. Positive charges are critical to the peptide's adsorption, whereas the presence of hydrophobic Trp8 and Trp9 leads to its deeper insertion. Combining the hypothetical relationships between the membrane disordering and the antimicrobial and hemolytical activities with the simulated results, we designed three new IL-analogous peptides: IL-K7 (Pro7 → Lys), IL-F89 (Trp8 and Trp9 → Phe), and IL-K7F89 (Pro7 → Lys; Trp8 and Trp9 → Phe). The hemolytic activity of IL-F89 is considerably lower than that of IL, whereas the antimicrobial activity of IL-K7 is greatly enhanced. In particular, the de novo peptide IL-K7F89 exhibits higher antimicrobial activity against Escherichia coli; its hemolytic activity decreased to only 10% of that of IL. Our simulated and experimental results correlated well. This approach—coupling MD simulations with experimental design—is a useful strategy toward the rational design of AMPs for potential therapeutic use.     

Lipid Gymnastics: Evidence of Complete Acyl Chain Reversal in Oxidized Phospholipids from Molecular Simulations

In oxidative environments, biomembranes contain oxidized lipids with short, polar acyl chains. Two stable lipid oxidation products are PoxnoPC and PazePC. PoxnoPC has a carbonyl group, and PazePC has an anionic carboxyl group pendant at the end of the short, oxidized acyl chain. We have used MD simulations to explore the possibility of complete chain reversal in OXPLs in POPC-OXPL mixtures. The polar AZ chain of PazePC undergoes chain reversal without compromising the lipid bilayer integrity at concentrations up to 25% OXPL, and the carboxyl group points into the aqueous phase. Counterintuitively, the perturbation of overall membrane structural and dynamic properties is stronger for PoxnoPC than for PazePC. This is because of the overall condensing and ordering effect of sodium ions bound strongly to the lipids in the PazePC simulations. The reorientation of AZ chain is similar for two different lipid force fields. This work provides the first molecular evidence of the “extended lipid conformation” in phospholipid membranes. The chain reversal of PazePC lipids decorates the membrane interface with reactive, negatively charged functional groups. Such chain reversal is likely to exert a profound influence on the structure and dynamics of biological membranes, and on membrane-associated biological processes.     

Partitioning of 2,6-Bis(1H-Benzimidazol-2-yl)pyridine fluorophore into a phospholipid bilayer: complementary use of fluorescence quenching studies and molecular dynamics simulations

Successful use of fluorescence sensing in elucidating the biophysical properties of lipid membranes requires knowledge of the distribution and location of an emitting molecule in the bilayer. We report here that 2,6-bis(1H-benzimidazol-2-yl)pyridine (BBP), which is almost non-fluorescent in aqueous solutions, reveals a strong emission enhancement in a hydrophobic environment of a phospholipid bilayer, making it interesting for fluorescence probing of water content in a lipid membrane. Comparing the fluorescence behavior of BBP in a wide variety of solvents with those in phospholipid vesicles, we suggest that the hydrogen bonding interactions between a BBP fluorophore and water molecules play a crucial role in the observed “light switch effect”. Therefore, the loss of water-induced fluorescence quenching inside a membrane are thought to be due to deep penetration of BBP into the hydrophobic, water-free region of a bilayer. Characterized by strong quenching by transition metal ions in solution, BBP also demonstrated significant shielding from the action of the quencher in the presence of phospholipid vesicles. We used the increase in fluorescence intensity, measured upon titration of probe molecules with lipid vesicles, to estimate the partition constant and the Gibbs free energy (ΔG) of transfer of BBP from aqueous buffer into a membrane. Partitioning BBP revealed strongly favorable ΔG, which depends only slightly on the lipid composition of a bilayer, varying in a range from − 6.5 to − 7.0 kcal/mol. To elucidate the binding interactions of the probe with a membrane on the molecular level, a distribution and favorable location of BBP in a POPC bilayer were modeled via atomistic molecular dynamics (MD) simulations using two different approaches: (i) free, diffusion-driven partitioning of the probe molecules into a bilayer and (ii) constrained umbrella sampling of a penetration profile of the dye molecule across a bilayer. Both of these MD approaches agreed with regard to the preferred location of a BBP fluorophore within the interfacial region of a bilayer, located between the hydrocarbon acyl tails and the initial portion of the lipid headgroups. MD simulations also revealed restricted permeability of water molecules into this region of a POPC bilayer, determining the strong fluorescence enhancement observed experimentally for the membrane-partitioned form of BBP.     

Effect of sodium chloride on a lipid bilayer

Electrostatic interactions govern structural and dynamical properties of membranes and can vary considerably with the composition of the aqueous buffer. We studied the influence of sodium chloride on a pure POPC lipid bilayer by fluorescence correlation spectroscopy experiments and molecular dynamics simulations. Increasing sodium chloride concentration was found to decrease the self-diffusion of POPC lipids within the bilayer. Self-diffusion coefficients calculated from the 100 ns simulations agree with those measured on a millisecond timescale, suggesting that most of the relaxation processes relevant for lipid diffusion are faster than the simulation timescale. As the dominant effect, the molecular dynamics simulations revealed a tight binding of sodium ions to the carbonyl oxygens of on average three lipids leading to larger complexes with reduced mobility. Additionally, the bilayer thickens by approximately 2 A, which increases the order parameter of the fatty acyl chains. Sodium binding alters the electrostatic potential, which is largely compensated by a changed polarization of the aqueous medium and a lipid dipole reorientation.     

Lipid hydration and mobility: an interplay between fluorescence solvent relaxation experiments and molecular dynamics simulations

Fluorescence solvent relaxation experiments are based on the characterization of time-dependent shifts in the fluorescence emission of a chromophore, yielding polarity and viscosity information about the chromophore’s immediate environment. A chromophore applied to a phospholipid bilayer at a well-defined location (with respect to the z-axis of the bilayer) allows monitoring of the hydration and mobility of the probed segment of the lipid molecules. Specifically, time-resolved fluorescence experiments, fluorescence quenching data and molecular dynamic (MD) simulations show that 6-lauroyl-2-dimethylaminonaphthalene (Laurdan) probes the hydration and mobility of the sn-1 acyl groups in a phosphatidylcholine bilayer. The time-dependent fluorescence shift (TDFS) of Laurdan provides information on headgroup compression and expansion induced by the addition of different amounts of cationic lipids to phosphatidylcholine bilayers. Those changes were predicted by previous MD simulations. Addition of truncated oxidized phospholipids leads to increased mobility and hydration at the sn-1 acyl level. This experimental finding can be explained by MD simulations, which indicate that the truncated chains of the oxidized lipid molecules are looping back into aqueous phase, hence creating voids below the glycerol level. Fluorescence solvent relaxation experiments are also useful in understanding salt effects on the structure and dynamics of lipid bilayers. For example, such experiments demonstrate that large anions increase hydration and mobility at the sn-1 acyl level of phosphatidylcholine bilayers, an observation which could not be explained by standard MD simulations. If polarizability is introduced into the applied force field, however, MD simulations show that big soft polarizable anions are able to interact with the hydrophilic/hydrophobic interface of the lipid bilayer, penetrating to the level probed by Laurdan, and that they expand and destabilize the bilayer making it more hydrated and mobile.     

Binding and insertion of alpha-helical anti-microbial peptides in POPC bilayers studied by molecular dynamics simulations

We have performed molecular dynamics simulations of the interactions of two alpha-helical anti-microbial peptides, magainin2 and its synthetic analog of MSI-78, with palmitoyl-oleoyl-phosphatidylcholine (POPC) lipid bilayers. We used various initial positions and orientations of the peptide with respect to the lipid bilayer, including a surface-bound state parallel to the interface, a trans-membrane state, and a partially inserted state. Our 20 ns long simulations show that both magainin2 and MSI-78 are most stable in the lipid environment, with the peptide destabilized to different extents in both aqueous and lipid/water interfacial environments. We found that there are strong specific interactions between the lysine residues of the peptides and the lipid head-group regions. MSI-78, owing to its large number of lysines, shows better binding characteristics and overall stability when compared to magainin2. We also find that both peptides destabilize the bilayer environment, as observed by the increase in lipid tail disorder and the induction of local curvature on the lipid head-groups by the peptides. From all the simulations, we conclude that the hydrogen bonding interactions between the lysines of the peptides and the oxygens of the polar lipid head-groups are the strongest and determine the overall peptide binding characteristics to the lipids.     

Molecular Dynamics Simulations of PIP2 and PIP3 in Lipid Bilayers: Determination of Ring Orientation, and the Effects of Surface Roughness on a Poisson-Boltzmann Description

Molecular dynamics (MD) simulations of phosphatidylinositol (4,5)-bisphosphate (PIP₂) and phosphatidylinositol (3,4,5)-trisphosphate (PIP₃) in 1-palmitoyl 2-oleoyl phosphatidylcholine (POPC) bilayers indicate that the inositol rings are tilted ∼40° with respect to the bilayer surface, as compared with 17° for the P-N vector of POPC. Multiple minima were obtained for the ring twist (analogous to roll for an airplane). The phosphates at position 1 of PIP₂ and PIP₃ are within an Ångström of the plane formed by the phosphates of POPC; lipids in the surrounding shell are depressed by 0.5-0.8 Å, but otherwise the phosphoinositides do not substantially perturb the bilayer. Finite size artifacts for ion distributions are apparent for systems of ∼26 waters/lipid, but, based on simulations with a fourfold increase of the aqueous phase, the phosphoinositide positions and orientations do not show significant size effects. Electrostatic potentials evaluated from Poisson-Boltzmann (PB) calculations show a strong dependence of potential height and ring orientation, with the maxima on the −25 mV surfaces (17.1 ± 0.1 Å for PIP₂ and 19.4 ± 0.3 Å for PIP₃) occurring near the most populated orientations from MD. These surfaces are well above the background height of 10 Å estimated for negatively charged cell membranes, as would be expected for lipids involved in cellular signaling. PB calculations on microscopically flat bilayers yield similar maxima as the MD-based (microscopically rough) systems, but show less fine structure and do not clearly indicate the most probable regions. Electrostatic free energies of interaction with pentalysine are also similar for the rough and flat systems. These results support the utility of a rigid/flat bilayer model for PB-based studies of PIP₂ and PIP₃ as long as the orientations are judiciously chosen.     

Fusogenic Alzheimer's peptide fragment Abeta (29-42) in interaction with lipid bilayers: secondary structure, dynamics, and specific interaction with phosphatidyl ethanolamine polar heads as revealed by solid-state NMR

The interaction of the native Alzheimer's peptide C-terminal fragment Abeta (29-42), and two mutants (G33A and G37A) with neutral lipid bilayers made of POPC and POPE in a 9:1 molar ratio was investigated by solid-state NMR. This fragment and the lipid composition were selected because they represent the minimum requirement for the fusogenic activity of the Alzheimer's peptide. The chemical shifts of alanine methyl isotropic carbon were determined by MAS NMR, and they clearly demonstrated that the major form of the peptide equilibrated in membrane is not in a helical conformation. (2)H NMR, performed with acyl chain deuterated POPC, demonstrated that there is no perturbation of the acyl chain's dynamics and of the lipid phase transition temperature. (2)H NMR, performed with alanine methyl-deuterated peptide demonstrated that the peptide itself has a limited mobility below and above the lipid phase transition temperature (molecular order parameter equal to 0.94). MAS (31)P NMR revealed a specific interaction with POPE polar head as seen by the enhancement of POPE phosphorus nuclei T(2) relaxation. All these results are in favor of a beta-sheet oligomeric association of the peptide at the bilayer interface, preferentially recruiting phosphatidyl ethanolamine polar heads.     

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.     

Vstx1, a modifier of Kv channel gating, localizes to the interfacial region of lipid bilayers

VSTx1 is a tarantula venom toxin which binds to the archaebacterial voltage-gated potassium channel KvAP. VSTx1 is thought to access the voltage sensor domain of the channel via the lipid bilayer phase. In order to understand its mode of action and implications for the mechanism of channel activation, it is important to characterize the interactions of VSTx1 with lipid bilayers. Molecular dynamics (MD) simulations (for a total simulation time in excess of 0.2 micros) have been used to explore VSTx1 localization and interactions with zwitterionic (POPC) and with anionic (POPE/POPG) lipid bilayers. In particular, three series of MD simulations have been used to explore the net drift of VSTx1 relative to the center of a bilayer, starting from different locations of the toxin. The preferred location of the toxin is at the membrane/water interface. Although there are differences between POPC and POPE/POPG bilayers, in both cases the toxin forms favorable interactions at the interface, maximizing H-bonding to lipid headgroups and to water molecules while retaining interactions with the hydrophobic core of the bilayer. A 30 ns unrestrained simulation reveals dynamic partitioning of VSTx1 into the interface of a POPC bilayer. The preferential location of VSTx1 at the interface is discussed in the context of Kv channel gating models and provides support for a mode of action in which the toxin interacts with the Kv voltage sensor "paddle" formed by the S3 and S4 helices.     

The effect of fructan on membrane lipid organization and dynamics in the dry state

Fructans are a group of fructose-based oligo- and polysaccharides, which appear to be involved in membrane preservation during dehydration by interacting with the membrane lipids. To get further understanding of the protective mechanism, the consequences of the fructan-membrane lipid interaction for the molecular organization and dynamics in the dry state were studied. POPC and DMPC were investigated in the dry state by (2)H, (31)P NMR, and Fourier transform infrared spectroscopy using two types of fructan and dextran. The order-disorder transition temperature of dry POPC was reduced by 70 degrees C in the presence of fructan. Fructan increased the mobility of the acyl chains, but immobilized the lipid headgroup region. Most likely, fructans insert between the headgroups of lipids, thereby spacing the acyl chains. This results in a much lower phase transition temperature. The headgroup is immobilized by the interaction with fructan. The location of the interaction with the lipid headgroup is different for the inulin-type fructan compared to the levan-type fructan, since inulin shows interaction with the lipid phosphate group, whereas levan does not. Dextran did not influence the phase transition temperature of dry POPC showing that reduction of this temperature is not a general property of polysaccharides.     

Implicit solvent model estimates of the stability of model structures of the alamethicin channel

Alamethicin is a hydrophobic helical peptide of 20 residues, which oligomerizes to form ion-conducting channels in membranes. The behavior of an intact alamethicin channel in POPC bilayers was recently studied, using 2 ns molecular dynamics (MD) simulations of a model hexameric channel. These simulations produced numerous conformations of the channel. In the present study, we used 11 of these channel conformations and carried out continuum-solvent model calculations, similar to those used for the monomers in our previous studies, to investigate the energetics of the channel inside the lipid bilayer. Our results suggest that, out of the 11 channel conformations produced by the MD simulations, only four are stable inside the lipid bilayer, with water-to-membrane free energies of transfer ranging from ~–6 to ~–10 kcal/mol. Analysis of the results suggests two causes for the apparent instability of the remainder of the structures inside the lipid bilayer, both resulting from the desolvation of channel polar groups (i.e. their transfer from the aqueous phase into the bilayer). The first is specific, uncompensated backbone hydrogen bonds, which exist in the region of the channel exposed to the hydrocarbon of the lipid bilayer. The second is exposure of intra-pore water molecules to the surrounding lipid. Thus, the association of these structures with the membrane involves a large electrostatic desolvation free-energy penalty. The apparent conflict between continuum-solvent and MD calculations, and its significance for the interpretation of membrane proteins simulations, are discussed.     

31P and 2H relaxation studies of helix VII and the cytoplasmic helix of the human cannabinoid receptors utilizing solid-state NMR techniques

Cannabinoid receptors are G-protein-coupled receptors comprised of seven transmembrane helices. We hypothesized that the extended helix of the receptor interacts differently with POPC bilayers due to the differing distribution of charged amino acid residues. To test this, hCB1(T377-E416) and hCB2(K278-H316) peptides were studied with 31P and 2H solid-state NMR spectroscopy by incorporating them into 1-palmitoyl-2-oleoyl-sn-glycerophosphocholine bilayers. Lipid affinities of the 40- and 39-residue peptides were analyzed on the basis of 31P and 2H spectral line shapes, order parameters, and T1 relaxation measurements of the POPC bilayers. Lipid headgroup perturbations were noticed in the 31P NMR spectra in the lipid/peptide mixtures when compared with the pure lipids. 2H order parameters were calculated from the quadrupolar splitting of the de-Paked 2H NMR spectra. At the top of the acyl chain, pure lipids had an average S(CD) approximately = 0.20, whereas S(CD) approximately = 0.16 and S(CD) approximately = 0.18 were found in the presence of hCB1(T377-E416) and hCB2(K278-H316), respectively. S(CD) values decreased in the central part of the acyl chains when compared to the pure POPC lipids, indicating a change in the dynamic properties of the lipid membrane in the presence of the cannabinoid peptides. R(1Z) vs S2(CD) plots exhibited a linear dependency with and without the peptides, with an increase in slope upon addition of the peptides to the POPC, indicating that the dynamics of the lipid bilayer is dominated by fast axially symmetric motion. This study provides insights into the interaction of cannabinoid peptides with the membrane bilayer by investigating the headgroup and acyl chain dynamics.     

2-Arachidonoyl-sn-glycero-3-phosphoinositol: a possible natural ligand for GPR55

GPR55 is a G protein-coupled receptor. Recently, we obtained evidence that lysophosphatidylinositol (LPI) is a possible endogenous ligand for GPR55. However, no information is currently available concerning the biological activities of the individual molecular species of LPI. Furthermore, little is known concerning the levels as well as the molecular species of LPI in mammalian tissues. In this study, we first examined whether LPI is present in rat brain. We found that rat brain contains 37.5 nmol/g tissue of LPI; the most predominant fatty acyl moiety is stearic acid (50.5%) followed by arachidonic acid (22.1%). We next compared the biological activities of various molecular species of LPI and related molecules using HEK293 cells expressing GPR55. We found that the level of biological activity of the 2-arachidonoyl species is markedly higher than those of others. These results strongly suggest that the 2-arachidonoyl species of LPI is the true natural ligand for GPR55.     

Effect of salt on the interactions of antimicrobial peptides with zwitterionic lipid bilayers

The effect of salt on the binding of the antimicrobial peptide magainin to POPC lipid bilayers is studied by 40-50 ns molecular dynamics simulations of a POPC bilayer in the presence of different concentrations of Na+ and Cl− ions, corresponding to effective concentrations of 0, 100, 150, 200, 250 and 300 millimolar NaCl, with and without a single molecule of antimicrobial peptide magainin. Simulations without magainin showed that increasing salt concentration leads to the decrease in the area per lipid, a decrease in the head group tilt of the lipids, as well as increased order of lipid tails, in agreement with other recent simulations. Simulations with magainin show that peptide binding to the lipids is stronger at lower concentrations of salt. The peptides disorder the lipids in their immediate vicinity, but this effect is diminished as the salt concentration increases. Our studies indicate that while 50 ns simulations give information on peptide hydrogen bonding and lipid tail ordering that is insensitive to the initial peptide orientation, this run time is not sufficient to equilibrate the peptide position and orientation within the bilayer.     

Molecular dynamics simulation of neuropeptide B and neuropeptide W in the dipalmitoylphosphatidylcholine membrane bilayer

GPR7 and GPR8 are recently deorphanized G-protein-coupled receptors that are implicated in the regulation of neuroendocrine function, feeding behavior, and energy homeostasis. Neuropeptide B (NPB) and neuropeptide W (NPW) are two membrane-bound hypothalamic peptides, which specifically antagonize GPR7 and GPR8. Despite years of research, an accurate estimation of structure and molecular recognition of these neuropeptide systems still remains elusive. Herein, we investigated the structure, orientation, and interaction of NPB and NPW in a dipalmitoylphosphatidylcholine bilayer using long-range molecular dynamics (MD) simulation. During 30-ns simulation, membrane-embedded helical axes of NPB and NPW tilted 30 and 15°, respectively, from the membrane normal in order to overcome possible hydrophobic mismatch with the lipid bilayer. The calculation of various structural parameters indicated that NPW is more rigid and compact as compared to NPB. Qualitatively, the peptides exhibited flexible N-terminal (residues 1–12) and rigid C-terminal α-helical parts (residues 13–21), confirming previous NMR data. A strong electrostatic attraction between C-termini and headgroup atoms caused translocation of the peptides towards lower leaflet of the bilayer. The stabilizing hydrogen bonds (H-bonds) between phosphate groups and Trp1, Lys3, and Arg15 of the peptides played important roles for membrane anchoring. MD simulations of Alanine (Ala) mutants revealed that WYK->Ala variant of NPB/NPW lacked crucial H-bond interactions with phospholipid headgroups and also caused severe misfolding in NPB. Altogether, the knowledge of preferred structural fold and interaction of neuropeptides within the membrane bilayer will be useful to develop synthetic agonist or antagonist peptides for GPR7 and GPR8.     

Membrane model for the G-protein-coupled receptor rhodopsin: hydrophobic interface and dynamical structure

Rhodopsin is the only member of the pharmacologically important superfamily of G-protein-coupled receptors with a known structure at atomic resolution. A molecular dynamics model of rhodopsin in a POPC phospholipid bilayer was simulated for 15 ns, revealing a conformation significantly different from the recent crystal structures. The structure of the bilayer compared with a protein-free POPC control indicated hydrophobic matching with the nonpolar interface of the receptor, in agreement with deuterium NMR experiments. A new generalized molecular surface method, based on a three-dimensional Voronoi cell construction for atoms with different radii, was developed to quantify cross-sectional area profiles for the protein, lipid acyl chains and headgroups, and water. Thus, it was possible to investigate the bilayer deformation due to curvature of the individual lipid monolayers. Moreover, the generalized molecular surface derived hydrophobic interface allowed benchmarking of the hydropathy sequence analysis, an important structural genomics tool. Five water molecules diffused into internal hydration sites during the simulation, yielding a total of 12 internal waters. The cytoplasmic loops and the C-terminal tail, containing the G-protein recognition and protein sorting sequences, exhibited a high mobility, in marked contrast to the extracellular and transmembrane domains. The proposed functional coupling of the highly conserved ERY motif to the lipid-water interface via the cytoplasmic loops provides insight into lipid effects on G-protein-coupled receptor activation in terms of a flexible surface model, involving the spontaneous monolayer curvature.