Crystals in Minutes: Instant On‐Site Microcrystallisation of Various Flavours of the CYP102A1 (P450BM3) Haem Domain

Despite CYP102A1 (P450BM3) representing one of the most extensively researched metalloenzymes, crystallisation of its haem domain upon modification can be a challenge. Crystal structures are indispensable for the efficient structure‐based design of P450BM3 as a biocatalyst. The abietane diterpenoid derivative N‐abietoyl‐l ‐tryptophan (AbiATrp) is an outstanding crystallisation accelerator for the wild‐type P450BM3 haem domain, with visible crystals forming within 2 hours and diffracting to a near‐atomic resolution of 1.22 Å. Using these crystals as seeds in a cross‐microseeding approach, an assortment of P450BM3 haem domain crystal structures, containing previously uncrystallisable decoy molecules and diverse artificial metalloporphyrins binding various ligand molecules, as well as heavily tagged haem‐domain variants, could be determined. Some of the structures reported herein could be used as models of different stages of the P450BM3 catalytic cycle.     

Investigating the Active Oxidants Involved in Cytochrome P450 Catalyzed Sulfoxidation Reactions

Cytochrome P450 (CYP) heme-thiolate monooxygenases catalyze the hydroxylation of the C−H bonds of organic molecules. This reaction is initiated by a ferryl-oxo heme radical cation (Cpd I). These enzymes can also catalyze sulfoxidation reactions and the ferric-hydroperoxy complex (Cpd 0) and the Fe(III)-H₂O₂ complex have been proposed as alternative oxidants for this transformation. To investigate this, the oxidation of 4-alkylthiobenzoic acids and 4-methoxybenzoic acid by the CYP199A4 enzyme from Rhodopseudomonas palustris HaA2 was compared using both monooxygenase and peroxygenase pathways. By examining mutants at the mechanistically important, conserved acid alcohol-pair (D251N, T252A and T252E) the relative amounts of the reactive intermediates that would form in these reactions were disturbed. Substrate binding and X-ray crystal structures helped to understand changes in the activity and enabled an attempt to evaluate whether multiple oxidants can participate in these reactions. In peroxygenase reactions the T252E mutant had higher activity towards sulfoxidation than O-demethylation but in the monooxygenase reactions with the WT enzyme the activity of both reactions was similar. The peroxygenase activity of the T252A mutant was greater for sulfoxidation reactions than the WT enzyme, which is the reverse of the activity changes observed for O-demethylation. The monooxygenase activity and coupling efficiency of sulfoxidation and oxidative demethylation were reduced by similar degrees with the T252A mutant. These observations infer that while Cpd I is required for O-dealkylation, another oxidant may contribute to sulfoxidation. Based on the activity of the CYP199A4 mutants it is proposed that this is the Fe(III)-H₂O₂ complex which would be more abundant in the peroxide-driven reactions.     

Characterizing the Intermediates Compound I and II in the Cytochrome P450 Catalytic Cycle with Nonlinear X‐ray Spectroscopy: A Simulation Study

Cytochrome P450 enzymes are an important family of biocatalysts that oxidize chemically inert C?H bonds. There are many unresolved questions regarding the catalytic reaction intermediates, in particular P450 Compound I (Cpd‐I) and II (Cpd‐II). By using simple molecular models, we simulate various X‐ray spectroscopy signals, including X‐ray absorption near‐edge structure (XANES), resonant inelastic X‐ray scattering (RIXS), and stimulated X‐ray Raman spectroscopy (SXRS) of the low‐ and high‐spin states of Cpd‐I and II. Characteristic peak patterns are presented and connected to the corresponding electronic structures. These X‐ray spectroscopy techniques are complementary to more conventional infrared and optical spectroscopy and they help to elucidate the evolving electronic structures of transient species along the reaction path.     

QM/MM study of the active species of the human cytochrome P450 3A4, and the influence thereof of the multiple substrate binding

