Binding of the protein Raf to the active form of Ras promotes activation of the MAP kinase signaling pathway, triggering cell growth and differentiation. Raf/Arg89 in the center of the binding interface plays an important role determining Ras–Raf binding affinity. We have investigated experimentally and computationally the Raf-R89K mutation, which abolishes signaling in vivo. The binding to [γ-35S]GTP-Ras of a fusion protein between the Raf-binding domain (RBD) of Raf and GST was reduced at least 175-fold by the mutation, corresponding to a standard binding free energy decrease of at least 3.0 kcal/mol. To compute this free energy and obtain insights into the microscopic interactions favoring binding, we performed alchemical simulations of the RBD, both complexed to Ras and free in solution, in which residue 89 is gradually mutated from Arg into Lys. The simulations give a standard binding free energy decrease of 2.9 ± 1.9 kcal/mol, in agreement with experiment. The use of numerous runs with three different force fields allows insights into the sources of uncertainty in the free energy and its components. The binding decreases partly because of a 7 kcal/mol higher cost to desolvate Lys upon binding, compared to Arg, due to better solvent interactions with the more concentrated Lys charge in the unbound state. This effect is expected to be general, contributing to the lower propensity of Lys to participate in protein–protein interfaces. Large contributions to the free energy change also arise from electrostatic interactions with groups up to 8 Å away, namely residues 37–41 in the conserved effector domain of Ras (including 4 kcal/mol from Ser39 which loses a bifurcated hydrogen bond to Arg89), the conserved Lys84 and Lys87 of Raf, and 2–3 specific water molecules. This analysis will provide insights into the large experimental database of Ras–Raf mutations.

The first crystal structure of an inorganic pyrophosphatase (S-PPase) from an archaebacterium, the thermophile Sulfolobus acidocaldarius, has been solved by molecular replacement and refined to an R-factor of 19.7% at 2.7 Å. S-PPase is a D3 homohexameric protein with one Mg2+ per active site in a position similar to, but not identical with, the first activating metal in mesophilic pyrophosphatases (PPase). In mesophilic PPases, Asp65, Asp70, and Asp102 coordinate the Mg2+, while only Asp65 and Asp102 do in S-PPase, and the Mg2+ moves by 0.7 Å. S-PPase may therefore be deactivated at low temperature by mispositioning a key metal ion.

The monomer S-PPase structure is very similar to that of Thermus thermophilus (T-PPase) and Escherichia coli (E-PPase), root-mean-square deviations around 1 Å/Cα. But the hexamer structures of S- and T-PPase are more tightly packed and more similar to each other than they are to that of E-PPase, as shown by the increase in surface area buried upon oligomerization. In T-PPase, Arg116 creates an interlocking ionic network to both twofold and threefold related monomers; S-PPase has hydrophilic interactions to threefold related monomers absent in both E- and T-PPase. In addition, the thermostable PPases have about 7% more hydrogen bonds per monomer than E-PPase, and, especially in S-PPase, additional ionic interactions anchor the C-terminus to the rest of the protein. Thermostability in PPases is thus due to subtle improvements in both monomer and oligomer interactions.

α1-Microglobulin (α1m) is an electrophoretically heterogeneous plasma protein. It belongs to the lipocalin superfamily, a group of proteins with a three-dimensional (3D) structure that forms an internal hydrophobic ligand-binding pocket. α1m carries a covalently linked unidentified chromophore that gives the protein a characteristic brown color and extremely heterogeneous optical properties. Twenty-one different colored tryptic peptides corresponding to residues 88–94, 118–121, and 122–134 of human α1m were purified. In these peptides, the side chains of Lys92, Lys118, and Lys130 carried size heterogeneous, covalently attached, unidentified chromophores with molecular masses between 122 and 282 atomic mass units (amu). In addition, a previously unknown uncolored lipophilic 282 amu compound was found strongly, but noncovalently associated with the colored peptides. Uncolored tryptic peptides containing the same Lys residues were also purified. These peptides did not carry any additional mass (i.e., chromophore) suggesting that only a fraction of the Lys92, Lys118, and Lys130 are modified. The results can explain the size, charge, and optical heterogeneity of α1m. A 3D model of α1m, based on the structure of rat epididymal retinoic acid-binding protein (ERABP), suggests that Lys92, Lys118, and Lys130 are semiburied near the entrance of the lipocalin pocket. This was supported by the fluorescence spectra of α1m under native and denatured conditions, which indicated that the chromophores are buried, or semiburied, in the interior of the protein. In human plasma, approximately 50% of α1m is complex bound to IgA. Only the free α1m carried colored groups, whereas α1m linked to IgA was uncolored.

