Robert J. O'Reilly

Model for the exceptional reactivity of peroxiredoxins 2 and 3 with hydrogen peroxide: a kinetic and computational study.

Peroxiredoxins (Prx) are thiol peroxidases that exhibit exceptionally high reactivity toward peroxides, but the chemical basis for this is not well understood. We present strong experimental evidence that two highly conserved arginine residues play a vital role in this activity of human Prx2 and Prx3. Point mutation of either ArgI or ArgII (in Prx3 Arg-123 and Arg-146, which are ∼3–4 Å or ∼6–7 Å away from the active site peroxidative cysteine (Cp), respectively) in each case resulted in a 5 orders of magnitude loss in reactivity. A further 2 orders of magnitude decrease in the second-order rate constant was observed for the double arginine mutants of both isoforms, suggesting a cooperative function for these residues. Detailed ab initio theoretical calculations carried out with the high level G4 procedure suggest strong catalytic effects of H-bond-donating functional groups to the Cp sulfur and the reactive and leaving oxygens of the peroxide in a cooperative manner. Using a guanidinium cation in the calculations to mimic the functional group of arginine, we were able to locate two transition structures that indicate rate enhancements consistent with our experimentally observed rate constants. Our results provide strong evidence for a vital role of ArgI in activating the peroxide that also involves H-bonding to ArgII. This mechanism could explain the exceptional reactivity of peroxiredoxins toward H2O2 and may have wider implications for protein thiol reactivity toward peroxides.

An assessment of theoretical procedures for π-conjugation stabilisation energies in enones

We introduce a representative database of 22 α,β- to β,γ-enecarbonyl isomerisation energies (to be known as the EIE22 data-set). Accurate reaction energies are obtained at the complete basis-set limit CCSD(T) level by means of the high-level W1-F12 thermochemical protocol. The isomerisation reactions involve a migration of one double bond that breaks the conjugated π-system. The considered enecarbonyls involve a range of common functional groups (e.g., Me, NH2, OMe, F, and CN). Apart from π-conjugation effects, the chemical environments are largely conserved on the two sides of the reactions and therefore the EIE22 data-set allows us to assess the performance of a variety of density functional theory (DFT) procedures for the calculation of π-conjugation stabilisation energies in enecarbonyls. We find that, with few exceptions (M05-2X, M06-2X, BMK, and BH&HLYP), all the conventional DFT procedures attain root mean square deviations (RMSDs) between 5.0 and 11.7 kJ mol−1. The range-separated and double-hybrid DFT procedures, on the other hand, show good performance with RMSDs below the ‘chemical accuracy’ threshold. We also examine the performance of composite and standard ab initio procedures. Of these, SCS-MP2 offers the best performance-to-computational cost ratio with an RMSD of 0.8 kJ mol−1.

Can DFT and ab initio methods describe all aspects of the potential energy surface of cycloreversion reactions

Abstract We introduce a representative benchmark database of 20 cycloreversion reaction energies obtained by means of the high-level W1 thermochemical protocol. We use these benchmark values to assess the performance of a variety of contemporary DFT, double-hybrid DFT (DHDFT), standard ab initio, and compound thermochemistry methods. We show that this set of reaction energies provides an extremely challenging test for nearly all of the considered DFT and DHDFT methods. For example, about 80% of the considered functionals result in root-mean-square deviations (RMSDs) above 10 kJ mol−1. The best DFT and DHDFT procedures are ωB97X and DSD-PBEP86-D3, with RMSDs of 4.7 and 7.9 kJ mol−1, respectively. Coupled with the fact that the barrier heights for these reactions also pose a significant challenge for many DFT methods, this work shows that only a handful of functionals can quantitatively describe all aspects of the potential energy surface of this important class of reactions. In addition, this work shows that London dispersion effects are particularly large for this class of reactions. For example, empirical D3 dispersion corrections reduce the RMSDs for the DFT and DHDFT procedures by amounts ranging from 3.5 (PBE and B2K-PLYP) to 22.0 (BLYP) kJ mol−1. GRAPHICAL ABSTRACT

Ab Initio Investigation of the Fragmentation of 5,5-Diamino-Substituted 1,4,2-Oxathiazoles

The mechanism for the fragmentation of 5,5-diamino-1,4,2-oxathiazole derivatives has been studied at the CCSD(T)/6-311+G(3df,2p)//MP2/6-31+G(2df,p) level of theory. The calculations suggest that the fragmentation occurs via a stepwise process involving the formation of polar intermediates that lie in shallow potential wells. We find a large thermodynamic driving force for fragmentation, which together with a weakening of the C-S bond through electron donation by the amino substituents provides the impetus for a low-barrier fragmentation.