Charged or aromatic anchor residue dependence of transmembrane peptide tilt

Classifying tetraspanins: A universal system for numbering residues and a proposal for naming structural motifs and subfamilies

All tetraspanins have four transmembrane domains (TMs). The large extracellular loop (LEL) that connects the third and fourth TMs contains multiple secondary structures together with the family's signature Cys-Cys-Gly motif. These intriguing membrane proteins are involved in diverse and incompletely understood cellular processes including cell adhesion, tissue differentiation, immune cell maturation and host-parasite interactions. Here we present a classification system that accurately describes the position of each amino acid within its primary sequence based on both sequence and topological conservation of the TMs and LEL. This builds on the numbering systems that have been used in the G protein-coupled receptor (GPCR) field for nearly three decades and which have aided the understanding of GPCR structure/activity relationships and ligand interactions. The high-resolution structures of the tetraspanins CD81, CD9, CD53 and Tspan15 were used to validate the structural relevance of our new tetraspanin classification system. Modelling of all tetraspanin LELs highlighted flexibility in LEL disulfide bonding across the family and suggests that the structural arrangement of tetraspanin LELs is more complex than previously thought. We therefore propose a new subfamily naming system that addresses this added complexity and facilitates the systematic classification of human tetraspanins, shedding light on all structural motifs within the family. We anticipate that our universal tetraspanin classification system will enable progress in defining how sequence and structure inform function.

p85S6K sustains synaptic GluA1 to ameliorate cognitive deficits in Alzheimer's disease

BACKGROUND

Ribosomal protein S6 kinase 1 (S6K1) is a serine-threonine kinase that has two main isoforms: p70S6K (70-kDa isoform) and p85S6K (85-kDa isoform). p70S6K, with its upstream mammalian target of rapamycin (mTOR), has been shown to be involved in learning and memory and participate in the pathophysiology of Alzheimer's disease (AD). However, the function of p85S6K has long been neglected due to its high similarity to p70S6k. The role of p85S6K in learning and memory is still largely unknown.

METHODS

We fractionated the postsynaptic densities to illustrate the differential distribution of p85S6K and p70S6K. Coimmunoprecipitation was performed to unveil interactions between p85S6K and the GluA1 subunit of AMPA receptor. The roles of p85S6K in synaptic targeting of GluA1 and learning and memory were evaluated by specific knockdown or overexpression of p85S6K followed by a broad range of methodologies including immunofluorescence, Western blot, in situ proximity ligation assay, morphological staining and behavioral examination. Further, the expression level of p85S6K was measured in brains from AD patients and AD model mice.

RESULTS

p85S6K, but not p70S6K, was enriched in the postsynaptic densities. Moreover, knockdown of p85S6K resulted in defective spatial and recognition memory. In addition, p85S6K could interact with the GluA1 subunit of AMPA receptor through synapse-associated protein 97 and A-kinase anchoring protein 79/150. Mechanistic studies demonstrated that p85S6K could directly phosphorylate GluA1 at Ser845 and increase the amount of GluA1 in synapses, thus sustaining synaptic function and spine densities. Moreover, p85S6K was found to be specifically decreased in the synaptosomal compartment in the brains of AD patients and AD mice. Overexpression of p85S6K ameliorated the synaptic deficits and cognitive impairment in transgenic AD model mice.

CONCLUSIONS

These results strongly imply a significant role for p85S6K in maintaining synaptic and cognitive function by interacting with GluA1. The findings provide an insight into the rational targeting of p85S6K as a therapeutic potential for AD.

Solid-State NMR: Methods for Biological Solids

In the last two decades, solid-state nuclear magnetic resonance (ssNMR) spectroscopy has transformed from a spectroscopic technique investigating small molecules and industrial polymers to a potent tool decrypting structure and underlying dynamics of complex biological systems, such as membrane proteins, fibrils, and assemblies, in near-physiological environments and temperatures. This transformation can be ascribed to improvements in hardware design, sample preparation, pulsed methods, isotope labeling strategies, resolution, and sensitivity. The fundamental engagement between nuclear spins and radio-frequency pulses in the presence of a strong static magnetic field is identical between solution and ssNMR, but the experimental procedures vastly differ because of the absence of molecular tumbling in solids. This review discusses routinely employed state-of-the-art static and MAS pulsed NMR methods relevant for biological samples with rotational correlation times exceeding 100's of nanoseconds. Recent developments in signal filtering approaches, proton methodologies, and multiple acquisition techniques to boost sensitivity and speed up data acquisition at fast MAS are also discussed. Several examples of protein structures (globular, membrane, fibrils, and assemblies) solved with ssNMR spectroscopy have been considered. We also discuss integrated approaches to structurally characterize challenging biological systems and some newly emanating subdisciplines in ssNMR spectroscopy.

Membrane electrostatics sensed by tryptophan anchors in hydrophobic model peptides depends on non-aromatic interfacial amino acids: implications in hydrophobic mismatch

WALPs are synthetic α-helical membrane-spanning peptides that constitute a well-studied system for exploring hydrophobic mismatch. These peptides represent a simplified consensus motif for transmembrane domains of intrinsic membrane proteins due to their hydrophobic core of alternating leucine and alanine flanked by membrane-anchoring aromatic tryptophan residues. Although the modulation of mismatch responses in WALPs by tryptophan anchors has been reported earlier, there have been limited attempts to utilize the intrinsic tryptophan fluorescence of this class of peptides in mismatch sensors. We have previously shown, utilizing the red edge excitation shift (REES) approach, that interfacial WALP tryptophan residues in fluid phase bilayers experience a dynamically constrained membrane microenvironment. Interestingly, emerging reports suggest the involvement of non-aromatic interfacially localized residues in modulating local structure and dynamics in WALP analogs. In this backdrop, we have explored the effect of interfacial amino acids, such as lysine (in KWALPs) and glycine (in GWALPs), on the tryptophan microenvironment of WALP analogs in zwitterionic and negatively charged membranes. We show that interfacial tryptophans in KWALP and GWALP experience a more restricted microenvironment, as reflected in the substantial increase in magnitude of REES and apparent rotational correlation time, relative to those in WALP in zwitterionic membranes. Interestingly, in contrast to WALP, the tryptophan anchors in KWALP and GWALP appear insensitive to the presence of negatively charged lipids in the membrane. These results reveal a subtle interplay between non-aromatic flanking residues in transmembrane helices and negatively charged lipids at the membrane interface, which could modulate the membrane microenvironment experienced by interfacially localized tryptophan residues. Since interfacial tryptophans are known to influence mismatch responses in WALPs, our results highlight the possibility of utilizing the fluorescence signatures of tryptophans in membrane proteins or model peptides such as WALP as markers for assessing protein responses to hydrophobic mismatch. More importantly, these results constitute one of the first reports on the influence of lipid headgroup charge in fine-tuning hydrophobic mismatch in membrane bilayers, thereby enriching the existing framework of hydrophobic mismatch.

Illuminating Disorder Induced by Glu in a Stable Arg-Anchored Transmembrane Helix

Membrane proteins are vital for biological function and are complex to study. Even in model peptide-lipid systems, the combined influence or interaction of pairs of chemical groups still is not well understood. Disordered proteins, whether in solution or near lipid membranes, are an emerging paradigm for the initiation and control of biological function. The disorder can involve molecular orientation as well as molecular folding. This paper reports an astonishing induction of disorder when one Glu residue is introduced into a highly stable 23-residue transmembrane helix. The parent helix is anchored by a single Arg residue, tilted at a well-defined angle with respect to the DOPC bilayer normal and undergoes rapid cone precession. When Glu is introduced two residues away from Arg, with 200° (or 160°) radial separation, the helix properties change radically to exhibit a multiplicity of three or more disordered states. The helix characteristics have been monitored by deuterium (²H) NMR spectroscopy as functions of the pH and lipid bilayer composition. The disordered multistate behavior of the (Glu, Arg)-containing helix varies with the lipid bilayer thickness and pH. The results highlight a fundamental induction of protein multistate properties by a single Glu residue in a lipid membrane environment.

Lipid-Dependent Titration of Glutamic Acid at a Bilayer Membrane Interface

