Dioxygen Binding in the Active Site of Histone Demethylase JMJD2A and the Role of the Protein Environment

JMJD2A catalyses the demethylation of di‐ and trimethylated lysine residues in histone tails and is a target for the development of new anticancer medicines. Mechanistic details of demethylation are yet to be elucidated and are important for the understanding of epigenetic processes. We have evaluated the initial step of histone demethylation by JMJD2A and demonstrate the dramatic effect of the protein environment upon oxygen binding using quantum mechanics/molecular mechanics (QM/MM) calculations. The changes in electronic structure have been studied for possible spin states and different conformations of O₂, using a combination of quantum and classical simulations. O₂ binding to this histone demethylase is computed to occur preferentially as an end‐on superoxo radical bound to a high‐spin ferric centre, yielding an overall quintet ground state. The favourability of binding is strongly influenced by the surrounding protein: we have quantified this effect using an energy decomposition scheme into electrostatic and dispersion contributions. His182 and the methylated lysine assist while Glu184 and the oxoglutarate cofactor are deleterious for O₂ binding. Charge separation in the superoxo‐intermediate benefits from the electrostatic stabilization provided by the surrounding residues, stabilizing the binding process significantly. This work demonstrates the importance of the extended protein environment in oxygen binding, and the role of energy decomposition in understanding the physical origin of binding/recognition.     

Data-mining the diaryl(thio)urea conformational landscape: Understanding the contrasting behavior of ureas and thioureas with quantum chemistry

The conformations adopted by urea and thiourea functional groups influence catalysis and binding. We combine data-mining with quantum chemical calculations to understand the differences in conformational behavior for these two important structural motifs. We developed a Python tool to automate the compilation of X-ray structural information and perform conformational clustering and visualization, based on SMILES input. While diarylureas have an overwhelming preference for the anti,anti-conformer, diarylthioureas adopt a mixture of anti,anti- and anti,syn-conformers. Computations show the anti,anti-thiourea conformer is destabilized by out-of-plane rotations which avoid a steric clash with the sulfur atom. These conformational preferences were studied computationally under a variety of conditions, and apart from in the gas-phase, a preference for anti,anti-ureas was found. Consistent with experiments, this preference increases in more polar environments. Quantitative predicted ratios are sensitive to the computational treatment of solvation effects, with COSMO-RS giving more realistic amounts of the anti,anti-conformer in THF and DMSO.     

Stereostructure assignment of flexible five-membered rings by GIAO 13C NMR calculations: prediction of the stereochemistry of elatenyne

The stereochemistry of conformationally mobile five-membered rings is often hard to assign from NMR data, and [2,2']bifuranyl systems are even more challenging. GIAO (13)C NMR chemical shifts have been calculated for a series of [2,2']bifuranyl and pyranopyran species, taking into account their conformational flexibility using weighted averages of the data for all low energy conformers. We show that calculation of (13)C NMR chemical shifts using the geometries obtained using molecular mechanics greatly reduces the computational expense without a significant loss of accuracy, even in this demanding system. The results were sufficiently accurate to distinguish not only the pyran and furanyl isomers but also between all the diastereoisomeric forms. As a result of this validation, we predict the stereochemistry for the recently proposed revised structure of the natural product elatenyne, which contains a [2,2']bifuranyl core.     

Metal-Free Arylation of Benzothiophenes at C4 by Activation as their Benzothiophene S-Oxides

Benzothiophenes, activated by oxidation to the corresponding S-oxides, undergo C−H/C−H-type coupling with phenols to give C4 arylation products. While an electron-withdrawing group at C3 of the benzothiophene is important, the process operates without a directing group and a metal catalyst, thus rendering it compatible with sensitive functionalities—e.g. halides and formyl groups. Quantum chemical calculations suggest a formal stepwise mechanism involving heterolytic cleavage of an aryloxysulfur species to give a π-complex of the corresponding benzothiophene and a phenoxonium cation. Subsequent addition of the phenoxonium cation to the C4 position of the benzothiophene is favored over the addition to C3; Fukui functions predict that the major regioisomer is formed at the more electron-rich position between C3 and C4. Varied selective manipulation of the benzothiophene products showcase the synthetic utility of the metal-free arylation process.     

Enantioselectivity in the boron aldol reactions of methyl ketones

DFT computed transition states quantitatively explain the surprising stereochemical outcome of unsubstituted enolborinates in diastereoselective and enantioselective boron aldol reactions.