Hydrophobicity scales

Hydrophobicity scales are values that define the relative hydrophobicity or hydrophilicity of amino acid residues. The more positive the value, the more hydrophobic are the amino acids located in that region of the protein. These scales are commonly used to predict the transmembrane alpha-helices of membrane proteins. When consecutively measuring amino acids of a protein, changes in value indicate attraction of specific protein regions towards the hydrophobic region inside lipid bilayer.

The hydrophobic or hydrophilic character of a compound or amino acid is its hydropathic character,^[1] hydropathicity, or hydropathy.

The hydrophobic effect represents the tendency of water to exclude non-polar molecules. The effect originates from the disruption of highly dynamic hydrogen bonds between molecules of liquid water. Polar chemical groups, such as OH group in methanol do not cause the hydrophobic effect. However, a pure hydrocarbon molecule, for example hexane, cannot accept or donate hydrogen bonds to water. Introduction of hexane into water causes disruption of the hydrogen bonding network between water molecules. The hydrogen bonds are partially reconstructed by building a water "cage" around the hexane molecule, similar to that in clathrate hydrates formed at lower temperatures. The mobility of water molecules in the "cage" (or solvation shell) is strongly restricted. This leads to significant losses in translational and rotational entropy of water molecules and makes the process unfavorable in terms of free energy of the system.^[2][3][4][5] In terms of thermodynamics, the hydrophobic effect is the free energy change of water surrounding a solute.^[6] A positive free energy change of the surrounding solvent indicates hydrophobicity, whereas a negative free energy change implies hydrophilicity. In this way, the hydrophobic effect not only can be localized but also decomposed into enthalpic and entropic contributions.

A number of different hydrophobicity scales have been developed.^{[3][1][7][8][9]} The Expasy Protscale website lists a total of 22 hydrophobicity scales.^[10]

There are clear differences between the four scales shown in the table.^[11] Both the second and fourth scales place cysteine as the most hydrophobic residue, unlike the other two scales. This difference is due to the different methods used to measure hydrophobicity. The method used to obtain the Janin and Rose et al. scales was to examine proteins with known 3-D structures and define the hydrophobic character as the tendency for a residue to be found inside of a protein rather than on its surface.^[12][13] Since cysteine forms disulfide bonds that must occur inside a globular structure, cysteine is ranked as the most hydrophobic. The first and third scales are derived from the physiochemical properties of the amino acid side chains. These scales result mainly from inspection of the amino acid structures.^[14][1] Biswas et al., divided the scales based on the method used to obtain the scale into five different categories.^[3]

The most common method of measuring amino acid hydrophobicity is partitioning between two immiscible liquid phases. Different organic solvents are most widely used to mimic the protein interior. However, organic solvents are slightly miscible with water and the characteristics of both phases change making it difficult to obtain pure hydrophobicity scale.^[3] Nozaki and Tanford proposed the first major hydrophobicity scale for nine amino acids.^[15] Ethanol and dioxane are used as the organic solvents and the free energy of transfer of each amino acid was calculated. Non liquid phases can also be used with partitioning methods such as micellar phases and vapor phases. Two scales have been developed using micellar phases.^[16][17] Fendler et al. measured the partitioning of 14 radiolabeled amino acids using sodium dodecyl sulfate (SDS) micelles. Also, amino acid side chain affinity for water was measured using vapor phases.^[14] Vapor phases represent the simplest non polar phases, because it has no interaction with the solute.^[18] The hydration potential and its correlation to the appearance of amino acids on the surface of proteins was studied by Wolfenden. Aqueous and polymer phases were used in the development of a novel partitioning scale.^[19] Partitioning methods have many drawbacks. First, it is difficult to mimic the protein interior.^[20][21] In addition, the role of self solvation makes using free amino acids very difficult. Moreover, hydrogen bonds that are lost in the transfer to organic solvents are not reformed but often in the interior of protein.^[22]

Hydrophobicity scales can also be obtained by calculating the solvent accessible surface areas for amino acid residues in the expended polypeptide chain^[22] or in alpha-helix and multiplying the surface areas by the empirical solvation parameters for the corresponding types of atoms.^[3] A differential solvent accessible surface area hydrophobicity scale based on proteins as compacted networks near a critical point, due to self-organization by evolution, was constructed based on asymptotic power-law (self-similar) behavior.^[23][24] This scale is based on a bioinformatic survey of 5526 high-resolution structures from the Protein Data Bank. This differential scale has two comparative advantages: (1) it is especially useful for treating changes in water-protein interactions that are too small to be accessible to conventional force-field calculations, and (2) for homologous structures, it can yield correlations with changes in properties from mutations in the amino acid sequences alone, without determining corresponding structural changes, either in vitro or in vivo.

Reversed phase liquid chromatography (RPLC) is the most important chromatographic method for measuring solute hydrophobicity.^[3][25] The non polar stationary phase mimics biological membranes. Peptide usage has many advantages because partition is not extended by the terminal charges in RPLC. Also, secondary structures formation is avoided by using short sequence peptides. Derivatization of amino acids is necessary to ease its partition into a C18 bonded phase. Another scale had been developed in 1971 and used peptide retention on hydrophilic gel.^[26] 1-butanol and pyridine were used as the mobile phase in this particular scale and glycine was used as the reference value. Pliska and his coworkers^[27] used thin layer chromatography to relate mobility values of free amino acids to their hydrophobicities. About a decade ago, another hydrophilicity scale was published, this scale used normal phase liquid chromatography and showed the retention of 121 peptides on an amide-80 column.^[28] The absolute values and relative rankings of hydrophobicity determined by chromatographic methods can be affected by a number of parameters. These parameters include the silica surface area and pore diameter, the choice and pH of aqueous buffer, temperature and the bonding density of stationary phase chains.^[3]

