Metal-Dependent Activity of Fe and Ni Acireductone Dioxygenases: How Two Electrons Reroute the Catalytic Pathway

Dioxygenase Chemistry in Nucleophilic Aldehyde Deformylations Utilizing Dicopper O2-Derived Peroxide Complexes

The chemistry of copper-dioxygen complexes is relevant to copper enzymes in biology as well as in (ligand)Cu-O2 (or Cu2-O2) species utilized in oxidative transformations. For overall energy considerations, as applicable in chemical synthesis, it is beneficial to have an appropriate atom economy; both O-atoms of O2(g) are transferred to the product(s). However, examples of such dioxygenase-type chemistry are extremely rare or not well documented. Herein, we report on nucleophilic oxidative aldehyde deformylation reactivity by the peroxo-dicopper(II) species [Cu2II(BPMPO-)(O22-)]1+ {BPMPO-H = 2,6-bis{[(bis(2-pyridylmethyl)amino]methyl}-4-methylphenol)} and [Cu2II(XYLO-)(O22-)]1+ (XYLO- = a BPMPO- analogue possessing bis(2-{2-pyridyl}ethyl)amine chelating arms). Their dicopper(I) precursors are dioxygenase catalysts. The O2(g)-derived peroxo-dicopper(II) intermediates react rapidly with aldehydes like 2-phenylpropionaldehyde (2-PPA) and cyclohexanecarboxaldehyde (CCA) in 2-methyltetrahydrofuran at -90 °C. Warming to room temperature (RT) followed by workup results in good yields of formate (HC(O)O-) along with ketones (acetophenone or cyclohexanone). Mechanistic investigation shows that [Cu2II(BPMPO-)(O22-)]1+ species initially reacts reversibly with the aldehydes to form detectable dicopper(II) peroxyhemiacetal intermediates, for which optical titrations provide the Keq (at -90 °C) of 73.6 × 102 M-1 (2-PPA) and 10.4 × 102 M-1 (CCA). In the reaction of [Cu2II(XYLO-)(O22-)]1+ with 2-PPA, product complexes characterized by single-crystal X-ray crystallography are the anticipated dicopper(I) complex, [Cu2I(XYLO-)]1+ plus a mixed-valent Cu(I)Cu(II)-formate species. Formate was further identified and confirmed by 1H NMR spectroscopy and electrospray ionization mass spectrometry (ESI-MS) analysis. Using 18O2(g)-isotope labeling the reaction produced a high yield of 18-O incorporated acetophenone as well as formate. The overall results signify that true dioxygenase reactions have occurred, supported by a thorough mechanistic investigation.

Unraveling the Pivotal Roles of Various Metal Ion Centers in the Catalysis of Quercetin 2,4-Dioxygenases

Quercetin 2,4-dioxygenase (QueD) with various transition metal ion co-factors shows great differences, but the internal reasons have not been illustrated in detail. In order to explore the effects of metal ion centers on the catalytic reactivity of QueD, we calculated and compared the minimum energy crossing point (MECP) of dioxygen from the relatively stable triplet state to the active singlet state under different conditions by using the DFT method. It was found that the metal ions play a more important role in the activation of dioxygen compared with the substrate and the protein environment. Simultaneously, the catalytic reactions of the bacterial QueDs containing six different transition metal ions were studied by the QM/MM approach, and we finally obtained the reactivity sequence of metal ions, Ni2+ > Co2+ > Zn2+ > Mn2+ > Fe2+ > Cu2+, which is basically consistent with the previous experimental results. Our calculation results indicate that metal ions act as Lewis acids in the reaction to stabilize the substrate anion and the subsequent superoxo and peroxo species in the reaction, and promote the proton coupled electron transfer (PCET) process. Furthermore, the coordination tendencies of transition metal ion centers also have important effects on the catalytic cycle. These findings have general implications on metalloenzymes, which can expand our understanding on how various metal ions play their key role in modulating catalytic reactivity.

Mononuclear nickel(ii)-flavonolate complexes of tetradentate tripodal 4N ligands as structural and functional models for quercetin 2,4-dioxygenase: structures, spectra, redox and dioxygenase activity

