Determination of intrinsic hydrophilicity/hydrophobicity of amino acid side chains in peptides in the absence of nearest-neighbor or conformational effects

Understanding the hydrophilicity/hydrophobicity of amino acid side chains in peptides/proteins is one the most important aspects of biology. Though many hydrophilicity/hydrophobicity scales have been generated, an "intrinsic" scale has yet to be achieved. "Intrinsic" implies the maximum possible hydrophilicity/hydrophobicity of side chains in the absence of nearest-neighbor or conformational effects that would decrease the full expression of the side-chain hydrophilicity/hydrophobicity when the side chain is in a polypeptide chain. Such a scale is the fundamental starting point for determining the parameters that affect side-chain hydrophobicity and for quantifying such effects in peptides and proteins. A 10-residue peptide sequence, Ac-X-G-A-K-G-A-G-V-G-L-amide, was designed to enable the determination of the intrinsic values, where position X was substituted by all 20 naturally occurring amino acids and norvaline, norleucine, and ornithine. The coefficients were determined by reversed-phase high-performance liquid chromatography using six different mobile phase conditions involving different pH values (2, 5, and 7), ion-pairing reagents, and the presence and absence of different salts. The results show that the intrinsic hydrophilicity/hydrophobicity of amino acid side chains in peptides (proteins) is independent of pH, buffer conditions, or whether C(8) or C(18) reversed-phase columns were used for 17 side chains (Gly, Ala, Cys, Pro, Val, nVal, Leu, nLeu, Ile, Met, Tyr, Phe, Trp, Ser, Thr, Asn, and Gln) and dependent on pH and buffer conditions, including the type of salt or ion-pairing reagent for potentially charged side chains (Orn, Lys, His, Arg, Asp, and Glu).     

Hydrophobicity regained.

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

Side chain burial and hydrophobic core packing in protein folding transition states

A critical step in the folding pathway of globular proteins is the formation of a tightly packed hydrophobic core. Several mutational studies have addressed the question of whether tight packing interactions are present during the rate-limiting step of folding. In some of these investigations, substituted side chains have been assumed to form native-like interactions in the transition state when the folding rates of mutant proteins correlate with their native-state stabilities. Alternatively, it has been argued that side chains participate in nonspecific hydrophobic collapse when the folding rates of mutant proteins correlate with side-chain hydrophobicity. In a reanalysis of published data, we have found that folding rates often correlate similarly well, or poorly, with both native-state stability and side-chain hydrophobicity, and it is therefore not possible to select an appropriate transition state model based on these one-parameter correlations. We show that this ambiguity can be resolved using a two-parameter model in which side chain burial and the formation of all other native-like interactions can occur asynchronously. Notably, the model agrees well with experimental data, even for positions where the one-parameter correlations are poor. We find that many side chains experience a previously unrecognized type of transition state environment in which specific, native-like interactions are formed, but hydrophobic burial dominates. Implications of these results to the design and analysis of protein folding studies are discussed.     

Contribution of amino acid substitutions at two different interior positions to the conformational stability of human lysozyme.

To elucidate correlative relationships between structural change and thermodynamic stability in proteins, a series of mutant human lysozymes modified at two buried positions (Ile56 and Ile59) were examined. Their thermodynamic parameters of denaturation and crystal structures were studied by calorimetry and X-ray crystallography. The mutants at positions 56 and 59 exhibited different responses to a series of amino acid substitutions. The changes in stability due to substitutions showed a linear correlation with changes in hydrophobicity of substituted residues, having different slopes at each mutation site. However, the stability of each mutant was found to be represented by a unique equation involving physical properties calculated from mutant structures. By fitting present and previous stability data for mutant human lysozymes substituted at various positions to the equation, the magnitudes of the hydrophobicity of a carbon atom and the hydrophobicity of nitrogen and neutral oxygen atoms were found to be 0.178 and -0.013 kJ/mol.A(2), respectively. It was also found that the contribution of a hydrogen bond with a length of 3.0 A to protein stability was 5.1 kJ/mol and the entropy loss of newly introduction of a water molecules was 7.8 kJ/mol.     

The role of position a in determining the stability and oligomerization state of alpha-helical coiled coils: 20 amino acid stability coefficients in the hydrophobic core of proteins

