Molecular modeling of the amphipathic helices of the plasma apolipoproteins.

In this paper we propose a classification of the amphipathic helical repeats occurring in the plasma apolipoprotein sequences. It is based upon the calculation of the molecular hydrophobicity potential around the helical segments. The repeats were identified using a new autocorrelation matrix, based upon similarities of hydrophobic and hydrophilic properties of the amino acid residues within the apolipoprotein sequences. The helices were constructed by molecular modeling, the molecular hydrophobicity potential was calculated, and isopotential contour lines drawn around the helices yielded a three-dimensional visualization of the hydrophobicity potential. Two classes of apolipoproteins could be differentiated by comparing the hydrophobic angles obtained by projection of the isopotential contour lines on a plane perpendicular to the long axis of the helix. The isopotential contour lines around apo AI, AIV, and E are more hydrophilic than hydrophobic, whereas they are of similar intensity for apo AII, CI, and CIII. In both cases discoidal lipid-protein complexes are generated, with the amphipathic helices around the edge of the lipid core. The long axis of the helices is oriented parallel to the phospholipid acyl chains and the hydrophilic side of the helix toward the aqueous phase. As a result of the differences in hydrophobicity potential, the contact between the hydrophobic side of the helices and the phospholipid acyl chains is larger for apo AII, CI, and CIII than for the other apolipoproteins. This might account for the greater stability of the discoidal complexes generated between phospholipids and these apoproteins.     

Conformational analysis of lipid-associating proteins in a lipid environment

Two major types of helical structures have been identified in lipid-associating proteins, being either amphipathic or transmembrane domains. A conformational analysis was carried out to characterize some of the properties of these helices. These calculations were performed both on isolated helices and in a lipid environment. According to the results of this analysis, the orientation of the line joining the hydrophobic and hydrophilic centers of the helix seems to determine the orientation of the helix at the lipid/water interface. The calculation of this parameter should be useful to discriminate between an amphipathic helix, parallel to the interface and a transmembrane helix orientated perpendicularly. The membrane-spanning helices are completely immersed in the phospholipid bilayer and their length corresponds to about the thickness of the hydrophobic core of the DPPC bilayer. The energy of interaction, expressed per phospholipid is significantly higher for the transmembrane compared to the amphipathic helices. For the membrane-spanning helices the mean energy of interaction is higher than the interaction energy between two phospholipids, while it is lower for most amphipathic helices. This might account for the stability of these protein-anchoring domains. This computer modeling approach should usefully complement the statistical analysis carried out on these helices, based on their hydrophobicity and hydrophobic moment. It represents a more refined analysis of the domains identified by the prediction techniques and stress the functional character of lipid-associating domains in membrane proteins as well as in soluble plasma lipoproteins.     

Calculation of the three-dimensional structure of Saccharomyces cerevisiae cytochrome b inserted in a lipid matrix

Cytochrome b is an integral membrane protein, which forms the core of the ubiquinol-cytochrome c oxidoreductase (cytochrome bc1) complex. A computer-aided three-dimensional modeling procedure was carried out in four steps. First, the candidate hydrophobic helices were searched for throughout the protein primary sequence by a computer procedure based upon the method of Eisenberg; second, a secondary helical structure was imposed to the transmembrane peptides; third, the helical segments at a lipid-water interface were oriented, and finally the possible interactions between helices with similar properties were investigated. This procedure enabled the identification of nine hydrophobic segments, of which eight are membrane-spanning helices while one has amphipathic properties. Three hydrophilic receptor-binding domains were also identified. Based upon their hydrophobicity profiles, the transmembrane helices could be associated in pairs inside the lipid bilayer. In our folding model proposed for cytochrome b, all mutation sites are not only located on the same side of the membrane but are also in close proximity in the three-dimensional structure. Inhibitor resistance mutational sites which were recently characterized (di Rago, J.-P., and Colson, A.-M. (1988) J. Biol. Chem. 263, 12564-12570) have been located on this model. Moreover, the receptor-binding domains and the mutation sites are close neighbors in the three-dimensional spatial representation.     

