Catalytic mechanism of Escherichia coli glycinamide ribonucleotide transformylase probed by site-directed mutagenesis and pH-dependent studies.

Site-directed mutagenesis followed by studies of the pH dependence of the kinetic parameters of the mutants has been used to probe the role of the active site residues and loops in catalysis by glycinamide ribonucleotide transformylase (EC 2.1.2.2). The analysis of the mutants of the strictly conserved active site residues, His108 and Asp144, revealed that His108 acts in a salt bridge with Asp144 as a general acid catalyst with a pK(a) value of 9.7. Asp144 also plays a key role in the preparation of the active site geometry for catalysis. The rate-limiting step in the pH range of 6-10 appears to be the catalytic steps involving tetrahedral intermediates, supported by the observation of a pL (L being H or D)-independent solvent deuterium isotope effect of 2. The ionization of the amino group of glycinamide ribonucleotide both as a free and as a bound form dominates the kinetic behavior at low pH. The analysis of a mutation, H121Q, within the loop spanning amino acids 111-131 suggests the closure of the loop is involved in the binding of the substrate. The kinetic behavior parallels pH effects revealed by a series of X-ray crystallographic structures of the apoenzyme and inhibitor-bound enzyme [Su, Y., Yamashita, M. M., Greasley, S. E. , Mullen, C. A., Shim, J. H., Jennings, P. A., Benkovic, S. J., and Wilson, I. A. (1998) J. Mol. Biol. 281, 485-499], permitting a more exact formulation of the probable catalytic mechanism.     

The pH dependence of stability of the activation helix and the catalytic site of GART

We have predicted the free energy of unfolding for the pH-dependent helix-coil transition of the activation helix of GART using continuum electrostatic calculations and structural modeling. We have assigned the contributions of each element of secondary structure and of each ionizable residue, within and in the vicinity of the activation helix, to the stability of several fragments of GART that participate in the formation of the catalytic site. We demonstrate that the interaction of His121-His132 contributes 2.2 kcal/mol to the ionization free energy between pH 0 and approximately 6. We also show that the ionization state of a network of five histidines, His108, His119, His121, His132 and His137, and two aspartic acids Asp141 and Asp144, contributes approximately 12 kcal/mol to the stability of the catalytic site of GART, out of a total stability of 16 kcal/mol of the whole enzyme. These interactions are important for the formation of the catalytic site of GART.     

Native-state conformational dynamics of GART: a regulatory pH-dependent coil-helix transition examined by electrostatic calculations

Glycinamide ribonucleotide transformylase (GART) undergoes a pH-dependent coil-helix transition with pK(a) approximately 7. An alpha-helix is formed at high pH spanning 8 residues of a 21-residue-long loop, comprising the segment Thr120-His121-Arg122-Gln123-Ala124-Leu125-Glu126-Asn127. To understand the electrostatic nature of this loop-helix, called the activation loop-helix, which leads to the formation and stability of the alpha-helix, pK(a) values of all ionizable residues of GART have been calculated, using Poisson-Boltzmann electrostatic calculations and crystallographic data. Crystallographic structures of high and low pH E70A GART have been used in our analysis. Low pK(a) values of 5.3, 5.3, 3.9, 1.7, and 4.7 have been calculated for five functionally important histidines, His108, His119, His121, His132, and His137, respectively, using the high pH E70A GART structure. Ten theoretical single and double mutants of the high pH E70A structure have been constructed to identify pairwise interactions of ionizable residues, which have aided in elucidating the multiplicity of electrostatic interactions of the activation loop-helix, and the impact of the activation helix on the catalytic site. Based on our pK(a) calculations and structural data, we propose that: (1) His121 forms a molecular switch for the coil-helix transition of the activation helix, depending on its protonation state; (2) a strong electrostatic interaction between His132 and His121 is observed, which can be of stabilizing or destabilizing nature for the activation helix, depending on the relative orientation and protonation states of the rings of His121 and His132; (3) electrostatic interactions involving His119 and Arg122 play a role in the stability of the activation helix; and (4) the activation helix contains the helix-promoting sequence Arg122-Gln123-Ala124-Leu125-Glu126, but its alignment relative to the N and C termini of the helix is not optimal, and is possibly of a destabilizing nature. Finally, we provide electrostatic evidence that the formation and closure of the activation helix create a hydrophobic environment for catalytic-site residue His108, to facilitate catalysis.     

