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

Simulation of oriented NMR spectra: Combining molecular dynamics and chemical shift tensor calculations

Solid state NMR is widely used to study the orientation and other structural features of proteins and peptides in lipid bilayers. Using data obtained by PISEMA (Polarization Inversion Spin Exchange at Magic Angle) experiments, periodic spectral patterns arise from well-aligned α-helical molecules. Significant problems in the interpretation of PISEMA spectra may arise for systems that do not form perfectly defined secondary structures, like α-helices, or the signal pattern is disturbed by molecular motion. Here, we present a new method that combines molecular dynamics simulation with tensorial orientational constraints (MDOC) and chemical shift tensor calculations for the simulation and interpretation of PISEMA-like spectra. The calculations include the spectra arising from non α-helical molecules and molecules with non-uniform intrinsic mobility. In a first step, dipolar or quadrupolar interaction tensors drive molecular rotations and reorientations to obtain the proper mean values as observed in corresponding NMR experiments. In a second step, the coordinate snapshots of the MDOC simulations are geometry optimized with the isotropic 15 N chemical shifts as constraints using Bond Polarization Theory (BPT) to provide reliable 15 N CS tensor data. The averaged dipolar 1 H-15 N couplings and the δzz tensor components can then be combined to simulate PISEMA patterns. We apply this method to the ß-helical peptide gramicidin A (gA) and demonstrate that this method enables the assignment of most PISEMA resonances. In addition, MDOC simulations provide local order parameters for the calculated sites. These local order parameters reveal large differences in backbone mobility between L- and D-amino acids of gA.

Implicit model to capture electrostatic features of membrane environment

Membrane protein structure prediction and design are challenging due to the complexity of capturing the interactions in the lipid layer, such as those arising from electrostatics. Accurately capturing electrostatic energies in the low-dielectric membrane often requires expensive Poisson-Boltzmann calculations that are not scalable for membrane protein structure prediction and design. In this work, we have developed a fast-to-compute implicit energy function that considers the realistic characteristics of different lipid bilayers, making design calculations tractable. This method captures the impact of the lipid head group using a mean-field-based approach and uses a depth-dependent dielectric constant to characterize the membrane environment. This energy function Franklin2023 (F23) is built upon Franklin2019 (F19), which is based on experimentally derived hydrophobicity scales in the membrane bilayer. We evaluated the performance of F23 on five different tests probing (1) protein orientation in the bilayer, (2) stability, and (3) sequence recovery. Relative to F19, F23 has improved the calculation of the tilt angle of membrane proteins for 90% of WALP peptides, 15% of TM-peptides, and 25% of the adsorbed peptides. The performances for stability and design tests were equivalent for F19 and F23. The speed and calibration of the implicit model will help F23 access biophysical phenomena at long time and length scales and accelerate the membrane protein design pipeline.

Alteration of Average Thickness of Lipid Bilayer by Membrane-Deforming Inclusions

Thickness of lipid bilayer membranes is a key physical parameter determining membrane permeability and stability with respect to formation of through pores. Most membrane inclusions or impurities like amphipathic peptides, transmembrane peptides, lipid inclusions of a different molecular shape, lipid domains, and protein-lipid domains, locally deform the membrane. The detailed structure of the locally deformed region of the membrane is a kind of "fingerprint" for the inclusion type. However, most experimental methods allow determining only averaged parameters of membranes with incorporated inclusions, thus preventing the direct obtaining of the characteristics of the inclusion. Here we developed a model that allows the obtaining of characteristic parameters of three types of membrane inclusions (amphipathic peptides, transmembrane peptides, monolayer lipid patches) from experimentally observable dependencies of the average thickness of lipid bilayer on the surface concentration of the inclusions. In the case of amphipathic peptides, the model provided the peptide parameters that were in qualitative agreement with the available experimental data.

Implicit model to capture electrostatic features of membrane environment

