Hydrophobic residues can identify native protein structures

Evaluation of protein structures needs a trustworthy potential function. Although several knowledge‐based potential functions exist, the impact of different types of amino acids in the scoring functions has not been studied yet. Previously, we have reported the importance of nonlocal interactions in scoring function (based on Delaunay tessellation) in discrimination of native structures. Then, we have questioned the structural impact of hydrophobic amino acids in protein fold recognition. Therefore, a Hydrophobic Reduced Model (HRM) was designed to reduce protein structure of FS (Full Structure) into RS (Reduced Structure). RS is considered as a reduced structure of only seven hydrophobic amino acids (L, V, F, I, A, W, Y) and all their interactions. The presented model was evaluated via four different performance metrics including the number of correctly identified natives, the Z‐score of the native energy, the RMSD of the minimum score, and the Pearson correlation coefficient between the energy and the model quality. Results indicated that only nonlocal interactions between hydrophobic amino acids could be sufficient and accurate enough for protein fold recognition. Interestingly, the results of HRM is significantly close to the model that considers all amino acids (20‐amino acid model) to discriminate the native structure of the proteins on eleven decoy sets. This indicates that the power of knowledge‐based potential functions in protein fold recognition is mostly due to hydrophobic interactions. Hence, we suggest combining a different well‐designed scoring function for non‐hydrophobic interactions with HRM to achieve better performance in fold recognition.     

Geometric cooperativity and anticooperativity of three-body interactions in native proteins

Characterizing multibody interactions of hydrophobic, polar, and ionizable residues in protein is important for understanding the stability of protein structures. We introduce a geometric model for quantifying 3-body interactions in native proteins. With this model, empirical propensity values for many types of 3-body interactions can be reliably estimated from a database of native protein structures, despite the overwhelming presence of pairwise contacts. In addition, we define a nonadditive coefficient that characterizes cooperativity and anticooperativity of residue interactions in native proteins by measuring the deviation of 3-body interactions from 3 independent pairwise interactions. It compares the 3-body propensity value from what would be expected if only pairwise interactions were considered, and highlights the distinction of propensity and cooperativity of 3-body interaction. Based on the geometric model, and what can be inferred from statistical analysis of such a model, we find that hydrophobic interactions and hydrogen-bonding interactions make nonadditive contributions to protein stability, but the nonadditive nature depends on whether such interactions are located in the protein interior or on the protein surface. When located in the interior, many hydrophobic interactions such as those involving alkyl residues are anticooperative. Salt-bridge and regular hydrogen-bonding interactions, such as those involving ionizable residues and polar residues, are cooperative. When located on the protein surface, these salt-bridge and regular hydrogen-bonding interactions are anticooperative, and hydrophobic interactions involving alkyl residues become cooperative. We show with examples that incorporating 3-body interactions improves discrimination of protein native structures against decoy conformations. In addition, analysis of cooperative 3-body interaction may reveal spatial motifs that can suggest specific protein functions.     

Trapping the on-pathway folding intermediate of Im7 at equilibrium