Human glycinamide ribonucleotide transformylase: active site mutants as mechanistic probes

Human glycinamide ribonucleotide transformylase (GART) (EC 2.1.2.2) is a validated target for cancer chemotherapy, but mechanistic studies of this therapeutically important enzyme are limited. Site-directed mutagenesis, initial velocity studies, pH-rate studies, and substrate binding studies have been employed to probe the role of the strictly conserved active site residues, N106, H108, and D144, and the semiconserved K170 in substrate binding and catalysis. Only two conservative substitutions, N106Q and K170R, resulted in catalytically active enzymes, and these active mutant enzymes gave pH-rate profiles and a steady-state kinetic mechanism essentially identical to those of the native enzyme. All inactive mutants were able to bind both substrates, ruling out disrupted formation of the ternary complex as the source of inactivity. Differences between human and Escherichia coli GART, previously used as a model for the human enzyme, were evident.     

Active-site mapping and site-specific mutagenesis of glycinamide ribonucleotide transformylase from Escherichia coli

The affinity reagent N10-(bromoacetyl)-5,8-dideazafolate has previously been shown to inactivate glycinamide ribonucleotide transformylase (EC 2.1.2.2) from Escherichia coli in an active-site-directed manner with a 1:1 stoichiometry [Inglese et al. (1990) Biochemistry 29, 1436-1443]. After a series of mild proteolytic digestions, the dideazafolate label was localized to an active-site peptide attached by an ester linkage to the highly conserved residue Asp 144. Subsequent site-specific mutagenesis of Asp 144 to Asn 144 resulted in a catalytically inactive enzyme that retained the ability to bind substrates and inhibitors. The Asn 144 mutant could be further labeled with the affinity reagent in an active-site-directed stoichiometric fashion; however, the site of modification in this case was His 119. These results imply that Asp 144 may function as a general base within the catalytic center of the transformylase and is in close proximity to His 119 in the folded protein.     

Crystal structure of glycinamide ribonucleotide transformylase from Escherichia coli at 3.0 A resolution. A target enzyme for chemotherapy.

The atomic structure of glycinamide ribonucleotide transformylase, an essential enzyme in purine biosynthesis, has been determined at 3.0 A resolution. The last three C-terminal residues and a sequence stretch of 18 residues (residues 113 to 130) are not visible in the electron density map. The enzyme forms a dimer in the crystal structure. Each monomer is divided into two domains, which are connected by a central mainly parallel seven-stranded beta-sheet. The N-terminal domain contains a Rossmann type mononucleotide fold with a phosphate ion bound to the C-terminal end of the first beta-strand. A long narrow cleft stretches from the phosphate to a conserved aspartic acid, Asp144, which has been suggested as an active-site residue. The cleft is lined by a cluster of residues, which are conserved between bacterial, yeast, avian and human enzymes, and likely represents the binding pocket and active site of the enzyme. GAR Tfase binds a reduced folate cofactor and glycinamide ribonucleotide for the catalysis of one of the initial steps in purine biosynthesis. Folate analogs and multi-substrate inhibitors of the enzyme have antineoplastic effects and the structure determination of the unliganded enzyme and enzyme-inhibitor complexes will aid the development of anti-cancer drugs.     

The apo and ternary complex structures of a chemotherapeutic target: human glycinamide ribonucleotide transformylase