Cytochrome P450 3A4 is involved in the metabolism of 50% of all swallowed drugs. The enzyme functions by means of a high-valent iron-oxo species, called compound I (Cpd I), which is formed after entrance of the substrate to the active site. We explored the features of Cpd I using hybrid quantum mechanical/molecular mechanical calculations on various models that are either substrate-free or containing one and two molecules of diazepam as a substrate. M?ssbauer parameters of Cpd I were computed. Our major finding shows that without the substrate, Cpd I tends to elongate its Fe-S bond, localize the radical on the sulfur, and form hydrogen bonds with A305 and T309, which may hypothetically lead to Cpd I consumption by H-abstraction. However, the positioning of diazepam close to Cpd I, as enforced by the effector molecule, was found to strengthen the NH...S interactions of the conserved I443 and G444 residues with the proximal cysteinate ligand. These interactions are known to stabilize the Fe-S bond, and as such, the presence of the substrate leads to a shorter Fe-S bond and it prevents the localization of the radical on the sulfur. This diazepam-Cpd I stabilization was manifested in the 1W0E conformer. The effector substrate did not influence Cpd I directly but rather by positioning the active substrate close to Cpd I, thus displacing the hydrogen bonds with A305 and T309, and thereby giving preference to substrate oxidation. It is hypothesized that these effects on Cpd I, promoted by the restrained substrate, may be behind the special metabolic behavior observed in cases of multiple substrate binding (also called cooperative binding). This restraint constitutes a mechanism whereby substrates stabilize Cpd I sufficiently long to affect monooxygenation by P450s at the expense of Cpd I destruction by the protein residues.     

An evaluation of molecular models of the cytochrome P450 Streptomyces griseolus enzymes P450SU1 and P450SU2

Summary P450SU1 and P450SU2 are herbicide-inducible bacterial cytochrome P450 enzymes from Streptomyces griseolus. They have two of the highest sequence identities to camphor hydroxylase (P450cam from Pseudomonas putida), the cytochrome P450 with the first known crystal structure. We have built several models of these two proteins to investigate the variability in the structures that can occur from using different modeling protocols. We looked at variability due to alignment methods, backbone loop conformations and refinement methods. We have constructed two models for each protein using two alignment algorithms, and then an additional model using an identical alignment but different loop conformations for both buried and surface loops. The alignments used to build the models were created using the Needleman-Wunsch method, adapted for multiple sequences, and a manual method that utilized both a dotmatrix search matrix and the Needleman-Wunsch method. After constructing the initial models, several energy minimization methods were used to explore the variability in the final models caused by the choice of minimization techniques. Features of cytochrome P450cam and the cytochrome P450 superfamily, such as the ferredoxin binding site, the heme binding site and the substrate binding site were used to evaluate the validity of the models. Although the final structures were very similar between the models with different alignments, active-site residues were found to be dependent on the conformations of buried loops and early stages of energy minimization. We show which regions of the active site are the most dependent on the particular methods used, and which parts of the structures seem to be independent of the methods.     

Active species of horseradish peroxidase (HRP) and cytochrome P450: two electronic chameleons

The active site of HRP Compound I (Cpd I) is modeled using hybrid density functional theory (UB3LYP). The effects of neighboring amino acids and of environmental polarity are included. The low-lying states have porphyrin radical cationic species (Por(*)(+)). However, since the Por(*)(+) species is a very good electron acceptor, other species, which can be either the ligand or side chain amino acid residues, may participate in electron donation to the Por(*)(+) moiety, thereby making Cpd I behave like a chemical chameleon. Thus, this behavior that was noted before for Cpd I of P450 is apparently much more wide ranging than initially appreciated. Since chemical chameleonic behavior property was found to be expressed not only in the properties of Cpd I itself, but also in its reactivity, the roots of this phenomenon are generalized. A comparative discussion of Cpd I species follows for the enzymes HRP, CcP, APX, CAT (catalase), and P450.     

Cytochrome P450 compound I in the plane wave pseudopotential framework: GGA electronic and geometric structure of thiolate‐ligated iron(IV)–oxo porphyrin

