Modeling of peptides and proteins in a membrane environment: I. A solvation model mimicking a lipid bilayer

A theoretical solvation model of peptides and proteins that mimics the heterogeneous membrane-water system was proposed. Our approach is based on the combined use of atomic parameters of solvation for water and hydrocarbons, which approximates the hydrated polar groups and acyl chains of lipids, respectively. This model was tested in simulations of several peptides: a nonpolar 20-mer polyleucine, a hydrophobic peptide with terminal polar groups, and a strongly amphiphilic peptide. The conformational space of the peptides in the presence of the membrane was studied by the Monte Carlo method. Unlike a polar solvent and vacuum, the membrane-like environment was shown to stabilize the α-helical conformation: low-energy structures have a helicity index of 100% in all cases. At the same time, the energetically most favorable orientations of the peptides relative to the membrane depend on their hydrophobic properties: nonpolar polyleucine is entirely immersed in the bilayer and the hydrophobic peptide with polar groups at the termini adopts a transbilayer orientation, whereas the amphiphilic peptide lies at the interface parallel to the membrane plane. The results of the simulations agree well with the available experimental data for these systems. In the following communications of this series, we plan to describe applications of the solvation model to membrane-bound proteins and peptides with biologically important functional activities.     

Modeling of Peptides in Implicit Membrane-Mimetic Media

Abstract

Lipid bilayer plays a crucial role in folding of membrane peptides and their stabilization in the membrane-bound state. Correct treatment of the media effects is thus essential for realistic simulations of peptides in bilayers. Previously (Volynsky et al., 1999), we proposed an efficient solvation model which mimics heterogeneous membrane-water system. The model is based on combined employment of atomic solvation parameters for water and hydrocarbon, which approximate hydrated headgroups and acyl chains of lipids, respectively. In this study, the model is employed in non-restrained Monte Carlo simulations of several peptides: totally apolar 20-residue poly-L-Leu, hydrophobic peptide with polar edges, and strongly amphiphilic pep-tide. The principal goals are: to explore energy landscape of these peptides in membrane; to characterize the structures of low-energy states and their orientations with respect to the bilayer. Simulations were performed starting from different structures (unordered or helical) and orientations. It was found that the membrane environment significantly promotes an α-helical conformation for all the peptides, while their energetically favourable orientations are quite different. Thus, poly-Leu was immobilized inside the membrane, the hydrophobic peptide with polar termini adapted transbilayer orientation, whereas the amphiphilic peptide stayed on the lipid-water interface in peripherial orientation. Energy barriers between different states were characterized. The computational results were compared with the experimental structural data.     

A solvent model for simulations of peptides in bilayers. I. Membrane-promoting alpha-helix formation

We describe an efficient solvation model for proteins. In this model atomic solvation parameters imitating the hydrocarbon core of a membrane, water, and weak polar solvent (octanol) were developed. An optimal number of solvation parameters was chosen based on analysis of atomic hydrophobicities and fitting experimental free energies of gas-cyclohexane, gas-water, and octanol-water transfer for amino acids. The solvation energy term incorporated into the ECEPP/2 potential energy function was tested in Monte Carlo simulations of a number of small peptides with known energies of bilayer-water and octanol-water transfer. The calculated properties were shown to agree reasonably well with the experimental data. Furthermore, the solvation model was used to assess membrane-promoting alpha-helix formation. To accomplish this, all-atom models of 20-residue homopolypeptides-poly-Leu, poly-Val, poly-Ile, and poly-Gly in initial random coil conformation-were subjected to nonrestrained Monte Carlo conformational search in vacuo and with the solvation terms mimicking the water and hydrophobic parts of the bilayer. All the peptides demonstrated their largest helix-forming tendencies in a nonpolar environment, where the lowest-energy conformers of poly-Leu, Val, Ile revealed 100, 95, and 80% of alpha-helical content, respectively. Energetic and conformational properties of Gly in all environments were shown to be different from those observed for residues with hydrophobic side chains. Applications of the solvation model to simulations of peptides and proteins in the presence of membrane, along with limitations of the approach, are discussed.     

Site-Specific Peptide Probes Detect Buried Water in a Lipid Membrane

