N–O bond cleavage mechanism(s) in nitrous oxide reductase

Coordination chemistry of the CuZ site in nitrous oxide reductase and its synthetic mimics

Atmospheric nitrous oxide (N2O) has garnered significant attention recently due to its dual roles as an ozone depletion agent and a potent greenhouse gas. Anthropogenic N2O emissions occur primarily through agricultural disruption of nitrogen homeostasis causing N2O to build up in the atmosphere. The enzyme responsible for N2O fixation within the geochemical nitrogen cycle is nitrous oxide reductase (N2OR), which catalyzes 2H+/2e- reduction of N2O to N2 and H2O at a tetranuclear active site, CuZ. In this review, the coordination chemistry of CuZ is reviewed. Recent advances in the understanding of biological CuZ coordination chemistry is discussed, as are significant breakthroughs in synthetic modeling of CuZ that have emerged in recent years. The latter topic includes both structurally faithful, synthetic [Cu4(µ4-S)] clusters that are able to reduce N2O, as well as dicopper motifs that shed light on reaction pathways available to the critical CuI-CuIV cluster edge of CuZ. Collectively, these advances in metalloenzyme studies and synthetic model systems provide meaningful knowledge about the physiologically relevant coordination chemistry of CuZ but also open new questions that will pose challenges in the near future.

Monophenolase and catecholase activity of Aspergillus oryzae catechol oxidase: insights from hybrid QM/MM calculations

Catechol oxidase from Aspergillus oryzae (AoCO4) can not only catalyze oxidation of o-diphenols to o-quinones, but can also catalyze monooxygenation of small phenolics. To gain insight into the catecholase and monophenolase activities of AoCO4, the reaction mechanism of catechol oxidation was investigated by means of hybrid quantum mechanical/molecular mechanical (QM/MM) calculations. The oxy-form of AoCO4 was found to be a μ-η2:η2 side-on peroxo dicopper(ii) complex, which can undergo a proton coupled electron transfer from the substrate rather than a proton transfer from the nearby Ser302 residue to generate a hydroperoxide. The μ-1,1-OOH Cu2(i,ii) complex is thermodynamically more stable than the μ-η1:η2 hydroperoxide. Moreover, the cleavage of the O-O bond in the μ-1,1-OOH Cu2(i,ii) intermediate has a much lower barrier than that in the μ-η1:η2 hydroperoxide species. In both cases, the O-O bond cleavage is the rate-limiting step, generating the reactive (μ-O˙)(μ-OH) dicopper(ii) complex. In addition, our results demonstrated that the oxidation of catechol to quinone is much more preferred than the hydroxylation reaction. These findings may provide useful information for understanding the reactivity of the Cu2O2 active site of coupled binuclear copper enzymes.

Biological and Bioinspired Inorganic N-N Bond-Forming Reactions

The metallobiochemistry underlying the formation of the inorganic N-N-bond-containing molecules nitrous oxide (N2O), dinitrogen (N2), and hydrazine (N2H4) is essential to the lifestyles of diverse organisms. Similar reactions hold promise as means to use N-based fuels as alternative carbon-free energy sources. This review discusses research efforts to understand the mechanisms underlying biological N-N bond formation in primary metabolism and how the associated reactions are tied to energy transduction and organismal survival. These efforts comprise studies of both natural and engineered metalloenzymes as well as synthetic model complexes.

Theory Demonstrated a "Coupled" Mechanism for O2 Activation and Substrate Hydroxylation by Binuclear Copper Monooxygenases

Multiscale simulations have been performed to address the longstanding issue of "dioxygen activation" by the binuclear copper monooxygenases (PHM and DβM), which have been traditionally classified as "noncoupled" binuclear copper enzymes. Our QM/MM calculations rule out that Cu_M(II)-O₂^• is an active species for H-abstraction from the substrate. In contrast, Cu_M(II)-O₂^• would abstract an H atom from the cosubstrate ascorbate to form a Cu_M(II)-OOH intermediate in PHM and DβM. Consistent with the recently reported structural features of DβM, the umbrella sampling shows that the "open" conformation of the Cu_M(II)-OOH intermediate could readily transform into the "closed" conformation in PHM, in which we located a mixed-valent μ-hydroperoxodicopper(I,II) intermediate, (μ-OOH)Cu(I)Cu(II). The subsequent O-O cleavage and OH moiety migration to Cu_H generate the unexpected species (μ-O^•)(μ-OH)Cu(II)Cu(II), which is revealed to be the reactive intermediate responsible for substrate hydroxylation. We also demonstrate that the flexible Met ligand is favorable for O-O cleavage reactions, while the replacement of Met with the strongly bound His ligand would inhibit the O-O cleavage reactivity. As such, the study not only demonstrates a "coupled" mechanism for O₂ activation by binuclear copper monooxygenases but also deciphers the full catalytic cycle of PHM and DβM in accord with the available experimental data. These findings of O₂ activation and substrate hydroxylation by binuclear copper monooxygenases could expand our understanding of the reactivities of the synthetic monocopper complexes.

