Coordination chemistry at the molybdenum site of sulfite oxidase: redox-induced structural changes in the cysteine 207 to serine mutant

The redox chemistry of the molybdenum site of the C207S mutant of recombinant human sulfite oxidase has been studied via potentiometric titrations employing both electron paramagnetic resonance (EPR) spectroscopy and X-ray absorption spectroscopy (XAS) as probes of the active site structure. In earlier EXAFS studies, oxidized Cys207Ser enzyme has been shown to possess a novel tri-oxo active site, in which Ser207 does not appear to be a ligand to Mo [George, G. N.; Garrett, R. M.; Prince, R. C.; Rajagopalan, K. V. J. Am. Chem. Soc. 1996, 118, 8588-8592]. Redox titrations show that the active site is modified under reducing conditions to a mono-oxo Mo(IV) species, probably with Ser207 ligated to the metal. The Mo(IV) species can be reoxidized to a mono-oxo Mo(V) species still coordinated to Ser207, which in turn can be further reoxidized to yield the initial tri-oxo Mo(VI) structure with loss of Ser207 ligation.     

Molybdenum site structure of Escherichia coli YedY, a novel bacterial oxidoreductase

We report a structural characterization using X-ray absorption spectroscopy of the molybdenum site of Escherichia coli YedY, a novel oxidoreductase related to be the sulfite oxidase family of molybdenum enzymes. We find that the enzyme can exist in Mo(V) and Mo(IV) oxidation states but cannot be readily oxidized to the Mo(VI) form. Mo(V) YedY has molybdenum coordination similar to that of sulfite oxidase, with one Mo═O at 1.71 ?, three Mo-S at 2.39 ?, and one Mo-OH at 2.09 ?, which elongates to 2.20 ? upon reduction to Mo(IV), indicating Mo-OH(2) coordination. The Mo(V) enzyme also possesses a long Mo-O coordination at 2.64 ?, which may be due to oxygen coordination by Asn-45 O(δ), with Mo-O(δ) approximately trans to the Mo═O group. A comparison with sulfite oxidase indicates that YedY possesses a much more uniform Mo-S coordination, with a maximum permitted deviation of less than 0.05 ?. Our results indicate that the YedY active site shows considerable similarity to but also important differences from that of reduced forms of sulfite oxidase.     

High-resolution EXAFS of the active site of human sulfite oxidase: comparison with density functional theory and X-ray crystallographic results

Much of our knowledge about molybdenum enzymes has originated from EXAFS spectroscopy. This technique provides excellent bond-length accuracy but has only limited bond-length resolution. We have used EXAFS spectroscopy with an extended data range in an attempt to improve bond-length resolution for the molybdenum enzyme sulfite oxidase. The Mo site of sulfite oxidase has two oxygen and three Mo-S ligands (two from cofactor dithiolene plus a cysteine). For the oxidized (Mo(VI)) enzyme, we find that the three Mo-S bond lengths are very similar (within 0.05 A) at 2.41 A, as are the Mo=O ligands at 1.72 A. Density functional theory shows that this is consistent with the proposed active-site structure. The reduced (Mo(IV)) enzyme shows two Mo-S bond lengths at 2.35 A and one at 2.41 A (assigned to cofactor dithiolene and cysteine, respectively, from DFT), together with one Mo=O at 1.72 A and one Mo-OH(2) at 2.30 A.     

The active site of arsenite oxidase from Alcaligenes faecalis

Arsenite oxidase, a member of the DMSO reductase family of molybdenum enzymes, has two molecules of guanosine dinucleotide molybdenum cofactor coordinating the molybdenum at the active site. X-ray absorption spectroscopy indicates that the Mo-S bonds shorten from 2.47 to 2.37 A upon reduction with the physiological substrate. It also indicates the presence of an oxo ligand at 1.70 A in both oxidized and reduced forms of the enzyme, together with a short, 1.83 A, Mo-O bond in the oxidized form that is lost upon reduction. Resonance Raman spectroscopy indicates that the two pterin dithiolene moieties have different aromaticities, with one, the Q-pterin, having a more discrete dithiolate structure while the other, the P-pterin, has considerable pi-delocalization. Our results indicate that the structure of arsenite oxidase is intermediate between that seen in other molybdenum enzymes, in which one ligand to the metal is provided by the polypeptide (serine, cysteine, or selenocysteine), and tungsten enzymes that lack a peptide ligand.     