Transmembrane peptides contain polar residues in the interior of the membrane, which may alter the electrostatic environment and favor hydration in the otherwise nonpolar environment of the membrane core. Here, we demonstrate a general, nonperturbative strategy to probe hydration of the peptide backbone at specific depths within the bilayer using a combination of site-specific isotope labels, ultrafast two-dimensional infrared spectroscopy, and spectral modeling based on molecular dynamics simulations. Our results show that the amphiphilic pH-low insertion peptide supports a highly heterogeneous environment, with significant backbone hydration of nonpolar residues neighboring charged residues. For example, a leucine residue located as far as 1 nm into the hydrophobic bulk reports hydrogen-bonded populations as high as ∼20%. These findings indicate that the polar nature of these residues may facilitate the transport of water molecules into the hydrophobic core of the membrane.     

Evaluation of a fast implicit solvent model for molecular dynamics simulations.

A solvation term based on the solvent accessible surface area (SASA) is combined with the CHARMM polar hydrogen force field for the efficient simulation of peptides and small proteins in aqueous solution. Only two atomic solvation parameters are used: one is negative for favoring the direct solvation of polar groups and the other positive for taking into account the hydrophobic effect on apolar groups. To approximate the water screening effects on the intrasolute electrostatic interactions, a distance-dependent dielectric function is used and ionic side chains are neutralized. The use of an analytical approximation of the SASA renders the model extremely efficient (i.e., only about 50% slower than in vacuo simulations). The limitations and range of applicability of the SASA model are assessed by simulations of proteins and structured peptides. For the latter, the present study and results reported elsewhere show that with the SASA model it is possible to sample a significant amount of folding/unfolding transitions, which permit the study of the thermodynamics and kinetics of folding at an atomic level of detail.     

Explicit and implicit water simulations of a beta-hairpin peptide.

The conformational properties of a beta-hairpin peptide (YITNSDGTWT) were studied by using both explicit and implicit water simulations. The conformational space of the peptide was scanned by using a restricted hydrogen-bonding search method. The search method used generated the conformational space with enough diversity and good representation of beta-hairpin structures. By using a total surface area-based treatment of hydrophobic interactions, implicit water simulations failed to discriminate between experimental beta-hairpin structures from the rest of the conformers present in the authors' conformation library. However, with inclusion of vibrational free energy and accounting separately for polar and nonpolar surface areas, the nuclear magnetic resonance structure was ranked successfully as the most stable conformation. There is a loose correlation between the conformational energies by the continuum model and the conformational energies by explicit water simulation for conformers with similar structures. However, in terms of solvation energy, both approaches have a much better correlation. By using proper treatment of surface effect (partition of the surface area into polar and nonpolar areas) and including vibrational free-energy contribution, the continuum models should be reliable. Furthermore, the authors found that, for this peptide, beta-hairpin structures have large vibrational entropy that contributes decisively to the stability of folded beta-hairpin structures. Proteins 1999;37:73-87.     

Effects of total hydrophobicity and length of the hydrophobic domain of a signal peptide on in vitro translocation efficiency.

The hydrophobic domain of the signal peptide of OmpF-Lpp, a model secretory protein, was systematically engineered so as to be composed of different lengths of polyleucine residues or polymers with alternate leucine and alanine residues, and the effects of the length and nature of the hydrophobic stretch on the rate of in vitro translocation were studied using everted membrane vesicles of Escherichia coli. The translocation reaction exhibited high substrate specificity as to the number of hydrophobic residues. The results suggest that the hydrophobic domain is recognized specifically by a component(s) of the secretory machinery rather than nonspecifically by the hydrophobic region of the membrane. The in vitro translocation thus demonstrated required SecA and ATP and was markedly enhanced upon imposition of the proton motive force, as in the case of secretory proteins possessing a natural signal peptide. The highest translocation rate was obtained with the octamer in the case of polyleucine-containing signal peptides, whereas it was the decamer in the case of ones containing both leucine and alanine. These results suggest that the total hydrophobicity of the hydrophobic region of the signal peptides is an important determinant of the substrate specificity.     

The leader peptides from bacteriorhodopsin and halorhodopsin are potential membrane-spanning amphipathic helices

We show that the N-terminal leader peptides from the bacterial membrane proteins bacteriorhodopsin and halorhodopsin can be expected to form amphipathic alpha-helics with a highly hydrophobic nonpolar face and a narrow, negatively charged polar face. This finding is discussed in terms of a model for the integration of these proteins into the bacterial membrane.     

Control of the transmembrane orientation and interhelical interactions within membranes by hydrophobic helix length

