A new method for identification of protein (sub)families in a set of proteins based on hydropathy distribution in proteins

Structural similarity among proteins is reflected in the distribution of hydropathicity along the amino acids in the protein sequence. Similarities in the hydropathy distributions are obvious for homologous proteins within a protein family. They also were observed for proteins with related structures, even when sequence similarities were undetectable. Here we present a novel method that employs the hydropathy distribution in proteins for identification of (sub)families in a set of (homologous) proteins. We represent proteins as points in a generalized hydropathy space, represented by vectors of specifically defined features. The features are derived from hydropathy of the individual amino acids. Projection of this space onto principal axes reveals groups of proteins with related hydropathy distributions. The groups identified correspond well to families of structurally and functionally related proteins. We found that this method accurately identifies protein families in a set of proteins, or subfamilies in a set of homologous proteins. Our results show that protein families can be identified by the analysis of hydropathy distribution, without the need for sequence alignment.     

A two-component nucleation model of protein hydrophobicity

A fundamental characteristic of soluble globular protein structure is a hydrophobic core and protein exterior comprised predominantly of hydrophilic residues. This distribution of amino acid residue hydrophobicity, from protein interior to exterior, has recently been profiled with the use of hydrophobic moments. The calculations enable comparison of the radial hydrophobicity distribution of different proteins and had revealed two features common to 30 proteins of diverse size and structure. One, a global feature, is the overall shape of the second-order ellipsoidal hydrophobic moment. The second, a specific feature, is a quasi-invariant hydrophobic-ratio of distances. Both features are dependent upon the rates of increase, from protein interior to exterior, of the accumulated numbers of hydrophobic and hydrophilic amino acid residues. These rates can be simulated simply with a two-component nucleation model of protein hydrophobicity. The model provides insight into the origin of the shape of the observed hydrophobic moment profiles and of the observed range of hydrophobic ratios. Consistent with observation, it is shown that a relatively wide range of hydrophobic and hydrophilic rates of increase yield a relatively narrow range of hydrophobic ratios. Furthermore, the model identifies one factor, the decrease in residue density with increasing distance from the protein interior, that is critical in providing the range of values that is comparable with the observed range.     

Intrinsic disorder in protein sense‐antisense recognition

In addition to the well‐established sense‐antisense complementarity abundantly present in the nucleic acid world and serving as a basic principle of the specific double‐helical structure of DNA, production of mRNA, and genetic code‐based biosynthesis of proteins, sense‐antisense complementarity is also present in proteins, where sense and antisense peptides were shown to interact with each other with increased probability. In nucleic acids, sense‐antisense complementarity is achieved via the Watson‐Crick complementarity of the base pairs or nucleotide pairing. In proteins, the complementarity between sense and antisense peptides depends on a specific hydropathic pattern, where codons for hydrophilic and hydrophobic amino acids in a sense peptide are complemented by the codons for hydrophobic and hydrophilic amino acids in its antisense counterpart. We are showing here that in addition to this pattern of the complementary hydrophobicity, sense and antisense peptides are characterized by the complementary order‐disorder patterns and show complementarity in sequence distribution of their disorder‐based interaction sites. We also discuss how this order‐disorder complementarity can be related to protein evolution.     

A simple method for displaying the hydropathic character of a protein

A computer program that progressively evaluates the hydrophilicity and hydrophobicity of a protein along its amino acid sequence has been devised. For this purpose, a hydropathy scale has been composed wherein the hydrophilic and hydrophobic properties of each of the 20 amino acid side-chains is taken into consideration. The scale is based on an amalgam of experimental observations derived from the literature. The program uses a moving-segment approach that continuously determines the average hydropathy within a segment of predetermined length as it advances through the sequence. The consecutive scores are plotted from the amino to the carboxy terminus. At the same time, a midpoint line is printed that corresponds to the grand average of the hydropathy of the amino acid compositions found in most of the sequenced proteins. In the case of soluble, globular proteins there is a remarkable correspondence between the interior portions of their sequence and the regions appearing on the hydrophobic side of the midpoint line, as well as the exterior portions and the regions on the hydrophilic side. The correlation was demonstrated by comparisons between the plotted values and known structures determined by crystallography. In the case of membrane-bound proteins, the portions of their sequences that are located within the lipid bilayer are also clearly delineated by large uninterrupted areas on the hydrophobic side of the midpoint line. As such, the membrane-spanning segments of these proteins can be identified by this procedure. Although the method is not unique and embodies principles that have long been appreciated, its simplicity and its graphic nature make it a very useful tool for the evaluation of protein structures.     

