Thermostability of membrane protein helix-helix interaction elucidated by statistical analysis

A prerequisite for the survival of (micro)organisms at high temperatures is an adaptation of protein stability to extreme environmental conditions. In contrast to soluble proteins, where many factors have already been identified, the mechanisms by which the thermostability of membrane proteins is enhanced are almost unknown. The hydrophobic membrane environment constrains possible stabilizing factors for transmembrane domains, so that a difference might be expected between soluble and membrane proteins. Here we present sequence analysis of predicted transmembrane helices of the genomes from eight thermophilic and 12 mesophilic organisms. A comparison of the amino acid compositions indicates that more polar residues can be found in the transmembrane helices of thermophilic organisms. Particularly, the amino acids aspartic acid and glutamic acid replace the corresponding amides. Cysteine residues are found to be significantly decreased by about 70% in thermophilic membrane domains suggesting a non-specific function of most cysteine residues in transmembrane domains of mesophilic organisms. By a pair-motif analysis of the two sets of transmembrane helices, we found that the small residues glycine and serine contribute more to transmembrane helix-helix interactions in thermophilic organisms. This may result in a tighter packing of the helices allowing more hydrogen bond formation.     

Extension of the fragment method to calculate amino acid zwitterion and side chain partition coefficients

The fragment method of calculating partition coefficients (P) has been extended to include the common amino acids (AAs). The results indicate that polar and charged side chains influence the hydrophobicity of atoms in the side chain in a predictable manner. Field effects, as evidenced through polar proximity factors and bond factors, need to be considered for accurate estimation of transfer phenomena. The calculated log P and delta G degree ' values of the 20 AAs agree well with the observed values. Pro calculates to be more hydrophilic than the observed log P. Hydrophobicity scales for peptide side chain residues are compared and evaluated in terms of suitability. Calculated pi values for nonpolar side chain residues agree well with the observed values; calculated values for uncharged polar side chain residues deviate by about 0.6 log units except for Gln and Cys; and polar side chain residues with charged side chains calculate as too hydrophilic. Reasons for the differences are explored. We also suggest that tightly bound water to polar moieties in amino acids and peptides may be transferred into the octanol phase during partitioning experiments. A quantitative methodology is presented which characterizes the thermodynamic partitioning of groups and individual atoms in amino acids and proteins.     

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.     

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.     

A new amphipathy scale I. Determination of the scale from molecular dynamics data

Two new amphipathy scales elaborated from molecular dynamics data are presented. Their applications contribute for the identification of the hydrophobic or hydrophilic regions in proteins solely from the primary structure. The new amphipathy coefficients (AC) reflect the side chain/solvent molecules configurational energies. A polar (water) and an apolar solvent, CCl4, were used resulting in the two ACwater and ACCCl4 scales. These solvents were chosen to simulate the aqueous phases and the transmembrane ambients of cellular membranes where the membrane proteins act. The new amphipathy scales were compared with some previous scales determined by different methods, which were also compared between them, indicating more than 90% of the correlation coefficients are less than 0.9: the scales are strictly dependent on the methodologies used in their determination. The ACCCl4 scale is related with the size of side chain amino acids while ACwater is related with the hydrophobicity of side chain amino acids. The quality of the scales was confirmed by an example of application where ACwater was able to identify correctly the transmembrane, hydrophobic regions of a membrane protein. These results also indicate that water is an important factor responsible for the tertiary structure of membrane proteins.     

Analysis and prediction of helix-helix interactions in membrane channels and transporters