The ionization properties of protein side chains in lipid-bilayer membranes will differ from the canonical values of side chains exposed to an aqueous solution. While the propensities of positively charged side chains of His, Lys, and Arg to release a proton in lipid membranes have been rather well characterized, the propensity for a negatively charged Glu side chain to receive a proton and achieve the neutral state in a bilayer membrane has been less well characterized. Indeed, the ionization of the glutamic acid side chain has been predicted to depend on its depth of burial in a lipid membrane but has been difficult to verify experimentally. To address the issue, we incorporated an interfacial Glu residue at position 4 of a distinct 23-residue transmembrane helix and used ²H NMR to examine the helix properties as a function of pH. We observe that the helix tilt and azimuthal rotation vary little with pH, but the extent of helix unraveling near residues 3 and 4 changes as the Glu residue E4 titrates. Remarkably, the ²H quadrupolar splitting for the side chain of alanine A3 responds to pH with an apparent pK _a of 4.8 in 1,2-dilauroyl-sn-glycero-3-phosphocholine (DLPC) and 6.3 in 1,2-dimyristoyl-sn-glycero-3-phosphatidylcholine (DMPC), but is unchanged up to pH 8.0 in 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) in the presence of residue E4. With bilayers composed of alkali-stable ether-linked lipids, the side chain of A3 responds to pH with an apparent pK _a of 11.0 in the ether analogue of DOPC. These results suggest that the depth dependence of Glu ionization in lipid-bilayer membranes may be steeper than previously predicted or envisioned.

Examination of pH dependency and orientation differences of membrane spanning alpha helices carrying a single or pair of buried histidine residues

We have employed the peptide framework of GWALP23 (acetyl-GGALWLALALALALALALWLAGA-amide) to examine the orientation, dynamics and pH dependence of peptides having buried single or pairs of histidine residues. When residue L8 is substituted to yield GWALP23-H⁸, acetyl-GGALWLAH⁸ALALALALALWLAGA-amide, the deuterium NMR spectra of ²H-labeled core alanine residues reveal a helix that occupies a single transmembrane orientation in DLPC, or in DMPC at low pH, yet shows multiple states at higher pH or in bilayers of DOPC. Moreover, a single histidine at position 8 or 16 in the GWALP23 framework is sensitive to pH. Titration points are observed near pH 3.5 for the deprotonation of H8 in lipid bilayers of DLPC or DMPC, and for H16 in DOPC. When residues L8 and L16 both are substituted to yield GWALP23-H^8,16, the ²H NMR spectra show, interestingly, no titration dependence from pH 2-8, yet bilayer thickness-dependent orientation differences. The helix with H8 and H16 is found to adopt a transmembrane orientation in thin bilayers of DLPC, a combination of transmembrane and surface orientations in DMPC, and then a complete transition to a surface bound orientation in the thicker DPoPC and DOPC lipid bilayers. In the surface orientations, alanine A7 no longer fits within the core helix. These results along with previous studies with different locations of histidine residues suggest that lipid hydrophobic thickness is a first determinant and pH a second determinant for the helical orientation, along with possible side-chain snorkeling, when the His residues are incorporated into the hydrophobic region of a lipid membrane-associated helix.

Comparing Interfacial Trp, Interfacial His and pH Dependence for the Anchoring of Tilted Transmembrane Helical Peptides

Charged and aromatic amino acid residues, being enriched toward the terminals of membrane-spanning helices in membrane proteins, help to stabilize particular transmembrane orientations. Among them, histidine is aromatic and can be positively charge at low pH. To enable investigations of the underlying protein-lipid interactions, we have examined the effects of single or pairs of interfacial histidine residues using the constructive low-dynamic GWALP23 (acetyl-GG²ALW⁵LALALALALALALW¹⁹LAG²²A-amide) peptide framework by incorporating individual or paired histidines at locations 2, 5, 19 or 22. Analysis of helix orientation by means of solid-state ²H NMR spectra of labeled alanine residues reveals marked differences with H^2,22 compared to W^2,22. Nevertheless, the properties of membrane-spanning H^2,22WALP23 helices show little pH dependence and are similar to those having Gly, Arg or Lys at positions 2 and 22. The presence of H5 or H19 influences the helix rotational preference but not the tilt magnitude. H5 affects the helical integrity, as residue 7 unwinds from the core helix; yet once again the helix orientation and dynamic properties show little sensitivity to pH. The overall results reveal that the detailed properties of transmembrane helices depend upon the precise locations of interfacial histidine residues.

Influence of interfacial tryptophan residues on an arginine-flanked transmembrane helix

The transmembrane helices of membrane proteins often are flanked by interfacial charged or aromatic residues that potentially help to anchor the membrane-spanning protein. For isolated single-span helices, the interfacial residues may be especially important for stabilizing particular tilted transmembrane orientations. The peptide RWALP23 (acetyl-GR²ALW(LA)₆LWLAR²²A-amide) has been employed to investigate the interplay between interfacial arginines and tryptophans. Here we replace the tryptophans of RWALP23 with A5 and A19, to investigate arginines alone with respect to helix fraying and orientation in varying lipid bilayers. Deuterated alanines incorporated into the central sequence allow the orientation and stability of the core helix to be assessed by means of solid -state ²H NMR in bilayers of DOPC, DMPC and DLPC. The helix tilt from the bilayer normal is found to increase slightly when R2 and R22 are present, and increases still further when the tryptophans W5 and W19 are replaced by alanines. The extent of helix dynamic averaging remains low in all cases. The preferred helix azimuthal rotation is essentially constant for all of the helices in each of the lipid membranes considered here. The alanines located outside of the core region of the peptide are sensitive to helical integrity. The new alanines, A5 and A19, therefore, provide new information about the length of the core helix and the onset of unraveling of the terminals. Residue A19 remains essentially on the central helix in each lipid membrane, while residues A3, A5 and A21 deviate from the core helix to an extent that depends on the membrane thickness. Differential unraveling of the two ends to expose peptide backbone groups for hydrogen bonding therefore acts together with specific interfacial side chains to stabilize a transmembrane helix.

Influence of Lipid Saturation, Hydrophobic Length and Cholesterol on Double-Arginine-Containing Helical Peptides in Bilayer Membranes

Membrane proteins are essential for many cell processes yet are more difficult to investigate than soluble proteins. Charged residues often contribute significantly to membrane protein function. Model peptides such as GWALP23 (acetyl-GGALW⁵ LAL⁸ LALALAL¹⁶ ALW¹⁹ LAGA-amide) can be used to characterize the influence of specific residues on transmembrane protein domains. We have substituted R8 and R16 in GWALP23 in place of L8 and L16, equidistant from the peptide center, and incorporated specific ² H-labeled alanine residues within the central sequence for detection by solid-state ² H NMR spectroscopy. The resulting pattern of [² H]Ala quadrupolar splitting (Δν_q ) magnitudes indicates the core helix for R^8,16 GWALP23 is significantly tilted to give a similar transmembrane orientation in thinner bilayers with either saturated C12:0 or C14:0 acyl chains (1,2-dilauroyl-sn-glycero-3-phosphocholine (DLPC) or 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC)) or unsaturated C16:1 Δ9 cis acyl chains. In bilayers of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC; C18:1 Δ9 cis) multiple orientations are indicated, whereas in longer, unsaturated 1,2-dieicosenoyl-sn-glycero-3-phosphocholine (DEiPC; C20:1 Δ11 cis) bilayers, the R^8,16 GWALP23 helix adopts primarily a surface orientation. The inclusion of 10-20 mol % cholesterol in DOPC bilayers drives more of the R^8,16 GWALP23 helix population to the membrane surface, thereby allowing both charged arginines access to the interfacial lipid head groups. The results suggest that hydrophobic thickness and cholesterol content are more important than lipid saturation for the arginine peptide dynamics and helix orientation in lipid membranes.

Breaking the Backbone: Central Arginine Residues Induce Membrane Exit and Helix Distortions within a Dynamic Membrane Peptide