We describe here a systematic investigation into the role of position a in the hydrophobic core of a model coiled-coil protein in determining coiled-coil stability and oligomerization state. We employed a model coiled coil that allowed the formation of an extended three-stranded trimeric oligomerization state for some of the analogs; however, due to the presence of a Cys-Gly-Gly linker, unfolding occurred from the same two-stranded monomeric oligomerization state for all of the analogs. Denaturation from a two-stranded state allowed us to measure the relative contribution of 20 different amino acid side chains to coiled-coil stability from chemical denaturation profiles. In addition, the relative hydrophobicity of the substituted amino acid side chains was assessed by reversed-phase high-performance liquid chromatography and found to correlate very highly (R = 0.95) with coiled-coil stability. We also determined the effect of position a in specifying the oligomerization state using ultracentrifugation as well as high-performance size-exclusion chromatography. We found that nine of the analogs populated one oligomerization state exclusively at peptide concentrations of 50 microM under benign buffer conditions. The Leu-, Tyr-, Gln-, and His-substituted analogs were found to be exclusively three-stranded trimers, while the Asn-, Lys-, Orn-, Arg-, and Trp-substituted analogs formed exclusively two-stranded monomers. Modeling results for the Leu-substituted analog showed that a three-stranded oligomerization state is preferred due to increased side-chain burial, while a two-stranded oligomerization state was observed for the Trp analog due to unfavorable cavity formation in the three-stranded state.     

Solvent effects on protein association and protein folding

Solvent effects on the thermodynamics of two processes--folding of proteins and association between proteins--are examined in detail. A complete inventory of the multitude of solvent effects may be obtained by employing the concept of conditional solvation free energy. This theoretical tool allows for the isolation of specific side-chain effects from the entire protein and for the study of its contribution to the overall free energy change in small model compounds. Some numerical examples are presented, and ways of estimating other cases, for which no relevant experimental data are available, are suggested. Our findings lead to the conclusion that the currently used hydrophobicity scales, based on partition coefficients between water and an organic solvent, are inadequate measures of the contribution of side chains being transferred from water to the interior of the protein. We have also tentatively concluded that correlation between hydrophilic functional groups might be more important than correlations between hydrophobic side chains.     

Contribution of intra- and intermolecular hydrogen bonds to the conformational stability of human lysozyme(,).

In globular proteins, there are intermolecular hydrogen bonds between protein and water molecules, and between water molecules, which are bound with the proteins, in addition to intramolecular hydrogen bonds. To estimate the contribution of these hydrogen bonds to the conformational stability of a protein, the thermodynamic parameters for denaturation and the crystal structures of five Thr to Val and five Thr to Ala mutant human lysozymes were determined. The denaturation Gibbs energy (DeltaG) of Thr to Val and Thr to Ala mutant proteins was changed from 4.0 to -5.6 kJ/mol and from 1.6 to -6.3 kJ/mol, respectively, compared with that of the wild-type protein. The contribution of hydrogen bonds to the stability (DeltaDeltaG(HB)) of the Thr and other mutant human lysozymes previously reported was extracted from the observed stability changes (DeltaDeltaG) with correction for changes in hydrophobicity and side chain conformational entropy between the wild-type and mutant structures. The estimation of the DeltaDeltaG(HB) values of all mutant proteins after removal of hydrogen bonds, including protein-water hydrogen bonds, indicates a favorable contribution of the intra- and intermolecular hydrogen bonds to the protein stability. The net contribution of an intramolecular hydrogen bond (DeltaG(HB[pp])), an intermolecular one between protein and ordered water molecules (DeltaG(HB[pw])), and an intermolecular one between ordered water molecules (DeltaG(HB[ww])) could be estimated to be 8. 5, 5.2, and 5.0 kJ/mol, respectively, for a 3 A long hydrogen bond. This result shows the different contributions to protein stability of intra- and intermolecular hydrogen bonds. The entropic cost due to the introduction of a water molecule (DeltaG(H)()2(O)) could be also estimated to be about 8 kJ/mol.     

Energetics, stability, and prediction of transmembrane helices

We show that the peptide backbone of an alpha-helix places a severe thermodynamic constraint on transmembrane (TM) stability. Neglect of this constraint by commonly used hydrophobicity scales underlies the notorious uncertainty of TM helix prediction by sliding-window hydropathy plots of membrane protein (MP) amino acid sequences. We find that an experiment-based whole-residue hydropathy scale (WW scale), which includes the backbone constraint, identifies TM helices of membrane proteins with an accuracy greater than 99 %. Furthermore, it correctly predicts the minimum hydrophobicity required for stable single-helix TM insertion observed in Escherichia coli. In order to improve membrane protein topology prediction further, we introduce the augmented WW (aWW) scale, which accounts for the energetics of salt-bridge formation. An important issue for genomic analysis is the ability of the hydropathy plot method to distinguish membrane from soluble proteins. We find that the method falsely predicts 17 to 43 % of a set of soluble proteins to be MPs, depending upon the hydropathy scale used.     