Three new nickel(ii)-flavonolate complexes of the type [Ni(L)(fla)](ClO4) 1-3, where L is the tripodal 4N ligand tris(pyrid-2-ylmethyl)amine (tpa, L1) or (pyrid-2-ylmethyl)bis(6-methylpyrid-2-ylmethyl)amine (6-Me2-tpa, L2) or tris(N-Et-benzimidazol-2-ylmethyl)amine (Et-ntb, L3), have been isolated as functional models for Ni(ii)-containing quercetin 2,4-dioxygenase. Single crystal X-ray structures of 1 and 3 reveal that Ni(ii) is involved in π-back bonding with flavonolate (fla-), as evident from enhancement in C[double bond, length as m-dash]O bond length upon coordination [H(fla), 1.232(3); 1, 1.245(7); 3, 1.262(8) Å]. More asymmetric chelation of fla- in 3 than in 1 [Δd = (Ni-Ocarbonyl - Ni-Oenolate): 1, 0.126; 3, 0.182 Å] corresponds to lower π-delocalization in 3 with electron-releasing N-Et substituent. The optimized structures of 1-3 and their geometrical isomers have been computed by DFT methods. The HOMO and LUMO, both localized on Ni(ii)-bound fla-, are highly conjugated bonding π- and antibonding π*-orbitals respectively. They are located higher in energy than the Ni(ii)-based MOs (HOMO-1, dx2-y2; HOMO-2/6, dz2), revealing that the Ni(ii)-bound fla- rather than Ni(ii) would undergo oxidation upon exposure to dioxygen. The results of computational studies, in combination with spectral and electrochemical studies, support the involvement of redox-inactive Ni(ii) in π-back bonding with fla-, tuning the π-delocalization in fla- and hence its activation. Upon exposure to dioxygen, all the flavonolate adducts in DMF solution decompose to produce CO and depside, which then is hydrolyzed to give the corresponding acids at 70 °C. The highest rate of dioxygenase reactivity of 3 (kO2: 3 (29.10 ± 0.16) > 1 (16.67 ± 0.70) > 2 (1.81 ± 0.04 × 10-1 M-1 s-1)), determined by monitoring the disappearance of the LMCT band in the range 440-450 nm, is ascribed to the electron-releasing N-Et substituent on bzim ring, which decreases the π-delocalization in fla- and enhances its activation.

QM/MM Molecular Dynamics Simulations Revealed Catalytic Mechanism of Urease

Urease catalyzes the hydrolysis of urea to form ammonia and carbamate, inducing an overall pH increase that affects both human health and agriculture. Inhibition, mutagenesis, and kinetic studies have provided insights into its enzymatic role, but there have been debates on the substrate binding mode as well as the reaction mechanism. In the present study, we report quatum mechanics-only (QM-only) and quantum mechanics/molecular mechanics molecular dynamics (QM/MM MD) calculations on urease that mainly investigate the binding mode of urea and the mechanism of the urease-catalyzed hydrolysis reaction. Comparison between the experimental data and our QM(GFN2-xTB)/MM metadynamics results demonstrates that urea hydrolysis via a complex with bidentate-bound urea is much more favorable than via that with monodentate-bound urea for both nucleophilic attack and the subsequent proton transfer steps. We also indicate that the bidentate coordination of urea fits the active site with a closed conformation of the mobile flap and can facilitate the stabilization of transition states and intermediates by forming multiple hydrogen bonds with certain active site residues.

Contrasting Effects of Inhibitors Li+ and Be2+ on Catalytic Cycle of Glycogen Synthase Kinase-3β

