Characterizing Residue-Bilayer Interactions Using Gramicidin A as a Scaffold and Tryptophan Substitutions as Probes

This document is the Accepted Manuscript version of a Published Work that appeared in final form in Journal of Chemical Theory and Computation, copyright © American Chemical Society after peer review and technical editing by the publisher. To access the final edited and published work see http://doi.org/10.1021/acs.jctc.7b00400.Previous experiments have shown that the lifetime of a gramicidin A dimer channel (which forms from two non-conducting monomers) in a lipid bilayer is modulated by mutations of the tryptophan (Trp) residues at the bilayer-water interface. We explore this further using extensive molecular dynamics simulations of various gA dimer and monomer mutants at the Trp positions in phosphatidylcholine bilayers with different tail lengths. gA interactions with the surrounding bilayer are strongly modulated by mutating these Trp residues. There are three principal effects: eliminating residue hydrogen bonding ability (i.e., reducing the channel-monolayer coupling strength) reduces the extent of the bilayer deformation caused by the assembled dimeric channel; a residue’s size and geometry affects its orientation, leading to different hydrogen bonding partners; and increasing a residue’s hydrophobicity increases the depth of gA monomer insertion relative to the bilayer center, thereby increasing the lipid bending frustration

Lipid bilayer thickness determines cholesterol's location in model membranes

Cholesterol is an essential biomolecule of animal cell membranes, and an important precursor for the biosynthesis of certain hormones and vitamins. It is also thought to play a key role in cell signaling processes associated with functional plasma membrane microdomains (domains enriched in cholesterol), commonly referred to as rafts. In all of these diverse biological phenomena, the transverse location of cholesterol in the membrane is almost certainly an important structural feature. Using a combination of neutron scattering and solid-state 2H NMR, we have determined the location and orientation of cholesterol in phosphatidylcholine (PC) model membranes having fatty acids of different lengths and degrees of unsaturation. The data establish that cholesterol reorients rapidly about the bilayer normal in all the membranes studied, but is tilted and forced to span the bilayer midplane in the very thin bilayers. The possibility that cholesterol lies flat in the middle of bilayers, including those made from PC lipids containing polyunsaturated fatty acids (PUFAs), is ruled out. These results support the notion that hydrophobic thickness is the primary determinant of cholesterol's location in membranes

Comparing Interfacial Trp, Interfacial His and pH Dependence for the Anchoring of Tilted Transmembrane Helical Peptides

Charged and aromatic amino acid residues, being enriched toward the terminals of membrane-spanning helices in membrane proteins, help to stabilize particular transmembrane orientations. Among them, histidine is aromatic and can be positively charge at low pH. To enable investigations of the underlying protein-lipid interactions, we have examined the effects of single or pairs of interfacial histidine residues using the constructive low-dynamic GWALP23 (acetyl-GG2ALW5LALALALALALALW19LAG22A-amide) peptide framework by incorporating individual or paired histidines at locations 2, 5, 19 or 22. Analysis of helix orientation by means of solid-state 2H NMR spectra of labeled alanine residues reveals marked differences with H2,22 compared to W2,22. Nevertheless, the properties of membrane-spanning H2,22WALP23 helices show little pH dependence and are similar to those having Gly, Arg or Lys at positions 2 and 22. The presence of H5 or H19 influences the helix rotational preference but not the tilt magnitude. H5 affects the helical integrity, as residue 7 unwinds from the core helix; yet once again the helix orientation and dynamic properties show little sensitivity to pH. The overall results reveal that the detailed properties of transmembrane helices depend upon the precise locations of interfacial histidine residues

A combined experimental and theoretical study of ion solvation in liquid N-methylacetamide

Most current biomolecular simulations are based on potential energy functions that treat the electrostatic energy as a sum of pairwise Coulombic interactions between effective fixed atomic charges. This approximation, in which many-body induced polarization effects are included in an average way, is expected to be satisfactory for a wide range of systems, but less accurate for processes involving the transfer and partition of ions among heterogeneous environments. The limitations of these potential energy functions are perhaps most obvious in studies of ion permeation through membrane channels. In many cases, the pore is so narrow that the permeating ion must shed most of its surrounding water molecules and the large energetic loss due to dehydration must be compensated by coordination with protein atoms. Interactions of cations with protein backbone carbonyl oxygens, in particular, play a critical role in several important biological channels. As a first step toward meeting the challenge of developing an accurate explicit accounting for induced polarization effects, the present work combines experiments and computation to characterize the interactions of alkali and halide ions with N-methylacetamide chosen to represent the peptide bond. From solubility measurements, we extract the solvation free energies of KCl and NaCl in liquid N-methylacetamide. Polarizable models based on the Drude oscillator are then developed and compared with available experimental and ab initio data. The good agreement for a range of structural and thermodynamic properties in the gas and condensed phases suggests that the polarizable models provide an accurate representation of ion−amide interactions in biological systems

Sensitivity of single membrane-spanning alpha-helical peptides to hydrophobic mismatch with a lipid bilayer: Effects on backbone structure, orientation, and extent of membrane incorporation