Mode of assembly of amphipathic helical segments in model high-density lipoproteins

The structure of discoidal apo A-I-phospholipid complexes, representing the metabolic precursors of mature high-density lipoprotein particles, was studied by a combination of both a theoretical and an experimental approach. The secondary structure of the complex was determined by circular dichroic measurements, while the relative orientation of the apo A-I helical segments and of the phospholipid acyl chains was determined by ATR infrared measurements. Fluorescence energy transfer between the tryptophan residues of apo A-I and fluorescent phospholipid probes yielded an estimation of the relative topography of the lipid and apolipoprotein components in discoidal and spherical particles. The theoretical approach consisted of the identification of the helical segments in various apo A-I species. These segments were then oriented at a lipid/water interface by minimization of their hydrophobic and hydrophilic transfer energies. The calculation of the hydrophobicity profiles along the axis of the helices leads to the identification of specific interactions between pairs of helices. The helices were further assembled together with the phospholipids by computer modelling, enabling an estimation of the dimensions of the complex. The combination of the experimental and theoretical results yielded a model for discoidal apolipoprotein-phospholipid complexes, in which the amphipathic helical segments are oriented along the edges of the discs. Such a model can be extended to the conversion of these complexes into mature spherical HDL, through the formation of a cholesteryl ester core.     

Local hydrophobicity stabilizes secondary structures in proteins

The probability of occurrence of helix and β-sheet residues in 47 globular proteins was determined as a function of local hydrophobicity, which was defined by the sum of the Nozaki-Tanford transfer free energies at two nearest-neighbors on both sides of the amino acid sequence. In general, hydrophilic amino acids favor neither helix nor β-sheet formations when neighbor residues are also hydrophilic but favor helix formation at higher local hydrophobicity. On the other hand, some hydrophobic amino acids such as Met, Leu, and Ile favor helix formation when neighbor residues are hydrophilic. None of the hydrophobic amino acids favor β-sheet formation with hydrophilic neighbors, but most of them strongly favor β-sheet formation at high local hydrophobicity. When the average of 20 amino acids is taken, both helix and β-sheet residue probabilities are higher at higher local hydrophobicity, although the increase is steeper for β-sheets. Therefore, β-sheet formation is more influenced by local hydrophobicity than helix formation. Generally, helices are nearer the surface and tend to have hydrophilic and hydrophobic faces at opposite sides. The tendency of alternating regions of hydrophilic and hydrophobic residues in a helical sequence was revealed by calculating the correlation of the Nozaki-Tanford values. Such amphipathic helices may be important in protein–protein and protein–lipid interactions and in forming hydrophilic channels in the membrane. The choice of 30 nonhomologous proteins as the data set did not alter the above results.     

Application of three-dimensional molecular hydrophobicity potential to the analysis of spatial organization of membrane domains in proteins: I. Hydrophobic properties of transmembrane segments of Na+, K+-ATPase

A new computer-aided molecular modeling approach based on the concept of three-dimensional (3D) molecular hydrophobicity potential has been developed to calculate the spatial organization of intramembrane domains in proteins. The method has been tested by calculating the arrangement of membrane-spanning segments in the photoreaction center ofRhodopseudomonas viridis and comparing the results obtained with those derived from the X-ray data. We have applied this computational procedure to the analysis of interhelical packing in membrane moiety of Na⁺, K⁺-ATPase. The work consists of three parts. In Part I, 3D distributions of electrostatic and molecular hydrophobicity potentials on the surfaces of transmembrane helical peptides were computed and visualized. The hydrophobic and electrostatic properties of helices are discussed from the point of view of their possible arrangement within the protein molecule. Interlocation of helical segments connected with short extramembrane loops found by means of optimization of their hydrophobic/hydrophilic contacts is considered in Part II. The most probable 3D model of packing of helical peptides in the membrane domain of Na⁺, K⁺-ATPase is discussed in the final part of the work.     

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.     

