Residues In Protein Structures Research Articles

Conserved residue clustering and protein structure prediction.

Protein residues that are critical for structure and function are expected to be conserved throughout evolution. Here, we investigate the extent to which these conserved residues are clustered in three-dimensional protein structures. In 92% of the proteins in a data set of 79 proteins, the most conserved positions in multiple sequence alignments are significantly more clustered than randomly selected sets of positions. The comparison to random subsets is not necessarily appropriate, however, because the signal could be the result of differences in the amino acid composition of sets of conserved residues compared to random subsets (hydrophobic residues tend to be close together in the protein core), or differences in sequence separation of the residues in the different sets. In order to overcome these limits, we compare the degree of clustering of the conserved positions on the native structure and on alternative conformations generated by the de novo structure prediction method Rosetta. For 65% of the 79 proteins, the conserved residues are significantly more clustered in the native structure than in the alternative conformations, indicating that the clustering of conserved residues in protein structures goes beyond that expected purely from sequence locality and composition effects. The differences in the spatial distribution of conserved residues can be utilized in de novo protein structure prediction: We find that for 79% of the proteins, selection of the Rosetta generated conformations with the greatest clustering of the conserved residues significantly enriches the fraction of close-to-native structures.

Geometry of interaction of the histidine ring with other planar and basic residues.

Among the aromatic residues in protein structures, histidine (His) is unique, as it can exist in the neutral or positively charged form at the physiological pH. As such, it can interact with other aromatic residues as well as form hydrogen bonds with polar and charged (both negative and positive) residues. We have analyzed the geometry of interaction of His residues with nine other planar side chains containing aromatic (residues Phe, Tyr, Trp, and His), carboxylate (Asp and Glu), carboxamide (Asn and Gln) and guanidinium (Arg) groups in 432 polypeptide chains. With the exception of the aspartic (Asp) and glutamic (Glu) acid side-chains, all other residues prefer to interact in a face-to-face or offset-face-stacked orientation with the His ring. Such a geometry is different from the edge-to-face relative orientation normally associated with the aromatic-aromatic interaction. His-His pair prefers to interact in a face-to-face orientation; however, when both the residues bind the same metal ion, the interplanar angle is close to 90 degrees. The occurrence of different interactions (including the nonconventional N-H...pi and C-H...pi hydrogen bonds) have been correlated with the relative orientations between the interacting residues. Several structural motifs, mostly involved in binding metal ions, have been identified by considering the cases where His residues are in contact with four other planar moieties. About 10% of His residues used here are also found in sequence patterns in PROSITE database. There are examples of the amino end of the Lys side chain interacting with His residues in such a way that it is located on an arc around a ring nitrogen atom.

Influence of Inter-Residue Interactions in Protein Structures

During the process of protein folding, the amino acid residues along the polypeptide chain interact with each other in a cooperative manner to form the stable native structure. The knowledge about interactions between amino acid residues in protein structures is very helpful to understand the mechanism of protein folding and stability. In this review, I classify the inter-residue interactions into short, medium and long range based on a simple geometric approach. The features of these interactions in different structural classes and folding types of globular proteins have been delineated. Further, the development of a new parameter, long-range order, based on inter-residue interactions and its application for predicting the protein folding rates are highlighted. The information gained from the studies on inter-residue interactions provides valuable insights for understanding protein folding and de novo protein design.

Protein packing: dependence on protein size, secondary structure and amino acid composition

We have used the occluded surface algorithm to estimate the packing of both buried and exposed amino acid residues in protein structures. This method works equally well for buried residues and solvent-exposed residues in contrast to the commonly used Voronoi method that works directly only on buried residues. The atomic packing of individual globular proteins may vary significantly from the average packing of a large data set of globular proteins. Here, we demonstrate that these variations in protein packing are due to a complex combination of protein size, secondary structure composition and amino acid composition. Differences in protein packing are conserved in protein families of similar structure despite significant sequence differences. This conclusion indicates that quality assessments of packing in protein structures should include a consideration of various parameters including the packing of known homologous proteins. Also, modeling of protein structures based on homologous templates should take into account the packing of the template protein structure.

Discovery of local packing motifs in protein structures

