Ab initio dynamics of the cytochrome P450 hydroxylation reaction

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.

Dynamical Origin of Rebound versus Dissociation Selectivity during Fe-Oxo-Mediated C-H Functionalization Reactions

The mechanism of catalytic C-H functionalization of alkanes by Fe-oxo complexes is often suggested to involve a hydrogen atom transfer (HAT) step with the formation of a radical-pair intermediate followed by diverging pathways for radical rebound, dissociation, or desaturation. Recently, we showed that in some Fe-oxo reactions, the radical pair is a nonstatistical-type intermediate and dynamic effects control rebound versus dissociation pathway selectivity. However, the effect of the solvent cage on the stability and lifetime of the radical-pair intermediate has never been analyzed. Moreover, because of the extreme complexity of motion that occurs during dynamics trajectories, the underlying physical origin of pathway selectivity has not yet been determined. For the reaction between [(TQA_Cl)FeIVO]+ and cyclohexane, here, we report explicit solvent trajectories and machine learning analysis on transition-state sampled features (e.g., vibrational, velocity, and geometric) that identified the transferring hydrogen atom kinetic energy as the most important factor controlling rebound versus nonrebound dynamics trajectories, which provides an explanation for our previously proposed dynamic matching effect in fast rebound trajectories that bypass the radical-pair intermediate. Manual control of the reaction trajectories confirmed the importance of this feature and provides a mechanism to enhance or diminish selectivity for the rebound pathway. This led to a general catalyst design principle and proof-of-principle catalyst design that showcases how to control rebound versus dissociation reaction pathway selectivity.

Molecular Dynamics of Iron Porphyrin-Catalyzed C-H Hydroxylation of Ethylbenzene

Quasi-classical molecular dynamics (MD) simulations were carried out to study the mechanism of iron porphyrin-catalyzed hydroxylation of ethylbenzene. The hydrogen atom abstraction from ethylbenzene by iron-oxo species is the rate-determining step, which generates the radical pair of iron-hydroxo species and the benzylic radical. In the subsequent radical rebound step, the iron-hydroxo species and benzylic radical recombine to form the hydroxylated product, which is barrierless on the doublet energy surface. In the gas-phase quasi-classical MD study on the doublet energy surface, 45% of the reactive trajectories lead directly to the hydroxylated product, and this increases to 56% in implicit solvent model simulations. The percentage of reactive trajectories leading to the separated radical pair is 98-100% on high-spin (quartet/sextet) energy surfaces. The low-spin state reactivity dominates in the hydroxylation of ethylbenzene, which is dynamically both concerted and stepwise, since the time gap between C-H bond cleavage and C-O bond formation ranges from 41 to 619 fs. By contrast, the high-spin state catalysis is an energetically stepwise process, which has a negligible contribution to the formation of hydroxylation products.

Direct Dynamics Trajectories Demonstrate Dynamic Matching and Nonstatistical Radical Pair Intermediates during Fe-Oxo-Mediated C-H Functionalization Reactions

The generally proposed mechanism for the reaction between non-heme Fe-oxo complexes and alkane C-H bonds involves a hydrogen atom transfer (HAT) reaction step with a radical pair intermediate that then has competitive radical rebound, dissociation, or desaturation pathways. Here, we report density functional theory-based quasiclassical direct dynamics trajectories that examine post-HAT reaction dynamics. Trajectories revealed that the radical pair intermediate can be a nonstatistical type intermediate without complete internal vibrational redistribution and post-HAT selectivity is generally determined by dynamic effects. Fast rebound trajectories occur through dynamic matching between the rotational motion of the newly formed Fe-OH bond and collision with the alkane radical, and all of this occurs through a nonsynchronous dynamically concerted process that circumvents the radical pair intermediate structure. For radical pair dissociation, trajectories proceeded to the radical pair intermediate for a very brief time, followed by complete dissociation. These trajectories provide a new viewpoint and model to understand the inherent reaction pathway selectivity for non-heme Fe-oxo-mediated C-H functionalization reactions.

Mechanism of Melatonin Metabolism by CYP1A1: What Determines the Bifurcation Pathways of Hydroxylation versus Deformylation?

