Pinholin S21 mutations induce structural topology and conformational changes

The bacteriophage infection cycle is terminated at a predefined time to release the progeny virions via a robust lytic system composed of holin, endolysin, and spanin proteins. Holin is the timekeeper of this process. Pinholin S²¹ is a prototype holin of phage Φ21, which determines the timing of host cell lysis through the coordinated efforts of pinholin and antipinholin. However, mutations in pinholin and antipinholin play a significant role in modulating the timing of lysis depending on adverse or favorable growth conditions. Earlier studies have shown that single point mutations of pinholin S²¹ alter the cell lysis timing, a proxy for pinholin function as lysis is also dependent on other lytic proteins. In this study, continuous wave electron paramagnetic resonance (CW-EPR) power saturation and double electron-electron resonance (DEER) spectroscopic techniques were used to directly probe the effects of mutations on the structure and conformational changes of pinholin S²¹ that correlate with pinholin function. DEER and CW-EPR power saturation data clearly demonstrate that increased hydrophilicity induced by residue mutations accelerate the externalization of antipinholin transmembrane domain 1 (TMD1), while increased hydrophobicity prevents the externalization of TMD1. This altered hydrophobicity is potentially accelerating or delaying the activation of pinholin S²¹. It was also found that mutations can influence intra- or intermolecular interactions in this system, which contribute to the activation of pinholin and modulate the cell lysis timing. This could be a novel approach to analyze the mutational effects on other holin systems, as well as any other membrane protein in which mutation directly leads to structural and conformational changes.     

Mutational analysis of the S21 pinholin

Lambdoid phage 21 has the prototype pinholin‐SAR endolysin lysis system, which is widely distributed among phages. Its prototype pinholin, S²¹68, triggers at an allele‐specific time to form small, heptameric lesions, or pinholes, in the cytoplasmic membrane, thus initiating lysis. S²¹68 has two transmembrane domains, TMD1 and TMD2. Only TMD2 is required for the formation of pinholes, whereas TMD1 acts as an inhibitor of TMD2 and must be externalized to the periplasm in the lytic pathway. Previously we provided evidence that S²¹68 first accumulates as inactive dimers with both transmembrane domains embedded in the bilayer. Here we analyse an extensive collection of S²¹ mutants to identify residues and domains critical to the function and regulation of the pinholin. Evidence is presented indicating that, within the inactive dimer, TMD1 acts in trans as an inhibitor of the lethal function of TMD2. A wide range of phenotypes, from absolute lysis defectives to accelerated lysis triggering, are observed for mutations mapping to each topological domain. The pattern of phenotypes allows the generation of a model for the structure of the inactive dimer. The model identifies the faces of the two transmembrane domains involved in intramolecular and intermolecular interactions, as well as interaction with the lipid.     

Probing the local secondary structure of bacteriophage S21 pinholin membrane protein using electron spin echo envelope modulation spectroscopy

There have recently been advances in methods for detecting local secondary structures of membrane protein using electron paramagnetic resonance (EPR). A three pulsed electron spin echo envelope modulation (ESEEM) approach was used to determine the local helical secondary structure of the small hole forming membrane protein, S²¹ pinholin. This ESEEM approach uses a combination of site-directed spin labeling and ²H-labeled side chains. Pinholin S²¹ is responsible for the permeabilization of the inner cytosolic membrane of double stranded DNA bacteriophage host cells. In this study, we report on the overall global helical structure using circular dichroism (CD) spectroscopy for the active form and the negative-dominant inactive mutant form of S²¹ pinholin. The local helical secondary structure was confirmed for both transmembrane domains (TMDs) for the active and inactive S²¹ pinholin using the ESEEM spectroscopic technique. Comparison of the ESEEM normalized frequency domain intensity for each transmembrane domain gives an insight into the α-helical folding nature of these domains as opposed to a π or 3₁₀-helix which have been observed in other channel forming proteins.     

Structure and transmembrane topology of Slc11a1 TMD1-5 in lipid membranes