Application of three-dimensional molecular hydrophobicity potential to the analysis of spatial organization of membrane domains in proteins: I. Hydrophobic properties of transmembrane segments of Na+, K+-ATPase

A new computer-aided molecular modeling approach based on the concept of three-dimensional (3D) molecular hydrophobicity potential has been developed to calculate the spatial organization of intramembrane domains in proteins. The method has been tested by calculating the arrangement of membrane-spanning segments in the photoreaction center ofRhodopseudomonas viridis and comparing the results obtained with those derived from the X-ray data. We have applied this computational procedure to the analysis of interhelical packing in membrane moiety of Na⁺, K⁺-ATPase. The work consists of three parts. In Part I, 3D distributions of electrostatic and molecular hydrophobicity potentials on the surfaces of transmembrane helical peptides were computed and visualized. The hydrophobic and electrostatic properties of helices are discussed from the point of view of their possible arrangement within the protein molecule. Interlocation of helical segments connected with short extramembrane loops found by means of optimization of their hydrophobic/hydrophilic contacts is considered in Part II. The most probable 3D model of packing of helical peptides in the membrane domain of Na⁺, K⁺-ATPase is discussed in the final part of the work.     

Understanding peptide interactions with the lipid bilayer: a guide to membrane protein engineering

In recent years, studies of the interactions of designed hydrophobic and amphipathic polypeptides with biological membranes have progressed considerably. In-plane and transmembrane helical domains have been engineered as well as sequences that exhibit dynamic distributions of different topologies. These sequences not only help our understanding of the thermodynamic interaction contributions that determine peptide orientation, but also exhibit interesting biological functions as autonomous units. In addition, helices are considered to be folding intermediates of multi-spanning membrane proteins. The specificity and thermodynamic stability of transmembrane helix-helix interactions or the assembly of beta sheets at the membrane surface have, therefore, become a focus of ongoing research.     

Lipid-Protein Interactions Are a Unique Property and Defining Feature of G Protein-Coupled Receptors

G protein-coupled receptors (GPCRs) are membrane-bound proteins that depend on their lipid environment to carry out their physiological function. Combined efforts from many theoretical and experimental studies on the lipid-protein interaction profile of several GPCRs hint at an intricate relationship of these receptors with their surrounding membrane environment, with several lipids emerging as particularly important. Using coarse-grained molecular dynamics simulations, we explore the lipid-protein interaction profiles of 28 different GPCRs, spanning different levels of classification and conformational states and totaling to 1 ms of simulation time. We find a close relationship with lipids for all GPCRs simulated, in particular, cholesterol and phosphatidylinositol phosphate (PIP) lipids, but the number, location, and estimated strength of these interactions is dependent on the specific GPCR as well as its conformational state. Although both cholesterol and PIP lipids bind specifically to GPCRs, they utilize distinct mechanisms. Interactions with PIP lipids are mediated by charge-charge interactions with intracellular loop residues and stabilized by one or both of the transmembrane helices linked by the loop. Interactions with cholesterol, on the other hand, are mediated by a hydrophobic environment, usually made up of residues from more than one helix, capable of accommodating its ring structure and stabilized by interactions with aromatic and charged/polar residues. Cholesterol binding to GPCRs occurs in a small number of sites, some of which (like the binding site on the extracellular side of transmembrane 6/7) are shared among many class A GPCRs. Combined with a thorough investigation of the local membrane structure, our results provide a detailed picture of GPCR-lipid interactions. Additionally, we provide an accompanying website to interactively explore the lipid-protein interaction profile of all GPCRs simulated to facilitate analysis and comparison of our data.     

Membrane interactions of synthetic peptides corresponding to amphipathic helical segments of the human immunodeficiency virus type-1 envelope glycoprotein.

