BDE261: a comprehensive set of high-level theoretical bond dissociation enthalpies

We have used the high-level W1w protocol to compile a comprehensive collection of 261 bond dissociation enthalpies (BDEs) for bonds connecting hydrogen, first-row and second-row p-block elements. Together they cover 45 bond types, and we term this the BDE261 set. We have used these benchmark values to assess the performance of computationally less demanding theoretical procedures, including density functional theory (DFT), double-hybrid DFT (DHDFT), and high-level composite procedures. We find that the M06-2X (DFT), ROB2-PLYP and DuT-D3 (DHDFT), and G3X(MP2)-RAD and G4(MP2)-6X (composite) procedures yield absolute BDEs with satisfactory to excellent accuracy. Overall, we recommend G4(MP2)-6X as an accurate and relatively cost-effective procedure for the direct computation of BDEs. One important finding is that the deviations for DFT and (especially) DHDFT procedures are often quite systematic. This allows an alternative approach to obtaining accurate absolute BDEs, namely, to evaluate accurate relative BDEs (RBDEs) using a computationally less demanding procedure, and to use these RBDEs in combination with appropriate and accurate reference BDEs to give accurate absolute BDEs. We recommend DuT-D3 for this purpose. For a still less computationally demanding approach, we introduce the deviation from additivity of the RBDE (DARBDE), and demonstrate that the combination of lower-level DARBDEs for larger systems and higher-level (W1w) reference RBDEs and BDEs for small systems can be utilized to obtain improved RBDEs for multiply substituted systems at low cost.

Hydrogen abstraction by chlorine atom from amino acids: remarkable influence of polar effects on regioselectivity

Quantum chemistry computations have been used to investigate hydrogen-atom abstraction by chlorine atom from protonated and N-acetylated amino acids. The results are consistent with the decreased reactivity at the backbone α-carbon and adjacent side-chain positions that is observed experimentally. The individual effects of NH(3)(+), COOH, and NHAc substituents have been examined and reveal important insights. The NH(3)(+) group in isolation is found to be deactivating at the α-position, while the acetamido group is activating. For the COOH group, polar effects lead to a contrathermodynamic deactivation of the thermodynamically most favorable α-abstraction. In the N-acetylamino acid, the α-position is deactivated by the combined inductive effect of the substituents and the presence of an early transition structure, again overriding the greater thermodynamic stability of the α-centered radical product. Deactivation of the α-, β-, and γ-positions results in a peculiar stability for amino acids and peptides and their derivatives with respect to radical degradation.

Obtaining Good Performance With Triple-ζ-Type Basis Sets in Double-Hybrid Density Functional Theory Procedures

A variety of combinations of B-LYP-based double-hybrid density functional theory (DHDFT) procedures and basis sets have been examined. A general observation is that the optimal combination of exchange contributions is in the proximity of 30% Becke 1988 (B88) exchange and 70% Hartree-Fock (HF) exchange, while for the correlation contributions, the use of independently optimized spin-component-scaled Møller-Plesset second-order perturbation theory (SCS-MP2) parameters (MP2OS and MP2SS) is beneficial. The triple-ζ Dunning aug'-cc-pVTZ+d and Pople 6-311+G(3df,2p)+d basis sets are found to be cost-effective for DHDFT methods. As a result, we have formulated the DuT-D3 DHDFT procedure, which employs the aug'-cc-pVTZ+d basis set and includes 30% B88 and 70% HF exchange energies, 59% LYP, 47% MP2OS, and 36% MP2SS correlation energies, and a D3 dispersion correction with the parameters s6 = 0.5, sr,6 = 1.569, and s8 = 0.35. Likewise, the PoT-D3 DHDFT procedure was formulated with the 6-311+G(3df,2p)+d basis set and has 32% B88 and 68% HF exchange energies, 63% LYP, 46% MP2OS, and 27% MP2SS correlation energies, and the D3 parameters s6 = 0.5, sr,6 = 1.569, and s8 = 0.30. Testing using the large E3 set of 740 energies demonstrates the robustness of these methods. Further comparisons show that the performance of these methods, particularly DuT-D3, compares favorably with the previously reported DSD-B-LYP and DSD-B-LYP-D3 methods used in conjunction with quadruple-ζ aug'-pc3+d and aug'-def2-QZVP basis sets but at lower computational expense. The previously reported ωB97X-(LP)/6-311++G(3df,3pd) procedure also performs very well. Our findings highlight the cost-effectiveness of appropriate- and moderate-sized triple-ζ basis sets in the application of DHDFT procedures.

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 (C(p)), 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 C(p) 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 H(2)O(2) and may have wider implications for protein thiol reactivity toward peroxides.