Spectroscopic studies of molybdenum and tungsten enzymes

Molybdenum and tungsten are the only second and third-row transition elements with a known function in living systems. Molybdenum fulfills functional roles in enzyme systems in almost all living creatures, from bacteria through plants to invertebrates and mammals, while tungsten takes the place of molybdenum in some prokaryotes, especially the hyperthermophilic archaea. The enzymes contain the metal bound by an unusual sulfur-containing cofactor. Despite possessing common structural elements, the enzymes are remarkable in the range of different chemical reactions that are catalyzed, although almost all are two-electron oxidation–reduction reactions in which an oxygen atom is transferred to or from the molybdenum. The functional roles filled by molybdenum enzymes are equally diverse; for example, they play essential roles in microbial respiration, in the uptake of nitrogen in green plants, in controlling insect eye color, and in human health. Spectroscopic studies, in particular electron paramagnetic resonance and X-ray absorption spectroscopy, have played an essential role in our understanding of the active site structures and catalytic mechanisms of the molybdenum and tungsten enzymes. This review summarizes the role spectroscopy has played in the state of our knowledge of the molybdenum and tungsten enzymes, with particular regard to structural information on the molybdenum sites.     

Structure of the active site of sulfite dehydrogenase from Starkeya novella

In this paper, we report the results of molybdenum K-edge X-ray absorption studies performed on the oxidized and reduced active sites of the sulfite dehydrogenase from Starkeya novella. Our results provide the first direct structural information on the active site of the oxidized form of this enzyme and confirm the conclusions derived from protein crystallography that the molybdenum coordination is analogous to that of the sulfite oxidases. The molybdenum atom of the oxidized enzyme is bound by two Mo=O ligands at 1.73 A and three thiolate Mo-S ligands at 2.42 A, whereas the reduced enzyme has one oxo at 1.74 A, one long oxygen at 2.19 A (characteristic of Mo-OH2), and three Mo-S ligands at 2.40 A.     

Rational engineering of a DNA glycosylase specific for an unnatural cytosine:pyrene base pair

A novel site-specific cytosine DNA glycosylase has been rationally engineered from the active site scaffold of the DNA repair enzyme uracil DNA glycosylase (UDG). UDG, which operates by a nucleotide flipping mechanism, was first converted into a sequence nonspecific cytosine DNA glycosylase (CDG) by altering the base-specific hydrogen bond donor-acceptor groups in the active site. A second mutation that renders UDG defective in nucleotide flipping was then introduced, and the double mutant was rescued using a substrate with a "preflipped" cytosine base. Substrate-assisted flipping was engineered by incorporation of an unnatural pyrene nucleotide wedge (Y) into the DNA strand opposite to the target cytosine. This new enzyme, CYDG, can be used to target cleavage of specific cytosine residues in the context of a C/Y base pair in any DNA fragment.     

Origins of enhanced proton transport in the Y7F mutant of human carbonic anhydrase II