The human and simian immunodeficiency virus envelope glycoproteins, which mediate virus-induced cell fusion, contain two putative amphipathic helical segments with large helical hydrophobic moments near their carboxyl-terminal ends. In an attempt to elucidate the biological role of these amphipathic helical segments, we have synthesized peptides corresponding to residues 768-788 and 826-854 of HIV-1/WMJ-22 gp160. Circular dichroism studies of the peptides showed that the alpha helicity of the peptides increased with the addition of dimyristoyl phosphatidylcholine (DMPC) indicating that the peptides form lipid-associating amphipathic helixes. The peptides solubilized turbid suspensions of DMPC vesicles, and electron microscopy of peptide-DMPC mixtures revealed the formation of discoidal complexes, suggesting that the peptides bind to and perturb lipid bilayers. The peptides were found to lyse lipid vesicles and caused carboxyfluorescein leakage from dye-entrapped egg phosphatidylcholine liposomes. The peptides also lysed human erythrocytes and were found to be toxic to cell cultures. At subtoxic concentrations, the peptides effectively inhibited the fusion of CD4+ cells infected with recombinant vaccinia virus expressing human immunodeficiency virus (HIV)-1 envelope proteins. Based on these results, and reported studies on the mutational analysis of HIV envelope proteins, we suggest that the amphipathic helical segments near the carboxyl terminus of HIV envelope proteins may play a role in lysis of HIV-infected cells and also may modulate the extent of cell fusion observed during HIV infection of CD4+ cells.     

Interaction of 18-residue peptides derived from amphipathic helical segments of globular proteins with model membranes

We investigated the interaction of six 18-residue peptides derived from amphipathic helical segments of globular proteins with model membranes. The net charge of the peptides at neutral pH varies from −1 to +6. Circular dichroism spectra indicate that peptides with a high net positive charge tend to fold into a helical conformation in the presence of negatively charged lipid vesicles. In helical conformation, their average hydrophobic moment and hydrophobicity would render them surface-active. The composition of amino acids on the polar face of the helix in the peptides is considerably different. The peptides show variations in their ability to permeabilise zwitterionic and anionic lipid vesicles. Whereas increased net positive charge favours greater permeabilisation, the distribution of charged residues in the polar face also plays a role in determining membrane activity. The distribution of amino acids in the polar face of the helix in the peptides that were investigated do not fall into the canonical classes described. Amphipathic helices, which are part of proteins, with a pattern of amino acid distribution different from those observed in class L, A and others, could help in providing newer insights into peptide-membrane interactions.     

Study on the packing geometry, stoichiometry, and membrane interaction of three analogs related to a pore-forming small globular protein

A de novo designed pore-forming small globular protein (SGP) with antitumor activity consists of four helices: 3 basic amphipathic helices composed of Leu and Lys surrounding a central hydrophobic helix composed of oligoalanine. These helices are connected by a beta-turn-forming sequence and two beta-turn-unfavorable ones (S. Lee, T. Kiyota, T. Kunitake, E. Matsumoto, S. Yamashita, K. Anzai, and G. Sugihara Biochemistry 1997, Vol. 36, pp. 3782-3791). In the present work, we designed and synthesized three new SGP analogs in order to study the stoichiometric packing geometry and stability of SGP. The replacement of alanines in the central helix of SGP with leucines (SGP-L), which make the helix much larger in size and more hydrophobic, resulted in an equilibrium of monomeric-trimeric structure. The replacement of some Lys residues by Glu residues in the hydrophilic regions of the amphipathic helices (SGP-E) led to a decrease in helical content and the formation of an equilibrium of monomeric-trimeric structure. The alteration of beta-turn regions with Gly residues, which makes these regions flexible (SGP-G), established an equilibrium of monomeric-dimeric states in buffer. The hydrophobic alpha-helix of SGP-L penetrated into the lipid bilayers in a manner that stabilized model membranes and biomembranes, whereas the central helices of SGP-G and -E destabilized them by forming channels. SGP and its analogs may be a useful model to study the role of the hydrophobic and hydrophilic regions in the formation of monomer-oligomer of proteins and to better understand the insertion of membrane targeting proteins into biomembranes.     

Orientational preferences of neighboring helices can drive ER insertion of a marginally hydrophobic transmembrane helix

