Substrate-assisted catalysis (SAC) is the process by which a functional group in a substrate contributes to catalysis by an enzyme. SAC has been demonstrated for representatives of three major enzyme classes: serine proteases, GTPases, and type II restriction endonucleases, as well as lysozyme and hexose-1-phosphate uridylyltransferase. Moreover, structure-based predictions of SAC have been made for many additional enzymes. Examples of SAC include both naturally occurring enzymes such as type II restriction endonucleases as well as engineered enzymes including serine proteases. In the latter case, a functional group from a substrate can substitute for a catalytic residue replaced by site-directed mutagenesis. From a protein engineering perspective, SAC provides a strategy for drastically changing enzyme substrate specificity or even the reaction catalyzed. From a biological viewpoint, SAC contributes significantly to the activity of some enzymes and may represent a functional intermediate in the evolution of catalysis. This review focuses on advances in engineering enzyme specificity and activity by SAC, together with the biological significance of this phenomenon.

Whereas previously we have successfully utilized the folding funnels concept to rationalize binding mechanisms (Ma B, Kumar S, Tsai CJ, Nussinov R, 1999, Protein Eng 12:713–720) and to describe binding (Tsai CJ, Kumar S, Ma B, Nussinov R, 1999, Protein Sci 8:1181–1190), here we further extend the concept of folding funnels, illustrating its utility in explaining enzyme pathways, multimolecular associations, and allostery. This extension is based on the recognition that funnels are not stationary; rather, they are dynamic, depending on the physical or binding conditions (Tsai CJ, Ma B, Nussinov R, 1999, Proc Natl Acad Sci USA 96:9970–9972). Different binding states change the surrounding environment of proteins. The changed environment is in turn expressed in shifted energy landscapes, with different shapes and distributions of populations of conformers. Hence, the function of a protein and its properties are not only decided by the static folded three-dimensional structure; they are determined by the distribution of its conformational substates, and in particular, by the redistributions of the populations under different environments. That is, protein function derives from its dynamic energy landscape, caused by changes in its surroundings.

Inspection of high resolution three-dimensional (3D) structures from the protein database reveals an increasing number of cis-Xaa-Pro and cis-Xaa-Yaa peptide bonds. However, we are still far from being able to predict whether these bonds will remain cis upon single-site substitution of Pro or Yaa and/or cleavage of a peptide bond close to it in the sequence. We have chosen oxidized Escherichia coli thioredoxin (Trx), a member of the Trx superfamily with a single α/β domain and cis P76 to determine the effect of single-site substitution and/or cleavage on this isomer. Standard two-dimensional (2D) NMR analysis were performed on cleaved Trx (1–73/74–108) and its P76A variant. Analysis of the NOE connectivities indicates remarkable similarity between the secondary and supersecondary structure of the noncovalent complexes and Trx. Analysis of the 2D version of the HCCH-TOCSY and HMQC-NOESY-HMQC and 13C-filtered HMQC-NOESY spectra of cleaved Trx with uniformly 13C-labeled I75 and P76 shows surprising conservation of both cis P76 and packing of I75 against W31. A similar NMR analysis of its P76A variant provides no evidence for cis A76 and shows only subtle local changes in both the packing of I75 and the interstrand connectivities between its most protected hydrophobic strands (β2 and β4). Indeed, a molecular simulation model for the trans P76A variant of Trx shows only subtle local changes around the substitution site. In conclusion, cleavage of R73 is insufficient to provoke cis/trans isomerization of P76, but cleavage and single-site substitution (P76A) favors the trans isomer.

