The interrelationships of side-chain and main-chain conformations in proteins

The accurate determination of a large number of protein structures by X-ray crystallography makes it possible to conduct a reliable statistical analysis of the distribution of the main-chain and side-chain conformational angles, how these are dependent on residue type, adjacent residue in the sequence, secondary structure, residue-residue interactions and location at the polypeptide chain termini. The interrelationship between the main-chain (phi, psi) and side-chain (chi 1) torsion angles leads to a classification of amino acid residues that simplify the folding alphabet considerably and can be a guide to the design of new proteins or mutational studies. Analyses of residues occurring with disallowed main-chain conformation or with multiple conformations shed some light on why some residues are less favoured in thermophiles.     

Electronic and vibrational circular dichroism of aromatic amino acids by density functional theory

Structures of model compounds mimicking aromatic amino acid residues in proteins are optimized by density functional theory (DFT), assuming that the main-chain conformation was a random coil. Excitation energies and dipole and rotational strengths for the optimized structures were calculated based on time-dependent DFT (TD-DFT). The electronic circular dichroism (ECD) bands of the models were significantly affected by side-chain conformations. Hydration models of the aromatic residues were also subjected to TD-DFT calculations, and the ECD bands of these models were found to be highly perturbed by the hydration of the main-chain amide groups. In addition to calculating the random-coil conformation, we also performed TD-DFT calculations of the aromatic residue models, assuming that the main-chain conformation was an alpha-helix or beta-strand. As expected, the overall feature of the ECD bands was also perturbed by the main-chain conformations. Moreover, vibrational circular dichroism (VCD) spectra of the hydration models in a random-coil structure were simulated by DFT, which showed that the VCD spectra are more sensitive to the side-chain conformations than the ECD spectra. The present results show that analyses combining ECD and VCD spectroscopy and using DFT calculations can elucidate the main- and side-chain conformations of aromatic residues in proteins.     

Validation of NMR side-chain conformations by packing calculations

The precision and accuracy of protein structures determined by nuclear magnetic resonance (NMR) spectroscopy depend on the completeness of input experimental data set. Typically, rather than a single structure, an ensemble of up to 20 equally representative conformers is generated and routinely deposited in the Protein Database. There are substantially more experimentally derived restraints available to define the main-chain coordinates than those of the side chains. Consequently, the side-chain conformations among the conformers are more variable and less well defined than those of the backbone. Even when a side chain is determined with high precision and is found to adopt very similar orientations among all the conformers in the ensemble, it is possible that its orientation might still be incorrect. Thus, it would be helpful if there were a method to assess independently the side-chain orientations determined by NMR. Recently, homology modeling by side-chain packing algorithms has been shown to be successful in predicting the side-chain conformations of the buried residues for a protein when the main-chain coordinates and sequence information are given. Since the main-chain coordinates determined by NMR are consistently more reliable than those of the side-chains, we have applied the side-chain packing algorithms to predict side-chain conformations that are compatible with the NMR-derived backbone. Using four test cases where the NMR solution structures and the X-ray crystal structure of the same protein are available, we demonstrate that the side-chain packing method can provide independent validation for the side-chain conformations of NMR structures. Comparison of the side-chain conformations derived by side-chain packing prediction and by NMR spectroscopy demonstrates that when there is agreement between the NMR model and the predicted model, on average 78% of the time the X-ray structure also concurs. While the side-chain packing method can confirm the reliable residue conformations in NMR models, more importantly, it can also identify the questionable residue conformations with an accuracy of 60%. This validation method can serve to increase the confidence level for potential users of structural models determined by NMR.     

Analysis of the relationship between side-chain conformation and secondary structure in globular proteins

The relationship between the preferred side-chain dihedral angles and the secondary structure of a residue was examined. The structures of 61 proteins solved to a resolution of 2.0 A (1 A = 0.1 nm) or better were analysed using a relational database to store the information. The strongest feature observed was that the chi 1 distribution for most side-chains in an alpha-helix showed an absence of the g- conformation and a shift towards the t conformation when compared to the non-alpha/beta structures. The exceptions to this tendency were for short polar side-chains that form hydrogen bonds with the main-chain which prefer g+. Shifts in the chi 1 preferences for residues in the beta-sheet were observed. Other side-chain dihedral angles (chi 2, chi 3, chi 4) were found to be influenced by the main-chain. This paper presents more accurate distributions for the side-chain dihedral angles which were obtained from the increased number of proteins determined to high resolution. The means and standard deviations for chi 1 and chi 2 angles are presented for all residues according to the secondary structure of the main-chain. The means and standard deviations are given for the most popular conformations for side-chains in which chi 3 and chi 4 rotations affect the position of C atoms.     