α-helical integral membrane proteins critically depend on the correct insertion of their transmembrane α helices into the lipid bilayer for proper folding, yet a surprisingly large fraction of the transmembrane α helices in multispanning integral membrane proteins are not sufficiently hydrophobic to insert into the target membrane by themselves. How can such marginally hydrophobic segments nevertheless form transmembrane helices in the folded structure? Here, we show that a transmembrane helix with a strong orientational preference (N(cyt)-C(lum) or N(lum)-C(cyt)) can both increase and decrease the hydrophobicity threshold for membrane insertion of a neighboring, marginally hydrophobic helix. This effect helps explain the "missing hydrophobicity" in polytopic membrane proteins.     

Stoichiometry of lipid interactions with transmembrane proteins--Deduced from the 3D structures

The stoichiometry of the first shell of lipids interacting with a transmembrane protein is defined operationally by the population of spin-labeled lipid chains whose motion is restricted directly by the protein. Interaction stoichiometries have been determined experimentally for a wide range of alpha-helical integral membrane proteins by using spin-label ESR spectroscopy. Here, we determine the spatially defined number of first-shell lipids at the hydrophobic perimeter of integral membrane proteins whose 3D structure has been determined by X-ray crystallography and lipid-protein interactions characterized by spin-labeling. Molecular modeling is used to build a single shell of lipids surrounding transmembrane structures derived from the PDB. Constrained energy optimization of the protein-lipid assemblies is performed by molecular mechanics. For relatively small proteins (up to 7-12 transmembrane helices), the geometrical first shell corresponds to that defined experimentally by perturbation of the lipid-chain dynamics. For larger, multi-subunit alpha-helical proteins, the lipids perturbed directly by the protein may either exceed or be less in number than those that can be accommodated at the intramembranous perimeter. In these latter cases, the motionally restricted spin-labeled lipids can be augmented by intercalation, or can correspond to a specific subpopulation at the protein interface, respectively. For monomeric beta-barrel proteins, the geometrical lipid stoichiometry corresponds to that determined from lipid mobility for a 22-stranded barrel, but fewer lipids are motionally restricted than can be accommodated around an eight-stranded barrel. Deviations from the geometrical first shell, in the beta-barrel case, are for the smaller protein with a highly curved barrel.     

Empirical lipid propensities of amino acid residues in multispan alpha helical membrane proteins

Characterizing the interactions between amino acid residues and lipid molecules is important for understanding the assembly of transmembrane helices and for studying membrane protein folding. In this study we develop TMLIP (TransMembrane helix-LIPid), an empirically derived propensity of individual residue types to face lipid membrane based on statistical analysis of high-resolution structures of membrane proteins. Lipid accessibilities of amino acid residues within the transmembrane (TM) region of 29 structures of helical membrane proteins are studied with a spherical probe of radius of 1.9 A. Our results show that there are characteristic preferences for residues to face the headgroup region and the hydrocarbon core region of lipid membrane. Amino acid residues Lys, Arg, Trp, Phe, and Leu are often found exposed at the headgroup regions of the membrane, where they have high propensity to face phospholipid headgroups and glycerol backbones. In the hydrocarbon core region, the strongest preference for interacting with lipids is observed for Ile, Leu, Phe and Val. Small and polar amino acid residues are usually buried inside helical bundles and are strongly lipophobic. There is a strong correlation between various hydrophobicity scales and the propensity of a given residue to face the lipids in the hydrocarbon region of the bilayer. Our data suggest a possibly significant contribution of the lipophobic effect to the folding of membrane proteins. This study shows that membrane proteins have exceedingly apolar exteriors rather than highly polar interiors. Prediction of lipid-facing surfaces of boundary helices using TMLIP1 results in a 54% accuracy, which is significantly better than random (25% accuracy). We also compare performance of TMLIP with another lipid propensity scale, kPROT, and with several hydrophobicity scales using hydrophobic moment analysis.     

Construction of protocellular structures under simulated primitive earth conditions