The crystal structures of four active site-directed thrombin inhibitors, 1–4, in a complex with human α-thrombin have been determined and refined at up to 2.0 Å resolution using X-ray crystallography. These compounds belong to a structurally novel family of inhibitors based on a 2,3-disubstituted benzo[b]thiophene structure. Compared to traditional active-site directed inhibitors, the X-ray crystal structures of these complexes reveal a novel binding mode. Unexpectedly, the lipophilic benzo[b]thiophene nucleus of the inhibitor appears to bind in the S1 specificity pocket. At the same time, the basic amine of the C-3 side chain of the inhibitor interacts with the mostly hydrophobic proximal, S2, and distal, S3, binding sites. The second, basic amine side chain at C-2 was found to point away from the active site, occupying a location between the S1 and S′1 sites. Together, the aromatic rings of the C-2 and C-3 side chains sandwich the indole ring of Trp60D contained in the thrombin S2 insertion loop defined by the sequence “Tyr-Pro-Pro-Trp.” [The thrombin residue numbering used in this study is equivalent to that reported for chymotrypsinogen (Hartley BS, Shotton DM, 1971, The enzymes, vol. 3. New York: Academic Press. pp 323–373).] In contrast to the binding mode of more classical thrombin inhibitors (d-Phe-Pro-Arg-H, NAPAP, Argatroban), this novel class of benzo[b]thiophene derivatives does not engage in hydrogen bond formation with Gly216 of the thrombin active site. A detailed analysis of the three-dimensional structures not only provides a clearer understanding of the interaction of these agents with thrombin, but forms a foundation for rational structure-based drug design. The use of the data from this study has led to the design of derivatives that are up to 2,900-fold more potent than the screening hit 1.

The dihydrolipoamide succinyltransferase (E2o) component of the α-ketoglutarate dehydrogenase complex catalyzes the transfer of a succinyl group from the S-succinyldihydrolipoyl moiety to coenzyme A. E2o is normally a 24-mer, but is found as a trimer when E2o is expressed with a C-terminal [His]6 tag. The crystal structure of the trimeric form of the catalytic domain (CD) of the Escherichia coli E2o has been solved to 3.0 Å resolution using the Molecular Replacement method. The refined model contains an intact trimer in the asymmetric unit and has an R-factor of 0.257 (Rfree = 0.286) for 18,699 reflections between 10.0 and 3.0 Å resolution. The core of tE2oCD (residues 187–396) superimposes onto that of the cubic E2oCD with an RMS difference of 0.4 Å for all main-chain atoms. The C-terminal end of tE2oCD (residues 397–404) rotates by an average of 37° compared to cubic E2oCD, disrupting the normal twofold interface. Despite the alteration of quaternary structure, the active site of tE2oCD shows no significant differences from that of the cubic E2oCD, although several side chains in the active site are more ordered in the trimeric form of E2oCD. Analysis of the available sequence data suggests that the majority of E2 components have active sites that resemble that of E. coli E2oCD. The remaining E2 components can be divided into three groups based on active-site sequence similarity. Analysis of the surface properties of both crystal forms of E. coli E2oCD suggests key residues that may be involved in the protein–protein contacts that occur between the catalytic and lipoyl domains of E2o.

The X-ray crystallographic structures of two mutants (K206Q and H207E) of the N-lobe of human transferrin (hTF/2N) have been determined to high resolution (1.8 and 2.0 Å, respectively). Both mutant proteins bind iron with greater affinity than native hTF/2N. The structures of the K206Q and H207E mutants show interactions (both H-bonding and electrostatic) that stabilize the interaction of Lys296 in the closed conformation, thereby stabilizing the iron bound forms.

As an alternative method to study the heterotropic mechanism of Escherichia coli aspartate transcarbamoylase, a series of nucleotide analogs were used. These nucleotide analogs have the advantage over site-specific mutagenesis experiments in that interactions between the backbone of the protein and the nucleotide could be evaluated in terms of their importance for function. The ATP analogs purine 5′-triphosphate (PTP), 6-chloropurine 5′-triphosphate (Cl-PTP), 6-mercaptopurine 5′-triphosphate (SH-PTP), 6-methylpurine 5′-triphosphate (Me-PTP), and 1-methyladenosine 5′-triphosphate (Me-ATP) were partially synthesized from their corresponding nucleosides. Kinetic analysis was performed on the wild-type enzyme in the presence of these ATP analogs along with GTP, ITP, and XTP. PTP, Cl-PTP, and SH-PTP each activate the enzyme at subsaturating concentrations of l-aspartate and saturating concentrations of carbamoyl phosphate, but not to the same extent as does ATP. These experiments suggest that the interaction between N6-amino group of ATP and the backbone of the regulatory chain is important for orienting the nucleotide and inducing the displacements of the regulatory chain backbone necessary for initiation of the regulatory response. Me-PTP and Me-ATP also activate the enzyme, but in a more complex fashion, which suggests differential binding at the two sites within each regulatory dimer. The purine nucleotides GTP, ITP, and XTP each inhibit the enzyme but to a lesser extent than CTP. The influence of deoxy and dideoxynucleotides on the activity of the enzyme was also investigated. These experiments suggest that the 2′ and 3′ ribose hydroxyl groups are not of significant importance for binding and orientation of the nucleotide in the regulatory binding site. 2′-dCTP inhibits the enzyme to the same extent as CTP, indicating that the interactions of the enzyme to the O2-carbonyl of CTP are critical for CTP binding, inhibition, and the ability of the enzyme to discriminate between ATP and CTP. Examination of the electrostatic surface potential of the nucleotides and the regulatory chain suggest that the complimentary electrostatic interactions between the nucleotides and the regulatory chain are important for binding and orientation of the nucleotide necessary to induce the local conformational changes that propagate the heterotropic effect.