Membrane protein structure prediction and design are challenging due to the complexity of capturing the interactions in the lipid layer, such as those arising from electrostatics. Accurately capturing electrostatic energies in the low-dielectric membrane often requires expensive Poisson-Boltzmann calculations that are not scalable for membrane protein structure prediction and design. In this work, we have developed a fast-to-compute implicit energy function that considers the realistic characteristics of different lipid bilayers, making design calculations tractable. This method captures the impact of the lipid head group using a mean-field-based approach and uses a depth-dependent dielectric constant to characterize the membrane environment. This energy function Franklin2023 (F23) is built upon Franklin2019 (F19), which is based on experimentally derived hydrophobicity scales in the membrane bilayer. We evaluated the performance of F23 on five different tests probing (1) protein orientation in the bilayer, (2) stability, and (3) sequence recovery. Relative to F19, F23 has improved the calculation of the tilt angle of membrane proteins for 90% of WALP peptides, 15% of TM-peptides, and 25% of the adsorbed peptides. The performances for stability and design tests were equivalent for F19 and F23. The speed and calibration of the implicit model will help F23 access biophysical phenomena at long time and length scales and accelerate the membrane protein design pipeline.

Effect of lipid saturation on the topology and oligomeric state of helical membrane polypeptides

Natural liquid crystalline membranes are made up of many different lipids carrying a mixture of saturated and unsaturated fatty acyl chains. Whereas in the past considerable attention has been paid to cholesterol content, the phospholipid head groups and the membrane surface charge the detailed fatty acyl composition was often considered less important. However, recent investigations indicate that the detailed fatty acyl chain composition has pronounced effects on the oligomerization of the transmembrane helical anchoring domains of the MHC II receptor or the membrane alignment of the cationic antimicrobial peptide PGLa. In contrast the antimicrobial peptides magainin 2 and alamethicin are less susceptible to lipid saturation. Using histidine-rich LAH4 designer peptides the high energetic contributions of lipid saturation in stabilizing transmembrane helical alignments are quantitatively evaluated. These observations can have important implications for the biological regulation of membrane proteins and should be taken into considerations during biophysical or structural experiments.

Azide- and diazirine-modified membrane lipids: Physicochemistry and applicability to study peptide/lipid interactions via cross-linking/mass spectrometry

Although the incorporation of photo-activatable lipids into membranes potentially opens new avenues for studying interactions with peptides and proteins, the question of whether azide- or diazirine-modified lipids are suitable for such studies remains controversial. We have recently shown that diazirine-modified lipids can indeed form cross-links to membrane peptides after UV activation and that these cross-links can be precisely determined in their position by mass spectrometry (MS). However, we also observed an unexpected backfolding of the lipid's diazirine-containing stearoyl chain to the membrane interface challenging the potential application of this modified lipid for future cross-linking (XL)-MS studies of protein/lipid interactions. In this work, we compared an azide- (AzidoPC) and a diazirine-modified (DiazPC) membrane lipid regarding their self-assembly properties, their mixing behavior with saturated bilayer-forming phospholipids, and their reactivity upon UV activation using differential scanning calorimetry (DSC), dynamic light scattering (DLS), small-angle X-ray scattering (SAXS), transmission electron microscopy (TEM), and MS. Mixtures of both modified lipids with DMPC were further used for photo-chemically induced XL experiments with a transmembrane model peptide (KLAW23) to elucidate similarities and differences between the azide and the diazirine moiety. We showed that both photo-reactive lipids can be used to study lipid/peptide and lipid/protein interactions. The AzidoPC proved easier to handle, whereas the DiazPC had fewer degradation products and a higher cross-linking yield. However, the problem of backfolding occurs in both lipids; thus, it seems to be a general phenomenon.

A Diazirine-Modified Membrane Lipid to Study Peptide/Lipid Interactions - Chances and Challenges

