Rieske dearomatizing dioxygenases utilize a Rieske iron-sulfur cluster and a mononuclear Fe(II) located 15 Å across a subunit boundary to catalyze O2-dependent formation of cis-dihydrodiol products from aromatic substrates. During catalysis, O2 binds to the Fe(II) while the substrate bind nearby. Single turnover reactions have shown that one electron from each metal center is required for catalysis. This finding suggested that the reactive intermediate is Fe(III)-(H)peroxo or HO-Fe(V)=O formed by O-O bond scission. Surprisingly, several kinetic phases were observed during the single turnover Rieske cluster oxidation. Here, the Rieske cluster oxidation and product formation steps of a single turnover of benzoate 1,2-dioxygenase are investigated using benzoate and three fluorinated analogs. It is shown that the rate constant for product formation correlates with the reciprocal relaxation time of only the fastest kinetic phase (RRT-1) for each substrate, suggesting that the slower phases are not mechanistically relevant. RRT-1 is strongly dependent on substrate type, suggesting a role for substrate in electron transfer from the Rieske cluster to the mononuclear iron site. This insight, together with the substrate and O2 concentration dependencies of RRT-1, indicates that a reactive species is formed after substrate and O2 binding, but before electron transfer from the Rieske cluster. Computational studies show that RRT-1 is correlated with the electron density at the substrate carbon closest to the Fe(II), consistent with initial electrophilic attack by an Fe(III)-superoxo intermediate. The resulting Fe(III)-peroxo-aryl radical species would then readily accept an electron from the Rieske cluster to complete the cis-dihydroxylation reaction.

Studies of the kinetics of polymerization of ε-caprolactone (CL) by salen-aluminum catalysts comprising ligands with similar steric profiles but different electron donating characteristics (R = OMe, Br, or NO2) were performed using high initial monomer concentrations (2 M < [CL]0 < 2.6 M) in toluene-d8 at temperatures ranging from 20 to 90 °C. Saturation behavior was observed, enabling determination of monomer equilibrium constants (Keq) and catalytic rate constants (k2) as a function of R and temperature. While Keq varied only slightly with the electron donating properties of R (Hammett ρ = +0.16(8)), k2 showed a more significant dependence reflected by ρ = +1.4(1). Thermodynamic parameters ΔG° (associated with Keq) and ΔG(‡) (associated with k2) were determined, with the former being ∼0 kcal/mol for all catalysts and the latter exhibiting the trend R = OMe > Br > NO2. Density functional theory (DFT) calculations were performed to characterize mechanistic pathways at a microscopic level of detail. Lowest energy transition-state structures feature incipient bonding of the nucleophile to the lactone carbonyl that is approaching the metal ion, but a distinct CL adduct is not an energy minimum on the reaction pathway, arguing against Keq being associated with coordination of monomer according to the typical coordination-insertion mechanism. An alternative hypothesis is presented associating Keq with "nonproductive" coordination of substrate in a manner that inhibits the polymerization reaction at high substrate concentrations.

Aluminum alkoxide complexes (2) of salen ligands with a three-carbon linker and para substituents having variable electron-withdrawing capabilities (X = NO2, Br, OMe) were prepared, and the kinetics of their ring-opening polymerization (ROP) of ε-caprolactone (CL) were investigated as a function of temperature, with the aim of drawing comparisons to similar systems with two-carbon linkers investigated previously (1). While 1 and 2 exhibit saturation kinetics and similar dependences of their ROP rates on substituents X (invariant Keq, similar Hammett ρ = +1.4(1) and 1.2(1) for k2, respectively), ROP by 2 was significantly faster than for 1. Theoretical calculations confirm that, while the reactant structures differ, the transition state geometries are quite similar, and by analyzing the energetics of the involved distortions accompanying the structural changes, a significant contribution to the basis for the rate differences was identified. Using this knowledge, a simplified computational method for evaluating ligand structural influences on cyclic ester ROP rates is proposed that may have utility for future catalyst design.

We report an NMR chemical shift study of conformationally challenging seven-membered lactones (1–11); computed and experimental data sets are compared. The computations involved full conformational analysis of each lactone, Boltzmann-weighted averaging of the chemical shifts across all conformers, and linear correction of the computed chemical shifts. DFT geometry optimizations [M06-2X/6-31+G(d,p)] and GIAO NMR chemical shift calculations [B3LYP/6-311+G(2d,p)] provide the computed chemical shifts. The corrected-mean absolute error (CMAE), the average of the differences between the computed and experimental chemical shifts for each of the eleven lactones, is encouragingly small (0.02–0.08 ppm for 1H or 0.8–2.2 ppm for 13C). Three pairs of cis vs. trans diastereomeric lactones were used to assess the ability of the method to distinguish between stereoisomers. The experimental shifts were compared with the computed shifts for each of the two possible isomers. We introduce the use of a “match ratio”—the ratio of the larger (worse fit) to the smaller (better fit) CMAE. A greater match ratio value indicates better distinguishing ability. The match ratios are larger for proton data [2.4–4.0 (ave = 3.2)] than for carbon [1.1–2.3 (ave = 1.6)], indicating that the former provide a better basis for discriminating these diastereomers.

Theoretical analysis of the mechanism of the intramolecular hexadehydro-Diels–Alder (HDDA) reaction, validated against prior and newly measured kinetic data for a number of different tethered yne-diynes, indicates that the reaction proceeds in a highly asynchronous fashion. The rate-determining step is bond formation at the alkyne termini nearest the tether, which involves a transition-state structure exhibiting substantial diradical character. Whether the reaction then continues to close the remaining bond in a concerted fashion or in a stepwise fashion (i.e., with an intervening intermediate) depends on the substituents at the remaining terminal alkyne positions. Computational modeling of the HDDA reaction is complicated by the significant diradical character that arises along the reaction coordinate, which leads to instabilities in both restricted singlet Kohn-Sham density functional theory (DFT) and coupled-cluster theory based on a Hartree-Fock reference wave function. A consistent picture emerges, however, from comparison of broken-symmetry DFT calculations and second-order perturbation theory based on complete-active-space self-consistent-field (CASPT2) calculations.