We have recently constructed hyperactive human deoxyribonuclease I (DNase I) variants that digest double-stranded DNA more efficiently under physiological saline conditions by introducing positively charged amino acids at eight positions that can interact favorably with the negatively charged DNA phosphates. In this study, we present data from supercoiled DNA nicking, linear DNA digestion, and hyperchromicity assays that distinguish two classes of DNase I hyperactive variants based upon their activity dependence on Ca2+. Class A variants are highly dependent upon Ca2+, having up to 300-fold lower activity in the presence of Mg2+ alone compared to that in the presence of Mg2+ and Ca2+, and include Q9R, H44K, and T205K, in addition to wild-type DNase I. In contrast, the catalytic activity of Class B variants, which comprise the E13R, T14K, N74K, S75K, and N110R hyperactive variants, is relatively Ca2+ independent. A significant proportion of this difference in Ca2+-dependent activity can be attributed to one of the two structural calcium binding sites in DNase I. Compared to wild-type, the removal of Ca2+ binding site 2 by alanine replacements at Asp99, Asp107, and Glu112 decreased activity up to 26-fold in the presence of Mg2+ and Ca2+, but had no effect in the presence of Mg2+ alone. We propose that the rate-enhancing effect of Ca2+ binding at site 2 can be replaced by favorable electrostatic interactions created by proximal positively charged amino acid substitutions such as those found in the Class B variants, thus reducing the dependence on Ca2+.

The thioredoxin action upon the 2-oxoacid dehydrogenase complexes is investigated by using different thioredoxins, both wild-type and mutated. The attacking cysteine residue of thioredoxin is established to be essential for the thioredoxin-dependent activation of the complexes. Mutation of the buried cysteine residue to serine is not crucial for the activation, but prevents inhibition of the complexes, exhibited by the Clamydomonas reinhardtii thioredoxin m disulfide. Site-directed mutagenesis of D26, W31, F/W12, and Y/A70 (the Escherichia coli thioredoxin numbering is employed for all the thioredoxins studied) indicates that both the active site and remote residues of thioredoxin are involved in its interplay with the 2-oxoacid dehydrogenase complexes. Sequences of 11 thioredoxin species tested biochemically are aligned. The thioredoxin residues at the contact between the α3/310 and α1 helices, the length of the α1 helix and the charges in the α2-β3 and β4-β5 linkers are found to correlate with the protein influence on the 2-oxoacid dehydrogenase complexes (the secondary structural elements of thioredoxin are defined according to Eklund H et al., 1991, Proteins 11:13–28). The distribution of the charges on the surface of the thioredoxin molecules is analyzed. The analysis reveals the species specific polarization of the thioredoxin active site surroundings, which corresponds to the efficiency of the thioredoxin interplay with the 2-oxoacid dehydrogenase systems. The most effective mitochondrial thioredoxin is characterized by the strongest polarization of this area and the highest value of the electrostatic dipole vector of the molecule. Not only the magnitude, but also the orientation of the dipole vector show correlation with the thioredoxin action. The dipole direction is found to be significantly influenced by the charges of the residues 13/14, 51, and 83/85, which distinguish the activating and inhibiting thioredoxin disulfides.

Thymidylate synthase (TS) converts dUMP to dTMP by reductive methylation, where 5,10-methylenetetrahydrofolate is the source of both the methylene group and reducing equivalents. The mechanism of this reaction has been extensively studied, mainly using the enzyme from Escherichia coli. Bacillus subtilis contains two genes for TSs, ThyA and ThyB. The ThyB enzyme is very similar to other bacterial TSs, but the ThyA enzyme is quite different, both in sequence and activity. In ThyA TS, the active site histidine is replaced by valine. In addition, the B. subtilis enzyme has a 2.4-fold greater kcat than the E. coli enzyme. The structure of B. subtilis thymidylate synthase in a ternary complex with 5-fluoro-dUMP and 5,10-methylenetetrahydrofolate has been determined to 2.5 Å resolution. Overall, the structure of B. subtilis TS (ThyA) is similar to that of the E. coli enzyme. However, there are significant differences in the structures of two loops, the dimer interface and the details of the active site. The effects of the replacement of histidine by valine and a serine to glutamine substitution in the active site area, and the addition of a loop over the carboxy terminus may account for the differences in kcat found between the two enzymes.

