Strain in protein structures as viewed through nonrotameric side chains: II. effects upon ligand binding.

The relation between the spatial positioning of nonrotameric residues and ligands was studied in 112 tertiary structures of protein-ligand complexes with a crystallographic resolution of 相似文献   

The substrate specificity of protein kinases which phosphorylate the alpha subunit of eukaryotic initiation factor 2.

The alpha subunit of eukaryotic protein synthesis initiation factor (eIF-2 alpha) is phosphorylated at a single serine residue (Ser51) by two distinct and well-characterized protein kinase, the haem-controlled repressor (HCR) and the double-stranded RNA-activated inhibitor (dsI). The sequence adjacent to Ser51 is rich in basic residues (Ser51-Arg-Arg-Arg-Ile-Arg) suggesting that they may be important in the substrate specificity of the two kinases, as is the case for several other protein kinases. A number of proteins and synthetic peptides containing clusters of basic residues were tested as substrates for HCR and dsI. Both kinases were able to phosphorylate histones and protamines ar multiple sites as judged by two-dimensional mapping of the tryptic phosphopeptides. These data also showed that the specificities of the two kinases were different from one another and from the specificities of two other protein kinases which recognise basic residues, cAMP-dependent protein kinase and protein kinase C. In histones, HCR phosphorylated only serine residues while dsI phosphorylated serine and threonine. Based on phosphoamino acid analyses and gel filtration of tryptic fragments, dsI was capable of phosphorylating both 'sites' in clupeine Y1 and salmine A1, whereas HCR acted only on the N-terminal cluster of serines in these protamines. The specificities of HCR and dsI were further studied using synthetic peptides with differing configurations of basic residues. Both kinases phosphorylated peptides containing C-terminal clusters of arginines on the 'target' serine residue, provided that they were present at positions +3 and/or +4 relative to Ser51. However, peptides containing only N-terminal basic residues were poor and very poor substrates for dsI and HCR, respectively. These findings are consistent with the disposition of basic residues near the phosphorylation site in eIF-2 alpha and show that the specificities of HCR and dsI differ from other protein kinases whose specificities have been studied.     

Linkers of secondary structures in proteins.

Linkers that connect repeating secondary structures fall into conformational classes based on distance and main-chain torsion clustering. A data set of 300 unique protein chains with low pairwise sequence identity was clustered into only a few groups representing the preferred motifs. The linkers of two to eight residues for the nonredundant data set are designated H-Ln-H, H-Ln-E, E-Ln-H, E-Ln-E, where n is the length, H stands for alpha-helices, and E for beta-strands. Most of the clusters identified here corroborate earlier findings. However, 19 new clusters are identified in this paper, with many of them having seven and eight residue linkers. In our first analysis, the secondary structures flanking the linkers are both interacting and noninteracting and there is no precise angle of orientation between them. A second analysis was performed on a set of proteins with restricted orientations for the flanking elements, namely, mainly alpha class of proteins with orthogonal architecture. Two definite clusters are identified, one corresponding to linkers of orthogonal helices and the other to linkers of antiparallel helices. Loops forming binding sites or involved in catalytic activity are important determinants of the function of proteins. Although the structural conservation of the residues around the catalytic triad of serine proteases has been studied widely, there has not been a systematic analysis of the conformation of the loops that contain them. Residues of the catalytic triad reside in the linkers of beta-strands, with varying lengths of more than eight residues. Here, we analyze the structural conservation of such linkers by superposition, and observe a conserved structural feature of the linkers incorporating each of the three residues of the catalytic triad.     

Automatic recognition of hydrophobic clusters and their correlation with protein folding units.

A method is described to objectively identify hydrophobic clusters in proteins of known structure. Clusters are found by examining a protein for compact groupings of side chains. Compact clusters contain seven or more residues, have an average of 65% hydrophobic residues, and usually occur in protein interiors. Although smaller clusters contain only side-chain moieties, larger clusters enclose significant portions of the peptide backbone in regular secondary structure. These clusters agree well with hydrophobic regions assigned by more intuitive methods and many larger clusters correlate with protein domains. These results are in striking contrast with the clustering algorithm of J. Heringa and P. Argos (1991, J Mol Biol 220:151-171). That method finds that clusters located on a protein's surface are not especially hydrophobic and average only 3-4 residues in size. Hydrophobic clusters can be correlated with experimental evidence on early folding intermediates. This correlation is optimized when clusters with less than nine hydrophobic residues are removed from the data set. This suggests that hydrophobic clusters are important in the folding process only if they have enough hydrophobic residues.     

