Catalytic Cycle of Multicopper Oxidases Studied by Combined Quantum- and Molecular-Mechanical Free-Energy Perturbation Methods

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, isomerisation 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

Reorganization Energy for Internal Electron Transfer in Multicopper Oxidases.

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

Stability mechanisms of a thermophilic laccase probed by molecular dynamics.

Laccases are highly stable, industrially important enzymes capable of oxidizing a large range of substrates. Causes for their stability are, as for other proteins, poorly understood. In this work, multiple-seed molecular dynamics (MD) was applied to a Trametes versicolor laccase in response to variable ionic strengths, temperatures, and glycosylation status. Near-physiological conditions provided excellent agreement with the crystal structure (average RMSD ∼0.92 Å) and residual agreement with experimental B-factors. The persistence of backbone hydrogen bonds was identified as a key descriptor of structural response to environment, whereas solvent-accessibility, radius of gyration, and fluctuations were only locally relevant. Backbone hydrogen bonds decreased systematically with temperature in all simulations (∼9 per 50 K), probing structural changes associated with enthalpy-entropy compensation. Approaching T opt (∼350 K) from 300 K, this change correlated with a beginning "unzipping" of critical β-sheets. 0 M ionic strength triggered partial denucleation of the C-terminal (known experimentally to be sensitive) at 400 K, suggesting a general salt stabilization effect. In contrast, F(-) (but not Cl(-)) specifically impaired secondary structure by formation of strong hydrogen bonds with backbone NH, providing a mechanism for experimentally observed small anion destabilization, potentially remedied by site-directed mutagenesis at critical intrusion sites. N-glycosylation was found to support structural integrity by increasing persistent backbone hydrogen bonds by ∼4 across simulations, mainly via prevention of F(-) intrusion. Hydrogen-bond loss in distinct loop regions and ends of critical β-sheets suggest potential strategies for laboratory optimization of these industrially important enzymes

Effect of metal Ions (Ni2+, Cu2+ and Zn2+) and water coordination on the structure of L-phenylalanine, L-tyrosine, L-tryptophan and their zwitterionic forms

Methods of quantum chemistry have been applied to double-charged complexes involving the transition metals Ni2+, Cu2+ and Zn2+ with the aromatic amino acids (AAA) phenylalanine, tyrosine and tryptophan. The effect of hydration on the relative stability and geometry of the individual species studied has been evaluated within the supermolecule approach. The interaction enthalpies, entropies and Gibbs energies of nine complexes Phe•M, Tyr•M, Trp•M, (M = Ni2+, Cu2+ and Zn2+) were determined at the Becke3LYP density functional level of theory. Of the transition metals studied the bivalent copper cation forms the strongest complexes with AAAs. For Ni2+and Cu2+ the most stable species are the NO coordinated cations in the AAA metal complexes, Zn2+cation prefers a binding to the aromatic part of the AAA (complex II). Some complexes energetically unfavored in the gas-phase are stabilized upon microsolvation

Theoretical Description of the Structure and Magnetic Properties of Nitroxide−Cu(II)−Nitroxide Spin Triads by Means of Multiconfigurational Ab Initio Calculations

The structural, electronic and magnetic properties of two different models of the heterospin polymer chain complexes of Cu2+ hexafluoroacetylacetonate with two pyrazole-substituted nitronyl nitroxides Cu(hfac)2LR have been studied by means of multiconfigurational perturbation theory based on a CASSCF (complete active space self-consistent field) wave function, i.e. the CASPT2 method. Our calculations reveal the presence of two minima in the electronic energy curve along the Cu−OL bond, separated by only 6 kcal/mol, and corresponding to the X-ray structures of the CuO6 centers in Cu(hfac)2LPr at 115 and 293 K, respectively. The two energetic minima are characterized by a different electronic structure, thus giving rise to a different three-spin exchange coupling and explaining the thermally induced spin transitions in this family of compounds. The concomitant variations in the magnetic properties, i.e. g factors and magnetic moments μeff(T) were calculated and compared with the experimental data of Cu(hfac)2LPr. Even if the correspondence is only qualitative, our calculations provide a convincing explanation of the observed magnetic peculiarities. In particular, at low temperatures, the predicted ground-state is 2Au, well separated from the 2Ag, 4Au states and therefore exclusively populated. Its calculated g factors, g∥ = 1.848, g⊥ = 1.965, 1.974, qualitatively correspond to the observed g N−O•−Cu−•O−N< for these linear spin triads is challenged by our calculations, pointing instead to a more important role of the end-standing NO in the exchange interactions with Cu(II)