Comparison of "Polarization inversion with spin exchange at magic angle" and "geometric analysis of labeled alanines" methods for transmembrane helix alignment

Exploring the trimerization process of a transmembrane helix with an ionizable residue by molecular dynamics simulations: a case study of transmembrane domain 5 of LMP-1

The oligomerization of membrane proteins is an important biological process that plays a critical role in the initialization of membrane protein receptor signaling. Unveiling how transmembrane segments oligomerize is critical for understanding the mechanism of membrane receptor signaling activation. Owing to the complicated membrane environment and the extraordinary dynamic properties of the ionizable residues in the transmembrane segment, it is extremely challenging to thoroughly understand the oligomerization process of the transmembrane domain. In this study, transmembrane domain 5 (TMD5) of latent membrane protein-1 from Epstein-Barr virus was used as a prototype model to investigate the trimerization process of the transmembrane segment with ionizable residues. The trimerization process of TMD5 was rebuilt and investigated via conventional molecular dynamics simulations and constant-pH molecular dynamics simulations. When TMD5s approached each other, the tilting angles of the TMD5 monomer decreased. TMD5s formed stable trimers until two interacting sites (D150s and Q139s) along each transmembrane helix were created to lock the TMD5s. The pK_a values of D150 shifted toward neutral states in the membrane environment. When TMD5s were monomers, the pK_a shift of D150 was mainly influenced by its microenvironment in the lipid bilayer. When TMD5s were moving close to each other, protein-protein interactions became the main contributing factor for the pK_a shift of D150s. Overall, this work elucidates the behavior of the TMD5 helix and the pK_a shift of ionizable residue D150 in the process of TMD5 oligomerization. This study may provide insight into the development of agents for targeting the oligomerization of membrane proteins.

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.

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.

Peptide-bicelle interaction: Following variations in size and morphology by a combined NMR-SAXS approach

Changes in membrane properties occurring upon protein interaction are key questions in understanding membrane protein function. To report on the occurring size and shape variation we present here a combined NMR-SAXS method performed under physiological conditions using the same samples, enabling determination of a global parameter, the hydration radius (r_H) and estimating the bicelle shape. We use zwitterionic (DMPC/DHPC) and negatively charged (DMPC/DHPC/DMPG) bicelles and investigate the interaction with model transmembrane and surface active peptides (KALP23 and melittin). ¹H NMR measurements based mostly on the translational diffusion coefficient D determination are used to characterize cmc values of DHPC micelles under the investigated conditions, to describe DHPC distribution with exact determination of the q (long chain/short chain) lipid ratio, to estimate aggregation numbers and effective r_H values. The scattering curve is used to fit a lenticular core-shell model enabling us to describe the bicelle shape in terms of ellipsoidal axis length parameters. For all studied systems formation of oblate ellipsoids is found. Even though the r_G/r_H ratio would be an elegant way to characterize shape variations, we show that changes occurring upon peptide-bicelle interaction in the "effective" size and in the measure on the anisometry - morphology - of the objects can be described by using r_H and the simplistic ellipsoidal core-shell model. While the influence of the transmembrane KALP peptide is significant, effects upon addition of surface active melittin peptide seem negligible. This synergy of techniques under controlled conditions can provide information about bicellar shape modulation occurring during peptide-bicelle interactions.

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.

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.

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.

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.

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.

Alanine scan and (2)H NMR analysis of the membrane-active peptide BP100 point to a distinct carpet mechanism of action