Although incorporation of photo‐activatable lipids into membranes potentially opens up novel avenues for investigating interactions with proteins, the question of whether diazirine‐modified lipids are suitable for such studies, remains under debate. Focusing on the potential for studying lipid/peptide interactions by cross‐linking mass spectrometry (XL‐MS), we developed a diazirine‐modified lipid (DiazPC), and examined its behaviour in membranes incorporating the model α‐helical peptide LAVA20. We observed an unexpected backfolding of the diazirine‐containing stearoyl chain of the lipid. This surprising behaviour challenges the potential application of DiazPC for future XL‐MS studies of peptide and protein/lipid interactions. The observations made for DiazPC most likely represent a general phenomenon for any type of membrane lipids with a polar moiety incorporated into the alkyl chain. Our finding is therefore of importance for future protein/lipid interaction studies relying on modified lipid probes.

A theoretical assessment of structure determination of multi-span membrane proteins by oriented sample solid-state NMR spectroscopy

Oriented sample solid state NMR (OS-ssNMR) spectroscopy allows direct determination of the structure and topology of membrane proteins reconstituted into aligned lipid bilayers. While OS-ssNMR theoretically has no upper size limit, its application to multi-span membrane proteins has not been established since most studies have been restricted to single or dual span proteins and peptides. Here, we present a critical assessment of the application of this method to multi-span membrane proteins. We used molecular dynamics simulations to back-calculate [¹⁵N-¹H] separated local field (SLF) spectra from a G protein-coupled receptor (GPCR) and show that fully resolved spectra can be obtained theoretically for a multi-span membrane protein with currently achievable resonance linewidths.

Structure, Topology, and Dynamics of Membrane-Inserted Polypeptides and Lipids by Solid-State NMR Spectroscopy: Investigations of the Transmembrane Domains of the DQ Beta-1 Subunit of the MHC II Receptor and of the COP I Protein p24

MHC class II receptors carry important function in adaptive immunity and their malfunctioning is associated with diabetes type I, chronic inflammatory diseases and other autoimmune diseases. The protein assembles from the DQ alpha-1 and DQ beta-1 subunits where the transmembrane domains of these type I membrane proteins have been shown to be involved in homo- and heterodimer formation. Furthermore, the DQ alpha 1 chain carries a sequence motif that has been first identified in the context of p24, a protein involved in the formation of COPI vesicles of the intracellular transport machinery, to specifically interact with sphingomyelin-C18 (SM-C18). Here we investigated the membrane interactions and dynamics of DQ beta-1 in liquid crystalline POPC phospholipid bilayers by oriented ¹⁵N solid-state NMR spectroscopy. The ¹⁵N resonances are indicative of a helical tilt angle of the membrane anchor sequence around 20°. Two populations can be distinguished by their differential dynamics probably corresponding the DQ beta-1 mono- and homodimer. Whereas, this equilibrium is hardly affected by the addition of 5 mole% SM-C18 a single population is visible in DMPC lipid bilayers suggesting that the lipid saturation is an important parameter. Furthermore, the DQ alpha-1, DQ beta-1 and p24 transmembrane helical domains were reconstituted into POPC or POPC/SM-C18 lipid bilayers where the fatty acyl chain of either the phosphatidylcholine or of the sphingolipid have been deuterated. Interestingly in the presence of both sphingolipid and polypeptide a strong decrease in the innermost membrane order of the POPC palmitoyl chain is observed, an effect that is strongest for DQ beta-1. In contrast, for the first time the polypeptide interactions were monitored by deuteration of the stearoyl chain of SM-C18. The resulting ²H solid-state NMR spectra show an increase in order for p24 and DQ alpha-1 which both carry the SM recognition motif. Thereby the data are suggestive that SM-C18 together with the transmembrane domains form structures imposing positive curvature strain on the surrounding POPC lipids. This effect is attenuated when SM-C18 is recognized by the protein.

Solid-State NMR Investigations of the MHC II Transmembrane Domains: Topological Equilibria and Lipid Interactions

