Investigation of membrane penetration depth and interactions of the amino-terminal domain of huntingtin: refined analysis by tryptophan fluorescence measurement

The membrane-association properties of the amino-terminal domain of huntingtin are accompanied by subcellular redistribution of the protein in cellular compartments. In this study we used tryptophan substitution of amino-acid residues at different positions of the huntingtin 1–17 domain (Htt17) to precisely determine, for the first time, the depth of penetration of the peptides within the lipid bilayer. Initially, secondary structure preferences and membrane association properties were quantitatively determined for several membrane lipid compositions; they were found to be closely related to those of the natural peptide, indicating that changes in the sequence had little effect on these characteristics of the domain. The tryptophan-substituted peptides became inserted into the membranes’ interfacial region, with average tryptophan positions between 7.5 and 11 Å from the bilayer center, in agreement with in-plane orientation of the peptide. Participation of the very-amino terminus of the peptide in the membrane-association process was demonstrated. The results not only revealed the occurrence of association intermediates when the huntingtin 1–17 anchoring sequence became inserted into the membrane but also suggest the formation of aggregates and/or oligomers during membrane association. When inserted, the F11W site was of crucial importance in lipid anchoring and stabilization of the whole peptide, whereas the terminal residues are located close to the membrane surface. The carboxy-terminal tryptophan (F17W), which also constitutes the site of the polyglutamine extension in the natural domain, was found closest to the aqueous environment, accompanied with the highest aqueous quenching constants. These results were used to propose a refined model of lipid interactions of the huntingtin 1–17 domain.     

Concentration effect of cimetidine with POPC bilayer: a molecular dynamics simulation study

All-atom molecular dynamics simulations have been performed on cimetidine in the presence of a palmitoyloleoylphosphatidylcholine (POPC) bilayer. The free energy profile of a single cimetidine molecule passing across POPC bilayer displays a minimum at the interface of bilayer and water. Ten cimetidine molecules were inserted into POPC bilayer to obtain an 8 mol % drug model, and molecular dynamics results showed that cimetidine molecules reside at the polar region of POPC bilayer with sulphur atoms directing to the hydrophobic region. By comparing the one drug model with 8 mol % drug model, one can see that the central barrier to cross the membrane increases while the free energy in bulk water decreases, indicating that the ability of cimetidine passing across the POPC bilayer weakens at increased concentration. In addition, the free energy minimum shifts closer to the hydrophobic core. Our results indicate that with the increased drug concentration, it is more difficult for cimetidine to enter and pass across POPC bilayer.     

Huntington’s Disease: The Past,Present, and Future Search for Disease Modifiers

Huntington’s disease (HD) is an autosomal dominant genetic disorder that specifically causes neurodegeneration of striatal neurons, resulting in a triad of symptoms that includes emotional, cognitive, and motor disturbances. The HD mutation causes a polyglutamine repeat expansion within the N-terminal of the huntingtin (Htt) protein. This expansion causes aggregate formation within the cytosol and nucleus due to the presence of misfolded mutant Htt, as well as altered interactions with Htt’s multiple binding partners, and changes in post-translational Htt modifications. The present review charts efforts toward a therapy that delays age of onset or slows symptom progression in patients affected by HD, as there is currently no effective treatment. Although silencing Htt expression appears promising as a disease modifying treatment, it should be attempted with caution in light of Htt’s essential roles in neural maintenance and development. Other therapeutic targets include those that boost aggregate dissolution, target excitotoxicity and metabolic issues, and supplement growth factors.     

The structure of the CD3 ζζ transmembrane dimer in POPC and raft-like lipid bilayer: A molecular dynamics study