We present a language for describing structural patterns of residues in protein structures and a method for the discovery of such patterns that recur in a set of protein structures. The patterns impose restrictions on the spatial position of each residue, their order along the amino acid chain, and which amino acids are allowed in each position. Unlike other methods for comparing sets of protein structures, our method is not based on the use of pairwise structure comparisons which is often time consuming and can produce inconsistent results. Instead, the method simultaneously takes into account information from all structures in the search for conserved structure patterns which are potential structure motifs. The method is based on describing the spatial neighborhoods of each residue in each structure as a string and applying a sequence pattern discovery method to find patterns common to subsets of these strings. Finally it is checked whether the similarities between the neighborhood strings correspond to spatially similar substructures. We apply the method to analyze sets of very disparate proteins from the four different protein families: serine proteases, cuprodoxins, cysteine proteinases, and ferredoxins. The motifs found by the method correspond well to the site and motif information given in the annotation of these proteins in PDB, Swiss-Prot, and PROSITE. Furthermore, the motifs are confirmed by using the motif data to constrain the structural alignment of the proteins obtained with the program SAP. This gave the best superposition/alignment of the proteins given the motif assignment. Proteins 1999;34:206–219. © 1999 Wiley-Liss, Inc.

Discovery of local packing motifs in protein structures

We present a language for describing structural patterns of residues in protein structures and a method for the discovery of such patterns that recur in a set of protein structures. The patterns impose restrictions on the spatial position of each residue, their order along the amino acid chain, and which amino acids are allowed in each position. Unlike other methods for comparing sets of protein structures, our method is not based on the use of pairwise structure comparisons which is often time consuming and can produce inconsistent results. Instead, the method simultaneously takes into account information from all structures in the search for conserved structure patterns which are potential structure motifs. The method is based on describing the spatial neighborhoods of each residue in each structure as a string and applying a sequence pattern discovery method to find patterns common to subsets of these strings. Finally it is checked whether the similarities between the neighborhood strings correspond to spatially similar substructures. We apply the method to analyze sets of very disparate proteins from the four different protein families: serine proteases, cuprodoxins, cysteine proteinases, and ferredoxins. The motifs found by the method correspond well to the site and motif information given in the annotation of these proteins in PDB, Swiss-Prot, and PROSITE. Furthermore, the motifs are confirmed by using the motif data to constrain the structural alignment of the proteins obtained with the program SAP. This gave the best superposition/alignment of the proteins given the motif assignment.

Guanidinated casein hydrolysate stimulates pancreatic secretagogue release by direct action to the intestine in rats.

Previously, we observed that pancreatic exocrine secretion was strongly enhanced after duodenal administration of guanidinated casein peptic hydrolysate (HGC) in rats with chronic bile-pancreatic juice (BPJ) diversion. Using a perifusion system of isolated dispersed rat intestinal mucosal cells, we investigated whether this phenomenon depends on a direct effect of protein on the intestine to release pancreatic secretagogues, such as cholecystokinin (CCK). Amylase release from isolated pancreatic acini was used as an assay to measure CCK or other pancreatic secretagogues in the effluents. Mucosal effluent with HGC stimulated amylase release from acini with or without soybean trypsin inhibitor (SBTI). Perifused effluent with low-concentrated SBTI did not stimulate amylase release. These results indicate that HGC stimulates release of pancreatic secretagogue from the intestinal mucosal cells independent of remaining trypsin activity in the isolated mucosa. Effluents with intact casein, its peptic digest, and homoarginine, which is a unique amino acid contained in HGC, were unable to stimulate amylase release from acini. Effluent with a high concentration of SBTI, which is rich in arginine residues, stimulated amylase release, but not with the same tryptic inhibitory effect of lima bean trypsin inhibitor, which is poor in arginine residues. These findings suggest that guanidyl residues in protein structure are responsible for release of pancreatic secretagogues from isolated intestinal mucosal cells. Finally, the increment of amylase release from pancreatic acini in response to the perifused effluent with HGC was eliminated in acini treated with a potent CCK antagonist, FK480. We conclude that in rats with BPJ diversion, HGC stimulates CCK release from the intestine by direct action on intestinal mucosa.

Ideal architecture of residue packing and its observation in protein structures.

A simple model of sphere packing has been investigated as an ideal model for long-range interactions for the packing of non-bonded residues in protein structures. By superposing all residues, the geometry of packing around a central residue is investigated. It is found that all residues conform almost perfectly to this lattice model for sphere packing when a radius of 6.5 A is used to define non-bonded (virtual) interacting residues. Side-chain positions with respect to sequential backbone segments are relatively regular as well. This lattice can readily be used in conformation simulations to reduce the conformational space.

