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

Sensing membrane voltage by reorientation of dipolar transmembrane peptides

Membrane voltage plays a vital role in the behavior and functions of the lipid bilayer membrane. For instance, it regulates the exchange of molecules across the membrane through transmembrane proteins such as ion channels. In this paper, we study the membrane voltage-sensing mechanism, which entails the reorientation of α-helices with a change in the membrane voltage. We consider a helix having a large electrical macrodipole embedded in a lipid bilayer as a model system. We performed extensive molecular dynamics simulations to study the effect of variation of membrane voltage on the tilt angle of peptides and ascertain the optimal parameters for designing such a voltage-sensing peptide. A theoretical model for the system is also developed to investigate the interplay of competing effects of hydrophobic mismatch and dipole-electric field coupling on the tilt of the peptide and further explore the parameter space. This work opens the possibility for the design and fabrication of artificial dipolar membrane voltage-sensing elements for biomedical applications.

Site-Specific Fluorescent Labeling of Hemagglutinin-Specific Antigen Binding Fragment through Amine Chemistry Revealed by Mass Spectrometry

To capture the structure of assembled hemagglutinin (HA) nanoparticles at single-particle resolution, HA-specific antigen binding fragments (Fabs) were labeled by fluorescent (FLR) dyes as probes to highlight the HA trimers displayed on the assembled tetravalent HA nanoparticles for a qualitative localization microscopic study. The FLR dyes were conjugated to the Fabs through N-hydroxysuccinimide (NHS) ester mediated amine coupling chemistry. The labeling profile, including labeling ratio, distribution, and site-specific labeling occupancy, can affect the imaging results and introduce inconsistency. To evaluate the labeling profile so as to evaluate the labeling efficiency, a combination of intact mass measurement by MALDI-MS and peptide mapping through LC-MS/MS was implemented. At the intact molecular level, the labeling ratio and distribution were determined. Through peptide mapping, the labeled residues were identified and the corresponding site-specific labeling occupancy was measured. A systematic comparative investigation of four different FLR-labeled 1H01-Fabs (generated from H1 strain HA specific mAb 1H01) allowed accurate profiling of the labeling pattern. The data indicate that the labeling was site-specific and semiquantitative. This warrants the consistency of single-particle fluorescent imaging experiments and allows a further imaging characterization of the single nanoparticles.

Membrane fluidity, composition, and charge affect the activity and selectivity of the AMP ascaphin-8

Ascaphins are cationic antimicrobial peptides that have been shown to have potential in the treatment of infectious diseases caused by multidrug-resistant pathogens (MDR). However, to date, their principal molecular target and mechanism of action are unknown. Results from peptide prediction software and molecular dynamics simulations confirmed that ascaphin-8 is an alpha-helical peptide. For the first time, the peptide was described as membranotrophic using biophysical approaches including calcein liposome leakage, Laurdan general polarization, and dynamic light scattering. Ascaphin-8's activity and selectivity were modulated by rearranging the spatial distribution of lysine (Var-K5), aspartic acid (Var-D4) residues, or substitution of phenylalanine with tyrosine (Var-Y). The parental peptide and its variants presented high affinity toward the bacterial membrane model (≤2 μM), but lost activity in sterol-enriched membranes (mammal and fungal models, with cholesterol and ergosterol, respectively). The peptide-induced pore size was estimated to be >20 nm in the bacterial model, with no difference among peptides. The same pattern was observed in membrane fluidity (general polarization) assays, where all peptides reduced membrane fluidity of the bacterial model but not in the models containing sterols. The peptides also showed high activity toward MDR bacteria. Moreover, peptide sensitivity of the artificial membrane models compared with pathogenic bacterial isolates were in good agreement.

Spiers Memorial Lecture: Analysis and de novo design of membrane-interactive peptides

Membrane-peptide interactions play critical roles in many cellular and organismic functions, including protection from infection, remodeling of membranes, signaling, and ion transport. Peptides interact with membranes in a variety of ways: some associate with membrane surfaces in either intrinsically disordered conformations or well-defined secondary structures. Peptides with sufficient hydrophobicity can also insert vertically as transmembrane monomers, and many associate further into membrane-spanning helical bundles. Indeed, some peptides progress through each of these stages in the process of forming oligomeric bundles. In each case, the structure of the peptide and the membrane represent a delicate balance between peptide-membrane and peptide-peptide interactions. We will review this literature from the perspective of several biologically important systems, including antimicrobial peptides and their mimics, α-synuclein, receptor tyrosine kinases, and ion channels. We also discuss the use of de novo design to construct models to test our understanding of the underlying principles and to provide useful leads for pharmaceutical intervention of diseases.

