A fast method for the quantitative estimation of the distribution of hydrophobic and hydrophilic segments in alpha-helices of membrane proteins

The work presents a fast quantitative approach for estimating the orientations of hydrophilic and hydrophobic regions in the helical wheels of membrane-spanning alpha-helices of transmembrane proteins. The common hydropathy analysis provides an estimate of the integral hydrophobicity in a moving window which scans an amino acid sequence. The new parameter, orientation hydrophobicity, is based on the estimate of hydrophobicity of the angular segment that scans the helical wheel of a given amino acid sequence. The corresponding procedure involves the treatment of transmembrane helices as cylinders with equal surface elements for each amino acid residue. The orientation hydrophobicity, P(phi), phi = 0-360 degrees, of a helical cylinder is given as a sum of hydrophobicities of individual amino acids which are taken as the S-shaped functions of the angle between the centre of amino acid surface element and the centre of the segment. Non-zero contribution to P(phi) comes only from the amino acids belonging to the angular segment for a given angle phi. The size of the angular segment is related to the size of the channel pore. The amplitudes of amino acid S-functions are calibrated in the way that their maximum values (reached when the amino acid is completely exposed into the pore) are equal to the corresponding hydropathy index in the selected scale (here taken as Goldman-Engelman-Steitz hydropathy scale). The given procedure is applied in the studies of three ionic channels with well characterized three-dimensional structures where the channel pore is formed by a bundle of alpha-helices: cholera toxin B, nicotinic acetylcholine homopentameric alpha7 receptor, and phospholamban. The estimated maximum of hydrophilic properties at the helical wheels are in a good agreement with the spatial orientations of alpha-helices in the corresponding channel pores.     

Hydrophobicity and prediction of the secondary structure of membrane proteins and peptides.

Reliability of the hydropathy method to predict the formation of membrane-spanning alpha-helices by integral membrane proteins and peptides whose structure is known from X-ray crystallography is analysed. It is shown that Kyte-Doolittle hydropathy plots do not predict accurately 22 transmembrane alpha-helices in the reaction centres (RC) of the photosynthetic bacteria Rhodopseudomonas viridis and Rhodobacter sphaeroides (R-26). The accuracy of prediction for these proteins was improved using an optimised Kyte-Doolittle hydrophobicity scale. However, this hydrophobicity scale did not improve the predictions for the alphabeta-peptides of the B800-850 (LH2) complexes of the photosynthetic bacteria Rhodopseudomonas acidophila and Rhodospirillum molischianum, which were excluded from the optimisation procedure. The best and worst predictions of membrane-spanning alpha-helices for the RC proteins and LH2 peptides, respectively, were obtained with a propensity scale (PRC) calculated from the amino acid sequences and X-ray data for the RC proteins. A propensity scale (PLH) obtained using the amino acid sequences and X-ray data for the alphabeta-peptides of the LH2 complexes did not give an acceptable prediction of the transmembrane segments in the LH2 peptides; moreover, it markedly contradicted the PRC scale. Amino acids have been concluded to have no significant preference to localisation in transmembrane segments. Therefore, the predictive ability of the hydropathy methodology appears to be limited: the number of transmembrane segments can be correctly calculated for the best case only, and the lengths and positions of membrane-spanning alpha-helices in a protein amino acid sequence can not be predicted exactly.     

What Amino Acid Properties Affect Protein Evolution?

We studied 10 protein-coding mitochondrial genes from 19 mammalian species to evaluate the effects of 10 amino acid properties on the evolution of the genetic code, the amino acid composition of proteins, and the pattern of nonsynonymous substitutions. The 10 amino acid properties studied are the chemical composition of the side chain, two polarity measures, hydropathy, isoelectric point, volume, aromaticity, aliphaticity, hydrogenation, and hydroxythiolation. The genetic code appears to have evolved toward minimizing polarity and hydropathy but not the other seven properties. This can be explained by our finding that the presumably primitive amino acids differed much only in polarity and hydropathy, but little in the other properties. Only the chemical composition (C) and isoelectric point (IE) appear to have affected the amino acid composition of the proteins studied, that is, these proteins tend to have more amino acids with typical C and IE values, so that nonsynonymous mutations tend to result in small differences in C and IE. All properties, except for hydroxythiolation, affect the rate of nonsynonymous substitution, with the observed amino acid changes having only small differences in these properties, relative to the spectrum of all possible nonsynonymous mutations. Received: 2 January 1998 / Accepted: 25 April 1998     