How do side chains orient globally in protein structures?

An angle Omega is defined to serve as a metric for global side-chain orientations, which reflects the orientation of the side chain relative to the radial vector from the center of the protein to an amino acid. The side-chain orientations of buried residues exhibit characteristically different orientations than do exposed residues, in both monomeric and dimeric structures. Overall, buried side chains point mostly inward, whereas surface side chains tend to point outward from the surface. This difference in behavior also correlates well with the residue hydrophobicity; so a global side-chain orientation can be viewed as a direct structural manifestation of hydrophobicity. When various solvent-accessible layers are considered, the behavior is relatively continuous between centrally located and exposed residues. In the case of interfacial residues between subunits, there are statistically significant differences between exposed residues and interface residues for ALA, ARG, ASN, ASP, GLU, HIS, LYS, THR, VAL, MET, PRO, and overall the interface residues have an increased tendency to point inward. Presumably, these substantial differences in orientations of side chains may be a manifestation of hydrophobic forces.     

Contribution of polar groups in the interior of a protein to the conformational stability

It has been generally believed that polar residues are usually located on the surface of protein structures. However, there are many polar groups in the interior of the structures in reality. To evaluate the contribution of such buried polar groups to the conformational stability of a protein, nonpolar to polar mutations (L8T, A9S, A32S, I56T, I59T, I59S, A92S, V93T, A96S, V99T, and V100T) in the interior of a human lysozyme were examined. The thermodynamic parameters for denaturation were determined using a differential scanning calorimeter, and the crystal structures were analyzed by X-ray crystallography. If a polar group had a heavy energy cost to be buried, a mutant protein would be remarkably destabilized. However, the stability (Delta G) of the Ala to Ser and Val to Thr mutant human lysozymes was comparable to that of the wild-type protein, suggesting a low-energy penalty of buried polar groups. The structural analysis showed that all polar side chains introduced in the mutant proteins were able to find their hydrogen bond partners, which are ubiquitous in protein structures. The empirical structure-based calculation of stability change (Delta Delta G) [Takano et al. (1999) Biochemistry 38, 12698--12708] revealed that the mutant proteins decreased the hydrophobic effect contributing to the stability (Delta G(HP)), but this destabilization was recovered by the hydrogen bonds newly introduced. The present study shows the favorable contribution of polar groups with hydrogen bonds in the interior of protein molecules to the conformational stability.     

Monitoring aromatic picosecond to nanosecond dynamics in proteins via 13C relaxation: expanding perturbation mapping of the rigidifying core mutation, V54A, in eglin c

Long-range effects, such as allostery, have evolved in proteins as a means of regulating function via communication between distal sites. An NMR-based perturbation mapping approach was used to more completely probe the dynamic response of the core mutation V54A in the protein eglin c by monitoring changes in picosecond to nanosecond aromatic side-chain dynamics and H/D exchange stabilities. Previous side-chain dynamics studies on this mutant were limited to methyl-bearing residues, most of which were found to rigidify on the picosecond to nanosecond time scale in the form of a contiguous "network". Here, high precision (13)C relaxation data from 13 aromatic side chains were acquired by applying canonical relaxation experiments to a newly developed carbon labeling scheme [Teilum et al. (2006) J. Am. Chem. Soc. 128, 2506-2507]. The fitting of model-free parameters yielded S (2) variability which is intermediate with respect to backbone and methyl-bearing side-chain variability and tau e values that are approximately 1 ns. Inclusion of the aromatic dynamic response results in an expanded network of dynamically coupled residues, with some aromatics showing increases in flexibility, which partially offsets the rigidification in methyl side chains. Using amide hydrogen exchange, dynamic propagation on a slower time scale was probed in response to the V54A perturbation. Surprisingly, regional stabilization (slowed exchange) 10-12 A from the site of mutation was observed despite a global destabilization of 1.5 kcal x mol (-1). Furthermore, this unlikely pocket of stabilized residues colocalizes with increases in aromatic flexibility on the faster time scale. Because the converse is also true (destabilized residues colocalize with rigidification on the fast time scale), a plausible entropy-driven mechanism is discussed for relating colocalization of opposing dynamic trends on vastly different time scales.     

Side-chain flexibility in proteins upon ligand binding

