Role of large hydrophobic residues in proteins

Large Hydrophobic Residues (LHR) such as phenylalanine, isoleucine, leucine, methionine and valine play an important rolein protein structure and activity. We describe the role of LHR in complete set of protein sequences in 15 different species.That is the distribution of LHR in different proteins of different species is reported. It is observed that the proteins prefer tohave 27% of large hydrophobic residues in total and all along the sequence. It is also observed that proteins accumulate moreLHR in its active sites. A window analysis on these protein sequences shows that the 27% of LHR is more frequent atwindow length of 45 amino acids. The influenza virus and P. falciparum show a random distribution of LHR in its proteinscompared to other model organisms.     

Strong hydrophobic nature of cysteine residues in proteins.

The differences between disulfide-bonding cystine (Cys_SS) and free cysteine (Cys_SH) residues were examined by analyzing the statistical distribution of both types of residue in proteins of known structure. Surprisingly, Cys_SH residues display stronger hydrophobicity than Cys_SS residues. A detailed survey of atoms which come into contact with the sulfhydryl group (sulfur atom) of Cys_SH revealed those atoms are essentially the same in number and variety as those of the methyl group of isoleucine, but are quite different to those of the hydroxyl group of serine. Moreover, the relationships among amino acids were also determined using the 3D-profile table of known protein structures. Cys_SH was located in the hydrophobic cluster, along with residues such as Met, Trp and Tyr, and was clearly separated from Ser and Thr in the polar cluster. These results imply that free cysteines behave as strongly hydrophobic, and not hydrophilic, residues in proteins.     

Native-like mean structure in the unfolded ensemble of small proteins

The nature of the unfolded state plays a great role in our understanding of proteins. However, accurately studying the unfolded state with computer simulation is difficult, due to its complexity and the great deal of sampling required. Using a supercluster of over 10,000 processors we have performed close to 800 micros of molecular dynamics simulation in atomistic detail of the folded and unfolded states of three polypeptides from a range of structural classes: the all-alpha villin headpiece molecule, the beta hairpin tryptophan zipper, and a designed alpha-beta zinc finger mimic. A comparison between the folded and the unfolded ensembles reveals that, even though virtually none of the individual members of the unfolded ensemble exhibits native-like features, the mean unfolded structure (averaged over the entire unfolded ensemble) has a native-like geometry. This suggests several novel implications for protein folding and structure prediction as well as new interpretations for experiments which find structure in ensemble-averaged measurements.     

Contribution of the free energy of mixing of hydrophobic side chains to the stability of the tertiary structure of proteins

The driving force for folding of polypeptide chains into their threedimensional compact units has been designated as being hydrophobic and a measure of the hydrophobic character of the constituent amino acids has been determined by relative solubility measurements. It has been found however that the hydrophobic character of a protein is not sufficient to account for the complete stabilization of the tertiary structure of proteins. It is suggested that if the free energy of mixing of the hydrophobic side chains in the interior of the protein is added to the free energy of desolvation, i.e. the hydrophobic free energy, then the total free energy of mixing and desolvation can account for the known stability of the tertiary structure of proteins.     

Conferment of folding ability to a naturally unfolded apocytochrome c through introduction of hydrophobic amino acid residues

Hyperthermophilic Aquifex aeolicus cytochrome c(555) (AA c(555)) exceptionally folds even in the apo state, unlike general cytochromes c including mesophilic Pseudomonas aeruginosa cytochrome c(551) (PA c(551)), which is structurally homologous to AA c(555) in the holo state. Here we hypothesized that the exceptional apo AA c(555) folding can be attributed to nine hydrophobic amino acid residues and proved this using a PA c(551) variant (denoted as PA-nh) carrying the nine hydrophobic residues at structurally corresponding positions. Circular dichroism experiments showed that the apo PA-nh variant became folded, unlike the wild-type apo PA c(551), and exhibited much higher stability than the wild type. Another difference between the holo forms of AA c(555) and PA c(551) is the existence of an extra helix in the former. Introduction of the amino acid residues forming the extra helix of AA c(555) into the PA-nh variant did not significantly affect its folding ability in the apo state. Therefore, the nine hydrophobic residues introduced into the apo PA-nh variant were enough to confer the folding ability. PA c(551) represents the first example of the conversion of an intrinsically unfolded apocytochrome c into an autonomously folded one, which was revealed by means of a protein engineering method without heme. Although heme is generally considered to be a trigger of apocytochrome c folding, the present results demonstrate a new heme-independent folding mechanism.     

Solution structure of a llama single-domain antibody with hydrophobic residues typical of the VH/VL interface

