Satisfaction of hydrogen-bonding potential influences the conservation of polar sidechains

Although polar amino acids tend to be found on the surface of proteins due to their hydrophilic nature, their important roles within the core of proteins are now becoming better recognized. It has long been understood that a significant number of mainchain functions will not achieve hydrogen bond satisfaction through the formation of secondary structures; in these circumstances, it is generally buried polar residues that provide hydrogen bond satisfaction. Here, we describe an analysis of the hydrogen-bonding of polar amino acids in a set of structurally aligned protein families. This allows us not only to calculate the conservation of each polar residue but also to assess whether conservation is correlated with the hydrogen-bonding potential of polar sidechains. We show that those polar sidechains whose hydrogen-bonding potential is satisfied tend to be more conserved than their unsatisfied or nonhydrogen-bonded counterparts, particularly when buried. Interestingly, these buried and satisfied polar residues are significantly more conserved than buried hydrophobic residues. Forming hydrogen bonds to mainchain amide atoms also influences conservation, with those satisfied buried polar residues that form two hydrogen bonds to mainchain amides being significantly more conserved than those that form only one or none. These results indicate that buried polar residues whose hydrogen-bonding potential is satisfied are likely to have important roles in maintaining protein structure.     

Tolerance to the substitution of buried apolar residues by charged residues in the homologous protein structures

Occurrence and accommodation of charged amino acid residues in proteins that are structurally equivalent to buried non-polar residues in homologues have been investigated. Using a dataset of 1,852 homologous pairs of crystal structures of proteins available at 2A or better resolution, 14,024 examples of apolar residues in the structurally conserved regions replaced by charged residues in homologues have been identified. Out of 2,530 cases of buried apolar residues, 1,677 of the equivalent charged residues in homologues are exposed and the rest of the charged residues are buried. These drastic substitutions are most often observed in homologous protein pairs with low sequence identity (<30%) and in large protein domains (>300 residues). Such buried charged residues in the large proteins are often located in the interface of sub-domains or in the interface of structural repeats, Beyond 7A of residue depth of buried apolar residues, or less than 4% of solvent accessibility, almost all the substituting charged residues are buried. It is also observed that acidic sidechains have higher preference to get buried than the positively charged residues. There is a preference for buried charged residues to get accommodated in the interior by forming hydrogen bonds with another sidechain than the main chain. The sidechains interacting with a buried charged residue are most often located in the structurally conserved regions of the alignment. About 50% of the observations involving hydrogen bond between buried charged sidechain and another sidechain correspond to salt bridges. Among the buried charged residues interacting with the main chain, positively charged sidechains form hydrogen bonds commonly with main chain carbonyls while the negatively charged residues are accommodated by hydrogen bonding with the main chain amides. These carbonyls and amides are usually located in the loops that are structurally variable among homologous proteins.     

Investigation of the PDZ domain ligand binding site using chemically modified peptides

Several chemically modified analogues to a tightly binding ligand for the second PDZ domain of MAGI-3 were synthesized and evaluated for their ability to compete with native peptide ligands. N-methyl scanning of the ligand backbone amides revealed the energetically important hydrogen bonds between the ligand backbone and the PDZ domain. Analogues to the ligand's conserved threonine/serine(-2) residue, involved in a side chain to side chain hydrogen bond with a conserved histidine in the PDZ domain, revealed that the interaction is highly sensitive to the steric structure around the hydroxyl group of this residue. Analogues of the ligand carboxy terminus revealed that the full hydrogen bond network of the GLGF loop is important in ligand binding.     

Patterns in ionizable side chain interactions in protein structures