The conformational changes during the photocycle of the photoactive yellow protein have been the subject of many recent studies. Spectroscopic measurements have shown that the photocycle also occurs in a crystalline environment, and this has been the basis for subsequent Laue diffraction and cryocrystallographic studies. These studies have shown that conformational changes during the photocycle are limited to the chromophore and its immediate environment. However, spectroscopic studies suggest the presence of large conformational changes in the protein. Here, we address this apparent discrepancy in two ways. First, we obtain a description of large concerted motions in the ground state of the yellow protein from NMR data and theoretical calculations. Second, we describe the high-resolution structure of the yellow protein crystallized in a different space group. The structure of the yellow protein differs significantly between the two crystal forms. We show that these differences can be used to obtain a description of the flexibility of the protein that is consistent with the motions observed in solution.

Parvalbumins constitute a class of calcium-binding proteins characterized by the presence of several helix-loop-helix (EF-hand) motifs. In a previous study (Revett SP, King G, Shabanowitz J, Hunt DF, Hartman KL, Laue TM, Nelson DJ, 1997, Protein Sci 7:2397–2408), we presented the sequence of the major parvalbumin isoform from the silver hake (Merluccius bilinearis) and presented spectroscopic and structural information on the excised “EF-hand” portion of the protein. In this study, the X-ray crystal structure of the silver hake major parvalbumin has been determined to high resolution, in the frozen state, using the molecular replacement method with the carp parvalbumin structure as a starting model. The crystals are orthorhombic, space group C2221, with a = 75.7 Å, b = 80.7 Å, and c = 42.1 Å. Data were collected from a single crystal grown in 15% glycerol, which served as a cryoprotectant for flash freezing at −188 °C. The structure refined to a conventional R-value of 21% (free R 25%) for observed reflections in the range 8 to 1.65 Å [I > 2σ(I)]. The refined model includes an acetylated amino terminus, 108 residues (characteristic of a β parvalbumin lineage), 2 calcium ions, and 114 water molecules per protein molecule. The resulting structure was used in molecular dynamics (MD) simulations focused primarily on the dynamics of the ligands coordinating the Ca2+ ions in the CD and EF sites. MD simulations were performed on both the fully Ca2+ loaded protein and on a Ca2+ deficient variant, with Ca2+ only in the CD site. There was substantial agreement between the MD and X-ray results in addressing the issue of mobility of key residues in the calcium-binding sites, especially with regard to the side chain of Ser55 in the CD site and Asp92 in the EF site.

Turkey ovomucoid third domain (OMTKY3) is a canonical inhibitor of serine proteinases. Upon complex formation, the inhibitors fully exposed P1 residue becomes fully buried in the preformed cavity of the enzyme. All 20 P1 variants of OMTKY3 have been obtained by recombinant DNA technology and their equilibrium association constants have been measured with six serine proteinases. To rationalize the trends observed in this data set, high resolution crystal structures have been determined for OMTKY3 P1 variants in complex with the bacterial serine proteinase, Streptomyces griseus proteinase B (SGPB). Four high resolution complex structures are being reported in this paper; the three β-branched variants, Ile18I, Val18I, and Thr18I, determined to 2.1, 1.6, and 1.7 Å resolution, respectively, and the structure of the Ser18I variant complex, determined to 1.9 Å resolution. Models of the Cys18I, Hse18I, and Ape18I variant complexes are also discussed. The β-branched side chains are not complementary to the shape of the S1 binding pocket in SGPB, in contrast to that of the wild-type γ-branched P1 residue for OMTKY3, Leu18I. χ1 angles of approximately 40° are imposed on the side chains of Ile18I, Val18I, and Thr18I within the S1 pocket. Dihedral angles of +60°, −60°, or 180° are more commonly observed but 40° is not unfavorable for the β-branched side chains. Thr18I Oγ1 also forms a hydrogen bond with Ser195 Oγ in this orientation. The Ser18I side chain adopts two alternate conformations within the S1 pocket of SGPB, suggesting that the side chain is not stable in either conformation.

