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.     

A probabilistic method to correlate ion pairs with protein thermostability

Recent developments in research on the stability of proteins - specifically, comparisons of the ion pairs of homologous structures - show that ion pairs potentially contribute to the thermostability of proteins. This study proposes a probabilistic Bayesian statistical method to efficiently predict the thermostability of proteins based on the properties of ion pairs. The experimental results suggest that the numbers, types and bond lengths of ion pairs can be used to predict with high accuracy (up to 80%) the thermostability of functionally similar proteins. The predictions have high precision (99%), especially for hyperthermophilic proteins. Results for proteins with differing functions also indicate that the number of ion pairs is related to the thermostability of proteins, and that predictions of thermostability can also be made for proteins with different functions.     

Fluctuations in ion pairs and their stabilities in proteins

This report investigates the effect of systemic protein conformational flexibility on the contribution of ion pairs to protein stability. Toward this goal, we use all NMR conformer ensembles in the Protein Data Bank (1) that contain at least 40 conformers, (2) whose functional form is monomeric, (3) that are nonredundant, and (4) that are large enough. We find 11 proteins adhering to these criteria. Within these proteins, we identify 22 ion pairs that are close enough to be classified as salt bridges. These are identified in the high-resolution crystal structures of the respective proteins or in the minimized average structures (if the crystal structures are unavailable) or, if both are unavailable, in the "most representative" conformer of each of the ensembles. We next calculate the electrostatic contribution of each such ion pair in each of the conformers in the ensembles. This results in a comprehensive study of 1,201 ion pairs, which allows us to look for consistent trends in their electrostatic contributions to protein stability in large sets of conformers. We find that the contributions of ion pairs vary considerably among the conformers of each protein. The vast majority of the ion pairs interconvert between being stabilizing and destabilizing to the structure at least once in the ensembles. These fluctuations reflect the variabilities in the location of the ion pairing residues and in the geometric orientation of these residues, both with respect to each other, and with respect to other charged groups in the remainder of the protein. The higher crystallographic B-factors for the respective side-chains are consistent with these fluctuations. The major conclusion from this study is that salt bridges observed in crystal structure may break, and new salt bridges may be formed. Hence, the overall stabilizing (or, destabilizing) contribution of an ion pair is conformer population dependent.     

Comparison of sequence-based and structure-based phylogenetic trees of homologous proteins: Inferences on protein evolution

Several studies based on the known three-dimensional (3-D) structures of proteins show that two homologous proteins with insignificant sequence similarity could adopt a common fold and may perform same or similar biochemical functions. Hence, it is appropriate to use similarities in 3-D structure of proteins rather than the amino acid sequence similarities in modelling evolution of distantly related proteins. Here we present an assessment of using 3-D structures in modelling evolution of homologous proteins. Using a dataset of 108 protein domain families of known structures with at least 10 members per family we present a comparison of extent of structural and sequence dissimilarities among pairs of proteins which are inputs into the construction of phylogenetic trees. We find that correlation between the structure-based dissimilarity measures and the sequence-based dissimilarity measures is usually good if the sequence similarity among the homologues is about 30% or more. For protein families with low sequence similarity among the members, the correlation coefficient between the sequence-based and the structure-based dissimilarities are poor. In these cases the structure-based dendrogram clusters proteins with most similar biochemical functional properties better than the sequence-similarity based dendrogram. In multi-domain protein families and disulphide-rich protein families the correlation coefficient for the match of sequence-based and structure-based dissimilarity (SDM) measures can be poor though the sequence identity could be higher than 30%. Hence it is suggested that protein evolution is best modelled using 3-D structures if the sequence similarities (SSM) of the homologues are very low.     

同义密码子的反常蛋白质二级结构偏好性

统计分析了 119种人蛋白质和 92种大肠杆菌蛋白质的mRNA序列和蛋白质二级结构的关系 .从二肽频数出发 ,研究了同义密码子使用对蛋白质二级结构的影响 ,证明其影响在 10 %到 2 0 %的量级 .对于人和大肠杆菌 ,在 90 %置信水平上 ,4 0 0对二肽中分别有 79对和 6 0对 ,在 95 %置信水平上 ,分别有 4 5对和 36对二肽的相应密码子二联体具有不同于氨基酸的反常二级结构偏好性 ,并且这种反常不能归因于随机涨落     

Crystal structure of highly thermostable glycerol kinase from a hyperthermophilic archaeon in a dimeric form