One-electron reduction of N-chlorinated and N-brominated species is a source of radicals and bromine atom formation

Hypochlorous (HOCl) and hypobromous (HOBr) acids are strong bactericidal oxidants that are generated by the human immune system but are implicated in the development of many human inflammatory diseases (e.g., atherosclerosis, asthma). These oxidants react readily with sulfur- and nitrogen-containing nucleophiles, with the latter generating N-halogenated species (e.g., chloramines/bromamines (RR'NX; X = Cl, Br)) as initial products. Redox-active metal ions and superoxide radicals (O(2)(•-)) can reduce N-halogenated species to nitrogen- and carbon-centered radicals. N-Halogenated species and O(2)(•-) are generated simultaneously at sites of inflammation, but the significance of their interactions remains unclear. In the present study, rate constants for the reduction of N-halogenated amines, amides, and imides to model potential biological substrates have been determined. Hydrated electrons reduce these species with k(2) > 10(9) M(-1) s(-1), whereas O(2)(•-) reduced only N-halogenated imides with complex kinetics indicative of chain reactions. For N-bromoimides, heterolytic cleavage of the N-Br bond yielded bromine atoms (Br(•)), whereas for other substrates, N-centered radicals and Cl(-)/Br(-) were produced. High-level quantum chemical procedures have been used to calculate gas-phase electron affinities and aqueous solution reduction potentials. The effects of substituents on the electron affinities of aminyl, amidyl, and imidyl radicals are rationalized on the basis of differential effects on the stabilities of the radicals and anions. The calculated reduction potentials are consistent with the experimental observations, with Br(•) production predicted for N-bromosuccinimide, while halide ion formation is predicted in all other cases. These data suggest that interaction of N-halogenated species with O(2)(•-) may produce deleterious N-centered radicals and Br(•).

O - H Bond Dissociation Energies

High-level composite, ab initio and density functional theory (DFT) procedures have been employed to study O–H bond dissociation energies (BDEs), as well as radical stabilization energies (RSEs) in the oxygen-centred radicals that are formed in the dissociation of the O–H bonds. Benchmark values are provided by Wn results up to W3.2 and W4.x. We are able to recommend revised BDE values for FO–H (415.6 ± 3 kJ mol–1), MeC(O)O–H (459.8 ± 6 kJ mol–1) and CF3CH2O–H (461.9 ± 6 kJ mol–1) on the basis of high-level calculations. We find that Gn-type procedures are generally reliable and cost-effective, and that some contemporary functionals and double-hybrid DFT procedures also provide adequate O–H BDEs/RSEs. We note that the variations in the O–H BDEs are associated with variations in the stabilities of not only the radicals but also the closed-shell precursor molecules. Most substituents destabilize both species, with σ-electron-withdrawing groups having larger destabilizing effects, while π-electron acceptors are stabilizing. Although there is little correlation between the stabilizing/destabilizing effects of the substituents and the RSEs, we present some general patterns in the RSEs that emerge from the present study.

A computational study of methanol-to-hydrocarbon conversion — Towards the design of a low-barrier process

Computational quantum chemistry has been employed to examine the production of ethylene with methanol-to-hydrocarbon (MTH) processes via a carbon pool mechanism. We find that the M05-2X functional performs well for the types of reactions that are involved. The methylation reactions of the aromatic cocatalyst are the most energy-demanding steps in the process. For the subsequent production of C2H4, we have identified a low-energy pathway that involves multiple methyl shifts, followed by concerted deprotonation and C2H4 elimination. The substitutions of the Al and Si atoms in the participating Si–OH–Al moiety of zeolite catalysts with Ga and Ge do not lead to lower barriers for the methylation reactions, nor does the use of a more electron-rich aromatic cocatalyst. However, we find that the use of two cocatalysts, a nucleophile and an aromatic carbon pool, can provide an overall low-energy pathway for the MTH process.