We examined the effect of the length of the hydrophobic core of Lys-flanked poly(Leu) peptides on their behavior when inserted into model membranes. Peptide structure and membrane location were assessed by the fluorescence emission lambdamax of a Trp residue in the center of the peptide sequence, the quenching of Trp fluorescence by nitroxide-labeled lipids (parallax analysis), and circular dichroism. Peptides in which the hydrophobic core varied in length from 11 to 23 residues were found to be largely alpha-helical when inserted into the bilayer. In dioleoylphosphatidylcholine (diC18:1PC) bilayers, a peptide with a 19-residue hydrophobic core exhibited highly blue-shifted fluorescence, an indication of Trp location in a nonpolar environment, and quenching localized the Trp to the bilayer center, an indication of transmembrane structure. A peptide with an 11-residue hydrophobic core exhibited emission that was red-shifted, suggesting a more polar Trp environment, and quenching showed the Trp was significantly displaced from the bilayer center, indicating that this peptide formed a nontransmembranous structure. A peptide with a 23-residue hydrophobic core gave somewhat red-shifted fluorescence, but quenching demonstrated the Trp was still close to the bilayer center, consistent with a transmembrane structure. Analogous behavior was observed when the behavior of individual peptides was examined in model membranes with various bilayer widths. Other experiments demonstrated that in diC18:1PC bilayers the dilution of the membrane concentration of the peptide with a 23-residue hydrophobic core resulted in a blue shift of fluorescence, suggesting the red-shifted fluorescence at higher peptide concentrations was due to helix oligomerization. The intermolecular self-quenching of rhodamine observed when the peptide was rhodamine-labeled, and the concentration dependence of self-quenching, supported this conclusion. These studies indicate that the mismatch between helix length and bilayer width can control membrane location, orientation, and helix-helix interactions, and thus may mismatch control both membrane protein folding and the interactions between membrane proteins.     

Positions of polar amino acids alter interactions between transmembrane segments and detergents

α-Helical transmembrane (TM) segments in membrane proteins are comprised primarily of hydrophobic amino acids that accommodate insertion from water into the nonpolar membrane bilayer. In many such segments, however, polar residues are also present for structural or functional reasons. These latter residues impair the local favorable acyl interactions required for solvation by hydrophobic media such as phospholipids in native bilayers or detergents used for in vitro characterization. Using a series of Lys-tagged designed TM-like peptides (typified by KK-YAAAIAAIAWAIAAIAAAIAA-KKK) in which single-Asn residue substitutions (from Ile or Ala) were made successively from the center of the hydrophobic region toward the C-terminus, we demonstrate that polar residues strongly alter the nature of the interaction between TM segments and the solvating detergent. Through the application of sodium dodecyl sulfate-polyacrylamide gel electrophoresis, circular dichroism spectroscopy, and tryptophan fluorescence, we observed drastic differences in the structures of the detergent-peptide complexes that contain relatively minor sequence differences. For example, the blue shift of the Trp fluorescence (indicating local detergent solvation at this location) differs by as much as ~10 nm depending upon the position of a single Asn substitution in an otherwise identical segment. The overall results suggest that polar point mutations occurring in a biological membrane will elicit comparable effects, placing a significant refolding burden on the local protein structure and potentially leading to disease states through altered protein--lipid interactions in membrane proteins.     

3D Hydrophobic Moment Vectors as a Tool to Characterize the Surface Polarity of Amphiphilic Peptides

The interaction of membranes with peptides and proteins is largely determined by their amphiphilic character. Hydrophobic moments of helical segments are commonly derived from their two-dimensional helical wheel projections, and the same is true for β-sheets. However, to the best of our knowledge, there exists no method to describe structures in three dimensions or molecules with irregular shape. Here, we define the hydrophobic moment of a molecule as a vector in three dimensions by evaluating the surface distribution of all hydrophilic and lipophilic regions over any given shape. The electrostatic potential on the molecular surface is calculated based on the atomic point charges. The resulting hydrophobic moment vector is specific for the instantaneous conformation, and it takes into account all structural characteristics of the molecule, e.g., partial unfolding, bending, and side-chain torsion angles. Extended all-atom molecular dynamics simulations are then used to calculate the equilibrium hydrophobic moments for two antimicrobial peptides, gramicidin S and PGLa, under different conditions. We show that their effective hydrophobic moment vectors reflect the distribution of polar and nonpolar patches on the molecular surface and the calculated electrostatic surface potential. A comparison of simulations in solution and in lipid membranes shows how the peptides undergo internal conformational rearrangement upon binding to the bilayer surface. A good correlation with solid-state NMR data indicates that the hydrophobic moment vector can be used to predict the membrane binding geometry of peptides. This method is available as a web application on http://www.ibg.kit.edu/HM/.     