Structural studies by 1H NMR of a prototypic alpha-helical peptide (LYQELQKLTQTLK) and homologs in trifluoroethanol/water and on sodium dodecyl sulfate micelles.

The 1H NMR-determined structure and dynamics of a synthetic, amphiphilic alpha-helical peptide, PH-1.0 (LYQELQKLTQTLK), and several homologs were compared in 50% trifluoroethanol-d2 (TFE-d2)/H20 and in sodium dodecyl-d25 sulfate (SDS-d25) micelles. The peptides were designed to test the influence on secondary structure of placement of favored and disfavored residues relative to a "longitudinal, hydrophobic strip-of-helix" defined by the repeating leucines. PH-1.0 was highly ordered as an alpha-helix in 50% TFE-d2/H20 and in SDS-d25 micelles. Homologs PH-1.1, in which L1 was replaced by T, and PH-1,4, in which L12 was replaced by T. were found to be partially helical in both media. Calculated structures in SDS-d25 revealed that the helix of PH-1.1 was slightly disordered at the N-terminus, but that of PH-1.4 was completely disordered at the C-terminus. Examination of distributions of hydrophobic residues in protein structures revealed that, when [symbol: see text] = LIVFM and [symbol: see text] = nonLIVFM, the pattern [symbol: see text] is favored and [symbol: see text] is disfavored in alpha-helices. Several analogs of PH-1.0 incorporating these patterns were studied. Peptide PH-1.12 (LYQELQKLLQTLK) retained alpha-helical structure in both 50% TFE-d2/H20 and in SDS-d25 micelles. However, although PH-1.13 (LYQELQKLTLTLK) was fully helical in 50% TFE, it was helical only through residue 6 in SDS micelles. Two homologs containing an additional loop of the helix and repeats of favored (PH-5.0, NYLQTLLETLKTLLQK) or suppressed LL patterns (PH-5.11, NYLQTLETLKLTQK) gave similar results, i.e. the latter peptide was helical only through residue 6 in SDS micelles. The degree of local order in these SDS micelle-adsorbed peptides correlates to placement of hydrophobic residues in motifs which are favored or disfavored in proteins in general and in alpha-helices specifically.

The substitution of proline 35 by alanine in Rhodobacter capsulatus cytochrome c2 affects the overall protein stability but not the alkaline transition.

It was shown by Koshy et al. [1990, Proc. Natl Acad. Sci USA, 87, 8697-8701; 1994, Biochem. J., 299, 347-350] that the substitution of proline 30 by alanine (P30A) of Drosophila melanogaster and rat cytochromes c exhibited decreased stabilities in both the heme iron-methionine sulfur (Fe-S) bond and overall protein conformation. Now we have found that the stability properties of the equivalent mutant of Rhodobacter capsulatus cytochrome c2 (P35A) are somewhat different. Based on optical and NMR spectroscopies, the Rb.capsulatus P35A alkaline transition (pKalk) was found to be unchanged with respect to the wild type, suggesting that the mutation in Rb.capsulatus cytochrome c2 has little effect on the stability of the Fe-S bond. However, Rb.capsulatus conformational stability was found to be decreased by 1.6 kcal/mol in the oxidized state. The difference in the stability properties of the equivalent proline to alanine substitutions in various species underscores the importance of studying mutations in more than one species before drawing generalizations about the role of conserved residues in protein structure and function.

Measuring Residue Association in Protein Structures Possible Implications for Protein Folding