Plasma membrane lipids significantly affect assembly and activity of many signaling networks. The present work is aimed at analyzing, by molecular dynamics simulations, the structure and dynamics of the CD3 ζζ dimer in palmitoyl-oleoyl-phosphatidylcholine bilayer (POPC) and in POPC/cholesterol/sphingomyelin bilayer, which resembles the raft membrane microdomain supposed to be the site of the signal transducing machinery. Both POPC and raft-like environment produce significant alterations in structure and flexibility of the CD3 ζζ with respect to nuclear magnetic resonance (NMR) model: the dimer is more compact, its secondary structure is slightly less ordered, the arrangement of the Asp6 pair, which is important for binding to the Arg residue in the alpha chain of the T cell receptor (TCR), is stabilized by water molecules. Different interactions of charged residues with lipids at the lipid–cytoplasm boundary occur when the two environments are compared. Furthermore, in contrast to what is observed in POPC, in the raft-like environment correlated motions between transmembrane and cytoplasmic regions are observed. Altogether the data suggest that when the TCR complex resides in the raft domains, the CD3 ζζ dimer assumes a specific conformation probably necessary to the correct signal transduction.     

Molecular dynamics simulation of lipid reorientation at bilayer edges

Understanding cellular membrane processes is critical for the study of events such as viral entry, neurotransmitter exocytosis, and immune activation. Supported lipid bilayers are commonly used to model these membrane processes experimentally. Despite the relative simplicity of such a system, many important structural and dynamic parameters are not experimentally observable with current techniques. Computational approaches allow the development of a high-resolution model of bilayer processes. We have performed molecular dynamics simulations of dimyristoylphosphatidylcholine (DMPC) bilayers to model the creation of bilayer gaps—a common process in bilayer patterning—and to analyze their structure and dynamics. We propose a model for gap formation in which the bilayer edges form metastable micelle-like structures on a nanosecond timescale. Molecules near edges structurally resemble lipids in ungapped bilayers but undergo small-scale motions more rapidly. These data suggest that lipids may undergo rapid local rearrangements during membrane fusion, facilitating the formation of fusion intermediates thought key to the infection cycle of viruses such as influenza, Ebola, and HIV.     

Structure, topology, and tilt of cell-signaling peptides containing nuclear localization sequences in membrane bilayers determined by solid-state NMR and molecular dynamics simulation studies

Cell-signaling peptides have been extensively used to transport functional molecules across the plasma membrane into living cells. These peptides consist of a hydrophobic sequence and a cationic nuclear localization sequence (NLS). It has been assumed that the hydrophobic region penetrates the hydrophobic lipid bilayer and delivers the NLS inside the cell. To better understand the transport mechanism of these peptides, in this study, we investigated the structure, orientation, tilt of the peptide relative to the bilayer normal, and the membrane interaction of two cell-signaling peptides, SA and SKP. Results from CD and solid-state NMR experiments combined with molecular dynamics simulations suggest that the hydrophobic region is helical and has a transmembrane orientation with the helical axis tilted away from the bilayer normal. The influence of the hydrophobic mismatch, between the hydrophobic length of the peptide and the hydrophobic thickness of the bilayer, on the tilt angle of the peptides was investigated using thicker POPC and thinner DMPC bilayers. NMR experiments showed that the hydrophobic domain of each peptide has a tilt angle of 15 +/- 3 degrees in POPC, whereas in DMPC, 25 +/- 3 degree and 30 +/- 3 degree tilts were observed for SA and SKP peptides, respectively. These results are in good agreement with molecular dynamics simulations, which predict a tilt angle of 13.3 degrees (SA in POPC), 16.4 degrees (SKP in POPC), 22.3 degrees (SA in DMPC), and 31.7 degrees (SKP in DMPC). These results and simulations on the hydrophobic fragment of SA or SKP suggest that the tilt of helices increases with a decrease in bilayer thickness without changing the phase, order, and structure of the lipid bilayers.     

Structure and Topology of the Huntingtin 1–17 Membrane Anchor by a Combined Solution and Solid-State NMR Approach

