Protein psychrophilicity: role of residual structural properties in adaptation of proteins to low temperatures

In order to investigate the structural distribution responsible for protein psychrophilicity, a systematic comparative analysis of 13 pairs of psychrophilic and mesophilic proteins is reported. Three kinds of residue structural states such as exposed, intermediate and buried were considered for analyzing the structural patterns of single amino acids and amino acids in different groups. The statistical test revealed that higher frequency in exposed state of Ala, higher frequency in intermediate state of His, lower frequency in buried state of Lys, lower frequency in exposed state of Gln, higher frequency in exposed state and in intermediate state of Thr, higher frequency in exposed and intermediate state of tiny and small amino acids groups could be critical factors related with protein psychrophilicity. Such structure-based differences of residual properties would help to develop a strategy for designing psychrophilic proteins.     

Protein thermostability: structure-based difference of residual properties between thermophilic and mesophilic proteins

Structure-based differences of residual properties between 20 pairs of thermophilic and mesophilic proteins were statistically analyzed to elucidate the factors governing protein thermostability. This study analyzed the distributions of outer residues, inner residues, flexible residues, rigid residues, hydrogen bonds, salt bridges, cation–pi interactions, and disulfide bonds in each protein in terms of residual structural states, which were determined as five kinds of states under residual packing value. Their structural patterns found in thermophilic protein groups were compared with those of mesophilic protein groups for showing distinctive difference of residual properties. The results of statistical tests (t-test) revealed that flexible residues in fully-exposed state and boundary state, salt bridges in exposed state, and hydrogen bonds in well-buried state could be critical factors related with protein thermostability. Such structure-based differences of residual properties would help to develop a strategy for enhancing protein thermostability.     

嗜压蛋白氨基酸残基溶剂可及性对其稳定性的影响

【目的】研究嗜压微生物中蛋白质在不同结构区域的特点对了解其稳定的结构基础及设计新型嗜压酶具有重要意义。【方法】利用43对嗜压和非嗜压同源蛋白的晶体结构信息,将氨基酸所处结构分成3种:分子表面、中间区及内核区,统计了不同结构状态中氨基酸的差异。【结果】统计分析结果表明,嗜压蛋白具有明显的溶剂可及性结构特点:其分子表面Cys、Asp、Asn和Lys含量显著高于非嗜压蛋白,而Pro和Arg则相反;中间区域Ile、Met含量显著高,而Trp则相反;内核区Cys、Ile含量显著高,而Ala则相反。同时,非嗜压氨基酸在嗜压蛋白分子表面及内核区含量均显著高于非嗜压蛋白。【结论】嗜压和非嗜压蛋白在分子表面差异最为明显,将是改造嗜压酶的首选区域;同时,需要统计更多的样本,对氨基酸压力不对称指数进行修订。     

Dependence of alpha-helical and beta-sheet amino acid propensities on the overall protein fold type

ABSTRACT: BACKGROUND: A large number of studies have been carried out to obtain amino acid propensities for alpha- helices and beta-sheets. The obtained propensities for alpha-helices are consistent with each other, and the pair-wise correlation coefficient is frequently high. On the other hand, the beta-sheet propensities obtained by several studies differed significantly, indicating that the context significantly affects beta-sheet propensity. RESULTS: We calculated amino acid propensities for alpha-helices and beta-sheets for 39 and 24 protein folds, respectively, and addressed whether they correlate with the fold. The propensities were also calculated for exposed and buried sites, respectively. Results showed that alpha-helix propensities do not differ significantly by fold, but beta-sheet propensities are diverse and depend on the fold. The propensities calculated for exposed sites and buried sites are similar for alpha-helix, but such is not the case for the beta-sheet propensities. We also found some fold dependence on amino acid frequency in beta-strands. Folds with a high Ser, Thr and Asn content at exposed sites in beta-strands tend to have a low Leu, Ile, Glu, Lys and Arg content (correlation coefficient = 0.90) and to have flat beta-sheets. At buried sites in beta-strands, the content of Tyr, Trp, Gln and Ser correlates negatively with the content of Val, Ile and Leu (correlation coefficient = 0.93). "All-beta" proteins tend to have a higher content of Tyr, Trp, Gln and Ser, whereas alpha/beta proteins tend to have a higher content of Val, Ile and Leu. CONCLUSIONS: The alpha-helix propensities are similar for all folds and for exposed and buried residues. However, beta-sheet propensities calculated for exposed residues differ from those for buried residues, indicating that the exposed-residue fraction is one of the major factors governing amino acid composition in beta-strands. Furthermore, the correlations we detected suggest that amino acid composition is related to folding properties such as the twist of a beta-strand or association between two beta sheets.     