Structural details of urea binding to barnase: a molecular dynamics analysis.

BACKGROUND: The molecular mechanism of urea-induced protein unfolding has not been established. It is generally thought that denaturation results from the stabilizing interactions of urea with portions of the protein that are buried in the native state and become exposed upon unfolding of the protein. RESULTS: We have performed molecular dynamics simulations of barnase (a 110 amino acid RNase from Bacillus amyloliquefaciens) with explicit water and urea molecules at 300 K and 360 K. The native conformation was unaffected in the 300 K simulations at neutral and low pH. Two of the three runs at 360 K and low pH showed some denaturation, with partial unfolding of the hydrophobic core 2. The first solvation shell has a much higher density of urea molecules (water/urea ratio ranging from 2.07 to 2.73) than the bulk (water/urea ratio of 4.56). About one half of the first-shell urea molecules are involved in hydrogen bonds with polar or charged groups on the barnase surface, and between 15% and 18% of the first-shell urea molecules participate in multiple hydrogen bonds with barnase. The more stably bound urea molecules tend to be in crevices or pockets on the barnase surface. CONCLUSIONS: The simulation results indicate that an aqueous urea solution solvates the surface of a polypeptide chain more favorably than pure water. Urea molecules interact more favorably with nonpolar groups of the protein than water does, and the presence of urea improves the interactions of water molecules with the hydrophilic groups of the protein. The results suggest that urea denaturation involves effects on both nonpolar and polar groups of proteins.     

Membrane adsorption, folding, insertion and translocation of synthetic trans-membrane peptides

Spontaneous membrane adsorption, folding and insertion of the synthetic WALP16 and KALP16 peptides was studied by computer simulations starting from completely extended conformations. The peptides were simulated using an unmodified all-atom force field in combination with an efficient Monte Carlo sampling algorithm. The membrane is represented implicitly as a hydrophobic zone inside a continuum solvent modelled using the generalized Born theory of solvation. The method was previously parameterized to match insertion energies of hydrophobic side chain analogs into cyclohexane and no parameters were optimized for the present simulations. Both peptides rapidly precipitate out of bulk solution and adsorb to the membrane surface. Interfacial folding into a helical conformation is followed by membrane insertion. Both the peptide conformations and their location in the membrane are strongly temperature dependent. The temperature dependent behaviour can be summarized by fitting to a four-state model, separating the system into folded and unfolded conformers, which are either inserted into the membrane or located at the interfaces. As the temperature is lowered the dominant peptide conformation of the system changes from unfolded surface bound configurations to folded surface bound states. Folded trans-membrane conformers represent the dominant configuration at low temperatures. The analysis allows direct estimates of the free energies of peptide folding and membrane insertion. In the case of WALP the quality of the fit is excellent and the thermodynamic behaviour is in good agreement with expected theoretical consideration. For KALP the fit is more problematic due to the large solvation energies of the charged lysine residues.     

Membrane adsorption,folding, insertion and translocation of synthetic trans-membrane peptides

Spontaneous membrane adsorption, folding and insertion of the synthetic WALP16 and KALP16 peptides was studied by computer simulations starting from completely extended conformations. The peptides were simulated using an unmodified all-atom force field in combination with an efficient Monte Carlo sampling algorithm. The membrane is represented implicitly as a hydrophobic zone inside a continuum solvent modelled using the generalized Born theory of solvation. The method was previously parameterized to match insertion energies of hydrophobic side chain analogs into cyclohexane and no parameters were optimized for the present simulations. Both peptides rapidly precipitate out of bulk solution and adsorb to the membrane surface. Interfacial folding into a helical conformation is followed by membrane insertion. Both the peptide conformations and their location in the membrane are strongly temperature dependent. The temperature dependent behaviour can be summarized by fitting to a four-state model, separating the system into folded and unfolded conformers, which are either inserted into the membrane or located at the interfaces. As the temperature is lowered the dominant peptide conformation of the system changes from unfolded surface bound configurations to folded surface bound states. Folded trans-membrane conformers represent the dominant configuration at low temperatures. The analysis allows direct estimates of the free energies of peptide folding and membrane insertion. In the case of WALP the quality of the fit is excellent and the thermodynamic behaviour is in good agreement with expected theoretical consideration. For KALP the fit is more problematic due to the large solvation energies of the charged lysine residues.     