Membrane proteins span a large variety of different functions such as cell-surface receptors, redox proteins, ion channels, and transporters. Proteins with functional pores show different characteristics of helix-helix packing as other helical membrane proteins. We found that the helix-helix contacts of 13 nonhomologous high-resolution structures of membrane channels and transporters are mainly accomplished by weakly polar amino acids (G > S > T > F) that preferably create contacts every fourth residue, typical for right-handed helix crossings. There is a strong correlation between the now available biological hydrophobicity scale and the propensities of the weakly polar and hydrophobic residues to be buried at helix-helix interfaces or to be exposed to the lipids in membrane channels and transporters. The polar residues, however, make no major contribution towards the packing of their transmembrane helices, and are therefore subsumed to be primarily exposed to the polar milieu during the folding process. The contact formation of membrane channels and transporters is therefore ruled by the solubility of the residues, which we suppose to be the driving force for the assembly of their transmembrane helices. By contrast, in 14 nonhomologous high-resolution structures of other membrane protein coils, also large and polar amino acids (D > S > M > Q) create characteristic contacts every 3.5th residues, which is a signature for left-handed helix crossings. Accordingly, it seems that dependent on the function, different concepts of folding and stabilization are realized for helical membrane proteins. Using a sequence-based matrix prediction method these differences are exploited to improve the prediction of buried and exposed residues of transmembrane helices significantly. When the sequence motifs typical for membrane channels and transporters were applied for the prediction of helix-helix contacts the quality of prediction rises by 16% to an average value of 76%, compared to the same approach when only single amino acid positions are taken into account.     

Helix kinks are equally prevalent in soluble and membrane proteins

Helix kinks are a common feature of α‐helical membrane proteins, but are thought to be rare in soluble proteins. In this study we find that kinks are a feature of long α‐helices in both soluble and membrane proteins, rather than just transmembrane α‐helices. The apparent rarity of kinks in soluble proteins is due to the relative infrequency of long helices (≥20 residues) in these proteins. We compare length‐matched sets of soluble and membrane helices, and find that the frequency of kinks, the role of Proline, the patterns of other amino acid around kinks (allowing for the expected differences in amino acid distributions between the two types of protein), and the effects of hydrogen bonds are the same for the two types of helices. In both types of protein, helices that contain Proline in the second and subsequent turns are very frequently kinked. However, there are a sizeable proportion of kinked helices that do not contain a Proline in either their sequence or sequence homolog. Moreover, we observe that in soluble proteins, kinked helices have a structural preference in that they typically point into the solvent. Proteins 2014; 82:1960–1970. © 2014 The Authors. Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.     

Missense mutations in transmembrane domains of proteins: phenotypic propensity of polar residues for human disease

Previous experiments on the cystic fibrosis transmembrane conductance regulator suggested that non-native polar residues within membrane domains can compromise protein structure/function. However, depending on context, replacement of a native residue by a non-native residue can result either in genetic disease or in benign effects (e.g., polymorphisms). Knowledge of missense mutations that frequently cause protein malfunction and subsequent disease can accordingly reveal information as to the impact of these residues in local protein environments. We exploited this concept by performing a statistical comparison of disease-causing mutations in protein membrane-spanning domains versus soluble domains. Using the Human Gene Mutation Database of 240 proteins (including 80 membrane proteins) associated with human disease, we compared the relative phenotypic propensity to cause disease of the 20 naturally occurring amino acids when removed from-or inserted into-native protein sequences. We found that in transmembrane domains (TMDs), mutations involving polar residues, and ionizable residues in particular (notably arginine), are more often associated with protein malfunction than soluble proteins. To further test the hypothesis that interhelical cross-links formed by membrane-embedded polar residues stabilize TMDs, we compared the occurrence of such residues in the TMDs of mesophilic and thermophilic prokaryotes. Results showed a significantly higher proportion of ionizable residues in thermophilic organisms, reinforcing the notion that membrane-embedded electrostatic interactions play critical roles in TMD stability.     

TIM23-mediated insertion of transmembrane α-helices into the mitochondrial inner membrane

While overall hydrophobicity is generally recognized as the main characteristic of transmembrane (TM) α-helices, the only membrane system for which there are detailed quantitative data on how different amino acids contribute to the overall efficiency of membrane insertion is the endoplasmic reticulum (ER) of eukaryotic cells. Here, we provide comparable data for TIM23-mediated membrane protein insertion into the inner mitochondrial membrane of yeast cells. We find that hydrophobicity and the location of polar and aromatic residues are strong determinants of membrane insertion. These results parallel what has been found previously for the ER. However, we see striking differences between the effects elicited by charged residues flanking the TM segments when comparing the mitochondrial inner membrane and the ER, pointing to an unanticipated difference between the two insertion systems.     

Important inter-residue contacts for enhancing the thermal stability of thermophilic proteins