We have developed experimental approaches for the construction of protocellular structures under simulated primitive earth conditions and studied their formation and characteristics. Three types of envelopes; protein envelopes, lipid envelopes, and lipid-protein envelopes are considered as candidates for protocellular structures. Simple protein envelopes and lipid envelopes are presumed to have originated at an early stage of chemical evolution, interaction mutually and then evolved into more complex envelopes composed of both lipids and proteins.Three kinds of protein envelopes were constructedin situ from amino acids under simulated primitive earth conditions such as a fresh water tide pool, a warm sea, and a submarine hydrothermal vent. One protein envelope was formed from a mixture of amino acid amides at 80 °C using multiple hydration-dehydration cycles. Marigranules, protein envelope structures, were produced from mixtures of glycine and acidic, basic and aromatic amino acids at 105 °C in a modified sea medium enriched with essential transition elements. Thermostable microspheres were also formed from a mixture of glycine, alanine, valine, and aspartic acid at 250 °C and above. The microspheres did not form at lower temperatures and consist of silicates and peptide-like polymers containing imide bonds and amino acid residues enriched in valine. Amphiphilic proteins with molecular weights of 2000 were necessary for the formation of the protein envelopes.Stable lipid envelopes were formed from different dialkyl phospholipids and fatty acids.Large, stable, lipid-protein envelopes were formed from egg lecithin and the solubilized marigranules. Polycations such as polylysine and polyhistidine, or basic proteins such as lysozyme and cytochromec also stabilized lipid-protein envelopes.     

Initial structural and dynamic characterization of the M2 protein transmembrane and amphipathic helices in lipid bilayers

Amphipathic helices in membrane proteins that interact with the hydrophobic/hydrophilic interface of the lipid bilayer have been difficult to structurally characterize. Here, the backbone structure and orientation of an amphipathic helix in the full-length M2 protein from influenza A virus has been characterized. The protein has been studied in hydrated DMPC/DMPG lipid bilayers above the gel to liquid-crystalline phase transition temperature by solid-state NMR spectroscopy. Characteristic PISA (Polar Index Slant Angle) wheels reflecting helical wheels have been observed in uniformly aligned bilayer preparations of both uniformly 15N labeled and amino acid specific labeled M2 samples. Hydrogen/deuterium exchange studies have shown the very slow exchange of some residues in the amphipathic helix and more rapid exchange for the transmembrane helix. These latter results clearly suggest the presence of an aqueous pore. A variation in exchange rate about the transmembrane helical axis provides additional support for this claim and suggests that motions occur about the helical axes in this tetramer to expose the entire backbone to the pore.     

Computational studies of colicin insertion into membranes: the closed state

Colicins are water-soluble toxins that, upon interaction with membranes, undergo a conformational change, insert, and form pores in them. Pore formation activity is localized in a bundle of 10 α-helices named the pore-forming domain (PFD). There is evidence that colicins attach to the membrane via a hydrophobic hairpin embedded in the core of the PFD. Two main models have been suggested for the membrane-bound state: penknife and umbrella, differing in regard to the orientation of the hydrophobic hairpin with respect to the membrane. The arrangement of the amphipathic helices has been described as either a compact three-dimensional structure or a two-dimensional array of loosely interacting helices on the membrane surface. Using molecular dynamics simulations with an implicit membrane model, we studied the structure and stability of the conformations proposed earlier for four colicins. We find that colicins are initially driven towards the membrane by electrostatic interactions between basic residues and the negatively charged membrane surface. They do not have a unique binding orientation, but in the predominant orientations the central hydrophobic hairpin is parallel to the membrane. In the inserted state, the estimated free energy tends to be lower for the compact arrangements of the amphipathic helix, but the more expanded ones are in better agreement with experimental distance distributions. The difference in energy between penknife and umbrella conformations is small enough for equilibrium to exist between them. Elongation of the hydrophobic hairpin helices and membrane thinning were found unable to produce stabilization of the transmembrane configuration of the hydrophobic hairpin.