The short membrane-active peptide BP100 [KKLFKKILKYL-NH2] is known as an effective antimicrobial and cell penetrating agent. For a functional alanine scan each of the 11 amino acids was replaced with deuterated Ala-d3, one at a time. MIC assays showed that a substitution of Lys did not affect the antimicrobial activity, but it decreased when a hydrophobic residue was replaced. In most cases, a reduction in hydrophobicity led to a decrease in hemolysis, and some peptide analogues had an improved therapeutic index. Circular dichroism showed that BP100 folds as an amphiphilic α-helix in a bilayer. Its alignment was determined from (2)H NMR in oriented membranes of different composition. The azimuthal rotation angle was the same under all conditions, but the average helix tilt angle and the dynamical behavior of the peptide varied in a systematic manner. In POPC/POPG bilayers, with a negative spontaneous curvature, the peptide was found to lie flat on the bilayer surface, and with little wobble. In DMPC/DMPG, with a positive spontaneous curvature, BP100 at higher concentrations became tilted obliquely into the membrane, with the uncharged C-terminus inserted more deeply into the lipid bilayer, experiencing significant fluctuations in tilt angle. In DMPC/DMPG/lyso-MPC, with a pronounced positive spontaneous curvature, the helix tilted even further and became even more mobile. The 11-mer BP100 is obviously too short to form transmembrane pores. We conclude that BP100 operates via a carpet mechanism, whereby the C-terminus gets inserted into the hydrophobic core of the bilayer, which leads to membrane perturbation and induces transient permeability.

Membrane environment modulates the pKa values of transmembrane helices

In this work, we apply the recently developed constant pH molecular dynamics technique to study protonation equilibria of titratable side chains in the context of simple transmembrane (TM) helices and explore the effect of pH on their configurations in membrane bilayers. We observe that, despite a significant shift toward neutral states, considerable population of different side chains stay in the charged state that give rise to pK_a values around 9.6 for Asp and Glu and 4.5 to 6 for His and Lys side chains, respectively. These charged states are highly stabilized by favorable interactions between head groups, water molecules, and the charged side chains that are facilitated by substantial changes in the configuration of the peptides. The pH dependent configurations and the measured pK_a values are in good agreement with relatively recent solid state NMR measurements. Our results presented here demonstrate that all-atom constant pH molecular dynamics can be applied to membrane proteins and peptides to obtain reliable pK_a values and pH dependent behavior for these systems.

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.

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.

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.

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.

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.

Determining the orientation and localization of membrane-bound peptides

Many naturally occurring bioactive peptides bind to biological membranes. Studying and elucidating the mode of interaction is often an essential step to understand their molecular and biological functions. To obtain the complete orientation and immersion depth of such compounds in the membrane or a membrane-mimetic system, a number of methods are available, which are separated in this review into four main classes: solution NMR, solid-state NMR, EPR and other methods. Solution NMR methods include the Nuclear Overhauser Effect (NOE) between peptide and membrane signals, residual dipolar couplings and the use of paramagnetic probes, either within the membrane-mimetic or in the solvent. The vast array of solid state NMR methods to study membrane-bound peptide orientation and localization includes the anisotropic chemical shift, PISA wheels, dipolar waves, the GALA, MAOS and REDOR methods and again the use of paramagnetic additives on relaxation rates. Paramagnetic additives, with their effect on spectral linewidths, have also been used in EPR spectroscopy. Additionally, the orientation of a peptide within a membrane can be obtained by the anisotropic hyperfine tensor of a rigidly attached nitroxide label. Besides these magnetic resonance techniques a series of other methods to probe the orientation of peptides in membranes has been developed, consisting of fluorescence-, infrared- and oriented circular dichroism spectroscopy, colorimetry, interface-sensitive X-ray and neutron scattering and Quartz crystal microbalance.

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.

Solid-state NMR ensemble dynamics as a mediator between experiment and simulation

Solid-state NMR (SSNMR) is a powerful technique to describe the orientations of membrane proteins and peptides in their native membrane bilayer environments. The deuterium ((2)H) quadrupolar splitting (DQS), one of the SSNMR observables, has been used to characterize the orientations of various single-pass transmembrane (TM) helices using a semistatic rigid-body model such as the geometric analysis of labeled alanine (GALA) method. However, dynamic information of these TM helices, which could be related to important biological function, can be missing or misinterpreted with the semistatic model. We have investigated the orientation of WALP23 in an implicit membrane of dimyristoylglycerophosphocholine by determining an ensemble of structures using multiple conformer models with a DQS restraint potential. When a single conformer is used, the resulting helix orientation (tilt angle (τ) of 5.6 ± 3.2° and rotation angle (ρ) of 141.8 ± 40.6°) is similar to that determined by the GALA method. However, as the number of conformers is increased, the tilt angles of WALP23 ensemble structures become larger (26.9 ± 6.7°), which agrees well with previous molecular dynamics simulation results. In addition, the ensemble structure distribution shows excellent agreement with the two-dimensional free energy surface as a function of WALP23's τ and ρ. These results demonstrate that SSNMR ensemble dynamics provides a means to extract orientational and dynamic information of TM helices from their SSNMR observables and to explain the discrepancy between molecular dynamics simulation and GALA-based interpretation of DQS data.