The four-helical protein Im7 folds via a rapidly formed on-pathway intermediate (k(UI)=3000 s(-1) at pH 7.0, 10 degrees C) that contains three (helices I, II and IV) of the four native alpha-helices. The relatively slow (k(IN)=300 s(-1)) conversion of this intermediate into the native structure is driven by the folding and docking of the six residue helix III onto the developing hydrophobic core. Here, we describe the structural properties of four Im7* variants designed to trap the protein in the intermediate state by disrupting the stabilising interactions formed between helix III and the rest of the protein structure. In two of these variants (I54A and L53AI54A), hydrophobic residues within helix III have been mutated to alanine, whilst in the other two mutants the sequence encompassing the native helix III was replaced by a glycine linker, three (H3G3) or six (H3G6) residues in length. All four variants were shown to be monomeric, as judged by analytical ultracentrifugation, and highly helical as measured by far-UV CD. In addition, all the variants denature co-operatively and have a stability (DeltaG(UF)) and buried hydrophobic surface area (M(UF)) similar to those of the on-pathway kinetic intermediate. Structural characterisation of these variants using 1-anilino-8-napthalene sulphonic acid (ANS) binding, near-UV CD and 1D (1)H NMR demonstrate further that the trapped intermediate ensemble is highly structured with little exposed hydrophobic surface area. Interestingly, however, the structural properties of the variants I54A and L53AI54A differ in detail from those of H3G3 and H3G6. In particular, the single tryptophan residue, located near the end of helix IV, and distant from helix III, is in a distinct environment in the two sets of mutants as judged by fluorescence, near-UV CD and the sensitivity of tryptophan fluorescence to iodide quenching. Overall, the results confirm previous kinetic analysis that demonstrated the hierarchical folding of Im7 via an on-pathway intermediate, and show that this species is a highly helical ensemble with a well-formed hydrophobic core. By contrast with the native state, however, the intermediate ensemble is flexible enough to change in response to mutation, its structural properties being tailored by residues in the sequence encompassing the native helix III.     

Visualization of the nature of protein folding by a study of a distance constraint approach in two-dimensional models

A distance constraint approach is applied to two-dimensional models of proteins in order to visualize the nature of protein folding and to examine the relative roles of different ranges of interaction. Three different native structures (I, II, and III) are considered; they have two different kinds of residues, viz., hydrophobic and hydrophilic, and different sequences of these residues. We examine how the distance constraint approach functions in the prediction of protein folding when we know the sequence of the residues, the (fixed) bond lengths, the mean distances between residues i and i + 2, and i and i + 3, and the mean distances for hydrophobic–hydrophobic, hydrophobic–hydrophilic, and hydrophilic–hydrophilic contacts between residues i and i + j, where j ≥ 4. This approach involves optimization of an object function with respect to 98 variables and is not free of the multiple-minimum problem. The optimization is always terminated if the chain is entangled and/or the segments (residues) are packed too compactly to move. In order to escape from such situations and to take the excluded-volume effect into account, a Monte Carlo method is used after the optimization is trapped in local minima. Success in the prediction of folding is found to depend on the starting conformations and on the native conformations. Fair success is obtained in predicting the helix-like structure in protein I and the overall structure of protein III, but not the β-like structures of proteins I and II. Insofar as the prediction of the structure of protein III is reasonable, it appears that some sequences of residues produce greater constraints on their conformations than others, if one considers only the hydrophobic and hydrophilic nature of the residues. These results imply that, in the folding of real proteins in three dimensions, the competition for hydrophobic (and hydrophilic) residues for inside (outside) positions in the molecule probably constitutes a necessary but not a sufficient condition to form and stabilize the native structure. The failure to predict the structure of protein II, and part of that of protein I, suggests that there are two types of long-range interactions. One (which we considered here) is nonspecific (i.e., is defined only in terms of contacts between residues of the same or different polarity) and acts at any stage of protein folding; the other (which we did not consider here) is a specific interaction between residues in pairs and contributes only when the residues in the specific pair take on the native conformation. Presumably, incorporation of such specific long-range interactions, together with the nonspecific ones, is necessary for successful protein folding, using the distance constraint approach.     

Polar residues in the protein core of Escherichia coli thioredoxin are important for fold specificity

Most globular proteins contain a core of hydrophobic residues that are inaccessible to solvent in the folded state. In general, polar residues in the core are thermodynamically unfavorable except when they are able to form intramolecular hydrogen bonds. Compared to hydrophobic interactions, polar interactions are more directional in character and may aid in fold specificity. In a survey of 263 globular protein structures, we found a strong positive correlation between the number of polar residues at core positions and protein size. To probe the importance of buried polar residues, we experimentally tested the effects of hydrophobic mutations at the five polar core residues in Escherichia coli thioredoxin. Proteins with single hydrophobic mutations (D26I, C32A, C35A, T66L, and T77V) all have cooperative unfolding transitions like the wild type (wt), as determined by chemical denaturation. Relative to wt, D26I is more stable while the other point mutants are less stable. The combined 5-fold mutant protein (IAALV) is less stable than wt and has an unfolding transition that is substantially less cooperative than that of wt. NMR spectra as well as amide deuterium exchange indicate that IAALV is likely sampling a number of low-energy structures in the folded state, suggesting that polar residues in the core are important for specifying a well-folded native structure.     