The three-dimensional structure of a llama single-domain antibody BrucD4-4 was established by use of solution NMR spectroscopy. BrucD4-4 has Val, Gly, Leu, and Trp residues at positions 37, 44, 45, and 47, which are considered to be a hallmark to distinguish llama VH from V(H)H fragments at the germline level. In contrast to the murine and human VHs, BrucD4-4 has sufficient solubility, is monomeric in solution, and displays high-quality NMR spectra characteristic of well-structured proteins. Amide proton/deuterium exchange and the (15)N relaxation data showed that BrucD4-4 has a classic protein structure with a well-packed core and comparatively mobile surface loops. The three-dimensional architecture of BrucD4-4 is analogous to that of VHs from murine and human F(v)s and camelid V(H)Hs with two pleated beta-sheets formed by four and five beta-strands. A canonical and undistorted beta-barrel exposes a number of hydrophobic residues into the solvent on the surface of the three-dimensional structure. The eight-residue H3 loop folds over the side chain of Val37 similarly to that in llama V(H)Hs; however, this interaction may be transient due to the H3 conformational flexibility. Overall, the surface characteristics of BrucD4-4 with respect to hydrophobicity appear to lie between the human VH domain from Fv Pot and the llama V(H)H fragment HC-V, which may explain its enhanced solubility allowing NMR structural analysis.     

Synthesis and beta-conformation of copolypeptides with alternating hydrophilic and hydrophobic residues

Copolypeptides with alternating hydrophilic and hydrophobic residues were prepared, and their ability to form beta-structures in aqueous solutions was investigated by circular dichroism. Optically pure samples of poly (Lys-Leu-Lys-Leu) and poly (Leu-Glu-Leu-Glu), obtained via the 2-hydroxyphenyl esters, undergo a coil-to-beta transition in presence of salt. The beta-structures obtained under identical conditions with partially racemized samples of poly (Leu-Lys)Np and poly (Leu-Glu)Np, prepared by polycondensation of the corresponding dipeptide p-nitrophenyl esters, appear to be less regular. Non-alternating poly (Gly-Lys-Leu-Lys-Leu) does not form beta-structures in presence of NaCl as does alternating poly (Lys-Leu-Lys-Leu) indicating that the amino acid sequence can dramatically change the tendency to form beta-structures.     

Contribution of secondary structure to the heat capacity and enthalpy distribution of the unfolded state in proteins.

We have recently shown that one can construct the enthalpy distribution for protein molecules from experimental knowledge of the temperature dependence of the heat capacity. For many proteins the enthalpy distribution evaluated at the midpoint of the denaturation transition (corresponding to the maximum in the heat capacity vs temperature curve) is broad and biphasic, indicating two different populations of molecules (native and unfolded) with distinctly different enthalpies. At temperatures above the denaturation point, the heat capacity for the unfolded state in many proteins is quite large and using the analysis just mentioned, we obtain a gaussian-like enthalpy distribution that is very broad. A large value of the heat capacity indicates that there are structural changes going on in the unfolded state above the transition temperature. In the present paper we investigate the origin of this large heat capacity by considering the presence of changing amounts of secondary structure (specifically, alpha-helix) in the unfolded state. For this purpose we use the empirical estimates of the Zimm-Bragg sigma and s factors for all of the native amino acids in water as determined by Scheraga and co-workers. Using myoglobin as an example, we calculate probability profiles and distribution functions for the total number of helix states in the specific-sequence molecule. Given the partition function for the specific-sequence molecule, we can then calculate a set of enthalpy moments for the molecule from which we obtain a good estimate of the enthalpy distribution in the unfolded state. This distribution turns out to be quite narrow when compared with the distribution obtained from the raw heat capacity data. We conclude that there must be other major structural changes (backbone and solvent) that are not accounted for by the inclusion of alpha-helix in the unfolded state.     

Stabilization of secondary structure elements by specific combinations of hydrophilic and hydrophobic amino acid residues is more important for proteins encoded by GC-poor genes

Stabilization of secondary structure elements by specific combinations of hydrophobic and hydrophilic amino acids has been studied by the way of analysis of pentapeptide fragments from twelve partial bacterial proteomes. PDB files describing structures of proteins from species with extremely high and low genomic GC-content, as well as with average G + C were included in the study. Amino acid residues in 78,009 pentapeptides from alpha helices, beta strands and coil regions were classified into hydrophobic and hydrophilic ones. The common propensity scale for 32 possible combinations of hydrophobic and hydrophilic amino acid residues in pentapeptide has been created: specific pentapeptides for helix, sheet and coil were described. The usage of pentapeptides preferably forming alpha helices is decreasing in alpha helices of partial bacterial proteomes with the increase of the average genomic GC-content in first and second codon positions. The usage of pentapeptides preferably forming beta strands is increasing in coil regions and in helices of partial bacterial proteomes with the growth of the average genomic GC-content in first and second codon positions. Due to these circumstances the probability of coil-sheet and helix-sheet transitions should be increased in proteins encoded by GC-rich genes making them prone to form amyloid in certain conditions. Possible causes of the described fact that importance of alpha helix and coil stabilization by specific combinations of hydrophobic and hydrophilic amino acids is growing with the decrease of genomic GC-content have been discussed.     