Response of GWALP transmembrane peptides to changes in the tryptophan anchor positions

While the interfacial partitioning of charged or aromatic anchor residues may determine the preferred orientations of transmembrane peptide helices, the dependence of helix orientation on anchor residue position is not well understood. When anchor residue locations are changed systematically, some adaptations of the peptide-lipid interactions may be required to compensate for the altered interfacial interactions. Recently, we have developed a novel transmembrane peptide, termed GW(5,19)ALP23 (acetyl-GGALW(5)LALALALALALALW(19)LAGA-ethanolamide), which proves to be a well-behaved sequence for an orderly investigation of protein-lipid interactions. Its roughly symmetric nature allows for shifting the anchoring Trp residues by one Leu-Ala pair inward (GW(7,17)ALP23) or outward (GW(3,21)ALP23), thus providing fine adjustments of the formal distance between the tryptophan residues. With no other obvious anchoring features present, we postulate that the inter-Trp distance may be crucial for aspects of the peptide-lipid interaction. Importantly, the amino acid composition is identical for each of the resulting related GWALP23 sequences, and the radial separation between the pairs of Trp residues on each side of the transmembrane α-helix remains similar. Here we address the adaptation of the aforementioned peptides to the varying Trp locations by means of solid-state (2)H nuclear magnetic resonance experiments in varying lipid bilayer membrane environments. All of the GW(x,y)ALP23 sequence isomers adopt transmembrane orientations in DOPC, DMPC, and DLPC environments, even when the Trp residues are quite closely spaced, in GW(7,17)ALP23. Furthermore, the dynamics for each peptide isomer are less extensive than for peptides possessing additional interfacial Trp residues. The helical secondary structure is maintained more strongly within the Trp-flanked core region than outside of the Trp boundaries. Deuterium-labeled tryptophan indole rings in the GW(x,y)ALP23 peptides provide additional insights into the behavior of the Trp side chains. A Trp side chain near the C-terminus adopts a different orientation and undergoes somewhat faster dynamics than a corresponding Trp side chain located an equivalent distance from the N-terminus. In contrast, as the inter-Trp distance changes, the variations among the average orientations of the Trp indole rings at either terminus are systematic yet fairly small. We conclude that subtle adjustments to the peptide tilt, and to the N- and C-terminal Trp side chain torsion angles, permit the GW(x,y)ALP23 peptides to maintain preferred transmembrane orientations while adapting to lipid bilayers with differing hydrophobic thicknesses.

An ensemble dynamics approach to decipher solid-state NMR observables of membrane proteins

Solid-state NMR (SSNMR) is an invaluable tool for determining orientations of membrane proteins and peptides in lipid bilayers. Such orientational descriptions provide essential information about membrane protein functions. However, when a semi-static single conformer model is used to interpret various SSNMR observables, important dynamics information can be missing, and, sometimes, even orientational information can be misinterpreted. In addition, over the last decade, molecular dynamics (MD) simulation and semi-static SSNMR interpretation have shown certain levels of discrepancies in terms of transmembrane helix orientation and dynamics. Dynamic fitting models have recently been proposed to resolve these discrepancies by taking into account transmembrane helix whole body motions using additional parameters. As an alternative approach, we have developed SSNMR ensemble dynamics (SSNMR-ED) using multiple conformer models, which generates an ensemble of structures that satisfies the experimental observables without any fitting parameters. In this review, various computational methods for determining transmembrane helix orientations are discussed, and the distributions of VpuTM (from HIV-1) and WALP23 (a synthetic peptide) orientations from SSNMR-ED simulations are compared with those from MD simulations and semi-static/dynamic fitting models. Such comparisons illustrate that SSNMR-ED can be used as a general means to extract both membrane protein structure and dynamics from the SSNMR measurements. This article is part of a Special Issue entitled: Membrane protein structure and function.