The coordination sphere of iron-sulfur clusters: lessons from site-directed mutagenesis experiments

 Cysteine is the ubiquitous ligand of iron-sulfur clusters in proteins, although chemical models have indicated that functional groups other than thiolates can coordinate iron in iron-sulfur compounds. Only a small number of naturally occurring examples of hydroxyl, histidinyl or carboxyl coordination have been clearly established but many others are suspected. Quite a few site-directed mutagenesis experiments have been aimed at replacing the cysteine ligands of iron-sulfur centers by other amino acids in various systems. The available data set shows that substituting one ligand, even by another functional residue, is very often destabilizing enough to impair cluster assembly; in some cases, the apoprotein cannot even be detected. One for one replacements have been demonstrated, but they have been so far almost exclusively confined to clusters with no more than one or two iron atoms. In contrast, changes of the cluster nuclearity or recruitment of free cysteine residues seem preferred ways for proteins containing larger clusters to cope with removal of a ligand, rather than using coordinating amino acids bearing different chemical functions. Furthermore, the possibility of replacing cysteines by other residues as ligands in iron-sulfur proteins does not uniquely depend on the ability of the cluster to accept other kinds of coordination than cysteinate; other factors such as the local flexibility of the polypeptide chain, the accessibility of the solvent and the electronic distribution on the active centers may also play a prominent role.     

All-atom contact potential approach to protein thermostability analysis

In this paper we use all-atom potential energy to define and analyze the inter-residue contacts in mesophilic and thermophilic proteins. Fifteen families of proteins are selected and each family has two representative proteins with greatly different preferred environmental temperatures. We find that both the number and energy of the contacts defined in this way show stronger correlations with the preferred temperatures of proteins than other factors used before. We also find that the charged-polar and charged-nonpolar residue contacts not only have larger contact numbers but also have lower single contact energies. Furthermore, the most important is that most of the thermophilic proteins have more charged-polar and charged-nonpolar residue contacts than their mesophilic counterparts. This suggests that they may play an important role in the thermostability of proteins, except usual charged-charged and nonpolar-nonpolar residue contacts. Charged residues may exert their profound influence by forming contacts not only with other charged residues but also with polar or nonpolar residues, thus further increasing the strength of contact network and then the thermostability of proteins.     

Computation of non-covalent interactions in TNF proteins and interleukins

The roles played by the non-covalent interactions have been investigated for a set of six TNF proteins and nine Interleukins. The stabilizing residues have been identified by a consensus approach using the concepts of available surface area, medium and long-range interactions and conservation of amino acid residues. The cation-pi interactions have been computed based on a geometric approach such as distance and energy criteria. We identified an average of 1 energetically significant cation-pi interactions in every 94 residues in TNF proteins and 1 in every 62 residues in Interleukins. In TNF proteins, the cationic groups Lys preferred to be in helix while Arg preferred to be in strand regions while in Interleukins the Arg residues preferred to be in helix and Lys preferred to be in strand regions. From the available surface area calculations, we found that, almost all the cation and pi residues in TNF proteins and Interleukins were either in buried or partially buried regions and none of them in the exposed regions. Medium and long-range interactions were predominant in both TNF proteins and Interleukins. It was observed that the percentage of stabilizing centers were more in TNF proteins as compared to the Interleukins, while the percentage of conserved residues were more in Interleukins than in TNF proteins. In the stabilizing residues Lys was observed to be a stabilizing residue in both TNF proteins and Interleukins. Among the aromatic group, Phe was seen to be a stabilizing residue in both TNF and Interleukins. We suggest that this study on the computation of cation-pi interactions in TNF proteins and Interleukins would be very helpful in further understanding the structure, stability and functional similarity of these proteins.     

Alternative nucleophilic residues in intein catalysis of protein splicing