Molecular dynamics calculations provide a method by which the dynamic properties of molecules can be explored over timescales and at a level of detail that cannot be obtained experimentally from NMR or X-ray analyses. Recent work (Philippopoulos M, Mandel AM, Palmer AG III, Lim C, 1997, Proteins 28:481–493) has indicated that the accuracy of these simulations is high, as measured by the correspondence of parameters extracted from these calculations to those determined through experimental means. Here, we investigate the dynamic behavior of the Src homology 3 (SH3) domain of hematopoietic cell kinase (Hck) via 15N backbone relaxation NMR studies and a set of four independent 4 ns solvated molecular dynamics calculations. We also find that molecular dynamics simulations accurately reproduce fast motion dynamics as estimated from generalized order parameter (S2) analysis for regions of the protein that have experimentally well-defined coordinates (i.e., stable secondary structural elements). However, for regions where the coordinates are not well defined, as indicated by high local root-mean-square deviations among NMR-determined structural family members or high B-factors/low electron density in X-ray crystallography determined structures, the parameters calculated from a short to moderate length (less than 5–10 ns) molecular dynamics trajectory are dependent on the particular coordinates chosen as a starting point for the simulation.

A divalent metal ion, such as Mn2+, is required for the catalytic reaction and allosteric regulation of pig heart NAD-dependent isocitrate dehydrogenase. The enzyme is irreversibly inactivated and cleaved by Fe2+ in the presence of O2 and ascorbate at pH 7.0. Mn2+ prevents both inactivation and cleavage. Nucleotide ligands, such as NAD, NADPH, and ADP, neither prevent nor promote inactivation or cleavage of the enzyme by Fe2+. The NAD-specific isocitrate dehydrogenase is composed of three distinct subunits in the ratio 2α:1β:1γ. The results indicate that the oxidative inactivation and cleavage are specific and involve the 40 kDa α subunit of the enzyme. A pair of major peptides is generated during Fe2+ inactivation: 29.5 + 10.5 kDa, as determined by SDS-PAGE. Amino-terminal sequencing reveals that these peptides arise by cleavage of the Val262–His263 bond of the α subunit. No fragments are produced when enzyme is incubated with Fe2+ and ascorbate under denaturing conditions in the presence of 6 M urea, indicating that the native structure is required for the specific cleavage. These results suggest that His263 of the α subunit may be a ligand of the divalent metal ion needed for the reaction catalyzed by isocitrate dehydrogenase. Isocitrate enhances the inactivation of enzyme caused by Fe2+ in the presence of oxygen, but prevents the cleavage, suggesting that inactivation occurs by a different mechanism when metal ion is bound to the enzyme in the presence of isocitrate: oxidation of cysteine may be responsible for the rapid inactivation in this case. Affinity cleavage caused by Fe2+ implicates α as the catalytic subunit of the multisubunit porcine NAD-dependent isocitrate dehydrogenase.

An analysis of the folding of the 94 residue tenth fibronectin type III (fnIII) domain of human fibronectin (FNfn10) is presented. Use of guanidine isothiocyanate as a denaturant allows us to obtain equilibrium and kinetic data across a broad range of denaturant concentrations that are unavailable in guanidine hydrochloride. Equilibrium unfolding experiments show that FNfn10 is significantly more stable than has been reported previously. Comparison of equilibrium and kinetic parameters reveals the presence of an intermediate that accumulates at low denaturant concentrations. This is the first demonstration of three-state folding kinetics for a fnIII domain. We have previously shown that a homologous domain from human tenascin (TNfn3) folds by a two-state mechanism, but this does not necessarily indicate that the two proteins fold by different folding pathways.