In a selected set of 44 high-resolution, non-homologous protein structures, the intramolecular hydrogen bonds or salt bridges formed by ionizable amino acid side chains were identified and analyzed. The analysis was based on the investigation of several properties of the involved residues such as their solvent exposure, their belonging to a certain secondary structural element, and their position relative to the N- and C-termini of their respective structural element. It was observed that two-thirds of the interactions made by basic or acidic side chains are hydrogen bonds to polar uncharged groups. In particular, the majority (78%) of the hydrogen bonds between ionizable side chains and main chain polar groups (sch:mch bonds) involved at least one buried atom, and in 42% of the cases both interacting atoms were buried. In α-helices, the sch:mch bonds observed in the proximity of the C- and N-termini show a clear preference for acidic and basic side chains, respectively. This appears to be due to the partial charges of peptide group atoms at the termini of α-helices, which establish energetically favorable electrostatic interactions with side chain carrying opposite charge, at distances even greater than 4.5 Å. The sch:mch interactions involving ionizable side chains that belong either to β-strands or to the central part of α-helices are based almost exclusively on basic residues. This results from the presence of main chain carbonyl oxygen atoms in the protein core which have unsatisfied hydrogen bonding capabilities.     

Genetic selection reveals the role of a buried, conserved polar residue

The burial of nonpolar surface area is known to enhance markedly the conformational stability of proteins. The contribution from the burial of polar surface area is less clear. Here, we report on the tolerance to substitution of Ser75 of bovine pancreatic ribonuclease (RNase A), a residue that has the unusual attributes of being buried, conserved, and polar. To identify variants that retain biological function, we used a genetic selection based on the intrinsic cytotoxicity of ribonucleolytic activity. Cell growth at 30 degrees C, 37 degrees C, and 44 degrees C correlated with residue size, indicating that the primary attribute of Ser75 is its small size. The side-chain hydroxyl group of Ser75 forms a hydrogen bond with a main-chain nitrogen. The conformational stability of the S75A variant, which lacks this hydrogen bond, was diminished by DeltaDeltaG = 2.5 kcal/mol. Threonine, which can reinstate this hydrogen bond, provided a catalytically active RNase A variant at higher temperatures than did some smaller residues (including aspartate), indicating that a secondary attribute of Ser75 is the ability of its uncharged side chain to accept a hydrogen bond. These results provide insight on the imperatives for the conservation of a buried polar residue.     

Hydrogen bonds between short polar side chains and peptide backbone: prevalence in proteins and effects on helix-forming propensities

A survey of 322 proteins showed that the short polar (SP) side chains of four residues, Thr, Ser, Asp, and Asn, have a very strong tendency to form hydrogen bonds with neighboring backbone amides. Specifically, 32% of Thr, 29% of Ser, 26% of Asp, and 19% of Asn engage in such hydrogen bonds. When an SP residue caps the N terminal of a helix, the contribution to helix stability by a hydrogen bond with the amide of the N3 or N2 residue is well established. When an SP residue is in the middle of a helix, the side chain is unlikely to form hydrogen bonds with neighboring backbone amides for steric and geometric reasons. In essence the SP side chain competes with the backbone carbonyl for the same hydrogen-bonding partner (i.e., the backbone amide) and thus SP residues tend to break backbone carbonyl-amide hydrogen bonds. The proposition that this is the origin for the low propensities of SP residues in the middle of alpha helices (relative to those of nonpolar residues) was tested. The combined effects of restricting side-chain rotamer conformations (documented by Creamer and Rose, Proc Acad Sci USA, 1992;89:5937-5941; Proteins, 1994;19:85-97) and excluding side- chain to backbone hydrogen bonds by the helix were quantitatively analyzed. These were found to correlate strongly with four experimentally determined scales of helix-forming propensities. The correlation coefficients ranged from 0.72 to 0.87, which are comparable to those found for nonpolar residues (for which only the loss of side-chain conformational entropy needs to be considered).     

The transmembrane 7-alpha-bundle of rhodopsin: distance geometry calculations with hydrogen bonding constraints.

