Comparative analysis of the orientation of transmembrane peptides using solid-state <Superscript>2</Superscript>H- and <Superscript>15</Superscript>N-NMR: mobility matters

Many solid-state nuclear magnetic resonance (NMR) approaches for membrane proteins rely on orientation-dependent parameters, from which the alignment of peptide segments in the lipid bilayer can be calculated. Molecules embedded in liquid-crystalline membranes, such as monomeric helices, are highly mobile, leading to partial averaging of the measured NMR parameters. These dynamic effects need to be taken into account to avoid misinterpretation of NMR data. Here, we compare two common NMR approaches: ²H-NMR quadrupolar waves, and separated local field ¹⁵N–¹H polarization inversion spin exchange at magic angle (PISEMA) spectra, in order to identify their strengths and drawbacks for correctly determining the orientation and mobility of α-helical transmembrane peptides. We first analyzed the model peptide WLP23 in oriented dimyristoylphosphatidylcholine (DMPC) membranes and then contrasted it with published data on GWALP23 in dilauroylphosphatidylcholine (DLPC). We only obtained consistent tilt angles from the two methods when taking dynamics into account. Interestingly, the two related peptides differ fundamentally in their mobility. Although both helices adopt the same tilt in their respective bilayers (~20°), WLP23 undergoes extensive fluctuations in its azimuthal rotation angle, whereas GWALP23 is much less dynamic. Both alternative NMR methods are suitable for characterizing orientation and dynamics, yet they can be optimally used to address different aspects. PISEMA spectra immediately reveal the presence of large-amplitude rotational fluctuations, which are not directly seen by ²H-NMR. On the other hand, PISEMA was unable to define the azimuthal rotation angle in the case of the highly dynamic WLP23, though the helix tilt could still be determined, irrespective of any dynamics parameters.     

Orientation and Dynamics of Peptides in Membranes Calculated from H-NMR Data

Solid-state ²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 ²H-NMR data on two representative α-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 ²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 ²H-NMR data are fitted to extract peptide orientation angles, an explicit representation of the peptide rigid-body angular fluctuations should be included.     

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.     

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 α-helix model. For WALP23 inserted in dioleoylphosphatidylcholine (DOPC), an estimated tilt angle of the helix with respect to the bilayer normal of 24° ± 5° 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 ²H NMR experiments, supporting recent suggestions that the relatively long timescale of ²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 ²H NMR experiments and molecular-dynamics simulations.     

Molecular dynamics simulations of the anchoring and tilting of the lung-surfactant peptide SP-B1-25 in palmitic acid monolayers

We have performed molecular dynamics simulations of multiple copies of the lung-surfactant peptide SP-B1-25 in a palmitic acid (PA) monolayer. SP-B1-25 is a shorter version of lung-surfactant protein B, an important component of lung surfactant. Up to 30 ns simulations of 20 wt % SP-B1-25 in the PA monolayers were performed with different surface areas of PA, extents of PA ionization, and various initial configurations of the peptides. Starting with initial peptide orientation perpendicular to the monolayer, the predicted final tilt angles average 54 degrees approximately 62 degrees with respect to the monolayer normal, similar to those measured experimentally by Lee et al. (Biophysical Journal. 2001. Synchrotron x-ray study of lung surfactant-specific protein SP-B in lipid monolayers. 81:572-585). In their final conformations, hydrogen-bond analysis and amino acid mutation studies show that the peptides are anchored by hydrogen bond interactions between the cationic residues Arg-12 and Arg-17 and the hydrogen bond acceptors of the ionized PA headgroup, and the tilt angle is affected by the interactions of Tyr-7 and Gln-19 with the PA headgroup. Our work indicates that the factors controlling orientation of small peptides in lipid layers can now be uncovered through molecular dynamics simulations.     

Interpretation of H-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 ²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.     

Structure and dynamics of an amphiphilic peptide in a lipid bilayer: a molecular dynamics study

