The Activation of Electrophile, Nucleophile and Leaving Group during the Reaction Catalysed by pI258 Arsenate Reductase

The reduction of arsenate to arsenite by pI258 arsenate reductase (ArsC) combines a nucleophilic displacement reaction with a unique intramolecular disulfide cascade. Within this reaction mechanism, the oxidative equivalents are translocated from the active site to the surface of ArsC. The first reaction step in the reduction of arsenate by pI258 ArsC consists of a nucleophilic displacement reaction carried out by Cys10 on dianionic arsenate. The second step involves the nucleophilic attack of Cys82 on the Cys10–arseno intermediate formed during the first reaction step. The onset of the second step is studied here by using quantum chemical calculations in a density functional theory context. The optimised geometry of the Cys10–arseno adduct in the ArsC catalytic site (sequence motif: Cys10–Thr11–Gly12–Asn13–Ser14–Cys15–Arg16–Ser17) forms the starting point for all subsequent calculations. Thermodynamic data and a hard and soft acids and bases (HSAB) reactivity analysis show a preferential nucleophilic attack on a monoanionic Cys10–arseno adduct, which is stabilised by Ser17. The P‐loop active site of pI258 ArsC activates first a hydroxy group and subsequently arsenite as the leaving group, as is clear from an increase in the calculated nucleofugality of these groups upon going from the gas phase to the solvent phase to the enzymatic environment. Furthermore, the enzymatic environment stabilises the thiolate form of the nucleophile Cys82 by 3.3 pH units through the presence of the eight‐residue α helix flanked by Cys82 and Cys89 (redox helix) and through a hydrogen bond with Thr11. The importance of Thr11 in the pKa regulation of Cys82 was confirmed by the observed decrease in the kcat value of the Thr11Ala mutant as compared to that of wild‐type ArsC. During the final reaction step, Cys89 is activated as a nucleophile by structural alterations of the redox helix that functions as a pKa control switch for Cys89; this final step is necessary to expose a Cys82–Cys89 disulfide.

[1]  Eaton E Lattman,et al.  Experimental pK(a) values of buried residues: analysis with continuum methods and role of water penetration. , 2002, Biophysical journal.

[2]  K. Morokuma,et al.  Effects of the protein environment on the structure and energetics of active sites of metalloenzymes. ONIOM study of methane monooxygenase and ribonucleotide reductase. , 2002, Journal of the American Chemical Society.

[3]  Giovanni Scalmani,et al.  New developments in the polarizable continuum model for quantum mechanical and classical calculations on molecules in solution , 2002 .

[4]  R. Parr,et al.  Hardness, softness, and the fukui function in the electronic theory of metals and catalysis. , 1985, Proceedings of the National Academy of Sciences of the United States of America.

[5]  L. Kiessling,et al.  Reactivity of a 2‘-Thio Nucleotide Analog , 1996 .

[6]  Rudolph Willem,et al.  Arsenate reductase from S. aureus plasmid pI258 is a phosphatase drafted for redox duty , 2001, Nature Structural Biology.

[7]  J. Steyaert,et al.  A Nucleophile Activation Dyad in Ribonucleases , 2002, The Journal of Biological Chemistry.

[8]  Goedele Roos,et al.  A Computational and Conceptual DFT Study on the Michaelis Complex of pI258 Arsenate Reductase. Structural Aspects and Activation of the Electrophile and Nucleophile , 2004 .

[9]  Jacopo Tomasi,et al.  Fast Evaluation of Geometries and Properties of Excited Molecules in Solution: A Tamm-Dancoff Model with Application to 4-Dimethylaminobenzonitrile , 2000 .

[10]  B. Edwards,et al.  Insights into the structure, solvation, and mechanism of ArsC arsenate reductase, a novel arsenic detoxification enzyme. , 2001, Structure.

[11]  R. Parr Density-functional theory of atoms and molecules , 1989 .

[12]  G. N. Ramachandhan Need for nonplanar peptide units in polypeptide chains. , 1968 .

[13]  J. Tomasi,et al.  Electrostatic interaction of a solute with a continuum. A direct utilizaion of AB initio molecular potentials for the prevision of solvent effects , 1981 .

[14]  E. Fauman,et al.  Crystal Structure of the Catalytic Domain of the Human Cell Cycle Control Phosphatase, Cdc25A , 1998, Cell.

[15]  K. Morokuma,et al.  ONIOM: A Multilayered Integrated MO + MM Method for Geometry Optimizations and Single Point Energy Predictions. A Test for Diels−Alder Reactions and Pt(P(t-Bu)3)2 + H2 Oxidative Addition , 1996 .

[16]  A. Warshel,et al.  What are the dielectric “constants” of proteins and how to validate electrostatic models? , 2001, Proteins.

[17]  Kenneth B. Wiberg,et al.  Application of the pople-santry-segal CNDO method to the cyclopropylcarbinyl and cyclobutyl cation and to bicyclobutane , 1968 .

[18]  Thom Vreven,et al.  Geometry optimization with QM/MM, ONIOM, and other combined methods. I. Microiterations and constraints , 2003, J. Comput. Chem..

[19]  K. Morokuma,et al.  Prediction of the dissociation energy of hexaphenylethane using the ONIOM(MO :MO:MO) method , 2002 .

[20]  P. Geerlings,et al.  Conceptual density functional theory. , 2003, Chemical reviews.

[21]  T. Koopmans,et al.  Über die Zuordnung von Wellenfunktionen und Eigenwerten zu den Einzelnen Elektronen Eines Atoms , 1934 .

[22]  P. Geerlings,et al.  MECHANISM OF 2 + 1 CYCLOADDITIONS OF HYDROGEN ISOCYANIDE TO ALKYNES : MOLECULAR ORBITAL AND DENSITY FUNCTIONAL THEORY STUDY , 1999 .