Transmembrane domains of membrane proteins sometimes contain conserved charged or ionizable residues which may be essential for protein function and regulation. This work examines the molecular interactions of single Arg residues within a highly dynamic transmembrane peptide helix. To this end, we have modified the GW^4,20ALP23 (acetyl-GGAW⁴(AL)₇AW²⁰AGA-amide) model peptide framework to incorporate Arg residues near the center of the peptide. Peptide helix formation, orientation and dynamics were analyzed by means of solid-state NMR spectroscopy to monitor specific ²H- or ¹⁵N-labeled residues. GW^4,20ALP23 itself adopts a tilted orientation within lipid bilayer membranes. Nevertheless, the GW^4,20ALP23 helix exhibits moderate to high dynamic averaging of NMR observables, such as ²H quadrupolar splittings or ¹⁵N-¹H dipolar couplings, due to competition between the interfacial Trp residues on opposing helix faces. Here we examine how the helix dynamics are impacted by the introduction of a single Arg residue at position 12 or 14. Residue R14 restricts the helix to low dynamic averaging and a well-defined tilt that varies inversely with the lipid bilayer thickness. To compensate for the dominance of R14, the competing Trp residues cause partial unwinding of the helix at the C-terminal. By contrast, R¹²GW^4,20ALP23 exits the DOPC bilayer to an interfacial surface-bound location. Interestingly, multiple orientations are exhibited by a single residue, Ala-9. Quadrupolar splittings generated by ²H-labeled residues A3, A5, A7, and A9 do not fit to the α-helical quadrupolar wave plot defined by residues A11, A13, A15, A17, A19, and A21. The discontinuity at residue A9 implicates a helical swivel distortion and an apparent 3₁₀-helix involving the N-terminal residues preceding A11. These molecular features suggest that, while arginine residues are prominent factors controlling transmembrane helix dynamics, the influence of interfacial tryptophan residues cannot be ignored.

Membrane Bending Moduli of Coexisting Liquid Phases Containing Transmembrane Peptide

A number of highly curved membranes in vivo, such as epithelial cell microvilli, have the relatively high sphingolipid content associated with "raft-like" composition. Given the much lower bending energy measured for bilayers with "nonraft" low sphingomyelin and low cholesterol content, observing high curvature for presumably more rigid compositions seems counterintuitive. To understand this behavior, we measured membrane rigidity by fluctuation analysis of giant unilamellar vesicles. We found that including a transmembrane helical GWALP peptide increases the membrane bending modulus of the liquid-disordered (Ld) phase. We observed this increase at both low-cholesterol fraction and higher, more physiological cholesterol fraction. We find that simplified, commonly used Ld and liquid-ordered (Lo) phases are not representative of those that coexist. When Ld and Lo phases coexist, GWALP peptide favors the Ld phase with a partition coefficient of 3-10 depending on mixture composition. In model membranes at high cholesterol fractions, Ld phases with GWALP have greater bending moduli than the Lo phase that would coexist.

Transmembrane Helix Integrity versus Fraying To Expose Hydrogen Bonds at a Membrane-Water Interface

Transmembrane helices dominate the landscape for many membrane proteins. Often flanked by interfacial aromatic residues, these transmembrane helices also contain loops and interhelix segments, which could help in stabilizing a transmembrane orientation. Using ²H nuclear magnetic resonance spectroscopy to monitor bilayer-incorporated model GWALP23 family peptides, we address systematically the issue of helix fraying in relation to the dynamics and orientation of highly similar individual transmembrane helices. We inserted aromatic (Phe, Trp, Tyr, and His) or non-aromatic residues (Ala and Gly) into positions 4 and 5 adjacent to a core transmembrane helix to examine the side-chain dependency of the transmembrane orientation, dynamics, and helix integrity (extent and location of unraveling). Incorporation of [²H]alanine labels enables one to assess the helicity of the core sequence and the peptide termini. For most of the helices, we observed substantial unwinding involving at least three residues at both ends. For the unique case of histidine at positions 4 and 5, an extended N-terminal unwinding was observed up to residue 7. For further investigation of the onset of fraying, we employed A^4,5GWALP23 with ²H labels at residues 4 and 5 and found that the number of terminal residues involved in the unwinding depends on bilayer thicknesses and helps to govern the helix dynamics. The combined results enable us to compare and contrast the extent of fraying for each related helix, as reflected by the deviation of experimental ²H quadrupolar splitting magnitudes of juxta-terminal alanines A3 and A21 from those represented by an ideal helix geometry.

Actionable Activating Oncogenic ERBB2/HER2 Transmembrane and Juxtamembrane Domain Mutations

Deregulated HER2 is a target of many approved cancer drugs. We analyzed 111,176 patient tumors and identified recurrent mutations in HER2 transmembrane domain (TMD) and juxtamembrane domain (JMD) that include G660D, R678Q, E693K, and Q709L. Using a saturation mutagenesis screen and testing of patient-derived mutations we found several activating TMD and JMD mutations. Structural modeling and analysis showed that the TMD/JMD mutations function by improving the active dimer interface or stabilizing an activating conformation. Further, we found that HER2 G660D employed asymmetric kinase dimerization for activation and signaling. Importantly, anti-HER2 antibodies and small-molecule kinase inhibitors blocked the activity of TMD/JMD mutants. Consistent with this, a G660D germline mutant lung cancer patient showed remarkable clinical response to HER2 blockade.

Wavelength-Selective Fluorescence of a Model Transmembrane Peptide: Constrained Dynamics of Interfacial Tryptophan Anchors

WALPs are prototypical, α-helical transmembrane peptides that represent a consensus sequence for transmembrane segments of integral membrane proteins and serve as excellent models for exploring peptide-lipid interactions and hydrophobic mismatch in membranes. Importantly, the WALP peptides are in direct contact with the lipids. They consist of a central stretch of alternating hydrophobic alanine and leucine residues capped at both ends by tryptophans. In this work, we employ wavelength-selective fluorescence approaches to explore the intrinsic fluorescence of tryptophan residues in WALP23 in 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) membranes. Our results show that the four tryptophan residues in WALP23 exhibit an average red edge excitation shift (REES) of 6 nm, implying their localization at the membrane interface, characterized by a restricted microenvironment. This result is supported by fluorescence anisotropy and lifetime measurements as a function of wavelength displayed by WALP23 tryptophans in POPC membranes. These results provide a new approach based on intrinsic fluorescence of interfacial tryptophans to address protein-lipid interaction and hydrophobic mismatch.

Control of Transmembrane Helix Dynamics by Interfacial Tryptophan Residues

Transmembrane protein domains often contain interfacial aromatic residues, which may play a role in the insertion and stability of membrane helices. Residues such as Trp or Tyr, therefore, are often found situated at the lipid-water interface. We have examined the extent to which the precise radial locations of interfacial Trp residues may influence peptide helix orientation and dynamics. To address these questions, we have modified the GW5,19ALP23 (acetyl-GGALW5(LA)6LW19LAGA-[ethanol]amide) model peptide framework to relocate the Trp residues. Peptide orientation and dynamics were analyzed by means of solid-state nuclear magnetic resonance (NMR) spectroscopy to monitor specific 2H- and 15N-labeled residues. GW5,19ALP23 adopts a defined, tilted orientation within lipid bilayer membranes with minimal evidence of motional averaging of NMR observables, such as 2H quadrupolar or 15N-1H dipolar splittings. Here, we examine how peptide dynamics are impacted by relocating the interfacial Trp (W) residues on both ends and opposing faces of the helix, for example by a 100° rotation on the helical wheel for positions 4 and 20. In contrast to GW5,19ALP23, the modified GW4,20ALP23 helix experiences more extensive motional averaging of the NMR observables in several lipid bilayers of different thickness. Individual and combined Gaussian analyses of the 2H and 15N NMR signals confirm that the extent of dynamic averaging, particularly rotational "slippage" about the helix axis, is strongly coupled to the radial distribution of the interfacial Trp residues as well as the bilayer thickness. Additional 2H labels on alanines A3 and A21 reveal partial fraying of the helix ends. Even within the context of partial unwinding, the locations of particular Trp residues around the helix axis are prominent factors for determining transmembrane helix orientation and dynamics within the lipid membrane environment.

Membrane-Induced p Ka Shifts in wt-pHLIP and Its L16H Variant