We propose a number of distance measures between residues in protein structures based on average, minimum and maximum distances of all atom (backbone and side-chain) coordinates or with respect to side-chain atom coordinates only. The d1-distance (D1-distance) refers to the average distance between side-chain (backbone and side-chain) atoms of a residue pair in a given structure. The dm-distance (Dm-distance) refers to the minimum distance between side-chain atoms (non-trivial minimum distance between all atoms of a residue pair). For each distance measure, averaging and normalizing over representative protein structures, association values and closeness orderings for all amino acid types are determined. The expected associations of side-chain interactions between oppositely charged residues, among hydrophobic residues and of cysteine with cysteine are confirmed. Several surprising associations are observed relative to (1) the aromatic residues tyrosine and tryptophan, but not phenylalanine; (2) multiple histidine residues; (3) asymmetries of arginine versus lysine, aspartate versus glutamate, alanine versus glycine, and asparagine versus glutamine; (4) absence of correlations of α-carbon distances with side-chain distances. The all atoms D1 -distance attractions are dominated by steric relationships, with glycine and alanine significantly close to all amino acids, whereas large residues are under-associated with all residue types. In contrast, for the closeness ordering corresponding to the minimum side-chain dm -distance, glycine and alanine are among the least associated. However, in the d1-distance alanine is significantly close to all hydrophobic residues with the exception of tryptophan. The dm-distance preferences display a pervasive attraction for tyrosine by almost all residue types, the prominence of tyrosine and tryptophan in cation-aromatic interactions, and the versatility of histidine in functionality. The principal findings suggest a new perspective on the early and intermediate stages of protein folding.

Polyphenols, astringency and proline-rich proteins

Recent, NMR and precipitation, studies of molecular recognition of proline-rich proteins and peptides by plant polyphenols are described and rationalized. The action of polysaccharides and caseins in the moderation of the astringent response, which is engendered by polyphenols present in foodstuffs and beverages, is described. The possible influence of plant cell wall glycoproteins on the process of lignification is discussed in the light of the observed affinity of phenolic substrates for prolyl residues in protein structures.

Geometry of interaction of metal ions with histidine residues in protein structures

An analysis of the geometry and the orientation of metal ions bound to histidine residues in proteins is presented. Cations are found to lie in the imidazole plane along the lone pair on the nitrogen atom. Out of the two tautomeric forms of the imidazole ring, the NE2-protonated form is normally preferred. However, when bound to a metal ion the ND1-protonated form is predominant and NE2 is the ligand atom. When the metal coordination is through ND1, steric interactions shift the side chain torsional angle, chi 2 from its preferred value of 90 or 270 degrees. The orientation of histidine residues is usually stabilized through hydrogen bonding; ND1-protonated form of a helical residue can form a hydrogen bond with the carbonyl oxygen atom in the preceding turn of the helix. A considerable number of ligands are found in helices and beta-sheets. A helical residue bound to a heme group is usually found near the C-terminus of the helix. Two ligand groups four residues apart in a helix, or two residues apart in a beta-strand are used in many proteins to bind metal ions.

Geometry of interaction of metal ions with sulfur-containing ligands in protein structures.

An analysis of the geometry of binding of metal ions by cysteine and methionine residues in protein structures has been made by using the Protein Data Bank. Metal ions have a distinct model of binding to each of these residues, and this is independent of the nature of the metal center or the type of protein. Metal ions tend to approach the sulfur of Met roughly 38 degrees from the perpendicular to the plane through atoms C gamma-S delta-C epsilon. For Cys, the approach direction is such that the M...S gamma-C beta-C alpha torsional angle is about +/- 90 or 180 degrees. The side-chain conformation of the cysteine residue is affected by the presence of the metal ion; there is a shift from the g+ conformation toward g- and mainly t conformations. When two Cys residues at positions i-3 and i bind to the same metal center, there appears to be some restriction on the geometry of metal binding by the residue i; for such a residue chi 1 and M...S gamma-C beta-C alpha angles are likely to be around 60 degrees and 270 degrees, respectively. Met and Cys residues coordinating to a metal ion are usually from coil or turn regions of the protein structure.

One- and two-dimensional NMR spectral analysis of the consequences of single amino acid replacements in proteins.

The traditional approach of using homologous sequences to elucidate the role of specific amino acid residues in protein structure and function becomes more meaningful as the number of differences is minimized, with the limit being alteration of a single residue. For small proteins in solution, NMR spectroscopy offers a means of obtaining detailed information about each residue and its response to a given change in the protein sequence. Extraction of this information has been aided by recent progress in spectrometer technology (higher magnetic fields, more sensitive signal detection, more sophisticated computers) and experimental strategies (new NMR pulse sequences including multiple-quantum and two-dimensional NMR methods). The set of avian ovomucoid third domains, which consists of the third domain proper plus a short leader (connecting peptide) and has a maximum of 56 amino acid residues, offers an attractive system for developing experimental methods for investigating sequence-structure and structure-function relationships in proteins. Our NMR results provide examples of sequence effects on pKa' values, average conformation, and internal motion of amino acid side chains.