The crystal structure of glycerol kinase from the hyperthermophilic archaeon Thermococcus kodakaraensis (Tk-GK) in a dimeric form was determined at a resolution of 2.4 A. This is the first crystal structure of a hyperthermophilic glycerol kinase. The overall structure of the Tk-GK dimer is very similar to that of the Escherichia coli glycerol kinase (Ec-GK) dimer. However, two dimers of Ec-GK can associate into a tetramer with a twofold axis, whereas those of Tk-GK cannot. This may be the reason why Tk-GK is not inhibited by fructose 1,6-bisphosphate, because the fructose 1,6-bisphosphate binding site is produced only when a tetrameric structure is formed. Differential scanning calorimetry analyses indicate that Tk-GK is a highly thermostable protein with a melting temperature (T(m)) of 105.4 degrees C for the major transition. This value is higher than that of Ec-GK by 34.1 degrees C. Comparison of the crystal structures of Tk-GK and Ec-GK indicate that there is a marked difference in the number of ion pairs in the alpha16 helix. Four ion pairs, termed IP1-IP4, are formed in this helix in the Tk-GK structure. To examine whether these ion pairs contribute to the stabilization of Tk-GK, four Tk-GK and four Ec-GK derivatives with reciprocal mutations at the IP1-IP4 sites were constructed. The determination of their stabilities indicates that the removal of each ion pair does not affect the stability of Tk-GK significantly, whereas the mutations designed to introduce one of these ion pairs stabilize or destabilize Ec-GK considerably. These results suggest that the ion pairs in the alpha16 helix contribute to the stabilization of Tk-GK in a cooperative manner.     

Ion pairs in alpha helices

A survey of 47 globular proteins was made to determine the probability of occurrence of ion pairs separated by 1,2,3,... and 8 residues in the alpha helices. As a control, the probability of occurrence of like charged pairs was also determined. The survey showed that ion pairs of the type i,i +/- 3 and i, i +/- 4 are the most predominant. Such a preference was not observed for like charged pairs. The observed frequency of ion pairs is significantly greater than their expected frequency. The normalized frequencies of occurrence of the ion pairs were also found to increase generally with the helix length. These results indicate that the ion pairs may contribute to the stability of solvent-exposed alpha helices. Since the stabilization of protein secondary structure enhances the stability of protein tertiary structure, these results may throw light on the mechanism of protein folding.     

Stereochemical modeling of disulfide bridges. Criteria for introduction into proteins by site-directed mutagenesis

A computer modeling procedure for assessing the stereochemical suitability of pairs of residues in proteins as potential sites for introduction of cystine disulfide crosslinks has been developed. Residue pairs with C alpha-C alpha distances of less than or equal to 6.5 A and C beta-C beta distances of less than or equal to 4.5 A are chosen for geometrical fixation of S atoms using the program MODIP. The stereochemistry of the modeled disulfides is evaluated using limits for the structural parameters of the various torsion angles and S-S bond length in the disulfide bridge. The ability of the procedure to correctly model disulfides has been checked with examples of cystine peptides of known crystal structures and 103 disulfide bridges from 25 available protein crystal structures determined at less than or equal to 2 A resolution. An analysis of results on three proteins with engineered disulfides, T4 lysozyme, dihydrofolate reductase and subtilisin, is presented. Two positions for the introduction of 'stereochemically optimal' disulfides are identified in subtilisin.     

Role of water molecules and ion pairs in Dps and related ferritin-like structures

A comparative study of water molecules and ion pairs in 11 Dps protein structures has been carried out. The invariant and common water molecules, the conserved residues interacting with them and the conserved ion pairs have been analyzed. Certain water molecules found on the interfaces between subunits are highly conserved and may be implicated in flexibility or continuing association of the subunits of the structure. It is possible that the water molecules, ion pairs and the special case of a water mediated charged network through a single water molecule are involved in maintaining the stability of the protein.     

Terminal ion pairs stabilize the second beta-hairpin of the B1 domain of protein G