Characteristic lengths of transmembrane peptides controlling their tilt and lateral distribution between membrane domains

Lipids and proteins of plasma membranes of eukaryotic cells are supposed to form protein-lipid domains, characterized by a different molecular order, bilayer thickness, and elastic parameters. Several mechanisms of preferable distribution of transmembrane proteins to the ordered or disordered membrane domains have been revealed. The mismatch between the length of the protein transmembrane domain and hydrophobic thickness of the lipid bilayer is considered to be an important driving force of protein lateral sorting. Utilizing the continuum theory of elasticity, we analyzed optimal configurations and preferable membrane domains for single-pass transmembrane peptides of various hydrophobic lengths and effective molecular shapes. We obtained that short transmembrane peptides stand perpendicularly to the membrane plane. The exceedance of a certain characteristic length leads to the tilt of the peptide. This length depends on the bilayer thickness. Thus, in the membrane with coexisting ordered (thicker) and disordered (thinner) phases tilting of the peptide in each phase is governed by its individual characteristic length. The lateral distribution of the peptides between ordered and disordered membrane domains is shown to be described by two additional characteristic lengths. The exceedance of the smaller one drives the peptide towards a more ordered and thicker membrane, while the exceedance of the larger characteristic length switches the preferable membrane domain from ordered and thicker to disordered and thinner. Thus, membrane proteins with long enough transmembrane domains are predicted to accumulate in the thinner disordered membrane as compared to the thicker ordered bilayer. For hourglass-like and barrel-like shaped transmembrane peptides the specific regime of sorting was obtained: the peptides distributed almost equally between the phases in a wide range of peptide lengths. This finding allowed explaining the experimental data on lateral distribution of transmembrane peptide tLAT.

Photophysical Properties of BADAN Revealed in the Study of GGBP Structural Transitions

The fluorescent dye BADAN (6-bromoacetyl-2-dimetylaminonaphtalene) is widely used in various fields of life sciences, however, the photophysical properties of BADAN are not fully understood. The study of the spectral properties of BADAN attached to a number of mutant forms of GGBP, as well as changes in its spectral characteristics during structural changes in proteins, allowed to shed light on the photophysical properties of BADAN. It was shown that spectral properties of BADAN are determined by at least one non-fluorescent and two fluorescent isomers with overlapping absorbing bands. It was found that BADAN fluorescence is determined by the unsolvated "PICT" (planar intramolecular charge transfer state) and solvated "TICT" (twisted intramolecular charge transfer state) excited states. While "TICT" state can be formed both as a result of the "PICT" state solvation and as a result of light absorption by the solvated ground state of the dye. BADAN fluorescence linked to GGBP/H152C apoform is quenched by Trp 183, but this effect is inhibited by glucose intercalation. New details of the changes in the spectral characteristics of BADAN during the unfolding of the protein apo and holoforms have been obtained.

Diverse Scientific Benchmarks for Implicit Membrane Energy Functions

Energy functions are fundamental to biomolecular modeling. Their success depends on robust physical formalisms, efficient optimization, and high-resolution data for training and validation. Over the past 20 years, progress in each area has advanced soluble protein energy functions. Yet, energy functions for membrane proteins lag behind due to sparse and low-quality data, leading to overfit tools. To overcome this challenge, we assembled a suite of 12 tests on independent data sets varying in size, diversity, and resolution. The tests probe an energy function's ability to capture membrane protein orientation, stability, sequence, and structure. Here, we present the tests and use the franklin2019 energy function to demonstrate them. We then identify areas for energy function improvement and discuss potential future integration with machine-learning-based optimization methods. The tests are available through the Rosetta Benchmark Server (https://benchmark.graylab.jhu.edu/) and GitHub (https://github.com/rfalford12/Implicit-Membrane-Energy-Function-Benchmark).