嗜热与嗜常温微生物的蛋白质氨基酸组成比较

嗜热微生物的嗜热特性与其蛋白质的高度热稳定性紧密相关。为了探索嗜热蛋白质的热稳定机制,比较嗜热和嗜常温微生物的蛋白质在氨基酸组成上的差别,收集110对分别来自嗜热和嗜常温微生物的同源蛋白质序列,比较两组蛋白质各种氨基酸含量以及疏水性氨基酸组成、疏水性指数和荷电氨基酸组成的差别,结果两者在多种氨基酸含量上存在微小但统计学上显著的差别,嗜热蛋白质比嗜常温蛋白质具有较高的平均疏水性和荷电氨基酸组成。对两组蛋白质的“脂肪族氨基酸指数”进行分析,证明嗜热蛋白质之所以具有较高的脂肪族氨基酸指数是由于其亮氨酸含量较高,与影响该指数的其它几种氨基酸无关;从而认为该指数的意义值得怀疑。通过对大量同源嗜热蛋白质和嗜常温蛋白质氨基酸组成的比较,能够揭示一些有关蛋白质热稳定性的普遍规律。     

The role of leucine and isoleucine in tuning the hydropathy of class A GPCRs

Leucine and Isoleucine are two amino acids that differ only by the positioning of one methyl group. This small difference can have important consequences in α-helices, as the β-branching of Ile results in helix destabilization. We set out to investigate whether there are general trends for the occurrences of Leu and Ile residues in the structures and sequences of class A GPCRs (G protein-coupled receptors). GPCRs are integral membrane proteins in which α-helices span the plasma membrane seven times and which play a crucial role in signal transmission. We found that Leu side chains are generally more exposed at the protein surface than Ile side chains. We explored whether this difference might be attributed to different functions of the two amino acids and tested if Leu tunes the hydrophobicity of the transmembrane domain based on the Wimley-White whole-residue hydrophobicity scales. Leu content decreases the variation in hydropathy between receptors and correlates with the non-Leu receptor hydropathy. Both measures indicate that hydropathy is tuned by Leu. To test this idea further, we generated protein sequences with random amino acid compositions using a simple numerical model, in which hydropathy was tuned by adjusting the number of Leu residues. The model was able to replicate the observations made with class A GPCR sequences. We speculate that the hydropathy of transmembrane domains of class A GPCRs is tuned by Leu (and to some lesser degree by Lys and Val) to facilitate correct insertion into membranes and/or to stably anchor the receptors within membranes.     

Protein map: an alignment-free sequence comparison method based on various properties of amino acids

In this paper, we propose a new protein map which incorporates with various properties of amino acids. As a powerful tool for protein classification, this new protein map both considers phylogenetic factors arising from amino acid mutations and provides computational efficiency for the huge amount of data. The ten amino acid physico-chemical properties (the chemical composition of the side chain, two polarity measures, hydropathy, isoelectric point, volume, aromaticity, aliphaticity, hydrogenation, and hydroxythiolation) are utilized according to their relative importance. Moreover, during the course of calculation of genetic distances between pairs of proteins, this approach does not require any alignment of sequences. Therefore, the proposed model is easier and quicker in handling protein sequences than multiple alignment methods, and gives protein classification greater evolutionary significance at the amino acid sequence level.     

Physicochemical characterization and functional analysis of some snake venom toxin proteins and related non-toxin proteins of other chordates

Snake venom contains a diverse array of proteins and polypeptides. Cytotoxins and short neurotoxins are non-enzymatic polypeptide components of snake venom. The three-dimensional structure of cytotoxin and short neurotoxin resembles a three finger appearance of three-finger protein super family. Different family members of three-finger protein super family are employed in diverse biological functions. In this work we analyzed the cytotoxin, short neurotoxin and related non-toxin proteins of other chordates in terms of functional analysis, amino acid compositional (%) profile, number of amino acids, molecular weight, theoretical isoelectric point (pI), number of positively charged and negatively charged amino acid residues, instability index and grand average of hydropathy with the help of different bioinformatical tools. Among all interesting results, profile of amino acid composition (%) depicts that all sequences contain a conserved cysteine amount but differential amount of different amino acid residues which have a family specific pattern. Involvement in different biological functions is one of the driving forces which contribute the vivid amino acid composition profile of these proteins. Different biological system dependent adaptation gives the birth of enriched bio-molecules. Understanding of physicochemical properties of these proteins will help to generate medicinally important therapeutic molecules for betterment of human lives.     