Computational simulation of the statistical properties of unfolded proteins

A simple Monte Carlo method was used to generate ensembles of simulated polypeptide conformations that are restricted only by steric repulsion. The models used for these simulations were based on the sequences of four real proteins, ranging in size from 26 to 268 amino acid residues, and included all non-hydrogen atoms. Two sets of calculations were performed, one that included only intra-residue steric repulsion terms and those between adjacent residues, and one that included repulsion terms between all possible atom pairs, so as to explicitly account for the excluded volume effect. Excluded volume was found to increase the average radius of gyration of the chains by 20-40%, with the expansion factor increasing with chain length. Contrary to recent suggestions, however, the excluded volume effect did not greatly restrict the distribution of dihedral angles or favor native-like topologies. The average dimensions of the ensembles calculated with excluded volume were consistent with those measured experimentally for unfolded proteins of similar sizes under denaturing conditions, without introducing any adjustable scaling factor. The simulations also reproduced experimentally determined effective concentrations for the formation of disulfide bonds in reduced and unfolded proteins. The statistically generated ensembles included significant numbers of conformations that were nearly as compact as the corresponding native proteins, as well as many that were as accessible to solvent as a fully extended chain. On the other hand, conformations with as much buried surface area as the native proteins were very rare, as were highly extended conformations. These results suggest that the overall properties of unfolded proteins can be usefully described by a random coil model and that an unfolded polypeptide can undergo significant collapse while losing only a relatively small fraction of its conformational entropy.     

Induction of changes in the secondary structure of globular proteins by a hydrophobic surface

Circular dichroism, ellipsometry and radiolabeling techniques were employed to study the induction of changes in the secondary structure of BSA, myoglobin and cytochrome C by a hydrophobic surface. The results showed that adsorbed protein molecules lose their ordered native structure in the initial stage of adsorption and the structure appears to be a random or disordered conformation. Protein molecules adsorbed in later stages adopt a more ordered secondary structure ( helix and structure). The changes of secondary structure of globular proteins induced by a hydrophobic surface can be explained by the steric interaction between adsorbed proteins as well as by hydrophobic interactions during the adsorption process. In addition, there is obviously an intermediate stage in which the protein molecules are mainly in the structure, indicating that for certain proteins, the structure may be a more stable secondary structure than helix on the hydrophobic surface. Correspondence to: S.-F. Sui     

Three-dimensional structure and active site of three hydrophobic molecule-binding proteins with significant amino acid sequence similarity.

We review our work on bovine and human retinol-binding protein (RBP), bovine beta lactoglobulin (BLG), and bovine odorant-binding protein (OBP). These three proteins share a sequence similarity high enough to justify the proposal that their three-dimensional structure ought to be quite similar, and they also share the function of similar or even identical hydrophobic ligand binding, although with a very different degree of specificity. Thus they constitute an ideal system to exhaustively explore the question of three-dimensional structure prediction from sequence similarity and the related question of binding site prediction for similar ligands. We have used x-ray diffraction techniques on single crystals of human and bovine RBP, bovine milk BLG, and bovine nasal mucosa OBP to investigate this problem. The results of these crystallographic studies indicate that to the level of resolution so far attained, the three-dimensional structure of these three proteins is reasonably predicted from the sequence similarity. The fold is the same and structural differences are rather subtle. Finally, we present experimental evidence that the binding sites of RBP, BLG, and OBP are in different regions of the molecules. Thus, it appears that although sequence alignment has correctly predicted the protein fold, it has incorrectly predicted the hydrophobic ligand-binding sites.     

Detection of secondary structure elements in proteins by hydrophobic cluster analysis.

Hydrophobic cluster analysis (HCA) is a protein sequence comparison method based on alpha-helical representations of the sequences where the size, shape and orientation of the clusters of hydrophobic residues are primarily compared. The effectiveness of HCA has been suggested to originate from its potential ability to focus on the residues forming the hydrophobic core of globular proteins. We have addressed the robustness of the bidimensional representation used for HCA in its ability to detect the regular secondary structure elements of proteins. Various parameters have been studied such as those governing cluster size and limits, the hydrophobic residues constituting the clusters as well as the potential shift of the cluster positions with respect to the position of the regular secondary structure elements. The following results have been found to support the alpha-helical bidimensional representation used in HCA: (i) there is a positive correlation (clearly above background noise) between the hydrophobic clusters and the regular secondary structure elements in proteins; (ii) the hydrophobic clusters are centred on the regular secondary structure elements; (iii) the pitch of the helical representation which gives the best correspondence is that of an alpha-helix. The correspondence between hydrophobic clusters and regular secondary structure elements suggests a way to implement variable gap penalties during the automatic alignment of protein sequences.     