[23]  J. Gázquez The Hard and Soft Acids and Bases Principle , 1997 .

[24]  Ulf Ryde,et al.  Comparison of methods for deriving atomic charges from the electrostatic potential and moments , 1998 .

[25]  A. Fersht,et al.  Histidine residues at the N- and C-termini of alpha-helices: perturbed pKas and protein stability. , 1992, Biochemistry.

[26]  R. Parr,et al.  Density-functional theory of the electronic structure of molecules. , 1995, Annual review of physical chemistry.

[27]  Jacopo Tomasi,et al.  Continuum solvation models: A new approach to the problem of solute’s charge distribution and cavity boundaries , 1997 .

[28]  K. Morokuma,et al.  A NEW ONIOM IMPLEMENTATION IN GAUSSIAN98. PART I. THE CALCULATION OF ENERGIES, GRADIENTS, VIBRATIONAL FREQUENCIES AND ELECTRIC FIELD DERIVATIVES , 1999 .

[29]  Goedele Roos,et al.  Origin of the pKa Perturbation of N-Terminal Cysteine in α- and 310-Helices: A Computational DFT Study , 2006 .

[30]  Ralph G. Pearson,et al.  Chemical Hardness: PEARSON:CHEM.HARDNESS O-BK , 1997 .

[31]  R. Zahler Enzyme Structure and Mechanism , 1979, The Yale Journal of Biology and Medicine.

[32]  K. Hallenga,et al.  Solution structure of the low‐molecular‐weight protein tyrosine phosphatase from Tritrichomonas foetus reveals a flexible phosphate binding loop , 2005, Protein science : a publication of the Protein Society.

[33]  J. Steyaert,et al.  Functional assessment of “in vivo” and “in silico” mutations in the guanine binding site of RNase T1: A DFT study , 2004 .

[34]  L. Curtiss,et al.  Intermolecular interactions from a natural bond orbital, donor-acceptor viewpoint , 1988 .

[35]  L. Wyns,et al.  Kinetics and active site dynamics of Staphylococcusaureus arsenate reductase , 2001, JBIC Journal of Biological Inorganic Chemistry.

[36]  James S. M. Anderson,et al.  Perturbative perspectives on the chemical reaction prediction problem , 2005 .

[37]  V. Barone,et al.  Recent advances in density functional methods , 1995 .

[38]  L. Wyns,et al.  Purification of an oxidation-sensitive enzyme, pI258 arsenate reductase from Staphylococcus aureus. , 2003, Journal of chromatography. B, Analytical technologies in the biomedical and life sciences.

[39]  Gilles Klopman,et al.  Chemical reactivity and reaction paths , 1974 .

[40]  Paul Geerlings,et al.  Leaving group activation by aromatic stacking: an alternative to general acid catalysis. , 2004, Journal of molecular biology.

[41]  James S. M. Anderson,et al.  Indices for predicting the quality of leaving groups. , 2005, Physical chemistry chemical physics : PCCP.

[42]  H. Berendsen,et al.  The α-helix dipole and the properties of proteins , 1978, Nature.

[43]  C. Breneman,et al.  Determining atom‐centered monopoles from molecular electrostatic potentials. The need for high sampling density in formamide conformational analysis , 1990 .

[44]  Carla Mattos,et al.  Structural mechanism of oxidative regulation of the phosphatase Cdc25B via an intramolecular disulfide bond. , 2005, Biochemistry.

[45]  P. Geerlings,et al.  HSAB principle: Applications of its global and local forms in organic chemistry , 2000 .

[46]  P A Kollman,et al.  Electrostatic potentials of proteins. 1. Carboxypeptidase A. , 1976, Journal of the American Chemical Society.

[47]  Rudolph Willem,et al.  All intermediates of the arsenate reductase mechanism, including an intramolecular dynamic disulfide cascade , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[48]  Robert G. Parr,et al.  Density functional approach to the frontier-electron theory of chemical reactivity , 1984 .

[49]  K D Watenpaugh,et al.  Crystal structure of the catalytic subunit of Cdc25B required for G2/M phase transition of the cell cycle. , 1999, Journal of molecular biology.

[50]  Goedele Roos,et al.  A Computational and Conceptual DFT Study of the Reactivity of Anionic Compounds: Implications for Enzymatic Catalysis , 2003 .

[51]  B. Edwards,et al.  Arginine 60 in the ArsC arsenate reductase of E. coli plasmid R773 determines the chemical nature of the bound As(III) product , 2004, Protein science : a publication of the Protein Society.

[52]  E. Fauman,et al.  A ligand‐induced conformational change in the yersinia protein tyrosine phosphatase , 1995, Protein science : a publication of the Protein Society.

[53]  Henry Chermette,et al.  Chemical reactivity indexes in density functional theory , 1999 .

[54]  G. A. Jeffrey,et al.  An Introduction to Hydrogen Bonding , 1997 .

[55]  Donald Bashford,et al.  Stabilization of Charges and Protonation States in the Active Site of the Protein Tyrosine Phosphatases: A Computational Study† , 2000 .

[56]  P. Geerlings,et al.  Conceptual and computational DFT in the study of aromaticity. , 2001, Chemical reviews.

[57]  F. Momany,et al.  Ab initio studies of molecular geometries. 27. Optimized molecular structures and conformational analysis of N.alpha.-acetyl-N-methylalaninamide and comparison with peptide crystal data and empirical calculations , 1983 .

[58]  J. Murray,et al.  Comparison of quantum chemical parameters and Hammett constants in correlating pK(a) values of substituted anilines. , 2001, The Journal of organic chemistry.