We describe the results of a procedure for maximizing the number of sequences that can be reliably linked to a protein of known three-dimensional structure. Unlike other methods, which try to increase sensitivity through the use of fold recognition software, we only use conventional sequence alignment tools, but apply them in a manner that significantly increases the number of relationships detected. We analyzed 11 genomes and found that, depending on the genome, between 23 and 32% of the ORFs had significant matches to proteins of known structure. In all cases, the aligned region consisted of either >100 residues or >50% of the smaller sequence. Slightly higher percentages could be attained if smaller motifs were also included. This is significantly higher than most previously reported methods, even those that have a fold-recognition component. We survey the biochemical and structural characteristics of the most frequently occurring proteins, and discuss the extent to which alignment methods can realistically assign function to gene products.

Effective inhibitors of matrix metalloproteinases (MMPs), a family of connective tissue-degrading enzymes, could be useful for the treatment of diseases such as cancer, multiple sclerosis, and arthritis. Many of the known MMP inhibitors are derived from peptide substrates, with high potency in vitro but little selectivity among MMPs and poor bioavailability. We have discovered nonpeptidic MMP inhibitors with improved properties, and report here the crystal structures of human stromelysin-1 catalytic domain (SCD) complexed with four of these inhibitors. The structures were determined and refined at resolutions ranging from 1.64 to 2.0 Å. Each inhibitor binds in the active site of SCD such that a bulky diphenyl piperidine moiety penetrates a deep, predominantly hydrophobic S′1 pocket. The active site structure of the SCD is similar in all four inhibitor complexes, but differs substantially from the peptide hydroxamate complex, which has a smaller side chain bound in the S′1 pocket. The largest differences occur in the loop forming the “top” of this pocket. The occupation of these nonpeptidic inhibitors in the S′1 pocket provides a structural basis to explain their selectivity among MMPs. An analysis of the unique binding mode predicts structural modifications to design improved MMP inhibitors.

Recent studies have pointed out the important role of local water structures in protein conformational stability. Here, we present an accurate and computationally effective way to estimate the free energy contribution of the simplest water structure motif—the water bridge. Based on the combination of empirical parameters for accessible protein surface area and the explicit consideration of all possible water bridges with the protein, we introduce an improved protein solvation model. We find that accounting for water bridge formation in our model is essential to understand the conformational behavior of polypeptides in water. The model formulation, in fact, does not depend on the polypeptide nature of the solute and is therefore applicable to other flexible biomolecules (i.e., DNAs, RNAs, polysaccharides, etc.).

A statistical analysis of a representative data set of 169 known protein structures was used to analyze the specificity of residue interactions between spatial neighboring strands in β-sheets. Pairwise potentials were derived from the frequency of residue pairs in nearest contact, second nearest and third nearest contacts across neighboring β-strands compared to the expected frequency of residue pairs in a random model. A pseudo-energy function based on these statistical pairwise potentials recognized native β-sheets among possible alternative pairings. The native pairing was found within the three lowest energies in 73% of the cases in the training data set and in 63% of β-sheets in a test data set of 67 proteins, which were not part of the training set. The energy function was also used to detect tripeptides, which occur frequently in β-sheets of native proteins. The majority of native partners of tripeptides were distributed in a low energy range. Self-correcting distance geometry (SECODG) calculations using distance constraints sets derived from possible low energy pairing of β-strands uniquely identified the native pairing of the β-sheet in pancreatic trypsin inhibitor (BPTI). These results will be useful for predicting the structure of proteins from their amino acid sequence as well as for the design of proteins containing β-sheets.

