Importance of H-abstraction in the final step of nitrosoalkane formation in the mechanism-based inactivation of cytochrome P450 by amine-containing drugs

Mechanism of the Oxidative Ring-Closure Reaction during Gliotoxin Biosynthesis by Cytochrome P450 GliF

During gliotoxin biosynthesis in fungi, the cytochrome P450 GliF enzyme catalyzes an unusual C-N ring-closure step while also an aromatic ring is hydroxylated in the same reaction cycle, which may have relevance to drug synthesis reactions in biotechnology. However, as the details of the reaction mechanism are still controversial, no applications have been developed yet. To resolve the mechanism of gliotoxin biosynthesis and gain insight into the steps leading to ring-closure, we ran a combination of molecular dynamics and density functional theory calculations on the structure and reactivity of P450 GliF and tested a range of possible reaction mechanisms, pathways and models. The calculations show that, rather than hydrogen atom transfer from the substrate to Compound I, an initial proton transfer transition state is followed by a fast electron transfer en route to the radical intermediate, and hence a non-synchronous hydrogen atom abstraction takes place. The radical intermediate then reacts by OH rebound to the aromatic ring to form a biradical in the substrate that, through ring-closure between the radical centers, gives gliotoxin products. Interestingly, the structure and energetics of the reaction mechanisms appear little affected by the addition of polar groups to the model and hence we predict that the reaction can be catalyzed by other P450 isozymes that also bind the same substrate. Alternative pathways, such as a pathway starting with an electrophilic attack on the arene to form an epoxide, are high in energy and are ruled out.

Energy Decomposition Analysis of the Nature of Coordination Bonding at the Heme Iron Center in Cytochrome P450 Inhibition

Drug compounds or their metabolic intermediates (MIs) sometimes inhibit the function of cytochrome P450 enzymes (P450s) by forming a coordination bond with the Fe(III) heme or Fe(II) heme of P450s. Such inhibition is one of the major causes of drug-drug interactions (DDIs), a subject of longstanding academic and practical interest. However, such coordination bonding is not fully understood at the quantum mechanical level, thus hampering rational improvement of the accuracy of DDI-related predictions. In this work, we employed density functional theory (DFT) and the generalized Kohn-Sham energy decomposition analysis (GKS-EDA) scheme to investigate the nature of the coordination bonding formed in the reversible and quasi-irreversible inhibition of P450s. The GKS-EDA results highlighted a previously unrecognized role of the electron correlation effect in P450 inhibition. The correlation effect tends to be larger in Fe(II) complexes of MI-type inhibitors and is particularly prominent for the nitrosoalkane ligand. An additional natural bond orbital (NBO) analysis provided insight into the relative significance of the σ donation and π backdonation effects in various heme-inhibitor complexes.

Theoretical Investigation on the Elusive Reaction Mechanism of Spirooxindole Formation Mediated by Cytochrome P450s: A Nascent Feasible Charge-Shift C-O Bond Makes a Difference

Spirooxindoles are pivotal biofunctional groups widely distributed in natural products and clinic drugs. However, construction of such subtle chiral skeletons is a long-standing challenge to both organic and bioengineering scientists. The knowledge of enzymatic spirooxindole formation in nature may inspire rational design of new catalysts. To this end, we presented a theoretical investigation on the elusive mechanism of the spiro-ring formation at the 3-position of oxindole mediated by cytochrome P450 enzymes (P450). Our calculated results demonstrated that the electrophilic attack of CpdI, the active species of P450, to the substrate, shows regioselectivity, i.e., the attack at the C9 position forms a tetrahedral intermediate involving an unusual feasible charge-shift C9^δ+-O^δ- bond, while the attack at the C1 position forms an epoxide intermediate. The predominant route is the first route with the charge-shift bonding intermediate due to holding a relatively lower barrier by >5 kcal mol^-1 than the epoxide route, which fits the experimental observations. Such a delocalized charge-shift bond facilitates the formation of a spiro-ring mainly through elongation of the C1-C9 bond to eliminate the aromatization of the tricyclic beta-carboline. Our theoretical results shed profound mechanistic insights for the first time into the elusive spirooxindole formation mediated by P450s.