Protein Structure Prediction and Design in a Biologically Realistic Implicit Membrane

Protein design is a powerful tool for elucidating mechanisms of function and engineering new therapeutics and nanotechnologies. Although soluble protein design has advanced, membrane protein design remains challenging because of difficulties in modeling the lipid bilayer. In this work, we developed an implicit approach that captures the anisotropic structure, shape of water-filled pores, and nanoscale dimensions of membranes with different lipid compositions. The model improves performance in computational benchmarks against experimental targets, including prediction of protein orientations in the bilayer, ΔΔG calculations, native structure discrimination, and native sequence recovery. When applied to de novo protein design, this approach designs sequences with an amino acid distribution near the native amino acid distribution in membrane proteins, overcoming a critical flaw in previous membrane models that were prone to generating leucine-rich designs. Furthermore, the proteins designed in the new membrane model exhibit native-like features including interfacial aromatic side chains, hydrophobic lengths compatible with bilayer thickness, and polar pores. Our method advances high-resolution membrane protein structure prediction and design toward tackling key biological questions and engineering challenges.

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.

A novel pH-dependent membrane peptide that binds to EphA2 and inhibits cell migration

Misregulation of the signaling axis formed by the receptor tyrosine kinase (RTK) EphA2 and its ligand, ephrinA1, causes aberrant cell-cell contacts that contribute to metastasis. Solid tumors are characterized by an acidic extracellular medium. We intend to take advantage of this tumor feature to design new molecules that specifically target tumors. We created a novel pH-dependent transmembrane peptide, TYPE7, by altering the sequence of the transmembrane domain of EphA2. TYPE7 is highly soluble and interacts with the surface of lipid membranes at neutral pH, while acidity triggers transmembrane insertion. TYPE7 binds to endogenous EphA2 and reduces Akt phosphorylation and cell migration as effectively as ephrinA1. Interestingly, we found large differences in juxtamembrane tyrosine phosphorylation and the extent of EphA2 clustering when comparing TYPE7 with activation by ephrinA1. This work shows that it is possible to design new pH-triggered membrane peptides to activate RTK and gain insights on its activation mechanism.

Refining amino acid hydrophobicity for dynamics simulation of membrane proteins

Coarse-grained (CG) models have been successful in simulating the chemical properties of lipid bilayers, but accurate treatment of membrane proteins and lipid-protein molecular interactions remains a challenge. The CgProt force field, original developed with the multiscale coarse graining method, is assessed by comparing the potentials of mean force for sidechain insertion in a DOPC bilayer to results reported for atomistic molecular dynamics simulations. Reassignment of select CG sidechain sites from the apolar to polar site type was found to improve the attractive interfacial behavior of tyrosine, phenylalanine and asparagine as well as charged lysine and arginine residues. The solvation energy at membrane depths of 0, 1.3 and 1.7 nm correlates with experimental partition coefficients in aqueous mixtures of cyclohexane, octanol and POPC, respectively, for sidechain analogs and Wimley-White peptides. These experimental values serve as important anchor points in choosing between alternate CG models based on their observed permeation profiles, particularly for Arg, Lys and Gln residues where the all-atom OPLS solvation energy does not agree well with experiment. Available partitioning data was also used to reparameterize the representation of the peptide backbone, which needed to be made less attractive for the bilayer hydrophobic core region. The newly developed force field, CgProt 2.4, correctly predicts the global energy minimum in the potentials of mean force for insertion of the uncharged membrane-associated peptides LS3 and WALP23. CgProt will find application in studies of lipid-protein interactions and the conformational properties of diverse membrane protein systems.

Conformational Changes in the Epidermal Growth Factor Receptor: Role of the Transmembrane Domain Investigated by Coarse-Grained MetaDynamics Free Energy Calculations