The very amino-terminal domain of the huntingtin protein is directly located upstream of the protein’s polyglutamine tract, plays a decisive role in several important properties of this large protein and in the development of Huntington’s disease. This huntingtin 1–17 domain is on the one hand known to markedly increase polyglutamine aggregation rates and on the other hand has been shown to be involved in cellular membrane interactions. Here, we determined the high-resolution structure of huntingtin 1–17 in dodecyl phosphocholine micelles and the topology of its helical domain in oriented phosphatidylcholine bilayers. Using two-dimensional solution NMR spectroscopy the low-energy conformations of the polypeptide were identified in the presence of dodecyl phosphocholine detergent micelles. In a next step a set of four solid-state NMR angular restraints was obtained from huntingtin 1–17 labeled with ¹⁵N and ²H at selected sites. Of the micellar ensemble of helical conformations only a limited set agrees in quantitative detail with the solid-state angular restraints of huntingtin 1–17 obtained in supported planar lipid bilayers. Thereby, the solid-state NMR data were used to further refine the domain structure in phospholipid bilayers. At the same time its membrane topology was determined and different motional regimes of this membrane-associated domain were explored. The pronounced structural transitions of huntingtin 1–17 upon membrane-association result in a α-helical conformation from K6 to F17, i.e., up to the very start of the polyglutamine tract. This amphipathic helix is aligned nearly parallel to the membrane surface (tilt angle ∼77°) and is characterized by a hydrophobic ridge on one side and an alternation of cationic and anionic residues that run along the hydrophilic face of the helix. This arrangement facilitates electrostatic interactions between huntingtin 1–17 domains and possibly with the proximal polyglutamine tract.     

NMR-Based Simulation Studies of Pf1 Coat Protein in Explicit Membranes

As time- and ensemble-averaged measures, NMR observables contain information about both protein structure and dynamics. This work represents a computational study to extract such information for membrane proteins from orientation-dependent NMR observables: solid-state NMR chemical shift anisotropy and dipolar coupling, and solution NMR residual dipolar coupling. We have performed NMR-restrained molecular dynamics simulations to refine the structure of the membrane-bound form of Pf1 coat protein in explicit lipid bilayers using the recently measured chemical shift anisotropy, dipolar coupling, and residual dipolar coupling data. From the simulations, we have characterized detailed protein-lipid interactions and explored the dynamics. All simulations are stable and the NMR restraints are well satisfied. The C-terminal transmembrane (TM) domain of Pf1 finds its optimal position in the membrane quickly (within 6 ns), illustrating efficient solvation of TM domains in explicit bilayer environments. Such rapid convergence also leads to well-converged interaction patterns between the TM helix and the membrane, which clearly show the interactions of interfacial membrane-anchoring residues with the lipids. For the N-terminal periplasmic helix of Pf1, we identify a stable, albeit dynamic, helix orientation parallel to the membrane surface that satisfies the amphiphatic nature of the helix in an explicit lipid bilayer. Such detailed information cannot be obtained solely from NMR observables. Therefore, the present simulations illustrate the usefulness of NMR-restrained MD refinement of membrane protein structure in explicit membranes.     

Polyglutamine expansion in huntingtin increases its insertion into lipid bilayers

An expanded polyglutamine (Q) tract (>37Q) in huntingtin (htt) causes Huntington disease. Htt associates with membranes and polyglutamine expansion in htt may alter membrane function in Huntington disease through a mechanism that is not known. Here we used differential scanning calorimetry to examine the effects of polyQ expansion in htt on its insertion into lipid bilayers. We prepared synthetic lipid vesicles composed of phosphatidylcholine and phosphatidylethanolamine and tested interactions of htt amino acids 1-89 with 20Q, 32Q or 53Q with the vesicles. GST-htt1-89 with 53Q inserted into synthetic lipid vesicles significantly more than GST-htt1-89 with 20Q or 32Q. We speculate that by inserting more into cell membranes, mutant huntingtin could increase disorder within the lipid bilayer and thereby disturb cellular membrane function.     

Dynamics properties of membrane proteins in native cell membranes revealed by solid-state NMR spectroscopy