The pH (low) insertion peptides (pHLIPs) is a family of peptides that are able to insert into a lipid bilayer at acidic pH. The molecular mechanism of pHLIPs insertion, folding, and stability in the membrane at low pH is based on multiple protonation events, which are challenging to study at the molecular level. More specifically, the relation between the experimental p K of insertion (p K^exp) of pHLIPs and the p K_a of the key residues is yet to be clarified. We carried out a computational study, complemented with new experimental data, and established the influence of (de)protonation of titrable residues on the stability of the peptide membrane-inserted state. Constant-pH molecular dynamics simulations were employed to calculate the p K_a values of these residues along the membrane normal. In the wt-pHLIP, we identified Asp14 as the key residue for the stability of the membrane-inserted state, and its p K_a value is strongly correlated with the experimental p K^exp measured in thermodynamics studies. Also, in order to narrow down the pH range at which pHLIP is stable in the membrane, we designed a new pHLIP variant, L16H, where Leu in the 16th position was replaced by a titrable His residue. Our results showed that the L16H variant undergoes two transitions. The calculated p K_a and experimentally observed p K^exp values are in good agreement. Two distinct p K^exp values delimit a pH range where the L16H peptide is stably inserted in the membrane, while, outside this range, the membrane-inserted state is destabilized and the peptide exits from the bilayer. pHLIP peptides have been successfully used to target cancer cells for the delivery of diagnostics and therapeutic agents to acidic tumors. The fine-tuning of the stability of the pHLIP inserted state and its restriction to a narrow well-defined pH range might allow the design of new peptides, able to discriminate between tissues with different extracellular pH values.

Helix formation and stability in membranes

In this article we review current understanding of basic principles for the folding of membrane proteins, focusing on the more abundant alpha-helical class. Membrane proteins, vital to many biological functions and implicated in numerous diseases, fold into their active conformations in the complex environment of the cell bilayer membrane. While many membrane proteins rely on the translocon and chaperone proteins to fold correctly, others can achieve their functional form in the absence of any translation apparatus or other aides. Nevertheless, the spontaneous folding process is not well understood at the molecular level. Recent findings suggest that helix fraying and loop formation may be important for overall structure, dynamics and regulation of function. Several types of membrane helices with ionizable amino acids change their topology with pH. Additionally we note that some peptides, including many that are rich in arginine, and a particular analogue of gramicidin, are able passively to translocate across cell membranes. The findings indicate that a final protein structure in a lipid-bilayer membrane is sequence-based, with lipids contributing to stability and regulation. While much progress has been made toward understanding the folding process for alpha-helical membrane proteins, it remains a work in progress. This article is part of a Special Issue entitled: Emergence of Complex Behavior in Biomembranes edited by Marjorie Longo.

Influence of glutamic acid residues and pH on the properties of transmembrane helices

Negatively charged side chains are important for the function of particular ion channels and certain other membrane proteins. To investigate the influence of single glutamic acid side chains on helices that span lipid-bilayer membranes, we have employed GWALP23 (acetyl-GGALW⁵LALALALALALALW¹⁹LAGA-amide) as a favorable host peptide framework. We substituted individual Leu residues with Glu residues (L12E or L14E or L16E) and incorporated specific ²H-labeled alanine residues within the core helical region or near the ends of the sequence. Solid-state ²H NMR spectra reveal little change for the core labels in GWALP23-E12, -E14 and -E16 over a pH range of 4 to 12.5, with the spectra being broader for samples in DOPC compared to DLPC bilayers. The spectra for samples with deuterium labels near the helix ends on alanines 3 and 21 show modest pH-dependent changes in the extent of unwinding of the helix terminals in DLPC and DOPC bilayers. The combined results indicate minor overall responses of these transmembrane helices to changes in pH, with the most buried residue E12 showing no pH dependence. While the Glu residues E14 and E16 may have high pK_a values in the lipid bilayer environment, it is also possible that a paucity of helix response is masking the pK_a values. Interestingly, when E16 is present, spectral changes at high pH report significant local unwinding of the core helix. Our results are consistent with the expectation that buried carboxyl groups aggressively hold their protons and/or waters of hydration.

Influence of High pH and Cholesterol on Single Arginine-Containing Transmembrane Peptide Helices

An essential component of mammalian cells, cholesterol exerts significant influence on the physical properties of the cell membrane and in turn its constituents, including membrane proteins. Although sparse, polar amino acid residues are highly conserved in membrane proteins and play pivotal roles in determining specific structural and functional properties. To improve our understanding of particular polar residues in the membrane environment, we have examined two specific “guest” Arg residues within a well-defined and deuterium-labeled “host” framework provided by the transmembrane helical peptide GWALP23 (acetyl-GGALWLALALALALALALWLAGA-amide). Solid-state ²H nuclear magnetic resonance (NMR) spectra from aligned bilayer membrane samples effectively report changes in the host helix properties because of the incorporation of the guest residues. The focus of this work is two-pronged. First, GWALP23-R14 was examined over a pH range of 2–13 in 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) ester- or ether-linked bilayer membranes. Our findings indicate that the Arg guanidinium side chain remains charged over this entire range, in agreement with numerous molecular dynamics simulations. Second, GWALP23-R12 and GWALP23-R14 peptides were characterized in DOPC bilayers with varying cholesterol content. Our findings suggest that 10 or 20% cholesterol content has minimal impact on the orientation of the R14 peptide. Although the NMR signals are broader and weaker in the presence of 20% cholesterol, the deuterium quadrupolar splittings for [²H]Ala residues in GWALP23-R14 change very little. Conversely, cholesterol appears to modulate the multistate behavior of GWALP23-R12 and to favor a major interfacial state for the helix, bound at the bilayer surface. These results indicate a conditional sensitivity of a complex multistate transmembrane Arg-containing peptide helix to the presence of cholesterol.

Impairment on the lateral mobility induced by structural changes underlies the functional deficiency of the lupus-associated polymorphism FcγRIIB-T232

Xu et al. show that the lupus-associated polymorphism FcγRIIB-T232 has structural changes of the TM domain that reduces lateral mobility and inhibitory functions.

FcγRIIB functions to suppress the activation of immune cells. A single-nucleotide polymorphism in the transmembrane (TM) domain of FcγRIIB, FcγRIIB-T232, is associated with lupus. In this study, we investigated the pathogenic mechanism of FcγRIIB-T232 at both functional and structural levels. Our results showed that FcγRIIB-T232 exhibited significantly reduced lateral mobility compared with FcγRIIB-I232 and was significantly less enriched into the microclusters of immune complexes (ICs) after stimulation. However, if sufficient responding time is given for FcγRIIB-T232 to diffuse and interact with the ICs, FcγRIIB-T232 can restore its inhibitory function. Moreover, substituting the FcγRIIB-T232 TM domain with that of a fast floating CD86 molecule restored both the rapid mobility and the inhibitory function, which further corroborated the importance of fast mobility for FcγRIIB to function. Mechanistically, the crippled lateral mobility of FcγRIIB-T232 can be explained by the structural changes of the TM domain. Both atomistic simulations and nuclear magnetic resonance measurement indicated that the TM helix of FcγRIIB-T232 exhibited a more inclined orientation than that of FcγRIIB-I232, thus resulting in a longer region embedded in the membrane. Therefore, we conclude that the single-residue polymorphism T232 enforces the inclination of the TM domain and thereby reduces the lateral mobility and inhibitory functions of FcγRIIB.

Ionization Properties of Histidine Residues in the Lipid Bilayer Membrane Environment

We address the critically important ionization properties of histidine side chains of membrane proteins, when exposed directly to lipid acyl chains within lipid bilayer membranes. The problem is important for addressing general principles that may underlie membrane protein function. To this end, we have employed a favorable host peptide framework provided by GWALP23 (acetyl-GGALW(5)LALALALALALALW(19)LAGA-amide). We inserted His residues into position 12 or 14 of GWALP23 (replacing either Leu(12) or Leu(14)) and incorporated specific [(2)H]Ala labels within the helical core sequence. Solid-state (2)H NMR spectra report the folding and orientation of the core sequence, revealing marked differences in the histidine-containing transmembrane helix behavior between acidic and neutral pH conditions. At neutral pH, the GWALP23-H12 and GWALP23-H14 helices exhibit well defined tilted transmembrane orientations in dioleoylphosphatidylcholine (DOPC)and dilauroylphosphatidylcholine (DLPC) bilayer membranes. Under acidic conditions, when His(12) is protonated and charged, the GWALP23-H12 helix exhibits a major population that moves to the DOPC bilayer surface and a minor population that occupies multiple transmembrane states. The response to protonation of His(14) is an increase in helix tilt, but GWALP23-H14 remains in a transmembrane orientation. The results suggest pKa values of less than 3 for His(12) and about 3-5 for His(14) in DOPC membranes. In the thinner DLPC bilayers, with increased water access, the helices are less responsive to changes in pH. The combined results enable us to compare the ionization properties of lipid-exposed His, Lys, and Arg side chains in lipid bilayer membranes.