Reactive Metabolites from Thiazole-Containing Drugs: Quantum Chemical Insights into Biotransformation and Toxicity

Drugs containing thiazole and aminothiazole groups are known to generate reactive metabolites (RMs) catalyzed by cytochrome P450s (CYPs). These RMs can covalently modify essential cellular macromolecules and lead to toxicity and induce idiosyncratic adverse drug reactions. Molecular docking and quantum chemical hybrid DFT study were carried out to explore the molecular mechanisms involved in the biotransformation of thiazole (TZ) and aminothiazole (ATZ) groups leading to RM epoxide, S-oxide, N-oxide, and oxaziridine. The energy barrier required for the epoxidation is 13.63 kcal/mol, that is lower than that of S-oxidation, N-oxidation, and oxaziridine formation (14.56, 17.90, and 20.20, kcal/mol respectively). The presence of the amino group in ATZ further facilitates all the metabolic pathways, for example, the barrier for the epoxidation reaction is reduced by ∼2.5 kcal/mol. Some of the RMs/their isomers are highly electrophilic and tend to form covalent bonds with nucleophilic amino acids, finally leading to the formation of metabolic intermediate complexes (MICs). The energy profiles of these competitive pathways have also been explored.

Mechanism-Based Inactivation of Cytochrome P450 Enzymes: Computational Insights

Mechanism-based inactivation (MBI) refers to the metabolic bioactivation of a xenobiotic by cytochrome P450s to a highly reactive intermediate which subsequently binds to the enzyme and leads to the quasi-irreversible or irreversible inhibition. Xenobiotics, mainly drugs with specific functional units, are the major sources of MBI. Two possible consequences of MBI by medicinal compounds are drug-drug interaction and severe toxicity that are observed and highlighted by clinical experiments. Today almost all of these latent functional groups (e.g., thiophene, furan, alkylamines, etc.) are known, and their features and mechanisms of action, owing to the vast experimental and theoretical studies, are determined. In the past decade, molecular modeling techniques, mostly density functional theory, have revealed the most feasible mechanism that a drug undergoes by P450 enzymes to generate a highly reactive intermediate. In this review, we provide a comprehensive and detailed picture of computational advances toward the elucidation of the activation mechanisms of various known groups with MBI activity. To this aim, we briefly describe the computational concepts to carry out and analyze the mechanistic investigations, and then, we summarize the studies on compounds with known inhibition activity including thiophene, furan, alkylamines, terminal acetylene, etc. This study can be reference literature for both theoretical and experimental (bio)chemists in several different fields including rational drug design, the process of toxicity prevention, and the discovery of novel inhibitors and catalysts.

Coupled electron and proton transfer in the piperidine drug metabolism pathway by the active species of cytochromes P450

KS-DFT and MSDFT studies reveal a novel CEPT step that triggers ring contraction of piperidines by P450.

Theoretical study on the mechanism of N- and α-carbon oxidation of lapatinib catalyzed by cytochrome P450 monooxygenase

Lapatinib, an orally active dual tyrosine kinase inhibitor, is efficacious in combination therapy with capecitabine for advanced metastatic breast cancer. Despite its importance, it has been associated with hepatotoxicity observed in clinical trials and postmarketing surveillance. The mechanisms of hepatotoxicity at the chemical and cellular levels may link to drug metabolism. In this study, the N- and α-carbon oxidation processes of lapatinib catalyzed by CYP3A4 were explored by density functional theory method. The calculation results show that oxidation of C₆ is the primary metabolic process and carboxylic acid is the main metabolic product. Both hydroxylation of C₈ and subsequent formation of primary amines are feasible. However, it is not easy for the primary amines to form active metabolites nitroso, which indicates that there are other paths for the production of nitroso. Carboxylic acid is not the main metabolite of N₇ oxidation because of higher hydrolysis energy barrier of intermediate nitrone. It is worthy to study subsequent N-hydroxylation and its downstream reaction, which may be the main pathway for the formation of nitroso. These results lay the foundation for drug design and optimization.

Suicide Inhibition of Cytochrome P450 Enzymes by Cyclopropylamines via a Ring-Opening Mechanism: Proton-Coupled Electron Transfer Makes a Difference