Protein-splicing inteins are widespread in nature and have found many applications in protein research and engineering. The mechanism of protein splicing typically requires a nucleophilic amino acid residue at both position 1 (first residue of intein) and position +1 (first residue after intein), however it was not clear whether or how the three different nucleophilic residues (Cys, Ser, and Thr) would work differently at these two positions. To use intein in a target protein of interest, one needs to choose an intein insertion site to have a nucleophilic residue at position +1, therefore it is desirable to know what nucleophilic residue(s) are preferred by different inteins. In this study we began with a statistical analysis of known inteins, which showed an unequal distribution of the three nucleophilic residues at positions 1 and +1, and then subjected six different mini-inteins to site-directed mutagenesis to systematically test the functionality of the three nucleophilic residues at the two positions. At position 1, most natural inteins had Cys and none had Thr. When the Cys at position 1 of the six inteins was mutated to Ser and Thr, the splicing activity was abolished in all except one case. At position +1, Cys and Ser were nearly equally abundant in natural inteins, and they were found to be functionally interchangeable in the six inteins of this study. When the two positions were studied as 1/+1 combination, the Cys/Ser combination was abundant in natural inteins, whereas the Ser/Cys combination was conspicuously absent. Similarly, all of the six inteins of this study spliced with the Cys/Ser combination, whereas none spliced with the Ser/Cys combination. These findings have interesting implications on the mechanism of splicing and the selection of intein insertion sites, and they also produced two rare mini-inteins that could splice with Thr at position +1.     

Analysis of homodimeric protein interfaces by graph-spectral methods

The quaternary structures impart structural and functional credibility to proteins. In a multi-subunit protein, it is important to understand the factors that drive the association or dissociation of the subunits. It is a well known fact that both hydrophobic and charged interactions contribute to the stability of the protein interface. The interface residues are also known to be highly conserved. Though they are buried in the oligomer, these residues are either exposed or partially exposed in the monomer. It is felt that a systematic and objective method of identifying interface clusters and their analysis can significantly contribute to the identification of a residue or a collection of residues important for oligomerization. Recently, we have applied the techniques of graph-spectral methods to a variety of problems related to protein structure and folding. A major advantage of this methodology is that the problem is viewed from a global protein topology point of view rather than localized regions of the protein structure. In the present investigation, we have applied the methods of graph-spectral analysis to identify side chain clusters at the interface and the centers of these clusters in a set of homodimeric proteins. These clusters are analyzed in terms of properties such as amino acid composition, accessibility to solvent and conservation of residues. Interesting results such as participation of charged and aromatic residues like arginine, glutamic acid, histidine, phenylalanine and tyrosine, consistent with earlier investigations, have emerged from these analyses. Important additional information is that the residues involved are a part of a cluster(s) and that they are sequentially distant residues which have come closer to each other in the three-dimensional structure of the protein. These residues can easily be detected using our graph-spectral algorithm. This method has also been used to identify important residues ('hot spots') in dimerization and also to detect dimerization sites on the monomer. The residues predicted using the present algorithm have correlated well with the experiments indicating the efficacy of this method in predicting residues involved in dimer stability.     

Recognizing protein folds by cluster distance geometry

Cluster distance geometry is a recent generalization of distance geometry whereby protein structures can be described at even lower levels of detail than one point per residue. With improvements in the clustering technique, protein conformations can be summarized in terms of alternative contact patterns between clusters, where each cluster contains four sequentially adjacent amino acid residues. A very simple potential function involving 210 adjustable parameters can be determined that favors the native contacts of 31 small, monomeric proteins over their respective sets of nonnative contacts. This potential then favors the native contacts for 174 small, monomeric proteins that have low sequence identity with any of the training set. A broader search finds 698 small protein chains from the Protein Data Bank where the native contacts are preferred over all alternatives, even though they have low sequence identity with the training set. This amounts to a highly predictive method for ab initio protein folding at low spatial resolution.     

Structural and energetic differences between insertions and substitutions in staphylococcal nuclease.

In a previous study, the small protein staphylococcal nuclease was shown to readily accommodate single alanine and glycine insertions, with average losses in stability comparable to substitutions at the same sites (PROT. 7:299-305, 1990). To more fully explore this unexpected adaptability to changes in residue spacing, 2 double amino acid insertions (alanyl-glycine, glycyl-glycine) and 3 additional single amino acid insertions with dissimilar side chains (proline, leucine, and glutamine) were constructed at 10 of the sites previously studied. At 8 of these sites, the type of amino acid side chain on the inserted residue significantly influenced the stability of the mutant protein. However, at 9 of the 10 sites, the double insertions were found to be no more destabilizing than the single alanine or glycine insertions. In contrast, double substitution mutations of staphylococcal nuclease, which replace two adjacent residues with alanine, do not show this striking degree of non-additivity. A comparison of the effects of single glutamine and single glycine insertions with alanyl-glycine insertions indicates that insertion of alanine into the peptide backbone is, on average, less destabilizing than appending the equivalent atoms onto the side chain of a glycine insertion. To explain their very different energetic effects, we propose that, unlike most substitutions, the inserted residue(s) must induce lateral displacements of the polypeptide chain, forcing the folded conformation away from that of wild type. The resulting obligatory shifts in the positioning of residues flanking the insertion generate a large number of degrees of freedom around which the mutant structure can relax.(ABSTRACT TRUNCATED AT 250 WORDS)     

