Catalytic mechanism of thioltransferase

To evaluate potential catalytic mechanism for thioltransferase thiol-disulfide exchange reactions, seven pig liver mutants were constructed by site-directed mutagenesis. All the expressed enzymes, including wild-type and mutants with the exception of the inactive mutant, ETT-C22S, were variably inhibited by iodoacetamide, and similar results were obtained when these enzymes were preincubated with GSH. However, when preincubated with S-sulfocysteine or hydroxyethyl disulfide, the activity of the enzymes was totally or partially protected against inhibition by iodoacetamide, with the exception of the mutants, ETT-C25S and ETT-C25A. When simultaneously pretreated with GSH and S-sulfocysteine, all enzymes were highly protected. Isoelectric focusing analysis of the above preincubation mixtures showed that different enzyme-substrate intermediates occurred. Using radioactively labeled substrates, [U-14C]cystine and [glycine-2-3H] GSH, enzyme-substrate intermediates were detected. These data indicate that reduced thioltransferase reacts first with disulfide substrates, then with a thiol substrate, e.g. GSH. The formation of either enzyme-substrate mixed disulfide or protein intramolecular disulfide protected the enzyme from inactivation by iodoacetamide. Based on the experimental results, alternative methods of the catalytic mechanism for thioltransferases are proposed.     

GSH-dependent peroxidase activity of the rice (Oryza sativa) glutaredoxin,a thioltransferase

Glutaredoxin (Grx) is a 12-kDa thioltransferase that reduces disulfide bonds of other proteins and maintains the redox potential of cells. In addition to its oxidoreductase activity, we report here that a rice Grx (OsGrx) can also function as a GSH-dependent peroxidase. Because of this antioxidant activity, OsGrx protects glutamine synthetase from oxidative damage. Individually replacing the conserved Cys residues in OsGrx with Ser shows that Cys(23), but not Cys(26), is essential for the thioltransferase and GSH-dependent peroxidase activities. Kinetic characterization of OsGrx reveals that the maximal catalytic efficiency (V(max)/K(m)) is obtained with cumene hydroperoxide rather than H(2)O(2) or t-butyl hydroperoxide.     

Enzyme I of the phosphoenolpyruvate:sugar phosphotransferase system. In vitro intragenic complementation: the roles of Arg126 in phosphoryl transfer and the C-terminal domain in dimerization

Enzyme I mutants of the Salmonella typhimurium phosphoenolpyruvate:sugar phosphotransferase system (PTS), which show in vitro intragenic complementation, have been identified as Arg126Cys (strain SB1690 ptsI34), Gly356Ser (strain SB1681 ptsI16), and Arg375Cys (strain SB1476 ptsI17). The mutation Arg126Cys is in the N-terminal HPr-binding domain, and complements Gly356Ser and Arg375Cys enzyme I mutations located in the C-terminal phosphoenolpyruvate(PEP)-binding domain. Complementation results in the formation of unstable heterodimers. None of the mutations alters the K(m) for HPr, which is phosphorylated by enzyme I. Arg126 is a conserved residue; the Arg126Cys mutation gives a V(max) of 0.04% wild-type, establishing a role in phosphoryl transfer. The Gly356Ser and Arg375Cys mutations reduce enzyme I V(max) to 4 and 2%, respectively, and for both, the PEP K(m) is increased from 0.1 to 3 mM. It is concluded that this activity was from the monomer, rather than the dimer normally found in assays of wild-type. In the presence of Arg126Cys enzyme, V(max) for Gly356Ser and Arg375Cys enzymes I increased 6- and 2-fold, respectively; the K(m) for PEP decreased to <10 microM, but the K(m) became dependent upon the stability of the heterodimer in the assay. Gly356 is conserved in enzyme I and pyruvate phosphate dikinase, which is a homologue of enzyme I, and this residue is part of a conserved sequence in the subunit interaction site. Gly356Ser mutation impairs enzyme I dimerization. The mutation Arg375Cys also impairs dimerization, but the equivalent residue in pyruvate phosphate dikinase is not associated with the subunit interaction site. A 37 000 Da, C-terminal domain of enzyme I has been expressed and purified; it dimerizes and complements Gly356Ser and Arg375Cys enzymes I proving that the association/dissociation properties of enzyme I are a function of the C-terminal domain.     