Ligand binding may involve a wide range of structural changes in the receptor protein, from hinge movement of entire domains to small side-chain rearrangements in the binding pocket residues. The analysis of side chain flexibility gives insights valuable to improve docking algorithms and can provide an index of amino-acid side-chain flexibility potentially useful in molecular biology and protein engineering studies. In this study we analyzed side-chain rearrangements upon ligand binding. We constructed two non-redundant databases (980 and 353 entries) of "paired" protein structures in complexed (holo-protein) and uncomplexed (apo-protein) forms from the PDB macromolecular structural database. The number and identity of binding pocket residues that undergo side-chain conformational changes were determined. We show that, in general, only a small number of residues in the pocket undergo such changes (e.g., approximately 85% of cases show changes in three residues or less). The flexibility scale has the following order: Lys > Arg, Gln, Met > Glu, Ile, Leu > Asn, Thr, Val, Tyr, Ser, His, Asp > Cys, Trp, Phe; thus, Lys side chains in binding pockets flex 25 times more often then do the Phe side chains. Normalizing for the number of flexible dihedral bonds in each amino acid attenuates the scale somewhat, however, the clear trend of large, polar amino acids being more flexible in the pocket than aromatic ones remains. We found no correlation between backbone movement of a residue upon ligand binding and the flexibility of its side chain. These results are relevant to 1. Reduction of search space in docking algorithms by inclusion of side-chain flexibility for a limited number of binding pocket residues; and 2. Utilization of the amino acid flexibility scale in protein engineering studies to alter the flexibility of binding pockets.     

Dynamics of a de novo designed three-helix bundle protein studied by 15N, 13C, and 2H NMR relaxation methods

Understanding how the amino acid sequence of a polypeptide chain specifies a unique, functional three-dimensional structure remains an important goal, especially in the context of the emerging discipline of de novo protein design. Alpha3D is a single chain protein of 73 amino acids resulting from a de novo design effort. Previous solution nuclear magnetic resonance studies of alpha3D confirm that the protein adopts the designed structure of a three-helix bundle. Furthermore, alpha3D has been previously shown to possess all of the major thermodynamic and structural characteristics of natural proteins, though it shares no sequence homology to any protein sequence in the database. In this work, the backbone and side-chain dynamics of alpha3D were investigated using 15N, 13C, and 2H nuclear magnetic resonance relaxation methods with the aim of assessing the character of the internal motions of this native-like protein of de novo design. At the backbone level, both 15N and 13C(alpha) relaxation studies indicate highly restrictive motion on the picosecond to nanosecond time scale in the alpha-helical regions of alpha3D, with increasing mobility at the ends of the alpha-helices and in the two loop regions. This is largely consistent with what is seen in proteins of natural origin. Overall, the view provided by both 2H and 13C methyl relaxation methods suggest that the side chains of alpha3D are more dynamic compared to natural proteins. Regions of relative flexibility bound clusters of rigid methyl-bearing side-chain groups that are interspersed with aromatic and beta-branched amino acids. The time scale of motions associated with methyl-bearing side chains of alpha3D are significantly longer than that seen in natural proteins. These results indicate that the strategies underlying the design of alpha3D have largely, but not completely, captured both the structural and dynamic character of natural proteins.     

Electronic properties of the amino acid side chains contribute to the structural preferences in protein folding.

A database of 118 non-redundant proteins was examined to determine the preferences of amino acids for secondary structures: alpha-helix, beta-strand and coil conformations. To better understand how the physicochemical properties of amino acid side chains might influence protein folding, several new scales have been suggested for quantifying the electronic effects of amino acids. These include the pKa at the amino group, localized effect substituent constants (esigma), and a composite of these two scales (epsilon). Amino acids were also classified into 5 categories on the basis of their electronic properties: O (strong electron donor), U (weak donor), Z (ambivalent), B (weak electron acceptor), and X (strong acceptor). Certain categories of amino acid appeared to be critical for particular conformations, e.g., O and U-type residues for alpha-helix formation. Pairwise analysis of the database according to these categories revealed significant context effects in the structural preferences. In general, the propensity of an amino acid for a particular conformation was related to the electronic features of the side chain. Linear regression analyses revealed that the electronic properties of amino acids contributed about as much to the folding preferences as hydrophobicity, which is a well-established determinant of protein folding. A theoretical model has been proposed to explain how the electronic properties of the side chain groups might influence folding along the peptide backbone.     

Protein stability and electrostatic interactions between solvent exposed charged side chains