Protein-DNA interactions: structural, thermodynamic and clustering patterns of conserved residues in DNA-binding proteins

Amino acid residues, which play important roles in protein function, are often conserved. Here, we analyze thermodynamic and structural data of protein-DNA interactions to explore a relationship between free energy, sequence conservation and structural cooperativity. We observe that the most stabilizing residues or putative hotspots are those which occur as clusters of conserved residues. The higher packing density of the clusters and available experimental thermodynamic data of mutations suggest cooperativity between conserved residues in the clusters. Conserved singlets contribute to the stability of protein-DNA complexes to a lesser extent. We also analyze structural features of conserved residues and their clusters and examine their role in identifying DNA-binding sites. We show that about half of the observed conserved residue clusters are in the interface with the DNA, which could be identified from their amino acid composition; whereas the remaining clusters are at the protein-protein or protein-ligand interface, or embedded in the structural scaffolds. In protein-protein interfaces, conserved residues are highly correlated with experimental residue hotspots, contributing dominantly and often cooperatively to the stability of protein-protein complexes. Overall, the conservation patterns of the stabilizing residues in DNA-binding proteins also highlight the significance of clustering as compared to single residue conservation.     

The four terminal components of the complement system are C-mannosylated on multiple tryptophan residues.

C-Mannosylation is a unique form of protein glycosylation, involving the C-glycosidic attachment of a mannosyl residue to the indole moiety of Trp. In the two examples found so far, human RNase 2 and interleukin-12, only the first Trp in the recognition motif WXXW is specifically C-mannosylated. To establish the generality of protein C-mannosylation, and to learn more about its mechanism, the terminal components of the human complement system (C6, C7, C8,and C9), which contain multiple and complex recognition motifs, were examined. Together with C5b they form the cytolytic agent, the membrane attack complex. These are the first proteins that are C-mannosylated on more than one Trp residue as follows: six in C6, four in C7, C8alpha, and C8beta, and two in C9. Thus, from the 113 Trp residues in the complete membrane attack complex, 50 were found to undergo C-mannosylation. The other important finding is that in C6, C7, C8, and C9 Trp residues without a second Trp (or another aromatic residue) at the +3 position can be C-mannosylated. This shows that they must contain an additional C-mannosylation signal. Whether this is encoded in the primary or tertiary structure is presently unknown. Finally, all modified Trp residues are part of the highly conserved core of the thrombospondin type 1 repeats present in these proteins. Since this module has been found in a large number of other proteins, the results suggest further candidates for C-mannosylation.     

Binding of basic peptides to acidic lipids in membranes: effects of inserting alanine(s) between the basic residues.

We studied the binding of peptides containing five basic residues to membranes containing acidic lipids. The peptides have five arginine or lysine residues and zero, one, or two alanines between the basic groups. The vesicles were formed from mixtures of a zwitterionic lipid, phosphatidylcholine, and an acidic lipid, either phosphatidylserine or phosphatidylglycerol. Measuring the binding using equilibrium dialysis, ultrafiltration, and electrophoretic mobility techniques, we found that all peptides bind to the membranes with a sigmoidal dependence on the mole fraction of acidic lipid. The sigmoidal dependence (Hill coefficient greater than 1 or apparent cooperativity) is due to both electrostatics and reduction of dimensionality and can be described by a simple model that combines Gouy-Chapman-Stern theory with mass action formalism. The adjustable parameter in this model is the microscopic association constant k between a basic residue and an acidic lipid (1 less than k less than 10 M-1). The addition of alanine residues decreases the affinity of the peptides for the membranes; two alanines inserted between the basic residues reduces k 2-fold. Equivalently, the affinity of the peptide for the membrane decreases 10-fold, probably due to a combination of local electrostatic effects and the increased loss of entropy that may occur when the more massive alanine-containing peptides bind to the membrane. The arginine peptides bind more strongly than the lysine peptides: k for an arginine residue is 2-fold higher than for a lysine residue. Our results imply that a cluster of arginine and lysine residues with interspersed electrically neutral amino acids can bind a significant fraction of a cytoplasmic protein to the plasma membrane if the cluster contains more than five basic residues.     

