Computational and Mutational Analysis of Human Glutaredoxin (Thioltransferase): Probing the Molecular Basis of the Low pK sub(a) 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 sub(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 sub(a) of its Cys22 thiol moiety and the pK sub(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 sub(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 sub(a) of Cys22. The variants K19Q and K19L were generated by molecular modeling, and the pK sub(a) value for Cys22 was calculated for each variant. For K19Q, the predicted Cys22 pK sub(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 sub(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 sub(a) of the Cys22 thiol. Additional contributions to stabilization of the Cys22 thiolate are lik