Determination of partial unfolding enthalpy for lysozyme upon interaction with dodecyltrimethylammonium bromide using an extended solvation model

The interactions of dodecyltrimethylammonium bromides (DTABs) with hen egg lysozyme have been investigated at pH = 7.0 and 27 degrees C in phosphate buffer by isothermal titration calorimetry. DTAB interacts endothermically and activate lysozyme. The endothermicity of the lysozyme-DTAB interaction is in marked contrast to the exothermic interactions between sodium dodecyl sulphate (SDS) and lysozyme which have been attributed to specific binding between the anionic sulphate head groups and cationic amino acid residues. The enthalpies of interaction between the cationic surfactant (DTAB) and lysozyme are dominated by the endothermic unfolding of the native structure followed by an exothermic solvation of the lysozyme-DTAB complex by the addition of extra DTAB. A new direct calorimetric method to follow protein denaturation, and the effect of surfactants on the stability of proteins was introduced. The extended solvation model was used to reproduce the enthalpies of lysozyme-DTAB interaction over the whole range of DTAB concentrations. The solvation parameters recovered from the new equation, attributed to the structural change of lysozyme and its biological activity. At low concentrations of DTAB, the binding is mainly electrostatic, with some simultaneous interaction of the hydrophobic tail with nearby hydrophobic patches on the lysozyme. These initial interactions presumably cause some protein unfolding and expose additional hydrophobic sites. The DTAB-induced denaturation enthalpy of lysozyme is 86.46 +/- 0.02 kJ mol(-1).     

Inter-residue and solvent-residue interactions in proteins: a statistical study on experimental structures

A large set of protein structures resolved by X-ray or NMR techniques has been extracted from the Protein Data Bank and analyzed using statistical methods. In particular, we investigate the interactions between side chains and the interactions between solvent and side chains, pointing out on the possibility of including the solvent as part of a knowledge-based potential. The solvent-residue contacts are accounted for on the basis of the Voronoi's polyhedron analysis. Our investigation confirms the importance of hydrophobic residues in determining the protein stability. We observe that in general hydrophobic-hydrophobic interactions and, more specifically, aromatic-aromatic contacts tend to be increasingly distally separated in the primary sequence of proteins, thus connecting distinct secondary structure elements. A simple relation expressing the dependence of the protein free energy by the number of residues is proposed. Such a relation includes both the residue-residue and the solvent-residue contributions. The former is dominant for large size proteins, whereas for small sizes (number of residues less than 100) the two terms are comparable. Gapless threading experiments show that the solvent-residue knowledge-based potential yields a significant contribution with respect to discriminating the native structure of proteins. Such contribution is important especially for proteins of small size and is similar to that given by the most favorable residue-residue knowledge-based potential referring to hydrophobic-hydrophobic interactions such as isoleucine-leucine. In general, the inclusion of the solvent-residue interaction produces a relevant increase of the free energy gap between the native structures and decoys.     

Use of hexadecyl Fractosil as a hydrophobic carrier for adsorptive immobilization of proteins

Fractosil, a porous form of silica, has been used for the preparation of a hydrophobically derivatized carrier for protein immobilization. Interaction of a number of arbitrarily chosen proteins with hexadecyl-substituted Fractosil has been investigated. Binding of proteins was found to take place with retention of their native properties. Glutamate dehydrogenase, used as a model allosteric protein, was found to retain its catalytic and allosteric properties upon binding to the adsorbent in the form of suspension or column. Positive cooperative interactions for binding of bovine serum albumin and glutamate dehydrogenase to the matrix were observed. These findings are discussed in terms of hydrophobic interactions occurring between various residues of the protein molecules and the hydrophobic ligands in addition to those interactions which may occur with the unsubstituted gel. Results presented on immobilized glutamate dehydrogenase, trypsin, alpha-chymotrypsin, alpha-amylase, and amyloglucosidase clearly indicate possible potential of the support for continuous catalytic transformations.     