Reconstruction of protein conformations from estimated positions of the C alpha coordinates.

Protein C alpha coordinates are used to accurately reconstruct complete protein backbones and side-chain directions. This work employs potentials of mean force to align semirigid peptide groups around the axes that connect successive C alpha atoms. The algorithm works well for all residue types and secondary structure classes and is stable for imprecise C alpha coordinates. Tests on known protein structures show that root mean square errors in predicted main-chain and C beta coordinates are usually less than 0.3 A. These results are significantly more accurate than can be obtained from competing approaches, such as modeling of backbone conformations from structurally homologous fragments.     

Deterministic features of side-chain main-chain hydrogen bonds in globular protein structures

A total of 19 835 polar residues from a data set of 250 non-homologous and highly resolved protein crystal structures were used to identify side-chain main-chain (SC-MC) hydrogen bonds. The ratio of the number of SC-MC hydrogen bonds to the total number of polar residues is close to 1:2, indicating the ubiquitous nature of such hydrogen bonds. Close to 56% of the SC-MC hydrogen bonds are local involving side-chain acceptor/donor ('i') and a main-chain donor/acceptor within the window i-5 to i+5. These short-range hydrogen bonds form well defined conformational motifs characterized by specific combinations of backbone and side-chain torsion angles. (a) The Ser/Thr residues show the greatest preference in forming intra-helical hydrogen bonds between the atoms O(gamma)(i) and O(i-4). More than half the examples of such hydrogen bonds are found at the middle of alpha-helices rather than at their ends. The most favoured motif of these examples is alpha(R)alpha(R)alpha(R)alpha(R)(g(-)). (b) These residues also show great preference to form hydrogen bonds between O(gamma)(i) and O(i-3), which are closely related to the previous type and though intra-helical, these hydrogen bonds are more often found at the C-termini of helices than at the middle. The motif represented by alpha(R)alpha(R)alpha(R)alpha(R)(g(+)) is most preferred in these cases. (c) The Ser, Thr and Glu are the most frequently found residues participating in intra-residue hydrogen bonds (between the side-chain and main-chain of the same residue) which are characterized by specific motifs of the form beta(g(+)) for Ser/Thr residues and alpha(R)(g(-)g(+)t) for Glu/Gln. (d) The side-chain acceptor atoms of Asn/Asp and Ser/Thr residues show high preference to form hydrogen bonds with acceptors two residues ahead in the chain, which are characterized by the motifs beta (tt')alphaR and beta(t)alpha(R), respectively. These hydrogen bonded segments, referred to as Asx turns, are known to provide stability to type I and type I' beta-turns. (e) Ser/Thr residues often form a combination of SC-MC hydrogen bonds, with the side-chain donor hydrogen bonded to the carbonyl oxygen of its own peptide backbone and the side-chain acceptor hydrogen bonded to an amide hydrogen three residues ahead in the sequence. Such motifs are quite often seen at the beginning of alpha-helices, which are characterized by the beta(g(+))alpha(R)alpha(R) motif. A remarkable majority of all these hydrogen bonds are buried from the protein surface, away from the surrounding solvent. This strongly indicates the possibility of side-chains playing the role of the backbone, in the protein interiors, to satisfy the potential hydrogen bonding sites and maintaining the network of hydrogen bonds which is crucial to the structure of the protein.     

Role of amino acid residues at turns in the conformational stability and folding of human lysozyme

To clarify the role of amino acid residues at turns in the conformational stability and folding of a globular protein, six mutant human lysozymes deleted or substituted at turn structures were investigated by calorimetry, GuHCl denaturation experiments, and X-ray crystal analysis. The thermodynamic properties of the mutant and wild-type human lysozymes were compared and discussed on the basis of their three-dimensional structures. For the deletion mutants, Delta47-48 and Delta101, the deleted residues are in turns on the surface and are absent in human alpha-lactalbumin, which is homologous to human lysozyme in amino acid sequence and tertiary structure. The stability of both mutants would be expected to increase due to a decrease in conformational entropy in the denatured state; however, both proteins were destabilized. The destabilizations were mainly caused by the disappearance of intramolecular hydrogen bonds. Each part deleted was recovered by the turn region like the alpha-lactalbumin structure, but there were differences in the main-chain conformation of the turn between each deletion mutant and alpha-lactalbumin even if the loop length was the same. For the point mutants, R50G, Q58G, H78G, and G37Q, the main-chain conformations of these substitution residues located in turns adopt a left-handed helical region in the wild-type structure. It is thought that the left-handed non-Gly residue has unfavorable conformational energy compared to the left-handed Gly residue. Q58G was stabilized, but the others had little effect on the stability. The structural analysis revealed that the turns could rearrange the main-chain conformation to accommodate the left-handed non-Gly residues. The present results indicate that turn structures are able to change their main-chain conformations, depending upon the side-chain features of amino acid residues on the turns. Furthermore, stopped-flow GuHCl denaturation experiments on the six mutants were performed. The effects of mutations on unfolding-refolding kinetics were significantly different among the mutant proteins. The deletion/substitutions in turns located in the alpha-domain of human lysozyme affected the refolding rate, indicating the contribution of turn structures to the folding of a globular protein.     

Structural and Functional Restraints on the Occurrence of Single Amino Acid Variations in Human Proteins

Human genetic variation is the incarnation of diverse evolutionary history, which reflects both selectively advantageous and selectively neutral change. In this study, we catalogue structural and functional features of proteins that restrain genetic variation leading to single amino acid substitutions. Our variation dataset is divided into three categories: i) Mendelian disease-related variants, ii) neutral polymorphisms and iii) cancer somatic mutations. We characterize structural environments of the amino acid variants by the following properties: i) side-chain solvent accessibility, ii) main-chain secondary structure, and iii) hydrogen bonds from a side chain to a main chain or other side chains. To address functional restraints, amino acid substitutions in proteins are examined to see whether they are located at functionally important sites involved in protein-protein interactions, protein-ligand interactions or catalytic activity of enzymes. We also measure the likelihood of amino acid substitutions and the degree of residue conservation where variants occur. We show that various types of variants are under different degrees of structural and functional restraints, which affect their occurrence in human proteome.     