The major histocompatibility complex class II (MHC II) membrane proteins are key players in the adaptive immune response. An aberrant function of these molecules is associated with a large number of autoimmune diseases such as diabetes type I and chronic inflammatory diseases. The MHC class II is assembled from DQ alpha 1 and DQ beta 1 which come together as a heterodimer through GXXXG-mediated protein-protein interactions and a highly specific protein-sphingomyelin-C18 interaction motif located on DQA1. This association can have important consequences in regulating the function of these membrane proteins. Here, we investigated the structure and topology of the DQA1 and DQB1 transmembrane helical domains by CD-, oriented ²H and ¹⁵N solid-state NMR spectroscopies. The spectra at peptide-to-lipid ratios of 0.5 to 2 mol% are indicative of a topological equilibrium involving a helix crossing the membrane with a tilt angle of about 20° and another transmembrane topology with around 30° tilt. The latter is probably representing a dimer. Furthermore, at the lowest peptide-to-lipid ratio, a third polypeptide population becomes obvious. Interestingly, the DQB1 and to a lesser extent the DQA1 transmembrane helical domains exhibit a strong fatty acyl chain disordering effect on the inner segments of the ²H-labelled palmitoyl chain of POPC bilayers. This phosphatidylcholine disordering requires the presence of sphingomyelin-C18 suggesting that the ensemble of transmembrane polypeptide and sphingolipid exerts positive curvature strain.

Amide I SFG Spectral Line Width Probes the Lipid-Peptide and Peptide-Peptide Interactions at Cell Membrane In Situ and in Real Time

The balance of lipid-peptide and peptide-peptide interactions at cell membrane is essential to a large variety of cellular processes. In this study, we have experimentally demonstrated for the first time that sum frequency generation vibrational spectroscopy can be used to probe the peptide-peptide and lipid-peptide interactions in cell membrane in situ and in real time by determination of the line width of amide I band of protein backbone. Using a "benchmark" model of α-helical WALP23, it is found that the dominated lipid-peptide interaction causes a narrow line width of the amide I band, whereas the peptide-peptide interaction can markedly broaden the line width. When WALP23 molecules insert into the lipid bilayer, a quite narrow line width of the amide I band is observed because of the lipid-peptide interaction. In contrast, when the peptide lies down on the bilayer surface, the line width of amide I band becomes very broad owing to the peptide-peptide interaction. In terms of the real-time change in the line width, the transition from peptide-peptide interaction to lipid-peptide interaction is monitored during the insertion of WALP23 into 1,2-dipalmitoyl- sn-glycero-3-phospho-(1'- rac-glycerol) (DPPG) lipid bilayer. The dephasing time of a pure α-helical WALP23 in 1-palmitoyl-2-oleoyl- sn-glycero-3-phospho-(1'- rac-glycerol) and DPPG bilayer is determined to be 2.2 and 0.64 ps, respectively. The peptide-peptide interaction can largely accelerate the dephasing time.

Structural Characterization of the Amyloid Precursor Protein Transmembrane Domain and Its γ-Cleavage Site

Alzheimer's disease is the most common form of dementia that affects about 50 million of sufferers worldwide. A major role for the initiation and progression of Alzheimer's disease has been associated with the amyloid β-peptide (Aβ), which is a protease cleavage product of the amyloid precursor protein. The amyloid precursor protein is an integral membrane protein with a single transmembrane domain. Here, we assessed the structural integrity of the transmembrane domain within oriented phosphatidylcholine lipid bilayers and determined the tilt angle distribution and dynamics of various subdomains using solid-state NMR and attenuated total reflectance Fourier transform infrared spectroscopies. Although the overall secondary structure of the transmembrane domain is α-helical, pronounced conformational and topological heterogeneities were observed for the γ- and, to a lesser extent, the ζ-cleavage site, with pronounced implications for the production of Aβ and related peptides, the development of the disease, and pharmaceutical innovation.

Azide-Modified Membrane Lipids: Synthesis, Properties, and Reactivity