Ionic lithium shows rare effectiveness for treating bipolar disorder and is a potential drug for neurodegenerative diseases. Unfortunately, lithium suffers from significant drawbacks, mainly a narrow therapeutic window. Among the targets of lithium, glycogen synthase kinase 3β (GSK-3β) may be responsible for its therapeutic effects. The development of alternative, selective inhibitors of this kinase could prevent lithium side effects, but such efforts have met little success so far. An atomistic understanding of Li⁺ inhibition and the GSK-3β phosphorylation reaction would therefore facilitate the development of new drugs. In this study, we use extensive sampling of catalytic states with our mixed quantum-classical dynamics method QM/DMD and binding affinities from a competitive metal affinity (CMA) approach to expand the atomistic picture of Li⁺ GSK-3β inhibition. We compare Li⁺ action with Be²⁺ and find our results in agreement with in vitro kinetics studies. Ultimately, our simulations show that Li⁺ inhibition is driven by decreasing the phosphorylation reaction rate, rather than reducing catalytic turnover through tight binding to different GSK-3β states like Be²⁺ inhibition. The effect of these metals derive from electrostatic differences and especially their smaller atomic radii compared to the native Mg²⁺ and thus provide insight for the development of GSK-3β inhibitors based on other paradigms.

A Model for the Solution Structure of Human Fe(II)-Bound Acireductone Dioxygenase and Interactions with the Regulatory Domain of Matrix Metalloproteinase I (MMP-I)

The metalloenzyme acireductone dioxygenase (ARD) shows metal-dependent physical and enzymatic activities depending upon the metal bound in the active site. The Fe(II)-bound enzyme catalyzes the penultimate step of the methionine salvage pathway, converting 1,2-dihydroxy-5-(methylthio)pent-1-en-3-one (acireductone) into formate and the ketoacid precursor of methionine, 2-keto-4-thiomethyl-2-oxobutanoate, using O2 as the oxidant. If Ni(II) is bound, an off-pathway shunt occurs, producing 3-methylthiopropionate, formate, and carbon monoxide from the same acireductone substrate. The solution structure of the Fe(II)-bound human enzyme, HsARD, is described and compared with the structures of Ni-bound forms of the closely related mouse enzyme, MmARD. Potential rationales for the different reactivities of the two isoforms are discussed. The human enzyme has been found to regulate the activity of matrix metalloproteinase I (MMP-I), which is involved in tumor metastasis, by binding the cytoplasmic transmembrane tail peptide of MMP-I. Nuclear magnetic resonance titration of HsARD with the MMP-I tail peptide permits identification of the peptide binding site on HsARD, a cleft anterior to the metal binding site adjacent to a dynamic proline-rich loop.

A family of structural and functional models for the active site of a unique dioxygenase: Acireductone dioxygenase (ARD)

We report the synthesis and biomimetic activity of a family of model complexes with relevance to acireductone dioxygenase (ARD), an enzyme that displays dual function based on metal identity found in the methionine salvage pathway (MSP). Three complexes with related structural motifs were synthesized and characterized derived from phenolate, and pyridine N₄O Schiff-base ligands. They display pseudo-octahedral Ni(II)-N₄O ligand coordination with water at the sixth site, in close alignment to the structure in the resting state of ARD. The three featured complexes exhibit carbon‑carbon bond cleavage activation of lithium acetylacetonate, which was used as a model enzyme substrate. Computationally derived mechanistic routes for the observed reactivity consistent with experimental conditions are herein proposed. The mechanism suggests the possibility of Ni(II)-substrate interactions, followed by oxygen insertion. These results constitute only the third functional model system of ARD, in an attempt to further advance biomimetic contributions to the ongoing debate of ARD's unique metal mediated, regioselective oxidative cleavage.

Toxic and Physiological Metal Uptake and Release by Human Serum Transferrin

An atomistic understanding of metal transport in the human body is critical to anticipate the side effects of metal-based therapeutics and holds promise for new drugs and drug delivery designs. Human serum transferrin (hTF) is a central part of the transport processes because of its ubiquitous ferrying of physiological Fe(III) and other transition metals to tightly controlled parts of the body. There is an atomistic mechanism for the uptake process with Fe(III), but not for the release process, or for other metals. This study provides initial insight into these processes for a range of transition metals-Ti(IV), Co(III), Fe(III), Ga(III), Cr(III), Fe(II), Zn(II)-through fully atomistic, extensive quantum mechanical/discrete molecular dynamics sampling and provides, to our knowledge, a new technique we developed to calculate relative binding affinities between metal cations and the protein. It identifies protonation of Tyr188 as a trigger for metal release rather than protonation of Lys206 or Lys296. The study identifies the difficulty of metal release from hTF as potentially related to cytotoxicity. Simulations identify a few critical interactions that stabilize the metal binding site in a flexible, nuanced manner.