The nitrogenase enzyme of Klebsiella pneumoniae consists of two separable proteins, each with multiple subunits and one or more oxygen sensitive metallocenters. The wild-type nitrogenase proteins are stable to electrophoresis in high concentrations of urea under anaerobic conditions. Addition of Mg+2 and ADP greatly increases the stability of the smaller Fe protein (from <4 to >6 M for full unfolding), an effect directly analogous to stabilization in p21ras induced by Mg+2 and GDP. Stabilization by Mg+2 is slight for the holo MoFe protein (from ∼1.5 to ∼2.4 M) but more dramatic for the apo protein form of the MoFe protein accumulated by certain Fe protein (nifH gene) mutants. The potent product inhibitor of nitrogenase function, MgADP, increases stability of the MoFe protein more than Mg+2 alone, to ∼3.6 M, showing that nucleotides interact with the MoFe protein. Mutations of the nifM gene result in slower accumulation of less stable Fe protein, indicating that NifM is involved in correct folding of the Fe protein. Mutationally altered proteins are often difficult to purify for study because of their inherent instability, low expression level, or oxygen lability. Crude extracts of 11 different mutants of Fe protein (nifH gene) were examined by transverse urea gradient gels to rapidly screen for stabilizing interactions in the presence or absence of substrate or inhibitor analogs. Amino acid alterations D44N and R188C, at the interface of the dimer, in the vicinity of the nucleotide binding site(s), have significantly lower stability than the wild-type enzyme in the absence of Mg+2 but comparable stability in its presence, showing the importance of Mg+2 in the subunit interactions. Mutations N163S and E266K, in which residues normally involved in hydrogen bonding far from the active site were altered, are more labile than the wild-type even with Mg+2 added. Seven other mutants, though nonfunctional, did not appear altered in stability compared to the wild-type.

The rate of macromolecular surface formation in yeast iso-2 cytochrome c and its site-specific mutant, N52I iso-2, has been studied using a monoclonal antibody that recognizes a tertiary epitope including K58 and H39. The results indicate that epitope refolding occurs after fast folding but prior to slow folding, in contrast to horse cytochrome c where surface formation occurs early. The antibody-detected (ad) kinetic phase accompanying epitope formation has kad = 0.2 s−1 and is ∼40-fold slower than the fastest detectable event in the folding of yeast iso-2 cytochrome c (k2f ∼ 8 s−1), but occurs prior to the absorbance- and fluorescence-detected slow folding steps (k1a ∼ 0.06 s−1; k1b ∼ 0.09 s−1). N52I iso-2 cytochrome c exhibits similar kinetic behavior with respect to epitope formation. A detailed dissection of the mechanistic differences between the folding pathways of horse and yeast cytochromes c identifies possible reasons for the slow surface formation in the latter. Our results suggest that non-native ligation involving H33 or H39 during refolding may slow down the formation of the tertiary epitope in iso-2 cytochrome c. This study illustrates that surface formation can be coupled to early events in protein folding. Thus, the rate of macromolecular surface formation is fine tuned by the residues that make up the surface and the interactions they entertain during refolding.

MAP30 is a 30 kDa single-stranded, type-I ribosome inactivating protein (RIP) possessing anti-tumor and anti-HIV activities. It binds both ribosomal RNA and the HIV-1 long-terminal repeat DNA. To understand the structural basis for MAP30 activities, we undertook the study of MAP30 by solution NMR spectroscopy. We report nearly complete 1H, 13C, and 15N chemical shift assignments of its 263 amino acids. Based upon an analysis of secondary 13C chemical shifts, 3JHNHA coupling constants, hydrogen exchange data, and nuclear Overhauser effect patterns, we find that the secondary structure and β-sheet topology of MAP30 are very similar to those of the ricin A chain, a subunit of the well-known type-II RIP, even though two proteins display distinct activities. We therefore suggest that MAP30 and ricin A chain share a similar three-dimensional fold, and that the reported functional differences between two proteins arise primarily from differences in local three-dimensional structure and other structural properties such as surface electrostatic potentials.