The effects of terminal ion pairs on the stability of a beta-hairpin peptide corresponding to the C-terminal residues of the B1 domain of protein G were determined using thermal unfolding as monitored by nuclear magnetic resonance and circular dichroism spectroscopy. Molecular dynamics (MD) simulations were also performed to examine the effect of ion pairs on the structures. Eight peptides were studied including the wild type (G41) and the N-terminal modified sequences that had the first residue deleted (E42), replaced with a Lys (K41), or extended by an additional Gly (G40). Acetylated variants were made to examine the effect of removing the positive N-terminal charge on beta-hairpin stability. The rank in stability determined experimentally is K41 > E42 approximately G41 approximately G40 > Ac-K41 > Ac-E42 approximately Ac-G41 > Ac-G40. The Tm of the K41 peptide is 12 degrees C higher than G41, while the Tm values for the acetylated peptides are less than their unacetylated forms by more than 15 degrees C. NOE cross-peaks between side-chain methylene groups at the N- and C-termini and larger CalphaH shifts compared to random values are seen for K41. The addition of 20% methanol increases the stability in K41 and G41. The MD studies complement these results by showing that the charged N-terminus is important to stability. The type of ion pair observed varies with peptide, and when formed the simulations show that the ion pair can prevent fraying of the beta-strands through electrostatic and hydrophobic contacts. Therefore, introducing favorable electrostatic interactions at the N- and C-termini can substantially enhance beta-hairpin stability and help define the structure.     

Comparison of X‐ray and NMR structures: Is there a systematic difference in residue contacts between X‐ray‐ and NMR‐resolved protein structures?

We have compared structures of 78 proteins determined by both NMR and X-ray methods. It is shown that X-ray and NMR structures of the same protein have more differences than various X-ray structures obtained for the protein, and even more than various NMR structures of the protein. X-ray and NMR structures of 18 of these 78 proteins have obvious large-scale structural differences that seem to reflect a difference of crystal and solution structures. The other 60 pairs of structures have only small-scale differences comparable with differences between various X-ray or various NMR structures of a protein; we have analyzed these structures more attentively. One of the main differences between NMR and X-ray structures concerns the number of contacts per residue: (1) NMR structures presented in PDB have more contacts than X-ray structures at distances below 3.0 A and 4.5-6.5 A, and fewer contacts at distances of 3.0-4.5 A and 6.5-8.0 A; (2) this difference in the number of contacts is greater for internal residues than for external ones, and it is larger for beta-containing proteins than for all-alpha proteins. Another significant difference is that the main-chain hydrogen bonds identified in X-ray and NMR structures often differ. Their correlation is 69% only. However, analogous difference is found for refined and rerefined NMR structures, allowing us to suggest that the observed difference in interresidue contacts of X-ray and NMR structures of the same proteins is due mainly to a difference in mathematical treatment of experimental results.     

Ion pairs and their role in modulating stability of cold- and warm-active uracil DNA glycosylase

MD simulations and continuum electrostatics calculations have been used to study the observed differences in thermostability of cold- and warm-active uracil DNA glycosylase (UDG). The present study focuses on the role of ion pairs and how they affect the thermal stability of the two enzymes. Analysis of the MD generated structural ensembles show that cod UDG (cUDG) and human UDG (hUDG) have 11 and 12 ion pairs which are present in at least 30% of the conformations. The electrostatic contribution of the ion pairs, computed using continuum electrostatics, is slightly more favorable in cUDG at 298 K. This is primarily attributed to more optimized interactions between the ion pairs and nearby dipoles/charges in cUDG. More global salt bridges are found in hUDG and are more stabilizing when compared to cUDG, possibly playing a role in maintaining overall stability and reducing conformational entropy. Both enzymes contain one three-member ionic network, but the one found in hUDG is far more stabilizing. Our results also suggest that care should be taken when performing statistical analysis of crystal structures with respect to ion pairs, and that crystallization conditions must be carefully examined when performing such analysis.     

Current computer modeling cannot explain why two highly similar sequences fold into different structures

The remarkable recent creation of two proteins that fold into two completely different and stable structures, exhibit different functions, yet differ by only a few amino acids poses a conundrum to those hoping to understand how sequence encodes structure. Here, computer modeling uniquely allows the characterization of not only the native structure of each minimally different sequence but also systems in which each sequence was modeled onto the fold of the alternate sequence. The reasons for the different structural preferences of two pairs of highly similar sequences are explored by a combination of structure analyses, comparison of potential energies calculated from energy-minimized single structures and trajectories produced from molecular dynamics simulations, and application of a novel method for calculating free energy differences. The sensitivity of such analyses to the choice of force field is also explored. Many of the hypotheses proposed on the basis of the nuclear magnetic resonance model structures of the proteins with 95% identical sequences are supported. However, each level of analysis provides different predictions regarding which sequence-structure combination should be most favored, highlighting the fact that protein structure and stability result from a complex combination of interdependent factors.     

