Metabolic-intermediate complex formation with cytochrome P450: theoretical studies in elucidating the reaction pathway for the generation of reactive nitroso intermediate

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.

Regioselective oxidation of heterocyclic aromatic hydrocarbons catalyzed by cytochrome P450: A case study of carbazole

Recently there are increasing interests in accurately evaluating the health effects of heterocyclic PAHs. However, the activation mechanism and possible metabolites of heterocyclic PAHs catalyzed by human CYP1A1 is still elusive to a great extent. Here, leveraged to high level QM/MM calculations, the corresponding activation pathways of a representative heterocyclic PAHs, carbazole, were systematically explored. The first stage is electrophilic addition or hydrogen abstraction from N-H group. Electrophilic addition was evidenced to be more feasible and regioselectivity at C3 and C4 sites were identified. Correlations between energy barriers and key structural/electrostatic parameters reveal that O-Cα distance and Fe-O-Cα angle are the main origin for the catalytic regioselectivity. Electrophilic addition was determined as the rate-determining step and the subsequent possible reactions include epoxidation, NIH shift (the hydrogen migration from the site of hydroxylation to the adjacent carbon) and proton shuttle. The corresponding products are epoxides, ketones and hydroxylated carbazoles, respectively. The main metabolites (hydroxylated carbazoles) are estimated to be more toxic than carbazole. The regioselectivity of carbazole activated by CYP1A1 is different from the environmental processes (gas and aqueous phase). Collectively, these results will inform the in-depth understanding the metabolic processes of heterocyclic PAHs and aid the accurate evaluation of their health effects.

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.

Remdesivir: Mechanism of Metabolic Conversion from Prodrug to Drug

BACKGROUND

Remdesivir (GS-5734) has emerged as a promising drug during the challenging times of COVID-19 pandemic. Being a prodrug, it undergoes several metabolic reactions before converting to its active triphosphate metabolite. It is important to establish the atomic level details and explore the energy profile of the prodrug to drug conversion process.

METHODS

In this work, Density Functional Theory (DFT) calculations were performed to explore the entire metabolic path. Further, the potential energy surface (PES) diagram for the conversion of prodrug remdesivir to its active metabolite was established. The role of catalytic triad of Hint1 phosphoramidase enzyme in P-N bond hydrolysis was also studied on a model system using combined molecular docking and quantum mechanics approach.

RESULTS

The overall energy of reaction is 11.47 kcal/mol exergonic and the reaction proceeds through many steps requiring high activation energies. In the absence of a catalyst, the P-N bond breaking step requires 41.78 kcal/mol, which is reduced to 14.26 kcal/mol in a catalytic environment.

CONCLUSION

The metabolic pathways of model system of remdesivir (MSR) were completely explored completely and potential energy surface diagrams at two levels of theory, B3LYP/6-311++G(d, p) and B3LYP/6-31+G(d), were established and compared. The results highlight the importance of an additional water molecule in the metabolic reaction. The P-N bond cleavage step of the metabolic process requires the presence of an enzymatic environment.

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.

Importance of substituents in ring opening: a DFT study on a model reaction of thiazole to thioamide

Thiazole ring is an important active molecular skeleton of drugs. Thiazole in natural products and drugs are usually harmlessly eliminated. However, hepatotoxic reactions may occur due to the biological activation of thiazole to produce reactive thioamide. A typical example is hepatotoxic sudoxicam and safety meloxicam. The only structural difference between them is a methyl group on C5 position of thiazole in meloxicam. The molecular basis for the difference remains unknown and the bioactivation mechanism of the thiazole ring is still obscure. Quantum chemical calculations were performed to elucidate the activation mechanism of the thiazole ring under P450 catalysis, and the influence of the substituents on the activation pathways of thiazole ring was also studied. The calculated results show that the activation of thiazole is closely related to the substituents on the thiazole and spin state of Cpd I. The thiazole and substituted thiazole directly open the ring when catalyzed by doublet spin state Cpd I that catalyzed by the quartet spin state Cpd I can open the ring directly or indirectly, which is related to the substituents. Thiazoles modified with electron-donating substituents mainly undergo direct ring opening, while thiazoles modified with electron-withdrawing groups or weak electron-donating groups mainly undergo indirect ring-opening process accompanied by intermediate formation. The research results laid the foundation for the design of thiazole ring drugs, and also laid a theoretical foundation for the study of reducing the toxicity of thiazole ring drugs.

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.