This method use DNA recombinant technology and it gives an actual measurement of protein stability. In his detailed site-directed mutagenesis studies, Utani and his coworkers substituted 19 amino acids at Trp49 of the tryptophan synthase and he measured the free energy of unfolding. They found that the increased stability is directly proportional to increase in hydrophobicity up to a certain size limit. The main disadvantage of site-directed mutagenesis method is that not all the 20 naturally occurring amino acids can substitute a single residue in a protein. Moreover, these methods have cost problems and is useful only for measuring protein stability.^[3][29]

The hydrophobicity scales developed by physical property methods are based on the measurement of different physical properties. Examples include, partial molar heat capacity, transition temperature and surface tension. Physical methods are easy to use and flexible in terms of solute. The most popular hydrophobicity scale was developed by measuring surface tension values for the naturally occurring 20 amino acids in NaCl solution.^[30] The main drawbacks of surface tension measurements is that the broken hydrogen bonds and the neutralized charged groups remain at the solution air interface.^[3][1] Another physical property method involve measuring the solvation free energy.^[31] The solvation free energy is estimated as a product of an accessibility of an atom to the solvent and an atomic solvation parameter. Results indicate the solvation free energy lowers by an average of 1 Kcal/residue upon folding.^[3]

Palliser and Parry have examined about 100 scales and found that they can use them for locating B-strands on the surface of proteins.^[32] Hydrophobicity scales were also used to predict the preservation of the genetic code.^[33] Trinquier observed a new order of the bases that better reflect the conserved character of the genetic code.^[3] They believed new ordering of the bases was uracil-guanine-cystosine-adenine (UGCA) better reflected the conserved character of the genetic code compared to the commonly seen ordering UCAG.^[3]

The Wimley–White whole residue hydrophobicity scales are significant for two reasons. First, they include the contributions of the peptide bonds as well as the sidechains, providing absolute values. Second, they are based on direct, experimentally determined values for transfer free energies of polypeptides.

Two whole-residue hydrophobicity scales have been measured:

• One for the transfer of unfolded chains from water to the bilayer interface (referred to as the Wimley–White interfacial hydrophobicity scale).
• One for the transfer of unfolded chains into octanol, which is relevant to the hydrocarbon core of a bilayer.

The Stephen H. White website^[34] provides an example of whole residue hydrophobicity scales showing the free energy of transfer ΔG(kcal/mol) from water to POPC interface and to n-octanol.^[34] These two scales are then used together to make Whole residue hydropathy plots.^[34] The hydropathy plot constructed using ΔG_woct − ΔG_wif shows favorable peaks on the absolute scale that correspond to the known TM helices. Thus, the whole residue hydropathy plots illustrate why transmembrane segments prefer a transmembrane location rather than a surface one.^{[35][36][37][38]}

Most of the existing hydrophobicity scales are derived from the properties of amino acids in their free forms or as a part of a short peptide. Bandyopadhyay-Mehler hydrophobicity scale was based on partitioning of amino acids in the context of protein structure. Protein structure is a complex mosaic of various dielectric medium generated by arrangement of different amino acids. Hence, different parts of the protein structure most likely would behave as solvents with different dielectric values. For simplicity, each protein structure was considered as an immiscible mixture of two solvents, protein interior and protein exterior. The local environment around individual amino acid (termed as "micro-environment") was computed for both protein interior and protein exterior. The ratio gives the relative hydrophobicity scale for individual amino acids. Computation was trained on high resolution protein crystal structures. This quantitative descriptor for microenvironment was derived from the octanol-water partition coefficient, (known as Rekker's Fragmental Constants) widely used for pharmacophores. This scale well correlate with the existing methods, based on partitioning and free energy computations. Advantage of this scale is it is more realistic, as it is in the context of real protein structures.^[9]

In the field of engineering, the hydrophobicity (or dewetting ability) of a flat surface (e.g., a counter top in kitchen or a cooking pan) can be measured by the contact angle of water droplet. A University of Nebraska-Lincoln team recently devised a computational approach that can relate the molecular hydrophobicity scale of amino-acid chains to the contact angle of water nanodroplet.^[39] The team constructed planar networks composed of unified amino-acid side chains with native structure of the beta-sheet protein. Using molecular dynamics simulation, the team is able to measure the contact angle of water nanodroplet on the planar networks (caHydrophobicity).

On the other hand, previous studies show that the minimum of excess chemical potential of a hard-sphere solute with respect to that in the bulk exhibits a linear dependence on cosine value of contact angle.^[40] Based on the computed excess chemical potentials of the purely repulsive methane-sized Weeks–Chandler–Andersen solute with respect to that in the bulk, the extrapolated values of cosine value of contact angle are calculated(ccHydrophobicity), which can be used to quantify the hydrophobicity of amino acid side chains with complete wetting behaviors.

• Hydrophobic mismatch

• ProtScale (web-based tool for calculating hydropathy plots)
• NetSurfP - Secondary Structure and Surface accessibility predictor
• Whole residue hydrophobicity scale
• Membrane protein explorer