Hydrophobic organization of α-helix membrane bundle in bacteriorhodopsin

The hydrophobic organization of the intramembraneα-helical bundle in bacteriorhodopsin (BRh) was assessed based on a new approach to characterization of spatial hydrophobic properties of transmembrane (TM)α-helical peptides. The method employs two independent techniques: Monte Carlo simulations of nonpolar solvent around TM peptides and analysis of molecular hydrophobicity potential on their surfaces. The results obtained by the two methods agree with each other and permit precise hydrophobicity mapping of TM peptides. Superimposition of such data on the experimentally derived spatial model of the membrane moiety together with 2D maps of hydrophobic hydrophilic contacts provide considerable insight into the hydrophobic organization of BRh. The helix bundle is stabilized to a large extent by hydrophobic interactions between helices—neighbors in the sequence of BRh, by long-range interactions in helix pairs C-E, C-F, and C-G, and by nonpolar contracts between retinal and helices C, D, E, F. Unlike globular proteins, no polar contacts between residues distantly separated in the sequence of BRh were found in the bundle. One of the most striking results of this study is the finding that the hydrophobic organization of BRh is significantly different from those in bacterial photoreaction centers. Thus, TMα-helices in BRh expose their most nonpolar sides to the bilayer as well as to the neighboring helices and to the interior of the bundle. Some of them contact lipids with their relatively hydrophilic surfaces. No correlation was found between disposition of the most hydrophobic and the most variable sides of the TM helices.     

Conformational changes of peptides at solid/liquid interfaces: a Monte Carlo study

Monte Carlo simulations were performed to study the conformational changes of negatively charged model peptides dissolved in water adsorbed onto charged surfaces. 8-, 16-, and 20-residues peptides were used, each of them consisted of repeating diblock units of aspartic acid (ASP, polar amino acid) and isoleucine (ILE, nonpolar amino acid) residues. We found that a water patch was retained at the charged surface, separating the peptide from it. We believed that these water molecules were primarily responsible for giving a particular orientation to the peptide at the surface. Water did play a role to some extent in the structural stability of the 8-residues peptide. However, for higher chain lengths (16-residues and 20-residues), the intrinsic hydrogen-bonding network (or intrinsic structural stability) showed a predominant effect over hydrophobic dehydration for the stability of the peptide at the surface.     

A solvent model for simulations of peptides in bilayers. II. Membrane-spanning alpha-helices

We describe application of the implicit solvation model (see the first paper of this series), to Monte Carlo simulations of several peptides in bilayer- and water-mimetic environments, and in vacuum. The membrane-bound peptides chosen were transmembrane segments A and B of bacteriorhodopsin, the hydrophobic segment of surfactant lipoprotein, and magainin2. Their conformations in membrane-like media are known from the experiments. Also, molecular dynamics study of surfactant lipoprotein with different explicit solvents has been reported (Kovacs, H., A. E. Mark, J. Johansson, and W. F. van Gunsteren. 1995. J. Mol. Biol. 247:808-822). The principal goal of this work is to compare the results obtained in the framework of our solvation model with available experimental and computational data. The findings could be summarized as follows: 1) structural and energetic properties of studied molecules strongly depend on the solvent; membrane-mimetic media significantly promote formation of alpha-helices capable of traversing the bilayer, whereas a polar environment destabilizes alpha-helical conformation via reduction of solvent-exposed surface area and packing; 2) the structures calculated in a membrane-like environment agree with the experimental ones; 3) noticeable differences in conformation of surfactant lipoprotein assessed via Monte Carlo simulation with implicit solvent (this work) and molecular dynamics in explicit solvent were observed; 4) in vacuo simulations do not correctly reproduce protein-membrane interactions, and hence should be avoided in modeling membrane proteins.     

Protein–lipid interactions studied with designed transmembrane peptides: role of hydrophobic matching and interfacial anchoring (Review)