In the present work, we describe the synthesis and the temperature-dependent behavior of photoreactive membrane lipids as well as their capability to study peptide/lipid interactions. The modified phospholipids contain an azide group either in the middle part or at the end of an alkyl chain and also differ in the linkage (ester vs ether) of the second alkyl chain. The temperature-dependent aggregation behavior of the azidolipids was studied using differential scanning calorimetry (DSC), Fourier-transform infrared (FTIR) spectroscopy, and small-angle X-ray scattering (SAXS). Aggregate structures were visualized by stain and cryo transmission electron microscopy (TEM) and were further characterized by dynamic light scattering (DLS). We show that the position of the azide group and the type of linkage of the alkyl chain at the sn-2 position of the glycerol influences the type of aggregates formed as well as their long-term stability: P10AzSPC and r12AzSHPC show the formation of extrudable liposomes, which are stable in size during storage. In contrast, azidolipids that carry a terminal azido moiety either form extrudable liposomes, which show time-dependent vesicle fusion (P15AzPdPC), or self-assemble in large sheet-like, nonextrudable aggregates (r15AzPdHPC) where the lipid molecules are arranged in an interdigitated orientation at temperatures below T_m (L_βI phase). Finally, a P10AzSPC:DMPC mixture was used for photochemically induced cross-linking experiments with a transmembrane peptide (WAL-peptide) to demonstrate the applicability of the azidolipids for the analysis of peptide/lipid interactions. The efficiency of photo-cross-linking was monitored by attenuated total reflection infrared (ATR-IR) spectroscopy and mass spectrometry (MS).

NMR Structural Studies of Antimicrobial Peptides: LPcin Analogs

Lactophoricin (LPcin), a component of proteose peptone (113-135) isolated from bovine milk, is a cationic amphipathic antimicrobial peptide consisting of 23 amino acids. We designed a series of N- or C-terminal truncated variants, mutated analogs, and truncated mutated analogs using peptide-engineering techniques. Then, we selected three LPcin analogs of LPcin-C8 (LPcin-YK1), LPcin-T2WT6W (LPcin-YK2), and LPcin-T2WT6W-C8 (LPcin-YK3), which may have better antimicrobial activities than LPcin, and successfully expressed them in E. coli with high yield. We elucidated the 3D structures and topologies of the three LPcin analogs in membrane environments by conducting NMR structural studies. We investigated the purity of the LPcin analogs and the α-helical secondary structures by performing (1)H-(15)N 2D HSQC and HMQC-NOESY liquid-state NMR spectroscopy using protein-containing micelle samples. We measured the 3D structures and tilt angles in membranes by conducting (15)N 1D and 2D (1)H-(15)N SAMMY type solid-state NMR spectroscopy with an 800 MHz in-house-built (1)H-(15)N double-resonance solid-state NMR probe with a strip-shield coil, using protein-containing large bicelle samples aligned and confirmed by molecular-dynamics simulations. The three LPcin analogs were found to be curved α-helical structures, with tilt angles of 55-75° for normal membrane bilayers, and their enhanced activities may be correlated with these topologies.

Solid-state NMR methods for oriented membrane proteins

Oriented-sample solid-state NMR represents one of few experimental methods capable of characterising the membrane-bound conformation of proteins in the cell membrane. Since the technique was developed 25 years ago, the technique has been applied to study the structure of helix bundle membrane proteins and antimicrobial peptides, characterise protein-lipid interactions, and derive information on dynamics of the membrane anchoring of membrane proteins. We will review the major developments in various aspects of oriented-sample solid-state NMR, including sample-preparation methods, pulse sequences, theory required to interpret the experiments, perspectives for and guidelines to new experiments, and a number of representative applications.

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.

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

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

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

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

Analyzing the effects of hydrophobic mismatch on transmembrane α-helices using tryptophan fluorescence spectroscopy

