We have calculated the reorganisation 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 reorganisation 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 optimisations and QM calculations in vacuum, and discuss differences between the results obtained at various levels of theory.

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 pK a units by a change in the oxidation state of the other site.