The structural and topological analysis of membrane-associated polypeptides by oriented solid-state NMR spectroscopy: established concepts and novel developments

Solid-state NMR spectroscopy is a powerful technique for the investigation of membrane-associated peptides and proteins as well as their interactions with lipids, and a variety of conceptually different approaches have been developed for their study. The technique is unique in allowing for the high-resolution investigation of liquid disordered lipid bilayers representing well the characteristics of natural membranes. Whereas magic angle solid-state NMR spectroscopy follows approaches that are related to those developed for solution NMR spectroscopy the use of static uniaxially oriented samples results in angular constraints which also provide information for the detailed analysis of polypeptide structures. This review introduces this latter concept theoretically and provides a number of examples. Furthermore, ongoing developments combining solid-state NMR spectroscopy with information from solution NMR spectroscopy and molecular modelling as well as exploratory studies using dynamic nuclear polarization solid-state NMR will be presented.

Molecular dynamics of membrane peptides and proteins: principles and comparison to experimental data

Molecular dynamics (MD) simulation is a standard tool used to assess the motion of biomolecules at atomic resolution. It requires a so-called "force field" that allows the evaluation of an empirical energy from the 3D coordinates of the atoms in the system. In this chapter, the application of MD simulations to membrane proteins and peptides is described with a particular emphasis on the comparison of MD results to experimental data. Such a comparison can be used either for (1) validating the results of a simulation, (2) interpreting an experiment at the atomic level, or (3) calibrating the force field. This last step is particularly important for the use of MD as a predictive tool. As an illustration, a comparison of (2)H NMR experiments to MD simulations of a transmembrane peptide is presented and discussed.

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)(n)LWWA-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(2)ALW(LA)(6)LWLAX(22)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 (2)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(5) and Trp(19)) 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.

Interpretation of 2H-NMR experiments on the orientation of the transmembrane helix WALP23 by computer simulations

Orientation, dynamics, and packing of transmembrane helical peptides are important determinants of membrane protein structure, dynamics, and function. Because it is difficult to investigate these aspects by studying real membrane proteins, model transmembrane helical peptides are widely used. NMR experiments provide information on both orientation and dynamics of peptides, but they require that motional models be interpreted. Different motional models yield different interpretations of quadrupolar splittings (QS) in terms of helix orientation and dynamics. Here, we use coarse-grained (CG) molecular dynamics (MD) simulations to investigate the behavior of a well-known model transmembrane peptide, WALP23, under different hydrophobic matching/mismatching conditions. We compare experimental (2)H-NMR QS (directly measured in experiments), as well as helix tilt angle and azimuthal rotation (not directly measured), with CG MD simulation results. For QS, the agreement is significantly better than previously obtained with atomistic simulations, indicating that equilibrium sampling is more important than atomistic details for reproducing experimental QS. Calculations of helix orientation confirm that the interpretation of QS depends on the motional model used. Our simulations suggest that WALP23 can form dimers, which are more stable in an antiparallel arrangement. The origin of the preference for the antiparallel orientation lies not only in electrostatic interactions but also in better surface complementarity. In most cases, a mixture of monomers and antiparallel dimers provides better agreement with NMR data compared to the monomer and the parallel dimer. CG MD simulations allow predictions of helix orientation and dynamics and interpretation of QS data without requiring any assumption about the motional model.

Order parameters of a transmembrane helix in a fluid bilayer: case study of a WALP peptide