Juxta-terminal Helix Unwinding as a Stabilizing Factor to Modulate the Dynamics of Transmembrane Helices

Transmembrane helices of integral membrane proteins often are flanked by interfacial aromatic residues that can serve as anchors to aid the stabilization of a tilted transmembrane orientation. Yet, physical factors that govern the orientation or dynamic averaging of individual transmembrane helices are not well understood and have not been adequately explained. Using solid-state (2) H NMR spectroscopy to examine lipid bilayer-incorporated model peptides of the GWALP23 (acetyl-GGALW(LA)6 LWLAGA-amide) family, we observed substantial unwinding at the terminals of several tilted helices spanning the membranes of DLPC, DMPC, or DOPC lipid bilayers. The fraying of helix ends might be vital for defining the dynamics and orientations of transmembrane helices in lipid bilayer membranes.

C≡N stretching vibration of 5-cyanotryptophan as an infrared probe of protein local environment: what determines its frequency?

Recently it has been suggested that the C≡N stretching vibration of a tryptophan analog, 5-cyanotryptophan, could be used as an infrared probe of the local environment, especially the hydration status, of tryptophan residues in proteins. However, the factors that influence the frequency of this vibrational mode are not understood. To determine these factors, herein we carried out linear and nonlinear infrared measurements on the C≡N stretching vibration of the sidechain of 5-cyanotryptophan, 3-methyl-5-cyanoindole, in a series of protic and aprotic solvents. We found that while the C≡N stretching frequencies obtained in these solvents do not correlate well with any individual Kamlet-Taft solvent parameter, i.e., π* (polarizability), β (hydrogen bond accepting ability), and α (hydrogen bond donating ability), they do however, collapse on a straight line when plotted against σ = π* + β - α. This linear relationship provides a firm indication that both specific interactions, i.e., hydrogen-bonding interactions with the C≡N (through α) and indole N-H (through β) groups, and non-specific interactions with the molecule (through π*) work together to determine the C≡N stretching frequency, thus laying a quantitative framework for applying 5-cyanotryptophan to investigate the microscopic environment of proteins in a site-specific manner. Furthermore, two-dimensional and pump-probe infrared measurements revealed that a significant portion (∼31%) of the ground state bleach signal has a decay time constant of ∼12.3 ps, due to an additional vibrational relaxation channel, making it possible to use 5-cyanotryptophan to probe dynamics occurring on a timescale on the order of tens of picoseconds.

Dynamic regulation of lipid-protein interactions

We review the importance of helix motions for the function of several important categories of membrane proteins and for the properties of several model molecular systems. For voltage-gated potassium or sodium channels, sliding, tilting and/or rotational movements of the S4 helix accompanied by a swapping of cognate side-chain ion-pair interactions regulate the channel gating. In the seven-helix G protein-coupled receptors, exemplified by the rhodopsins, collective helix motions serve to activate the functional signaling. Peptides which initially associate with lipid-bilayer membrane surfaces may undergo dynamic transitions from surface-bound to tilted-transmembrane orientations, sometimes accompanied by changes in the molecularity, formation of a pore or, more generally, the activation of biological function. For single-span membrane proteins, such as the tyrosine kinases, an interplay between juxtamembrane and transmembrane domains is likely to be crucial for the regulation of dimer assembly that in turn is associated with the functional responses to external signals. Additionally, we note that experiments with designed single-span transmembrane helices offer fundamental insights into the molecular features that govern protein-lipid interactions. This article is part of a Special Issue entitled: Lipid-protein interactions.

3D hydrophobic moment vectors as a tool to characterize the surface polarity of amphiphilic peptides

The interaction of membranes with peptides and proteins is largely determined by their amphiphilic character. Hydrophobic moments of helical segments are commonly derived from their two-dimensional helical wheel projections, and the same is true for β-sheets. However, to the best of our knowledge, there exists no method to describe structures in three dimensions or molecules with irregular shape. Here, we define the hydrophobic moment of a molecule as a vector in three dimensions by evaluating the surface distribution of all hydrophilic and lipophilic regions over any given shape. The electrostatic potential on the molecular surface is calculated based on the atomic point charges. The resulting hydrophobic moment vector is specific for the instantaneous conformation, and it takes into account all structural characteristics of the molecule, e.g., partial unfolding, bending, and side-chain torsion angles. Extended all-atom molecular dynamics simulations are then used to calculate the equilibrium hydrophobic moments for two antimicrobial peptides, gramicidin S and PGLa, under different conditions. We show that their effective hydrophobic moment vectors reflect the distribution of polar and nonpolar patches on the molecular surface and the calculated electrostatic surface potential. A comparison of simulations in solution and in lipid membranes shows how the peptides undergo internal conformational rearrangement upon binding to the bilayer surface. A good correlation with solid-state NMR data indicates that the hydrophobic moment vector can be used to predict the membrane binding geometry of peptides. This method is available as a web application on http://www.ibg.kit.edu/HM/.

Solid-state nuclear magnetic resonance spectroscopy for membrane protein structure determination

Solid-state NMR (ssNMR) is a versatile technique that can provide high-resolution (sub-angstrom) structural data for integral membrane proteins embedded in native and model membrane environments. The methodologies for a priori structure determination have for the most part been developed using samples with crystalline and fibrous morphologies. However, the techniques are now being applied to large, polytopic membrane proteins including receptors, ion channels, and porins. ssNMR data may be used to annotate and refine existing structures in regions of the protein not fully resolved by crystallography (including ligand-binding sites and mobile solvent accessible loop regions). This review describes the spectroscopic experiments and data analysis methods (including assignment) used to generate high-resolution structural data for membrane proteins. We also consider the range of sample morphologies that are appropriate for study by this method.

Transferring the PRIMO Coarse-Grained Force Field to the Membrane Environment: Simulations of Membrane Proteins and Helix-Helix Association

An extension of the recently developed PRIMO coarse-grained force field to membrane environments, PRIMO-M, is described. The membrane environment is modeled with the heterogeneous dielectric generalized Born (HDGB) methodology that simply replaces the standard generalized Born model in PRIMO without further parametrization. The resulting model was validated by comparing amino acid insertion free energy profiles and application in molecular dynamics simulations of membrane proteins and membrane-interacting peptides. Membrane proteins with 148–661 amino acids show stable root-mean-squared-deviations (RMSD) between 2 and 4 Å for most systems. Transmembrane helical peptides maintain helical shape and exhibit tilt angles in good agreement with experimental or other simulation data. The association of two glycophorin A (GpA) helices was simulated using replica exchange molecular dynamics simulations yielding the correct dimer structure with a crossing angle in agreement with previous studies. Finally, conformational sampling of the influenza fusion peptide also generates structures in agreement with previous studies. Overall, these findings suggest that PRIMO-M can be used to study membrane bound peptides and proteins and validates the transferable nature of the PRIMO coarse-grained force field.

Comparisons of interfacial Phe, Tyr, and Trp residues as determinants of orientation and dynamics for GWALP transmembrane peptides

Aromatic amino acids often flank the transmembrane alpha helices of integral membrane proteins. By favoring locations within the membrane–water interface of the lipid bilayer, aromatic residues Trp, Tyr, and sometimes Phe may serve as anchors to help stabilize a transmembrane orientation. In this work, we compare the influence of interfacial Trp, Tyr, or Phe residues upon the properties of tilted helical transmembrane peptides. For such comparisons, it has been critical to start with no more than one interfacial aromatic residue near each end of a transmembrane helix, for example, that of GWALP23 (acetyl-GGALW⁵(LA)₆LW¹⁹LAGA-[ethanol]amide). To this end, we have employed ²H-labeled alanines and solid-state NMR spectroscopy to investigate the consequences of moving or replacing W5 or W19 in GWALP23 with selected Tyr, Phe, or Trp residues at the same or proximate locations. We find that GWALP23 peptides having F5, Y5, or W5 exhibit essentially the same average tilt and similar dynamics in bilayer membranes of 1,2-dilauroylphosphatidylcholine (DLPC) or 1,2-dioleoylphosphatidylcholine (DOPC). When double Tyr anchors are present, in Y^4,5GWALP23 the NMR observables are markedly more subject to dynamic averaging and at the same time are less responsive to the bilayer thickness. Decreased dynamics are nevertheless observed when ring hydrogen bonding is removed, such that F^4,5GWALP23 exhibits a similar extent of low dynamic averaging as GWALP23 itself. When F5 is the sole aromatic group in the N-interfacial region, the dynamic averaging is (only) slightly more extensive than with W5, Y5, or Y4 alone or with F4,5, yet it is much less than that observed for Y^4,5GWALP23. Interestingly, moving Y5 to Y4 or W19 to W18, while retaining only one hydrogen-bond-capable aromatic ring at each interface, maintains the low level of dynamic averaging but alters the helix azimuthal rotation. The rotation change is about 40° for Y4 regardless of whether the host lipid bilayer is DLPC or DOPC. The rotational change (Δρ) is more dramatic and more complex when W19 is moved to W18, as Δρ is about +90° in DLPC but about −60° in DOPC. Possible reasons for this curious lipid-dependent helix rotation could include not only the separation distances between flanking aromatic or hydrophobic residues but also the absolute location of the W19 indole ring. For the more usual cases, when the helix azimuthal rotation shows little dependence on the host bilayer identity, excepting W¹⁸GWALP23, the transmembrane helices adapt to different lipids primarily by changing the magnitude of their tilt. We conclude that, in the absence of other functional groups, interfacial aromatic residues determine the preferred orientations and dynamics of membrane-spanning peptides. The results furthermore suggest possibilities for rotational and dynamic control of membrane protein function.