The primary structure of pig liver thioltransferase

The complete amino acid sequence of pig liver thioltransferase has been determined. The homogeneous protein was cleaved by trypsin, chymotrypsin, Staphylococcus aureus V8 protease, and cyanogen bromide. The resulting peptides were purified by reversed-phase high performance liquid chromatography and ion-exchange fast protein liquid chromatography. Sequencing of the fragments was achieved with either automated Edman degradation or fast atom bombardment-mass spectrometry. Pig liver thioltransferase is a single polypeptide with 105 amino acid residues and an acetylated glutamine N terminus. The protein has 2 cysteine pairs with sequences of -Cys-Pro-Phe-Cys- and -Cys-Ile-Gly-Gly-Cys-, the first pair of which (Cys22 and Cys25) is located at the potential active site of the enzyme. The sequence of pig liver thioltransferase displays close homology (82%) with calf thymus glutaredoxin, suggesting that they belong to the same evolutionary family.     

Analyzing the catalytic mechanism of protein tyrosine phosphatase PtpB from Staphylococcus aureus through site-directed mutagenesis

Protein tyrosine phosphatase B (PtpB) from Staphylococcus aureus, MRSA 252, is a low molecular weight protein tyrosine phosphatase involved in its pathogenicity. PtpB has been modeled in silico and site-directed mutagenesis performed to ascertain the importance of active site residues Cys8, Arg14, Ser15 and Asp120 in its catalytic mechanism. Kinetic characterization of wild-type and the mutant PtpBs, C8S, R14A, S15T, S15A, D120A, D120E, D120N revealed the reaction mechanism followed by this LMWPTPase. The mutations caused major changes in the local environment resulting in significant decrease of its catalytic activity. Inhibition kinetics for the wild-type enzyme was performed with maleimide and maleimidobutyric acid.     

The region from phenylalanine-17 to phenylalanine-28 of a yeast mitochondrial ATPase inhibitor is essential for its ATPase inhibitory activity.

Mitochondrial ATP synthase (F(1)F(o)-ATPase) is regulated by an intrinsic ATPase inhibitor protein. In the present study, we investigated the structure-function relationship of the yeast ATPase inhibitor by amino acid replacement. A total of 22 mutants were isolated and characterized. Five mutants (F17S, R20G, R22G, E25A, and F28S) were entirely inactive, indicating that the residues, Phe17, Arg20, Arg22, Glu25, and Phe28, are essential for the ATPase inhibitory activity of the protein. The activity of 7 mutants (A23G, R30G, R32G, Q36G, L37G, L40S, and L44G) decreased, indicating that the residues, Ala23, Arg30, Arg32, Gln36, Leu37, Leu40, and Leu44, are also involved in the activity. Three mutants, V29G, K34Q, and K41Q, retained normal activity at pH 6.5, but were less active at pH 7.2, indicating that the residues, Val29, Lys34, and Lys41, are required for the protein's action at higher pH. The effects of 6 mutants (D26A, E35V, H39N, H39R, K46Q, and K49Q) were slight or undetectable, and the residues Asp26, Glu35, His39, Lys46, and Lys49 thus appear to be dispensable. The mutant E21A retained normal ATPase inhibitory activity but lacked pH-sensitivity. Competition experiments suggested that the 5 inactivated mutants (F17S, R20G, R22G, E25A, and F28S) could still bind to the inhibitory site on F(1)F(o)-ATPase. These results show that the region from the position 17 to 28 of the yeast inhibitor is the most important for its activity and is required for the inhibition of F(1), rather than binding to the enzyme.     