The epidermal growth factor receptor (EGFR) is a dimeric membrane protein that regulates key aspects of cellular function. Activation of the EGFR is linked to changes in the conformation of the transmembrane (TM) domain, brought about by changes in interactions of the TM helices of the membrane lipid bilayer. Using an advanced computational approach that combines Coarse-Grained molecular dynamics and well-tempered MetaDynamics (CG-MetaD), we characterize the large-scale motions of the TM helices, simulating multiple association and dissociation events between the helices in membrane, thus leading to a free energy landscape of the dimerization process. The lowest energy state of the TM domain is a right-handed dimer structure in which the TM helices interact through the N-terminal small-X₃-small sequence motif. In addition to this state, which is thought to correspond to the active form of the receptor, we have identified further low-energy states that allow us to integrate with a high level of detail a range of previous experimental observations. These conformations may lead to the active state via two possible activation pathways, which involve pivoting and rotational motions of the helices, respectively. Molecular dynamics also reveals correlation between the conformational changes of the TM domains and of the intracellular juxtamembrane domains, paving the way for a comprehensive understanding of EGFR signaling at the cell membrane.

Fluorescence quenching of (dimethylamino)naphthalene dyes Badan and Prodan by tryptophan in cytochromes P450 and micelles

Fluorescence of 2-(N,N-dimethylamino)-6-propionylnaphthalene dyes Badan and Prodan is quenched by tryptophan in Brij 58 micelles as well as in two cytochrome P450 proteins (CYP102, CYP119) with Badan covalently attached to a cysteine residue. Formation of nonemissive complexes between a dye molecule and tryptophan accounts for about 76% of the fluorescence intensity quenching in micelles, the rest is due to diffusive encounters. In the absence of tryptophan, fluorescence of Badan-labeled cytochromes decays with triexponential kinetics characterized by lifetimes of about 100 ps, 700–800 ps, and 3 ns. Site mutation of a histidine residue in the vicinity of the Badan label by tryptophan results in shortening of all three decay lifetimes. The relative amplitude of the fastest component increases at the expense of the two slower ones. The average quenching rate constants are 4.5 × 10⁸ s^–1 (CYP102) and 3.7 × 10⁸ s^–1 (CYP119), at 288 K. Cyclic voltammetry of Prodan in MeCN shows a reversible reduction peak at −1.85 V vs NHE that becomes chemically irreversible and shifts positively upon addition of water. A quasireversible reduction at −0.88 V was observed in an aqueous buffer (pH 7.3). The excited-state reduction potential of Prodan (and Badan) is estimated to vary from about +0.6 V (vs NHE) in polar aprotic media (MeCN) to approximately +1.6 V in water. Tryptophan quenching of Badan/Prodan fluorescence in CYPs and Brij 58 micelles is exergonic by ≤0.5 V and involves tryptophan oxidation by excited Badan/Prodan, coupled with a fast reaction between the reduced dye and water. Photoreduction is a new quenching mechanism for 2-(N,N-dimethylamino)-6-propionylnaphthalene dyes that are often used as solvatochromic polarity probes, FRET donors and acceptors, as well as reporters of solvation dynamics.

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

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

More than the sum of its parts: coarse-grained peptide-lipid interactions from a simple cross-parametrization

Interfacial systems are at the core of fascinating phenomena in many disciplines, such as biochemistry, soft-matter physics, and food science. However, the parametrization of accurate, reliable, and consistent coarse-grained (CG) models for systems at interfaces remains a challenging endeavor. In the present work, we explore to what extent two independently developed solvent-free CG models of peptides and lipids--of different mapping schemes, parametrization methods, target functions, and validation criteria--can be combined by only tuning the cross-interactions. Our results show that the cross-parametrization can reproduce a number of structural properties of membrane peptides (for example, tilt and hydrophobic mismatch), in agreement with existing peptide-lipid CG force fields. We find encouraging results for two challenging biophysical problems: (i) membrane pore formation mediated by the cooperative action of several antimicrobial peptides, and (ii) the insertion and folding of the helix-forming peptide WALP23 in the membrane.

Dual-fluorescence L-amino acid reports insertion and orientation of melittin peptide in cell membranes

Monitoring insertion and orientation of peptides in situ on cell membranes remains a challenge. To this end, we synthesized an l-amino acid (AFaa) containing a dual-fluorescence dye of the 3-hydroxyflavone family, as a side chain. In contrast to other labeling approaches using a flexible linker, the AFaa fluorophore, introduced by solid phase synthesis into desired position of a peptide, is attached closely to its backbone with well-defined orientation, and, therefore, could reflect its localization in the membrane. This concept was validated by replacing the leucine-9 (L9) and tryptophan-19 (W19) residues by AFaa in melittin, a well-studied membrane-active peptide. Due to high sensitivity of AFaa dual emission to the environment polarity, we detected a much deeper insertion of L9 peptide position into the bilayer, compared to the W19 position. Moreover, using fluorescence microscopy with a polarized light excitation, we found different orientation of AFaa at L9 and W19 positions of melittin in the bilayers of giant vesicles and cellular membranes. These results suggested that in the natural membranes, similarly to the model lipid bilayers, melittin is preferentially oriented parallel to the membrane surface. The developed amino acid and the proposed methodology will be of interest to study other membrane peptides.