Role of mutational bias and natural selection on genome-wide nucleotide bias in prokaryotic organisms

Correlations between genomic GC contents and amino acid frequencies were studied in the homologous sequences of 12 eubacterial genomes. Results show that amino acids encoded by GC-rich codons increases significantly with genomic GC contents, whereas opposite trend was observed in case of amino acids encoded by GC-poor codons. Further studies show all the amino acids do not change in the predicted direction according to their genomic GC pressure, suggesting that protein evolution is not entirely dictated by their nucleotide frequencies. Amino acid substitution matrix calculated among hydrophobic, amphipathic and hydrophilic amino acid groups' shows that amphipathic and hydrophilic amino acids are more frequently substituted by hydrophobic amino acids than from hydrophobic to hydrophilic or amphipathic amino acids. This indicates that nucleotide bias induces a directional changes in proteome composition in such a way that underwent strong changes in hydropathy values. In fact, significant increases in hydrophobicity values have also been observed with the increase of genomic GC contents. Correlations between GC contents and amino acid compositions in three different predicted protein secondary structures show that hydropathy values increases significantly with GC contents in aperiodic and helix structures whereas strand structure remains insensitive with the genomic GC levels. The relative importance of mutation and selection on the evolution of proteins have been discussed on the basis of these results.     

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.     

Aminoacyl-tRNA synthetase classes and groups in prokaryotes

Knowledge on the evolution of aminoacyl-tRNA synthetases is crucial to studies on the origins of life. The relationships between the different aminoacyl-tRNA synthetase specificities in prokaryotic organisms are studied in this work. We reconstructed the ancestor sequences and the phylogenetic relationships utilizing the Maximum Likelihood method. The results suggest that in class I the evolution of the N-terminal segment was strongly influenced by the amino acid hydropathy in both domains of prokaryotes. The results for the C-terminal segments of class I were different in the two domains, indicating that its evolution was strongly influenced by the specific types of tRNA modification in each domain. The class II groups in Archaea were more heterogeneous with respect to the hydropathy of amino acids, indicating the interference of other influences. In bacteria, the configuration was also complex but the overall consensual division in two groups was maintained, group IIa forming a single branch with the five hydroapathetic amino acid specificities and group IIb containing the specificities for the moderately hydrophobic together with the hydrophilic amino acids. It is indicated that the aminoacyl-tRNA synthetase in both domains were subjected to different selective forces in diverse parts of the proteins, resulting in complex phylogenetic patterns.     

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.     

Trends of amino acid usage in the proteins from the unicellular parasite Giardia lamblia

Correspondence analysis of amino acid frequencies was applied to 75 complete coding sequences from the unicellular parasite Giardia lamblia, and it was found that three major factors influence the variability of amino acidic composition of proteins. The first trend strongly correlated with (a) the cysteine content and (b) the mean weight of the amino acids used in each protein. The second trend correlated with the global levels of hydropathy and aromaticity of each protein. Both axes might be related with the defense of the parasite to oxygen free radicals. Finally, the third trend correlated with the expressivity of each gene, indicating that in G. lamblia highly expressed sequences display a tendency to preferentially use a subset of the total amino acids.     

Protein-based phylogenetic analysis by using hydropathy profile of amino acids

So far, various approaches for phylogenetic analysis have been developed. Almost all of them put stress on analyzing nucleic acid sequences or protein primary structures. In this paper, we take the physicochemical properties of amino acids into account and introduce the hydropathy profile of amino acids into phylogenetic analysis. We find that this introduction is effectual and our method may be used to complement phylogenetic analysis.     

Hydrophilicity of polar amino acid side-chains is markedly reduced by flanking peptide bonds