A molecular dynamics simulation of a simple model membrane system composed of a single amphiphilic helical peptide (ace-K2GL16K2A-amide) in a fully hydrated 1,2-dimyristoyl-sn-glycero-3-phosphocholine bilayer was performed for a total of 1060 ps. The secondary structure of the peptide and its stability were described in terms of average dihedral angles, phi and psi, and the C alpha torsion angles formed by backbone atoms; by the average translation per residue along the helix axis; and by the intramolecular peptide hydrogen bonds. The results indicated that residues 6 through 15 remain in a stable right-handed alpha-helical conformation, whereas both termini exhibit substantial fluctuations. A change in the backbone dihedral angles for residues 16 and 17 is accompanied by the loss of two intramolecular hydrogen bonds, leading to a local but long-lived disruption of the helix. The dynamics of the peptide was characterized in terms of local and global helix motions. The local motions of the N-H bond angles were described in terms of the autocorrelation functions of P2[cos thetaNH(t, t + tau)] and reflected the different degrees of local peptide order as well as a variation in time scale for local motions. The chi1 and chi2 dihedral angles of the leucine side chains underwent frequent transitions between potential minima. No connection between the side-chain positions and their mobility was observed, however. In contrast, the lysine side chains displayed little mobility during the simulation. The global peptide motions were characterized by the tilting and bending motions of the helix. Although the peptide was initially aligned parallel to the bilayer normal, during the simulation it was observed to tilt away from the normal, reaching an angle of approximately 25 degrees by the end of the simulation. In addition, a slight bend of the helix was detected. Finally, the solvation of the peptide backbone and side-chain atoms was also investigated.     

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₃- or CF₃-group attached to the C^αC^β 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 ¹⁵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.     

Orientation of a beta-hairpin antimicrobial peptide in lipid bilayers from two-dimensional dipolar chemical-shift correlation NMR

The orientation of a beta-sheet membrane peptide in lipid bilayers is determined, for the first time, using two-dimensional (2D) (15)N solid-state NMR. Retrocyclin-2 is a disulfide-stabilized cyclic beta-hairpin peptide with antibacterial and antiviral activities. We used 2D separated local field spectroscopy correlating (15)N-(1)H dipolar coupling with (15)N chemical shift to determine the orientation of multiply (15)N-labeled retrocyclin-2 in uniaxially aligned phosphocholine bilayers. Calculated 2D spectra exhibit characteristic resonance patterns that are sensitive to both the tilt of the beta-strand axis and the rotation of the beta-sheet plane from the bilayer normal and that yield resonance assignment without the need for singly labeled samples. Retrocyclin-2 adopts a transmembrane orientation in dilauroylphosphatidylcholine bilayers, with the strand axis tilted at 20 degrees +/- 10 degrees from the bilayer normal, but changes to a more in-plane orientation in thicker 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidyl-choline (POPC) bilayers with a tilt angle of 65 degrees +/- 15 degrees . These indicate that hydrophobic mismatch regulates the peptide orientation. The 2D spectra are sensitive not only to the peptide orientation but also to its backbone (phi, psi) angles. Neither a bent hairpin conformation, which is populated in solution, nor an ideal beta-hairpin with uniform (phi, psi) angles and coplanar strands, agrees with the experimental spectrum. Thus, membrane binding orders the retrocyclin conformation by reducing the beta-sheet curvature but does not make it ideal. (31)P NMR spectra of lipid bilayers with different compositions indicate that retrocyclin-2 selectively disrupts the orientational order of anionic membranes while leaving zwitteronic membranes intact. These structural results provide insights into the mechanism of action of this beta-hairpin antimicrobial peptide.     

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

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

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/¹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.     

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/¹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.     

Assessing the effects of time and spatial averaging in 15N chemical shift/15N-1H dipolar correlation solid state NMR experiments