Comparison of protein-protein interactions in serine protease-inhibitor and antibody-antigen complexes: implications for the protein docking problem

The protein-protein interaction energy of 12 nonhomologous serine protease-inhibitor and 15 antibody-antigen complexes is calculated using a molecular mechanics formalism and dissected in terms of the main-chain vs. side-chain contribution, nonrotameric side-chain contributions, and amino acid residue type involvement in the interface interaction. There are major differences in the interactions of the two types of protein-protein complex. Protease-inhibitor complexes interact predominantly through a main-chain-main-chain mechanism while antibody-antigen complexes interact predominantly through a side-chain-side-chain or a side-chain-main-chain mechanism. However, there is no simple correlation between the main-chain-main-chain interaction energy and the percentage of main-chain surface area buried on binding. The interaction energy is equally effected by the presence of nonrotameric side-chain conformations, which constitute approximately 20% of the interaction energy. The ability to reproduce the interface interaction energy of the crystal structure if original side-chain conformations are removed from the calculation is much greater in the protease-inhibitor complexes than the antibody-antigen complexes. The success of a rotameric model for protein-protein docking appears dependent on the extent of the main-chain-main-chain contribution to binding. Analysis of (1) residue type and (2) residue pair interactions at the interface show that antibody-antigen interactions are very restricted with over 70% of the antibody energy attributable to just six residue types (Tyr > Asp > Asn > Ser > Glu > Trp) in agreement with previous studies on residue propensity. However, it is found here that 50% of the antigen energy is attributable to just four residue types (Arg = Lys > Asn > Asp). On average just 12 residue pair interactions (6%) contribute over 40% of the favorable interaction energy in the antibody-antigen complexes, with charge-charge and charge/polar-tyrosine interactions being prominent. In contrast protease inhibitors use a diverse set of residue types and residue pair interactions.     

The differential effect on two hybrid proteins of deletion mutations within the hydrophobic region of the Escherichia coli OmpA signal peptide

Oligonucleotide-directed site-specific mutagenesis was used to systematically shorten the hydrophobic region within the signal peptide of the Escherichia coli outer membrane protein OmpA. DNA encoding the wild type and mutant OmpA signal peptides were then fused in frame to DNA encoding the mature regions of Staphylococcus aureus nuclease A and TEM beta-lactamase. The ability of these signal peptides to direct processing of the resulting hybrid proteins was dependent on both their length and the protein to which they were fused. Deletion of two or more residues progressively slowed processing of pro-OmpA-nuclease. By contrast, pro-OmpA-beta-lactamase was less sensitive to the length of the hydrophobic region than to the nature of the deleted residue(s). Deletion of an Ala residue tended to reduce processing efficiency of pro-OmpA-beta-lactamase, while deletion of an Ile residue, together with the Ala residue, resulted in improvement. The loss of either 3 or 4 residues abolished processing of both hybrids. These data indicate that both the length as well as the identity of residues in the hydrophobic region are important. The relative importance of these two factors depends on the mature region of the protein being secreted.     

Three-dimensional view of the surface motif associated with the P-loop structure: cis and trans cases of convergent evolution

Here we identify the determinants of the nucleotide-binding ability associated with the P-loop-containing proteins, inferring their functional importance from their structural convergence to a unique three- dimensional (3D) motif. (1) A new surface 3D pattern is identified for the P-loop nucleotide-binding region, which is more selective than the corresponding sequence pattern; (2) the signature displays one residue that we propose is the determinant for the guanine-binding ability (the residues aligned to ras D119; this residue is known to be important only in the G-proteins, we extend the prediction to all the other P-loop- containing proteins); and (3) two cases of convergent evolution at the molecular level are highlighted in the analysis of the active site: the positive charge aligned to ras K117 and the arginine residues aligned to the GAP arginine finger.The analysis of the residues conserved on protein surfaces allows one to identify new functional or evolutionary relationships among protein structures that would not be detectable by conventional sequence or structure comparison methods.     

Three-dimensional cluster analysis identifies interfaces and functional residue clusters in proteins