Structural and mechanistic investigations on Salmonella typhimurium acetate kinase (AckA): identification of a putative ligand binding pocket at the dimeric interface

Background

A large number of studies have been carried out to obtain amino acid propensities for ??-helices and ??-sheets. The obtained propensities for ??-helices are consistent with each other, and the pair-wise correlation coefficient is frequently high. On the other hand, the ??-sheet propensities obtained by several studies differed significantly, indicating that the context significantly affects ??-sheet propensity.

Results

We calculated amino acid propensities for ??-helices and ??-sheets for 39 and 24 protein folds, respectively, and addressed whether they correlate with the fold. The propensities were also calculated for exposed and buried sites, respectively. Results showed that ??-helix propensities do not differ significantly by fold, but ??-sheet propensities are diverse and depend on the fold. The propensities calculated for exposed sites and buried sites are similar for ??-helix, but such is not the case for the ??-sheet propensities. We also found some fold dependence on amino acid frequency in ??-strands. Folds with a high Ser, Thr and Asn content at exposed sites in ??-strands tend to have a low Leu, Ile, Glu, Lys and Arg content (correlation coefficient = ?0.90) and to have flat ??-sheets. At buried sites in ??-strands, the content of Tyr, Trp, Gln and Ser correlates negatively with the content of Val, Ile and Leu (correlation coefficient = ?0.93). "All-??" proteins tend to have a higher content of Tyr, Trp, Gln and Ser, whereas "??/??" proteins tend to have a higher content of Val, Ile and Leu.

Conclusions

The ??-helix propensities are similar for all folds and for exposed and buried residues. However, ??-sheet propensities calculated for exposed residues differ from those for buried residues, indicating that the exposed-residue fraction is one of the major factors governing amino acid composition in ??-strands. Furthermore, the correlations we detected suggest that amino acid composition is related to folding properties such as the twist of a ??-strand or association between two ?? sheets.     

Group II archaeal chaperonin recognition of partially folded human γD‐crystallin mutants

The features in partially folded intermediates that allow the group II chaperonins to distinguish partially folded from native states remain unclear. The archaeal group II chaperonin from Methanococcus Mauripaludis (Mm‐Cpn) assists the in vitro refolding of the well‐characterized β‐sheet lens protein human γD‐crystallin (HγD‐Crys). The domain interface and buried cores of this Greek key conformation include side chains, which might be exposed in partially folded intermediates. We sought to assess whether particular features buried in the native state, but absent from the native protein surface, might serve as recognition signals. The features tested were (a) paired aromatic side chains, (b) side chains in the interface between the duplicated domains of HγD‐Crys, and (c) side chains in the buried core which result in congenital cataract when substituted. We tested the Mm‐Cpn suppression of aggregation of these HγD‐Crys mutants upon dilution out of denaturant. Mm‐Cpn was capable of suppressing the off‐pathway aggregation of the three classes of mutants indicating that the buried residues were not recognition signals. In fact, Mm‐Cpn recognized the HγD‐Crys mutants better than (wild‐type) WT and refolded most mutant HγD‐Crys to levels twice that of WT HγD‐Crys. This presumably represents the increased population or longer lifetimes of the partially folded intermediates of the mutant proteins. The results suggest that Mm‐Cpn does not recognize the features of HγD‐Crys tested—paired aromatics, exposed domain interface, or destabilized core—but rather recognizes other features of the partially folded β‐sheet conformation that are absent or inaccessible in the native state of HγD‐Crys.     