In vivo formation and stability of engineered disulfide bonds in subtilisin

Computer modeling suggested that a disulfide bond could be built into Bacillus amyloliquefaciens subtilisin between positions 22 (wild-type, Thr) and 87 (Ser) or between positions 24 (Ser) and 87 (Ser). Single cysteines were introduced into this cysteine-free protease at positions 22, 24, or 87 by site-directed mutagenesis of the cloned subtilisin gene. The corresponding double-cysteine mutants were constructed, and recombinant plasmids were expressed in Bacillus subtilis. Double-cysteine mutant enzymes were secreted as efficiently as wild-type, and disulfide bonds were formed quantitatively in vivo. These disulfide bonds were introduced approximately 24 A away from the catalytic site and had no detectable effect on either the specific activities or the pH optima of the mutant enzymes. The equilibrium constants for the reduction of the mutant disulfide bonds by dithiothreitol were determined to be 82 +/- 22 and 20 +/- 5 for Cys22/Cys87 and Cys24/Cys87, respectively. Studies of autoproteolytic inactivation of wild-type subtilisin support a relationship between autolytic stability and conformational stability of the protein. The stabilities of Cys24/Cys87 and wild-type enzymes to autolysis were essentially the same; however, Cys22/Cys87 was actually less stable to autolysis. Reduction of the disulfide cross-bridge lowered the autolytic stability of both double-cysteine mutants relative to their disulfide forms. This correlates with a lowered autolytic stability for the Cys22 and Cys87 single-cysteine mutants, and the fact that an intramolecular hydrogen bond between the hydroxyl groups of Thr22 and Ser87 is likely to be disrupted in the Cys22 and Cys87 single-cysteine mutant proteins.     

Catalytic mechanism of the glutathione peroxidase-type tryparedoxin peroxidase of Trypanosoma brucei

Trypanosoma brucei, the causative agent of African sleeping sickness, encodes three nearly identical genes for cysteine-homologues of the selenocysteine-containing glutathione peroxidases. The enzymes, which are essential for the parasites, lack glutathione peroxidase activity but catalyse the trypanothione/Tpx (tryparedoxin)-dependent reduction of hydroperoxides. Cys47, Gln82 and Trp137 correspond to the selenocysteine, glutamine and tryptophan catalytic triad of the mammalian selenoenzymes. Site-directed mutagenesis revealed that Cys47 and Gln82 are essential. A glycine mutant of Trp137 had 13% of wild-type activity, which suggests that the aromatic residue may play a structural role but is not directly involved in catalysis. Cys95, which is conserved in related yeast and plant proteins but not in the mammalian selenoenzymes, proved to be essential as well. In contrast, replacement of the highly conserved Cys76 by a serine residue resulted in a fully active enzyme species and its role remains unknown. Thr50, proposed to stabilize the thiolate anion at Cys47, is also not essential for catalysis. Treatment of the C76S/C95S but not of the C47S/C76S double mutant with H2O2 induced formation of a sulfinic acid and covalent homodimers in accordance with Cys47 being the peroxidative active site thiol. In the wild-type peroxidase, these oxidations are prevented by formation of an intramolecular disulfide bridge between Cys47 and Cys95. As shown by MS, regeneration of the reduced enzyme by Tpx involves a transient mixed disulfide between Cys95 of the peroxidase and Cys40 of Tpx. The catalytic mechanism of the Tpx peroxidase resembles that of atypical 2-Cys-peroxiredoxins but is distinct from that of the selenoenzymes.     

The short‐chain oxidoreductase Q9HYA2 from Pseudomonas aeruginosa PAO1 contains an atypical catalytic center