Theoretical Studies of Nickel-Dependent Enzymes

The advancements of quantum chemical methods and computer power allow detailed mechanistic investigations of metalloenzymes. In particular, both quantum chemical cluster and combined QM/MM approaches have been used, which have been proven to successfully complement experimental studies. This review starts with a brief introduction of nickel-dependent enzymes and then summarizes theoretical studies on the reaction mechanisms of these enzymes, including NiFe hydrogenase, methyl-coenzyme M reductase, nickel CO dehydrogenase, acetyl CoA synthase, acireductone dioxygenase, quercetin 2,4-dioxygenase, urease, lactate racemase, and superoxide dismutase.

Uncovered Dynamic Coupling Resolves the Ambiguous Mechanism of Phenylalanine Hydroxylase Oxygen Binding

Phenylalanine hydroxylase (PAH) is an iron enzyme catalyzing the oxidation of l-Phe to l-Tyr during phenylalanine catabolism. Dysfunction of PAH leads to the debilitating condition phenylketonuria (PKU), which prompted research into the structure and function of PAH over the last 50 years. Despite intensive study, there is no consensus on the atomistic details of the mechanism of O₂ binding and splitting by wild-type (WT) PAH and how it varies with PKU-inducing mutations, Arg158Gln and Glu280Lys. We studied structures involved in a proposed mechanism for the WT and mutants using extensive mixed quantum-classical molecular dynamics simulations. Simulations reveal a previously unobserved dynamic coupling between the active site and the mutation sites, suggesting how they can affect the catalytic performance of PAH. Furthermore, the effect of the coupling on the PAH structure agrees with and expands our understanding of the experimentally observed differences in activity between the WT and mutants.

Revisiting the methionine salvage pathway and its paralogues

Methionine is essential for life. Its chemistry makes it fragile in the presence of oxygen. Aerobic living organisms have selected a salvage pathway (the MSP) that uses dioxygen to regenerate methionine, associated to a ratchet-like step that prevents methionine back degradation. Here, we describe the variation on this theme, developed across the tree of life. Oxygen appeared long after life had developed on Earth. The canonical MSP evolved from ancestors that used both predecessors of ribulose bisphosphate carboxylase oxygenase (RuBisCO) and methanethiol in intermediate steps. We document how these likely promiscuous pathways were also used to metabolize the omnipresent by-products of S-adenosylmethionine radical enzymes as well as the aromatic and isoprene skeleton of quinone electron acceptors.

Predictive methods for computational metalloenzyme redesign - a test case with carboxypeptidase A

Computational metalloenzyme design is a multi-scale problem. It requires treating the metal coordination quantum mechanically, extensive sampling of the protein backbone, and additionally accounting for the polarization of the active site by both the metal cation and the surrounding protein (a phenomenon called electrostatic preorganization). We bring together a combination of theoretical methods that jointly offer these desired qualities: QM/DMD for mixed quantum-classical dynamic sampling, quantum theory of atoms in molecules (QTAIM) for the assessment of electrostatic preorganization, and Density Functional Theory (DFT) for mechanistic studies. Within this suite of principally different methods, there are both complementarity of capabilities and cross-validation. Using these methods, predictions can be made regarding the relative activities of related enzymes, as we show on the native Zn²⁺-dependent carboxypeptidase A (CPA), and its mutant proteins, which are hypothesized to hydrolyze modified substrates. For the native CPA, we replicated the catalytic mechanism and the rate in close agreement with the experiment, giving validity to the QM/DMD predicted structure, the DFT mechanism, and the QTAIM assessment of catalytic activity. For most sequences of the modified substrate and tried CPA mutants, substantially worsened activity is predicted. However, for the substrate mutant that contains Asp instead of Phe at the C-terminus, one CPA mutant exhibits a reasonable activity, as predicted across the theoretical methods. CPA is a well-studied system, and here it serves as a testing ground for the offered methods.

To what extent do structural changes in catalytic metal sites affect enzyme function?