Topography of human plasma alpha1-acid glycoprotein.

Human plasma alpha1-acid glycoprotein, whose linear amino acid sequence has recently been elucidated (Schmid et al. (1973), Biochemistry 12, 2711), was further investigated with regard to its topography. Nitration of this protein and subsequent elucidation of the structures of the peptides containing modified tyrosine indicated that residues 27, 37, 78, 115, 127, and 157 are free, 50 and 91 are in an intermediate state, and 65, 74, 110, and 142 are buried. CD measurements between pH 10 and 12 demonstrated that the buried tyrosines are strongly hydrogen bonded and are probably responsible to a considerable extent for the stability of this protein. Of the three tryptophans of this protein, residue 122 proved to be partially reactive with Koshland reagent while the other two (25 and 160) were found to be unreactive. The state of the two disulfide bonds, established by differential reduction and alkylation with specific reagents, was shown to be of an intermediate type. Using carboxymethylation with bromoacetate at pH 7.0 for 8 days, the three histidines (97, 100, and 171) and methionine 111 could be shown to be in intermediate states. All lysines were treated with trinitrobenzenesulfonate and thus were assumed to be free. Of the 40 carboxylic groups, which were amidated with glycine methyl ester, 32 including the 14 sialyl residues were found to be free, six in an intermediate and the remaining two in a buried state. The present study describes the states of almost half of the amino acid residues of alpha1-acid glycoprotein, a knowledge important for the construction of a preliminary three-dimensional model of this conjugated protein.     

Hydration effect on low-frequency protein dynamics observed in simulated neutron scattering spectra

Hydration effects on protein dynamics were investigated by comparing the frequency dependence of the calculated neutron scattering spectra between full and minimal hydration states at temperatures between 100 and 300 K. The protein boson peak is observed in the frequency range 1-4 meV at 100 K in both states. The peak frequency in the minimal hydration state shifts to lower than that in the full hydration state. Protein motions with a frequency higher than 4 meV were shown to undergo almost harmonic motion in both states at all temperatures simulated, whereas those with a frequency lower than 1 meV dominate the total fluctuations above 220 K and contribute to the origin of the glass-like transition. At 300 K, the boson peak becomes buried in the quasielastic contributions in the full hydration state but is still observed in the minimal hydration state. The boson peak is observed when protein dynamics are trapped within a local minimum of its energy surface. Protein motions, which contribute to the boson peak, are distributed throughout the whole protein. The fine structure of the dynamics structure factor is expected to be detected by the experiment if a high resolution instrument (<∼20 μeV) is developed in the near future.     

Cold denaturation and heat denaturation of Streptomyces subtilisin inhibitor. 2. 1H NMR studies

Structural transitions of the protein Streptomyces subtilisin inhibitor (SSI) from the native state to the cold-denatured and heat-denatured states were studied by 1H NMR spectroscopy in the temperature range from -10 to 60 degrees C in the acidic pH range. Assignments of some of the 1H NMR signals of SSI in the cold-denatured and heat-denatured states were performed by a combined use of selective deuteration and site-directed mutagenesis. Throughout the pH range from 2.1 to 3.1, both transitions were cooperative and basically only three distinct spectra corresponding to structures in the cold-denatured, native, and heat-denatured states were detected. In the cold-denatured state, the side-chain signals of Met73, His106, at least one Val, and two Leu were observed at distinctly shifted positions from those for a random-coiled structure, suggesting the formation of a tertiary structure, while those of Met70, His43, and Ala2 were observed at positions for a random-coiled structure. This tertiary structure in the cold-denatured state is entirely different from that in the native state, as some amino acid residues exposed to the solvent in the native state (e.g., Met73, His106) are buried while those sequestered in the native state (e.g., His43) are exposed. In the heat-denatured state, however, most 1H NMR signals were observed at random-coiled positions, indicating that there is much less tertiary structure in the heat-denatured state than in the cold-denatured state. At pH values below 2.09, a structural transition was observed from the cold-denatured state to the heat-denatured state without passing through the native state.(ABSTRACT TRUNCATED AT 250 WORDS)     