Human carbonic anhydrase II (HCA II), among the fastest enzymes known, catalyzes the reversible hydration of CO 2 to HCO 3 (-). The rate-limiting step of this reaction is believed to be the formation of an intramolecular water wire and transfer of a proton across the active site cavity from a zinc-bound solvent to a proton shuttling residue (His64). X-ray crystallographic studies have shown this intramolecular water wire to be directly stabilized through hydrogen bonds via a small well-defined set of amino acids, namely, Tyr7, Asn62, Asn67, Thr199, and Thr200. Furthermore, X-ray crystallographic and kinetic studies have shown that the mutation of tyrosine 7 to phenylalanine, Y7F HCA II, has the effect of increasing the proton transfer rate by 7-fold in the dehydration direction of the enzyme reaction compared to wild-type (WT). This increase in the proton transfer rate is postulated to be linked to the formation of a more directional, less branched, water wire. To evaluate this proposal, molecular dynamics simulations have been employed to study water wire formation in both the WT and Y7F HCA II mutant. These studies reveal that the Y7F mutant enhances the probability of forming small water wires and significantly extends the water wire lifetime, which may account for the elevated proton transfer seen in the Y7F mutant. Correlation analysis of the enzyme and intramolecular water wire indicates that the Y7F mutant significantly alters the interaction of the active site waters with the enzyme while occupancy data of the water oxygens reveals that the Y7F mutant stabilizes the intramolecular water wire in a manner that maximizes smaller water wire formation. This increase in the number of smaller water wires is likely to elevate the catalytic turnover of an already very efficient enzyme.     

X-ray crystal structure and EPR spectra of "arsenite-inhibited" Desulfovibriogigas aldehyde dehydrogenase: a member of the xanthine oxidase family

X-ray crystallography has been used to determine the structure of arsenite-inhibited aldehyde dehydrogenase from Desulfovibrio gigas, a member of the xanthine oxidase family of mononuclear molybdenum enzymes. The structure shows an AsO3 moiety bound to the molybdenum atom of the active site through one of the oxygen atoms. A reduced sample of arsenite-inhibited aldehyde dehydrogenase has a Mo(V) signal that shows anisotropic hyperfine and quadrupole coupling to one arsenic atom. This signal has a strong resemblance with a previously reported signal for arsenite-inhibited xanthine oxidase.     

Unfavorable electrostatic and steric interactions in DNA polymerase β E295K mutant interfere with the enzyme's pathway

Mutations in DNA polymerase β (pol β) have been associated with approximately 30% of human tumors. The E295K mutation of pol β has been linked to gastric carcinoma via interference with base excision repair. To interpret the different behavior of E295K as compared to wild-type pol β in atomic and energetic detail, we resolve a binary crystal complex of E295K at 2.5 ? and apply transition path sampling (TPS) to delineate the closing pathway of the E295K pol β mutant. Conformational changes are important components in the enzymatic pathway that lead to and ready the enzyme for the chemical reaction. Our analyses show that the closing pathway of E295K mutant differs from the wild-type pol β in terms of the individual transition states along the pathway, associated energies, and the active site conformation in the final closed form of the mutant. In particular, the closed state of E295K has a more distorted active site than the active site in the wild-type pol β. In addition, the total energy barrier in the conformational closing pathway is 65 ± 11 kJ/mol, much higher than that estimated for both correct (e.g., G:C) and incorrect (e.g., G:A) wild-type pol β systems (42 ± 8 and 45 ± 7 kJ/mol, respectively). In particular, the rotation of Arg258 is the rate-limiting step in the conformational pathway of E295K due to unfavorable electrostatic and steric interactions. The distorted active site in the closed relative to open state and the high energy barrier in the conformational pathway may explain in part why the E295K mutant is observed to be inactive. Interestingly, however, following the closing of the thumb but prior to the rotation of Arg258, the E295K mutant complex has a similar energy level as compared to the wild-type pol β. This suggests that the E295K mutant may associate with DNA with similar affinity, but it may be hampered in continuing the process of chemistry. Supporting experimental data come from the observation that the catalytic activity of wild-type pol β is hampered when E295K is present: this may arise from the competition between E295K and wild-type enzyme for the DNA. These combined results suggest that the low insertion efficiency of E295K mutant as compared to wild-type pol β may be related to a closed form distorted by unfavorable electrostatic and steric interactions between Arg258 and other key residues. The active site is thus less competent for proceeding to the chemical reaction, which may also involve a higher reaction barrier than the wild-type or may not be possible in this mutant. Our analysis also suggests further experiments for other mutants to test the above hypothesis and dissect the roles of steric and electrostatic factors on enzyme behavior.     