The jigsaw puzzle model: search for conformational specificity in protein interiors

The jigsaw puzzle model postulates that the predominant factor relating primary sequence to three-dimensional fold lies in the stereospecific packing of interdigitating side-chains within densely packed protein interiors. An attempt has been made to check the validity of the model by means of a surface complementarity function. Out of a database of 100 highly resolved protein structures the contacts between buried hydrophobic residues (Leu, Ile, Val, Phe) and their neighbours have been categorized in terms of the extent of side-chain surface area involved in a contact (overlap) and their steric fit (Sm). The results show that the majority of contacts between a buried residue and its immediate neighbours (side-chains) are of high steric fit and in the case of extended overlap at least one of the angular parameters characterizing interresidue geometry to have pronounced deviation from a random distribution, estimated by chi(2). The calculations thus tend to support the "jigsaw puzzle" model in that 75-85% of the contacts involving hydrophobic residues are of high surface complementarity, which, coupled to high overlap, exercise fairly stringent constraints over the possible geometrical orientations between interacting residues. These constraints manifest in simple patterns in the distributions of orientational angles. Approximately 60-80% of the buried side-chain surface packs against neighbouring side-chains, the rest interacting with main-chain atoms. The latter partition of the surface maintains an equally high steric fit (relative to side-chain contacts) emphasizing a non-trivial though secondary role played by main-chain atoms in interior packing. The majority of this class of contacts, though of high complementarity, is of reduced overlap. All residues whether hydrophobic or polar/charged show similar surface complementarity measures upon burial, indicating comparable competence of all amino acids in packing effectively with their atomic environments. The specificity thus appears to be distributed over the entire network of contacts within proteins. The study concludes with a proposal to classify contacts as specific and non-specific (based on overlap and fit), with the former perhaps contributing more to the specificity between sequence and fold than the latter.     

Hydrophobic complementarity in protein-protein docking

Formation of hydrophobic contacts across a newly formed interface is energetically favorable. Based on this observation we developed a geometric-hydrophobic docking algorithm that estimates quantitatively the hydrophobic complementarity at protein-protein interfaces. Each molecule to be docked is represented as a grid of complex numbers, storing information regarding the shape of the molecule in the real part and information regarding the hydropathy of the surface in the imaginary part. The grid representations are correlated using fast Fourier transformations. The algorithm is used to compare the extent of hydrophobic complementarity in oligomers (represented by D2 tetramers) and in hetero-dimers of soluble proteins (complexes). We also test the implication of hydrophobic complementarity in distinguishing correct from false docking solutions. We find that hydrophobic complementarity at the interface exists in oligomers and in complexes, and in both groups the extent of such complementarity depends on the size of the interface. Thus, the non-polar portions of large interfaces are more often juxtaposed than non-polar portions of small interfaces. Next we find that hydrophobic complementarity helps to point out correct docking solutions. In oligomers it significantly improves the ranks of nearly correct reassembled and modeled tetramers. Combining geometric, electrostatic and hydrophobic complementarity for complexes gives excellent results, ranking a nearly correct solution < 10 for 5 of 23 tested systems, < 100 for 8 systems and < 1000 for 19 systems.     

A quantitative measure of error minimization in the genetic code

Summary We have calculated the average effect of changing a codon by a single base for all possible single-base changes in the genetic code and for changes in the first, second, and third codon positions separately. Such values were calculated for an amino acid's polar requirement, hydropathy, molecular volume, and isoelectric point. For each attribute the average effect of single-base changes was also calculated for a large number of randomly generated codes that retained the same level of redundancy as the natural code. Amino acids whose codons differed by a single base in the first and third codon positions were very similar with respect to polar requirement and hydropathy. The major differences between amino acids were specified by the second codon position. Codons with U in the second position are hydrophobic, whereas most codons with A in the second position are hydrophilic. This accounts for the observation of complementary hydropathy. Single-base changes in the natural code had a smaller average effect on polar requirement than all but 0.02% of random codes. This result is most easily explained by selection to minimize deleterious effects of translation errors during the early evolution of the code.     