Melatonin, a widely applied cosmetic active ingredient, has a variety of uses as a skin protector through antioxidant and anti-inflammatory functions as well as giving the body UV-induced defenses and immune system support. In the body, melatonin is synthesized from a tryptophan amino acid in a cascade of reactions, but as melatonin is toxic at high concentrations, it is metabolized in the human skin by the cytochrome P450 enzymes. The P450s are diverse heme-based mono-oxygenases that catalyze oxygen atom-transfer processes that trigger metabolism and detoxification reactions in the body. In the catalytic cycle of the P450s, a short-lived high-valent iron(IV)-oxo heme cation radical is formed that has been proposed to be the active oxidant. How and why it activates melatonin in the human body and what the origin of the product distributions is, are unknown. This encouraged us to do a detailed computational study on a typical human P450 isozyme, namely CYP1A1. We initially did a series of molecular dynamics simulations with substrate docked into several orientations. These simulations reveal a number of stable substrate-bound positions in the active site, which may lead to differences in substrate activation channels. Using tunneling analysis on the full protein structures, we show that two of the four binding conformations lead to open substrate-binding pockets. As a result, in these open pockets, the substrate is not tightly bound and can escape back into the solution. In the closed conformations, in contrast, the substrate is mainly oriented with the methoxy group pointing toward the heme, although under a different angle. We then created large quantum cluster models of the enzyme and focused on the chemical reaction mechanisms for melatonin activation, leading to competitive O-demethylation and C⁶-aromatic hydroxylation pathways. The calculations show that active site positioning determines the product distributions, but the bond that is activated is not necessarily closest to the heme in the enzyme-substrate complex. As such, the docking and molecular dynamics positioning of the substrate versus oxidant can give misleading predictions on product distributions. In particular, in quantum mechanics cluster model I, we observe that through a tight hydrogen bonding network, a preferential 6-hydroxylation of melatonin is obtained. However, O-demethylation becomes possible in alternative substrate-binding orientations that have the C⁶-aromatic ring position shielded. Finally, we investigated enzymatic and non-enzymatic O-demethylation processes and show that the hydrogen bonding network in the substrate-binding pocket can assist and perform this step prior to product release from the enzyme.

Theoretical analysis of pressure induced spin crossover phenomenon in a di-nuclear Fe(II) molecular complex

We have studied a Fe-based di-nuclear molecular complex having the chemical formula [{Fe(bpp)(NCS)₂}₂([Formula: see text]'-bipy)]·2MeOH (where bpp  =  [Formula: see text]-bis(pyrazol-3-yl) pyridine and [Formula: see text]'-bipy  =  [Formula: see text]'-bipyridine, 1) using density functional theory and model Hamiltonian approach. Our study provides insight to the pressure driven spin-crossover (SCO) phenomena observed experimentally in these systems. Upon increasing the pressure, the spin state of Fe(II) cation gradually changes from a high spin state (S  =2) to a low spin (LS) state (S  =0) accompanied by volume contraction. The gradual increase in pressure shrinks Fe-N bond length and also causes angular deviation of the FeN₆ octahedron leading to full conversion to the LS state without global structural phase transition. We have carried out exact diagonalization study of an effective single site Hamiltonian and confirmed the importance of intramolecular interaction for SCO phenomena. We have investigated the cooperativity of the observed SCO phenomena. We have also studied the effect of Co doping on the spin state of Fe and find that the spin state of Fe has a subtle dependency on the concentration of dopant atoms. Excess Co doping pave the way towards the possibility of an intermediate spin state for Fe and can give rise to a bistable spin transition process.

QM Calculations in ADMET Prediction

In recent years, there has been an increase in the application of quantum mechanics (QM) methods to describe properties related to the ADMET profile of small molecules. The application of these methods allows calculating useful descriptors and physiochemical properties contributing to ADMET prediction. Considering that QM methods are the only one that describe the electronic state of a molecules, such methods are particularly useful for studying the metabolism of drugs; furthermore, the introduction of mixed QM and molecular mechanics (QM/MM) is also increasing the understanding of drug interaction with cytochromes from a mechanistic point of view. Finally, combining the increase number of experimental data with machine learning algorithms and QM-derived descriptors allowed the creation of an end-user software capable of affecting the drug discovery process.

Mapping hole hopping escape routes in proteins

A recently proposed oxidative damage protection mechanism in proteins relies on hole hopping escape routes formed by redox-active amino acids. We present a computational tool to identify the dominant charge hopping pathways through these residues based on the mean residence times of the transferring charge along these hopping pathways. The residence times are estimated by combining a kinetic model with well-known rate expressions for the charge-transfer steps in the pathways. We identify the most rapid hole hopping escape routes in cytochrome P450 monooxygenase, cytochrome c peroxidase, and benzylsuccinate synthase (BSS). This theoretical analysis supports the existence of hole hopping chains as a mechanism capable of providing hole escape from protein catalytic sites on biologically relevant timescales. Furthermore, we find that pathways involving the [4Fe4S] cluster as the terminal hole acceptor in BSS are accessible on the millisecond timescale, suggesting a potential protective role of redox-active cofactors for preventing protein oxidative damage.