Enzyme dynamics and tunneling enhanced by compression in the hydrogen abstraction catalyzed by soybean lipoxygenase-1

A fully microscopical simulation of the rate-limiting hydrogen abstraction catalyzed by soybean lipoxygenase-1 (SLO-1) has been carried out. This enzyme exhibits the largest, and weakly temperature dependent, experimental H/D kinetic isotope effect (KIE) reported for a biological system. The theoretical model used here includes the complete enzyme with a solvation shell of water molecules, the Fe(III)-OH- cofactor, and the linoleic acid substrate. We have used a hybrid QM(PM3/d-SRP)/MM method to describe the potential energy surface of the whole system, and the ensemble-averaged variational transition-state theory with multidimensional tunneling (EA-VTST/MT) to calculate the rate constant and the primary KIE. The computational results show that the compression of the wild-type active site enzyme results in the huge contribution of tunneling (99%) to the rate of the hydrogen abstraction. Importantly, the active site becomes more flexible in the Ile553Ala mutant reactant complex simulation (for which a markedly temperature dependent KIE has been experimentally determined), thus justifying the proposed key role of the gating promoting mode in the reaction catalyzed by SLO-1. Finally, the results indicate that the calculated KIE for the wild-type enzyme has an important dependence on the barrier width.     

Hydrogen‐bonding interactions in the active site of a low molecular weight protein‐tyrosine phosphatase

Molecular dynamics simulation of the Michaelis complex, phospho‐enzyme intermediate, and the wild‐type and C12S mutant have been carried out to examine hydrogen‐bonding interactions in the active site of the bovine low molecular weight protein‐tyrosine phosphatase (BPTP). It was found that the S_γ atom of the nucleophilic residue Cys‐12 is ideally located at a position opposite from the phenylphosphate dianion for an inline nucleophilic substitution reaction. In addition, electrostatic and hydrogen‐bonding interactions from the backbone amide groups of the phosphate‐binding loop strongly stabilize the thiolate anion, making Cys‐12 ionized in the active site. In the phospho‐enzyme intermediate, three water molecules are found to form strong hydrogen bonds with the phosphate group. In addition, another water molecule can be identified to form bridging hydrogen bonds between the phosphate group and Asp‐129, which may act as the nucleophile in the subsequent phosphate hydrolysis reaction, with Asp‐129 serving as a general base. The structural difference at the active site between the wild‐type and C12S mutant has been examined. It was found that the alkoxide anion is significantly shifted toward one side of the phosphate binding loop, away from the optimal position enjoyed by the thiolate anion of the wild‐type enzyme in an S_N2 process. This, coupled with the high pK_a value of an alcoholic residue, makes the C12S mutant catalytically inactive. These molecular dynamics simulations provided details of hydrogen bonding interactions in the active site of BPTP, and a structural basis for further studies using combined quantum mechanical and molecular mechanical potential to model the entire dephosphorylation reaction by BPTP. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 1192–1203, 2000     

Spectroscopic and kinetic studies of PKU-inducing mutants of phenylalanine hydroxylase: Arg158Gln and Glu280Lys

Phenylalanine hydroxylase (PAH) is a tetrahydrobiopterin-dependent, nonheme iron enzyme that catalyzes the hydroxylation of L-Phe to L-Tyr in the rate-limiting step of phenylalanine catabolism. This reaction is tightly coupled in the wild-type enzyme to oxidation of the tetrahydropterin cofactor. Dysfunction of PAH activity in humans leads to the disease phenylketonuria (PKU). We have investigated two PKU-inducing mutants, Arg158Gln and Glu280Lys, using kinetic methods, magnetic circular dichrosim (MCD) spectroscopy, and X-ray absorption spectroscopy (XAS). Analysis of the products produced by the mutant enzymes shows that although both oxidize pterin at more than twice the rate of wild-type enzyme, these reactions are only approximately 20% coupled to production of L-Tyr. Previous MCD and XAS studies had demonstrated that the resting Fe(II) site is six-coordinate in the wild-type enzyme and converts to a five-coordinate site when both L-Phe and reduced pterin are present in the active site. Although the Arg158Gln mutant forms the five-coordinate site when both cosubstrates are bound, the Fe(II) site of the Glu280Lys mutant remains six-coordinate. These results provide insight into the PAH reaction and disease mechanism at a molecular level, indicating that the first step of the mechanism is formation of a peroxy-pterin species, which subsequently reacts with the Fe(II) site if the pterin is properly oriented for formation of an Fe-OO-pterin bridge and an open coordination position is available on the Fe(II).     