Time-dependent inhibition (TDI) of CYP1A2 by a CYP3A4-mediated reactive metabolite: proposal for a novel mechanism of irreversible TDI by a non-suicide substrate

1. The purpose of this study was to clarify the mechanism of DSP-1053 time-dependent inhibition (TDI) for CYP1A2. 2. DSP-1053 inhibited time- and concentration-dependently CYP1A2 activity in human liver microsomes even in a dilution assay. However, DSP-1053 was not metabolized by recombinant human CYP1A2. These findings indicate that the inhibitory effect of DSP-1053 on CYP1A2 does not follow a general mechanism-based inhibition (MBI) because it did not seem to be a suicide substrate. 3. In fact, CYP1A2 was not inhibited with DSP-1053 pre-incubation in recombinant human CYP1A2. On the other hand, CYP1A2 was potently inhibited after pre-incubation with DSP-1053 in a mixture of human recombinant CYP1A2 and CYP3A4. In addition, DSP-1053 TDI of CYP1A2 in human liver microsomes was drastically reduced not only by addition of a CYP3A4 inhibitor, but also by addition of potassium cyanide (KCN), which is a trapping agent for iminium ions. We also confirmed in this study that CYP1A2 suicide inhibition by DSP-1053 metabolites generated by CYP3A4 had only minimal role in DSP-1053 TDI of CYP1A2. 4. In conclusion, a possible mechanism for DSP-1053 TDI of CYP1A2 is that DSP-1053 iminium ion, which is generated by CYP3A4, departs from CYP3A4 without inhibiting it and covalently binds to CYP1A2.

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.

Toxicity Originating from Thiophene Containing Drugs: Exploring the Mechanism using Quantum Chemical Methods