Cell membranes provide an environment that is essential to the functions of membrane proteins. Cell membranes are mainly composed of proteins and highly diverse phospholipids. The influence of diverse lipid compositions of native cell membranes on the dynamics of the embedded membrane proteins has not been examined. Here we employ solid-state NMR to investigate the dynamics of E. coli Aquaporin Z (AqpZ) in its native inner cell membranes, and reveal the influence of diverse lipid compositions on the dynamics of AqpZ by comparing it in native cell membranes to that in POPC/POPG bilayers. We demonstrate that the dynamic rigidity of AqpZ generally conserves in both native cell membranes and POPC/POPG bilayers, due to its tightly packed tetrameric structure. In the gel and the liquid crystal phases of lipids, our experimental results show that AqpZ is more dynamic in native cell membranes than that in POPC/POPG bilayers. In addition, we observe that phase transitions of lipids in native membranes are less sensitive to temperature variations compared with that in POPC/POPG bilayers, which results in that the dynamics of AqpZ is less affected by the phase transitions of lipids in native cell membranes than that in POPC/POPG bilayers. This study provides new insights into the dynamics of membrane proteins in native cell membranes.     

Immobilization and aggregation of the antimicrobial peptide protegrin-1 in lipid bilayers investigated by solid-state NMR

The dynamics and aggregation of a beta-sheet antimicrobial peptide, protegrin-1 (PG-1), are investigated using solid-state NMR spectroscopy. Chemical shift anisotropies of F12 and V16 carbonyl carbons are uniaxially averaged in 1,2-dilauryl-sn-glycero-3-phosphatidylcholine (DLPC) bilayers but approach rigid-limit values in the thicker 1-palmitoyl-2-oleoyl-sn-glycerol-3-phosphatidylcholine (POPC) bilayers. The Calpha-Halpha dipolar coupling of L5 is scaled by a factor of 0.16 in DLPC bilayers but has a near-unity order parameter of 0.96 in POPC bilayers. The larger couplings of PG-1 in POPC bilayers indicate immobilization of the peptide, suggesting that PG-1 forms oligomeric aggregates at the biologically relevant bilayer thickness. Exchange NMR experiments on F12 (13)CO-labeled PG-1 show that the peptide undergoes slow reorientation with a correlation time of 0.7 +/- 0.2 s in POPC bilayers. This long correlation time suggests that in addition to aggregation, geometric constraints in the membrane may also contribute to PG-1 immobilization. The PG-1 aggregates contact both the surface and the hydrophobic center of the POPC bilayer, as determined by (1)H spin-diffusion measurements. Thus, solid-state NMR provides a wide range of information about the molecular details of membrane peptide immobilization and aggregation in lipid bilayers.     

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.     

Integrating Solid-State NMR and Computational Modeling to Investigate the Structure and Dynamics of Membrane-Associated Ghrelin

The peptide hormone ghrelin activates the growth hormone secretagogue receptor 1a, also known as the ghrelin receptor. This 28-residue peptide is acylated at Ser3 and is the only peptide hormone in the human body that is lipid-modified by an octanoyl group. Little is known about the structure and dynamics of membrane-associated ghrelin. We carried out solid-state NMR studies of ghrelin in lipid vesicles, followed by computational modeling of the peptide using Rosetta. Isotropic chemical shift data of isotopically labeled ghrelin provide information about the peptide’s secondary structure. Spin diffusion experiments indicate that ghrelin binds to membranes via its lipidated Ser3. Further, Phe4, as well as electrostatics involving the peptide’s positively charged residues and lipid polar headgroups, contribute to the binding energy. Other than the lipid anchor, ghrelin is highly flexible and mobile at the membrane surface. This observation is supported by our predicted model ensemble, which is in good agreement with experimentally determined chemical shifts. In the final ensemble of models, residues 8–17 form an α-helix, while residues 21–23 and 26–27 often adopt a polyproline II helical conformation. These helices appear to assist the peptide in forming an amphipathic conformation so that it can bind to the membrane.     

Transmembrane helix-helix interactions: comparative simulations of the glycophorin a dimer

The glycophorin helix dimer is a paradigm for the exploration of helix-helix interactions in integral membrane proteins. Two NMR structures of the dimer are known, one in a detergent micelle and one in a lipid bilayer. Multiple (4 x 50 ns) molecular dynamics simulations starting from each of the two NMR structures, with each structure in either a dodecyl phosphocholine (DPC) micelle or a dimyristoyl phosphatidylcholine (DMPC) bilayer, have been used to explore the conformational dynamics of the helix dimer. Analysis of the helix-helix interaction, mediated by the GxxxG sequence motif, suggests convergence of the simulations to a common model. This is closer to the NMR structure determined in a bilayer than to micelle structure. The stable dimer interface in the final simulation model is characterized by (i) Gly/Gly packing and (ii) Thr/Thr interhelix H-bonds. These results demonstrate the ability of extended molecular dynamics simulations in a lipid bilayer environment to refine membrane protein structures or models derived from experimental data obtained in protein/detergent micelles.     