Extraction of detergents from hydrophobic proteins with isopentanol: application to electrophoretic analysis of photosynthetic bacterial hydrophobic proteins

The method for extracting Triton X-100 used by I. H. Mather and C. B. Tampling [Anal. Biochem. 93, 139-142 (1979)], has been extended to other detergents of different charge and chemical nature. All the detergents tested can be extracted with isopentanol in conditions in which not more than 8% of hydrophobic or hydrophilic protein is lost from the water phase. The removal of detergent from reaction centers and light harvesting protein-pigment complexes of photosynthetic bacteria, eliminates the artifacts of oligomers when analyzed by sodium dodecyl sulfate-gel electrophoresis.     

Characterization of the small hydrophobic proteins associated with pulmonary surfactant

Lipid extracts of bovine pulmonary surfactant, which retain many of the biophysical characteristics of natural surfactant, contain approx. 98% lipid and 2% protein, as determined by amino acid analysis. Polyacrylamide/urea gel electrophoresis reveals that lipid extract surfactant contained a major apoprotein band with apparent Mr 3500 and minor apoprotein bands with apparent Mr 15,000 and 7000. After reduction, the 15 kDa band disappears and is replaced by a prominent band with apparent Mr = 5000. Reduction also results in a relative diminution of the 7 kDa band and a relative increase in the intensity of the 3.5-kDa band. Edman degradation reveals two major peptide sequences which have been designated surfactant-associated peptide (N-terminal Phe) and surfactant-associated peptide (N-terminal Leu) and a minor sequence designated surfactant-associated peptide (N-terminal Ile). The latter surfactant-associated peptide appears to be related to the N-terminal Leu peptide but lacks the terminal Leu. N-Terminal analysis by dansylation demonstrates that the 15 and 5 kDa (reduced) apoprotein species contain N-terminal Phe, Leu and Ile. The 3.5 and 7 kDa bands contain only N-terminal Leu and Ile. Chromatography of lipid extracts on silicic acid columns gives rise to fraction I, which contains protein and phosphatidylglycerol, and fraction II, which contains protein, phosphatidylglycerol and phosphatidylethanolamine. Fraction I was primarily composed of the 15-kDa apoproteins, while fraction II contained mainly the 3.5 and 7 kDa apoproteins. Both fractions exhibited biophysical activity after reconstitution with dipalmitoylphosphatidylcholine. These results indicate that lipid extracts contain an oligomer of 15 kDa containing surfactant-associated peptide (N-terminal Phe) and surfactant-associated peptides (N-terminal Leu or Ile) which interact through sulfhydryl and perhaps other bonds. Lipid extracts also contain 3.5 kDa monomers of surfactant-associated peptides with N-terminal Leu and N-terminal Ile which can dimerize through sulfhydryl and perhaps hydrophobic interactions.     

Fluorescence and the structure of proteins. XXI. Fluorescence of aminotyrosyl residues in peptides and helical proteins.

1. Five peptides containing tyrosine were converted to the 3-aminotyrosyl peptides by nitration with tetranitromethane and subseuqent reduction of the nitro groups to amino groups. The fluorescence of these aminotyrosyl residues was found to be quite similar to that of 3-aminotyrosine and it is concluded that the fluorescence is not sensitive to incorporation of the amino acid into the peptide chain. 2. Fluorescence of 3-aminotyrosine derivatives was sensitive, however, to the nature of the solvent; as the dielectric constant decreased, fluorescence was enhanced ten fold and the emission maximum shifted from the 350-370 nm value in aqueous solution to 320 nm. It is predicted that similar differences might be expected for exposed and buried aminotyrosyl residues in a protein. 3. Exposed tyrosyl residues on the helical protein tropomyosin and a helical segment of paramyosin were aminated in part (39% and 34% of the total tyrosyl residues, respectively). The fluorescence of the aminated tyrosyl residues on these proteins was similar to that of the aminotyrosyl peptides in an aqueous medium. Although the fluorescence efficiency of an aminotyrosyl residue was much lower than that of a tyrosyl residue, it was easy to distinguish the fluorescence of the aminotyrosyl residues (350-355 nm) on the protein from that arising from unmodified tyrosyl residues (305 nm).     

A molecular portrait of the response to unfolded proteins

Using DNA microarrays, 381 genes have been found to be induced in response to unfolded proteins. The identity of the previously characterized 208 of these, and further experiments, have revealed new details on the scope of the unfolded protein response and its connection to the degradation of proteins at the endoplasmic reticulum.