The moderate thermophilic eubacterium Alicyclobacillus (formerly Bacillus) acidocaldarius expresses a thermostable carboxylesterase (esterase 2) belonging to the hormone-sensitive lipase (HSL)-like group of the esterase/lipase family. Based on secondary structures predictions and a secondary structure-driven multiple sequence alignment with remote homologous protein of known three-dimensional (3D) structure, we previously hypothesized for this enzyme the α/β-hydrolase fold typical of several lipases and esterases and identified Ser155, Asp252, and His282 as the putative members of the catalytic triad. In this paper we report the construction of a 3D model for this enzyme based on the structure of mouse acetylcholinesterase complexed with fasciculin. The model reveals the topological organization of the fold corroborating our predictions. As regarding the active-site residues, Ser155, Asp252, and His282 are located close to each other at hydrogen bond distances. Their catalytic role was here probed by biochemical and mutagenic studies. Moreover, on the basis of the secondary structure-driven multiple sequence alignment and the 3D structural model, a residue supposed important for catalysis, Gly84, was mutated to Ser. The activity of the mutated enzyme was drastically reduced. We propose that Gly84 is part of a putative “oxyanion hole” involved in the stabilization of the transition state similar to the C group of the esterase/lipase family.

Cytochrome cH is the electron donor to the oxidase in methylotrophic bacteria. Its amino acid sequence suggests that it is a typical Class I cytochrome c, but some features of the sequence indicated that its structure might be of special interest. The structure of oxidized cytochrome cH has been solved to 2.0 Å resolution by X-ray diffraction. It has the classical tertiary structure of the Class 1 cytochromes c but bears a closer gross resemblance to mitochondrial cytochrome c than to the bacterial cytochrome c2. The left-hand side of the haem cleft is unique; in particular, it is highly hydrophobic, the usual water is absent, and the “conserved” Tyr67 is replaced by tryptophan. A number of features of the structure demonstrate that the usual hydrogen bonding network involving water in the haem channel is not essential and that other mechanisms may exist for modulation of redox potentials in this cytochrome.

The binding free energies of four inhibitors to bovine β-trypsin are calculated. The inhibitors use either ornithine, lysine, or arginine to bind to the S1 specificity site. The electrostatic contribution to binding free energy is calculated by solving the finite difference Poisson–Boltzmann equation, the contribution of nonpolar interactions is calculated using a free energy-surface area relationship and the loss of conformational entropy is estimated both for trypsin and ligand side chains. Binding free energy values are of a reasonable magnitude and the relative affinity of the four inhibitors for trypsin is correctly predicted. Electrostatic interactions are found to oppose binding in all cases. However, in the case of ornithine- and lysine-based inhibitors, the salt bridge formed between their charged group and the partially buried carboxylate of Asp189 is found to stabilize the complex. Our analysis reveals how the molecular architecture of the trypsin binding site results in highly specific recognition of substrates and inhibitors. Specifically, partially burying Asp189 in the inhibitor-free enzyme decreases the penalty for desolvation of this group upon complexation. Water molecules trapped in the binding interface further stabilize the buried ion pair, resulting in a favorable electrostatic contribution of the ion pair formed with ornithine and lysine side chains. Moreover, all side chains that form the trypsin specificity site are partially buried, and hence, relatively immobile in the inhibitor-free state, thus reducing the entropic cost of complexation. The implications of the results for the general problem of recognition and binding are considered. A novel finding in this regard is that like charged molecules can have electrostatic contributions to binding that are more favorable than oppositely charged molecules due to enhanced interactions with the solvent in the highly charged complex that is formed.