A 3D model of the transmembrane 7-alpha-bundle of rhodopsin-like G-protein-coupled receptors (GPCRs) was calculated using an iterative distance geometry refinement with an evolving system of hydrogen bonds, formed by intramembrane polar side chains in various proteins of the family and collectively applied as distance constraints. The alpha-bundle structure thus obtained provides H bonding of nearly all buried polar side chains simultaneously in the 410 GPCRs considered. Forty evolutionarily conserved GPCR residues form a single continuous domain, with an aliphatic "core" surrounded by six clusters of polar and aromatic side chains. The 7-alpha-bundle of a specific GPCR can be calculated using its own set of H bonds as distance constraints and the common "average" model to restrain positions of the helices. The bovine rhodopsin model thus determined is closely packed, but has a few small polar cavities, presumably filled by water, and has a binding pocket that is complementary to 11-cis (6-s-cis, 12-s-trans, C = N anti)-retinal or to all-trans-retinal, depending on conformations of the Lys296 and Trp265 side chains. A suggested mechanism of rhodopsin photoactivation, triggered by the cis-trans isomerization of retinal, involves rotations of Glu134, Tyr223, Trp265, Lys296, and Tyr306 side chains and rearrangement of their H bonds. The model is in agreement with published electron cryomicroscopy, mutagenesis, chemical modification, cross-linking, Fourier transform infrared spectroscopy, Raman spectroscopy, electron paramagnetic resonance spectroscopy, NMR, and optical spectroscopy data. The rhodopsin model and the published structure of bacteriorhodopsin have very similar retinal-binding pockets.     

The 2.3-A resolution structure of the maltose- or maltodextrin-binding protein, a primary receptor of bacterial active transport and chemotaxis

The three-dimensional structure of the maltose- or maltodextrin-binding protein (Mr = 40,622) with bound maltose has been obtained by crystallographic analysis at 2.8-A resolution. The structure, which has been partially refined at 2.3 A, is ellipsoidal with overall dimensions of 30 x 40 x 65 A and divided into two distinct globular domains by a deep groove. Although each domain is built from two peptide segments from the amino- and carboxyl-terminal halves, both domains exhibit similar supersecondary structure, consisting of a central beta-pleated sheet flanked on both sides with two or three parallel alpha-helices. The groove, which has a depth of 18 A and a base of about 9 x 18 A, contains the maltodextrin-binding site. We have previously observed the same general features in the well-refined structures of six other periplasmic receptors with specificities for L-arabinose, D-galactose/D-glucose, sulfate, phosphate, leucine/isoleucine/valine, and leucine. The bound maltose is buried in the groove and almost completely inaccessible to the bulk solvent. The groove is heavily populated by polar and aromatic groups many of which are involved in extensive hydrogen-bonding and van der Waals interactions with the maltose. All the disaccharide hydroxyl groups, which form a peripheral polar surface approximately in the plane of the sugar rings, are tied in a total of 11 direct hydrogen bonds with six charged side chains, one Trp side chain, and one peptide backbone NH, and five indirect hydrogen bonds via water molecules. The maltose is wedged between four aromatic side chains. The resulting stacking of these aromatic residues on the faces of the glucosyl units provides a majority of the van der Waals contacts in the complex. The nonreducing glucosyl unit of the maltose is involved in approximately twice as many hydrogen bonds and van der Waals contacts as the glucosyl unit at the reducing end. The binding protein-maltose complex shows the best example of the extensive use of polar and aromatic residues in binding oligosaccharides. The tertiary structure of the maltodextrin-binding protein, along with the results of genetic studies by a number of investigators, has also enabled us for the first time to map the different regions on the surface of the protein involved in the interactions with the membrane-bound protein components necessary for transport of and chemotaxis toward maltodextrins. These sites permit distinction of the "open cleft" (without bound sugar) and closed (with bound sugar) conformations of the binding protein by the chemotactic signal transducer with which the maltodextrin-binding protein interacts.(ABSTRACT TRUNCATED AT 250 WORDS)     

On the evolutionary conservation of hydrogen bonds made by buried polar amino acids: the hidden joists,braces and trusses of protein architecture

Background  

The hydrogen bond patterns between mainchain atoms in protein structures not only give rise to regular secondary structures but also satisfy mainchain hydrogen bond potential. However, not all mainchain atoms can be satisfied through hydrogen bond interactions that arise in regular secondary structures; in some locations sidechain-to-mainchain hydrogen bonds are required to provide polar group satisfaction. Buried polar residues that are hydrogen-bonded to mainchain amide atoms tend to be highly conserved within protein families, confirming that mainchain architecture is a critical restraint on the evolution of proteins. We have investigated the stabilizing roles of buried polar sidechains on the backbones of protein structures by performing an analysis of solvent inaccessible residues that are entirely conserved within protein families and superfamilies and hydrogen bonded to an equivalent mainchain atom in each family member.     