About half of known enzymatic reactions involve metals. Enzymes belonging to the same superfamily often evolve to catalyze different reactions on the same structural scaffold. The work presented here investigates how functional differentiation, within superfamilies that contain metalloenzymes, relates to structural changes at the catalytic metal site. In general, when the catalytic metal site is unchanged across the enzymes of a superfamily, the functional differentiation within the superfamily tends to be low and the mechanism conserved. Conversely, all types of structural changes in the metal binding site are observed for superfamilies with high functional differentiation. Overall, the catalytic role of the metal ions appears to be one of the most conserved features of the enzyme mechanism within metalloenzyme superfamilies. In particular, when the catalytic role of the metal ion does not involve a redox reaction (i.e. there is no exchange of electrons with the substrate), this role is almost always maintained even when the site undergoes significant structural changes. In these enzymes, functional diversification is most often associated with modifications in the surrounding protein matrix, which has changed so much that the enzyme chemistry is significantly altered. On the other hand, in more than 50% of the examples where the metal has a redox role in catalysis, changes at the metal site modify its catalytic role. Further, we find that there are no examples in our dataset where metal sites with a redox role are lost during evolution.

SYNOPSIS

In this paper we investigate how functional diversity within superfamilies of metalloenzymes relates to structural changes at the catalytic metal site. Evolution tends to strictly conserve the metal site. When changes occur, they do not modify the catalytic role of non-redox metals whereas they affect the role of redox-active metals.

On the Structure and Reaction Mechanism of Human Acireductone Dioxygenase

Acireductone dioxygenase (ARD) is an intriguing enzyme from the methionine salvage pathway that is capable of catalysing two different oxidation reactions with the same substrate depending on the type of the metal ion in the active site. To date, the structural information regarding the ARD-acireductone complex is limited and possible reaction mechanisms are still under debate. The results of joint experimental and computational studies undertaken to advance knowledge about ARD are reported. The crystal structure of an ARD from Homo sapiens was determined with selenomethionine. EPR spectroscopy suggested that binding acireductone triggers one protein residue to dissociate from Fe²⁺ , which allows NO (and presumably O₂ ) to bind directly to the metal. Mössbauer spectroscopic data (interpreted with the aid of DFT calculations) was consistent with bidentate binding of acireductone to Fe²⁺ and concomitant dissociation of His88 from the metal. Major features of Fe vibrational spectra obtained for the native enzyme and upon addition of acireductone were reproduced by QM/MM calculations for the proposed models. A computational (QM/MM) study of the reaction mechanisms suggests that Fe²⁺ promotes O-O bond homolysis, which elicits cleavage of the C1-C2 bond of the substrate. Higher M³⁺ /M²⁺ redox potentials of other divalent metals do not support this pathway, and instead the reaction proceeds similarly to the key reaction step in the quercetin 2,3-dioxygenase mechanism.

Insights into the dioxygen activation and catalytic mechanism of the nickel-containing quercetinase

The catalytic mechanism of Ni-QueD^FLA was elucidated by QM/MM calculations, and the different reactivities of nickel and iron were illuminated.

The Metal Drives the Chemistry: Dual Functions of Acireductone Dioxygenase

Acireductone dioxygenase (ARD) from the methionine salvage pathway (MSP) is a unique enzyme that exhibits dual chemistry determined solely by the identity of the divalent transition-metal ion (Fe²⁺ or Ni²⁺) in the active site. The Fe²⁺-containing isozyme catalyzes the on-pathway reaction using substrates 1,2-dihydroxy-3-keto-5-methylthiopent-1-ene (acireductone) and dioxygen to generate formate and the ketoacid precursor of methionine, 2-keto-4-methylthiobutyrate, whereas the Ni²⁺-containing isozyme catalyzes an off-pathway shunt with the same substrates, generating methylthiopropionate, carbon monoxide, and formate. The dual chemistry of ARD was originally discovered in the bacterium Klebsiella oxytoca, but it has recently been shown that mammalian ARD enzymes (mouse and human) are also capable of catalyzing metal-dependent dual chemistry in vitro. This is particularly interesting, since carbon monoxide, one of the products of off-pathway reaction, has been identified as an antiapoptotic molecule in mammals. In addition, several biochemical and genetic studies have indicated an inhibitory role of human ARD in cancer. This comprehensive review describes the biochemical and structural characterization of the ARD family, the proposed experimental and theoretical approaches to establishing mechanisms for the dual chemistry, insights into the mechanism based on comparison with structurally and functionally similar enzymes, and the applications of this research to the field of artificial metalloenzymes and synthetic biology.