Reverse micelles as a tool for probing solvent modulation of protein dynamics: Reverse micelle encapsulated hemoglobin

Hydration waters impact protein dynamics. Dissecting the interplay between hydration waters and dynamics requires a protein that manifests a broad range of dynamics. Proteins in reverse micelles (RMs) have promise as tools to achieve this objective because the water content can be manipulated. Hemoglobin is an appropriate tool with which to probe hydration effects. We describe both a protocol for hemoglobin encapsulation in reverse micelles and a facile method using PEG and cosolvents to manipulate water content. Hydration properties are probed using the water-sensitive fluorescence from Hb bound pyranine and covalently attached Badan. Protein dynamics are probed through ligand recombination traces derived from photodissociated carbonmonoxy hemoglobin on a log scale that exposes the potential role of both α and β solvent fluctuations in modulating protein dynamics. The results open the possibility of probing hydration level phenomena in this system using a combination of NMR and optical probes.

Dynamic Heterogeneous Dielectric Generalized Born (DHDGB): An implicit membrane model with a dynamically varying bilayer thickness

An extension to the heterogeneous dielectric generalized Born (HDGB) implicit membrane formalism is presented to allow dynamic membrane deformations in response to membrane-inserted biomolecules during molecular dynamic simulations. The flexible membrane is implemented through additional degrees of freedom that represent the membrane deformation at the contact points of a membrane-inserted solute with the membrane. The extra degrees of freedom determine the dielectric and non-polar solvation free energy profiles that are used to obtain the solvation free energy in the presence of the membrane and are used to calculate membrane deformation free energies according to an elastic membrane model. With the dynamic HDGB (DHDGB) model the membrane is able to deform in response to the insertion of charged molecules thereby avoiding the overestimation of insertion free energies with static membrane models. The DHDGB model also allows the membrane to respond to the insertion of membrane-spanning solutes with hydrophobic mismatch. The model is tested with the membrane insertion of amino acid side chain analogs, arginine-containing helices, the WALP23 peptide, and the gramicidin A channel.

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.

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.

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.

Lipid packing drives the segregation of transmembrane helices into disordered lipid domains in model membranes

Cell membranes are comprised of multicomponent lipid and protein mixtures that exhibit a complex partitioning behavior. Regions of structural and compositional heterogeneity play a major role in the sorting and self-assembly of proteins, and their clustering into higher-order oligomers. Here, we use computer simulations and optical microscopy to study the sorting of transmembrane helices into the liquid-disordered domains of phase-separated model membranes, irrespective of peptide-lipid hydrophobic mismatch. Free energy calculations show that the enthalpic contribution due to the packing of the lipids drives the lateral sorting of the helices. Hydrophobic mismatch regulates the clustering into either small dynamic or large static aggregates. These results reveal important molecular driving forces for the lateral organization and self-assembly of transmembrane helices in heterogeneous model membranes, with implications for the formation of functional protein complexes in real cells.

Influence of hydrophobic mismatch and amino acid composition on the lateral diffusion of transmembrane peptides

We investigated the effect of amino acid composition and hydrophobic length of alpha-helical transmembrane peptides and the role of electrostatic interactions on the lateral diffusion of the peptides in lipid membranes. Model peptides of varying length and composition, and either tryptophans or lysines as flanking residues, were synthesized. The peptides were labeled with the fluorescent label Alexa Fluor 488 and incorporated into phospholipid bilayers of different hydrophobic thickness and composition. Giant unilamellar vesicles were formed by electroformation, and the lateral diffusion of the transmembrane peptides (and lipids) was determined by fluorescence correlation spectroscopy. In addition, we performed coarse-grained molecular-dynamics simulations of single peptides of different hydrophobic lengths embedded in planar membranes of different thicknesses. Both the experimental and simulation results indicate that lateral diffusion is sensitive to membrane thickness between the peptides and surrounding lipids. We did not observe a difference in the lateral diffusion of the peptides with respect to the presence of tryptophans or lysines as flanking residues. The specific lipid headgroup composition of the membrane has a much less pronounced impact on the diffusion of the peptides than does the hydrophobic thickness.