On the role of electrostatic interactions in the design of protein-protein interfaces

Here, the methods of continuum electrostatics are used to investigate the contribution of electrostatic interactions to the binding of four protein-protein complexes; barnase-barstar, human growth hormone and its receptor, subtype N9 influenza virus neuraminidase and the NC41 antibody, the Ras binding domain (RBD) of kinase cRaf and a Ras homologue Rap1A. In two of the four complexes electrostatics are found to strongly oppose binding (hormone-receptor and neuraminidase-antibody complexes), in one case the net effect is close to zero (barnase-barstar) and in one case electrostatics provides a significant driving force favoring binding (RBD-Rap1A). In order to help understand the wide range of electrostatic contributions that were calculated, the electrostatic free energy was partitioned into contributions of individual charged and polar residues, salt bridges and networks involving salt bridges and hydrogen bonds. Although there is no one structural feature that accounts for the differences between the four interfaces, the extent to which the desolvation of buried charges is compensated by the formation of hydrogen bonds and ion pairs appears to be an important factor. Structural features that are correlated with contribution of an individual residue to stability are also discussed. These include partial burial of a charged group in the free monomer, the formation of networks involving charged and polar amino acids, and the formation of partially exposed ion-pairs. The total electrostatic contribution to binding is found to be inversely correlated with buried total and non-polar surface area. This suggests that different interfaces can be designed to exploit electrostatic and hydrophobic forces in very different ways.     

Perturbing the polar environment of Asp102 in trypsin: consequences of replacing conserved Ser214.

Much of the catalytic power of trypsin is derived from the unusual buried, charged side chain of Asp102. A polar cave provides the stabilization for maintaining the buried charge, and it features the conserved amino acid Ser214 adjacent to Asp102. Ser214 has been replaced with Ala, Glu, and Lys in rat anionic trypsin, and the consequences of these changes have been determined. Three-dimensional structures of the Glu and Lys variant trypsins reveal that the new 214 side chains are buried. The 2.2-A crystal structure (R = 0.150) of trypsin S214K shows that Lys214 occupies the position held by Ser214 and a buried water molecule in the buried polar cave. Lys214-N zeta is solvent inaccessible and is less than 5 A from the catalytic Asp102. The side chain of Glu214 (2.8 A, R = 0.168) in trypsin S214E shows two conformations. In the major one, the Glu carboxylate in S214E forms a hydrogen bond with Asp102. Analytical isoelectrofocusing results show that trypsin S214K has a significantly different isoelectric point than trypsin, corresponding to an additional positive charge. The kinetic parameter kcat demonstrates that, compared to trypsin, S214K has 1% of the catalytic activity on a tripeptide amide substrate and S214E is 44% as active. Electrostatic potential calculations provide corroboration of the charge on Lys214 and are consistent with the kinetic results, suggesting that the presence of Lys214 has disturbed the electrostatic potential of Asp102.     

The 1.7 A crystal structure of BPI: a study of how two dissimilar amino acid sequences can adopt the same fold

We have extended the resolution of the crystal structure of human bactericidal/permeability-increasing protein (BPI) to 1.7 A. BPI has two domains with the same fold, but with little sequence similarity. To understand the similarity in structure of the two domains, we compare the corresponding residue positions in the two domains by the method of 3D-1D profiles. A 3D-1D profile is a string formed by assigning each position in the 3D structure to one of 18 environment classes. The environment classes are defined by the local secondary structure, the area of the residue which is buried from solvent, and the fraction of the area buried by polar atoms. A structural alignment between the two BPI domains was used to compare the 3D-1D environments of structurally equivalent positions. Greater than 31% of the aligned positions have conserved 3D-1D environments, but only 13% have conserved residue identities. Analysis of the 3D-1D environmentally conserved positions helps to identify pairs of residues likely to be important in conserving the fold, regardless of the residue similarity. We find examples of 3D-1D environmentally conserved positions with dissimilar residues which nevertheless play similar structural roles. To generalize our findings, we analyzed four other proteins with similar structures yet dissimilar sequences. Together, these examples show that aligned pairs of dissimilar residues often share similar structural roles, stabilizing dissimilar sequences in the same fold.     