Revisiting “reverse hydrophobic effect”: Applicable only to coil mutations at the surface

In a seminal paper, Pakula and Sauer (Nature, 1990, 344, 363–364) demonstrated that the increase in side‐chain hydrophobicity has a reverse relationship with protein stability. We have addressed this problem with several examples of mutants that span at different locations in protein structure based on secondary structure and solvent accessibility. We confirmed that the stability change upon single coil mutation at exposed region is reversely correlated with hydrophobicity with a single exception. In addition, we found the existence of such relationship in partially buried coil mutants. The stability of exposed helical mutants is governed by conformational properties. In buried and partially buried helical and strand mutants properties reflecting hydrophobicity have direct relationship with stability, whereas an opposite relationship was obtained with entropy and flexibility. The structural analysis of partially buried/exposed mutants showed that the surrounding residues are important for the stability change upon mutation. These results provide insights to understand the general behavior for the stability of proteins upon amino acid substitutions. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 591–599, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Role of medium- and long-range interactions to the stability of the mutants of T4 lysozyme

Inter-residue interactions play an important role to the folding and stability of protein molecules. In this work, we analyze the role of medium- and long-range interactions to the stability of T4 lysozyme mutants. We found that, in buried mutations, the increase in long-range contacts upon mutations destabilizes the protein, whereas, in surface mutations, the increase in long-range contacts increases the stability, indicating the importance of surrounding polar residues to the stability of surface mutations. Further, the increase in medium-range contacts decreases the stability of buried and surface mutations and a direct relationship is observed between the increase of medium-range contacts and increase in stability for partially buried/exposed mutations. Moreover, the relationship between amino acid properties and stability of T4 lysozyme mutants at positions Ile3, Phe53, and Leu99 showed that the effect of medium- and long-range contacts is less for buried mutations and the inter-residue contacts have significant correlation with the stability of partially buried mutations.     

RVP-net: online prediction of real valued accessible surface area of proteins from single sequences

SUMMARY: RVP-net is an online program for the prediction of real valued solvent accessibility. All previous methods of accessible surface area (ASA) predictions classify amino acid residues into exposure states and named them buried or exposed based on different thresholds. Real values in some cases were generated by taking the mid points of these state thresholds. This is the first method, which provides a direct prediction of ASA without making exposure categories and achieves results better than 19% mean absolute error. To facilitate batch processing of several sequences, a standalone version of this tool is also provided. AVAILABILITY: Online predictions are available at http://www.netasa.org/rvp-net/. Standalone version of the program can be obtained from the corresponding author by E-mail request.     

Analysis of temperature factor distribution in high-resolution protein structures.

