The mechanism of human aromatase (CYP 19A1) revisited: DFT and QM/MM calculations support a compound I-mediated pathway for the aromatization process

Aromatase inhibitors for the treatment of breast cancer: An overview (2019-2023)

Aromatase inhibition is considered a legitimate approach for the treatment of ER-positive (ER+) breast cancer as it accounts for more than 70% of breast cancer cases. Aromatase inhibitor therapy has been demonstrated to be highly effective in decreasing tumour size, increasing survival rates, and lowering the chance of cancer recurrence. The present review deliberates the pathophysiology and the role of aromatase in estrogen biosynthesis. Estrogen biosynthesis, various androgens, and their function in the human body have also been discussed. The salient aspects of the aromatase active site, its mode of action, and AIs, along with their intended interactions with presently FDA-approved inhibitors, have been briefly discussed. It has been detailed how different reported AIs were designed, their SAR investigations, in silico analysis, and biological evaluations. Various AIs from multiple origins, such as synthetic and semi-synthetic, have also been discussed.

Roles of Ferric Peroxide Anion Intermediates (Fe3+O2 -, Compound 0) in Cytochrome P450 19A1 Steroid Aromatization and a Cytochrome P450 2B4 Secosteroid Oxidation Model

Cytochrome P450 (P450, CYP) 19A1 is the steroid aromatase, the enzyme responsible for the 3-step conversion of androgens (androstenedione or testosterone) to estrogens. The final step is C-C bond scission (removing the 19-oxo group as formic acid) that proceeds via a historically controversial reaction mechanism. The two competing mechanistic possibilities involve a ferric peroxide anion (Fe3+O2 -, Compound 0) and a perferryl oxy species (FeO3+, Compound I). One approach to discern the role of each species in the reaction is with the use of oxygen-18 labeling, i.e., from 18O2 and H2 18O of the reaction product formic acid. We applied this approach, using several technical improvements, to study the deformylation of 19-oxo-androstenedione by human P450 19A1 and of a model secosteroid, 3-oxodecaline-4-ene-10-carboxaldehyde (ODEC), by rabbit P450 2B4. Both aldehyde substrates were sensitive to non-enzymatic acid-catalyzed deformylation, yielding 19-norsteroids, and conditions were established to avoid issues with artifactual generation of formic acid. The Compound 0 reaction pathway predominated (i.e., Fe3+O2 -) in both P450 19A1 oxidation of 19-oxo-androstenedione and P450 2B4 oxidation of ODEC. The P450 19A1 results contrast with our prior conclusions (J. Am. Chem. Soc. 2014, 136, 15016-16025), attributed to several technical modifications.

Theoretical Study of Hydroxylation of α- and β-Pinene by a Cytochrome P450 Monooxygenase Model

Previous studies on biocatalytic transformations of pinenes by cytochrome P450 (CYP) enzymes reveal the formation of different oxygenated products from a single substrate due to the multistate reactivity of CYP and the many reactive sites in the pinene scaffold. Up until now, the detailed mechanism of these biocatalytic transformations of pinenes have not been reported. Hereby, we report a systematic theoretical study of the plausible hydrogen abstraction and hydroxylation reactions of α- and β-pinenes by CYP using the density functional theory (DFT) method. All DFT calculations in this study were based on B3LYP/LAN computational methodology using the Gaussian09 software. We used the B3LYP functional with corrections for dispersive forces, BSSE, and anharmonicity to study the mechanism and thermodynamic properties of these reactions using a bare model (without CYP) and a pinene-CYP model. According to the potential energy surface and Boltzmann distribution for radical conformers, the major reaction products of CYP-catalyzed hydrogen abstraction from β-pinene are the doublet trans (53.4%) and doublet cis (46.1%) radical conformer at delta site. The formation of doublet cis/trans hydroxylated products released a total Gibbs free energy of about 48 kcal/mol. As for alpha pinene, the most stable radicals were trans-doublet (86.4%) and cis-doublet (13.6%) at epsilon sites, and their hydroxylation products released a total of ~50 kcal/mol Gibbs free energy. Our results highlight the likely C-H abstraction and oxygen rebounding sites accounting for the multi-state of CYP (doublet, quartet, and sextet spin states) and the formation of different conformers due to the presence of cis/trans allylic hydrogen in α-pinene and β-pinene molecules.

A comprehensive insight into aldehyde deformylation: mechanistic implications from biology and chemistry

Aldehyde deformylation is an important reaction in biology, organic chemistry and inorganic chemistry and the process has been widely applied and utilized. For instance, in biology, the aldehyde deformylation reaction has wide differences in biological function, whereby cyanobacteria convert aldehydes into alkanes or alkenes, which are used as natural products for, e.g., defense mechanisms. By contrast, the cytochromes P450 catalyse the biosynthesis of hormones, such as estrogen, through an aldehyde deformylation reaction step. In organic chemistry, the aldehyde deformylation reaction is a common process for replacing functional groups on a molecule, and as such, many different synthetic methods and procedures have been reported that involve an aldehyde deformylation step. In bioinorganic chemistry, a variety of metal(iii)-peroxo complexes have been synthesized as biomimetic models and shown to react efficiently with aldehydes through deformylation reactions. This review paper provides an overview of the various aldehyde deformylation reactions in organic chemistry, biology and biomimetic model systems, and shows a broad range of different chemical reaction mechanisms for this process. Although a nucleophilic attack at the carbonyl centre is the consensus reaction mechanism, several examples of an alternative electrophilic reaction mechanism starting with hydrogen atom abstraction have been reported as well. There is still much to learn and to discover on aldehyde deformylation reactions, as deciphered in this review paper.