From NAD+ to Nickel Pincer Complex: A Significant Cofactor Evolution Presented by Lactate Racemase

Lactate racemase (LarA), a new nickel enzyme discovered recently, catalyzes the racemization between d- and l-lactates with a novel nickel pincer cofactor (Ni-PTTMN) derived from nicotinic acid. In this study, by using DFT and a 200-atom active-site model, LarA is revealed to employ a modified proton-coupled hydride-transfer mechanism in which a hydride is transferred to a cofactor pyridine carbon from the substrate α-carbon along with proton transfer from the substrate hydroxy group to a histidine, and then moved back from the opposite side. Tyr294 and Lys298 provide significant acceleration effects by orientating substrates and stabilizing the negative charge developing at the substrate hydroxy oxygen. The barrier was determined to be 12.0 kcal mol^-1 , which reveals enhanced racemase activity relative to the LarA reaction using NAD⁺ -like cofactors. Compared with NAD⁺ , Ni-PTTMN has a stronger hydride-addition reactivity in moderate and high environmental polarity and may fit perfectly the moderately polar active site of LarA.

Oxygen activation by mononuclear Mn, Co, and Ni centers in biology and synthetic complexes

The active sites of metalloenzymes that catalyze O₂-dependent reactions generally contain iron or copper ions. However, several enzymes are capable of activating O₂ at manganese or nickel centers instead, and a handful of dioxygenases exhibit activity when substituted with cobalt. This minireview summarizes the catalytic properties of oxygenases and oxidases with mononuclear Mn, Co, or Ni active sites, including oxalate-degrading oxidases, catechol dioxygenases, and quercetin dioxygenase. In addition, recent developments in the O₂ reactivity of synthetic Mn, Co, or Ni complexes are described, with an emphasis on the nature of reactive intermediates featuring superoxo-, peroxo-, or oxo-ligands. Collectively, the biochemical and synthetic studies discussed herein reveal the possibilities and limitations of O₂ activation at these three "overlooked" metals.

Metal-Dependent Function of a Mammalian Acireductone Dioxygenase

The two acireductone dioxygenase (ARD) isozymes from the methionine salvage pathway of Klebsiella oxytoca are the only known pair of naturally occurring metalloenzymes with distinct chemical and physical properties determined solely by the identity of the divalent transition metal ion (Fe(2+) or Ni(2+)) in the active site. We now show that this dual chemistry can also occur in mammals. ARD from Mus musculus (MmARD) was studied to relate the metal ion identity and three-dimensional structure to enzyme function. The iron-containing isozyme catalyzes the cleavage of 1,2-dihydroxy-3-keto-5-(thiomethyl)pent-1-ene (acireductone) by O2 to formate and the ketoacid precursor of methionine, which is the penultimate step in methionine salvage. The nickel-bound form of ARD catalyzes an off-pathway reaction resulting in formate, carbon monoxide (CO), and 3-(thiomethyl) propionate. Recombinant MmARD was expressed and purified to obtain a homogeneous enzyme with a single transition metal ion bound. The Fe(2+)-bound protein, which shows about 10-fold higher activity than that of others, catalyzes on-pathway chemistry, whereas the Ni(2+), Co(2+), or Mn(2+) forms exhibit off-pathway chemistry, as has been seen with ARD from Klebsiella. Thermal stability of the isozymes is strongly affected by the metal ion identity, with Ni(2+)-bound MmARD being the most stable, followed by Co(2+) and Fe(2+), and Mn(2+)-bound ARD being the least stable. Ni(2+)- and Co(2+)-bound MmARD were crystallized, and the structures of the two proteins found to be similar. Enzyme-ligand complexes provide insight into substrate binding, metal coordination, and the catalytic mechanism.

Metal-binding promiscuity in artificial metalloenzyme design