Energetics of charge-charge interactions between residues adjacent in sequence

Statistical analysis of the residue separation between a pair of ionizable side chains within 4 ? of each other was performed on a set of 1560 non-homologous PDB structures. We found that the frequency of pairs of like charges (i.e., pairs consisting of acidic residues Asp and Glu or pairs consisting of basic residues Arg and Lys) is two orders of magnitude lower than the pairs of oppositely charged residues (salt-bridges). We also found that for pairs of like charges the distribution is skewed dramatically towards short residue separation (<3). On the basis of these observations, we hypothesize that at short residue separation the repulsion between charges does not contribute much to the protein stability and the effects are largely dominated by the long range charge-charge interactions with other ionizable groups in the protein molecule. To test this hypothesis, we incorporated various pairs of charged residues at position 63 and 64 of ubiquitin and compared the stabilities of these variants. We also performed calculations of the expected changes in the charge-charge interactions. A very good correlation between experimental changes in the stability of ubiquitin variants, and changes in the energy of charge-charge interactions provides support for the hypothesis that a pair of ionizable residues next to each other in sequence modulates protein stability via long range charge-charge interactions with the rest of the protein.     

Comparative analyses of quaternary arrangements in homo‐oligomeric proteins in superfamilies: Functional implications

A comprehensive analysis of the quaternary features of distantly related homo‐oligomeric proteins is the focus of the current study. This study has been performed at the levels of quaternary state, symmetry, and quaternary structure. Quaternary state and quaternary structure refers to the number of subunits and spatial arrangements of subunits, respectively. Using a large dataset of available 3D structures of biologically relevant assemblies, we show that only 53% of the distantly related homo‐oligomeric proteins have the same quaternary state. Considering these homologous homo‐oligomers with the same quaternary state, conservation of quaternary structures is observed only in 38% of the pairs. In 36% of the pairs of distantly related homo‐oligomers with different quaternary states the larger assembly in a pair shows high structural similarity with the entire quaternary structure of the related protein with lower quaternary state and it is referred as “Russian doll effect.” The differences in quaternary state and structure have been suggested to contribute to the functional diversity. Detailed investigations show that even though the gross functions of many distantly related homo‐oligomers are the same, finer level differences in molecular functions are manifested by differences in quaternary states and structures. Comparison of structures of biological assemblies in distantly and closely related homo‐oligomeric proteins throughout the study differentiates the effects of sequence divergence on the quaternary structures and function. Knowledge inferred from this study can provide insights for improved protein structure classification and function prediction of homo‐oligomers. Proteins 2016; 84:1190–1202. © 2016 Wiley Periodicals, Inc.     

Electrostatic contributions to the stability of the GCN4 leucine zipper structure

Ion pairs are ubiquitous in X-ray structures of coiled coils, and mutagenesis of charged residues can result in large stability losses. By contrast, pK_a values determined by NMR in solution often predict only small contributions to stability from charge interactions. To help reconcile these results we used triple-resonance NMR to determine pK_a values for all groups that ionize between pH 1 and 13 in the 33 residue leucine zipper fragment, GCN4p. In addition to the native state we also determined comprehensive pK_a values for two models of the GCN4p denatured state: the protein in 6 M urea, and unfolded peptide fragments of the protein in water. Only residues that form ion pairs in multiple X-ray structures of GCN4p gave large pK_a differences between the native and denatured states. Moreover, electrostatic contributions to stability were not equivalent for oppositely charged partners in ion pairs, suggesting that the interactions between a charge and its environment are as important as those within the ion pair. The pH dependence of protein stability calculated from NMR-derived pK_a values agreed with the stability profile measured from equilibrium urea-unfolding experiments as a function of pH. The stability profile was also reproduced with structure-based continuum electrostatic calculations, although contributions to stability were overestimated at the extremes of pH. We consider potential sources of errors in the calculations, and how pK_a predictions could be improved. Our results show that although hydrophobic packing and hydrogen bonding have dominant roles, electrostatic interactions also make significant contributions to the stability of the coiled coil.     

Aromatic side-chain interactions in proteins. II. Near- and far-sequence Phe-X pairs