Biological membranes are characterized by a heterogeneous composition, which is not only manifested in the wide variety of their components, but also in aspects like the lateral organization, topology, and conformation of proteins and lipids. In bringing about the correct membrane structure, protein–lipid interactions can be expected to play a prominent role. The extent of hydrophobic matching between transmembrane protein segments and lipids potentially constitutes a versatile director of membrane organization, because a tendency to avoid hydrophobic mismatch could result in compensating adaptations such as tilt of the transmembrane segment or segregation into distinct domains. Also, interfacial interactions between lipid headgroups and the aromatic and charged residues that typically flank transmembrane domains may act as an organizing element. In this review, we discuss the numerous model studies that have systematically explored the influence of hydrophobic matching and interfacial anchoring on membrane structure. Designed peptides consisting of a polyleucine or polyleucine/alanine hydrophobic stretch, which is flanked on both sides by tryptophan or lysine residues, reflect the general layout of transmembrane protein segments. It is shown for phosphatidylcholine bilayers and for other model membranes that these peptides adapt a transmembrane topology without extensive peptide or lipid adaptations under conditions of hydrophobic matching, but that significant rearrangements can result from hydrophobic mismatch. Moreover, these effects depend on the nature of the flanking residues, implying a modulation of the mismatch response by interfacial interactions of the flanking residues. The implications of these model studies for the organization of biomembranes are discussed in the context of recent experiments with more complex systems.     

Protein-lipid interactions studied with designed transmembrane peptides: role of hydrophobic matching and interfacial anchoring

Biological membranes are characterized by a heterogeneous composition, which is not only manifested in the wide variety of their components, but also in aspects like the lateral organization, topology, and conformation of proteins and lipids. In bringing about the correct membrane structure, protein-lipid interactions can be expected to play a prominent role. The extent of hydrophobic matching between transmembrane protein segments and lipids potentially constitutes a versatile director of membrane organization, because a tendency to avoid hydrophobic mismatch could result in compensating adaptations such as tilt of the transmembrane segment or segregation into distinct domains. Also, interfacial interactions between lipid headgroups and the aromatic and charged residues that typically flank transmembrane domains may act as an organizing element. In this review, we discuss the numerous model studies that have systematically explored the influence of hydrophobic matching and interfacial anchoring on membrane structure. Designed peptides consisting of a polyleucine or polyleucine/alanine hydrophobic stretch, which is flanked on both sides by tryptophan or lysine residues, reflect the general layout of transmembrane protein segments. It is shown for phosphatidylcholine bilayers and for other model membranes that these peptides adapt a transmembrane topology without extensive peptide or lipid adaptations under conditions of hydrophobic matching, but that significant rearrangements can result from hydrophobic mismatch. Moreover, these effects depend on the nature of the flanking residues, implying a modulation of the mismatch response by interfacial interactions of the flanking residues. The implications of these model studies for the organization of biomembranes are discussed in the context of recent experiments with more complex systems.     

Amino-acid solvation structure in transmembrane helices from molecular dynamics simulations

Understanding the solvation of amino acids in biomembranes is an important step to better explain membrane protein folding. Several experimental studies have shown that polar residues are both common and important in transmembrane segments, which means they have to be solvated in the hydrophobic membrane, at least until helices have aggregated to form integral proteins. In this work, we have used computer simulations to unravel these interactions on the atomic level, and classify intramembrane solvation properties of amino acids. Simulations have been performed for systematic mutations in poly-Leu helices, including not only each amino acid type, but also every z-position in a model helix. Interestingly, many polar or charged residues do not desolvate completely, but rather retain hydration by snorkeling or pulling in water/headgroups--even to the extent where many of them exist in a microscopic polar environment, with hydration levels corresponding well to experimental hydrophobicity scales. This suggests that even for polar/charged residues a large part of solvation cost is due to entropy, not enthalpy loss. Both hydration level and hydrogen bonding exhibit clear position-dependence. Basic side chains cause much less membrane distortion than acidic, since they are able to form hydrogen bonds with carbonyl groups instead of water or headgroups. This preference is supported by sequence statistics, where basic residues have increased relative occurrence at carbonyl z-coordinates. Snorkeling effects and N-/C-terminal orientation bias are directly observed, which significantly reduces the effective thickness of the hydrophobic core. Aromatic side chains intercalate efficiently with lipid chains (improving Trp/Tyr anchoring to the interface) and Ser/Thr residues are stabilized by hydroxyl groups sharing hydrogen bonds to backbone oxygens.