Sensitivity and resolution enhancement of oriented solid-state NMR: application to membrane proteins

Oriented solid-state NMR (O-ssNMR) spectroscopy is a major technique for the high-resolution analysis of the structure and topology of transmembrane proteins in native-like environments. Unlike magic angle spinning (MAS) techniques, O-ssNMR spectroscopy requires membrane protein preparations that are uniformly oriented (mechanically or magnetically) so that anisotropic NMR parameters, such as dipolar and chemical shift interactions, can be measured to determine structure and orientation of membrane proteins in lipid bilayers. Traditional sample preparations involving mechanically aligned lipids often result in short relaxation times which broaden the (15)N resonances and encumber the manipulation of nuclear spin coherences. The introduction of lipid bicelles as membrane mimicking systems has changed this scenario, and the more favorable relaxation properties of membrane protein (15)N and (13)C resonances make it possible to develop new, more elaborate pulse sequences for higher spectral resolution and sensitivity. Here, we describe our recent progress in the optimization of O-ssNMR pulse sequences. We explain the theory behind these experiments, demonstrate their application to small and medium size proteins, and describe the technical details for setting up these new experiments on the new generation of NMR spectrometers.

Single tryptophan and tyrosine comparisons in the N-terminal and C-terminal interface regions of transmembrane GWALP peptides

Hydrophobic membrane-spanning helices often are flanked by interfacial aromatic or charged residues. In this paper, we compare the consequences of single Trp → Tyr substitutions at each interface for the properties of a defined transmembrane helix in the absence of charged residues. The choice of molecular framework is critical for these single-residue experiments because the presence of "too many" aromatic residues (more than one at either membrane-water interface) introduces excess dynamic averaging of solid state NMR observables. To this end, we compare the outcomes when changing W(5) or W(19), or both of them, to tyrosine in the well-characterized transmembrane peptide acetyl-GGALW(5)(LA)6LW(19)LAGA-amide ("GWALP23"). By means of solid-state (2)H and (15)N NMR experiments, we find that Y(19)GW(5)ALP23 displays similar magnitudes of peptide helix tilt as Y(5)GW(19)ALP23 and responds similarly to changes in bilayer thickness, from DLPC to DMPC to DOPC. The presence of Y(19) changes the azimuthal rotation angle ρ (about the helix axis) to a similar extent as Y(5), but in the opposite direction. When tyrosines are substituted for both tryptophans to yield GY(5,19)ALP23, the helix tilt angle is again of comparable magnitude, and furthermore, the preferred azimuthal rotation angle ρ is relatively unchanged from that of GW(5,19)ALP23. The extent of dynamic averaging increases marginally when Tyr replaces Trp. Yet, importantly, all members of the peptide family having single Tyr or Trp residues near each interface exhibit only moderate and not highly extensive dynamic averaging. The results provide important benchmarks for evaluating conformational and dynamic control of membrane protein function.

Importance of indole N-H hydrogen bonding in the organization and dynamics of gramicidin channels

The linear ion channel peptide gramicidin represents an excellent model for exploring the principles underlying membrane protein structure and function, especially with respect to tryptophan residues. The tryptophan residues in gramicidin channels are crucial for the structure and function of the channel. In order to test the importance of indole hydrogen bonding for the biophysical properties of gramicidin channels, we monitored the effect of N-methylation of gramicidin tryptophans, using a combination of steady state and time-resolved fluorescence approaches along with circular dichroism spectroscopy. We show here that in the absence of the hydrogen bonding ability of tryptophans, tetramethyltryptophan gramicidin (TM-gramicidin) is unable to maintain the single stranded, head-to-head dimeric channel conformation in membranes. Our results show that TM-gramicidin displays a red-shifted fluorescence emission maximum, lower red edge excitation shift (REES), and higher fluorescence intensity and lifetime, consistent with its nonchannel conformation. This is in agreement with the measured location (average depth) of the 1-methyltryptophans in TM-gramicidin using the parallax method. These results bring out the usefulness of 1-methyltryptophan as a fluorescent tool to examine the hydrogen bonding ability of tryptophans in proteins and peptides. We conclude that changes in the hydrogen bonding ability of tryptophans, along with coupled changes in peptide backbone structure induce the loss of single stranded β(6.3) helical dimer conformation. These results agree with earlier results from size-exclusion chromatography and single-channel measurements for TM-gramicidin, and confirm the importance of indole hydrogen bonding for the conformation and function of ion channels and membrane proteins.

Canonical azimuthal rotations and flanking residues constrain the orientation of transmembrane helices

In biological membranes the alignment of embedded proteins provides crucial structural information. The transmembrane (TM) parts have well-defined secondary structures, in most cases α-helices and their orientation is given by a tilt angle and an azimuthal rotation angle around the main axis. The tilt angle is readily visualized and has been found to be functionally relevant. However, there exist no general concepts on the corresponding azimuthal rotation. Here, we show that TM helices prefer discrete rotation angles. They arise from a combination of intrinsic properties of the helix geometry plus the influence of the position and type of flanking residues at both ends of the hydrophobic core. The helical geometry gives rise to canonical azimuthal angles for which the side chains of residues from the two ends of the TM helix tend to have maximum or minimum immersion within the membrane. This affects the preferential position of residues that fall near hydrophobic/polar interfaces of the membrane, depending on their hydrophobicity and capacity to form specific anchoring interactions. On this basis, we can explain the orientation and dynamics of TM helices and make accurate predictions, which correspond well to the experimental values of several model peptides (including dimers), and TM segments of polytopic membrane proteins.

Combination of ¹⁵N reverse labeling and afterglow spectroscopy for assigning membrane protein spectra by magic-angle-spinning solid-state NMR: application to the multidrug resistance protein EmrE

Magic-angle-spinning (MAS) solid-state NMR spectroscopy has emerged as a viable method to characterize membrane protein structure and dynamics. Nevertheless, the spectral resolution for uniformly labeled samples is often compromised by redundancy of the primary sequence and the presence of helical secondary structure that results in substantial resonance overlap. The ability to simplify the spectrum in order to obtain unambiguous site-specific assignments is a major bottleneck for structure determination. To address this problem, we used a combination of (15)N reverse labeling, afterglow spectroscopic techniques, and frequency-selective dephasing experiments that dramatically improved the ability to resolve peaks in crowded spectra. This was demonstrated using the polytopic membrane protein EmrE, an efflux pump involved in multidrug resistance. Residues preceding the (15)N reverse labeled amino acid were imaged using a 3D NCOCX afterglow experiment and those following were recorded using a frequency-selective dephasing experiment. Our approach reduced the spectral congestion and provided a sensitive way to obtain chemical shift assignments for a membrane protein where no high-resolution structure is available. This MAS methodology is widely applicable to the study of other polytopic membrane proteins in functional lipid bilayer environments.

Buried lysine, but not arginine, titrates and alters transmembrane helix tilt