A protein dissection study demonstrates that two specific hydrophobic clusters play a key role in stabilizing the core structure of the molten globule state of human alpha-lactalbumin

The molten globule state of alpha-lactalbumin (alpha LA) has served as a paradigm for understanding the role of these partially folded states in protein folding. We previously showed that a peptide construct consisting of the A and B helices (residues 1-38) cross-linked to the D- and C-terminal 3(10) helices (residues 101-120) of alpha LA is capable of folding to a stable molten globule-like state. Here, we report the study of three peptide constructs that are designed to investigate the contribution two short hydrophobic sequences located near the C-terminus of alpha LA make to the structure and stability of the alpha LA molten globule state. These regions of the protein have been shown to form stable non-native structures in isolation. The three peptide constructs contain residues 1-38 cross-linked to three separate C-terminal peptides via the native 28-111 disulfide bond. The C-terminal peptides consist of residues 101-114, 106-120, and 106-114. The results of CD, fluorescence, ANS binding, and urea denaturation experiments indicate that constructs that lack either of the hydrophobic sequences (residues 101-105 and 115-120) are significantly less structured. These results highlight the importance of long-range, mutually stabilizing interactions within the molten globule state of the protein. Proteins 2001;42:237-242.     

Monte Carlo simulations of the peptide recognition at the consensus binding site of the constant fragment of human immunoglobulin G: the energy landscape analysis of a hot spot at the intermolecular interface

Monte Carlo simulations of molecular recognition at the consensus binding site of the constant fragment (Fc) of human immunoglobulin G (Ig) protein have been performed to analyze structural and thermodynamic aspects of binding for the 13-residue cyclic peptide DCAWHLGELVWCT. The energy landscape analysis of a hot spot at the intermolecular interface using alanine scanning and equilibrium-simulated tempering dynamics with the simplified, knowledge-based energy function has enabled the role of the protein hot spot residues in providing the thermodynamic stability of the native structure to be determined. We have found that hydrophobic interactions between the peptide and the Met-252, Ile-253, His-433, and His-435 protein residues are critical to guarantee the thermodynamic stability of the crystallographic binding mode of the complex. Binding free energy calculations, using a molecular mechanics force field and a solvation energy model, combined with alanine scanning have been conducted to determine the energetic contribution of the protein hot spot residues in binding affinity. The conserved Asn-434, Ser-254, and Tyr-436 protein residues contribute significantly to the binding affinity of the peptide-protein complex, serving as an energetic hot spot at the intermolecular interface. The results suggest that evolutionary conserved hot spot protein residues at the intermolecular interface may be partitioned in fulfilling thermodynamic stability of the native binding mode and contributing to the binding affinity of the complex.     

Characterization of a partially denatured state of a protein by two-dimensional NMR: reduction of the hydrophobic interactions in ubiquitin

A stable, partially structured state of ubiquitin, the A-state, is formed at pH 2.0 in 60% methanol/40% water at 298 K. Detailed characterization of the structure of this state has been carried out by 2D NMR spectroscopy. Assignment of slowly exchanging amide resonances protected from the solvent in the native and A-state shows that gross structural reorganization of the protein has not occurred and that the A-state contains a subset of the interactions present in the native state (N-state). Vicinal coupling constants and NOESY data show the presence of the first two strands of the five-strand beta-sheet that is present in the native protein and part of the third beta-strand. The hydrophobic face of the beta-sheet in the A-state is covered by a partially structured alpha-helix, tentatively assigned to residues 24-34, that is considerably more flexible than the alpha-helix in the N-state. There is evidence for some fixed side-chain--side-chain interactions between these two units of structure. The turn-rich area of the protein, which contains seven reverse turns and a short piece of 3(10) helix, does not appear to be structured in the A-state and is approaching random coil.     