Proteins from thermophilic organisms exhibit high thermal stability, but have structures that are very similar to their mesophilic homologues. In order to gain insight into the basis of thermostability, we have analyzed the medium- and long-range contacts in mesophilic and thermophilic proteins of 16 different families. We found that the thermophiles prefer to have contacts between residues with hydrogen-bond-forming capability. Apart from hydrophobic contacts, more contacts are observed between polar and non-polar residues in thermophiles than mesophiles. Residue-wise analysis showed that Tyr has good contacts with several other residues, and Cys has considerably higher long-range contacts in thermophiles compared with mesophiles. Furthermore, the residues occurring in the range of 31-34 residues apart in the sequence contribute significant long-range contacts to the stability of thermophilic proteins.     

Differences in amino acids composition and coupling patterns between mesophilic and thermophilic proteins

Summary. Thermophilic proteins show substantially higher intrinsic thermal stability than their mesophilic counterparts. Amino acid composition is believed to alter the intrinsic stability of proteins. Several investigations and mutagenesis experiment have been carried out to understand the amino acid composition for the thermostability of proteins. This review presents some generalized features of amino acid composition found in thermophilic proteins, including an increase in residue hydrophobicity, a decrease in uncharged polar residues, an increase in charged residues, an increase in aromatic residues, certain amino acid coupling patterns and amino acid preferences for thermophilic proteins. The differences of amino acids composition between thermophilic and mesophilic proteins are related to some properties of amino acids. These features provide guidelines for engineering mesophilic protein to thermophilic protein. Authors’ addresses: Yuan-Jiang Pan, Institute of Chemical Biology and Pharmaceutical Chemistry, Zhejiang University, Zhejiang University Road 38, Hangzhou 310027, China; Wei-Fen Li, Microbiology Division, College of Animal Science, Zhejiang University, Hangzhou 310029, China     

A new amphipathy scale: I. Determination of the scale from molecular dynamics data

Two new amphipathy scales elaborated from molecular dynamics data are presented. Their applications contribute for the identification of the hydrophobic or hydrophilic regions in proteins solely from the primary structure. The new amphipathy coefficients (AC) reflect the side chain/solvent molecules configurational energies. A polar (water) and an apolar solvent, CCl₄, were used resulting in the two AC_water and AC_CCl4 scales. These solvents were chosen to simulate the aqueous phases and the transmembrane ambients of cellular membranes where the membrane proteins act. The new amphipathy scales were compared with some previous scales determined by different methods, which were also compared between them, indicating more than 90% of the correlation coefficients are less than 0.9: the scales are strictly dependent on the methodologies used in their determination. The AC_CCl4 scale is related with the size of side chain amino acids while AC_water is related with the hydrophobicity of side chain amino acids. The quality of the scales was confirmed by an example of application where AC_water was able to identify correctly the transmembrane, hydrophobic regions of a membrane protein. These results also indicate that water is an important factor responsible for the tertiary structure of membrane proteins.     

Hydrophobicity regained.

