Theoretical study of the electronic spectra of small molecules that incorporate analogues of the copper-cysteine bond

The copper-sulphur bond which binds cysteinate to the metal centre is a key factor in the spectroscopy of blue copper proteins. We present theoretical calculations describing the electronically excited states of small molecules, including CuSH, CuSCH_3, (CH_3)_2SCuSH, (imidazole)-CuSH and (imidazole)_2-CuSH, derived from the active site of blue copper proteins that contain the copper-sulphur bond in order to identify small molecular systems that have electronic structure that is analogous to the active site of the proteins. Both neutral and cationic forms are studied, since these represent the reduced and oxidised forms of the protein, respectively. For CuSH and CuSH^+, excitation energies from time-dependent density functional theory with the B97-1 exchange-correlation functional agree well with the available experimental data and multireference configuration interaction calculations. For the positive ions, the singly occupied molecular orbital is formed from an antibonding combination of a 3d orbital on copper and a 3pπ orbital on sulphur, which is analogous to the protein. This leads several of the molecules to have qualitatively similar electronic spectra to the proteins. For the neutral molecules, changes in the nature of the low lying virtual orbitals leads the predicted electronic spectra to vary substantially between the different molecules. In particular, addition of a ligand bonded directly to copper results in the low-lying excited states observed in CuSH and CuSCH_33 to be absent or shifted to higher energies

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

The turbulent burning velocity of methanol-air mixtures

Methanol is a sustainable and versatile alternative fuel for spark-ignition engines and other combustion applications. To characterize the combustion behavior of this fuel, a good understanding of the factors affecting its turbulent burning velocity is required. This paper presents experimental values of the turbulent burning velocity of methanol-air mixtures obtained in a fan-stirred bomb, for u'= 2-6 m/s, phi= 0.8-1.4, T=358 K and pressures up to 0.5 MPa. In combination with laminar burning velocity values previously obtained on the same rig, these measurements are used to provide better insight into the various factors affecting u_t of methanol, and to assess to what degree existing turbulent combustion models can reproduce experimental trends. It appeared that most models correctly accounted for the effects of turbulent rms velocity u'. With respect to the effects of phi and pressure, however, models accounting for flame stretch and instabilities, through the inclusion of model terms depending on thermodiffusive mixture properties and pressure, had a slight edge on simpler formulations