The ionization states of individual amino acid residues of membrane proteins are difficult to decipher or assign directly in the lipid-bilayer membrane environment. We address this issue for lysines and arginines in designed transmembrane helices. For lysines (but not arginines) at two locations within dioleoyl-phosphatidylcholine bilayer membranes, we measure pK(a) values below 7.0. We find that buried charged lysine, in fashion similar to arginine, will modulate helix orientation to maximize its own access to the aqueous interface or, if occluded by aromatic rings, may cause a transmembrane helix to exit the lipid bilayer. Interestingly, the influence of neutral lysine (vis-à-vis leucine) upon helix orientation also depends upon its aqueous access. Our results suggest that changes in the ionization states of particular residues will regulate membrane protein function and furthermore illustrate the subtle complexity of ionization behavior with respect to the detailed lipid and protein environment.

Three-dimensional stress field around a membrane protein: atomistic and coarse-grained simulation analysis of gramicidin A

Using both atomistic and coarse-grained (CG) models, we compute the three-dimensional stress field around a gramicidin A (gA) dimer in lipid bilayers that feature different degrees of negative hydrophobic mismatch. The general trends in the computed stress field are similar at the atomistic and CG levels, supporting the use of the CG model for analyzing the mechanical features of protein/lipid/water interfaces. The calculations reveal that the stress field near the protein-lipid interface exhibits a layered structure with both significant repulsive and attractive regions, with the magnitude of the stress reaching 1000 bar in certain regions. Analysis of density profiles and stress field distributions helps highlight the Trp residues at the protein/membrane/water interface as mechanical anchors, suggesting that similar analysis is useful for identifying tension sensors in other membrane proteins, especially membrane proteins involved in mechanosensation. This work fosters a connection between microscopic and continuum mechanics models for proteins in complex environments and makes it possible to test the validity of assumptions commonly made in continuum mechanics models for membrane mediated processes. For example, using the calculated stress field, we estimate the free energy of membrane deformation induced by the hydrophobic mismatch, and the results for regions beyond the annular lipids are in general consistent with relevant experimental data and previous theoretical estimates using elasticity theory. On the other hand, the assumptions of homogeneous material properties for the membrane and a bilayer thickness at the protein/lipid interface being independent of lipid type (e.g., tail length) appear to be oversimplified, highlighting the importance of annular lipids of membrane proteins. Finally, the stress field analysis makes it clear that the effect of even rather severe hydrophobic mismatch propagates to only about two to three lipid layers, thus putting a limit on the range of cooperativity between membrane proteins in crowded cellular membranes.

Properties of membrane-incorporated WALP peptides that are anchored on only one end

Peptides of the "WALP" family, acetyl-GWW(LA)(n)LWWA-[ethanol]amide, have proven to be opportune models for investigating lipid-peptide interactions. Because the average orientations and motional behavior of the N- and C-terminal Trp (W) residues differ, it is of interest to investigate how the positions of the tryptophans influence the properties of the membrane-incorporated peptides. To address this question, we synthesized acetyl-GGWW(LA)(n)-ethanolamide and acetyl-(AL)(n)WWG-ethanolamide, in which n = 4 or 8, which we designate as "N-anchored" and "C-anchored" peptides, respectively. Selected (2)H or (15)N labels were incorporated for solid-state nuclear magnetic resonance (NMR) spectroscopy. These peptides can be considered "half"-anchored WALP peptides, having only one pair of interfacial Trp residues near either the amino or the carboxyl terminus. The hydrophobic lengths of the (n = 8) peptides are similar to that of WALP23. These longer half-anchored WALP peptides incorporate into lipid bilayers as α-helices, as reflected in their circular dichroism spectra. Solid-state NMR experiments indicate that the longer peptide helices assume defined transmembrane orientations with small non-zero average tilt angles and moderate to high dynamic averaging in bilayer membranes of 1,2-dioleoylphosphatidylcholine, 1,2-dimyristoylphosphatidylcholine, and 1,2-dilauroylphosphatidylcholine. The intrinsically small apparent tilt angles suggest that interactions of aromatic residues with lipid headgroups may play an important role in determining the magnitude of the peptide tilt in the bilayer membrane. The shorter (n = 4) peptides, in stark contrast to the longer peptides, display NMR spectra that are characteristic of greatly reduced motional averaging, probably because of peptide aggregation in the bilayer environment, and CD spectra that are characteristic of β-structure.

Accommodation of a central arginine in a transmembrane peptide by changing the placement of anchor residues

Both Trp and Arg in transmembrane protein domains make important interactions with lipids at the membrane/water interface, but at different depths. Derivatives of the designed peptide GWALP23, acetyl-GGALW(5)LALALALALALALW(19)LAGA-amide, with single Trp anchors, have proven useful for characterizing such interactions. Indeed, previous work revealed quite different effects emanating from Arg substitutions at positions 12 and 14 within GWALP23, with the R12 peptide exhibiting multiple positions and orientations with respect to DOPC bilayer membranes (Vostrikov et al. J. Am. Chem. Soc. 2010, 132, 5803-5811). To gain further understanding of the multistate behavior, we moved the Trp "anchor" residues to more outer positions 3 and 21 in GWALP23 itself, and in the R12 and R14 derivatives. The locations and orientations of the peptides with respect to lipid bilayer membranes of differing thickness were investigated by means of solid-state (2)H NMR spectroscopy, using labeled alanines, and coarse-grained molecular dynamics simulations. Interestingly, relatively intense and narrow (2)H resonances from selected backbone C(α) deuterons were observed over quite narrow ranges of frequency and sample orientation. The backbone resonances reflect dynamic complexities and at the same time provide important contributions for the analysis of peptide transmembrane orientation. With the Trp(3,21) anchors relatively far from the peptide and bilayer center, the results indicate significantly large apparent tilt angles, for example, close to 30° for the new R12 and R14 peptides with respect to the bilayer normal of DLPC membranes. The R12 side chain indeed is "rescued" to a stable position, where it is accommodated within the transmembrane helix, when the Trp anchors are moved outward and to another face of the helix. At the same time, the R14 side chain of transmembrane GW(3,21)ALP23 also retains a stable favored position.

Mechanism of bacterial signal transduction revealed by molecular dynamics of Tsr dimers and trimers of dimers in lipid vesicles

Bacterial chemoreceptors provide an important model for understanding signalling processes. In the serine receptor Tsr from E. coli, a binding event in the periplasmic domain of the receptor dimer causes a shift in a single transmembrane helix of roughly 0.15 nm towards the cytoplasm. This small change is propagated through the ∼22 nm length of the receptor, causing downstream inhibition of the kinase CheA. This requires interactions within a trimer of receptor dimers. Additionally, the signal is amplified across a 53,000 nm² array of chemoreceptor proteins, including ∼5,200 receptor trimers-of-dimers, at the cell pole. Despite a wealth of experimental data on the system, including high resolution structures of individual domains and extensive mutagenesis data, it remains uncertain how information is communicated across the receptor from the binding event to the downstream effectors. We present a molecular model of the entire Tsr dimer, and examine its behaviour using coarse-grained molecular dynamics and elastic network modelling. We observe a large bending in dimer models between the linker domain HAMP and coiled-coil domains, which is supported by experimental data. Models of the trimer of dimers, built from the dimer models, are more constrained and likely represent the signalling state. Simulations of the models in a 70 nm diameter vesicle with a biologically realistic lipid mixture reveal specific lipid interactions and oligomerisation of the trimer of dimers. The results indicate a mechanism whereby small motions of a single helix can be amplified through HAMP domain packing, to initiate large changes in the whole receptor structure.

To understand cell signalling events requires a physical model of the structure and behaviour of the signalling proteins involved. The methyl-accepting chemoreceptor proteins direct bacterial movement towards food sources and away from toxins. Based on experimental data we have built structural models of the serine chemoreceptor (Tsr) as a dimer, which is incapable of activating the downstream kinase CheA, and as a trimer of dimers, which can activate CheA. We have performed molecular dynamics simulation to reveal the behaviour of these two forms in a planar lipid bilayer and in a 70 nm diameter lipid vesicle with a mixture of lipids mimicking the E. coli inner membrane. We show that in isolation the dimers undergo a bending movement around the central HAMP domain, whereas the trimer-of-dimers model does not. Comparison with published experimental data suggests that these bending motions are real, and that they occur in the trimer of dimers only in response to ligand binding. Drawing together these observations with studies showing that the signalling event involves small piston motions in the transmembrane helices suggests that the bending motion is frustrated in the unliganded trimer of dimers, and that ligand binding induces bending by repacking the HAMP interface.