Revisiting hydrophobic mismatch with free energy simulation studies of transmembrane helix tilt and rotation

Protein-lipid interaction and bilayer regulation of membrane protein functions are largely controlled by the hydrophobic match between the transmembrane (TM) domain of membrane proteins and the surrounding lipid bilayer. To systematically characterize responses of a TM helix and lipid adaptations to a hydrophobic mismatch, we have performed a total of 5.8-mus umbrella sampling simulations and calculated the potentials of mean force (PMFs) as a function of TM helix tilt angle under various mismatch conditions. Single-pass TM peptides called WALPn (n = 16, 19, 23, and 27) were used in two lipid bilayers with different hydrophobic thicknesses to consider hydrophobic mismatch caused by either the TM length or the bilayer thickness. In addition, different flanking residues, such as alanine, lysine, and arginine, instead of tryptophan in WALP23 were used to examine their influence. The PMFs, their decomposition, and trajectory analysis demonstrate that 1), tilting of a single-pass TM helix is the major response to a hydrophobic mismatch; 2), TM helix tilting up to approximately 10 degrees is inherent due to the intrinsic entropic contribution arising from helix precession around the membrane normal even under a negative mismatch; 3), the favorable helix-lipid interaction provides additional driving forces for TM helix tilting under a positive mismatch; 4), the minimum-PMF tilt angle is generally located where there is the hydrophobic match and little lipid perturbation; 5), TM helix rotation is dependent on the specific helix-lipid interaction; and 6), anchoring residues at the hydrophilic/hydrophobic interface can be an important determinant of TM helix orientation.

Monitoring membrane binding and insertion of peptides by two-color fluorescent label

Herein, we developed an approach for monitoring membrane binding and insertion of peptides using a fluorescent environment-sensitive label of the 3-hydroxyflavone family. For this purpose, we labeled the N-terminus of three synthetic peptides, melittin, magainin 2 and poly-l-lysine capable to interact with lipid membranes. Binding of these peptides to lipid vesicles induced a strong fluorescence increase, which enabled to quantify the peptide-membrane interaction. Moreover, the dual emission of the label in these peptides correlated well with the depth of its insertion measured by the parallax quenching method. Thus, in melittin and magainin 2, which show deep insertion of their N-terminus, the label presented a dual emission corresponding to a low polar environment, while the environment of the poly-l-lysine N-terminus was rather polar, consistent with its location close to the bilayer surface. Using spectral deconvolution to distinguish the non-hydrated label species from the hydrated ones and two photon fluorescence microscopy to determine the probe orientation in giant vesicles, we found that the non-hydrated species were vertically oriented in the bilayer and constituted the best indicators for evaluating the depth of the peptide N-terminus in membranes. Thus, this label constitutes an interesting new tool for monitoring membrane binding and insertion of peptides.

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.

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.

Profiling of dynamics in protein-lipid-water systems: a time-resolved fluorescence study of a model membrane protein with the label BADAN at specific membrane depths

Profiles of lipid-water bilayer dynamics were determined from picosecond time-resolved fluorescence spectra of membrane-embedded BADAN-labeled M13 coat protein. For this purpose, the protein was labeled at seven key positions. This places the label at well-defined locations from the water phase to the center of the hydrophobic acyl chain region of a phospholipid model membrane, providing us with a nanoscale ruler to map membranes. Analysis of the time-resolved fluorescence spectroscopic data provides the characteristic time constant for the twisting motion of the BADAN label, which is sensitive to the local flexibility of the protein–lipid environment. In addition, we obtain information about the mobility of water molecules at the membrane–water interface. The results provide an unprecedented nanoscale profiling of the dynamics and distribution of water in membrane systems. This information gives clear evidence that the actual barrier of membranes for ions and aqueous solvents is located at the region of carbonyl groups of the acyl chains.