Properties of polyproline II, a secondary structure element implicated in protein-protein interactions

The polyproline II (PPII) conformation of protein backbone is an important secondary structure type. It is unusual in that, due to steric constraints, its main-chain hydrogen-bond donors and acceptors cannot easily be satisfied. It is unable to make local hydrogen bonds, in a manner similar to that of alpha-helices, and it cannot easily satisfy the hydrogen-bonding potential of neighboring residues in polyproline conformation in a manner analogous to beta-strands. Here we describe an analysis of polyproline conformations using the HOMSTRAD database of structurally aligned proteins. This allows us not only to determine amino acid propensities from a much larger database than previously but also to investigate conservation of amino acids in polyproline conformations, and the conservation of the conformation itself. Although proline is common in polyproline helices, helices without proline represent 46% of the total. No other amino acid appears to be greatly preferred; glycine and aromatic amino acids have low propensities for PPII. Accordingly, the hydrogen-bonding potential of PPII main-chain is mainly satisfied by water molecules and by other parts of the main-chain. Side-chain to main-chain interactions are mostly nonlocal. Interestingly, the increased number of nonsatisfied H-bond donors and acceptors (as compared with alpha-helices and beta-strands) makes PPII conformers well suited to take part in protein-protein interactions.     

Mutations and selection in the generation of class II histocompatibility antigen polymorphism.

A comparison of seven human DR and DC class II histocompatibility antigen beta-chain amino acid sequences indicates that the allelic variation is of comparable magnitude within the DR and DC beta-chain genes. Silent and replacement nucleotide substitutions in six DR and DC beta-chain sequences, as well as in seven murine class II sequences (three I-A beta and four I-A alpha alleles) were analyzed. The results suggest that the mutation rates are of a comparable magnitude in the nucleotide sequences encoding the first and second external domains of the class II molecules. Nevertheless, the allelic amino acid replacements are predominantly located in the first domains. We conclude that a conservative selective pressure acts on the second domains, whereas in many positions in the first domains replacement substitutions are selectively neutral or maybe even favoured. Thus, the difference between the first and second domains as regards the number of amino acid replacements is mainly due to selection.     

Crystal structure of an Fab fragment in complex with a meningococcal serosubtype antigen and a protein G domain.