To investigate the contribution to protein stability of electrostatic interactions between charged surface residues, we have studied the effect of substituting three negatively charged solvent exposed residues with their side-chain amide analogs in bovine calbindin D9k--a small (Mr 8,500) globular protein of the calmodulin superfamily. The free energy of urea-induced unfolding for the wild-type and seven mutant proteins has been measured. The mutant proteins have increased stability towards unfolding relative to the wild-type. The experimental results correlate reasonably well with theoretically calculated relative free energies of unfolding and show that electrostatic interactions between charges on the surface of a protein can have significant effects on protein stability.     

Electronic Properties of the Amino Acid Side Chains Contribute to the Structural Preferences in Protein Folding

Abstract

A database of 118 non-redundant proteins was examined to determine the preferences of amino acids for secondary structures: α-helix, β-strand and coil conformations. To better understand how the physicochemical properties of amino acid side chains might influence protein folding, several new scales have been suggested for quantifying the electronic effects of amino acids. These include the pKa at the amino group, localized effect substituent constants (eσ), and a composite of these two scales (ε). Amino acids were also classified into 5 categories on the basis of their electronic properties: O (strong electron donor), U (weak donor), Z (ambivalent), B (weak electron acceptor), and X (strong acceptor). Certain categories of amino acid appeared to be critical for particular conformations, e.g., O and U-type residues for α-helix formation. Pairwise analysis of the database according to these categories revealed significant context effects in the structural preferences. In general, the propensity of an amino acid for a particular conformation was related to the electronic features of the side chain. Linear regression analyses revealed that the electronic properties of amino acids contributed about as much to the folding preferences as hydrophobicity, which is a well-established determinant of protein folding. A theoretical model has been proposed to explain how the electronic properties of the side chain groups might influence folding along the peptide backbone.     

Solvent-amino acid interaction energies in three-dimensional-lattice Monte Carlo simulations of a model 27-mer protein: Folding thermodynamics and kinetics

Amino acid residue-solvent interactions are required for lattice Monte Carlo simulations of model proteins in water. In this study, we propose an interaction-energy scale that is based on the interaction scale by Miyazawa and Jernigan. It permits systematic variation of the amino acid-solvent interactions by introducing a contrast parameter for the hydrophobicity, C(s), and a mean attraction parameter for the amino acids, omega. Changes in the interaction energies strongly affect many protein properties. We present an optimized energy parameter set for best representing realistic behavior typical for many proteins (fast folding and high cooperativity for single chains). Our optimal parameters feature a much weaker hydrophobicity contrast and mean attraction than does the original interaction scale. The proposed interaction scale is designed for calculating the behavior of proteins in bulk and at interfaces as a function of solvent characteristics, as well as protein size and sequence.     

Role of amino acid residues at turns in the conformational stability and folding of human lysozyme

To clarify the role of amino acid residues at turns in the conformational stability and folding of a globular protein, six mutant human lysozymes deleted or substituted at turn structures were investigated by calorimetry, GuHCl denaturation experiments, and X-ray crystal analysis. The thermodynamic properties of the mutant and wild-type human lysozymes were compared and discussed on the basis of their three-dimensional structures. For the deletion mutants, Delta47-48 and Delta101, the deleted residues are in turns on the surface and are absent in human alpha-lactalbumin, which is homologous to human lysozyme in amino acid sequence and tertiary structure. The stability of both mutants would be expected to increase due to a decrease in conformational entropy in the denatured state; however, both proteins were destabilized. The destabilizations were mainly caused by the disappearance of intramolecular hydrogen bonds. Each part deleted was recovered by the turn region like the alpha-lactalbumin structure, but there were differences in the main-chain conformation of the turn between each deletion mutant and alpha-lactalbumin even if the loop length was the same. For the point mutants, R50G, Q58G, H78G, and G37Q, the main-chain conformations of these substitution residues located in turns adopt a left-handed helical region in the wild-type structure. It is thought that the left-handed non-Gly residue has unfavorable conformational energy compared to the left-handed Gly residue. Q58G was stabilized, but the others had little effect on the stability. The structural analysis revealed that the turns could rearrange the main-chain conformation to accommodate the left-handed non-Gly residues. The present results indicate that turn structures are able to change their main-chain conformations, depending upon the side-chain features of amino acid residues on the turns. Furthermore, stopped-flow GuHCl denaturation experiments on the six mutants were performed. The effects of mutations on unfolding-refolding kinetics were significantly different among the mutant proteins. The deletion/substitutions in turns located in the alpha-domain of human lysozyme affected the refolding rate, indicating the contribution of turn structures to the folding of a globular protein.     

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

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