Quantum-Mechanical/Molecular-Mechanical Studies of CYP11A1-Catalyzed Biosynthesis of Pregnenolone from Cholesterol Reveal a C-C Bond Cleavage Reaction That Occurs by a Compound I-Mediated Electron Transfer

We explore here a long-standing mechanistic question by using quantum-mechanical/molecular-mechanical (QM/MM) methodology. The question concerns the mechanism of steroid hormone biosynthesis, whereby the P450 enzyme, CYP11A1, catalyzes the C20-C22 bond-cleavage in the 20,22-hydroxylated cholesterol, 20R,22R-DiOHCH, leading to pregnenolone, which is critical for the subsequent production of all steroid hormones. This is an unusual feat whereby the P450 enzyme breaks two O-H bonds and one C-C bond, while making two C═O bonds. How does the enzyme perform such a complex and highly energy-demanding reaction? Our computational results rule out the previously proposed Compound I (Cpd I) electrophilic attack mechanism via the formation of a peroxide intermediate as well as the H-abstraction-mediated C-C cleavage mechanism. Notably, oxygen-rebound cannot transpire, in spite of the fact that the classical active species, Cpd I, participates in the catalytic process. Our findings reveal a mechanism whereby C-C bond cleavage is mediated by an electron transfer from the C22-O^--deprotonated substrate to Cpd I. As such, our QM/MM calculations demonstrate that Cpd I acts as an electron sink that facilitates the C-C bond cleavage.

Density Functional Theory Study on the Demethylation Reaction between Methylamine, Dimethylamine, Trimethylamine, and Tamoxifen Catalyzed by a Fe(IV)-Oxo Porphyrin Complex

In this work, we studied computationally the N-demethylation reaction of methylamine, dimethylamine, and trimethylamine as archetypal examples of primary, secondary, and tertiary amines catalyzed by high-field low-spin Fe-containing enzymes such as cytochromes P450. Using DFT calculations, we found that the expected C-H hydroxylation process was achieved for trimethylamine. When dimethylamine and methylamine were studied, two different reaction mechanisms (C-H hydroxylation and a double hydrogen atom transfer) were computed to be energetically accessible and both are equally preferred. Both processes led to the formation of formaldehyde and the N-demethylated substrate. Finally, as an illustrative example, the relative contribution of the three primary oxidation routes of tamoxifen was rationalized through energetic barriers obtained from density functional calculations and docking experiments involving CYP3A4 and CYP2D6 isoforms. We found that the N-demethylation process was the intrinsically favored one, whereas other oxidation reactions required most likely preorganization imposed by the residues close to the active sites.

Combined Docking and Quantum Chemical Study on CYP-Mediated Metabolism of Estrogens in Man

Long-term exposure to estrogens seriously increases the incidence of various diseases including breast cancer. Experimental studies indicate that cytochrome P450 (CYP) enzymes catalyze the bioactivation of estrogens to catechols, which can exert their harmful effects via various routes. It has been shown that the 4-hydroxylation pathway of estrogens is the most malign, while 2-hydroxylation is considered a benign pathway. It is also known experimentally that with increasing unsaturation of ring B of estrogens the prevalence of the 4-hydroxylation pathway significantly increases. In this study, we used a combination of structural analysis, docking, and quantum chemical calculations at the B3LYP/6-311+G* level to investigate the factors that influence the regioselectivity of estrogen metabolism in man. We studied the structure of human estrogen metabolizing enzymes (CYP1A1, CYP1A2, CYP1B1, and CYP3A4) in complex with estrone using docking and investigated the susceptibility of estrone, equilin, and equilenin (which only differ in the unsaturation of ring B) to undergo 2- and 4-hydroxylation using several models of CYP enzymes (Compound I, methoxy, and phenoxy radical). We found that even the simplest models could account for the experimental difference between the 2- and 4- hydroxylation pathways and thus might be used for fast screening purposes. We also show that reactivity indices, specifically in this case the radical and nucleophilic condensed Fukui functions, also correctly predict the likeliness of estrogen derivatives to undergo 2- or 4-hydroxylation.

Theoretical Study of the Mechanism of Exemestane Hydroxylation Catalyzed by Human Aromatase Enzyme

Human aromatase (CYP19A1) aromatizes the androgens to form estrogens via a three-step oxidative process. The estrogens are necessary in humans, mainly in women, because of the role they play in sexual and reproductive development. However, these also are involved in the development and growth of hormone-dependent breast cancer. Therefore, inhibition of the enzyme aromatase, by means of drugs known as aromatase inhibitors, is the frontline therapy for these types of cancers. Exemestane is a suicidal third-generation inhibitor of aromatase, currently used in breast cancer treatment. In this study, the hydroxylation of exemestane catalyzed by aromatase has been studied by means of hybrid QM/MM methods. The Free Energy Perturbation calculations provided a free energy of activation for the hydrogen abstraction step (rate-limiting step) of 17 kcal/mol. The results reveal that the hydroxylation of exemestane is not the inhibition stage, suggesting a possible competitive mechanism between the inhibitor and the natural substrate androstenedione in the first catalytic subcycle of the enzyme. Furthermore, the analysis of the interaction energy for the substrate and the cofactor in the active site shows that the role of the enzymatic environment during this reaction consists of a transition state stabilization by means of electrostatic effects.

Recent application of calculations of metal complexes based on density functional theory

Recent application of density functional theory (DFT) for metal complexes is reviewed to show the achievements of DFT and the challenges for it, as well as the methods for selecting proper functionals.