Surface, subunit interfaces and interior of oligomeric proteins

The solvent-accessible surface area (As) of 23 oligomeric proteins is calculated using atomic co-ordinates from high-resolution and well-refined crystal structures. As is correlated with the protein molecular weight, and a power law predicts its value to within 5% on average. The accessible surface of the average oligomer is similar to that of monomeric proteins in its hydropathy and amino acid composition. The distribution of the 20 amino acid types between the protein surface and its interior is also the same as in monomers. Interfaces, i.e. surfaces involved in subunit contacts, differ from the rest of the subunit surface. They are enriched in hydrophobic side-chains, yet they contain a number of charged groups, especially from Arg residues, which are the most abundant residues at interfaces except for Leu. Buried Arg residues are involved in H-bonds between subunits. We counted H-bonds at interfaces and found that several have none, others have one H-bond per 200 A2 of interface area on average (1 A = 0.1 nm). A majority of interface H-bonds involve charged donor or acceptor groups, which should make their contribution to the free energy of dissociation significant, even when they are few. The smaller interfaces cover about 700 A2 of the subunit surface. The larger ones cover 3000 to 10,000 A2, up to 40% of the subunit surface area in catalase. The lower value corresponds to an estimate of the accessible surface area loss required for stabilizing subunit association through the hydrophobic effect alone. Oligomers with small interfaces have globular subunits with accessible surface areas similar to those of monomeric proteins. We suggest that these oligomers assemble from preformed monomers with little change in conformation. In oligomers with large interfaces, isolated subunits should be unstable given their excessively large accessible surface, and assembly is expected to require major structural changes.     

Sequence characteristics responsible for protein-protein interactions in the intrinsically disordered regions of caseins,amelogenins, and small heat-shock proteins

Milk caseins and dental amelogenins are intrinsically disordered proteins (IDPs) that associate with themselves and others. Paradoxically, they are also described as hydrophobic proteins, which is difficult to reconcile with a solvent-exposed conformation. We attempt to resolve this paradox. We show that caseins and amelogenins are not hydrophobic proteins but they are more hydrophobic than most IDPs. Remarkably, uncharged residues from different regions of these mature proteins have a nearly constant average hydropathy but these regions exhibit different charged residue frequencies. A novel sequence analysis method was developed to identify hydrophobic and order-promoting regions that would favor conformational collapse. We found that such regions were uncommon; most hydrophobic and order-promoting residues were adjacent to hydrophilic or disorder-promoting residues. A further reason why caseins and amelogenins do not collapse is their high proportion of disorder-promoting proline residues. We conclude that in these proteins the hydrophobic effect is not large enough to cause conformational collapse but it can contribute, along with polar interactions, to protein-protein interactions. This behaviour is similar to the interaction of the disordered N-terminal region of small heat-shock proteins with either themselves during oligomer formation or other, unfolding, proteins during chaperone action.     

Separation of membrane proteins according to their hydropathy by serial phase partitioning with Triton X-114

The detergent Triton X-114, because of its convenient cloud point temperature (22 °C), has been used extensively to extract membrane proteins and to separate them in two phases according to their hydropathy. The upper detergent-poor phase contains mostly hydrophilic proteins, whereas hydrophobic ones are found mainly in the lower detergent-rich phase. In this work, we developed a method to fractionate membrane proteins and estimate their hydropathy based on a series of cloud point partitions with Triton X-114. With this method, beetroot plasma membrane proteins were separated in different fractions according to their hydropathy, following the binomial distribution law as expected. This method revealed the presence of both hydrophilic and hydrophobic Ca²⁺-dependent protein kinases in those membranes. At least five distinct Ca²⁺-dependent kinases were observed in in-gel kinase activity assays. This separation procedure was also used as the first step in the purification of a hydrophobic 60-kDa kinase.     

Location of helical regions in tetrapyrrole-containing proteins by a helical hydrophobic moment analysis. Application to phytochrome