A new solid-state NMR-based strategy is established for the precise and efficient analysis of orientation and dynamics of transmembrane peptides in fluid bilayers. For this purpose, several dynamically averaged anisotropic constraints, including (13)C and (15)N chemical shift anisotropies and (13)C-(15)N dipolar couplings, were determined from two different triple-isotope-labeled WALP23 peptides ((2)H, (13)C, and (15)N) and combined with previously published quadrupolar splittings of the same peptide. Chemical shift anisotropy tensor orientations were determined with quantum chemistry. The complete set of experimental constraints was analyzed using a generalized, four-parameter dynamic model of the peptide motion, including tilt and rotation angle and two associated order parameters. A tilt angle of 21 degrees was determined for WALP23 in dimyristoylphosphatidylcholine, which is much larger than the tilt angle of 5.5 degrees previously determined from (2)H NMR experiments. This approach provided a realistic value for the tilt angle of WALP23 peptide in the presence of hydrophobic mismatch, and can be applied to any transmembrane helical peptide. The influence of the experimental data set on the solution space is discussed, as are potential sources of error.

Changes in transmembrane helix alignment by arginine residues revealed by solid-state NMR experiments and coarse-grained MD simulations

Independent experimental and computational approaches show agreement concerning arginine/membrane interactions when a single arginine is introduced at selected positions within the membrane-spanning region of acetyl-GGALW(5)LALALAL(12)AL(14)ALALW(19)LAGA-ethanolamide, designated GWALP23. Peptide sequence isomers having Arg in position 12 or position 14 display markedly different behaviors, as deduced by both solid-state NMR experiments and coarse-grained molecular dynamics (CG-MD) simulations. With respect to the membrane normal of DOPC or DPPC lipid bilayer membranes, GWALP23-R14 shows one major state whose apparent average tilt is approximately 10 degrees greater than that of GWALP23. The presence of R14 furthermore induces bilayer thinning and peptide displacement to "lift" the charged guanidinium toward the bilayer surface. By contrast, GWALP23-R12 exhibits multiple states that are in slow exchange on the NMR time scale, with CG-MD simulations indicating two distinct positions with different screw rotation angles in the membrane, along with an increased tendency to exit the lipid bilayer.

Sensitivity enhancement of separated local field experiments: application to membrane proteins

Separated local field (SLF) experiments have been used for almost three decades to obtain structural information in solid-state NMR. These experiments resolve chemical shift anisotropy (CSA) from dipole-dipole interactions (dipolar couplings, DC) in isolated spin systems. Both CSA and DC data can be converted into orientational constraints to elucidate the secondary structure and topology of membrane proteins in oriented lipid bilayers. Here, we propose a new suite of sensitivity enhanced SLF pulse sequences to measure CSA and DC for aligned membrane proteins and liquid crystalline molecules that will decrease the time needed for data acquisition. We demonstrate the efficacy of these new sensitivity enhanced experiments using both a single crystal of N-acetyl leucine and a single pass membrane protein sarcolipin reconstituted in aligned lipid bicelles. These results lay the groundwork for the routine application of this methodology for studying the structure and topology of membrane proteins.

Influence of proline upon the folding and geometry of the WALP19 transmembrane peptide

The orientations, geometries, and lipid interactions of designed transmembrane (TM) peptides have attracted significant experimental and theoretical interest. Because the amino acid proline will introduce a known discontinuity into an alpha helix, we have sought to measure the extent of helix kinking caused by a single proline within the isolated TM helical domain of WALP19. For this purpose, we synthesized acetyl-GWWLALALAP(10)ALALALWWA-ethanolamide and included pairs of deuterated alanines by using 60-100% Fmoc-l-Ala-d(4) at selected sequence positions. Solid-state deuterium ((2)H) magnetic resonance spectra from oriented, hydrated samples (1/40, peptide/lipid; using several lipids) reveal signals from many of the alanine backbone C(alpha) deuterons as well as the alanine side-chain C(beta) methyl groups, whereas signals from C(alpha) deuterons generally have not been observed for similar peptides without proline. It is conceivable that altered peptide dynamics may be responsible for the apparent "unmasking" of the backbone resonances in the presence of the proline. Data analysis using the geometric analysis of labeled alanines (GALA) method reveals that the peptide helix is distorted due to the presence of the proline. To provide additional data points for evaluating the segmental tilt angles of the two halves of the peptide, we substituted selected leucines with l-Ala-d(4). Using this approach, we were able to deduce that the apparent average tilt of the C-terminal increases from approximately 4 degrees to approximately 12 degrees when Pro(10) is introduced. The segment N-terminal to proline is more complex and possibly is more dynamically flexible; Leu to Ala mutations within the N-terminal segment alter the average orientations of alanines in both segments. Nevertheless, in DOPC, we could estimate an apparent kink angle of approximately 19 degrees . Together, the results suggest that the central proline influences not only the geometry but also the dynamics of the membrane-spanning peptide. The results make up an important basis for understanding the functional role of proline in several families of membrane proteins.