Drug metabolism of thiophene containing substrates by cytochrome P450s (CYP450) leads to toxic side effects, for example, nephrotoxicity (suprofen, ticlopidine), hepatotoxicity (tienilic acid), thrombotic thrombocytopenic purpura (clopidogrel), and aplastic anemia (ticlopidine). The origin of toxicity in these cases has been attributed to two different CYP450 mediated metabolic reactions: S-oxidation and epoxidation. In this work, the molecular level details of the bioinorganic chemistry associated with the generation of these competitive reactions are reported. Density functional theory was utilized (i) to explore the molecular mechanism for S-oxidation and epoxidation using the radical cationic center Cpd I [(iron(IV)-oxo-heme porphine system with SH(-) as the axial ligand, to mimic CYP450s] as the model oxidant, (ii) to establish the 3D structures of the reactants, transition states, and products on both the metabolic pathways, and (iii) to examine the potential energy (PE) profile for both the pathways to determine the energetically preferred toxic metabolite formation. The energy barrier required for S-oxidation was observed to be 14.75 kcal/mol as compared to that of the epoxidation reaction (13.23 kcal/mol) on the doublet PE surface of Cpd I. The formation of the epoxide metabolite was found to be highly exothermic (-23.24 kcal/mol), as compared to S-oxidation (-8.08 kcal/mol). Hence, on a relative scale the epoxidation process was observed to be thermodynamically and kinetically more favorable. The energy profiles associated with the reactions of the S-oxide and epoxide toxic metabolites were also explored. This study helps in understanding the CYP450-catalyzed toxic reactions of drugs containing the thiophene ring at the atomic level.

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.

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

The metabolism of amine-containing drugs by cytochrome P450 enzymes (P450s) is prone to form a nitrosoalkane metabolic intermediate (MI), which subsequently coordinates to the heme iron of a P450, to produce a metabolic-intermediate complex (MIC). This type of P450 inhibition, referred to as mechanism-based inactivation (MBI), presents a serious concern in drug discovery processes. We applied density functional theory (DFT) to the reaction between N-methylhydroxylamine (NMH) and the compound I reactive species of P450, in an effort to elucidate the mechanism of the putative final step of the MI formation in the alkylamine metabolism. Our DFT calculations show that H-abstraction from the hydroxyl group of NMH is the most favorable pathway via which the nitrosoalkane intermediate is produced spontaneously. H-abstraction from the N–H bond was slightly less favorable. In contrast, N-oxidation and H-abstraction from the C–H bond of the methyl group had much higher energy barriers. Hence, if the conversion of NMH to nitrosoalkane is catalyzed by a P450, the reaction should proceed preferentially via H-abstraction, either from the O–H bond or from the N–H bond. Our theoretical analysis of the interaction between the MI and pentacoordinate heme moieties provided further insights into the coordination bond in the MIC.

How is a metabolic intermediate formed in the mechanism-based inactivation of cytochrome P450 by using 1,1-dimethylhydrazine: hydrogen abstraction or nitrogen oxidation?

A precise understanding of the mechanism-based inactivation of cytochrome P450 enzymes (P450s) at the quantum mechanical level should allow more reliable predictions of drug-drug interactions than those currently available. Hydrazines are among the molecules that act as mechanism-based inactivators to terminate the function of P450s, which are essential heme enzymes responsible for drug metabolism in the human body. Despite its importance, the mechanism explaining how a metabolic intermediate (MI) is formed from hydrazine is not fully understood. We used density functional theory (DFT) calculations to compare four possible mechanisms underlying the reaction between 1,1-dimethylhydrazine (or unsymmetrical dimethylhydrazine, UDMH) and the reactive compound I (Cpd I) intermediate of P450. Our DFT calculations provided a clear view on how an aminonitrene-type MI is formed from UDMH. In the most favorable pathway, hydrogen is spontaneously abstracted from the N2 atom of UDMH by Cpd I, followed by a second hydrogen abstraction from the N2 atom by Cpd II. Nitrogen oxidation of nitrogen atoms and hydrogen abstraction from the C-H bond of the methyl group were found to be less favorable than the hydrogen abstraction from the N-H bond. We also found that the reaction of protonated UDMH with Cpd I is rather sluggish. The aminonitrene-type MI binds to the ferric heme more strongly than a water molecule. This is consistent with the notion that the catalytic cycle of P450 is impeded when such an MI is produced through the P450-catalyzed reaction.

Density functional study on the cytochrome-mediated S-oxidation: identification of crucial reactive intermediate on the metabolic path of thiazolidinediones

S-Oxidation is an important cytochrome P450 (CYP450)-catalyzed reaction, and the structural and energetic details of this process can only be studied by using quantum chemical methods. Thiazolidinedione (TZD) ring metabolism involving initial S-oxidation leads to the generation of reactive metabolites (RMs) and subsequent toxicity forcing the withdrawal of the glitazone class of drugs, thus, the study of the biochemical pathway of TZD ring metabolism is a subject of interest. The S-oxidation of the TZD ring and the formation of the isocyanate intermediate (ISC) was implicated as a possible pathway; however, there are several questions still unanswered in this biochemical pathway. The current study focuses on the CYP450-mediated S-oxidation, fate of the sulfoxide product (TZDSO), ring cleavage to ISC, and formation of nucleophilic adducts. The process of S-oxidation was explored by using Cpd I (iron(IV)-oxo porphyrin, to mimic CYP450) at TZVP/6-311+G(d) basis set. The barriers were calculated after incorporating dispersion and solvent corrections. The metabolic conversion from TZDSO to ISC (studied at B3LYP/6-311++G(2df,3pd)//B3LYP/6-31+G(d)) required a novel protonated intermediate, TZDSOH(+). The effect of higher basis sets (6-311+G(d,p), aug-cc-pvqz) on this conversion was studied. TZDSOH(+) was observed to be more reactive and thermodynamically accessible than ISC, indicating that TZDSOH(+) is the actual reactive intermediate leading to toxicity of the TZD class of compounds.