Several amino acid side-chain hydropathy scales have been devised on the basis of solubility and water/organic solvent partitioning data obtained with free amino acids or side-chain analogs. In nearly all cases, these scales are based upon the structure-additivity assumption; it has been assumed that the transfer free energies of the amino acid side-chains are the same in these model compounds as they are in a polypeptide. This assumption is probably wrong. In the present study, deviations from additivity for amino acid side-chains are demonstrated by comparing a theoretically derived scale, which N-acetylamino acid amides. The results show that the flanking peptide bonds dramatically reduce the hydrophilicity of the polar side-chains, with deviations up to several kilocalories (1 kcal = 4.184 kJ) for the charged side-chains at pH 7.0. Further calculation shows that these deviations are due to reductions of 40 to 85% in the unfavorable transfer free energy of the polar functional groups. In addition, proximity of the neighboring amide bonds in the parent molecule (N-acetylglycine amide) decreases the hydrophilicity of the -CONH-backbone unit by 36%. This decrease is expected to be twice as large for -CONH- units in the interior of a polypeptide backbone. The significance of these observations is: (1) valid hydropathy scales can be obtained only with model peptides; (2) deviations from additivity are expected in all solvent systems, including non-polar solvents that are thought to mimic the interior of a membrane; (3) the spontaneous insertion of polypeptides into membranes is likely to occur much more readily than has been previously thought. In order to estimate the free energy of transferring the side-chains and the polypeptide backbone from water to the interior of a lipid bilayer, the results of this study are used to construct a hydropathy scale based upon the partitioning of solutes between water and non-polar solvents. The validity of hydropathy scales that are based on criteria other than solubility and water/organic solvent partitioning data is also discussed.     

Detection of a genome-linked protein (VPg) of hepatitis A virus and its comparison with other picornaviral VPgs.

The nucleotide sequence corresponding to the P3 region of the hepatitis A virus (HAV) polyprotein genome was determined from cloned cDNA and translated into an amino acid sequence. Comparison of the amino acid sequences of the genome-linked proteins (VPgs) of other picornaviruses with the predicted amino acid sequence of HAV was used to locate the primary structure of a putative VPg within the genome of HAV. The sequence of HAV VPg, like those of other picornaviral VPg molecules, contains a tyrosine residue as a potential binding site for HAV RNA in position 3 from its N terminus. The potential cleavage sites to generate VPg from a putative HAV polyprotein are between glutamic acid and glycine at the N terminus and glutamic acid and serine or glutamine and serine at the C terminus. A synthetic peptide corresponding to 10 amino acids of the predicted C terminus of HAV VPg induced anti-peptide antibodies in rabbits when it was conjugated to thyroglobulin as a carrier. These antibodies were specific for the peptide and precipitated VPg, linked to HAV RNA, from purified HAV and from lysates of HAV-infected cells. The precipitation reaction was blocked by the synthetic peptide (free in solution or coupled to carrier proteins) and prevented by pretreatment of VPg RNA with protease. Thus, our predicted amino acid sequence is colinear with the nucleotide sequence of the VPg gene in the HAV genome. From our results we concluded that HAV has the typical organization of picornavirus genes in this part of its genome. Similarity among hydrophobicity patterns of amino acid sequences of different picornaviral VPgs was revealed in hydropathy plots. Thus, the VPg of HAV appears to be closely related to VPg1 and VPg2 of foot-and-mouth disease virus. In contrast, HAV VPg has a unique isoelectric point (pI = 7.15) among the picornavirus VPgs.     

A strong effect of AT mutational bias on amino acid usage in Buchnera is mitigated at high-expression genes

The advent of full genome sequences provides exceptionally rich data sets to explore molecular and evolutionary mechanisms that shape divergence among and within genomes. In this study, we use multivariate analysis to determine the processes driving genome-wide patterns of amino usage in the obligate endosymbiont Buchnera and its close free-living relative Escherichia coli. In the AT-rich Buchnera genome, the primary source of variation in amino acid usage differentiates high- and low-expression genes. Amino acids of high-expression Buchnera genes are generally less aromatic and use relatively GC-rich codons, suggesting that selection against aromatic amino acids and against amino acids with AT-rich codons is stronger in high-expression genes. Selection to maintain hydrophobic amino acids in integral membrane proteins is a primary factor driving protein evolution in E. coli but is a secondary factor in Buchnera. In E. coli, gene expression is a secondary force driving amino acid usage, and a correlation with tRNA abundance suggests that translational selection contributes to this effect. Although this and previous studies demonstrate that AT mutational bias and genetic drift influence amino acid usage in Buchnera, this genome-wide analysis argues that selection is sufficient to affect the amino acid content of proteins with different expression and hydropathy levels.     

Looking at structure, stability, and evolution of proteins through the principal eigenvector of contact matrices and hydrophobicity profiles