Hydrophobic matching between transmembrane protein segments and the lipid bilayer in which they are embedded is a significant factor in the behavior and orientation of such transmembrane segments. The condition of hydrophobic mismatch occurs when the hydrophobic thickness of a lipid bilayer is significantly different than the length of the membrane spanning segment of a protein, resulting in a mismatch. This mismatch can result in altered function of proteins as well as nonnative structural arrangements including effects on transmembrane α-helix tilt angles, oligomerization state, and/or the formation of non-transmembrane topographies. Here, a fluorescence-based protocol is described for testing model transmembrane α-helices and their sensitivity to hydrophobic mismatch by measuring the propensity of these helices to form non-transmembrane structures. Overall, good hydrophobic matching between the bilayer and transmembrane segments is an important factor that must be considered when designing membrane proteins or peptides.

Supramolecular structure of membrane-associated polypeptides by combining solid-state NMR and molecular dynamics simulations

Elemental biological functions such as molecular signal transduction are determined by the dynamic interplay between polypeptides and the membrane environment. Determining such supramolecular arrangements poses a significant challenge for classical structural biology methods. We introduce an iterative approach that combines magic-angle spinning solid-state NMR spectroscopy and atomistic molecular dynamics simulations for the determination of the structure and topology of membrane-bound systems with a resolution and level of accuracy difficult to obtain by either method alone. Our study focuses on the Shaker B ball peptide that is representative for rapid N-type inactivating domains of voltage-gated K(+) channels, associated with negatively charged lipid bilayers.

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.

The determinants of hydrophobic mismatch response for transmembrane helices

Hydrophobic mismatch arises from a difference in the hydrophobic thickness of a lipid membrane and a transmembrane protein segment, and is thought to play an important role in the folding, stability and function of membrane proteins. We have investigated the possible adaptations that lipid bilayers and transmembrane α-helices undergo in response to mismatch, using fully-atomistic molecular dynamics simulations totaling 1.4 μs. We have created 25 different tryptophan-alanine-leucine transmembrane α-helical peptide systems, each composed of a hydrophobic alanine-leucine stretch, flanked by 1-4 tryptophan side chains, as well as the β-helical peptide dimer, gramicidin A. Membrane responses to mismatch include changes in local bilayer thickness and lipid order, varying systematically with peptide length. Adding more flanking tryptophan side chains led to an increase in bilayer thinning for negatively mismatched peptides, though it was also associated with a spreading of the bilayer interface. Peptide tilting, bending and stretching were systematic, with tilting dominating the responses, with values of up to ~45° for the most positively mismatched peptides. Peptide responses were modulated by the number of tryptophan side chains due to their anchoring roles and distributions around the helices. Potential of mean force calculations for local membrane thickness changes, helix tilting, bending and stretching revealed that membrane deformation is the least energetically costly of all mismatch responses, except for positively mismatched peptides where helix tilting also contributes substantially. This comparison of energetic driving forces of mismatch responses allows for deeper insight into protein stability and conformational changes in lipid membranes.

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.

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.

Probing the lipid-protein interface using model transmembrane peptides with a covalently linked acyl chain

The aim of this study was to gain insight into how interactions between proteins and lipids in membranes are sensed at the protein-lipid interface. As a probe to analyze this interface, we used deuterium-labeled acyl chains that were covalently linked to a model transmembrane peptide. First, a perdeuterated palmitoyl chain was coupled to the Trp-flanked peptide WALP23 (Ac-CGWW(LA)(8)LWWA-NH(2)), and the deuterium NMR spectrum was analyzed in di-C18:1-phosphatidylcholine (PC) bilayers. We found that the chain order of this peptide-linked chain is rather similar to that of a noncovalently coupled perdeuterated palmitoyl chain, except that it exhibits a slightly lower order. Similar results were obtained when site-specific deuterium labels were used and when the palmitoyl chain was attached to the more-hydrophobic model peptide WLP23 (Ac-CGWWL(17)WWA-NH(2)) or to the Lys-flanked peptide KALP23 (Ac-CGKK(LA)(8)LKKA-NH(2)). The experiments showed that the order of both the peptide-linked chains and the noncovalently coupled palmitoyl chains in the phospholipid bilayer increases in the order KALP23 < WALP23 < WLP23. Furthermore, changes in the bulk lipid bilayer thickness caused by varying the lipid composition from di-C14:1-PC to di-C18:1-PC or by including cholesterol were sensed rather similarly by the covalently coupled chain and the noncovalently coupled palmitoyl chains. The results indicate that the properties of lipids adjacent to transmembrane peptides mostly reflect the properties of the surrounding lipid bilayer, and hence that (at least for the single-span model peptides used in this study) annular lipids do not play a highly specific role in protein-lipid interactions.

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.