The partly folded states of α-lactalbumin (α-LA) exposed to acid solution at pH 2.0 (A-state) or at neutral pH upon EDTA-mediated removal of the single protein-bound calcium ion (apo form) have been probed by limited proteolysis experiments. These states are nowadays commonly considered to be molten globules and thus protein-folding intermediates. Pepsin was used for proteolysis at acid pH, while proteinase K and chymotrypsin at neutral pH. The expectations were that these proteolytic probes would detect sites and/or chain regions in the partly folded states of α-LA sufficiently dynamic, or even unfolded, capable of binding and adaptation to the specific stereochemistry of the protease's active site. A time-course analysis of the proteolytic events revealed that the fast, initial proteolytic cuts of the 123-residue chain of α-LA in its A-state or apo form by the three proteases occur at the same chain region 39–54, the actual site(s) of cleavage depending upon the protease employed. This region in native α-LA encompasses the β-sheets of the protein. Subsequent cleavages occur mostly at chain regions 31–35 and 95–105. Four fragment species of α-LA have been isolated by reverse-phase high-performance liquid chromatography, and their conformational properties examined by circular dichroism and fluorescence emission spectroscopy. The single chain fragment 53–103, containing all the binding sites for calcium in native α-LA and cross-linked by two disulfide bridges, maintains in aqueous buffer and in the presence of calcium ions a folded structure characterized by the same content of α-helix of the corresponding chain segment in native α-LA. Evidence for some structure was also obtained for the two-chain species 1–40 and 104–123, as well as 1–31 and 105–123, both systems being covalently linked by two disulfide bonds. In contrast, the protein species given by fragment 1–34 connected to fragment 54–123 or 57–123 via four disulfide bridges adopts in solution a folded structure with the helical content expected for a native-like conformation. Of interest, the proteolytic fragment species herewith isolated correspond to the structural domains and subdomains of α-LA that can be identified by computational analysis of the three-dimensional structure of native α-LA (Siddiqui AS, Barton GI, 1995, Protein Sci 4:872–884). The fast, initial cleavages at the level of the β-sheet region of native α-LA indicate that this region is highly mobile or even unfolded in the α-LA molten globule(s), while the rest of the protein chain maintains sufficient structure and rigidity to prevent extensive proteolysis. The subsequent cleavages at chain segment 95–105 indicate that also this region is somewhat mobile in the A-state or apo form of the protein. It is concluded that the overall domain topology of native α-LA is maintained in acid or at neutral pH upon calcium depletion. Moreover, the molecular properties of the partly folded states of α-LA deduced here from proteolysis experiments do correlate with those derived from previous NMR and other physicochemical measurements.

The structure and dynamics of the fatty acid binding cavity in I-FABP (rat intestinal fatty acid binding protein) were analyzed. In the crystal structure of apo I-FABP, the probe occupied cavity volume and surface are 539 ± 8 Å3 and 428 Å2, respectively (1.4 Å probe). A total of 31 residues contact the cavity with their side chains. The side-chain cavity surface is partitioned according to the residue type as follows: 36–39% hydrophobic, 21–25% hydrophilic, and 37–43% neutral or ambivalent. Thus, the cavity surface is neither like a typical protein interior core, nor is like a typical protein external surface. All hydrophilic residues that contact the cavity—with the exception of Asp74—are clustered on the one side of the cavity. The cavity appears to expand its hydrophobic surface upon fatty acid binding on the side opposite to this hydrophilic patch. In holo I-FABP the fatty acid chain interactions with the hydrophilic side chains are mediated by water molecules. Molecular dynamics (MD) simulation of fully solvated apo I-FABP showed global conformational changes of I-FABP, which resulted in a large, but seemingly transient, exposure of the cavity to the external solvent. The packing density of the side chains lining the cavity, studied by Voronoi volumes, showed the presence of two distinctive small hydrophobic cores. The MD simulation predicts significant structural perturbations of the cavity on the subnanosecond time scale, which are capable of facilitating exchange of I-FABP internal water.

Human growth hormone (hGH) binds to its receptor (hGHr) in a three-body interaction: one molecule of the hormone and two identical monomers of the receptor form a trimer. Curiously, the hormone-receptor interactions in the trimer are not equivalent and the formation of the complex occurs in a specific kinetic order (Cunningham BC, Ultsch M, De Vos AM, Mulkerrin MG, Clauser KR, Wells JA, 1991, Science 254:821–825). In this paper, we model the recognition of hGH to the hGHr using shape complementarity of the three-dimensional structures and macromolecular docking to explore possible binding modes between the receptor and hormone. The method, reported previously (Hendrix DK, Kuntz ID, 1998, Pacific symposium on biocomputing 1998, pp 1234–1244), is based upon matching complementary-shaped strategic sites on the molecular surface. We modify the procedure to examine three-body systems. We find that the order of binding seen experimentally is also essential to our model. We explore the use of mutational data available for hGH to guide our model. In addition to docking hGH to the hGHr, we further test our methodology by successfully reproducing 16 macromolecular complexes from X-ray crystal structures, including enzyme-inhibitor, antibody-antigen, protein dimer, and protein-DNA complexes.