Glycinamide ribonucleotide transformylase (GART; 10-formyltetrahydrofolate:5'-phosphoribosylglycinamide formyltransferase, EC 2.1.2.2), an essential enzyme in de novo purine biosynthesis, has been a chemotherapeutic target for several decades. The three-dimensional structure of the GART domain from the human trifunctional enzyme has been solved by X-ray crystallography. Models of the apoenzyme, and a ternary complex with the 10-formyl-5,8-dideazafolate cosubstrate and a glycinamide ribonucleotide analogue, hydroxyacetamide ribonucleotide [alpha,beta-N-(hydroxyacetyl)-d-ribofuranosylamine], are reported to 2.2 and 2.07 A, respectively. The model of the apoenzyme represents the first structure of GART, from any source, with a completely unoccupied substrate and cosubstrate site, while the ternary complex is the first structure of the human GART domain that is bound at both the substrate and cosubstrate sites. A comparison of the two models therefore reveals subtle structural differences that reflect substrate and cosubstrate binding effects and implies roles for the invariant residues Gly 133, Gly 146, and His 137. Preactivation of the DDF formyl group appears to be key for catalysis, and structural flexibility of the active end of the substrate may facilitate nucleophilic attack. A change in pH, rather than folate binding, correlates with movement of the folate binding loop, whereas the phosphate binding loop position does not vary with pH. The electrostatic surface potentials of the human GART domain and Escherichia coli enzyme explain differences in the binding affinity of polyglutamylated folates, and these differences have implications to future chemotherapeutic agent design.     

Improvement in the efficiency of formyl transfer of a GAR transformylase hybrid enzyme

A hybrid glycinamide ribonucleotide transformylase was assembled from two protein domains that were treated as discrete modules. One module contained the ribonucleotide binding domain from the purN glycinamide ribonucleotide transformylase; the second module contained the catalytic machinery and the formyl tetrahydrofolate binding domain from the enzyme encoded by purU, formyl tetrahydrofolate hydrolase. The resultant enzyme showed 0.1% catalytic activity of the wild-type glycinamide ribonucleotide transformylase enzyme but had a formyl transfer efficiency of 10%. A combinatorial mutagenesis approach was used to improve the solubility and formyl transfer properties of the hybrid enzyme. The mutagenized hybrid glycinamide ribonucleotide transformylase was initially expressed as a fusion to the alpha-peptide of beta-galactosidase. Clones were selected for improvement in solubility by determining which clones were capable of alpha-complementation using a blue/white screen. One clone was further characterized and found to have an improved efficiency of transfer of the ribonucleotide increasing this property to >95%.     

Roles of active site residues in Pseudomonas aeruginosa phosphomannomutase/phosphoglucomutase

In Pseudomonas aeruginosa, the dual-specificity enzyme phosphomannomutase/phosphoglucomutase catalyzes the transfer of a phosphoryl group from serine 108 to the hydroxyl group at the 1-position of the substrate, either mannose 6-P or glucose 6-P. The enzyme must then catalyze transfer of the phosphoryl group on the 6-position of the substrate back to the enzyme. Each phosphoryl transfer is expected to require general acid-base catalysis, provided by amino acid residues at the enzyme active site. An extensive survey of the active site residues by site-directed mutagenesis failed to identify a single key residue that mediates the proton transfers. Mutagenesis of active site residues Arg20, Lys118, Arg247, His308, and His329 to residues that do not contain ionizable groups produced proteins for which V(max) was reduced to 4-12% of that of the wild type. The fact that no single residue decreased catalytic activity more significantly, and that several residues had similar effects on V(max), suggested that the ensemble of active site amino acids act by creating positive electrostatic potential, which serves to depress the pK of the substrate hydroxyl group so that it binds in ionized form at the active site. In this way, the necessity of positioning the reactive hydroxyl group near a specific amino acid residue is avoided, which may explain how the enzyme is able to promote catalysis of both phosphoryl transfers, even though the 1- and 6-positions do not occupy precisely the same position when the substrate binds in the two different orientations in the active site. When Ser108 is mutated, the enzyme retains a surprising amount of activity, which has led to the suggestion that an alternative residue becomes phosphorylated in the absence of Ser108. (31)P NMR spectra of the S108A protein confirm that it is phosphorylated. Although the S108A/H329N protein had no detectable catalytic activity, the (31)P NMR spectra were not consistent with a phosphohistidine residue.     

Measurement of glycinamide ribonucleotide synthetase activity in intact human fibroblasts and Chinese hamster ovary cells

An intact cell assay system based on Tween 80 permeabilization was used to investigate glycinamide ribonucleotide (GAR) synthetase activity in human fibroblasts and Chinese hamster ovary cells. Optimal conditions for the assay of the enzyme were determined with regards to ATP, MgCl2, NH4Cl and ribose-5'-phosphate concentrations as well as pH. Using the optimal assay conditions, the Vmax values as determined by Lineweaver-Burke double reciprocal plots were found to be 5.19 nmol GAR formed/5 X 10(5) cells/30 min for the fibroblasts and 13.4 nmol GAR formed/5 X 10(5) cells/30 min for the Chinese hamster ovary cells.     

A rationalization of the acidic pH dependence for stromelysin-1 (Matrix metalloproteinase-3) catalysis and inhibition

The pH dependence of matrix metalloproteinase (MMP) catalysis is described by a broad bell-shaped curve, indicating the involvement of two unspecified ionizable groups in proteolysis. Stromelysin-1 has a third pK(a) near 6, resulting in a uniquely sharp acidic catalytic optimum, which has recently been attributed to His(224). This suggests the presence of a critical, but unidentified, S1' substructure. Integrating biochemical characterizations of inhibitor-enzyme interactions with active site topography from corresponding crystal structures, we isolated contributions to the pH dependence of catalysis and inhibition of active site residues Glu(202) and His(224). The acidic pK(a) 5.6 is attributed to the Glu(202).zinc.H(2)O complex, consistent with a role for the invariant active site Glu as a general base in MMP catalysis. The His(224)-dependent substructure is identified as a tripeptide (Pro(221)-Leu(222)-Tyr(223)) that forms the substrate cleft lower wall. Substrate binding induces a beta-conformation in this sequence, which extends and anchors the larger beta-sheet of the enzyme. substrate complex and appears to be essential for productive substrate binding. Because the PXY tripeptide is strictly conserved among MMPs, this "beta-anchor" may represent a common motif required for macromolecular substrate hydrolysis. The striking acidic profile of stromelysin-1 defined by the combined ionization of Glu(202) and His(224) allows the design of highly selective inhibitors.     

Structural insights into the avian AICAR transformylase mechanism

ATIC encompasses both AICAR transformylase and IMP cyclohydrolase activities that are responsible for the catalysis of the penultimate and final steps of the purine de novo synthesis pathway. The formyl transfer reaction catalyzed by the AICAR Tfase domain is substantially more demanding than that catalyzed by the other folate-dependent enzyme of the purine biosynthesis pathway, GAR transformylase. Identification of the AICAR Tfase active site and key catalytic residues is essential to elucidate how the non-nucleophilic AICAR amino group is activated for formyl transfer. Hence, the crystal structure of dimeric avian ATIC was determined as a complex with the AICAR Tfase substrate AICAR, as well as with an IMP cyclohydrolase inhibitor, XMP, to 1.93 A resolution. AICAR is bound at the dimer interface of the transformylase domains and forms an extensive hydrogen bonding network with a multitude of active site residues. The crystal structure suggests that the conformation of the 4-carboxamide of AICAR is poised to increase the nucleophilicity of the C5 amine, while proton abstraction occurs via His(268) concomitant with formyl transfer. Lys(267) is likely to be involved in the stabilization of the anionic formyl transfer transition state and in subsequent protonation of the THF leaving group.     

Crystal structure of haloalkane dehalogenase LinB from Sphingomonas paucimobilis UT26 at 0.95 A resolution: dynamics of catalytic residues

We present the structure of LinB, a 33-kDa haloalkane dehalogenase from Sphingomonas paucimobilis UT26, at 0.95 A resolution. The data have allowed us to directly observe the anisotropic motions of the catalytic residues. In particular, the side-chain of the catalytic nucleophile, Asp108, displays a high degree of disorder. It has been modeled in two conformations, one similar to that observed previously (conformation A) and one strained (conformation B) that approached the catalytic base (His272). The strain in conformation B was mainly in the C(alpha)-C(beta)-C(gamma) angle (126 degrees ) that deviated by 13.4 degrees from the "ideal" bond angle of 112.6 degrees. On the basis of these observations, we propose a role for the charge state of the catalytic histidine in determining the geometry of the catalytic residues. We hypothesized that double-protonation of the catalytic base (His272) reduces the distance between the side-chain of this residue and that of the Asp108. The results of molecular dynamics simulations were consistent with the structural data showing that protonation of the His272 side-chain nitrogen atoms does indeed reduce the distance between the side-chains of the residues in question, although the simulations failed to demonstrate the same degree of strain in the Asp108 C(alpha)-C(beta)-C(gamma) angle. Instead, the changes in the molecular dynamics structures were distributed over several bond and dihedral angles. Quantum mechanics calculations on LinB with 1-chloro-2,2-dimethylpropane as a substrate were performed to determine which active site conformations and protonation states were most likely to result in catalysis. It was shown that His272 singly protonated at N(delta)(1) and Asp108 in conformation A gave the most exothermic reaction (DeltaH = -22 kcal/mol). With His272 doubly protonated at N(delta)(1) and N(epsilon)(2), the reactions were only slightly exothermic or were endothermic. In all calculations starting with Asp108 in conformation B, the Asp108 C(alpha)-C(beta)-C(gamma) angle changed during the reaction and the Asp108 moved to conformation A. The results presented here indicate that the positions of the catalytic residues and charge state of the catalytic base are important for determining reaction energetics in LinB.     

10-(2-benzoxazolcarbonyl)-5,10-dideaza-acyclic-5,6,7,8-tetrahydrofolic acid: a potential inhibitor of GAR transformylase and AICAR transformylase

The design and synthesis of 10-(2-benzoxazolcarbonyl)-DDACTHF (1) as an inhibitor of glycinamide ribonucleotide transformylase (GAR Tfase) and aminoimidazole carboxamide transformylase (AICAR Tfase) are reported. Ketone 1 and the corresponding alcohol 13 were evaluated for inhibition of GAR Tfase and AICAR Tfase and the former was found to be a potent inhibitor of recombinant human (rh) GAR Tfase (Ki=600 nM).     

Glycinamide ribonucleotide synthetase from Escherichia coli: cloning, overproduction, sequencing, isolation, and characterization

The purD gene of Escherichia coli encoding the enzyme glycinamide ribonucleotide (GAR) synthetase, which catalyzes the conversion of phosphoribosylamine (PRA), glycine, and MgATP to glycinamide ribonucleotide, MgADP, and Pi, has been cloned and sequenced. The protein, as deduced by the structural gene sequence, contains 430 amino acids and has a calculated Mr of 45,945. Construction of an overproducing strain behind a lambda pL promoter allowed a 4-fold purification of the protein to homogeneity. N-Terminal sequence analysis and comparison of the sequence with those of other GAR synthetases confirm the amino acid sequence deduced from the gene sequence. Initial velocity studies and product and dead-end inhibition studies are most consistent with a sequential ordered mechanism of substrate binding and product release in which PRA binds first followed by MgATP and then glycine; Pi leaves first, followed by loss of MgADP and finally GAR. Incubation of [18O]glycine, ATP, and PRA results in quantitative transfer of the 18O to Pi. GAR synthetase is very specific for its substrate glycine.     

pK values of the ionizable groups of proteins

We have used potentiometric titrations to measure the pK values of the ionizable groups of proteins in alanine pentapeptides with appropriately blocked termini. These pentapeptides provide an improved model for the pK values of the ionizable groups in proteins. Our pK values determined in 0.1 M KCl at 25 degrees C are: 3.67+/-0.03 (alpha-carboxyl), 3.67+/-0.04 (Asp), 4.25+/-0.05 (Glu), 6.54+/-0.04 (His), 8.00+/-0.03 (alpha-amino), 8.55+/-0.03 (Cys), 9.84+/-0.11 (Tyr), and 10.40+/-0.08 (Lys). The pK values of some groups differ from the Nozaki and Tanford (N & T) pK values often used in the literature: Asp (3.67 this work vs. 4.0 N & T); His (6.54 this work vs. 6.3 N & T); alpha-amino (8.00 this work vs. 7.5 N & T); Cys (8.55 this work vs. 9.5 N & T); and Tyr (9.84 this work vs. 9.6 N & T). Our pK values will be useful to those who study pK perturbations in folded and unfolded proteins, and to those who use theory to gain a better understanding of the factors that determine the pK values of the ionizable groups of proteins.     

Kinetic and mutational analyses of the major cytosolic exopolyphosphatase from Saccharomyces cerevisiae

Yeast exopolyphosphatase (scPPX) processively splits off the terminal phosphate group from linear polyphosphates longer than pyrophosphate. scPPX belongs to the DHH phosphoesterase superfamily and is evolutionarily close to the well characterized family II pyrophosphatase (PPase). Here, we used steady-state kinetic and binding measurements to elucidate the metal cofactor requirement for scPPX catalysis over the pH range 4.2-9.5. A single tight binding site for Mg(2+) (K(d) of 24 microm) was detected by equilibrium dialysis. Steady-state kinetic analysis of tripolyphosphate hydrolysis revealed a second site that binds Mg(2+) in the millimolar range and modulates substrate binding. This step requires two protonated and two deprotonated enzyme groups with pK(a) values of 5.0-5.3 and 7.6-8.2, respectively. The catalytic step requiring two deprotonated groups (pK(a) of 4.6 and 5.6) is modulated by ionization of a third group (pK(a) of 8.7). Conservative mutations of Asp(127), His(148), His(149) (conserved in scPPX and PPase), and Asn(35) (His in PPase) reduced activity by a factor of 600-5000. N35H and D127E substitutions reduced the Mg(2+) affinity of the tight binding site by 25-60-fold. Contrary to expectations, the N35H variant was unable to hydrolyze pyrophosphate, but markedly altered metal cofactor specificity, displaying higher catalytic activity with Co(2+) bound to the weak binding site versus the Mg(2+)- or Mn(2+)-bound enzyme. These results provide an initial step toward understanding the dynamics of scPPX catalysis and reveal significant functional differences between structurally similar scPPX and family II PPase.     

Conformationally restricted analogues designed for selective inhibition of GAR Tfase versus thymidylate synthase or dihydrofolate reductase

The synthesis and evaluation of a series of conformationally restricted analogues of 10-formyl-tetrahydrofolate as potential inhibitors of glycinamide ribonucleotide transformylase (GAR Tfase) or aminoimidazole carboxamide transformylase (AICAR Tfase) are reported.     

A quantum chemical study on the mechanism of glycinamide ribonucleotide transformylase inhibitor: 10-Formyl-5,8,10-trideazafolic acid

A density functional theory (DFT) study is presented on the reaction mechanism of glycinamide ribonucleotide (GAR) with 10-formyl-5,8,10-trideazafolic acid (10f-TDAF), which is an inhibitor designed for GAR transformylase (GAR Tfase). There are three different paths for this system and the results indicate that inhibitor 10f-TDAF can form a very stable intermediate with the substrate GAR or generate an imine bond with GAR by elimination of water. The results have verified the presumption from available experiments and implied that 10f-TDAF would be an important target for anti-neoplastic intervention.     

Interaction between the carboxyl groups of Asp127 and Asp129 in the active site of Escherichia coli phosphofructokinase.

The pH dependence of the enzymic properties of the phosphofructokinase from Escherichia coli was compared to those of two mutants in which one carboxyl group of the active site has been removed from either Asp127 or Asp129. All measurements of activity were made in the presence of allosteric activator ADP or GDP to eliminate any cooperative process. Asp129 is a crucial residue for the activity of phosphofructokinase since its conversion to Ser decreases the catalytic activity by 2-3 orders of magnitude in both the forward and reverse reactions, but the ionization of Asp129 is not directly related the pH dependence of phosphofructokinase activity. This pH dependence is however modified by the Asp129----Ser mutation, which decreases the pK of another residue, Asp127, by as much as pH of 1.5. The side chain of Asp127 has the catalytic role proposed earlier: its deprotonated form acts as a base in the forward reaction, and its protonated form acts as an acid in the reverse reaction. The protonated form of Asp127 is also required for the binding of fructose 1,6-bisphosphate. The electrostatic interaction between the carboxyl groups of Asp127 and Asp129 seems different in free phosphofructokinase to that in enzyme/substrate complexes, suggesting that a conformational change occurs upon substrate binding. The pH dependence of phosphofructokinase activity involves one other ionizable group with a pK of approximately 6 which does not belong to the side chains of Asp127 or Asp129.