Helical regions in many tetrapyrrole proteins are highly amphiphilic, one side interacting with a hydrophobic core and another side interacting with the polar solvent. The mean helical hydrophobic moment is a measure of amphiphilicity of a helix. Helical regions in myoglobin, the alpha and beta subunits of C-phycocyanin, and cytochrome c can be distinguished from nonhelical regions by use of a hydrophobic moment analysis. 24 of 27 (89%) of the helical regions in these proteins were located by this analysis. Calculations were also performed on chymotrypsin, ribonuclease, and papain, which do not possess as pronounced a hydrophobic core as the tetrapyrrole-containing proteins. Less than 50% of the helical regions were correctly located, indicating a lack of amphiphilicity in the helices of these proteins. The hydrophobic moment analysis was also used to predict helical regions in phytochrome, the ubiquitous photoreceptor in plants. Additionally, this analysis is used to quickly locate internal hydrophilic residues which may be functionally important. The distribution of hydrophobic moments from a random sequence was determined so that qualitative and to some extent quantitative comparisons between different amphiphilic helices may be made.     

Hydrophobic patches on protein subunit interfaces: Characteristics and prediction

Hydrophobic patches, defined as clusters of neighboring apolar atoms deemed accessible on a given protein surface, have been investigated on protein subunit interfaces. The data were taken from known tertiary structures of multimeric protein complexes. Amino acid composition and preference, patch size distribution, and patch contact complementarity across associating subunits were examined and compared with hydrophobic patches found on the solvent-accessible surface of the multimeric complexes. The largest or second largest patch on the accessible surface of the entire subunit was involved in multimeric interfaces in 90% of the cases. These results should prove useful for subunit design and engineering as well as for prediction of subunit interface regions. Proteins 28:333–343, 1997. © 1997 Wiley-Liss, Inc.     

CARd: Carbon distribution analysis program for protein sequences

Carbon distribution is responsible for stability and structure of proteins. Arrangement of carbon along the protein sequence is depends on how the amino acids are organized and is guided by mRNAs. An atomic level revision is important for understanding these codes. This will ultimately help in identification of disorders and suggest mutations. For this purpose a carbon distribution analysis program has been developed. This program captures the hydrophobic / hydrophilic / disordered regions in a protein. The program gives accurate results. The calculations are precise and sensitive to single amino acid resolution. This program is to help in mutational studies leading to protein stabilisation.     

Classification of 29 families of secondary transport proteins into a single structural class using hydropathy profile analysis

A classification scheme for membrane proteins is proposed that clusters families of proteins into structural classes based on hydropathy profile analysis. The averaged hydropathy profiles of protein families are taken as fingerprints of the 3D structure of the proteins and, therefore, are able to detect more distant evolutionary relationships than amino acid sequences. A procedure was developed in which hydropathy profile analysis is used initially as a filter in a BLAST search of the NCBI protein database. The strength of the procedure is demonstrated by the classification of 29 families of secondary transporters into a single structural class, termed ST[3]. An exhaustive search of the database revealed that the 29 families contain 568 unique sequences. The proteins are predominantly from prokaryotic origin and most of the characterized transporters in ST[3] transport organic and inorganic anions and a smaller number are Na(+)/H(+) antiporters. All modes of energy coupling (symport, antiport, uniport) are found in structural class ST[3]. The relevance of the classification for structure/function prediction of uncharacterised transporters in the class is discussed.     

Structure and dynamics of the fatty acid binding cavity in apo rat intestinal fatty acid binding protein.

The structure and dynamics of the fatty acid binding cavity in I-FABP (rat intestinal fatty acid binding protein) were analyzed. In the crystal structure of apo I-FABP, the probe occupied cavity volume and surface are 539+/-8 A3 and 428 A2, respectively (1.4 A probe). A total of 31 residues contact the cavity with their side chains. The side-chain cavity surface is partitioned according to the residue type as follows: 36-39% hydrophobic, 21-25% hydrophilic, and 37-43% neutral or ambivalent. Thus, the cavity surface is neither like a typical protein interior core, nor is like a typical protein external surface. All hydrophilic residues that contact the cavity-with the exception of Asp74-are clustered on the one side of the cavity. The cavity appears to expand its hydrophobic surface upon fatty acid binding on the side opposite to this hydrophilic patch. In holo I-FABP the fatty acid chain interactions with the hydrophilic side chains are mediated by water molecules. Molecular dynamics (MD) simulation of fully solvated apo I-FABP showed global conformational changes of I-FABP, which resulted in a large, but seemingly transient, exposure of the cavity to the external solvent. The packing density of the side chains lining the cavity, studied by Voronoi volumes, showed the presence of two distinctive small hydrophobic cores. The MD simulation predicts significant structural perturbations of the cavity on the subnanosecond time scale, which are capable of facilitating exchange of I-FABP internal water.     

