Search

Contains fulltext : 201415.pdf (publisher's version ) (Open Access)

We have used combined quantum mechanical and molecular mechanical free-energy perturbation methods in combination with explicit solvent simulations to study the reaction mechanism of the multicopper oxidases, in particular, the regeneration of the reduced state from the native intermediate. For 52 putative states of the trinuclear copper cluster, differing in the oxidation states of the copper ions and the protonation states of water- and O2-derived ligands, we have studied redox potentials, acidity constants, isomerization reactions, as well as water- and O2 binding reactions. Thereby, we can propose a full reaction mechanism of the multicopper oxidases with atomic detail. We also show that the two copper sites in the protein communicate so that redox potentials and acidity constants of one site are affected by up to 0.2 V or 3 pKa units by a change in the oxidation state of the other site.

The catalytic cycle of multicopper oxidases (MCOs) involves intramolecular electron transfer (IET) from the Cu-T1 copper ion, which is the primary site of the one-electron oxidations of the substrate, to the trinuclear copper cluster (TNC), which is the site of the four-electron reduction of dioxygen to water. In this study we report a detailed characterization of the kinetic and electrochemical properties of bilirubin oxidase (BOx) a member of the MCO family. The experimental results strongly indicate that under certain conditions, e.g. in alkaline solutions, the IET can be the rate-limiting step in the BOx catalytic cycle. The data also suggest that one of the catalytically relevant intermediates (most likely characterized by an intermediate oxidation state of the TNC) formed during the catalytic cycle of BOx has a redox potential close to 0.4 V, indicating an uphill IET process from the T1 copper site (0.7 V) to the Cu-T23. These suggestions are supported by calculations of the IET rate, based on the experimentally observed Gibbs free energy change and theoretical estimates of reorganization energy obtained by combined quantum and molecular mechanical (QM/MM) calculations.

The key step in the reaction mechanism of multicopper oxidases (MCOs)-the cleavage of the O-O bond in O-2-has been investigated using combined quantum mechanical and molecular mechanical (QM/MM) methods. This process represents a reaction pathway from the peroxy intermediate after it accepts one electron from the nearby type-1 Cu site to the experimentally-observed native intermediate, which is the only fully oxidised catalytically relevant state in MCOs. Scans of the QM(DFT)/MM potential energy surface have allowed us to obtain estimates of the activation energies. Furthermore, vacuum calculations on a smaller model of the active site have allowed us to estimate the entropy contributions to the barrier height and to obtain further insight into the reaction by comparing the small cluster model with the QM/MM model, which includes the entire protein. Owing to the complicated electronic structure of these low-spin exchange coupled systems, multireference quantum chemical calculations at the complete-active space second-order perturbation theory (CASPT2) were used in an attempt to benchmark the barrier heights obtained at the DFT(B3LYP) level. Our best estimate of the activation barrier is Delta G = 60-65 kJ mol(-1), in good agreement with the experimental barrier of similar to 55 kJ mol(-1), which can be inferred from the experimental rate constant of k > 350 s(-1). It has also been shown that the reaction involves protonation of the O-2 moiety before bond cleavage. The proton likely comes from a nearby carboxylate residue which was recently suggested by the experiments.

We have calculated the reorganization energy for the intramolecular electron transfer between the reduced type 1 copper site and the peroxy intermediate of the trinuclear cluster in the multicopper oxidase CueO. The calculations are performed at the combined quantum mechanics and molecular mechanics (QM/MM) level, based on molecular dynamics simulations with tailored potentials for the two copper sites. We obtain a reorganization energy of 91-133 kJ/mol, depending on the theoretical treatment. The two Cu sites contribute by 12 and 22 kJ/mol to this energy, whereas the solvent contribution is 34 kJ/mol. The rest comes from the protein, involving small contributions from many residues. We have also estimated the energy difference between the two electron-transfer states and show that the reduction of the peroxy intermediate is exergonic by 43-87 kJ/mol, depending on the theoretical method. Both the solvent and the protein contribute to this energy difference, especially charged residues close to the two Cu sites. We compare these estimates with energies obtained from QM/MM optimizations and QM calculations in a vacuum and discuss differences between the results obtained at various levels of theory.

Manganese superoxide dismutases (MnSODs) are enzymes that convert two molecules of the poisonous superoxide radical into molecular oxygen and hydrogen peroxide. During the reaction, the manganese ion cycles between the Mn2+ and Mn3+ oxidation states and accomplishes its enzymatic action in two half-cycles (corresponding to the oxidation and reduction of O-2(center dot-)). Despite many experimental and theoretical studies dealing with SODs, including quantum chemical active-site-model studies of numerous variants of the reaction mechanisms, several details of MnSOD enzymatic action are still unclear. In this study, we have modeled and compared four reaction pathways (one associative, one dissociative, and two second-sphere) in a protein environment using the QM/MM approach (combined quantum and molecular mechanics calculations) at the density functional theory level. The results were complemented by CASSCF/CASPT2/MM single-point energy calculations for the most plausible models to account properly for the multireference character of the various spin multiplets. The results indicate that the oxidation of O-2(center dot-) to O-2 most likely occurs by an associative mechanism following a two-state (quartet-octet) reaction profile. The barrier height is estimated to be less than 25 kJ.mol(-1). On the other hand, the conversion of O-2(center dot-) to H2O2 is likely to take place by a second-sphere mechanism, that is, without direct coordination of the superoxide radical to the manganese center. The reaction pathway involves the conical intersection of two quintet states, giving rise to an activation barrier of similar to 60 kJ.mol(-1). The calculations also indicate that the associative mechanism can represent a competitive pathway in the second half-reaction with the overall activation barrier being only slightly higher than the activation barrier in the second-sphere mechanism. The activation barriers along the proposed reaction pathways are in very good agreement with

We have developed a computational method that combines extended X-ray absorption fine structure (EXAFS) refinements with the integrated quantum mechanical and molecular mechanics (QM/MM) method. This method allows us to obtain a structure of a metal site inside a protein that is compatible with both EXAFS data and QM calculations (i.e., that is chemically reasonable). Thereby, the QM/MM calculations play the same role as MM in nearly all NMR and crystallographic refinements-EXAFS ensures that the metal-ligand distances are accurate and QM/MM fills in all the other structural data. We have used this method to show that a structure with a peroxide ion in the center of the trinuclear cluster fits experimental EXAFS data better than a structure with the peroxide ion on the side of the cluster for the peroxide adduct of multicopper oxidases.

We have studied the detailed reaction mechanism of iron and manganese superoxide dismutase with density functional calculations on realistic active-site models, with large basis sets and including solvation, zero-point, and thermal effects. The results indicate that the conversion of O-2(-) to O-2 follows an associative mechanism, with O-2 directly binding to the metal, followed by the protonation of the metal-bound hydroxide ion, and the dissociation of O-3(2). All these reaction steps are exergonic. Likewise, we suggest that the conversion of O-2(-) to H2O2 follows an at least a partly second-sphere pathway. There are small differences in the preferred oxidation and spin states, as well as in the geometries, of Fe and Mn, but these differences have little influence on the energetics, and therefore on the reaction mechanism of the two types of superoxide dismutases. For example, the two metals have very similar reduction potentials in the active-site models, although they differ by 0.7 V in water solution. The reaction mechanisms and spin states seem to have been designed to avoid spin conversions or to facilitate them by employing nearly degenerate spin states. (c) 2006 Wiley Periodicals, Inc.

The four-electron reduction of dioxygen to water in multicopper oxidases takes place in a trinuclear copper cluster, which is linked to a mononuclear blue copper site, where the substrates are oxidized. Recently, several intermediates in the catalytic cycle have been spectroscopically characterized, and two possible structural models have been suggested for both the peroxy and native intermediates, In this study, these spectroscopic results are complemented by hybrid quantum and molecular mechanical (QM/MM) calculations, taking advantage of recently available crystal structures with a full complement of copper ions. Thereby, we obtain optimized molecular structures for all of the Experimentally studied intermediates involved in the reductive cleavage of the O-2 molecule and energy profiles for individual reaction steps, This allows identification of the experimentally observed intermediates and further insight into the reaction mechanism that is probably relevant for the whole class of multicopper oxidases, We suggest that the peroxy intermediate contains an O-2(2-) ion, in which one oxygen atom bridges the type 2 copper ion and one of the type 3 copper ions, whereas the other one coordinates to the other type 3 copper ion, One-electron reduction of this intermediate triggers the cleavage of the O-O bond, which involves the uptake of a proton, The product of this cleavage is the observed native intermediate, which we suggest to contain a O-2 ion coordinated to all three of the copper ions in the center of the cluster.