The denaturant-induced (un)folding of apoflavodoxin from Azotobacter vinelandii has been followed at the residue level by NMR spectroscopy. NH groups of 21 residues of the protein could be followed in a series of 1H–15N heteronuclear single-quantum coherence spectra recorded at increasing concentrations of guanidinium hydrochloride despite the formation of protein aggregate. These NH groups are distributed throughout the whole apoflavodoxin structure. The midpoints of unfolding determined by NMR coincide with the one obtained by fluorescence emission spectroscopy. Both techniques give rise to unfolding curves with transition zones at significantly lower denaturant concentrations than the one obtained by circular dichroism spectroscopy. The NMR (un)folding data support a mechanism for apoflavodoxin folding in which a relatively stable intermediate is involved. Native apoflavodoxin is shown to cooperatively unfold to a molten globule-like state with extremely broadened NMR resonances. This initial unfolding step is slow on the NMR chemical shift timescale. The subsequent unfolding of the molten globule is faster on the NMR chemical shift timescale and the limited appearance of 1H–15N HSQC cross peaks of unfolded apoflavodoxin in the denaturant range studied indicates that it is noncooperative.

A method is presented that allows the calculation of the lifetimes of tryptophan residues on the basis of spectral and structural data. It is applied to two different proteins. The calcium binding protein from the sarcoplasm of the muscles of the sand worm Nereis diversicolor (NSCP) changes its conformation upon binding of Ca2+ or Mg2+. NSCP contains three tryptophan residues at position 4, 57, and 170, respectively. The fluorescence lifetimes of W57 are investigated in a mutant in which W4 and W170 have been replaced. The time resolved fluorescence properties of W57 are linked to its different microconformations, which were determined by a molecular dynamics simulation map. Together with the determination of the radiative rate constant from the wavelength of maximum intensity of the decay associated spectra, it was possible to determine an exponential relation between the nonradiative rate constant and the distance between the indole CE3 atom and the carbonyl carbon of the peptide bond reflecting a mechanism of electron transfer as the main determinant of the value for the nonradiative rate constant. This result allows the calculation of the fluorescence lifetimes from the protein structure and the spectra. This method was further tested for the tryptophan of Ha-ras p21 (W32) and for W43 of the Tet repressor, which resulted in acceptable values for the predicted lifetimes.

High mobility group (HMG) proteins 14 and 17 are nonhistone nuclear proteins that have been implicated in control of transcription and chromatin structure. To examine the posttranslational modifications of HMG-14 and -17 in vivo, HMG proteins were prepared from nuclear vs. cytosolic fractions of human K562 cells treated with 12-O-tetradecanoylphorbol 13-acetate (TPA) or okadaic acid (OA) and examined by electrospray mass spectrometry. Analysis of full-length masses demonstrated mono-, di-, and triphosphorylation of HMG-14 and mono- and diphosphorylation of HMG-17 from OA treated cells, whereas HMG-14 and -17 from TPA treated cells were monophosphorylated. Peptide mass and sequence analysis showed major and minor phosphorylation sites, respectively, at Ser24 and Ser28 in HMG-17, and Ser20 and Ser24 in HMG-14. These sites were found in the consensus sequence RRSARLSAK, within the nucleosomal binding domain of each protein. A third phosphorylation site in HMG-14 was located at either Ser6 or Ser7. Interestingly, the proportion of HMG-14 and -17 found in cytosolic pools increased significantly after 1 h of treatment compared to control cells and showed preferential phosphorylation compared with proteins from nuclear fractions. These results suggest that phosphorylation of HMG-14 and -7 interferes with nuclear localization mechanisms in a manner favoring release from nuclei.

Comparisons of protein sequence via cyclic training of Hidden Markov Models (HMMs) in conjunction with alignments of three-dimensional structure, using the Combinatorial Extension (CE) algorithm, reveal two putative EF-hand metal binding domains in acetylcholinesterase. Based on sequence similarity, putative EF-hands are also predicted for the neuroligin family of cell surface proteins. These predictions are supported by experimental evidence. In the acetylcholinesterase crystal structure from Torpedo californica, the first putative EF-hand region binds the Zn2+ found in the heavy metal replacement structure. Further, the interaction of neuroligin 1 with its cognate receptor neurexin depends on Ca2+. Thus, members of the α,β hydrolase fold family of proteins contain potential Ca2+ binding sites, which in some family members may be critical for heterologous cell associations.