We review and further develop an analytical model that describes how thermodynamic constraints on the stability of the native state influence protein evolution in a site-specific manner. To this end, we represent both protein sequences and protein structures as vectors: structures are represented by the principal eigenvector (PE) of the protein contact matrix, a quantity that resembles closely the effective connectivity of each site; sequences are represented through the "interactivity" of each amino acid type, using novel parameters that are correlated with hydropathy scales. These interactivity parameters are more strongly correlated than the other hydropathy scales that we examine with: (1) the change upon mutations of the unfolding free energy of proteins with two-states thermodynamics; (2) genomic properties as the genome-size and the genome-wide GC content; (3) the main eigenvectors of the substitution matrices. The evolutionary average of the interactivity vector correlates very strongly with the PE of a protein structure. Using this result, we derive an analytic expression for site-specific distributions of amino acids across protein families in the form of Boltzmann distributions whose "inverse temperature" is a function of the PE component. We show that our predictions are in agreement with site-specific amino acid distributions obtained from the Protein Data Bank, and we determine the mutational model that best fits the observed site-specific amino acid distributions. Interestingly, the optimal model almost minimizes the rate at which deleterious mutations are eliminated by natural selection.     

Spasmin and a Putative Spasmin Binding Protein(s) Isolated from Solubilized Spasmonemes1

ABSTRACT The amino acid composition and hydrophobicity scale (hydropathy) of calcium-binding proteins contained in the contractile spasmoneme of Carchesium polypinum was compared with other calcium-binding proteins from eukaryotes. Spasmins which may hind at most 4 calcium ions simultaneously and initiate spasmoneme contraction cooperatively belong to a super family of proteins including; centrin/caltractin and calmodulin. Based on chemical modification of tryptophan and methionine, these residues are involved in contraction but the spasmin proteins contain little or none of these amino acids. Based on this evidence, it is suggested that another, non-calcium binding protein(s) is involved in spasmoneme contraction.     

Plasmodium falciparum: gene structure and hydropathy profile of the major merozoite surface antigen (gp195) of the Uganda-Palo Alto isolate

The gene encoding the 195,000-Da major merozoite surface antigen (gp195) of the FUP (Uganda-Palo Alto) isolate of Plasmodium falciparum, a strain widely used for monkey vaccination experiments, has been cloned and sequenced. The translated amino acid sequence of the FUP gp195 protein is closely related to the sequences of corresponding proteins of the CAMP (Malaysia) and MAD-20 (Papua New Guinea) isolates and more distantly related to those of the Wellcome (West Africa) and K1 (Thailand) isolates, supporting the proposed allelic dimorphism of gp195 within the parasite population. The prevalence of dimorphic sequences within the gp195 protein suggests that many gp195 epitopes would be group-specific. Despite the extensive differences in amino acid sequence between gp195 proteins of these two groups, the hydropathy profiles of proteins representative of both groups are very similar. The conservation of overall secondary structure shown by the hydropathy profile comparison indicates that gp195 proteins of the various P. falciparum isolates are functionally equivalent. This information on the primary structure of the FUP gp195 protein will enable us to evaluate the possible roles of conserved, group-specific and variable epitopes in immunity to the blood stage of the malaria parasite.     

Characterization and functional expression in Escherichia coli of the sodium/proton/glutamate symport proteins of Bacillus stearothermophilus and Bacillus caldotenax

The genes encoding the Na+/H+/L-glutamate symport proteins of the thermophilic organisms Bacillus stearothermophilus (gltTBs) and Bacillus caldotenax (gltTBc) were cloned by complementation of Escherichia coli JC5412 for growth on glutamate as sole source of carbon, energy and nitrogen. The nucleotide sequences of the gltTBs and gltTBc genes were determined. In both cases the translated sequences corresponded with proteins of 421 amino acid residues (96.7% amino acid identity between GltTBs and GltTBc). Putative promoter, terminator and ribosome-binding-site sequences were found in the flanking regions. These expression signals were functional in E. coli. The hydropathy profiles indicate that the proteins are hydrophobic and could form 12 membrane-spanning regions. The Na+/H+ coupled L-glutamate symport proteins GltTBs and GltTBc are homologous to the strictly H+ coupled L-glutamate transport protein of E. coli K-12 (overall 57.2% identity). Functional expression of glutamate transport activity was demonstrated by uptake of glutamate in whole cells and membrane vesicles. In accordance with previous observations (de Vrij et al., 1989; Heyne et al., 1991), glutamate uptake was driven by the electrochemical gradients of sodium ions and protons.