The characteristic oxidation or reduction reaction mechanisms of short‐chain oxidoreductase (SCOR) enzymes involve a highly conserved Asp‐Ser‐Tyr‐Lys catalytic tetrad. The SCOR enzyme Q9HYA2 from the pathogenic bacterium Pseudomonas aeruginosa was recognized to possess an atypical catalytic tetrad composed of Lys118‐Ser146‐Thr159‐Arg163. Orthologs of Q9HYA2 containing the unusual catalytic tetrad along with conserved substrate and cofactor recognition residues were identified in 27 additional species, the majority of which are bacterial pathogens. However, this atypical catalytic tetrad was not represented within the Protein Data Bank. The crystal structures of unligated and NADPH‐complexed Q9HYA2 were determined at 2.3 Å resolution. Structural alignment to a polyketide ketoreductase (KR), a typical SCOR, demonstrated that Q9HYA2's Lys118, Ser146, and Arg163 superimposed upon the KR's catalytic Asp114, Ser144, and Lys161, respectively. However, only the backbone of Q9HYA2's Thr159 overlapped KR's catalytic Tyr157. The Thr159 hydroxyl in apo Q9HYA2 is poorly positioned for participating in catalysis. In the Q9HYA2–NADPH complex, the Thr159 side chain was modeled in two alternate rotamers, one of which is positioned to interact with other members of the tetrad and the bound cofactor. A chloride ion is bound at the position normally occupied by the catalytic tyrosine hydroxyl. The putative active site of Q9HYA2 contains a chemical moiety at each catalytically important position of a typical SCOR enzyme. This is the first observation of a SCOR protein with this alternate catalytic center that includes threonine replacing the catalytic tyrosine and an ion replacing the hydroxyl moiety of the catalytic tyrosine.     

Immunological characterization of thioltransferase from pig liver

Polyclonal antibodies against pig liver thioltransferase were raised in a New Zealand rabbit. These antibodies completely neutralized the thioltransferase activity of the homogeneous enzyme and that in the crude cytosolic homogenate at an equivalent titer. The antibodies also cross-reacted equally with calf thymus glutaredoxin and calf liver thioltransferase, but not with Escherichia coli thioredoxin, suggesting that thioltransferase and glutaredoxin from the same species are identical. Immunoblotting analysis of the cytosolic proteins from 14 different pig tissues revealed that most pig tissues contain a 12-kDa protein which reacts with these antibodies. This protein is found in greater abundance in stomach, small intestine, liver, skeletal muscle, kidney, heart, lung, and cerebral cortex, whereas retina, cerebellum, spleen, pancreas, and thymus have low levels of the protein. No reactive protein was detected in the lens. The tissue distribution of the protein was also determined by assay of the enzyme activity and was generally in good agreement with that obtained from the immunoblotting survey. Pig liver thioltransferase was cleaved by trypsin, chymotrypsin, Staphylococcus aureus V8 protease, and cyanogen bromide. The selected peptides purified by reversed phase high performance liquid chromatography or ion exchange fast protein liquid chromatography were subjected to reaction with the polyclonal antibodies against pig liver thioltransferase. Four antigenically reactive fragments were detected by dot-blotting analysis. These peptides are located in the first 30-amino acid residues from the NH2 terminus and the sequence from amino acid residues 39-67, indicating that the active site of the enzyme, Cys22 and Cys25, is located on one of the antigenic determinant domains.     

Identification and reactivity of the catalytic site of pig liver thioltransferase

The active site cysteine of pig liver thioltransferase was identified as Cys22. The kinetics of the reaction between Cys22 of the reduced enzyme and iodoacetic acid as a function of pH revealed that the active site sulfhydryl group had a pKa of 2.5. Incubation of reduced enzyme with [1-14C]cysteine prevented the inactivation of the enzyme by iodoacetic acid at pH 6.5, and no stable protein-cysteine disulfide was found when the enzyme was separated from excess [1-14C]cysteine, suggesting an intramolecular disulfide formation. The results suggested a reaction mechanism for thioltransferase. The thiolated Cys22 first initiates a nucleophilic attack on a disulfide substrate, resulting in the formation of an unstable mixed disulfide between Cys22 and the substrate. Subsequently, the sulfhydryl group at Cys25 is deprotonated as a result of micro-environmental changes within the active site domain, releasing the mixed disulfide and forming an intramolecular disulfide bond. Reduced glutathione, the second substrate, reduces the intramolecular disulfide forming a transient mixed disulfide which is then further reduced by glutathione to regenerate the reduced enzyme and form oxidized glutathione. The rate-limiting step for a typical reaction between a disulfide and reduced glutathione is proposed to be the reduction of the intramolecular disulfide form of the enzyme by reduced glutathione.     

