Solid-state NMR methods employing (2)H NMR and geometric analysis of labeled alanines (GALA) were used to study the structure and orientation of the transmembrane alpha-helical peptide acetyl-GWW(LA)(8)LWWA-amide (WALP23) in phosphatidylcholine (PC) bilayers of varying thickness. In all lipids the peptide was found to adopt a transmembrane alpha-helical conformation. A small tilt angle of 4.5 degrees was observed in di-18:1-PC, which has a hydrophobic bilayer thickness that approximately matches the hydrophobic length of the peptide. This tilt angle increased slightly but systematically with increasing positive mismatch to 8.2 degrees in di-C12:0-PC, the shortest lipid used. This small increase in tilt angle is insufficient to significantly change the effective hydrophobic length of the peptide and thereby to compensate for the increasing hydrophobic mismatch, suggesting that tilt of these peptides in a lipid bilayer is energetically unfavorable. The tilt and also the orientation around the peptide axis were found to be very similar to the values previously reported for a shorter WALP19 peptide (GWW(LA)(6)LWWA). As also observed in this previous study, the peptide rotates rapidly around the bilayer normal, but not around its helix axis. Here we show that these properties allow application of the GALA method not only to macroscopically aligned samples but also to randomly oriented samples, which has important practical advantages. A minimum of four labeled alanine residues in the hydrophobic transmembrane sequence was found to be required to obtain accurate tilt values using the GALA method.

The orientation of the transmembrane peptide WALP23 under small hydrophobic mismatch has been assessed through long-time-scale molecular dynamics simulations of hundreds of nanoseconds. Each simulation gives systematically large tilt angles (>30 degrees). In addition, the peptide visits various azimuthal rotations that mostly depend on the initial conditions and converge very slowly. In contrast, small tilt angles as well as a well-defined azimuthal rotation were suggested by recent solid-state 2H NMR studies on the same system. To optimally compare our simulations with NMR data, we concatenated the different trajectories in order to increase the sampling. The agreement with 2H NMR quadrupolar splittings is spectacularly better when these latter are back-calculated from the concatenated trajectory than from any individual simulation. From these ensembled-average quadrupolar splittings, we then applied the GALA method as described by Strandberg et al. (Biophys J. 2004, 86, 3709-3721), which basically derives the peptide orientation (tilt and azimuth) from the splittings. We find small tilt angles (6.5 degrees), whereas the real observed tilt in the concatenated trajectory presents a higher value (33.5 degrees). We thus propose that the small tilt angles estimated by the GALA method are the result of averaging effects, provided that the peptide visits many states of different azimuthal rotations. We discuss how to improve the method and suggest some other experiments to confirm this hypothesis. This work also highlights the need to run several and rather long trajectories in order to predict the peptide orientation from computer simulations.