Many pathogens present highly variable surface proteins to their host as a means of evading immune responses. The structure of a peptide antigen corresponding to the subtype P1.7 variant of the porin PorA from the human pathogen Neisseria meningitidis was determined by solution of the X-ray crystal structure of the ternary complex of the peptide (ANGGASGQVK) in complex with a Fab fragment and a domain from streptococcal protein G to 1.95 A resolution. The peptide adopted a beta-hairpin structure with a type I beta-turn between residues Gly4P and Gly7P, the conformation of the peptide being further stabilised by a pair of hydrogen bonds from the side-chain of Asn2P to main-chain atoms in Val9P. The antigen binding site within the Fab formed a distinct crevice lined by a high proportion of apolar amino acids. Recognition was supplemented by hydrogen bonds from heavy chain residues Thr50H, Asp95H, Leu97H and Tyr100H to main-chain and side-chain atoms in the peptide. Complementarity-determining region (CDR) 3 of the heavy chain was responsible for approximately 50 % of the buried surface area formed by peptide-Fab binding, with the remainder made up from CDRs 1 and 3 of the light chain and CDRs 1 and 2 of the heavy chain. Knowledge of the structures of variable surface antigens such as PorA is an essential prerequisite to a molecular understanding of antigenic variation and its implications for vaccine design.     

Canonical structures for the hypervariable regions of immunoglobulins

We have analysed the atomic structures of Fab and VL fragments of immunoglobulins to determine the relationship between their amino acid sequences and the three-dimensional structures of their antigen binding sites. We identify the relatively few residues that, through their packing, hydrogen bonding or the ability to assume unusual phi, psi or omega conformations, are primarily responsible for the main-chain conformations of the hypervariable regions. These residues are found to occur at sites within the hypervariable regions and in the conserved beta-sheet framework. Examination of the sequences of immunoglobulins of unknown structure shows that many have hypervariable regions that are similar in size to one of the known structures and contain identical residues at the sites responsible for the observed conformation. This implies that these hypervariable regions have conformations close to those in the known structures. For five of the hypervariable regions, the repertoire of conformations appears to be limited to a relatively small number of discrete structural classes. We call the commonly occurring main-chain conformations of the hypervariable regions "canonical structures". The accuracy of the analysis is being tested and refined by the prediction of immunoglobulin structures prior to their experimental determination.     

The 2 A resolution structure of the sulfate-binding protein involved in active transport in Salmonella typhimurium

The crystal structure of the liganded form of the sulfate-binding protein, an initial receptor for active transport of sulfate in Salmonella typhimurium, has been solved and refined at 2.0 A resolution (1 A = 0.1 nm). The final model, which consists of 2422 non-hydrogen atoms, one sulfate substrate and 143 water molecules, yields a crystallographic R-factor of 14.0% for 16,959 reflections between 8 and 2 A. The structure deviates from ideal bond lengths and angle distances by 0.015 A and 0.037 A, respectively. The protein is ellipsoid with overall dimensions of 35 A x 35 A x 65 A and consists of two similar globular domains. The two domains are linked by three distinct peptide segments, which though widely separated in the amino acid sequence, are in close proximity in the tertiary structure. As these connecting segments are located near the periphery of the molecule, they further serve as the base or a "boundary" of the deep cleft formed between the two domains. Despite the unusual interdomain connectivity, both domains have similar supersecondary structure consisting of a central five-stranded beta-pleated sheet sandwiched by alpha-helices on either side. The arrangement of the two domains gives rise to the ellipsoidal shape and to the cleft between the two domains wherein the sulfate substrate is found and completely engulfed. A discovery of considerable importance is that the sulfate substrate is tightly held in place primarily by seven hydrogen bonds, five of which are donated by main-chain peptide NH groups, another by a serine hydroxyl and the last by the indole NH moiety of a tryptophan side-chain; there are no positively charged residues, nor cations, nor water molecules within van der Waals' distance to the sulfate dianion. All the main-chain peptide units associated with the sulfate are in turn linked (via the peptide CO group) to arrays of hydrogen bonds. Three of these arrays are composed of alternating peptide units and hydrogen bonds within the solvent-exposed part of three alpha-helices and two are linked to a histidine and an arginine residue. The sulfate-binding protein bears strong similarity to the structures of four other periplasmic binding proteins solved in our laboratory which are specific for L-arabinose, D-galactose/D-glucose, leucine/isoleucine/valine and leucine. The similarity includes the ellipsoidal shape and the two globular domain structures, each domain consisting of a central beta-pleated sheet flanked by alpha-helices.(ABSTRACT TRUNCATED AT 400 WORDS)     