The extent of matching of membrane hydrophobic thickness with the hydrophobic length of transmembrane protein segments potentially constitutes a major director of membrane organization. Therefore, the extent of mismatch that can be compensated, and the types of membrane rearrangements that result, can provide valuable insight into membrane functionality. In the present study, a large family of synthetic peptides and lipids is used to investigate a range of mismatch situations. Peptide conformation, orientation, and extent of incorporation are assessed by infrared spectroscopy, tryptophan fluorescence, circular dichroism, and sucrose gradient centrifugation. It is shown that peptide backbone structure is not significantly affected by mismatch, even when the extent of mismatch is large. Instead, this study demonstrates that for tryptophan-flanked peptides the dominant response of a membrane to large mismatch is that the extent of incorporation is reduced, when the peptide is both too short and too long. With increasing mismatch, a smaller fraction of peptide is incorporated into the lipid bilayer, and a larger fraction is present in extramembranous aggregates. Relatively long peptides that remain incorporated in the bilayer have a small tilt angle with respect to the membrane normal. The observed effects depend on the nature of the flanking residues: long tryptophan-flanked peptides do not associate well with thin bilayers, while equisized lysine-flanked peptides associate completely, thus supporting the notion that tryptophan and lysine interact differently with membrane-water interfaces. The different properties that aromatic and charged flanking residues impart on transmembrane protein segments are discussed in relation to protein incorporation in biological systems

The effect of peptide/lipid hydrophobic mismatch on the phase behavior of model membranes mimicking the lipid composition in Escherichia coli membranes.

The effect of hydrophobic peptides on the lipid phase behavior of an aqueous dispersion of dioleoylphosphatidylethanolamine and dioleoylphosphatidylglycerol (7:3 molar ratio) was studied by (31)P NMR spectroscopy. The peptides (WALPn peptides, where n is the total number of amino acid residues) are designed as models for transmembrane parts of integral membrane proteins and consist of a hydrophobic sequence of alternating leucines and alanines, of variable length, that is flanked on both ends by tryptophans. The pure lipid dispersion was shown to undergo a lamellar-to-isotropic phase transition at approximately 60 degrees C. Small-angle x-ray scattering showed that at a lower water content a cubic phase belonging to the space group Pn3m is formed, suggesting also that the isotropic phase in the lipid dispersion represents a cubic liquid crystalline phase. It was found that the WALP peptides very efficiently promote formation of nonlamellar phases in this lipid system. At a peptide-to-lipid (P/L) molar ratio of 1:1000, the shortest peptide used, WALP16, lowered the lamellar-to-isotropic phase transition by approximately 15 degrees C. This effect was less for longer peptides. For all of the WALP peptides used, an increase in peptide concentration led to a further lowering of the phase transition temperature. At the highest P/L ratio (1:25) studied, WALP16 induced a reversed hexagonal liquid crystalline (H(II)) phase, while the longer peptides still promoted the formation of an isotropic phase. Peptides with a hydrophobic length larger than the bilayer thickness were found to be unable to inhibit formation of the isotropic phase. The results are discussed in terms of mismatch between the hydrophobic length of the peptide and the hydrophobic thickness of the lipid bilayer and its consequences for lipid-protein interactions in membranes

Influence of Lipid/Peptide Hydrophobic Mismatch on the Thickness of Diacylphosphatidylcholine Bilayers. A 2H NMR and ESR Study Using Designed Transmembrane Alpha-Helical Peptides and Gramicidin A

We have investigated the effect of a series of hydrophobic polypeptides (WALP peptides) on the mean hydrophobic thickness of (chain-perdeuterated) phosphatidylcholines (PCs) with different acyl chain length, using 2H NMR and ESR techniques. The WALP peptides are uncharged and consist of a sequence with variable length of alternating leucine and alanine, flanked on both sides by two tryptophans, and with the N- and C-termini blocked, e.g., FmAW2(LA)nW2AEtn. 2H NMR measurements showed that the shortest peptide with a total length of 16 amino acids (WALP16) causes an increase of 0.6 Å in bilayer thickness in di-C12-PC, a smaller increase in di-C14-PC, no effect in di-C16-PC, and a decrease of 0.4 Å in di-C18-PC, which was the largest decrease observed in any of the peptide/lipid systems. The longest peptide, WALP19, in di-C12-PC caused the largest increase in thickness of the series (+1.4 Å), which decreased again for longer lipids toward di-C18-PC, in which no effect was noticed. WALP17 displayed an influence intermediate between that of WALP16 and WALP19. Altogether, incorporation of the WALP peptides was found to result in small but very systematic changes in bilayer thickness and area per lipid molecule, depending on the difference in hydrophobic length between the peptide and the lipid bilayer in the liquid-crystalline phase. ESR measurements with spin-labeled lipid probes confirmed this result. Because thickness is expected to be influenced most at the lipids directly adjacent to the peptides, also the maximal adaptation of these first-shell lipids was estimated. The calculation was based on the assumption that there is little or no aggregation of the WALP peptides, as was supported by ESR, and that lipid exchange is rapid on the 2H NMR time scale. It was found that even the maximal possible changes in first-shell lipid length were relatively small and represented only a partial response to mismatch. The synthetic WALP peptides are structurally related to the gramicidin channel, which was therefore used for comparison. In most lipid systems, gramicidin proved to be a stronger perturber of bilayer thickness than WALP19, although its length should approximate that of the shorter WALP16. The effects of gramicidin and WALP peptides on bilayer thickness were evaluated with respect to previous 31P NMR studies on the effects of these peptides on macroscopic lipid phase behavior. Both approaches indicate that, in addition to the effective hydrophobic length, also the physical nature of the peptide surface is a modulator of lipid order