Simplicial edge representation of protein structures and alpha contact potential with confidence measure

Protein representation and potential function are two important ingredients for studying protein folding, equilibrium thermodynamics, and sequence design. We introduce a novel geometric representation of protein contact interactions using the edge simplices from the alpha shape of the protein structure. This representation can eliminate implausible neighbors that are not in physical contact, and can avoid spurious contact between two residues when a third residue is between them. We developed statistical alpha contact potential using an odds-ratio model. A studentized bootstrap method was then introduced to assess the 95% confidence intervals for each of the 210 propensity parameters. We found, with confidence, that there is significant long-range propensity (>30 residues apart) for hydrophobic interactions. We tested alpha contact potential for native structure discrimination using several sets of decoy structures, and found that it often performs comparably with atom-based potentials requiring many more parameters. We also show that accurate geometric representation is important, and that alpha contact potential has better performance than potential defined by cutoff distance between geometric centers of side chains. Hierarchical clustering of alpha contact potentials reveals natural grouping of residues. To explore the relationship between shape and physicochemical representations, we tested the minimum alphabet size necessary for native structure discrimination. We found that there is no significant difference in performance of discrimination when alphabet size varies from 7 to 20, if geometry is represented accurately by alpha simplicial edges. This result suggests that the geometry of packing plays an important role, but the specific residue types are often interchangeable.     

De novo design of the hydrophobic cores of proteins.

We have developed and experimentally tested a novel computational approach for the de novo design of hydrophobic cores. A pair of computer programs has been written, the first of which creates a "custom" rotamer library for potential hydrophobic residues, based on the backbone structure of the protein of interest. The second program uses a genetic algorithm to globally optimize for a low energy core sequence and structure, using the custom rotamer library as input. Success of the programs in predicting the sequences of native proteins indicates that they should be effective tools for protein design. Using these programs, we have designed and engineered several variants of the phage 434 cro protein, containing five, seven, or eight sequence changes in the hydrophobic core. As controls, we have produced a variant consisting of a randomly generated core with six sequence changes but equal volume relative to the native core and a variant with a "minimalist" core containing predominantly leucine residues. Two of the designs, including one with eight core sequence changes, have thermal stabilities comparable to the native protein, whereas the third design and the minimalist protein are significantly destabilized. The randomly designed control is completely unfolded under equivalent conditions. These results suggest that rational de novo design of hydrophobic cores is feasible, and stress the importance of specific packing interactions for the stability of proteins. A surprising aspect of the results is that all of the variants display highly cooperative thermal denaturation curves and reasonably dispersed NMR spectra. This suggests that the non-core residues of a protein play a significant role in determining the uniqueness of the folded structure.     

Solvent accessibility in native and isolated domain environments: general features and implications to interface predictability

A non-redundant database of 4536 structural domains, comprising more than 790,000 residues, has been used for the calculation of their solvent accessibility in the native protein environment and then in the isolated domain environment. Nearly 140,000 (18%) residues showed a change in accessible surface area in the above two conditions. General features of this change under these two circumstances have been pointed out. Propensities of these interfacing amino acid residues have been calculated and their variation for different secondary structure types has been analyzed. Actual amount of surface area lost by different secondary structures is higher in the case of helix and strands compared to coil and other conformations. Overall change in surface area in hydrophobic and uncharged residues is higher than that in charged residues. An attempt has been made to know the predictability of interface residues from sequence environments. This analysis and prediction results have significant implications towards determining interacting residues in proteins and for the prediction of protein-protein, protein-ligand, protein-DNA and similar interactions.     

Specific interactions of serpins in their native forms attenuate their conformational transitions