Functional analysis of CYP2D6.31 variant: homology modeling suggests possible disruption of redox partner interaction by Arg440His substitution

Cytochrome P450 2D6 (CYP2D6) is an important human drug-metabolizing enzyme that exhibits a marked genetic polymorphism. Numerous CYP2D6 alleles have been characterized at a functional level, although the consequences for expression and/or catalytic activity of a substantial number of rare variants remain to be investigated. One such allele, CYP2D6*31, is characterized by mutations encoding three amino acid substitutions: Arg296Cys, Arg440His and Ser486Thr. The identification of this allele in an individual with an apparent in vivo poor metabolizer phenotype prompted us to analyze the functional consequence of these substitutions on enzyme activity using yeast as a heterologous expression system. We demonstrated that the Arg440His substitution, alone or in combination with Arg296Cys and/or Ser486Thr, altered the respective kinetic parameters [Km (microM) and kcat (min(-1))] of debrisoquine 4-hydroxylation (wild-type, 25; 0.92; variants, 43-68; 0.05-0.11) and dextromethorphan O-demethylation (wild-type, 1; 4.72; variants, 12-23; 0.64-1.43), such that their specificity constants (kcat/Km) were decreased by more than 95% compared to those observed with the wild-type enzyme. The rates of oxidation of rac-metoprolol at single substrate concentrations of 40 and 400 microM were also markedly decreased by approximately 90% with each CYP2D6 variant containing the Arg440His substitution. These in vitro data confirm that the CYP2D6*31 allele encodes an enzyme with a severely impaired but residual catalytic activity and, furthermore, that the Arg440His exchange alone is the inactivating mutation. A homology model of CYP2D6 based on the crystal structure of rabbit CYP2C5 locates Arg440 on the proximal surface of the protein. Docking the structure of the FMN domain of human cytochrome P450 reductase to the CYP2D6 model suggests that Arg440 is a key member of a cluster of basic amino acid residues important for reductase binding.     

Identification of MAPKAP kinase 2 as a major enzyme responsible for the phosphorylation of the small mammalian heat shock proteins.

MAP kinase-activated protein kinase-2 (MAPKAP kinase-2) phosphorylates the serine residues in murine heat shock protein 25 (hsp25) and human heat shock protein 27 (hsp27) which are phosphorylated in vivo in response to growth factors and heat shock, namely Ser15 and Ser86 (hsp25) and Ser15, Ser78 and Ser82 (hsp27). Ser86 of hsp25 and the equivalent residue in hsp27 (Ser82) are phosphorylated preferentially in vitro. The small heat shock protein is present in rabbit skeletal muscle and hsp25 kinase activity in skeletal muscle extracts co-purifies with MAPKAP kinase-2 activity throughout the purification of the latter enzyme. These results suggest that MAPKAP kinase-2 is the enzyme responsible for the phosphorylation of these small heat shock proteins in mammalian cells.     

Enhancement of the thermal stability of pyroglutamyl peptidase I by introduction of an intersubunit disulfide bond.

From the comparison of the three-dimensional structure of mesophilic pyroglutamyl peptidase from Bacillus amyloliquefaciens and the thermophilic enzyme from Thermococcus litoralis, the intersubunit disulfide bond was estimated to be one of the factors for thermal stability. Since Ser185 was corresponded to Cys190 of the thermophilic enzyme by sequence alignment, the Ser185 residue was replaced with cysteine by site-directed mutagenesis. The S185C mutant enzyme appeared to form a disulfide bond, which was confirmed by SDS-PAGE with and without 2-mercaptoethanol. The mutant enzyme showed a catalytic efficiency equivalent to that of the wild-type enzyme for hydrolysis of a synthetic peptide substrate. However, the thermal stability of the S185C mutant was found to be 30 degrees C higher than that of wild-type. Thus the introduction of a disulfide bond enhanced thermal stability without changing the catalytic efficiency of the enzyme.     