Tilt and rotation angles of a transmembrane model peptide as studied by fluorescence spectroscopy

In this study the membrane orientation of a tryptophan-flanked model peptide, WALP23, was determined by using peptides that were labeled at different positions along the sequence with the environmentally sensitive fluorescent label BADAN. The fluorescence properties, reflecting the local polarity, were used to determine the tilt and rotation angles of the peptide based on an ideal alpha-helix model. For WALP23 inserted in dioleoylphosphatidylcholine (DOPC), an estimated tilt angle of the helix with respect to the bilayer normal of 24 degrees +/- 5 degrees was obtained. When the peptides were inserted into bilayers with different acyl chain lengths or containing different concentrations of cholesterol, small changes in tilt angle were observed as response to hydrophobic mismatch, whereas the rotation angle appeared to be independent of lipid composition. In all cases, the tilt angles were significantly larger than those previously determined from (2)H NMR experiments, supporting recent suggestions that the relatively long timescale of (2)H NMR measurements may result in an underestimation of tilt angles due to partial motional averaging. It is concluded that although the fluorescence technique has a rather low resolution and limited accuracy, it can be used to resolve the discrepancies observed between previous (2)H NMR experiments and molecular-dynamics simulations.

Orientation and dynamics of transmembrane peptides: the power of simple models

In this review we discuss recent insights obtained from well-characterized model systems into the factors that determine the orientation and tilt angles of transmembrane peptides in lipid bilayers. We will compare tilt angles of synthetic peptides with those of natural peptides and proteins, and we will discuss how tilt can be modulated by hydrophobic mismatch between the thickness of the bilayer and the length of the membrane spanning part of the peptide or protein. In particular, we will focus on results obtained on tryptophan-flanked model peptides (WALP peptides) as a case study to illustrate possible consequences of hydrophobic mismatch in molecular detail and to highlight the importance of peptide dynamics for the experimental determination of tilt angles. We will conclude with discussing some future prospects and challenges concerning the use of simple peptide/lipid model systems as a tool to understand membrane structure and function.

Novel free energy calculations to explore mechanisms and energetics of membrane protein structure and function

Understanding the delicate balance of forces governing helix or beta-hairpin interactions in transmembrane (TM) proteins is central to understanding membrane structure and function. These membrane constituent interactions play an essential role in determining the structure and function of membrane proteins, and protein interactions in membranes, and thus form the basis for many vital processes, including TM signaling, transport of ions and small molecules, energy transduction, and cell-cell recognition. "Why does a single-pass TM helix or beta-hairpin have specific orientations in membranes?" "What are the roles of hydrogen bonds, close packing, and helix-lipid or beta-hairpin-lipid interactions in helix or beta-hairpin associations in membranes?" "How do these interactions change the membrane structures?" "How do TM domains transmit signals across membranes?" These are important membrane biophysical questions that can be addressed by understanding the delicate balance of forces governing helix or beta-hairpin interactions with/in membranes. In this work, we summarize a series of helix/beta-hairpin restraint potentials that we have developed, and illustrate their applications that begin to address the complicated energetics and molecular mechanisms of these interactions at the atomic level by calculating the potentials of mean force (PMFs) along reaction coordinates relevant to helix/beta-hairpin motions in membranes and dissecting the total PMF into the contributions arising from physically important microscopic forces.