DFT investigation of the molybdenum cofactor in periplasmic nitrate reductases: structure of the Mo(V) EPR-active species

The periplasmic nitrate reductase NAP belongs to the DMSO reductase family that regroups molybdoenzymes housing a bis-molybdopterin cofactor as the active site. Several forms of the Mo(V) state, an intermediate redox state in the catalytic cycle of the enzyme, have been evidenced by EPR spectroscopy under various conditions, but their structure and catalytic relevance are not fully understood. On the basis of structural data available from the literature, we built several models that reproduce the first coordination sphere of the molybdenum cofactor and used DFT methods to make magneto-structural correlations on EPR-detected species. "High-g" states, which are the most abundant Mo(V) species, are characterized by a low-anisotropy g tensor and a high g(min) value. We assign this signature to a six-sulfur coordination sphere in a pseudotrigonal prismatic geometry with a partial disulfide bond. The "very high-g" species is well described with a sulfido ion as the sixth ligand. The "low-g" signal can be successfully associated to a Mo(V) sulfite-oxidase-type active site with only one pterin moiety coordinated to the molybdenum ion with an oxo or sulfido axial ligand. For all these species we investigate their catalytic activity using a thermodynamic point of view on the molybdenum coordination sphere. Beyond the periplasmic nitrate reductase case, this work provides useful magneto-structural correlations to characterize EPR-detected species in mononuclear molybdoenzymes.     

Structure-based redesign of GST A2-2 for enhanced catalytic efficiency with azathioprine

Glutathione transferase (GST) A2-2 is the most efficient human enzyme in the biotransformation of the prodrug azathioprine (Aza). The activation of Aza has therapeutic potential for possible use of GSTs in targeted enzyme-prodrug treatment of diseases. Based on the assumed catalytic mechanism and computational docking of Aza to the active site of the enzyme, active-site residues were selected for construction of focused mutant libraries, which were thereafter screened for Aza activity. Mutants with elevated Aza activity were identified, DNA sequenced, and the proteins purified. The two most active mutants showed up to 70-fold higher catalytic efficiency than the parental GST A2-2. The structure of the most active triple mutant (L107G/L108D/F222H) enzyme was determined by X-ray crystallography demonstrating significant changes in the topography of the active site facilitating productive binding of Aza as a substrate.     

Computational design of a thermostable mutant of cocaine esterase via molecular dynamics simulations

Cocaine esterase (CocE) has been known as the most efficient native enzyme for metabolizing naturally occurring cocaine. A major obstacle to the clinical application of CocE is the thermoinstability of native CocE with a half-life of only ～11 min at physiological temperature (37 °C). It is highly desirable to develop a thermostable mutant of CocE for therapeutic treatment of cocaine overdose and addiction. To establish a structure-thermostability relationship, we carried out molecular dynamics (MD) simulations at 400 K on wild-type CocE and previously known thermostable mutants, demonstrating that the thermostability of the active form of the enzyme correlates with the fluctuation (characterized as the root-mean square deviation and root-mean square fluctuation of atomic positions) of the catalytic residues (Y44, S117, Y118, H287, and D259) in the simulated enzyme. In light of the structure-thermostability correlation, further computational modelling including MD simulations at 400 K predicted that the active site structure of the L169K mutant should be more thermostable. The prediction has been confirmed by wet experimental tests showing that the active form of the L169K mutant had a half-life of 570 min at 37 °C, which is significantly longer than those of the wild-type and previously known thermostable mutants. The encouraging outcome suggests that the high-temperature MD simulations and the structure-thermostability relationship may be considered as a valuable tool for the computational design of thermostable mutants of an enzyme.     