Osteopetrosis mutation R444L causes endoplasmic reticulum retention and misprocessing of vacuolar H+-ATPase a3 subunit

Osteopetrosis is a genetic bone disease characterized by increased bone density and fragility. The R444L missense mutation in the human V-ATPase a3 subunit (TCIRG1) is one of several known mutations in a3 and other proteins that can cause this disease. The autosomal recessive R444L mutation results in a particularly malignant form of infantile osteopetrosis that is lethal in infancy, or early childhood. We have studied this mutation using the pMSCV retroviral vector system to integrate the cDNA construct for green fluorescent protein (GFP)-fused a3(R445L) mutant protein into the RAW 264.7 mouse osteoclast differentiation model. In comparison with wild-type a3, the mutant glycoprotein localized to the ER instead of lysosomes and its oligosaccharide moiety was misprocessed, suggesting inability of the core-glycosylated glycoprotein to traffic to the Golgi. Reduced steady-state expression of the mutant protein, in comparison with wild type, suggested that the former was being degraded, likely through the endoplasmic reticulum-associated degradation pathway. In differentiated osteoclasts, a3(R445L) was found to degrade at an increased rate over the course of osteoclastogenesis. Limited proteolysis studies suggested that the R445L mutation alters mouse a3 protein conformation. Together, these data suggest that Arg-445 plays a role in protein folding, or stability, and that infantile malignant osteopetrosis caused by the R444L mutation in the human V-ATPase a3 subunit is another member of the growing class of protein folding diseases. This may have implications for early-intervention treatment, using protein rescue strategies.

The Transmembrane Helix Tilt May Be Determined by the Balance between Precession Entropy and Lipid Perturbation

Hydrophobic helical peptides interact with lipid bilayers in various modes, determined by the match between the length of the helix’s hydrophobic core and the thickness of the hydrocarbon region of the bilayer. For example, long helices may tilt with respect to the membrane normal to bury their hydrophobic cores in the membrane, and the lipid bilayer may stretch to match the helix length. Recent molecular dynamics simulations and potential of mean force calculations have shown that some TM helices whose lengths are equal to, or even shorter than, the bilayer thickness may also tilt. The tilt is driven by a gain in the helix precession entropy, which compensates for the free energy penalty resulting from membrane deformation. Using this free energy balance, we derived theoretically an equation of state, describing the dependence of the tilt on the helix length and membrane thickness. To this end, we conducted coarse-grained Monte Carlo simulations of the interaction of helices of various lengths with lipid bilayers of various thicknesses, reproducing and expanding the previous molecular dynamics simulations. Insight from the simulations facilitated the derivation of the theoretical model. The tilt angles calculated using the theoretical model agree well with our simulations and with previous calculations and measurements.

Proline kink angle distributions for GWALP23 in lipid bilayers of different thicknesses

By using selected (2)H and (15)N labels, we have examined the influence of a central proline residue on the properties of a defined peptide that spans lipid bilayer membranes by solid-state nuclear magnetic resonance (NMR) spectroscopy. For this purpose, GWALP23 (acetyl-GGALW(5)LALALALALALALW(19)LAGA-ethanolamide) is a suitable model peptide that employs, for the purpose of interfacial anchoring, only one tryptophan residue on either end of a central α-helical core sequence. Because of its systematic behavior in lipid bilayer membranes of differing thicknesses [Vostrikov, V. V., et al. (2010) J. Biol. Chem. 285, 31723-31730], we utilize GWALP23 as a well-characterized framework for introducing guest residues within a transmembrane sequence; for example, a central proline yields acetyl-GGALW(5)LALALAP(12)ALALALW(19)LAGA-ethanolamide. We synthesized GWALP23-P12 with specifically placed (2)H and (15)N labels for solid-state NMR spectroscopy and examined the peptide orientation and segmental tilt in oriented DMPC lipid bilayer membranes using combined (2)H GALA and (15)N-(1)H high-resolution separated local field methods. In DMPC bilayer membranes, the peptide segments N-terminal and C-terminal to the proline are both tilted substantially with respect to the bilayer normal, by ~34 ± 5° and 29 ± 5°, respectively. While the tilt increases for both segments when proline is present, the range and extent of the individual segment motions are comparable to or smaller than those of the entire GWALP23 peptide in bilayer membranes. In DMPC, the proline induces a kink of ~30 ± 5°, with an apparent helix unwinding or "swivel" angle of ~70°. In DLPC and DOPC, on the basis of (2)H NMR data only, the kink angle and swivel angle probability distributions overlap those of DMPC, yet the most probable kink angle appears to be somewhat smaller than in DMPC. As has been described for GWALP23 itself, the C-terminal helix ends before Ala(21) in the phospholipids DMPC and DLPC yet remains intact through Ala(21) in DOPC. The dynamics of bilayer-incorporated, membrane-spanning GWALP23 and GWALP23-P12 are less extensive than those observed for WALP family peptides that have more than two interfacial Trp residues.

Tyrosine replacing tryptophan as an anchor in GWALP peptides

Synthetic model peptides have proven useful for examining fundamental peptide-lipid interactions. A frequently employed peptide design consists of a hydrophobic core of Leu-Ala residues with polar or aromatic amino acids flanking each side at the interfacial positions, which serve to "anchor" a specific transmembrane orientation. For example, WALP family peptides (acetyl-GWW(LA)(n)LWWA-[ethanol]amide), anchored by four Trp residues, have received particular attention in both experimental and theoretical studies. A recent modification proved successful in reducing the number of Trp anchors to only one near each end of the peptide. The resulting GWALP23 (acetyl-GGALW(5)(LA)(6)LW(19)LAGA-[ethanol]amide) displays reduced dynamics and greater sensitivity to lipid-peptide hydrophobic mismatch than traditional WALP peptides. We have further modified GWALP23 to incorporate a single tyrosine, replacing W(5) with Y(5). The resulting peptide, Y(5)GWALP23 (acetyl-GGALY(5)(LA)(6)LW(19)LAGA-amide), has a single Trp residue that is sensitive to fluorescence experiments. By incorporating specific (2)H and (15)N labels in the core sequence of Y(5)GWALP23, we were able to use solid-state NMR spectroscopy to examine the peptide orientation in hydrated lipid bilayer membranes. The peptide orients well in membranes and gives well-defined (2)H quadrupolar splittings and (15)N/(1)H dipolar couplings throughout the core helical sequence between the aromatic residues. The substitution of Y(5) for W(5) has remarkably little influence on the tilt or dynamics of GWALP23 in bilayer membranes of the phospholipids DOPC, DMPC, or DLPC. A second analogue of the peptide with one Trp and two Tyr anchors, Y(4,5)GWALP23, is generally less responsive to the bilayer thickness and exhibits lower apparent tilt angles with evidence of more extensive dynamics. In general, the peptide behavior with multiple Tyr anchors appears to be quite similar to the situation when multiple Trp anchors are present, as in the original WALP series of model peptides.

Hydrophobic mismatch of mobile transmembrane helices: Merging theory and experiments

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

On the combined analysis of ²H and ¹⁵N/¹H solid-state NMR data for determination of transmembrane peptide orientation and dynamics

The dynamics of membrane-spanning peptides have a strong affect on the solid-state NMR observables. We present a combined analysis of ²H-alanine quadrupolar splittings together with ¹⁵N/(1)H dipolar couplings and ¹⁵N chemical shifts, using two models to treat the dynamics, for the systematic evaluation of transmembrane peptides based on the GWALP23 sequence (acetyl-GGALW(LA)₆LWLAGA-amide). The results indicate that derivatives of GWALP23 incorporating diverse guest residues adopt a range of apparent tilt angles that span 5°-35° in lipid bilayer membranes. By comparing individual and combined analyses of specifically ²H- or ¹⁵N-labeled peptides incorporated in magnetically or mechanically aligned lipid bilayers, we examine the influence of data-set size/identity, and of explicitly modeled dynamics, on the deduced average orientations of the peptides. We conclude that peptides with small apparent tilt values (<∼10°) can be fitted by extensive families of solutions, which can be narrowed by incorporating additional ¹⁵N as well as ²H restraints. Conversely, peptides exhibiting larger tilt angles display more narrow distributions of tilt and rotation that can be fitted using smaller sets of experimental constraints or even with ²H or ¹⁵N data alone. Importantly, for peptides that tilt significantly more than 10° from the bilayer-normal, the contribution from rigid body dynamics can be approximated by a principal order parameter.