Solid state NMR analysis of peptides in membranes: Influence of dynamics and labeling scheme

The functional state of a membrane-active peptide is often defined by its conformation, molecular orientation, and its oligomeric state in the lipid bilayer. These "static" structural properties can be routinely studied by solid state NMR using isotope-labeled peptides. In the highly dynamic environment of a liquid crystalline biomembrane, however, the whole-body fluctuations of a peptide are also of paramount importance, although difficult to address and most often ignored. Yet it turns out that disregarding such motional averaging in calculating the molecular alignment from orientational NMR-constraints may give a misleading, if not false picture of the system. Here, we demonstrate that the reliability of a simplified static or an advanced dynamic data analysis depends critically on the choice of isotope labeling scheme used. Two distinctly different scenarios have to be considered. When the labels are placed on the side chains of a helical peptide (such as a CD(3)- or CF(3)-group attached to the C(alpha)C(beta) bond), their nuclear spin interaction tensors are very sensitive to motional averaging. If this effect is not properly accounted for, the helix tilt angle tends to be severely underestimated. At the same time, the analysis of labels in the side chains allows to extract valuable dynamical information about whole-body fluctuations of the peptide helix in the membrane. On the other hand, the alternative labeling scheme where (15)N-labels are accommodated within the peptide backbone, will yield nearly correct helix tilt angles, irrespective as to whether dynamics are taken into account or not.

Tilt and azimuthal angles of a transmembrane peptide: a comparison between molecular dynamics calculations and solid-state NMR data of sarcolipin in lipid membranes

We report molecular dynamics simulations in the explicit membrane environment of a small membrane-embedded protein, sarcolipin, which regulates the sarcoplasmic reticulum Ca-ATPase activity in both cardiac and skeletal muscle. In its monomeric form, we found that sarcolipin adopts a helical conformation, with a computed average tilt angle of 28 +/- 6 degrees and azymuthal angles of 66 +/- 22 degrees, in reasonable accord with angles determined experimentally (23 +/- 2 degrees and 50 +/- 4 degrees, respectively) using solid-state NMR with separated-local-field experiments. The effects of time and spatial averaging on both (15)N chemical shift anisotropy and (1)H/(15)N dipolar couplings have been analyzed using short-time averages of fast amide out-of-plane motions and following principal component dynamic trajectories. We found that it is possible to reproduce the regular oscillatory patterns observed for the anisotropic NMR parameters (i.e., PISA wheels) employing average amide vectors. This work highlights the role of molecular dynamics simulations as a tool for the analysis and interpretation of solid-state NMR data.

Orientation and dynamics of peptides in membranes calculated from 2H-NMR data

Solid-state (2)H-NMR is routinely used to determine the alignment of membrane-bound peptides. Here we demonstrate that it can also provide a quantitative measure of the fluctuations around the distinct molecular axes. Using several dynamic models with increasing complexity, we reanalyzed published (2)H-NMR data on two representative alpha-helical peptides: 1), the amphiphilic antimicrobial peptide PGLa, which permeabilizes membranes by going from a monomeric surface-bound to a dimeric tilted state and finally inserting as an oligomeric pore; and 2), the hydrophobic WALP23, which is a typical transmembrane segment, although previous analysis had yielded helix tilt angles much smaller than expected from hydrophobic mismatch and molecular dynamics simulations. Their (2)H-NMR data were deconvoluted in terms of the two main helix orientation angles (representing the time-averaged peptide tilt and azimuthal rotation), as well as the amplitudes of fluctuation about the corresponding molecular axes (providing the dynamic picture). The mobility of PGLa is found to be moderate and to correlate well with the respective oligomeric states. WALP23 fluctuates more vigorously, now in better agreement with the molecular dynamics simulations and mismatch predictions. The analysis demonstrates that when (2)H-NMR data are fitted to extract peptide orientation angles, an explicit representation of the peptide rigid-body angular fluctuations should be included.