Computational studies of the metal-binding site of the wild-type and the H46R mutant of the copper, zinc superoxide dismutase

Impairment of the Zn(II)-binding site of the copper, zinc superoxide dismutase (CuZnSOD) protein is involved in a number of hypotheses and explanations for the still unknown toxic gain of function mutant varieties of CuZnSOD that are associated with familial forms of amyotrophic lateral sclerosis (ALS). In this work, computational chemistry methods have been used for studying models of the metal-binding site of the ALS-linked H46R mutant of CuZnSOD and of the wild-type variety of the enzyme. By comparing the energy and electronic structure of these models, a plausible explanation for the effect of the H46R mutation on the zinc site is obtained. The computational study clarifies the role of the D124 and D125 residues for keeping the structural integrity of the Zn(II)-binding site, which was known to exist but its mechanism has not been explained. Earlier results suggest that the explanation for the impairment of the Zn(II)-site proposed in this work may be useful for understanding the mechanism of action of the ALS-linked mutations in CuZnSOD in general.     

Enhanced degradation of chemical warfare agents through molecular engineering of the phosphotriesterase active site

The bacterial phosphotriesterase has been utilized as a template for the evolution of improved enzymes for the catalytic decomposition of organophosphate nerve agents. A combinatorial library of active site mutants was constructed by randomizing residues His-254, His-257, and Leu-303. The collection of mutant proteins was screened for the ability to hydrolyze a chromogenic analogue of the most toxic stereoisomer of the chemical warfare agent, soman. The mutant H254G/H257W/L303T catalyzed the hydrolysis of the target substrate nearly 3 orders of magnitude faster than the wild-type enzyme. The X-ray crystal structure was solved in the presence and absence of diisopropyl methyl phosphonate. The mutant enzyme was ligated to an additional divalent cation at the active site that was displaced upon the binding of the substrate analogue inhibitor. These studies demonstrate that substantial changes in substrate specificity can be achieved by relatively minor changes to the primary amino acid sequence.     

X-ray crystal structure of arsenite-inhibited xanthine oxidase: μ-sulfido,μ-oxo double bridge between molybdenum and arsenic in the active site

Xanthine oxidoreductase is a molybdenum-containing enzyme that catalyzes the hydroxylation reaction of sp(2)-hybridized carbon centers of a variety of substrates, including purines, aldehydes, and other heterocyclic compounds. The complex of arsenite-inhibited xanthine oxidase has been characterized previously by UV-vis, electron paramagnetic resonance, and X-ray absorption spectroscopy (XAS), and the catalytically essential sulfido ligand of the square-pyrimidal molybdenum center has been suggested to be involved in arsenite binding through either a μ-sulfido,μ-oxo double bridge or a single μ-sulfido bridge. However, this is contrary to the crystallographically observed single μ-oxo bridge between molybdenum and arsenic in the desulfo form of aldehyde oxidoreductase from Desulfovibrio gigas (an enzyme closely related to xanthine oxidase), whose molybdenum center has an oxo ligand replacing the catalytically essential sulfur, as seen in the functional form of xanthine oxidase. Here we use X-ray crystallography to characterize the molybdenum center of arsenite-inhibited xanthine oxidase and solve the structures of the oxidized and reduced inhibition complexes at 1.82 and 2.11 ? resolution, respectively. We observe μ-sulfido,μ-oxo double bridges between molybdenum and arsenic in the active sites of both complexes. Arsenic is four-coordinate with a distorted trigonal-pyramidal geometry in the oxidized complex and three-coordinate with a distorted trigonal-planar geometry in the reduced complex. The doubly bridged binding mode is in agreement with previous XAS data indicating that the catalytically essential sulfur is also essential for the high affinity of reduced xanthine oxidoreductase for arsenite.