N-benzyl-N-cyclopropylamine (BCA) has been attracting great interests for decades for its partial suicide inactivation role to cytochrome P450 (P450) via a ring-opening mechanism besides acting as a role of normal substrates. Understanding the mechanism of such partial inactivation is vital to the clinical drug design. Thus, density functional theoretical (DFT) calculations were carried out on such P450-catalyzed reactions, not only on the metabolic pathway, but on the ring-opening inactivation one. Our theoretical results demonstrated that, in the metabolic pathway, besides the normal carbinolamine, an unexpected enamine was formed via the dual hydrogen abstraction (DHA) process, in which the competition between rotation of the H-abstracted substrate radical and the rotation of hydroxyl group of the protonated Cpd II moiety plays a significant role in product branch; In the inactivation pathway, the well-noted single electron transfer (SET) mechanism-involved process was invalidated for its high energy barrier, a proton-coupled electron transfer [PCET(ET)] mechanism plays a role. Our results are consistent with other related theoretical works on heteroatom-hydrogen (X-H, X = O, N) activation and revealed new features. The revealed mechanisms will play a positive role in relative drug design.

Detoxification of 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP) by cytochrome P450 enzymes: A theoretical investigation

Two types of detoxification routes, N-demethylation to form 4-phenyl-1,2,3,6-tetrahydropyridine (PTP) and aromatic hydroxylation to generate 4-(4'-hydroxyphenyl)-1-methyl-1,2,3,6-tetrahydropyridine (MPTP-OH), for 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP) mediated by Compound I (Cpd I) of cytochrome P450 are investigated theoretically using hybrid density functional calculations. Quantum chemical results reveal that for the N-demethylation, the initial C-H bond activation is achieved via a hydrogen atom transfer (HAT) mechanism. This is followed by a subsequent O-rebound to yield the carbinolamine intermediate. Due to the nature of pericyclic reaction, the generated carbinolamine decomposes in a non-enzymatic aqueous environment with the assistance of water molecules, forming amine and hydrated formaldehyde. For the aromatic hydroxylation, an initial addition of Cpd I to the substrate occurs mainly through a side-on approach with a subsequent proton shuttle to form the phenol product. A comparison of the energy barriers for both routes indicates that the N-demethylation (7.5/5.7kcal/mol for the quartet/doublet state in solvent) is thermodynamically more favorable than the aromatic hydroxylation process (14.9/14.8kcal/mol for the quartet/doublet state in solvent). This trend is in good agreement with the experimental product distribution, viz., the N-demethylation product PTP is more than the aromatic hydroxylation product MPTP-OH. Taken together, these observations not only enrich our knowledge on the mechanistic details of the N-dealkylation and the aromatic hydroxylation by P450s, but also provide certain insights into the metabolism of other analogous toxins.

Applications of density functional theory to iron-containing molecules of bioinorganic interest

The past decades have seen an explosive growth in the application of density functional theory (DFT) methods to molecular systems that are of interest in a variety of scientific fields. Owing to its balanced accuracy and efficiency, DFT plays particularly useful roles in the theoretical investigation of large molecules. Even for biological molecules such as proteins, DFT finds application in the form of, e.g., hybrid quantum mechanics and molecular mechanics (QM/MM), in which DFT may be used as a QM method to describe a higher prioritized region in the system, while a MM force field may be used to describe remaining atoms. Iron-containing molecules are particularly important targets of DFT calculations. From the viewpoint of chemistry, this is mainly because iron is abundant on earth, iron plays powerful (and often enigmatic) roles in enzyme catalysis, and iron thus has the great potential for biomimetic catalysis of chemically difficult transformations. In this paper, we present a brief overview of several recent applications of DFT to iron-containing non-heme synthetic complexes, heme-type cytochrome P450 enzymes, and non-heme iron enzymes, all of which are of particular interest in the field of bioinorganic chemistry. Emphasis will be placed on our own work.

Theoretical insights into the reductive metabolism of CCl4 by cytochrome P450 enzymes and the CCl4-dependent suicidal inactivation of P450

The one-electron reduction product, ˙CCl₃, irreversibly inactivates P450 via covalently binding to the meso-carbon, whereas the two successive one-electron reductions product, :CCl₂, reversibly inhibits P450 by coordinating to iron.