Lipidation catalyzed by protein prenyltransferases is essential for the biological function of a number of eukaryotic proteins, many of which are involved in signal transduction and vesicular traffic regulation. Sequence similarity searches reveal that the α-subunit of protein prenyltransferases (PTα) is a member of the tetratricopeptide repeat (TPR) superfamily. This finding makes the three-dimensional structure of the rat protein farnesyltransferase the first structural model of a TPR protein interacting with its protein partner. Structural comparison of the two TPR domains in protein farnesyltransferase and protein phosphatase 5 indicates that variation in TPR consensus residues may affect protein binding specificity through altering the overall shape of the TPR superhelix. A general approach to evolutionary analysis of proteins with repetitive sequence motifs has been developed and applied to the protein prenyltransferases and other TPR proteins. The results suggest that all members in PTα family originated from a common multirepeat ancestor, while the common ancestor of PTα and other members of TPR superfamily is likely to be a single repeat protein.

The heat stable inhibitor of cAMP-dependent protein kinase (PKIα) contains both a nuclear export signal (NES) and a high affinity inhibitory region that is essential for inhibition of the catalytic subunit of the kinase. These functions are sequentially independent. Two-dimensional NMR spectroscopy was performed on uniformly [15N]-labeled PKIα to examine its structure free in solution. Seventy out of 75 residues were identified, and examination of the CαH chemical shifts revealed two regions of upfield chemical shifts characteristic of α-helices. When PKIα was fragmented into two functionally distinct peptides for study at higher concentrations, no significant alterations in chemical shifts or secondary structure were observed. The first ordered region, identified in PKIα (1–25), contains an α-helix from residues 1–13. This helix extends by one turn the helix observed in the crystal structure of a PKIα (5–24) peptide bound to the catalytic subunit. The second region of well-defined secondary structure, residues 35–47, overlaps with the nuclear export signal in the PKIα (26–75) fragment. This secondary structure consists of a helix with a hydrophobic face comprised of Leu37, Leu41, and Leu44, followed by a flexible turn containing Ile46. These four residues are critical for nuclear export function. The remainder of the protein in solution appears relatively unstructured, and this lack of structure surrounding a few essential and well-defined signaling elements may be characteristic of a growing family of small regulatory proteins that interact with protein kinases.

Circular dichroism and NMR spectroscopy have been used to determine the structure of the low-density lipoprotein (LDL) receptor-binding peptide, comprising residues 130–152, of the human apolipoprotein E. This peptide has little persistent three-dimensional structure in solution, but when bound to micelles of dodecylphosphocholine (DPC) it adopts a predominantly α-helical structure. The three-dimensional structure of the DPC-bound peptide has been determined by using 1H-NMR spectroscopy: the structure derived from NOE-based distance constraints and restrained molecular dynamics is largely helical. The derived φ and ψ angle order parameters show that the helical structure is well defined but with some flexibility that causes the structures not to be superimposable over the full peptide length. Deuterium exchange experiments suggest that many peptide amide groups are readily accessible to the solvent, but those associated with hydrophobic residues exchange more slowly, and this helix is thus likely to be positioned on the surface of the DPC micelles. In this conformation the peptide has one hydrophobic face and two that are rich in basic amino acid side chains. The solvent-exposed face of the peptide contains residues previously shown to be involved in binding to the LDL receptor.

The structure of the trigonal crystal form of bovine β-lactoglobulin variant B at pH 7.1 has been determined by X-ray diffraction methods at a resolution of 2.22 Å and refined to values for R and Rfree of 0.239 and 0.286, respectively. By comparison with the structure of the trigonal crystal form of bovine β-lactoglobulin variant A at pH 7.1, which was determined previously [Qin BY et al., 1998, Biochemistry 37:14014–14023], the structural consequences of the sequence differences D64G and V118A of variants A and B, respectively, have been investigated. Only minor differences in the core calyx structure occur. In the vicinity of the mutation site D64G on loop CD (residues 61–67), there are small changes in main-chain conformation, whereas the substitution V118A on β-strand H is unaccompanied by changes in the surrounding structure, thereby creating a void volume and weakened hydrophobic interactions with a consequent loss of thermal stability relative to variant A. A conformational difference is found for the loop EF, implicated in the pH-dependent conformational change known as the Tanford transition, but it is not clear whether this reflects differences intrinsic to the variants in solution or differences in crystallization.