Are membrane proteins "inside-out" proteins?

One of the central paradigms of structural biology is that membrane proteins are "inside-out" proteins, in that they have a core of polar residues surrounded by apolar residues. This is the reverse of the characteristics found in water-soluble proteins. We have decided to test this paradigm, now that sufficient numbers of transmembrane alpha-helical structures are accessible to statistical analysis. We have analyzed the correlation between accessibility and hydrophobicity of both individual residues and complete helices. Our analyses reveal that hydrophobicity of residues in a transmembrane helical bundle does not correlate with any preferred location and that the hydrophilic vector of a helix is a poor indicator of the solvent exposed face of a helix. Neither polar nor hydrophobic residues show any bias for the exterior or the interior of a transmembrane domain. As a control, analysis of water-soluble helical bundles performed in a similar manner has yielded clear correlations between hydrophobicity and accessibility. We therefore conclude that, based on the data set used, membrane proteins as "inside-out" proteins is an unfounded notion, suggesting that packing of alpha-helices in membranes is better understood by maximization of van der Waal's forces, rather than by a general segregation of hydrophobicities driven by lipid exclusion.     

Protein patterning on silicon-based surface using background hydrophobic thin film

A new and convenient protein patterning method on silicon-based surface was developed for protein array by spin coating of hydrophobic thin film (CYTOP). Photolithographic lift-off process was used to display two-dimensional patterns of spatially hydrophilic region. The background hydrophobic thin film was used to suppress nonspecific protein binding, and the hydrophilic target protein binding region was chemically modified to introduce aldehyde group after removal of the photoresist layer. The difference in surface energy between the hydrophilic pattern and background hydrophobic film would induce easier covalent binding of proteins onto defined hydrophilic areas having physical and chemical constraints. Below 1 microg/ml of total protein concentration, the CYTOP hydrophobic film effectively suppressed nonspecific binding of the protein. During the process of protein patterning, inherent property of the hydrophobic thin film was not changed judging from static and dynamic contact angle survey. Quantitative analysis of the protein binding was demonstrated by streptavidin-biotin system.     

Biological implications of complementary hydropathy of amino acids

The principle of complementary hydropathy predicts that peptides coded for by opposing DNA strands will bind one another because highly hydrophilic amino acids will be complemented by hydrophobic ones and vice versa. This paper provides the chemical plausibility for such interactions. It is suggested that exons coding for interacting peptides were juxtaposed and co-evolved together. Present day genes are no longer thus arranged because of duplications and exon shuffling.     

In silico inference of inclusion membrane protein family in obligate intracellular parasites chlamydiae.

Chlamydiae are obligate intracellular pathogens that proliferate only inside a vacuole, called an inclusion. Chlamydial Inc proteins are known to be a major component of the inclusion membrane, but little is known about the gene number and function. The Inc proteins share very low sequence similarity but a similar hydropathy profile among them. Using the hydropathy profile, we computationally searched the open reading frames (ORFs) having a similar profile and predicted 90 and 36 ORFs (Inc-like ORFs) as candidates for Inc proteins in Chlamydia pneumoniae J138 and Chlamydia trachomatis serovar D, respectively. On the other hand, only a few Inc-like ORFs were found in organisms other than chlamydiae, suggesting that the Inc-like ORFs are specific to chlamydiae. Comparative genome analysis also revealed that the Inc-like ORFs have multiplied and diverged as paralogues and orthologues in the chlamydial genomes, and that some Inc-like ORFs lacked the N-terminal portion or encoded the split form. The data suggest that these gene products constitute a large protein family and may play an important role in chlamydial infection, growth and survival in the host cell.