Three-dimensional cluster analysis offers a method for the prediction of functional residue clusters in proteins. This method requires a representative structure and a multiple sequence alignment as input data. Individual residues are represented in terms of regional alignments that reflect both their structural environment and their evolutionary variation, as defined by the alignment of homologous sequences. From the overall (global) and the residue-specific (regional) alignments, we calculate the global and regional similarity matrices, containing scores for all pairwise sequence comparisons in the respective alignments. Comparing the matrices yields two scores for each residue. The regional conservation score (C(R)(x)) defines the conservation of each residue x and its neighbors in 3D space relative to the protein as a whole. The similarity deviation score (S(x)) detects residue clusters with sequence similarities that deviate from the similarities suggested by the full-length sequences. We evaluated 3D cluster analysis on a set of 35 families of proteins with available cocrystal structures, showing small ligand interfaces, nucleic acid interfaces and two types of protein-protein interfaces (transient and stable). We present two examples in detail: fructose-1,6-bisphosphate aldolase and the mitogen-activated protein kinase ERK2. We found that the regional conservation score (C(R)(x)) identifies functional residue clusters better than a scoring scheme that does not take 3D information into account. C(R)(x) is particularly useful for the prediction of poorly conserved, transient protein-protein interfaces. Many of the proteins studied contained residue clusters with elevated similarity deviation scores. These residue clusters correlate with specificity-conferring regions: 3D cluster analysis therefore represents an easily applied method for the prediction of functionally relevant spatial clusters of residues in proteins.     

Hot regions in protein--protein interactions: the organization and contribution of structurally conserved hot spot residues

Structurally conserved residues at protein-protein interfaces correlate with the experimental alanine-scanning hot spots. Here, we investigate the organization of these conserved, computational hot spots and their contribution to the stability of protein associations. We find that computational hot spots are not homogeneously distributed along the protein interfaces; rather they are clustered within locally tightly packed regions. Within the dense clusters, they form a network of interactions and consequently their contributions to the stability of the complex are cooperative; however the contributions of independent clusters are additive. This suggests that the binding free energy is not a simple summation of the single hot spot residue contributions. As expected, around the hot spot residues we observe moderately conserved residues, further highlighting the crucial role of the conserved interactions in the local densely packed environment. The conserved occurrence of these organizations suggests that they are advantageous for protein-protein associations. Interestingly, the total number of hydrogen bonds and salt bridges contributed by hot spots is as expected. Thus, H-bond forming residues may use a "hot spot for water exclusion" mechanism. Since conserved residues are located within highly packed regions, water molecules are easily removed upon binding, strengthening electrostatic contributions of charge-charge interactions. Hence, the picture that emerges is that protein-protein associations are optimized locally, with the clustered, networked, highly packed structurally conserved residues contributing dominantly and cooperatively to the stability of the complex. When addressing the crucial question of "what are the preferred ways of proteins to associate", these findings point toward a critical involvement of hot regions in protein-protein interactions.     

Determinants of strand register in antiparallel beta-sheets of proteins.

Antiparallel beta-sheets present two distinct environments to inter-strand residue pairs: beta(A,HB) sites have two backbone hydrogen bonds; whereas at beta(A,NHB) positions backbone hydrogen bonding is precluded. We used statistical methods to compare the frequencies of amino acid pairs at each site. Only approximately 10% of the 210 possible pairs showed occupancies that differed significantly between the two sites. Trends were clear in the preferred pairs, and these could be explained using stereochemical arguments. Cys-Cys, Aromatic-Pro, Thr-Thr, and Val-Val pairs all preferred the beta(A,NHB) site. In each case, the residues usually adopted sterically favored chi1 conformations, which facilitated intra-pair interactions: Cys-Cys pairs formed disulfide bonds; Thr-Thr pairs made hydrogen bonds; Aromatic-Pro and Val-Val pairs formed close van der Waals contacts. In contrast, to make intimate interactions at a beta(A,HB) site, one or both residues had to adopt less favored chi1 geometries. Nonetheless, pairs containing glycine and/or aromatic residues were favored at this site. Where glycine and aromatic side chains combined, the aromatic residue usually adopted the gauche conformation, which promoted novel aromatic ring-peptide interactions. This work provides rules that link protein sequence and tertiary structure, which will be useful in protein modeling, redesign, and de novo design. Our findings are discussed in light of previous analyses and experimental studies.