Structural evidence for a dominant role of nonpolar interactions in the binding of a transport/chemosensory receptor to its highly polar ligands.

The receptor, a maltose/maltooligosaccharide-binding protein, has been found to be an excellent system for the study of molecular recognition because its polar and nonpolar binding functions are segregated into two globular domains. The X-ray structures of the "closed" and "open" forms of the protein complexed with maltose and maltotetraitol have been determined. These sugars have approximately 3 times more accessible polar surface (from OH groups) than nonpolar surface (from small clusters of sugar ring CH bonds). In the closed structures, the oligosaccharides are buried in the groove between the two domains of the protein and bound by extensive hydrogen bonding interactions of the OH groups with the polar residues confined mostly in one domain and by nonpolar interactions of the CH clusters with four aromatic residues lodged in the other domain. Substantial contacts between the sugar hydroxyls and aromatic residues are also formed. In the open structures, the oligosaccharides are bound almost exclusively in the domain rich in aromatic residues. This finding, along with the analysis of buried surface area due to complex formations in the open and closed structures, supports a major role for nonpolar interactions in initial ligand binding even when the ligands have significantly greater potential for highly specific polar interactions.     

Crystallographic structure at 1.6-A resolution of the human adenovirus proteinase in a covalent complex with its 11-amino-acid peptide cofactor: insights on a new fold

The crystal structure of the human adenovirus proteinase (AVP), a cysteine proteinase covalently bound to its 11-amino-acid peptide cofactor pVIc, has been solved to 1.6-A resolution with a crystallographic R-factor of 0.136, R(free)=0.179. The fold of AVP-pVIc is new and the structural basis for it is described in detail. The polypeptide chain of AVP folds into two domains. One domain contains a five-strand beta-sheet with two peripheral alpha-helices; this region represents the hydrophobic core of the protein. A second domain contains the N terminus, several C-terminal alpha-helices, and a small peripheral anti-parallel beta-sheet. The domains interact through an extended polar interface. pVIc spans the two domains like a strap, its C-terminal portion forming a sixth strand on the beta-sheet. The active site is in a long, deep groove located between the two domains. Portions are structurally similar to the active site of the prototypical cysteine proteinase papain, especially some of the Calpha backbone atoms (r.m.s. deviation of 0.354 A for 12 Calpha atoms). The active-site nucleophile of AVP, the conserved Cys(122), was shown to have a pK(a) of 4.5, close to the pK(a) of 3.0 for the nucleophile of papain, suggesting that a similar ion pair arrangement with His(54) may be present in AVP-pVIc. The interactions between AVP and pVIc include 24 non-beta-strand hydrogen bonds, six beta-strand hydrogen bonds and one covalent bond. Of the 204 amino acid residues in AVP, 33 are conserved among the many serotypes of adenovirus, and these aid in forming the active site groove, are involved in substrate specificity or interact between secondary structure elements.     

Completely buried, non-ion-paired glutamic acid contributes favorably to the conformational stability of pyrrolidone carboxyl peptidases from hyperthermophiles

Pyrrolidone carboxyl peptidases (PCPs) from hyperthermophiles have a structurally conserved and completely buried Glu192 in the hydrophobic core; in contrast, the corresponding residue in the mesophile protein is a hydrophobic residue, Ile. Does the buried ionizable residue contribute to stabilization or destabilization of hyperthermophile PCPs? To elucidate the role of the buried glutamic acid in stabilizing PCP from hyperthermophiles, we constructed five Glu192 mutants of PCP-0SH (C142S/C188S, Cys-free double mutant of PCP) from Pyrococcus furiosus and examined their thermal and pH-induced unfolding and crystal structures and compared them with those of PCP-0SH. The stabilities of apolar (E192A/I/V) and polar (E192D/Q) mutants were less than PCP-0SH at acidic pH values. In the alkaline region, the mutant proteins, except for E192D, were more stable than PCP-0SH. The thermal stability data and theoretical calculations indicated an apparent pKa value > or = 7.3 for Glu192. Present results confirmed that the protonated Glu192 in PCP-0SH forms strong hydrogen bonds with the carbonyl oxygen and peptide nitrogen of Pro168. New intermolecular hydrogen bonds in the E --> A/D mutants were formed by a water molecule introduced into the cavity created around position 192, whereas the hydrogen bonds disappeared in the E --> I/V mutants. Structure-based empirical stability of mutant proteins was in good agreement with the experimental results. The results indicated that (1) completely buried Glu192 contributes to the stabilization of PCP-0SH because of the formation of strong intramolecular hydrogen bonds and (2) the hydrogen bonds by the nonionized and buried Glu can contribute more than the burial of hydrophobic groups to the conformational stability of proteins.     