Mechanism of NO–CO reaction over highly dispersed cuprous oxide on γ-alumina catalyst using a metal–support interfacial site in the presence of oxygen: similarities to and differences from biological systems

The NO–CO reaction mechanism over the Cu/γ-Al₂O₃ catalyst was elucidated using DFT and a cluster model.

The catalytic cycle of nitrous oxide reductase - The enzyme that catalyzes the last step of denitrification

The reduction of the potent greenhouse gas nitrous oxide requires a catalyst to overcome the large activation energy barrier of this reaction. Its biological decomposition to the inert dinitrogen can be accomplished by denitrifiers through nitrous oxide reductase, the enzyme that catalyzes the last step of the denitrification, a pathway of the biogeochemical nitrogen cycle. Nitrous oxide reductase is a multicopper enzyme containing a mixed valence CuA center that can accept electrons from small electron shuttle proteins, triggering electron flow to the catalytic sulfide-bridged tetranuclear copper "CuZ center". This enzyme has been isolated with its catalytic center in two forms, CuZ*(4Cu1S) and CuZ(4Cu2S), proven to be spectroscopic and structurally different. In the last decades, it has been a challenge to characterize the properties of this complex enzyme, due to the different oxidation states observed for each of its centers and the heterogeneity of its preparations. The substrate binding site in those two "CuZ center" forms and which is the active form of the enzyme is still a matter of debate. However, in the last years the application of different spectroscopies, together with theoretical calculations have been useful in answering these questions and in identifying intermediate species of the catalytic cycle. An overview of the spectroscopic, kinetics and structural properties of the two forms of the catalytic "CuZ center" is given here, together with the current knowledge on nitrous oxide reduction mechanism by nitrous oxide reductase and its intermediate species.

Spectroscopic Definition of the CuZ° Intermediate in Turnover of Nitrous Oxide Reductase and Molecular Insight into the Catalytic Mechanism

Spectroscopic methods and density functional theory (DFT) calculations are used to determine the geometric and electronic structure of Cu_Z°, an intermediate form of the Cu₄S active site of nitrous oxide reductase (N₂OR) that is observed in single turnover of fully reduced N₂OR with N₂O. Electron paramagnetic resonance (EPR), absorption, and magnetic circular dichroism (MCD) spectroscopies show that Cu_Z° is a 1-hole (i.e., 3Cu^ICu^II) state with spin density delocalized evenly over Cu_I and Cu_IV. Resonance Raman spectroscopy shows two Cu-S vibrations at 425 and 413 cm^-1, the latter with a -3 cm^-1 O¹⁸ solvent isotope shift. DFT calculations correlated to these spectral features show that Cu_Z° has a terminal hydroxide ligand coordinated to Cu_IV, stabilized by a hydrogen bond to a nearby lysine residue. Cu_Z° can be reduced via electron transfer from Cu_A using a physiologically relevant reductant. We obtain a lower limit on the rate of this intramolecular electron transfer (IET) that is >10⁴ faster than the unobserved IET in the resting state, showing that Cu_Z° is the catalytically relevant oxidized form of N₂OR. Terminal hydroxide coordination to Cu_IV in the Cu_Z° intermediate yields insight into the nature of N₂O binding and reduction, specifying a molecular mechanism in which N₂O coordinates in a μ-1,3 fashion to the fully reduced state, with hydrogen bonding from Lys397, and two electrons are transferred from the fully reduced μ₄S^2- bridged tetranuclear copper cluster to N₂O via a single Cu atom to accomplish N-O bond cleavage.

Theoretically exploring the key role of the Lys412 residue in the conversion of N2O to N2by nitrous oxide reductase from Achromobacter cycloclastes

The active Cu_Zcluster of NOR provides strong back-donation to coordinated N₂O and activates the O atom of the N₂O group facilitating H-bonding and protonationviathe Lys412 residue.

Computational electrochemistry: prediction of liquid-phase reduction potentials

This article reviews recent developments and applications in the area of computational electrochemistry. Our focus is on predicting the reduction potentials of electron transfer and other electrochemical reactions and half-reactions in both aqueous and nonaqueous solutions. Topics covered include various computational protocols that combine quantum mechanical electronic structure methods (such as density functional theory) with implicit-solvent models, explicit-solvent protocols that employ Monte Carlo or molecular dynamics simulations (for example, Car-Parrinello molecular dynamics using the grand canonical ensemble formalism), and the Marcus theory of electronic charge transfer. We also review computational approaches based on empirical relationships between molecular and electronic structure and electron transfer reactivity. The scope of the implicit-solvent protocols is emphasized, and the present status of the theory and future directions are outlined.