Structural and functional consequences of the replacement of proximal residues Cys(172) and Cys(192) in the large subunit of ribulose-1,5-bisphosphate carboxylase/oxygenase from Chlamydomonas reinhardtii

Proximal Cys(172) and Cys(192) in the large subunit of the photosynthetic enzyme Rubisco (ribulose-1,5-bisphosphate carboxylase/oxygenase; EC 4.1.1.39) are evolutionarily conserved among cyanobacteria, algae and higher plants. Mutation of Cys(172) has been shown to affect the redox properties of Rubisco in vitro and to delay the degradation of the enzyme in vivo under stress conditions. Here, we report the effect of the replacement of Cys(172) and Cys(192) by serine on the catalytic properties, thermostability and three-dimensional structure of Chlamydomonas reinhardtii Rubisco. The most striking effect of the C172S substitution was an 11% increase in the specificity factor when compared with the wild-type enzyme. The specificity factor of C192S Rubisco was not altered. The V(c) (V(max) for carboxylation) was similar to that of wild-type Rubisco in the case of the C172S enzyme, but approx. 30% lower for the C192S Rubisco. In contrast, the K(m) for CO(2) and O(2) was similar for C192S and wild-type enzymes, but distinctly higher (approximately double) for the C172S enzyme. C172S Rubisco showed a critical denaturation temperature approx. 2 degrees C lower than wild-type Rubisco and a distinctly higher denaturation rate at 55 degrees C, whereas C192S Rubisco was only slightly more sensitive to temperature denaturation than the wild-type enzyme. X-ray crystal structures reveal that the C172S mutation causes a shift of the main-chain backbone atoms of beta-strand 1 of the alpha/beta-barrel affecting a number of amino acid side chains. This may cause the exceptional catalytic features of C172S. In contrast, the C192S mutation does not produce similar structural perturbations.     

枯草蛋白酶E的突变体

Ser236位于横贯枯草蛋白酶E的α螺旋末端,远离催化活性中心,Ser236的突变不会对酶的活性产生大的影响。用定点突变的方法对枯草蛋白酶E的基因进行改造引入Ser236Cys,可能会形成分子间二硫键,有利于提高酶的稳定性。Ser236Cys变体酶(BP1)活性是野生型蛋白酶E的15倍,热稳定性提高3倍;进一步在其他位点引入突变的变体酶BU1(A1a15Asp/Gly20His/Ser236Cys)和BW1(Ser24His/Lys27Asp/Ser236Cys)活性都比野生型蛋白酶E低,但BW1的稳定性稍高于野生型蛋白酶E。     

Increasing the thermal stability of the water-soluble pyrroloquinoline quinone glucose dehydrogenase by single amino acid replacement

Based on the characterization of a PCR mutation of water-soluble glucose dehydrogenase possessing pyrroloquinoline quinone (PQQ), PQQGDH-B, Ser231Cys, we have constructed a series of Ser231 variants. The replacement of Ser231 to Cys, Met, Leu, Asp, Asn, His, or Lys resulted in an increase in thermal stability. Among these variants, Ser231Lys showed the highest level of thermal stability and also showed high catalytic activity. Considering that Ser231Lys showed more than an 8-fold increase in its half-life during the thermal inactivation at 55 degrees C compared with the wild-type enzyme, and also retained catalytic activity similar to a wild-type enzyme, the application of this mutant enzyme as a glucose sensor constituent may develop into a stable glucose sensor construction.     

An unusually low pK(a) for Cys282 in the active site of human muscle creatine kinase