Conformational preferences and pK(a) value of selenocysteine residue

The conformational preferences of the L-selenocysteine (Sec) dipeptides with selenol and selenolate groups (Ac-Sec-NHMe and Ac-Sec(-) -NHMe, respectively) and the apparent (i.e., macroscopic) pK(a) value of the Sec residue have been studied using the dispersion-corrected density functionals M06-2X and B2PLYP-D with the implicit solvation method in the gas phase and in water. In the gas phase, the backbone-to-backbone and/or side chain-to-backbone hydrogen bonds are found to contribute in stabilizing the most preferred conformations for the Sec and Sec(-) residues, as seen for the Cys and Cys(-) residues. However, the polyproline II-like conformations prevail over the conformations with the backbone-to-backbone hydrogen bonds in water because of the weakened hydrogen bonds by the favorable direct interactions between the backbone C?O and H?N groups and water molecules. The Sec and Sec(-) residues are found to adopt more various conformations than the Cys and Cys(-) residues in water, although the most preferred conformations of the neutral and/or anionic forms of the two residues are similar each other in the gas phase and in water. Using the statistically weighted free energies of the Sec and Sec(-) dipeptides in the gas phase and their solvation free energies, the pK(a) value of the Sec residue is estimated to be 5.47 at 25°C, which is in good agreement with the experimental value of 5.43 ± 0.02. It is found that the lower pK(a) value of the selenol side chain for the Sec residue by ～3 units than the thiol side chain for the Cys residue is ascribed to the higher gas-phase acidity of the Sec residue.     

Structure of actinidin,after refinement at 1.7 Å resolution

The structure of the sulphydryl protease, actinidin, after refinement at 1.7 Å resolution, is described. The positions of most of the 1666 atoms have been determined with an accuracy better than 0.1 Å; only two residues (219 and 220) and the side-chain of a third (87) cannot be seen. In addition, the model contains 272 solvent molecules, all taken as water, except one which may be an ammonium ion. Atomic B values give a good indication of the mobility of different parts of the structure. Actinidin has a double domain structure, with one domain mostly helical in its secondary structure, and the other domain built around a twisted β-sheet. The geometry of hydrogen bonds in helices, β-structure and turns has been analysed. All are significantly non-linear, with the angle N-?…O ~160 °. Carbonyl groups are tilted outwards from the axis of each helix, the tilting apparently unaffected by whether or not additional hydrogen bonds are made (e.g. to water or side-chain atoms). Each domain is folded round a substantial core of non-polar side-chains, but the interface between domains is mostly polar. Interactions across this interface involve a network of eight buried water molecules, the buried carboxyl and amino groups of Glu35, Glu50, Lys181 and Lys17, other polar side-chains and a few hydrophobic groups. One other internal charged side-chain, that of Glu52, is adjacent to a buried solvent molecule, probably an ammonium ion. Other side-chain environments are described. One proline residue has a cis configuration. The sulphydryl group is oxidized, probably to SO₂^?, with one oxygen atom clearly visible but the other somewhat less certain. The active site geometry is otherwise compatible with the mechanism proposed by Drenth et al. (1975,1976) for papain. The positions of the 272 solvent molecules are described. The best-ordered water molecules are those that are internal (total of 17), in surface pockets, or in the intermolecular contact regions. These generally form three or four hydrogen bonds, two to proton acceptors and one or two to proton donors. Other water molecules make water bridges on the surface, sometimes covering the exposed edges of non-polar groups. Intermolecular contacts involve few protein atoms, but many water molecules.     