Transmembrane helix orientation and dynamics: insights from ensemble dynamics with solid-state NMR observables

As the major component of membrane proteins, transmembrane helices embedded in anisotropic bilayer environments adopt preferential orientations that are characteristic or related to their functional states. Recent developments in solid-state nuclear magnetic resonance (SSNMR) spectroscopy have made it possible to measure NMR observables that can be used to determine such orientations in a native bilayer environment. A quasistatic single conformer model is frequently used to interpret the SSNMR observables, but important motional information can be missing or misinterpreted in the model. In this work, we have investigated the orientation of the single-pass transmembrane domain of viral protein "u" (VpuTM) from HIV-1 by determining an ensemble of structures using multiple conformer models based on the SSNMR ensemble dynamics technique. The resulting structure ensemble shows significantly larger orientational fluctuations while the ensemble-averaged orientation is compatible with the orientation based on the quasistatic model. This observation is further corroborated by comparison with the VpuTM orientation from comparative molecular dynamics simulations in explicit bilayer membranes. SSNMR ensemble dynamics not only reveals the importance of transmembrane helix dynamics in interpretation of SSNMR observables, but also provides a means to simultaneously extract both transmembrane helix orientation and dynamics information from the SSNMR measurements.

Orientational and motional narrowing of solid-state NMR lineshapes of uniaxially aligned membrane proteins

A unified theory for the NMR line shapes of aligned membrane proteins arising from uniaxial disorder (mosaic spread) and global rotational diffusion about the director axis is presented. A superoperator formalism allows one to take into account the effects of continuous radiofrequency irradiation and frequency offsets in the presence of dynamics. A general method based on the Stochastic Liouville Equation makes it possible to bridge the static and dynamic limits in a single model. Simulations of solid-state NMR spectra are performed for a uniform α helix by considering orientational disorder and diffusion of the helix as a whole relative to the alignment axis. The motional narrowing of the resonance lines is highly inhomogeneous and can be used as an additional angular restraint in structure calculations. Experimental solid-state NMR spectra of Pf1 coat protein support the conclusions of the theory for two limiting cases. The static disorder dominates the (15)N NMR spectra of Pf1 aligned on a phage, while fast uniaxial diffusion provides a line narrowing mechanism for the Pf1 protein reconstituted in magnetically aligned bicelles.

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.

Orientation and dynamics of synthetic transbilayer polypeptides containing GpATM dimerization motifs

Deuterium NMR spectroscopy was used to study how the positioning of a dimerization motif within a transbilayer polypeptide influences its orientation and dynamics in bilayers. Three polypeptide variants comprising glycophorin A transmembrane (GpATM) dimerization motifs incorporated into lysine-terminated poly-leucine-alanine helices were mixed into 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine multilamellar vesicles. The variants differed in orientation of the motif segment around the helix axis with respect to the peptide ends. Polypeptides were labeled with methyl-deuterated alanines at positions that were identically situated relative to the peptide ends (Ala-20 and Ala-22) and at two positions within the motif. An analysis of quadrupole splittings revealed similar tilts and orientations of the peptide ends for all three variants, suggesting that average orientations were dominated by interactions at the bilayer surface. For one variant, however, fast orientational fluctuations about the helix axis were significantly smaller. This may indicate some perturbation of peptide dynamics and conformation by interactions that are sensitive to the motif orientation relative to the peptide ends. For the variant that displayed distinct dynamics, one orientation consistent with observed splittings corresponded to the motif being situated such that its two glycines were particularly accessible to adjacent peptides.

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.

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.