The effect of time and spatial averaging on ¹⁵N chemical shift/¹H-¹⁵N dipolar correlation spectra, i.e., PISEMA spectra, of -helical membrane peptides and proteins is investigated. Three types of motion are considered: (a) Librational motion of the peptide planes in the -helix; (b) rotation of the helix about its long axis; and (c) wobble of the helix about a nominal tilt angle. A 2ns molecular dynamics simulation of helix D of bacteriorhodopsin is used to determine the effect of librational motion on the spectral parameters. For the time averaging, the rotation and wobble of this same helix are modelled by assuming either Gaussian motion about the respective angles or a uniform distribution of a given width. For the spatial averaging, regions of possible ¹⁵N chemical shift/¹H-¹⁵N dipolar splittings are computed for a distribution of rotations and/or tilt angles of the helix. The computed spectra show that under certain motional modes the ¹⁵N chemical shift/¹H-¹⁵N dipolar pairs for each of the residues do not form patterns which mimic helical wheel patterns. As a result, the unambiguous identification of helix tilt and helix rotation without any resonance assignments or on the basis of a single assignment may be difficult.     

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.     

Peptide insertion, positioning, and stabilization in a membrane: insight from an all-atom molecular dynamics simulation

Peptide insertion, positioning, and stabilization in a model membrane are probed via an all-atom molecular dynamics (MD) simulation. One peptide (WL5) is simulated in each leaflet of a solvated dimyristoylglycero-3-phosphate (DMPC) membrane. Within the first 5 ns, the peptides spontaneously insert into the membrane and then stabilize during the remaining 70 ns of simulation time. In both leaflets, the peptides localize to the membrane interface, and this localization is attributed to the formation of peptide-lipid hydrogen bonds. We show that the single tryptophan residue in each peptide contributes significantly to these hydrogen bonds; specifically, the nitrogen heteroatom of the indole ring plays a critical role. The tilt angles of the indole rings relative to the membrane normal in the upper and lower leaflets are approximately 26 degrees and 54 degrees , respectively. The tilt angles of the entire peptide chain are 62 degrees and 74 degrees . The membrane induces conformations of the peptide that are characteristic of beta-sheets, and the peptide enhances the lipid ordering in the membrane. Finally, the diffusion rate of the peptides in the membrane plane is calculated (based on experimental peptide concentrations) to be approximately 6 A(2)/ns, thus suggesting a 500 ns time scale for intermolecular interactions.     

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

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

Structural dynamics and topology of phospholamban in oriented lipid bilayers using multidimensional solid-state NMR

Phospholamban (PLN), a single-pass membrane protein, regulates heart muscle contraction and relaxation by reversible inhibition of the sarco(endo)plasmic reticulum Ca-ATPase (SERCA). Studies in detergent micelles and oriented lipid bilayers have shown that in its monomeric form PLN adopts a dynamic L shape (bent or T state) that is in conformational equilibrium with a more dynamic R state. In this paper, we use solid-state NMR on both uniformly and selectively labeled PLN to refine our initial studies, describing the topology and dynamics of PLN in oriented lipid bilayers. Two-dimensional PISEMA (polarization inversion spin exchange at the magic angle) experiments carried out in DOPC/DOPE mixed lipid bilayers reveal a tilt angle of the transmembrane domain with respect to the static magnetic field, of 21 +/- 2 degrees and, at the same time, map the rotation angle of the transmembrane domain with respect to the bilayer. PISEMA spectra obtained with selectively labeled samples show that the cytoplasmic domain of PLN is helical and makes an angle of 93 +/- 6 degrees with respect to the bilayer normal. In addition, using samples tilted by 90 degrees , we find that the transmembrane domain of PLN undergoes fast long-axial rotational diffusion about the bilayer normal with the cytoplasmic domain undergoing this motion and other complex dynamics, scaling the values of chemical shift anisotropy. While this dynamic was anticipated by previous solution NMR relaxation studies in micelles, these measurements in the anisotropic lipid environment reveal new dynamic and conformational features encoded in the free protein that might be crucial for SERCA recognition and subsequent inhibition.