Singlet-triplet gaps in large multireference systems: Spin-flip-driven alternatives for bioinorganic modeling

Quantum mechanical calculations suggest that lytic polysaccharide monooxygenases use a copper-oxyl, oxygen-rebound mechanism

Lytic polysaccharide monooxygenases (LPMOs) exhibit a mononuclear copper-containing active site and use dioxygen and a reducing agent to oxidatively cleave glycosidic linkages in polysaccharides. LPMOs represent a unique paradigm in carbohydrate turnover and exhibit synergy with hydrolytic enzymes in biomass depolymerization. To date, several features of copper binding to LPMOs have been elucidated, but the identity of the reactive oxygen species and the key steps in the oxidative mechanism have not been elucidated. Here, density functional theory calculations are used with an enzyme active site model to identify the reactive oxygen species and compare two hypothesized reaction pathways in LPMOs for hydrogen abstraction and polysaccharide hydroxylation; namely, a mechanism that employs a η(1)-superoxo intermediate, which abstracts a substrate hydrogen and a hydroperoxo species is responsible for substrate hydroxylation, and a mechanism wherein a copper-oxyl radical abstracts a hydrogen and subsequently hydroxylates the substrate via an oxygen-rebound mechanism. The results predict that oxygen binds end-on (η(1)) to copper, and that a copper-oxyl-mediated, oxygen-rebound mechanism is energetically preferred. The N-terminal histidine methylation is also examined, which is thought to modify the structure and reactivity of the enzyme. Density functional theory calculations suggest that this posttranslational modification has only a minor effect on the LPMO active site structure or reactivity for the examined steps. Overall, this study suggests the steps in the LPMO mechanism for oxidative cleavage of glycosidic bonds.

Transmission coefficients for chemical reactions with multiple states: role of quantum decoherence

Transition-state theory (TST) is a widely accepted paradigm for rationalizing the kinetics of chemical reactions involving one potential energy surface (PES). Multiple PES reaction rate constants can also be estimated within semiclassical approaches provided the hopping probability between the quantum states is taken into account when determining the transmission coefficient. In the Marcus theory of electron transfer, this hopping probability was historically calculated with models such as Landau-Zener theory. Although the hopping probability is intimately related to the question of the transition from the fully quantum to the semiclassical description, this issue is not adequately handled in physicochemical models commonly in use. In particular, quantum nuclear effects such as decoherence or dephasing are not present in the rate constant expressions. Retaining the convenient semiclassical picture, we include these effects through the introduction of a phenomenological quantum decoherence function. A simple modification to the usual TST rate constant expression is proposed: in addition to the electronic coupling, a characteristic decoherence time τ(dec) now also appears as a key parameter of the rate constant. This new parameter captures the idea that molecular systems, although intrinsically obeying quantum mechanical laws, behave semiclassically after a finite but nonzero amount of time (τ(dec)). This new degree of freedom allows a fresh look at the underlying physics of chemical reactions involving more than one quantum state. The ability of the proposed formula to describe the main physical lines of the phenomenon is confirmed by comparison with results obtained from density functional theory molecular dynamics simulations for a triplet to singlet transition within a copper dioxygen adduct relevant to the question of dioxygen activation by copper monooxygenases.

Study of the docking of competitive inhibitors at a model of tyrosinase active site: insights from joint broken-symmetry/Spin-Flip DFT computations and ELF topological analysis

Following our previous study (Piquemal et al., New J. Chem., 2003, 27, 909), we present here a DFT study of the inhibition of the Tyrosinase enzyme. Broken-symmetry DFT computations are supplemented with Spin-Flip TD-DFT calculations, which, for the first time, are applied to such a dicopper enzyme. The chosen biomimetic model encompasses a dioxygen molecule, two Cu(II) cations, and six imidazole rings. The docking energy of a natural substrate, namely phenolate, together with those of several inhibitor and non-inhibitor compounds, are reported and show the ability of the model to rank the most potent inhibitors in agreement with experimental data. With respect to broken-symmetry calculations, the Spin-Flip TD-DFT approach reinforces the possibility for theory to point out potent inhibitors: the need for the deprotonation of the substrates, natural or inhibitors, is now clearly established. Moreover, Electron Localization Function (ELF) topological analysis computations are used to deeply track the particular electronic distribution of the Cu-O-Cu three-center bonds involved in the enzymatic Cu(2)O(2) metallic core (Piquemal and Pilmé, J. Mol. Struct.: Theochem, 2006, 77, 764). It is shown that such bonds exhibit very resilient out-of-plane density expansions that play a key role in docking interactions: their 3D-orientation could be the topological electronic signature of oxygen activation within such systems.

Density functional theory for transition metals and transition metal chemistry

We introduce density functional theory and review recent progress in its application to transition metal chemistry. Topics covered include local, meta, hybrid, hybrid meta, and range-separated functionals, band theory, software, validation tests, and applications to spin states, magnetic exchange coupling, spectra, structure, reactivity, and catalysis, including molecules, clusters, nanoparticles, surfaces, and solids.

Long Distance Electron-Transfer Mechanism in Peptidylglycine α-Hydroxylating Monooxygenase: A Perfect Fitting for a Water Bridge

The active sites of copper enzymes have been the subject of many theoretical and experimental investigations from a number of years. Such studies have embraced topics devoted to the modeling of the first coordination sphere at the metallic cations up to the development of biomimetic, or bioinspired, catalytic systems. At least from the theoretical viewpoint, fewer efforts have been dedicated to elucidate how the two copper cations act concertedly in noncoupled dicopper enzymes such as peptidylglycine alpha-hydroxylating monooxygenase (PHM) and dopamine beta-monooxygenase (DbetaM). In these metalloenzymes, an electronic transfer is assumed between the two distant copper cations (11 A). Recent experimental results suggest that this transfer occurs through water molecules, a phenomenon which has been theoretically evidenced to be of high efficiency in the case of cytochrome b5 (Science, 2005, 310, 1311). In the present contribution dedicated to PHM, we overpass the common theoretical approaches dedicated to the electronic and geometrical structures of sites CuM or CuH restricted to their first coordination spheres and aim at directly comparing theoretical results to the experimentally measured activity of the PHM enzyme. To achieve this goal, molecular dynamics simulations were performed on wild-type and various mutants of PHM. More precisely, we provide an estimate of the electron-transfer efficiency between the CuM and CuH sites by means of such molecular dynamics simulations coupled to Marcus theory joined to the Beratan model to approximate the required coupling matrix elements. The theoretical results are compared to the kinetics measurements performed on wild and mutated PHM. The present work, the dynamic aspects of which are essential, accounts for the experimental results issued from mutagenesis. It supports the conclusion that an electronic transfer can occur between two copper(I) sites along a bridge involving a set of hydrogen and chemical bonds. Residue Gln170 is evidenced to be the keystone of this water-mediated pathway.

Mononuclear Cu-O2 complexes: geometries, spectroscopic properties, electronic structures, and reactivity

Using interwoven experimental and theoretical methods, detailed studies of several structurally defined 1:1 Cu-O 2 complexes have provided important fundamental chemical information useful for understanding the nature of intermediates involved in aerobic oxidations in synthetic and enzymatic copper-mediated catalysis. In particular, these studies have shed new light on the factors that influence the mode of O 2 coordination (end-on vs side-on) and the electronic structure, which can vary between Cu(II)-superoxo and Cu(III)-peroxo extremes.