The cytochromes P450 constitute a ubiquitous family of metalloenzymes, catalyzing manifold reactions of biological and synthetic importance via a thiolate‐ligated iron‐oxo (IV) porphyrin radical species denoted compound I (Cpd I). Experimental investigations have implicated this intermediate in a broad spectrum of biophysically interesting phenomena, further augmenting the importance of a Cpd I model system. Ab initio molecular dynamics, including Car–Parrinello and path integral methods, conjoin electronic structure theory with finite temperature simulation, affording tools most valuable to approach such enzymes. These methods are typically driven by density functional theory (DFT) in a plane‐wave pseudopotential framework; however, existing studies of Cpd I have been restricted to localized Gaussian basis sets. The appropriate choice of density functional and pseudopotential for such simulations is accordingly not obvious. To remedy this situation, a systematic benchmarking of thiolate‐ligated Cpd I is performed using several generalized‐gradient approximation (GGA) functionals in the Martins–Troullier and Vanderbilt ultrasoft pseudopotential schemes. The resultant electronic and structural parameters are compared to localized–basis DFT calculations using GGA and hybrid density functionals. The merits and demerits of each scheme are presented in the context of reproducing existing experimental and theoretical results for Cpd I. © 2013 Wiley Periodicals, Inc.     

Engineering Cytochrome P450BM3 Enzymes for Direct Nitration of Unsaturated Hydrocarbons

Applications of the peroxidase activity of cytochrome P450 enzymes in synthetic chemistry remain largely unexplored. We present herein a protein engineering strategy to increase cytochrome P450BM3 peroxidase activity for the direct nitration of aromatic compounds and terminal aryl-substituted olefins in the presence of a dual-functional small molecule (DFSM). Site-directed mutations of key active-site residues allowed the efficient regulation of steric effects to limit substrate access and, thus, a significant decrease in monooxygenation activity and increase in peroxidase activity. Nitration of several phenol and aniline compounds also yielded ortho- and para-nitration products with moderate-to-high total turnover numbers. Besides direct aromatic nitration by P450 variants using nitrite as a nitrating agent, we also demonstrated the use of the DFSM-facilitated P450 peroxidase system for the nitration of the vinyl group of styrene and its derivatives.     

Spectroscopic and Kinetic Evidence for the Crucial Role of Compound 0 in the P450cam‐Catalyzed Hydroxylation of Camphor by Hydrogen Peroxide

The hydroperoxo iron(III) intermediate P450_camFe^III–OOH, being the true Compound 0 (Cpd 0) involved in the natural catalytic cycle of P450_cam, could be transiently observed in the peroxo‐shunt oxidation of the substrate‐free enzyme by hydrogen peroxide under mild basic conditions and low temperature. The prolonged lifetime of Cpd 0 enabled us to kinetically examine the formation and reactivity of P450_camFe^III–OOH species as a function of varying reaction conditions, such as pH, and concentration of H₂O₂, camphor, and potassium ions. The mechanism of hydrogen peroxide binding to the substrate‐free form of P450_cam differs completely from that observed for other heme proteins possessing the distal histidine as a general acid–base catalyst and is mainly governed by the ability of H₂O₂ to undergo deprotonation at the hydroxo ligand coordinated to the iron(III) center under conditions of pH≥p${K{{{\rm P450}\hfill \atop {\rm a}\hfill}}}$ . Notably, no spectroscopic evidence for the formation of either Cpd I or Cpd II as products of heterolytic or homolytic O?O bond cleavage, respectively, in Cpd 0 could be observed under the selected reaction conditions. The kinetic data obtained from the reactivity studies involving (1R)‐camphor, provide, for the first time, experimental evidence for the catalytic activity of the P450Fe^III–OOH intermediate in the oxidation of the natural substrate of P450_cam.     

QM/MM study of mechanisms for compound I formation in the catalytic cycle of cytochrome P450cam