The temperature factors obtained from X-ray refinement of proteins at high resolution show large variations from one structure to another. However, the B-values expressed in units of standard deviation about their mean value (B'-factor) at the C alpha atoms show remarkably characteristic frequency distribution. In all of the 110 proteins examined in this study, the frequency distribution exhibited a bimodal distribution. The peaks in the B'-factor frequency distribution occur at -1.1 and 0.4 for a bin size of 0.5. The peak at lower temperature factor corresponds largely to buried residues, whereas the peak at larger value corresponds to exposed residues. The distribution could be accurately described as a superposition of two Gaussian functions. The parameters describing the distribution are therefore characteristic of protein structures. The frequency distribution for a given amino acid over all the proteins also shows a similar bimodal distribution, although the areas under the two Gaussians differ from one amino acid to another. The area under the frequency distribution curve for any interval in B'-factor represents the propensity of the amino acid to occur in that interval. This propensity is related both to the hydrophilicity/hydrophobicity of the residue and the tendency of the residue to impose a different degree of rigidity on the polypeptide chain. The frequency distribution of stretches of high B'-factors departs appreciably from that expected for a random distribution. The correlation in the B-values of sequentially proximal residues is probably responsible for the bimodal distribution.     

Packing-based difference of structural features between thermophilic and mesophilic proteins

Twenty pairs of thermophilic and mesophilic proteins were compared in terms of residue packing distribution to obtain structural features related to protein thermostability. Based on residue packing concept, structural features of residues such as residue packing distribution, inner/outer position, secondary structure and water solvation were investigated. The statistical tests revealed that higher frequency in well-packed state of residues, lower frequency in exposed state and higher frequency in well-packed state of inner positioned residues, and higher frequency in well-packed state of 3/10 helix residues could be general structural features thermophilic proteins have.     

Functional dynamics of the hydrophobic cleft in the N-domain of calmodulin

Molecular dynamics studies of the N-domain (amino acids 1-77; CaM(1-77)) of Ca2+-loaded calmodulin (CaM) show that a solvent exposed hydrophobic cleft in the crystal structure of CaM exhibits transitions from an exposed (open) to a buried (closed) state over a time scale of nanoseconds. As a consequence of burying the hydrophobic cleft, the R(g) of the protein is reduced by 1.5 A. Based on this prediction, x-ray scattering experiments were conducted on this domain over a range of concentrations. Models built from the scattering data show that the R(g) and general shape is consistent with the simulation studies of CaM(1-77). Based on these observations we postulate a model in which the conformation of CaM fluctuates between two different states that expose and bury this hydrophobic cleft. In aqueous solution the closed state dominates the population, while in the presence of peptides, the open state dominates. This inherent flexibility of CaM may be the key to its versatility in recognizing structurally distinct peptide sequences. This model conflicts with the currently accepted hypothesis based on observations in the crystal structure, where upon Ca2+ binding the hydrophobic cleft is exposed to solvent. We postulate that crystal packing forces stabilize the protein conformation toward the open configuration.     

Predicting surface exposure of amino acids from protein sequence

The amino acid residues on a protein surface play a key role in interaction with other molecules, determined many physical properties, and constrain the structure of the folded protein. A database of monomeric protein crystal structures was used to teach computer-simulated neural networks rules for predicting surface exposure from local sequence. These trained networks are able to correctly predict surface exposure for 72% of residues in a testing set using a binary model, (buried/exposed) and for 54% of residues using a ternary model (buried/intermediate/exposed). In the ternary model, only 11% of the exposed residues are predicted as buried and only 5% of the buried residues are predicted as exposed. Also, since the networks are able to predict exposure with a quantitative confidence estimate, it is possible to assign exposure for over half of the residues in a binary model with greater than 80% accuracy. Even more accurate predictions are obtained by making a consensus prediction of exposure for a homologous family. The effect of the local environment of an amino acid on its accessibility, though smaller than expected, is significant and accounts for the higher success rate of prediction than obtained with previously used criteria. In the absence of a three-dimensional structure, the ability to predict surface accessibility of amino acids directly from the sequence is a valuable tool in choosing sites of chemical modification or specific mutations and in studies of molecular interaction.     

Solvent accessibility study on tRNAPhe

In order to assess the solvent–solute association in the tRNA^Phe molecule, solvent accessibility calculations were carried out for its crystalline and completely extended states following the method of Lee and Richards. To do this, results from the calculations on model trinucleotide systems pApXpA with different bases at position X were used. In the folded form of the molecule, it was found that the oxygen atoms O(I) and O(II) of almost all the phosphate groups and the O(2′) atoms of the sugar rings situated throughout the backbone were highly exposed to the solvent. The amount of reduction found in the solvent accessibilities of the various atoms in going from the extended state to the folded state of the molecule indicates the kind of compactness of the tertiary structure in tRNA^Phe. The results give quantitative support to many characteristics of the tRNA molecule, such as loop sections, buried/exposed residues, hydrophobic interactions, etc., which were thought to be due to other factors.     

Prediction of protein mutant stability using classification and regression tool

Prediction of protein stability upon amino acid substitutions is an important problem in molecular biology and the solving of which would help for designing stable mutants. In this work, we have analyzed the stability of protein mutants using two different datasets of 1396 and 2204 mutants obtained from ProTherm database, respectively for free energy change due to thermal (DeltaDeltaG) and denaturant denaturations (DeltaDeltaG(H(2)O)). We have used a set of 48 physical, chemical energetic and conformational properties of amino acid residues and computed the difference of amino acid properties for each mutant in both sets of data. These differences in amino acid properties have been related to protein stability (DeltaDeltaG and DeltaDeltaG(H(2)O)) and are used to train with classification and regression tool for predicting the stability of protein mutants. Further, we have tested the method with 4 fold, 5 fold and 10 fold cross validation procedures. We found that the physical properties, shape and flexibility are important determinants of protein stability. The classification of mutants based on secondary structure (helix, strand, turn and coil) and solvent accessibility (buried, partially buried, partially exposed and exposed) distinguished the stabilizing/destabilizing mutants at an average accuracy of 81% and 80%, respectively for DeltaDeltaG and DeltaDeltaG(H(2)O). The correlation between the experimental and predicted stability change is 0.61 for DeltaDeltaG and 0.44 for DeltaDeltaG(H(2)O). Further, the free energy change due to the replacement of amino acid residue has been predicted within an average error of 1.08 kcal/mol and 1.37 kcal/mol for thermal and chemical denaturation, respectively. The relative importance of secondary structure and solvent accessibility, and the influence of the dataset on prediction of protein mutant stability have been discussed.     

The states of tyrosyl and tryptophyl residues in a protein proteinase inhibitor (Streptomyces subtilisin inhibitor

The states of tyrosyl and tryptophyl residues of a dimeric protein proteinase inhibitor, Streptomyces subtilisin inhibitor (Sato, S & Murao, S. (1973), Agric. Biol. Chem. 37, 1067) were studies by solvent perturbation difference spectroscopy with methanol, ethylene glycol, polyethylene glycol, and deuterium oxide as perturbants, and by spectrophotometric titration at alkaline pH. It appeared that all three tyrosyl residues per monomer of the inhibitor were exposed on the surface of the molecule, and their apparent pK values were estimated separately to be 9.58, 11.10, and 12.42. The single tryptophyl residue per monomer of the inhibitor appeared to be partially buried in the protein molecule.     

High polar content of long buried blocks of sequence in protein domains suggests selection against amyloidogenic non-polar sequences

Native protein structures achieve stability in part by burying hydrophobic side-chains. About 75% of all amino acid residues buried in protein interiors are non-polar. Buried residues are not uniformly distributed in protein sequences, but sometimes cluster as contiguous polypeptide stretches that run through the interior of protein domain structures. Such regions have an intrinsically high local sequence density of non-polar residues, creating a potential problem: local non-polar sequences also promote protein misfolding and aggregation into non-native structures such as the amyloid fibrils in Alzheimer's disease. Here we show that long buried blocks of sequence in protein domains of known structure have, on average, a lower content of non-polar amino acids (about 70%) than do isolated buried residues (about 80%). This trend is observed both in small and in large protein domains and is independent of secondary structure. Long, completely non-polar buried stretches containing many large side-chains are particularly avoided. Aspartate residues that are incorporated in long buried stretches were found to make fewer polar interactions than those in short stretches, hinting that they may be destabilizing to the native state. We suggest that evolutionary pressure is acting on non-native properties, causing buried polar residues to be placed at positions where they would break up aggregation-prone non-polar sequences, perhaps even at some cost to native state stability.