The reconstruction of side-chain conformations from protein Cα coordinates

A model of nine proteins including side-chain atoms have been built from the known C^α coordinates and amino acid sequences using a Monte Carlo Protein Building Annealing method. The Cartesian coordinates for the side-chain atoms were established with bond lengths and angles selected randomly from within previously determined ranges. A simulated annealing technique is used to generate some 300 structures with differing side-chain conformations. The atomic coordinates of the backbone atoms are fixed during the simulated annealing process. The coordinates of the side-chain atoms of 300 low energy conformations are averaged to obtain a mean structure that is minimized with the C^α atoms constrained to their position in the x-ray structure using the OPLS/AMBER force field with the GB/SA water model. The rms deviation of the main-chain atoms (without C^β) compared with the corresponding crystal structures is in the range 0.20–0.64 Å. The rms deviation of the side-chain atoms is between 1.72 and 2.71 Å and for all atoms is between 1.19 and 1.99 Å. The method is insensitive to random errors in the C^α positions and the computational requirement is modest. © 1997 John Wiley & Sons, Inc.     

Crystal structures of phenylalanyl-tRNA synthetase complexed with phenylalanine and a phenylalanyl-adenylate analogue

The crystal structures of Thermus thermophilus phenylalanyl-tRNA synthetase (PheRS) complexed with phenylalanine and phenylalaninyl-adenylate (PheOH-AMP), the synthetic analogue of phenylalanyl-adenylate, have been determined at 2.7A and 2.5A resolution, respectively. Both Phe and PheOH-AMP are engulfed in the active site cleft of the catalytic alpha-subunit of PheRS, and neither makes contact with the PheRS beta-subunit. The conformations and binding of Phe are almost identical in both complexes. The recognition of Phe by PheRS is achieved through a mixture of multiple van der Waals interactions and hydrogen bonds. The side-chain of the Phe substrate is sandwiched between the hydrophobic side-chains of Phealpha258 and Phealpha260 on one side, and the main-chain atoms of the two adjacent beta-strands on the other. The side-chains of Valalpha261 and Alaalpha314 form the back wall of the amino acid binding pocket. In addition, PheRS residues (Trpalpha149, Seralpha180, Hisalpha178, Argalpha204, Glnalpha218, and Glualpha220) form a total of seven hydrogen bonds with the main-chain atoms of Phe. The conformation of PheOH-AMP and the network of interactions of its AMP moiety with PheRS are reminiscent of the other class II synthetases. The structural similarity between PheRS and histidyl-tRNA synthetase extends to the amino acid binding site, which is normally unique for each enzyme. The complex structures suggest that the PheRS beta-subunit may affect the first step of the reaction (formation of phenylalanyl-adenylate) through the metal-mediated conserved alpha/beta-subunit interface. The modeling of tyrosine in the active site of PheRS revealed no apparent close contacts between tyrosine and the PheRS residues. This result implies that the proofreading mechanism against activated tyrosine, rather than direct recognition, may play the major role in the PheRS specificity.     

Nearest-neighbor effects and structural preferences in dipeptides are a function of the electronic properties of amino acid side-chains

The electronic properties of amino acid side-chains are emerging as an important factor in the preference for secondary structure in proteins. These properties have not been fully characterized, nor has their role in the behavior of peptides been explored in any detail. The present studies sought to evaluate several possibilities: 1) that hydrophilicity can be expressed solely in electronic terms, 2) that substituent effects of side-chains extend across the peptide bond, and (3) nearest-neighbor effects in dipeptides correlate with secondary structural preferences. Quantum mechanics (QM) calculations were used to define the electronic properties of individual amino acids and dipeptides. It was found that the hydrophilicity of an amino acid side-chain can be accurately represented as a function of the electron densities of its component atoms. In addition, the nature of an amino acid in the second position of a dipeptide affects the electronic properties (Mulliken populations and electron densities) of the main-chain atoms of the first residue. Certain electronic features of the dipeptides strongly correlated with propensity for secondary structure. Specifically, Mulliken population data at the Calpha atom and N atom predicted preference for alpha-helices versus coil and strand conformations, respectively. Analysis of dipeptides arrayed in either helical or extended structures revealed lengthening of main-chain bonds in the alpha-helical conformations. A thorough characterization of the electronic properties of amino acids and short peptide segments may provide a better understanding of the forces that determine secondary structure in proteins.     