In the catalytic cycle of cytochrome P450cam, after molecular oxygen binds as a ligand to the heme iron atom to yield a ferrous dioxygen complex, there are fast proton transfers that lead to the formation of the active species, Compound I (Cpd I), which are not well understood because they occur so rapidly. In the present work, the conversion of the ferric hydroperoxo complex (Cpd 0) to Cpd I has been investigated by combined quantum-mechanical/molecular-mechanical (QM/MM) calculations. The residues Asp(251) and Glu(366) are considered as proton sources. In mechanism I, a proton is transported to the distal oxygen atom of the hydroperoxo group via a hydrogen bonding network to form protonated Cpd 0 (prot-Cpd0: FeOOH(2)), followed by heterolytic O-O bond cleavage that generates Cpd I and water. Although a local minimum is found for prot-Cpd0 in the Glu(366) channel, it is very high in energy (more than 20 kcal/mol above Cpd 0) and the barriers for its decay are only 3-4 kcal/mol (both toward Cpd 0 and Cpd I). In mechanism II, an initial O-O bond cleavage followed by a concomitant proton and electron transfer yields Cpd I and water. The rate-limiting step in mechanism II is O-O cleavage with a barrier of about 13-14 kcal/mol. According to the QM/MM calculations, the favored low-energy pathway to Cpd I is provided by mechanism II in the Asp(251) channel. Cpd 0 and Cpd I are of similar energies, with a slight preference for Cpd I.     

Direct Hydroxylation of Benzene to Phenol by Cytochrome P450BM3 Triggered by Amino Acid Derivatives

The selective hydroxylation of benzene to phenol, without the formation of side products resulting from overoxidation, is catalyzed by cytochrome P450BM3 with the assistance of amino acid derivatives as decoy molecules. The catalytic turnover rate and the total turnover number reached 259 min⁻¹ P450BM3⁻¹ and 40 200 P450BM3⁻¹ when N‐heptyl‐l ‐proline modified with l ‐phenylalanine (C7‐l ‐Pro‐l ‐Phe) was used as the decoy molecule. This work shows that amino acid derivatives with a totally different structure from fatty acids can be used as decoy molecules for aromatic hydroxylation by wild‐type P450BM3. This method for non‐native substrate hydroxylation by wild‐type P450BM3 has the potential to expand the utility of P450BM3 for biotransformations.     

Compound I in heme thiolate enzymes: a comparative QM/MM study

This study directly compares the active species of heme enzymes, so-called Compound I (Cpd I), across the heme-thiolate enzyme family. Thus, sixty-four different Cpd I structures are calculated by hybrid quantum mechanical/molecular mechanical (QM/MM) methods using four different cysteine-ligated heme enzymes (P450(cam), the mutant P450(cam)-L358P, CPO and NOS) with varying QM region sizes in two multiplicities each. The overall result is that these Cpd I species are similar to each other with regard to many characteristic features. Hence, using the more stable CPO Cpd I as a model for P450 Cpd I in experiments should be a reasonable approach. However, systematic differences were also observed, and it is shown that NOS stands out in most comparisons. By analyzing the electrical field generated by the enzyme on the QM region, one can see that (a) the protein exerts a large influence and modifies all the Cpd I species compared with the gas-phase situation and (b) in NOS this field is approximately planar to the heme plane, whereas it is approximately perpendicular in the other enzymes, explaining the deviating results on NOS. The calculations on the P450(cam) mutant L358P show that the effects of removing the hydrogen bond between the heme sulfur and L358 are small at the Cpd I stage. Finally, Mossbauer parameters are calculated for the different Cpd I species, enabling future comparisons with experiments. These results are discussed in the broader context of recent findings of Cpd I species that exhibit large variations in the electronic structure due to the presence of the substrate.     

Gauging the relative oxidative powers of compound I, ferric-hydroperoxide, and the ferric-hydrogen peroxide species of cytochrome P450 toward C-H hydroxylation of a radical clock substrate

Density functional calculations were performed in response to the controversies regarding the identity of the oxidant species in cytochrome P450. The calculations were used to gauge the relative C-H hydroxylation reactivity of three potential oxidant species of the enzyme, the high-valent oxo-iron species Compound I (Cpd I), the ferric hydroperoxide Compound 0 (Cpd 0), and the ferric-hydrogen peroxide complex Fe(H(2)O(2)). The results for the hydroxylation of a radical probe substrate, 1, show the following trends: (a) Cpd I is the most reactive species; in its presence the other two reagents will be silent. (b) In the absence of Cpd I, substrate oxidation by Cpd 0 and Fe(H(2)O(2)) will take place via a stepwise mechanism that involves initial O-O homolysis followed by H-abstraction from 1. (c) Cpd 0 will undergo mostly porphyrin hydroxylation and only approximately 15% of substrate oxidation producing mostly the rearranged alcohol, 3 (Scheme 2). (d) Fe(H(2)O(2)) will generate mostly free hydrogen peroxide (uncoupling). A small fraction will perform substrate oxidation and lead mostly to 3. Reactivity probes for these reagents are kinetic isotope effect (KIE) and the product ratio of unrearranged to rearranged alcohols, [2/3]. Thus, for substrate oxidation by Cpd 0 or Fe(H(2)O(2)) KIE will be small, approximately 2, while Cpd I will have large KIE values. Typically both Cpd 0 and Fe(H(2)O(2)) will lead to a [2/3] ratio < 1, while Cpd I will lead to ratios > 1. In addition, the product isotope effect (KIE(2)/KIE(3) not equal 1) is expected from the reactivity of Cpd I.     