Enzyme informatics

Over the last 50 years, sequencing, structural biology and bioinformatics have completely revolutionised biomolecular science, with millions of sequences and tens of thousands of three dimensional structures becoming available. The bioinformatics of enzymes is well served by, mostly free, online databases. BRENDA describes the chemistry, substrate specificity, kinetics, preparation and biological sources of enzymes, while KEGG is valuable for understanding enzymes and metabolic pathways. EzCatDB, SFLD and MACiE are key repositories for data on the chemical mechanisms by which enzymes operate. At the current rate of genome sequencing and manual annotation, human curation will never finish the functional annotation of the ever-expanding list of known enzymes. Hence there is an increasing need for automated annotation, though it is not yet widespread for enzyme data. In contrast, functional ontologies such as the Gene Ontology already profit from automation. Despite our growing understanding of enzyme structure and dynamics, we are only beginning to be able to design novel enzymes. One can now begin to trace the functional evolution of enzymes using phylogenetics. The ability of enzymes to perform secondary functions, albeit relatively inefficiently, gives clues as to how enzyme function evolves. Substrate promiscuity in enzymes is one example of imperfect specificity in protein-ligand interactions. Similarly, most drugs bind to more than one protein target. This may sometimes result in helpful polypharmacology as a drug modulates plural targets, but also often leads to adverse side-effects. Many chemoinformatics approaches can be used to model the interactions between druglike molecules and proteins in silico. We can even use quantum chemical techniques like DFT and QM/MM to compute the structural and energetic course of enzyme catalysed chemical reaction mechanisms, including a full description of bond making and breaking.

[Cu2O]2+ active site formation in Cu-ZSM-5: geometric and electronic structure requirements for N2O activation

Understanding the formation mechanism of the [Cu2O](2+) active site in Cu-ZSM-5 is important for the design of efficient catalysts to selectively convert methane to methanol and related value-added chemicals and for N2O decomposition. Spectroscopically validated DFT calculations are used here to evaluate the thermodynamic and kinetic requirements for formation of [Cu2O](2+) active sites from the reaction between binuclear Cu(I) sites and N2O in the 10-membered rings Cu-ZSM-5. Thermodynamically, the most stable Cu(I) center prefers bidentate coordination with a close to linear bite angle. This binuclear Cu(I) site reacts with N2O to generate the experimentally observed [Cu2O](2+) site. Kinetically, the reaction coordinate was evaluated for two representative binuclear Cu(I) sites. When the Cu-Cu distance is sufficiently short (<4.2 Å), N2O can bind in a "bridged" μ-1,1-O fashion and the oxo-transfer reaction is calculated to proceed with a low activation energy barrier (2 kcal/mol). This is in good agreement with the experimental Ea for N2O activation (2.5 ± 0.5 kcal/mol). However, when the Cu-Cu distance is long (>5.0 Å), N2O binds in a "terminal" η(1)-O fashion to a single Cu(I) site of the dimer and the resulting E(a) for N2O activation is significantly higher (16 kcal/mol). Therefore, bridging N2O between two Cu(I) centers is necessary for its efficient two-electron activation in [Cu2O](2+) active site formation. In nature, this N2O reduction reaction is catalyzed by a tetranuclear CuZ cluster that has a higher E(a). The lower E(a) for Cu-ZSM-5 is attributed to the larger thermodynamic driving force resulting from formation of strong Cu(II)-oxo bonds in the ZSM-5 framework.

Computational modeling of the direct hydride transfer mechanism for the MAO catalyzed oxidation of phenethylamine and benzylamine: ONIOM (QM/QM) calculations

Monoamine oxidases are two isozymic flavoenzymes which are the important targets for drugs used in the treatment of depression, Parkinson and Alzheimer's diseases. The catalytic reaction taking place between the cofactor FAD and amine substrate is still not completely understood. Herein we employed quantum chemical methods on the recently proposed direct hydride transfer mechanism including full active site residues of MAO isoforms in the calculations. Activation free energy barriers of direct hydride transfer mechanism for MAO-A and MAO-B were calculated by ONIOM (our own n-layered integrated molecular orbital + molecular mechanics) method with QM/QM (quantum mechanics:quantum mechanics) approach employing several density functional theory functionals, B3LYP, WB97XD, CAM-B3LYP and M06-2X, for the high layer. The formation of very recently proposed αC-flavin N5 adduct inside the enzyme has been investigated. ONIOM (M06-2X/6-31+G(d,p):PM6) results revealed that such an adduct may form only in MAO-B suggesting slightly different hydride transfer mechanisms for MAO-A and MAO-B.