A widespread practice is to use free energies of transfer between organic solvents and water (delta G0transfer to define hydrophobicity scales for the amino acid side chains. A comparison of four delta G0transfer scales reveals that the values for hydrogen-bonding side chains are highly dependent on the non-aqueous environment. This property of polar side chains violates the assumptions underlying the paradigm of equating delta G0transfer with hydrophobicity or even with a generic solvation energy that is directly relevant to protein stability and ligand binding energetics. This simple regaining of the original concept of hydrophobicity reveals a flaw in approaches that use delta G0transfer values to derive generic estimates of the energetics of the burial of polar groups, and allows the introduction of a "pure" hydrophobicity scale for the amino acid residues.     

Snorkeling preferences foster an amino acid composition bias in transmembrane helices

By analyzing transmembrane (TM) helices in known structures, we find that some polar amino acids are more frequent at the N terminus than at the C terminus. We propose the asymmetry occurs because most polar amino acids are better able to snorkel their polar atoms away from the membrane core at the N terminus than at the C terminus. Two findings lead us to this proposition: (1) side-chain conformations are influenced strongly by the N or C-terminal position of the amino acid in the bilayer, and (2) the favored snorkeling direction of an amino acid correlates well with its N to C-terminal composition bias. Our results suggest that TM helix predictions should incorporate an N to C-terminal composition bias, that rotamer preferences of TM side-chains are position-dependent, and that the ability to snorkel influences the evolutionary selection of amino acids for the helix N and C termini.     

Study of the molecular organization of visual rhodopsin in photoreceptor membranes by limited proteolysis

Proteolysis of rhodopsin in disc membranes of right-side out orientation by thermolysin, papain and St. aureus V8 protease allowed to identify two highly exposed regions of polypeptide chain located on the cytoplasmic membrane surface: carboxyl terminal sequence 321-348 and the fragment 236-241. Incubation with chymotrypsin reveals the third site on the cytoplasmic surface, 146-147, accessible to proteolytic enzymes. Frozen-thawed membranes comprise a mixture of vesicles with normal and inverted orientation. Both thermolytic and chymotryptic digests of rhodopsin in these membranes contain the polypeptide which represents the amino terminal sequence lacking the first 30 amino acid residues. Thus at least 30 amino acids from the N-terminus must protrude into the intradiscal space. One additional site was located on the intradiscal surface: papain digests rhodopsin in the inverted membranes at the position 186-187. Localization of the proteolytic cleavage sites allowed to propose a model for rhodopsin topography in disc membrane: the polypeptide chain traverses the bilayer thickness seven times; each of seven transmembrane segments containing approximately 40 amino acid residues includes a sequence of approximately 30 hydrophobic amino acids; which are probably in close contact with hydrocarbon matrix of the membrane. Hydrophobic sequences are terminated with fragments containing clusters of hydrophilic amino acids, possibly interacting with lipid polar head groups and orienting each segment in the bilayer.     

Patterns of temperature adaptation in proteins from Methanococcus and Bacillus

It has long been known that amino acid substitutions in proteins of organisms living at moderate and high temperatures (mesophiles and thermophiles, respectively) are not all symmetrical; for example, more aligned sites have lysine in mesophiles and arginine in thermophiles than have the opposite pattern. This is generally taken to indicate that certain amino acids are favored over others by selection at different temperatures. Previous comparisons of protein sequences from mesophiles and thermophiles have used relatively small numbers of sequences from a diverse array of species, meaning that only the most common amino acid substitutions could be examined and any taxon-specific patterns would be obscured. Here, we compare a large number of proteins between mesophiles and thermophiles in the archaeal genus Methanococcus and the bacterial genus Bacillus. Each genus exhibits dramatically asymmetrical substitution patterns for many pairs of amino acids. There are several pairs of amino acids for which one amino acid is favored in thermophilic Bacillus and the other is favored in thermophilic Methanococcus; this appears to result from the higher G + C content of the DNA of thermophilic Bacillus, a complication not seen in Methanococcus.     

Hydrophobic environment is a key factor for the stability of thermophilic proteins

The stability of thermophilic proteins has been viewed from different perspectives and there is yet no unified principle to understand this stability. It would be valuable to reveal the most important interactions for designing thermostable proteins for such applications as industrial protein engineering. In this work, we have systematically analyzed the importance of various interactions by computing different parameters such as surrounding hydrophobicity, inter‐residue interactions, ion‐pairs and hydrogen bonds. The importance of each interaction has been determined by its predicted relative contribution in thermophiles versus the same contribution in mesophilic homologues based on a dataset of 373 protein families. We predict that hydrophobic environment is the major factor for the stability of thermophilic proteins and found that 80% of thermophilic proteins analyzed showed higher hydrophobicity than their mesophilic counterparts. Ion pairs, hydrogen bonds, and interaction energy are also important and favored in 68%, 50%, and 62% of thermophilic proteins, respectively. Interestingly, thermophilic proteins with decreased hydrophobic environments display a greater number of hydrogen bonds and/or ion pairs. The systematic elimination of mesophilic proteins based on surrounding hydrophobicity, interaction energy, and ion pairs/hydrogen bonds, led to correctly identifying 95% of the thermophilic proteins in our analyses. Our analysis was also applied to another, more refined set of 102 thermophilic–mesophilic pairs, which again identified hydrophobicity as a dominant property in 71% of the thermophilic proteins. Further, the notion of surrounding hydrophobicity, which characterizes the hydrophobic behavior of residues in a protein environment, has been applied to the three‐dimensional structures of elongation factor‐Tu proteins and we found that the thermophilic proteins are enriched with a hydrophobic environment. The results obtained in this work highlight the importance of hydrophobicity as the dominating characteristic in the stability of thermophilic proteins, and we anticipate this will be useful in our attempts to engineering thermostable proteins. © Proteins 2013. © 2012 Wiley Periodicals, Inc.     

The role of leucine and isoleucine in tuning the hydropathy of class A GPCRs

Leucine and Isoleucine are two amino acids that differ only by the positioning of one methyl group. This small difference can have important consequences in α-helices, as the β-branching of Ile results in helix destabilization. We set out to investigate whether there are general trends for the occurrences of Leu and Ile residues in the structures and sequences of class A GPCRs (G protein-coupled receptors). GPCRs are integral membrane proteins in which α-helices span the plasma membrane seven times and which play a crucial role in signal transmission. We found that Leu side chains are generally more exposed at the protein surface than Ile side chains. We explored whether this difference might be attributed to different functions of the two amino acids and tested if Leu tunes the hydrophobicity of the transmembrane domain based on the Wimley-White whole-residue hydrophobicity scales. Leu content decreases the variation in hydropathy between receptors and correlates with the non-Leu receptor hydropathy. Both measures indicate that hydropathy is tuned by Leu. To test this idea further, we generated protein sequences with random amino acid compositions using a simple numerical model, in which hydropathy was tuned by adjusting the number of Leu residues. The model was able to replicate the observations made with class A GPCR sequences. We speculate that the hydropathy of transmembrane domains of class A GPCRs is tuned by Leu (and to some lesser degree by Lys and Val) to facilitate correct insertion into membranes and/or to stably anchor the receptors within membranes.     

Branched-chain amino acid metabolism in higher plants

Valine, leucine and isoleucine contain short branched carbohydrate residues responsible for their classification as branched-chain amino acids (BCAA). Among the proteinogenic amino acids, BCAA show the highest hydrophobicity and are accordingly the major constituents of transmembrane regions of membrane proteins. BCAA cannot be synthesized by humans and thus belong to the essential amino acids. In contrast, plants are able to synthesize these amino acids de novo and are an important source for these compounds in the human diet. However, BCAA cannot only be synthesized in plants, leucine and probably also valine and isoleucine can also be degraded. Many enzymes operating in turnover are found in mitochondria, while some catabolizing activities are located in peroxisomes. The breakdown of BCAA is physically separated from their biosynthesis in chloroplasts. Additionally, in the order of the Capparales, enzymes of the leucine metabolism seem to be evolutionary related to or may even participate in the methionine chain elongation pathway, the early part of the biosynthesis of aliphatic glucosinolates. In summary, in higher plants a complex network of pathways interferes with the homeostasis of Val, Leu and Ile.     

Protein surface amino acid compositions distinctively differ between thermophilic and mesophilic bacteria

One of the well-known observations of proteins from thermophilic bacteria is the bias of the amino acid composition in which charged residues are present in large numbers, and polar residues are scarce. On the other hand, it has been reported that the molecular surfaces of proteins are adapted to their subcellular locations, in terms of the amino acid composition. Thus, it would be reasonable to expect that the differences in the amino acid compositions between proteins of thermophilic and mesophilic bacteria would be much greater on the protein surface than in the interior. We performed systematic comparisons between proteins from thermophilic bacteria and mesophilic bacteria, in terms of the amino acid composition of the protein surface and the interior, as well as the entire amino acid chains, by using sequence information from the genome projects. The biased amino acid composition of thermophilic proteins was confirmed, and the differences from those of mesophilic proteins were most obvious in the compositions of the protein surface. In contrast to the surface composition, the interior composition was not distinctive between the thermophilic and mesophilic proteins. The frequency of the amino acid pairs that are closely located in the space was also analyzed to show the same trend of the single amino acid compositions. Interestingly, extracellular proteins from mesophilic bacteria showed an inverse trend against thermophilic proteins (i.e. a reduced number of charged residues and rich in polar residues). Nuclear proteins from eukaryotes, which are known to be abundant in positive charges, showed different compositions as a whole from the thermophiles. These results suggest that the bias of the amino acid composition of thermophilic proteins is due to the residues on the protein surfaces, which may be constrained by the extreme environment.