The importance of solute carrier family 11 (Slc11) in divalent metal-ion transport has been well established. The core domains TMD1-5 and TMD6-10 of the proteins were modeled as a symmetric but inversely orientated arrangement with respect to membrane normal. In this article, the structures and transmembrane topologies of TMD1-5 of Slc11a1 incorporated with phospholipids 1-palmitoyl-2-oleoyl- sn-glycero-3-phosphocholine (POPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-[phospho-rac-(1-glycerol)] (sodium salt) (POPG), and POPC/POPG (3:1) were explored using circular dichroism, fluorescence, and attenuated total reflection Fourier transform infrared (ATR-FTIR) spectroscopies. The segments TMD2-5 were inserted in lipid membranes mainly as an α-helix with orientations of helices along membrane normal. The tilt angles of the helices were in an order of TMD3 > TMD4 > TMD2 > TMD5 in these membranes. In contrast, TMD1 was partly inserted in membranes, leaving partial segment at membrane surface. The amount of the lipid component with negatively charged headgroups had an effect on both the helicity and orientation of the transmembrane domains (TMDs). Nevertheless, the helices maintained similar topologies in various membranes.     

Stability and membrane orientation of the fukutin transmembrane domain: a combined multiscale molecular dynamics and circular dichroism study

The N-terminal domain of fukutin-I has been implicated in the localization of the protein in the endoplasmic reticulum and Golgi Apparatus. It has been proposed to mediate this through its interaction with the thinner lipid bilayers found in these compartments. Here we have employed multiscale molecular dynamics simulations and circular dichroism spectroscopy to explore the structure, stability, and orientation of the short 36-residue N-terminus of fukutin-I (FK1TMD) in lipids with differing tail lengths. Our results show that FK1TMD adopts a stable helical conformation in phosphatidylcholine lipids when oriented with its principal axis perpendicular to the bilayer plane. The stability of the helix is largely insensitive to the lipid tail length, preventing hydrophobic mismatch by virtue of its mobility and ability to tilt within the lipid bilayers. This suggests that changes in FK1TMD tilt in response to bilayer properties may be implicated in the regulation of its trafficking. Coarse-grained simulations of the complex Golgi membrane suggest the N-terminal domain may induce the formation of microdomains in the surrounding membrane through its preferential interaction with 1,2-dipalmitoyl-sn-glycero-3-phosphatidylinositol 4,5-bisphosphate lipids.     

Cholesterol-induced protein sorting: an analysis of energetic feasibility

The mechanism(s) underlying the sorting of integral membrane proteins between the Golgi complex and the plasma membrane remain uncertain because no specific Golgi retention signal has been found. Moreover one can alter a protein's eventual localization simply by altering the length of its transmembrane domain (TMD). M. S. Bretscher and S. Munro (SCIENCE: 261:1280-1281, 1993) therefore proposed a physical sorting mechanism based on the hydrophobic match between the proteins' TMD and the bilayer thickness, in which cholesterol would regulate protein sorting by increasing the lipid bilayer thickness. In this model, Golgi proteins with short TMDs would be excluded from cholesterol-enriched domains (lipid rafts) that are incorporated into transport vesicles destined for the plasma membrane. Although attractive, this model remains unproven. We therefore evaluated the energetic feasibility of a cholesterol-dependent sorting process using the theory of elastic liquid crystal deformations. We show that the distribution of proteins between cholesterol-enriched and cholesterol-poor bilayer domains can be regulated by cholesterol-induced changes in the bilayer physical properties. Changes in bilayer thickness per se, however, have only a modest effect on sorting; the major effect arises because cholesterol changes also the bilayer material properties, which augments the energetic penalty for incorporating short TMDs into cholesterol-enriched domains. We conclude that cholesterol-induced changes in the bilayer physical properties allow for effective and accurate sorting which will be important generally for protein partitioning between different membrane domains.     

Charged or Aromatic Anchor Residue Dependence of Transmembrane Peptide Tilt

The membrane-spanning segments of integral membrane proteins often are flanked by aromatic or charged amino acid residues, which may “anchor” the transmembrane orientation. Single spanning transmembrane peptides such as those of the WALP family, acetyl-GWW(LA)_nLWWA-amide, furthermore adopt a moderate average tilt within lipid bilayer membranes. To understand the anchor residue dependence of the tilt, we introduce Leu-Ala “spacers” between paired anchors and in some cases replace the outer tryptophans. The resulting peptides, acetyl-GX²ALW(LA)₆LWLAX²²A-amide, have Trp, Lys, Arg, or Gly in the two X positions. The apparent average orientations of the core helical sequences were determined in oriented phosphatidylcholine bilayer membranes of varying thickness using solid-state ²H NMR spectroscopy. When X is Lys, Arg, or Gly, the direction of the tilt is essentially constant in different lipids and presumably is dictated by the tryptophans (Trp⁵ and Trp¹⁹) that flank the inner helical core. The Leu-Ala spacers are no longer helical. The magnitude of the apparent helix tilt furthermore scales nicely with the bilayer thickness except when X is Trp. When X is Trp, the direction of tilt is less well defined in each phosphatidylcholine bilayer and varies up to 70° among 1,2-dioleoyl-sn-glycero-3-phosphocholine, 1,2-dimyristoyl-sn-glycero-3-phosphocholine, and 1,2-dilauroyl-sn-glycero-3-phosphocholine bilayer membranes. Indeed, the X = Trp case parallels earlier observations in which WALP family peptides having multiple Trp anchors show little dependence of the apparent tilt magnitude on bilayer thickness. The results shed new light on the interactions of arginine, lysine, tryptophan, and even glycine at lipid bilayer membrane interfaces.     

Molecular dynamics simulations of the human ocular lens with age and cataract

The human ocular lens consists primarily of elongated, static fibers characterized by high stability and low turnover, which differ dramatically in their composition and properties from other biological membranes. Cholesterol (Chol) and sphingolipids (SL) are present at high concentrations, including saturated SLs, such as dihyrosphingomyelin (DHSM). Past molecular dynamics simulations demonstrated that the presence of DHSM and high Chol concentration contributes to higher order in lipid membranes. This current study simulated more complex models of human lens membranes. Models were developed representing physiological compositions in cataractous lenses aged 74 ± 6 years and in healthy lenses aged 22 ± 4, 41 ± 6, and 69 ± 3 years. With older age, Chol and ceramide concentrations increase and glycerophospholipid concentration decreases. With cataract, ceramide concentration increases and Chol and glycerophospholipid concentrations decrease. Surface area per lipid, deuterium order parameters (S_CD), sterol tilt angle, electron density profiles, bilayer thickness, chain interdigitation, two-dimensional radial distribution functions (2D-RDF), lipid clustering, and hydrogen bonding were calculated for all simulations. All systems exhibited low surface area per lipid and high bilayer thickness, indicative of strong vertical packing. S_CD parameters suggest similarly, with saturated tails in the hydrophobic core of the membrane having elevated order. Vertical packing and acyl tail order increased with both age and cataract condition. Lateral diffusion decreased with age and cataracts, with the older and cataractous models demonstrating increased long-range structure by the 2D-RDF analysis. In future work examining the membrane proteins of the lens, these models can serve as a physiologically accurate representation of the lens lipidome.     

Membrane interface composition drives the structure and the tilt of the single transmembrane helix protein PMP1: MD studies

PMP1, a regulatory subunit of the yeast plasma membrane H⁺-ATPase, is a single transmembrane helix protein. Its cytoplasmic C-terminus possesses several positively charged residues and interacts with phosphatidylserine lipids as shown through both ¹H- and ²H-NMR experiments. We used all-atom molecular dynamics simulations to obtain atomic-scale data on the effects of membrane interface lipid composition on PMP1 structure and tilt. PMP1 was embedded in two hydrated bilayers, differing in the composition of the interfacial region. The neutral bilayer is composed of POPC (1-palmitoyl-2-oleoyl-3-glycero-phosphatidylcholine) lipids and the negatively charged bilayer is composed of POPC and anionic POPS (1-palmitoyl-2-oleoyl-3-glycero-phosphatidylserine) lipids. Our results were consistent with NMR data obtained previously, such as a lipid sn-2 chain lying on the W28 aromatic ring and in the groove formed on one side of the PMP1 helix. In pure POPC, the transmembrane helix is two residues longer than the initial structure and the helix tilt remains constant at 6 ± 3°. By contrast, in mixed POPC-POPS, the initial helical structure of PMP1 is stable throughout the simulation time even though the C-terminal residues interact strongly with POPS headgroups, leading to a significant increase of the helix tilt within the membrane to 20 ± 5°.     

Membrane partitioning of peptide aggregates: coarse-grained molecular dynamics simulations

Abstract

Coarse-grained molecular dynamics (CGMD) simulation technique (MARTINI force field) is applied to monitor the aggregation of helical peptides representing the transmembrane sequence and its extension of bone marrow stromal cell antigen 2 (BST-2). One of the peptides is coupled with a protein transducing domain (PTD) of nine arginine residues (R9) at its N-terminal side as well as a peptide, pep11**, which has been shown to bind to human papilloma virus 16 (HPV16) E6 oncoprotein. A short hydrophobic stretch of the transmembrane domain (TMD) of BST-2 aggregates the fastest and inserts into a lipid membrane. An aggregate of R9-pep11** attaches to the membrane via simultaneous contact of many arginine residues. Monomers from the aggregates of the shortest of the hydrophobic TMDs dissolve into the opposing leaflet when the aggregate spans the bilayer. A ‘flipping’ of the individual monomeric peptides is not observed.

Communicated by Ramaswamy H. Sarma     

Lipid-anchored Synaptobrevin Provides Little or No Support for Exocytosis or Liposome Fusion

SNARE proteins catalyze many forms of biological membrane fusion, including Ca²⁺-triggered exocytosis. Although fusion mediated by SNAREs generally involves proteins anchored to each fusing membrane by a transmembrane domain (TMD), the role of TMDs remains unclear, and previous studies diverge on whether SNAREs can drive fusion without a TMD. This issue is important because it relates to the question of the structure and composition of the initial fusion pore, as well as the question of whether SNAREs mediate fusion solely by creating close proximity between two membranes versus a more active role in transmitting force to the membrane to deform and reorganize lipid bilayer structure. To test the role of membrane attachment, we generated four variants of the synaptic v-SNARE synaptobrevin-2 (syb2) anchored to the membrane by lipid instead of protein. These constructs were tested for functional efficacy in three different systems as follows: Ca²⁺-triggered dense core vesicle exocytosis, spontaneous synaptic vesicle exocytosis, and Ca²⁺-synaptotagmin-enhanced SNARE-mediated liposome fusion. Lipid-anchoring motifs harboring one or two lipid acylation sites completely failed to support fusion in any of these assays. Only the lipid-anchoring motif from cysteine string protein-α, which harbors many lipid acylation sites, provided support for fusion but at levels well below that achieved with wild type syb2. Thus, lipid-anchored syb2 provides little or no support for exocytosis, and anchoring syb2 to a membrane by a TMD greatly improves its function. The low activity seen with syb2-cysteine string protein-α may reflect a slower alternative mode of SNARE-mediated membrane fusion.     

Assembling an ion channel: ORF 3a from SARS‐CoV

Protein 3a is a 274 amino acid polytopic channel protein with three putative transmembrane domains (TMDs) encoded by severe acute respiratory syndrome corona virus (SARS‐CoV). Synthetic peptides corresponding to each of its three individual transmembrane domains (TMDs) are reconstituted into artificial lipid bilayers. Only TMD2 and TMD3 induce channel activity. Reconstitution of the peptides as TMD1 + TMD3 as well as TMD2 + TMD3 in a 1 : 1 mixture induces membrane activity for both mixtures. In a 1 : 1 : 1 mixture, channel like behavior is almost restored. Expression of full length 3a and reconstitution into artificial lipid bilayers reveal a weak cation selective (P_K ≈ 2 P_Cl) rectifying channel. In the presence of nonphysiological concentration of Ca‐ions the channel develops channel activity. © 2013 Wiley Periodicals, Inc. Biopolymers 99:628–635, 2013.     

Conformation and Trimer Association of the Transmembrane Domain of the Parainfluenza Virus Fusion Protein in Lipid Bilayers from Solid-State NMR: Insights into the Sequence Determinants of Trimer Structure and Fusion Activity

Enveloped viruses enter cells by using their fusion proteins to merge the virus lipid envelope and the cell membrane. While crystal structures of the water-soluble ectodomains of many viral fusion proteins have been determined, the structure and assembly of the C-terminal transmembrane domain (TMD) remains poorly understood. Here we use solid-state NMR to determine the backbone conformation and oligomeric structure of the TMD of the parainfluenza virus 5 fusion protein. ¹³C chemical shifts indicate that the central leucine-rich segment of the TMD is α-helical in POPC/cholesterol membranes and POPE membranes, while the Ile- and Val-rich termini shift to the β-strand conformation in the POPE membrane. Importantly, lipid mixing assays indicate that the TMD is more fusogenic in the POPE membrane than in the POPC/cholesterol membrane, indicating that the β-strand conformation is important for fusion by inducing membrane curvature. Incorporation of para-fluorinated Phe at three positions of the α-helical core allowed us to measure interhelical distances using ¹⁹F spin diffusion NMR. The data indicate that, at peptide:lipid molar ratios of ~ 1:15, the TMD forms a trimeric helical bundle with inter-helical distances of 8.2–8.4 Å for L493F and L504F and 10.5 Å for L500F. These data provide high-resolution evidence of trimer formation of a viral fusion protein TMD in phospholipid bilayers, and indicate that the parainfluenza virus 5 fusion protein TMD harbors two functions: the central α-helical core is the trimerization unit of the protein, while the two termini are responsible for inducing membrane curvature by transitioning to a β-sheet conformation.     

Atomistic-level study of the interactions between hIAPP protofibrils and membranes: Influence of pH and lipid composition

The pathology of type 2 diabetes mellitus is associated with the aggregation of human islet amyloid polypeptide (hIAPP) and aggregation-mediated membrane disruption. The interactions of hIAPP aggregates with lipid membrane, as well as the effects of pH and lipid composition at the atomic level, remain elusive. Herein, using molecular dynamics simulations, we investigate the interactions of hIAPP protofibrillar oligomers with lipids, and the membrane perturbation that they induce, when they are partially inserted in an anionic dipalmitoyl-phosphatidylglycerol (DPPG) membrane or a mixed dipalmitoyl-phosphatidylcholine (DPPC)/DPPG (7:3) lipid bilayer under acidic/neutral pH conditions. We observed that the tilt angles and insertion depths of the hIAPP protofibril are strongly correlated with the pH and lipid composition. At neutral pH, the tilt angle and insertion depth of hIAPP protofibrils at a DPPG bilayer reach ~52° and ~1.62 nm with respect to the membrane surface, while they become ~77° and ~1.75 nm at a mixed DPPC/DPPG membrane. The calculated tilt angle of hIAPP at DPPG membrane is consistent with a recent chiral sum frequency generation spectroscopic study. The acidic pH induces a smaller tilt angle of ~40° and a shallower insertion depth (~1.24 nm) of hIAPP at the DPPG membrane surface, mainly due to protonation of His18 near the turn region. These differences mainly result from a combination of distinct electrostatic, van der Waals, hydrogen bonding and salt-bridge interactions between hIAPP and lipid bilayers. The hIAPP-membrane interaction energy analysis reveals that besides charged residues K1, R11 and H18, aromatic residues Phe15 and Phe23 also exhibit strong interactions with lipid bilayers, revealing the crucial role of aromatic residues in stabilizing the membrane-bound hIAPP protofibrils. hIAPP-membrane interactions disturb the lipid ordering and the local bilayer thickness around the peptides. Our results provide atomic-level information of membrane interaction of hIAPP protofibrils, revealing pH-dependent and membrane-modulated hIAPP aggregation at the early stage.     

HIV-1 Vpu protein antagonizes innate restriction factor BST-2 via lipid-embedded helix-helix interactions

The Vpu protein of HIV-1 antagonizes BST-2 (tetherin), a broad spectrum effector of the innate immune response to viral infection, by an intermolecular interaction that maps genetically to the α-helical transmembrane domains (TMDs) of each protein. Here we utilize NMR spectroscopy to describe key features of the helix-helix pairing that underlies this interaction. The antagonism of BST-2 involves a sequence of three alanines and a tryptophan spaced at four residue intervals within the Vpu TMD helix. Responsiveness to Vpu involves bulky hydrophobic residues in the C-terminal region of the BST-2 TMD helix that likely fit between the alanines on the interactive face of Vpu. These aspects of Vpu and BST-2 form an anti-parallel, lipid-embedded helix-helix interface. Changes in human BST-2 that mimic sequences found in nonhuman primate orthologs unresponsive to Vpu change the tilt angle of the TMD in the lipid bilayer without abrogating its intrinsic ability to interact with Vpu. These data explain the mechanism by which HIV-1 evades a key aspect of innate immunity and the species specificity of Vpu using an anti-parallel helix-helix packing model.     

The pinholin of lambdoid phage 21: control of lysis by membrane depolarization

The phage 21 holin, S²¹, forms small membrane holes that depolarize the membrane and is designated as a pinholin, as opposed to large-hole-forming holins, like S^λ. Pinholins require secreted SAR endolysins, a pairing that may represent an intermediate in the evolution of canonical holin-endolysin systems.     

Effects of amantadine on the dynamics of membrane-bound influenza A M2 transmembrane peptide studied by NMR relaxation

The molecular motions of membrane proteins in liquid-crystalline lipid bilayers lie at the interface between motions in isotropic liquids and in solids. Specifically, membrane proteins can undergo whole-body uniaxial diffusion on the microsecond time scale. In this work, we investigate the ¹H rotating-frame spin-lattice relaxation (T _1ρ) caused by the uniaxial diffusion of the influenza A M2 transmembrane peptide (M2TMP), which forms a tetrameric proton channel in lipid bilayers. This uniaxial diffusion was proved before by ²H, ¹⁵N and ¹³C NMR lineshapes of M2TMP in DLPC bilayers. When bound to an inhibitor, amantadine, the protein exhibits significantly narrower linewidths at physiological temperature. We now investigate the origin of this line narrowing through temperature-dependent ¹H T _1ρ relaxation times in the absence and presence of amantadine. Analysis of the temperature dependence indicates that amantadine decreases the correlation time of motion from 2.8 ± 0.9 μs for the apo peptide to 0.89 ± 0.41 μs for the bound peptide at 313 K. Thus the line narrowing of the bound peptide is due to better avoidance of the NMR time scale and suppression of intermediate time scale broadening. The faster diffusion of the bound peptide is due to the higher attempt rate of motion, suggesting that amantadine creates better-packed and more cohesive helical bundles. Analysis of the temperature dependence of $ { \ln }\left( {T_{1\rho }^{ - 1} } \right) $ indicates that the activation energy of motion increased from 14.0 ± 4.0 kJ/mol for the apo peptide to 23.3 ± 6.2 kJ/mol for the bound peptide. This higher activation energy indicates that excess amantadine outside the protein channel in the lipid bilayer increases the membrane viscosity. Thus, the protein-bound amantadine speeds up the diffusion of the helical bundles while the excess amantadine in the bilayer increases the membrane viscosity.     

Molecular dynamics simulation of Bombolitin II in the dipalmitoylphosphatidylcholine membrane bilayer

The orientation behavior of Bombolitin II (BLT2) in the dipalmitoylphosphatidylcholine membrane bilayer was investigated by using molecular-dynamics simulation. During the 20-ns simulation, the BLT2 began to tilt and finally reached the angle of 51° from the membrane-normal. The structure of the peptide formed the amphipathic α-helical structure during the entire simulation time. The peptide tilts with its hydrophobic side faced to the hydrophobic core of the bilayer. We analyzed the mechanism of the tilting behavior of the peptide associated with the membrane in detail. The analysis showed that the hydrogen-bond interaction and the electrostatic interaction were found to exist between Lys¹² and a lipid molecule. These interactions are considered to work as an important factor in tilting the peptide to the membrane-normal.     

High-Resolution Orientation and Depth of Insertion of the Voltage-Sensing S4 Helix of a Potassium Channel in Lipid Bilayers

The opening and closing of voltage-gated potassium (Kv) channels are controlled by several conserved Arg residues in the S4 helix of the voltage-sensing domain. The interaction of these positively charged Arg residues with the lipid membrane has been of intense interest for understanding how membrane proteins fold to allow charged residues to insert into lipid bilayers against free-energy barriers. Using solid-state NMR, we have now determined the orientation and insertion depth of the S4 peptide of the KvAP channel in lipid bilayers. Two-dimensional ¹⁵N correlation experiments of macroscopically oriented S4 peptide in phospholipid bilayers revealed a tilt angle of 40° and two possible rotation angles differing by 180° around the helix axis. Remarkably, the tilt angle and one of the two rotation angles are identical to those of the S4 helix in the intact voltage-sensing domain, suggesting that interactions between the S4 segment and other helices of the voltage-sensing domain are not essential for the membrane topology of the S4 helix. ¹³C-³¹P distances between the S4 backbone and the lipid ³¹P indicate a ∼ 9 Å local thinning and 2 Å average thinning of the DMPC (1,2-dimyristoyl-sn-glycero-3-phosphochloline)/DMPG (1,2-dimyristoyl-sn-glycero-3-phosphatidylglycerol) bilayer, consistent with neutron diffraction data. Moreover, a short distance of 4.6 Å from the guanidinium C^ζ of the second Arg to ³¹P indicates the existence of guanidinium phosphate hydrogen bonding and salt bridges. These data suggest that the structure of the Kv gating helix is mainly determined by protein-lipid interactions instead of interhelical protein-protein interactions, and the S4 amino acid sequence encodes sufficient information for the membrane topology of this crucial gating helix.     

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.