An N-terminal Nuclear Export Signal Regulates Trafficking and Aggregation of Huntingtin (Htt) Protein Exon 1

Huntington disease is a dominantly inherited neurodegenerative condition caused by polyglutamine expansion in the N terminus of the huntingtin protein (Htt). The first 17 amino acids (N17) of Htt play a key role in regulating its toxicity and aggregation. Both nuclear export and cytoplasm retention functions have been ascribed to N17. We have determined that N17 acts as a nuclear export sequence (NES) within Htt exon and when fused to yellow fluorescent protein. We have defined amino acids within N17 that constitute the nuclear export sequence (NES). Mutation of any of the conserved residues increases nuclear accumulation of Htt exon 1. Nuclear export of Htt is sensitive to leptomycin B and is reduced by knockdown of exportin 1. In HEK293 cells, NES mutations decrease overall Htt aggregation but increase the fraction of cells with nuclear inclusions. In primary cultured neurons, NES mutations increase nuclear accumulation and increase overall aggregation. This work defines a bona fide nuclear export sequence within N17 and links it to effects on protein aggregation. This may help explain the important role of N17 in controlling Htt toxicity.     

Solid-state nuclear magnetic resonance relaxation studies of the interaction mechanism of antimicrobial peptides with phospholipid bilayer membranes

An 18-residue peptide, KWGAKIKIGAKIKIGAKI-NH(2) was designed to form amphiphilic beta-sheet structures when bound to lipid bilayers. The peptide possesses high antimicrobial activity when compared to naturally occurring linear antimicrobial peptides, most of which adopt an amphipathic alpha-helical conformation upon binding to the lipids. The perturbation of the bilayer by the peptide was studied by static (31)P and (2)H solid-state NMR spectroscopy using POPC and POPG/POPC (3/1) bilayer membranes with sn-1 chain perdeuterated POPC and POPG as the isotopic labels. (31)P NMR powder spectra exhibited two components for POPG/POPC bilayers upon addition of the peptide but only a slight change in the line shape for POPC bilayers, indicating that the peptide selectively disrupted the membrane structure consisting of POPG lipids. (2)H NMR powder spectra indicated a reduction in the lipid chain order for POPC bilayers and no significant change in the ordering for POPG/POPC bilayers upon association of the peptide with the bilayers, suggesting that the peptide acts as a surface peptide in POPG/POPC bilayers. Relaxation rates are more sensitive to the motions of the membranes over a large range of time scales. Longer (31)P longitudinal relaxation times for both POPG and POPC in the presence of the peptide indicated a direct interaction between the peptide and the POPG/POPC bilayer membranes. (31)P longitudinal relaxation studies also suggested that the peptide prefers to interact with the POPG phospholipids. However, inversion-recovery (2)H NMR spectroscopic experiments demonstrated a change in the relaxation rate of the lipid acyl chains for both the POPC membranes and the POPG/POPC membranes upon interaction with the peptide. Transverse relaxation studies indicated an increase in the spectral density of the collective membrane motion caused by the interaction between the peptide and the POPG/POPC membrane. The experimental results demonstrate significant dynamic changes in the membrane in the presence of the antimicrobial peptide and support a carpet mechanism for the disruption of the membranes by the antimicrobial peptide.     

Effects of a carane derivative local anesthetic on a phospholipid bilayer studied by molecular dynamics simulation