A method for determining the positions of polar hydrogens added to a protein structure that maximizes protein hydrogen bonding.

An automated method for the optimal placement of polar hydrogens in a protein structure is described. This method treats the polar, side chain hydrogens of lysine, serine, threonine, and tyrosine and the amino terminus of a protein. The program, called NETWORK, divides the potential hydrogen-bonding pairs of a protein into groups of interacting donors and acceptors. A search is conducted on each of the local groups to find an arrangement which forms the most hydrogen bonds. If two or more arrangements have the same number of hydrogen bonds, the arrangement with the shortest set of hydrogen bonds is selected. The polar hydrogens of the histidyl side chain are specifically treated, and the ionization state of this residue is allowed to change, if this change results in additional hydrogen bonds for the local group. The program will accept Protein Data Bank as well as Biosym-format coordinate files. Input and output routines can be easily modified to accept other coordinate file formats. The predictions from this method are compared to known hydrogen positions for bovine pancreatic trypsin inhibitor, insulin, RNase-A, and trypsin for which the neutron diffraction structures have been determined. The usefulness of this program is further demonstrated by a comparison of molecular dynamics simulations for the enzyme cytochrome P-450cam with and without using NETWORK.     

Structural and functional consequences of binding site mutations in transferrin: crystal structures of the Asp63Glu and Arg124Ala mutants of the N-lobe of human transferrin

Human transferrin is a serum protein whose function is to bind Fe(3+) with very high affinity and transport it to cells, for delivery by receptor-mediated endocytosis. Structurally, the transferrin molecule is folded into two globular lobes, representing its N-terminal and C-terminal halves, with each lobe possessing a high-affinity iron binding site, in a cleft between two domains. Central to function is a highly conserved set of iron ligands, including an aspartate residue (Asp63 in the N-lobe) that also hydrogen bonds between the two domains and an arginine residue (Arg124 in the N-lobe) that binds an iron-bound carbonate ion. To further probe the roles of these residues, we have determined the crystal structures of the D63E and R124A mutants of the N-terminal half-molecule of human transferrin. The structure of the D63E mutant, determined at 1.9 A resolution (R = 0.245, R(free) = 0.261), showed that the carboxyl group still binds to iron despite the larger size of the Glu side chain, with some slight rearrangement of the first turn of alpha-helix residues 63-72, to which it is attached. The structure of the R124A mutant, determined at 2.4 A resolution (R = 0.219, R(free) = 0.288), shows that the loss of the arginine side chain results in a 0.3 A displacement of the carbonate ion, and an accompanying movement of the iron atom. In both mutants, the iron coordination is changed slightly, the principal change being in each case a lengthening of the Fe-N(His249) bond. Both mutants also release iron more readily than the wild type, kinetically and in terms of acid lability of iron binding. We attribute this to more facile protonation of the synergistically bound carbonate ion, in the case of R124A, and to strain resulting from the accommodation of the larger Glu side chain, in the case of D63E. In both cases, the weakened Fe-N(His) bond may also contribute, consistent with protonation of the His ligand being an early intermediate step in iron release, following the protonation of the carbonate ion.     

Environment-specific amino acid substitution tables: tertiary templates and prediction of protein folds.

The local environment of an amino acid in a folded protein determines the acceptability of mutations at that position. In order to characterize and quantify these structural constraints, we have made a comparative analysis of families of homologous proteins. Residues in each structure are classified according to amino acid type, secondary structure, accessibility of the side chain, and existence of hydrogen bonds from the side chains. Analysis of the pattern of observed substitutions as a function of local environment shows that there are distinct patterns, especially for buried polar residues. The substitution data tables are available on diskette with Protein Science. Given the fold of a protein, one is able to predict sequences compatible with the fold (profiles or templates) and potentially to discriminate between a correctly folded and misfolded protein. Conversely, analysis of residue variation across a family of aligned sequences in terms of substitution profiles can allow prediction of secondary structure or tertiary environment.