This review presents recent examples of metal-binding promiscuity in protein scaffolds and highlights the effect of metal variation on catalytic functionality. Naturally evolved binding sites, as well as unnatural amino acids and cofactors can bind a diverse range of metals, including non-biological transition elements. Computational screening and rational design have been successfully used to create promiscuous binding-sites. Incorporation of non-native metals into proteins expands the catalytic range of transformations catalysed by enzymes and enhances their potential for application in chemicals synthesis.

Mysteries of metals in metalloenzymes

Natural metalloenzymes are often the most proficient catalysts in terms of their activity, selectivity, and ability to operate at mild conditions. However, metalloenzymes are occasionally surprising in their selection of catalytic metals, and in their responses to metal substitution. Indeed, from the isolated standpoint of producing the best catalyst, a chemist designing from first-principles would likely choose a different metal. For example, some enzymes employ a redox active metal where a simple Lewis acid is needed. Such are several hydrolases. In other cases, substitution of a non-native metal leads to radical improvements in reactivity. For example, histone deacetylase 8 naturally operates with Zn(2+) in the active site but becomes much more active with Fe(2+). For β-lactamases, the replacement of the native Zn(2+) with Ni(2+) was suggested to lead to higher activity as predicted computationally. There are also intriguing cases, such as Fe(2+)- and Mn(2+)-dependent ribonucleotide reductases and W(4+)- and Mo(4+)-dependent DMSO reductases, where organisms manage to circumvent the scarcity of one metal (e.g., Fe(2+)) by creating protein structures that utilize another metal (e.g., Mn(2+)) for the catalysis of the same reaction. Naturally, even though both metal forms are active, one of the metals is preferred in every-day life, and the other metal variant remains dormant until an emergency strikes in the cell. These examples lead to certain questions. When are catalytic metals selected purely for electronic or structural reasons, implying that enzymatic catalysis is optimized to its maximum? When are metal selections a manifestation of competing evolutionary pressures, where choices are dictated not just by catalytic efficiency but also by other factors in the cell? In other words, how can enzymes be improved as catalysts merely through the use of common biological building blocks available to cells? Addressing these questions is highly relevant to the enzyme design community, where the goal is to prepare maximally efficient quasi-natural enzymes for the catalysis of reactions that interest humankind. Due to competing evolutionary pressures, many natural enzymes may not have evolved to be ideal catalysts and can be improved for the isolated purpose of catalysis in vitro when the competing factors are removed. The goal of this Account is not to cover all the possible stories but rather to highlight how variable enzymatic catalysis can be. We want to bring up possible factors affecting the evolution of enzyme structure, and the large- and intermediate-scale structural and electronic effects that metals can induce in the protein, and most importantly, the opportunities for optimization of these enzymes for catalysis in vitro.

Cu,Zn-superoxide dismutase without Zn is folded but catalytically inactive

Amyotrophic lateral sclerosis has been linked to the gain of aberrant function of superoxide dismutase, Cu,Zn-SOD1 upon protein misfolding. The mechanism of SOD1 misfolding is thought to involve mutations leading to the loss of Zn, followed by protein unfolding and aggregation. We show that the removal of Zn from SOD1 may not lead to an immediate unfolding but immediately deactivates the enzyme through a combination of subtle structural and electronic effects. Using quantum mechanics/discrete molecular dynamics, we showed that both Zn-less wild-type (WT)-SOD1 and its D124N mutant that does not bind Zn have at least metastable folded states. In those states, the reduction potential of Cu increases, leading to the presence of detectable amounts of Cu(I) instead of Cu(II) in the active site, as confirmed experimentally. The Cu(I) protein cannot participate in the catalytic Cu(I)-Cu(II) cycle. However, even without the full reduction to Cu(I), the Cu site in the Zn-less variants of SOD1 is shown to be catalytically incompetent: unable to bind superoxide in a way comparable to the WT-SOD1. The changes are more radical and different in the D124N Zn-less mutant than in the Zn-less WT-SOD1, suggesting D124N being perhaps not the most adequate model for Zn-less SOD1. Overall, Zn in SOD1 appears to be influencing the Cu site directly by adjusting its reduction potential and geometry. Thus, the role of Zn in SOD1 is not just structural, as was previously thought; it is a vital part of the catalytic machinery.