We have collected all aromatic pairs (3152) involving an N-phenyl partner in a dataset of 593 proteins of the PDB: 728 of these pairs involve a partner residue less than 6 apart in the sequence. These near-sequence Phe-X pairs correspond to specific conformations that stabilize secondary structures, mainly alpha-helices when the residues are 1, 3, and 4 apart, and beta-strands when they are 2 apart in the sequence. These conformations are not spatially random and have been examined in detail. The remaining phenylalanine pairs (2424) are between partners more than 5 apart in the sequence. Of these far-sequence pairs, 34% of occurrences are in sheets. Next in frequencies are pairs that bridge a beta-strand to a helix (24%), followed by pairs that bridge a beta-strand to a random coiled structure (15%). Helix to helix pairs only constitute 12% of these far-sequence pairs. Analysis of the pairing frequency supports the hypothesis that aromatic interactions are late events of protein folding.     

Implication for buried polar contacts and ion pairs in hyperthermostable enzymes

Understanding the structural basis of thermostability and thermoactivity, and their interdependence, is central to the successful future exploitation of extremophilic enzymes in biotechnology. However, the structural basis of thermostability is still not fully characterized. Ionizable residues play essential roles in proteins, modulating protein stability, folding and function. The dominant roles of the buried polar contacts and ion pairs have been reviewed by distinguishing between the inside polar contacts and the total intramolecular polar contacts, and by evaluating their contribution as molecular determinants for protein stability using various protein structures from hyperthermophiles, thermophiles and mesophilic organisms. The analysis revealed that the remarkably increased number of internal polar contacts in a monomeric structure probably play a central role in enhancing the melting temperature value up to 120 degrees C for hyperthermophilic enzymes from the genus Pyrococcus. These results provide a promising contribution for improving the thermostability of enzymes by modulating buried polar contacts and ion pairs.     

Contribution of cation-pi interactions to the stability of protein-DNA complexes

Cation-pi interactions between an aromatic ring and a positive charge located above it have proven to be important in protein structures and biomolecule associations. Here, the role of these interactions at the interface of protein-DNA complexes is investigated, by means of ab initio quantum mechanics energy calculations and X-ray structure analyses. Ab initio energy calculations indicate that Na ions and DNA bases can form stable cation-pi complexes, whose binding strength strongly depends on the type of base, on the position of the Na ion, and whether the base is isolated or included in a double-stranded B-DNA. A survey of protein-DNA complex structures using appropriate geometrical criteria revealed cation-pi interactions in 71% of the complexes. More than half of the cation-pi pairs involve arginine residues, about one-third asparagine or glutamine residues that only carry a partial charge, and one-seventh lysine residues. The most frequently observed pair, which is also the most stable as monitored by ab initio energy calculations, is arginine- guanine. Arginine-adenine interactions are also favorable in general, although to a lesser extent, whereas those with thymine and cytosine are not. Our calculations show that the major contribution to cation-pi interactions with DNA bases is of electrostatic nature. These interactions often occur concomitantly with hydrogen bonds with adjacent bases; their strength is estimated to be from three to four times lower than that of hydrogen bonds. Finally, the role of cation-pi interactions in the stability and specificity of protein-DNA complexes is discussed.     

Tolerance to the substitution of buried apolar residues by charged residues in the homologous protein structures

Occurrence and accommodation of charged amino acid residues in proteins that are structurally equivalent to buried non-polar residues in homologues have been investigated. Using a dataset of 1,852 homologous pairs of crystal structures of proteins available at 2A or better resolution, 14,024 examples of apolar residues in the structurally conserved regions replaced by charged residues in homologues have been identified. Out of 2,530 cases of buried apolar residues, 1,677 of the equivalent charged residues in homologues are exposed and the rest of the charged residues are buried. These drastic substitutions are most often observed in homologous protein pairs with low sequence identity (<30%) and in large protein domains (>300 residues). Such buried charged residues in the large proteins are often located in the interface of sub-domains or in the interface of structural repeats, Beyond 7A of residue depth of buried apolar residues, or less than 4% of solvent accessibility, almost all the substituting charged residues are buried. It is also observed that acidic sidechains have higher preference to get buried than the positively charged residues. There is a preference for buried charged residues to get accommodated in the interior by forming hydrogen bonds with another sidechain than the main chain. The sidechains interacting with a buried charged residue are most often located in the structurally conserved regions of the alignment. About 50% of the observations involving hydrogen bond between buried charged sidechain and another sidechain correspond to salt bridges. Among the buried charged residues interacting with the main chain, positively charged sidechains form hydrogen bonds commonly with main chain carbonyls while the negatively charged residues are accommodated by hydrogen bonding with the main chain amides. These carbonyls and amides are usually located in the loops that are structurally variable among homologous proteins.