Statistical theory for protein combinatorial libraries. Packing interactions, backbone flexibility, and the sequence variability of a main-chain structure

Combinatorial experiments provide new ways to probe the determinants of protein folding and to identify novel folding amino acid sequences. These types of experiments, however, are complicated both by enormous conformational complexity and by large numbers of possible sequences. Therefore, a quantitative computational theory would be helpful in designing and interpreting these types of experiment. Here, we present and apply a statistically based, computational approach for identifying the properties of sequences compatible with a given main-chain structure. Protein side-chain conformations are included in an atom-based fashion. Calculations are performed for a variety of similar backbone structures to identify sequence properties that are robust with respect to minor changes in main-chain structure. Rather than specific sequences, the method yields the likelihood of each of the amino acids at preselected positions in a given protein structure. The theory may be used to quantify the characteristics of sequence space for a chosen structure without explicitly tabulating sequences. To account for hydrophobic effects, we introduce an environmental energy that it is consistent with other simple hydrophobicity scales and show that it is effective for side-chain modeling. We apply the method to calculate the identity probabilities of selected positions of the immunoglobulin light chain-binding domain of protein L, for which many variant folding sequences are available. The calculations compare favorably with the experimentally observed identity probabilities.     

The structural basis of receptor-binding by Escherichia coli associated with diarrhea and septicemia

GafD in Escherichia coli G (F17) fimbriae is associated with diarrheal disease, and the structure of the ligand-binding domain, GafD1-178, has been determined at 1.7A resolution in the presence of the receptor sugar N-acetyl-D-glucosamine. The overall fold is a beta-barrel jelly-roll fold. The ligand-binding site was identified and localized to the side of the molecule. Receptor binding is mediated by side-chain as well main-chain interactions. Ala43-Asn44, Ser116-Thr117 form the sugar acetamide specificity pocket, while Asp88 confers tight binding and Trp109 appears to position the ligand. There is a disulfide bond that rigidifies the acetamide specificity pocket. The three fimbrial lectins, GafD, FimH and PapG share similar beta-barrel folds but display different ligand-binding regions and disulfide-bond patterns. We suggest an evolutionary path for the evolution of the very diverse fimbrial lectins from a common ancestral fold.     

Statistical and molecular dynamics studies of buried waters in globular proteins

Buried solvent molecules are common in the core of globular proteins and contribute to structural stability. Folding necessitates the burial of polar backbone atoms in the protein core, whose hydrogen-bonding capacities should be satisfied on average. Whereas the residues in alpha-helices and beta-sheets form systematic main-chain hydrogen bonds, the residues in turns, coils and loops often contain polar atoms that fail to form intramolecular hydrogen bonds. The statistical analysis of 842 high resolution protein structures shows that well-resolved, internal water molecules preferentially reside near residues without alpha-helical and beta-sheet secondary structures. These buried waters most often form primary hydrogen bonds to main-chain atoms not involved in intramolecular hydrogen bonds, providing strong evidence that hydrating main-chain atoms is a key structural role of buried water molecules. Additionally, the average B-factor of protein atoms hydrogen-bonded to waters is smaller than that of protein atoms forming intramolecular hydrogen bonds, and the average B-factor of water molecules involved in primary hydrogen bonds with main-chain atoms is smaller than the average B-factor of water molecules involved in secondary hydrogen bonds to protein atoms that form concurrent intramolecular hydrogen bonds. To study the structural coupling between internal waters and buried polar atoms in detail we simulated the dynamics of wild-type FKBP12, in which a buried water, Wat137, forms one side-chain and multiple main-chain hydrogen bonds. We mutated E60, whose side-chain hydrogen bonds with Wat137, to Q, N, S or A, to modulate the multiplicity and geometry of hydrogen bonds to the water. Mutating E60 to a residue that is unable to form a hydrogen bond with Wat137 results in reorientation of the water molecule and leads to a structural readjustment of residues that are both near and distant to the water. We predict that the E60A mutation will result in a significantly reduced affinity of FKBP12 for its ligand FK506. The propensity of internal waters to hydrogen bond to buried polar atoms suggests that ordered water molecules may constitute fundamental structural components of proteins, particularly in regions where alpha-helical or beta-sheet secondary structure is not present.