Molecular dynamics (MD) simulations of two hydrated palmitoyloleoylphosphatidylcholine (POPC) bilayers each containing eight carane derivative (KP-23) local anesthetic (LA) molecules in neutral (POPC-LA) or protonated (POPC-LAH) forms were carried out to investigate the effect of KP-23 and its protonation on the bilayer. 3-ns trajectories were used for analyses. A pure POPC bilayer was employed as a reference system. In both POPC-LA and POPC-LAH systems a few KP-23 molecules intercalated into the bilayer and moved near the bilayer/water interface. They were located on the hydrophobic core side of the interface in the POPC-LA bilayer, but on the water phase side in the POPC-LAH bilayer. The order of the POPC chains was higher in the POPC-LA bilayer than in the pure POPC bilayer and was lower in the POPC-LAH bilayer. Interactions between polar groups of KP-23 and POPC or water were responsible for a lower hydration of POPC headgroups in POPC bilayers containing KP-23 than in the pure POPC bilayer. KP-23 molecules were found to form aggregates both in POPC-LA and POPC-LAH bilayers. Due to higher amphiphilicity of LAH, the LAH aggregate was more micelle-like and larger than the LA one. The results demonstrate the rapid timescales of the initial processes that take place at and near the bilayer interface as well as details of the atomic level interactions between local anesthetic and the lipid matrix of a cell membrane.     

Location, Structure, and Dynamics of the Synthetic Cannabinoid Ligand CP-55,940 in Lipid Bilayers

The widely used hydrophobic cannabinoid ligand CP-55,940 partitions with high efficiency into biomembranes. We studied the location, orientation, and dynamics of CP-55,940 in POPC bilayers by solid-state NMR. Chemical-shift perturbation of POPC protons from the aromatic ring-current effect, as well as ¹H NMR cross-relaxation rates, locate the hydroxyphenyl ring of the ligand near the lipid glycerol, carbonyls, and upper acyl-chain methylenes. Order parameters of the hydroxyphenyl ring determined by the ¹H-¹³C DIPSHIFT experiment indicate that the bond between the hydroxyphenyl and hydroxycyclohexyl rings is oriented perpendicular to the bilayer normal. ²H NMR order parameters of the nonyl tail are very low, indicating that the hydrophobic chain maintains a high level of conformational flexibility in the membrane. Lateral diffusion rates of CP-55,940 and POPC were measured by ¹H magic-angle spinning NMR with pulsed magnetic field gradients. The rate of CP-55,940 diffusion is comparable to the rate of lipid diffusion. The magnitude of cross-relaxation and diffusion rates suggests that associations between CP-55,940 and lipids are with lifetimes of a fraction of a microsecond. With its flexible hydrophobic tail, CP-55,940 may efficiently approach the binding site of the cannabinoid receptor from the lipid-water interface by lateral diffusion.     

Comparative Molecular Dynamics Simulation Studies of Protegrin-1 Monomer and Dimer in Two Different Lipid Bilayers

Antimicrobial peptides interact specifically with the membrane of a pathogen and kill the pathogen by releasing its cellular contents. Protegrin-1 (PG-1), a β-hairpin antimicrobial peptide, is known to exist as a transmembrane monomer in a 1,2-dilauroylphosphatidylcholine (DLPC) bilayer and shows concentration-dependent oligomerization in a 1-palmitoyl-2-oleoylphosphatidylcholine (POPC) bilayer. To examine its structure, dynamics, orientation, and interaction in membranes, we performed comparative molecular dynamics simulations of PG-1 monomer and dimer in DLPC and POPC bilayers for a total of 840 ns. The PG-1 monomer exhibits larger tilting in DLPC than in POPC due to a hydrophobic mismatch. PG-1 tilting is dependent on its rotation angle. The specific orientation of PG-1 in membranes is governed by the interactions of its aromatic residues with lipid headgroups. The calculated ¹⁵N and ¹³CO chemical shifts of Val¹⁶ in DLPC reveal that there are different sets of tilt and rotation angles that satisfy the experimental values reasonably, suggesting that more experiments are needed to determine its orientation. The dimer simulations show that the dimer interface is better preserved in POPC than in DLPC because POPC's greater hydrophobic thickness causes reduced flexibility of the C-terminal strands. Both monomer and dimer simulations show membrane thinning around PG-1, largely due to arginine-lipid interactions.