All phosphagen kinases contain a conserved cysteine residue which has been shown by crystallographic studies, on both creatine kinase and arginine kinase, to be located in the active site. There are conflicting reports as to whether this cysteine is essential for catalysis. In this study we have used site-directed mutagenesis to replace Cys282 of human muscle creatine kinase with serine and methionine. In addition, we have replaced Cys282, conserved across all creatine kinases, with alanine. No activity was found with the C282M mutant. The C282S mutant showed significant, albeit greatly reduced, activity in both the forward (creatine phosphorylation) and reverse (MgADP phosphorylation) reactions. The K(m) for creatine was increased approximately 10-fold, but the K(m) for phosphocreatine was relatively unaffected. The V and V/K pH-profiles for the wild-type enzyme were similar to those reported for rabbit muscle creatine kinase, the most widely studied creatine kinase isozyme. However, the V/K(creatine) profile for the C282S mutant was missing a pK of 5.4. This suggests that Cys282 exists as the thiolate anion, and is necessary for the optimal binding of creatine. The low pK of Cys282 was also determined spectrophotometrically and found to be 5.6 +/- 0.1. The S284A mutant was found to have reduced catalytic activity, as well as a 15-fold increase in K(m) for creatine. The pK(a) of Cys282 in this mutant was found to be 6.7 +/- 0.1, indicating that H-bonding to Ser284 is an important, but not the sole, factor contributing to the unusually low pK(a) of Cys282.     

Computational and mutational analysis of human glutaredoxin (thioltransferase): probing the molecular basis of the low pKa of cysteine 22 and its role in catalysis

Human glutaredoxin (GRx), also known as thioltransferase, is a 12 kDa thiol-disulfide oxidoreductase that is highly selective for reduction of glutathione-containing mixed disulfides. The apparent pK(a) for the active site Cys22 residue is approximately 3.5. Previously we observed that the catalytic enhancement by glutaredoxin could be ascribed fully to the difference between the pK(a) of its Cys22 thiol moiety and the pK(a) of the product thiol, each acting as a leaving group in the enzymatic and nonenzymatic reactions, respectively [Srinivasan et al. (1997), Biochemistry 36, 3199-3206]. Continuum electrostatic calculations suggest that the low pK(a) of Cys22 results primarily from stabilization of the thiolate anion by a specific ion-pairing with the positively charged Lys19 residue, although hydrogen bonding interactions with Thr21 also appear to contribute. Variants of Lys19 were considered to further assess the predicted role of Lys19 on the pK(a) of Cys22. The variants K19Q and K19L were generated by molecular modeling, and the pK(a) value for Cys22 was calculated for each variant. For K19Q, the predicted Cys22 pK(a) is 7.3, while the predicted value is 8.3 for K19L. The effects of the mutations on the interaction energy between the adducted glutathionyl moiety and GRx were roughly estimated from the van der Waals and electrostatic energies between the glutathionyl moiety and proximal protein residues in a mixed disulfide adduct of GRx and glutathione, i.e., the GRx-SSG intermediate. The values for the K19 mutants differed by only a small amount compared to those for the wild type enzyme intermediate. Together, the computational analysis predicted that the mutant enzymes would have markedly reduced catalytic rates while retaining the glutathionyl specificity displayed by the wild type enzyme. Accordingly, we constructed and characterized the K19L and K19Q mutants of two forms of the GRx enzyme. Each of the mutants retained glutathionyl specificity as predicted and displayed diminution in activity, but the decreases in activity were not to the extent predicted by the theoretical calculations. Changes in the respective Cys22-thiol pK(a) values of the mutant enzymes, as shown by pH profiles for iodoacetamide inactivation of the respective enzymes, clearly revealed that the K19-C22 ion pair cannot fully account for the low pK(a) of the Cys22 thiol. Additional contributions to stabilization of the Cys22 thiolate are likely donated by Thr21 and the N-terminal partial positive charge of the neighboring alpha-helix.