Enabling highly (R)-enantioselective epoxidation of styrene by engineering unique non-natural P450 peroxygenases

Unlike the excellent (S)-enantioselective epoxidation of styrene performed by natural styrene monooxygenases (ee > 99%), the (R)-enantioselective epoxidation of styrene has not yet achieved a comparable efficiency using natural or engineered oxidative enzymes. This report describes the H₂O₂-dependent (R)-enantioselective epoxidation of unfunctionalized styrene and its derivatives by site-mutated variants of a unique non-natural P450BM3 peroxygenase, working in tandem with a dual-functional small molecule (DFSM). The observed (R)-enantiomeric excess (ee) of styrene epoxidation is up to 99% with a turnover number (TON) of 918 by the best enantioselective mutant F87A/T268I/L181Q, while the best active mutant F87A/T268I/V78A/A184L (with 98% ee) gave a catalytic TON of 4350, representing the best activity of a P450 peroxygenase towards styrene epoxidation to date. Following this approach, a set of styrene derivatives, such as o-, m-, p-chlorostyrenes and fluorostyrenes, could also be epoxidized with modest to very good TONs (362–3480) and high (R)-enantioselectivities (95–99% ee). The semi-preparative scale synthesis of (R)-styrene oxide performed at 0 °C with high conversion, maintaining enantioselectivity, and moderate isolated yields, further suggests the potential application of the current P450 enzymatic system in styrene epoxidation. This study indicates that the synergistic use of protein engineering and an exogenous DFSM constitutes an efficient strategy to control the enantioselectivity of styrene epoxidation, thus substantially expanding the chemical scope of P450 enzymes as useful bio-oxidative catalysts.

H₂O₂-dependent epoxidation of unfunctionalized styrenes is achieved with high (R)-enantioselectivity and moderate to excellent TONs by combining site-mutated variants of cytochrome P450BM3 monooxygenase and a dual-functional small molecule (DFSM).     

New features in the catalytic cycle of cytochrome P450 during the formation of compound I from compound 0

Density functional theory (DFT) is applied to the dark section of the catalytic cycle of the enzyme cytochrome P450, namely, the formation of the active species, Compound I (Cpd I), from the ferric-hydroperoxide species (Cpd 0) by a protonation-assisted mechanism. The chosen 96-atom model includes the key functionalities deduced from experiment: Asp(251), Thr(252), Glu(366), and the water channels that relay the protons. The DFT model calculations show that (a) Cpd I is not formed spontaneously from Cpd 0 by direct protonation, nor is the process very exothermic. The process is virtually thermoneutral and involves a significant barrier such that formation of Cpd I is not facile on this route. (b) Along the protonation pathway, there exists an intermediate, a protonated Cpd 0, which is a potent oxidant since it is a ferric complex of water oxide. Preliminary quantum mechanical/molecular mechanical calculations confirm that Cpd 0 and Cpd I are of similar energy for the chosen model and that protonated Cpd 0 may exist as an unstable intermediate. The paper also addresses the essential role of Thr(252) as a hydrogen-bond acceptor (in accord with mutation studies of the OH group to OMe).     

Improved Cyclopropanation Activity of Histidine‐Ligated Cytochrome P450 Enables the Enantioselective Formal Synthesis of Levomilnacipran