Plasminogen activator inhibitor-1 (PAI-1) belongs to the serine protease inhibitor (serpin) protein superfamily. Serpins are unique in that their native forms are not the most thermodynamically stable conformation; instead, a more stable, latent conformation exists. During the transition to the latent form, the first strand of beta-sheet C (s1C) in the serpin is peeled away from the beta-sheet, and the reactive center loop (RCL) is inserted into beta-sheet A, rendering the serpin inactive. To elucidate the contribution of specific interactions in the metastable native form to the latency transition, we examined the effect of mutations at the s1C of PAI-1, specifically in positions P4' through P10'. Several mutations strengthened the interactions between these residues and the core protein, and slowed the transition of the protein from the metastable native form to the latent form. In particular, anchoring of the strand to the protein's hydrophobic core at the beginning (P4' site) and center of the strand (P8' site) greatly retarded the latency transition. Mutations that weakened the interactions at the s1C region facilitated the conformational conversion of the protein to the latent form. PAI-1's overall structural stability was largely unchanged by the mutations, as evaluated by urea-induced equilibrium unfolding monitored via fluorescence emission. Therefore, the mutations likely exerted their effects by modulating the height of the energy barrier from the native to the latent form. Our results show that interactions found only in the metastable native form of serpins are important structural features that attenuate folding of the proteins into their latent forms.     

Protein stability and surface electrostatics: a charged relationship

Engineering proteins to withstand a broad range of conditions continues to be a coveted objective, holding the potential to advance biomedicine, industry, and our understanding of disease. One way of achieving this goal lies in elucidating the underlying interactions that define protein stability. It has been shown that the hydrophobic effect, hydrogen bonding, and packing interactions between residues in the protein interior are dominant factors that define protein stability. The role of surface residues in protein stability has received much less attention. It has been believed that surface residues are not important for protein stability particularly because their interactions with the solvent should be similar in the native and unfolded states. In the case of surface charged residues, it was sometimes argued that solvent exposure meant that the high dielectric of the solvent will further decrease the strength of the charge-charge interactions. In this paper, we challenge the notion that the surface charged residues are not important for protein stability. We computationally redesigned sequences of five different proteins to optimize the surface charge-charge interactions. All redesigned proteins exhibited a significant increase in stability relative to their parent proteins, as experimentally determined by circular dichroism spectroscopy and differential scanning calorimetry. These results suggest that surface charge-charge interactions are important for protein stability and that rational optimization of charge-charge interactions on the protein surface can be a viable strategy for enhancing protein stability.     

Pressure versus heat-induced unfolding of ribonuclease A: the case of hydrophobic interactions within a chain-folding initiation site

To investigate the characteristics of the postulated carboxy terminal chain-folding initiation site in bovine pancreatic ribonuclease A (RNase A) (residues 106-118), important in the early stages of the folding pathway, we have engineered by site-directed mutagenesis a set of 14 predominantly conservative hydrophobic variants of the protein. The stability of each variant has been compared by pressure and temperature-induced unfolding, monitored by fourth derivative UV absorbance spectroscopy. Apparently simple two-state, reversible unfolding transitions are observed, suggesting that the disruption of tertiary structure of each protein at high pressure or temperature is strongly cooperative. Within the limits of the technique, we are unable to detect significant differences between the two processes of denaturation. Both steady-state kinetic parameters for the enzyme reaction and UV CD spectra of each RNase A variant indicate that truncation of hydrophobic side chains in this region has, in general, little or no effect on the native structure and function of the enzyme. Furthermore, the decreases in free energy of unfolding upon pressure and thermal denaturation of all the variants, particularly those modified at residues 106 and 108, suggest that the hydrophobic residues and side chain packing interactions of this region play an important role in maintaining the conformational stability of RNase A. We also demonstrate the potential of Tyr115 replacement by Trp as a non-destabilizing fluorescence probe of conformational changes local to the region.     

The folding of an enzyme. VI. The folding pathway of barnase: comparison with theoretical models.

The sequence of events in the refolding pathway of barnase has been analysed to search for general principles in protein folding. There appears to be a correlation between burying hydrophobic surface area and early folding events. All the regions that fold early interact extensively with the beta-sheet. These interactions involve predominantly hydrophobic interactions and the burial of very extensive hydrophobic areas in which multiple, close, hydrophobic-hydrophobic contacts are established around a central group of aliphatic residues. There is no burial of hydrophilic residues in these regions; those that are partly screened from the solvent make hydrogen bonds. All the regions or interactions that are made late in the folding pathway do not make extensive contacts with the beta-sheet. Their buried hydrophobic regions lack a central hydrophobic residue or residues around which other hydrophobic residues pack. Further, in some of these regions there is an extensive burial of hydrophilic residues. The results are consistent with one of the earlier events in protein folding being the local formation of native-like secondary structure elements driven by local hydrophobic surface burial. A possible candidate for an initiation site is a beta-hairpin between beta-strands 3 and 4 that is conserved in the microbial ribonuclease family. A comparison of structures in this family shows that those regions that can be superimposed, or have sequence homology, correspond to elements of structure that are formed and interact with each other early in the folding pathway, suggesting that some of these residues could be involved in directing the folding process. The data on barnase combined with results from other laboratories suggest the following tentative conclusions for the refolding of small monomeric proteins. (1) The refolding pathway is, at least in part, sequential and of compulsory order. (2) Secondary structure formation is driven by local hydrophobic surface burial and precedes the formation of most tertiary interactions. These elements are then stabilized and sometimes elongated by tertiary interactions. It is plausible that there are stop signals encoded in the linear sequence that prevent the elongation of isolated secondary structure elements in solution to a larger extent than is found in the folded protein. (3) Many tertiary interactions are not very constrained in the intermediate but become more and more defined as the hydrophobic cores consolidate, loop structures form and the configuration of surface residues takes place. The interactions between different elements of secondary structure are the last ones to be consolidated while the interactions within the secondary structure elements are consolidated earlier.(ABSTRACT TRUNCATED AT 400 WORDS)     

Oxidation of tryptophans in an interhelical hydrophobic cluster of myoglobin alters the thermodynamics of the denaturation transition

Model folding studies of sperm whale myoglobin have illustrated the presence of hydrophobic interfacial regions between elements of secondary structure. The specific oxidation of two tryptophan residues, in the A-H helix contact of sperm whale myoglobin, to the less hydrophobic oxindolylalanine residues is utilized to probe the contribution of hydrophobic packing density in this contact region. The acid denaturation of the modified protein is no longer a simple two-state process exhibiting the presence of stable intermediates. The relative stability of the intermediate is shown to be +5.3 kcal/mol less stable than native myoglobin. This value is consistent with the predicted relative stability, based upon electrostatic model calculations, of the docking of the A helix with a des-A helix myoglobin. The presence of stable intermediate structures in the denaturation pathway of the modified protein is consistent with the proposed role of hydrophobic interactions in damping structural fluctuations and statistical mechanical models of noncooperative protein unfolding. These results demonstrate the relationship between large-scale fluctuations and the frictional forces governing small-scale motions within the protein core.     

Identification of amino acids involved in protein structural uniqueness: implication for de novo protein design

Structural uniqueness is characteristic of native proteins and is essential to express their biological functions. The major factors that bring about the uniqueness are specific interactions between hydrophobic residues and their unique packing in the protein core. To find the origin of the uniqueness in their amino acid sequences, we analyzed the distribution of the side chain rotational isomers (rotamers) of hydrophobic amino acids in protein tertiary structures and derived deltaS(contact), the conformational-entropy changes of side chains by residue-residue contacts in each secondary structure. The deltaS(contact) values indicate distinct tendencies of the residue pairs to restrict side chain conformation by inter-residue contacts. Of the hydrophobic residues in alpha-helices, aliphatic residues (Leu, Val, Ile) strongly restrict the side chain conformations of each other. In beta-sheets, Met is most strongly restricted by contact with Ile, whereas Leu, Val and Ile are less affected by other residues in contact than those in alpha-helices. In designed and native protein variants, deltaS(contact) was found to correlate with the folding-unfolding cooperativity. Thus, it can be used as a specificity parameter for designing artificial proteins with a unique structure.