Engineering enzymes capable of modes of activation unprecedented in nature will increase the range of industrially important molecules that can be synthesized through biocatalysis. However, low activity for a new function is often a limitation in adopting enzymes for preparative‐scale synthesis, reaction with demanding substrates, or when a natural substrate is also present. By mutating the proximal ligand and other key active‐site residues of the cytochrome P450 enzyme from Bacillus megaterium (P450‐BM3), a highly active His‐ligated variant of P450‐BM3 that can be employed for the enantioselective synthesis of the levomilnacipran core was engineered. This enzyme, BM3‐Hstar, catalyzes the cyclopropanation of N,N‐diethyl‐2‐phenylacrylamide with an estimated initial rate of over 1000 turnovers per minute and can be used under aerobic conditions. Cyclopropanation activity is highly dependent on the electronic properties of the P450 proximal ligand, which can be used to tune this non‐natural enzyme activity.     

A Panel of Cytochrome P450 BM3 Variants to Produce Drug Metabolites and Diversify Lead Compounds

Herein we demonstrate that a small panel of variants of cytochrome P450 BM3 from Bacillus megaterium covers the breadth of reactivity of human P450s by producing 12 of 13 mammalian metabolites for two marketed drugs, verapamil and astemizole, and one research compound. The most active enzymes support preparation of individual metabolites for preclinical bioactivity and toxicology evaluations. Underscoring their potential utility in drug lead diversification, engineered P450 BM3 variants also produce novel metabolites by catalyzing reactions at carbon centers beyond those targeted by animal and human P450s. Production of a specific metabolite can be improved by directed evolution of the enzyme catalyst. Some variants are more active on the more hydrophobic parent drug than on its metabolites, which limits production of multiply‐hydroxylated species, a preference that appears to depend on the evolutionary history of the P450 variant.     

Dual‐Functional Small Molecules for Generating an Efficient Cytochrome P450BM3 Peroxygenase

We report a unique strategy for the development of a H₂O₂‐dependent cytochrome P450BM3 system, which catalyzes the monooxygenation of non‐native substrates with the assistance of dual‐functional small molecules (DFSMs), such as N‐(ω‐imidazolyl fatty acyl)‐l ‐amino acids. The acyl amino acid group of DFSM is responsible for bounding to enzyme as an anchoring group, while the imidazolyl group plays the role of general acid–base catalyst in the activation of H₂O₂. This system affords the best peroxygenase activity for the epoxidation of styrene, sulfoxidation of thioanisole, and hydroxylation of ethylbenzene among those P450–H₂O₂ system previously reported. This work provides the first example of the activation of the normally H₂O₂‐inert P450s through the introduction of an exogenous small molecule. This approach improves the potential use of P450s in organic synthesis as it avoids the expensive consumption of the reduced nicotinamide cofactor NAD(P)H and its dependent electron transport system. This introduces a promising approach for exploiting enzyme activity and function based on direct chemical intervention in the catalytic process.     

Interplay between two spin states determines the hydroxylation catalyzed by P450 monooxygenase from Trichoderma brevicompactum

Tri11 (now renamed as tri22) encoded cytochrome P₄₅₀ monooxygenase in Trichoderma brevicompactum catalyzes the C-4 C-H hydroxylation of 12, 13-epoxytrichothec-9-ene (EPT) to produce trichodermol in the biosynthetic pathway of trichodermin/harzianum A. The density functional theory (DFT)-quantum mechanics (QM) approach is applied to elucidate the hydroxylation of EPT by using a model active species of P₄₅₀ (Cpd I). The QM calculations were performed on the active site complex, to find out transition-state structure, intermediate, and product complexes for the two spin states at different potential energy surfaces. The two state reactivity rebound-free product formation resulted from the interplay of two spin states (doublet and quartet).     

The stereochemistry of fatty acid hydroxylation by cytochrome P450BM3

The stereochemical preference for the cytochrome P450_BM3-catalysed hydroxylation of tetradecanoic and pentadecanoic acids has been determined via comparison with authentic non-racemic standards utilising enantioselective HPLC. The sub-terminal hydroxylation of these fatty acids by P450_BM3 is highly selective for the formation of the R-alcohols. This is the same enantioselectivity as is seen for hexadecanoic acid oxidation but contrasts with a previous report of S-hydroxylation of pentadecanoic acid by P450_BM3.