Molecular recognition in (+)-alpha-pinene oxidation by cytochrome P450cam

Oxygenated derivatives of the monoterpene (+)-alpha-pinene are found in plant essential oils and used as fragrances and flavorings. (+)-alpha-Pinene is structurally related to (+)-camphor, the natural substrate of the heme monooxygenase cytochrome P450(cam) from Pseudomonas putida. The aim of the present work was to apply the current understanding of P450 substrate binding and catalysis to engineer P450(cam) for the selective oxidation of (+)-alpha-pinene. Consideration of the structures of (+)-camphor and (+)-alpha-pinene lead to active-site mutants containing combinations of the Y96F, F87A, F87L, F87W, and V247L mutations. All mutants showed greatly enhanced binding and rate of oxidation of (+)-alpha-pinene. Some mutants had tighter (+)-alpha-pinene binding than camphor binding by the wild-type. The most active was the Y96F/V247L mutant, with a (+)-alpha-pinene oxidation rate of 270 nmol (nmol of P450(cam))(-)(1) min(-)(1), which was 70% of the rate of camphor oxidation by wild-type P450(cam). Camphor is oxidized by wild-type P450(cam) exclusively to 5-exo-hydroxycamphor. If the gem dimethyl groups of (+)-alpha-pinene occupied similar positions to those found for camphor in the wild-type structure, (+)-cis-verbenol would be the dominant product. All P450(cam) enzymes studied gave (+)-cis-verbenol as the major product but with much reduced selectivity compared to camphor oxidation by the wild-type. (+)-Verbenone, (+)-myrtenol, and the (+)-alpha-pinene epoxides were among the minor products. The crystal structure of the Y96F/F87W/V247L mutant, the most selective of the P450(cam) mutants initially examined, was determined to provide further insight into P450(cam) substrate binding and catalysis. (+)-alpha-Pinene was bound in two orientations which were related by rotation of the molecule. One orientation was similar to that of camphor in the wild-type enzyme while the other was significantly different. Analysis of the enzyme/substrate contacts suggested rationalizations of the product distribution. In particular competition rather than cooperativity between the F87W and V247L mutations and substrate movement during catalysis were proposed to be major factors. The crystal structure lead to the introduction of the L244A mutation to increase the selectivity of pinene oxidation by further biasing the binding orientation toward that of camphor in the wild-type structure. The F87W/Y96F/L244A mutant gave 86% (+)-cis-verbenol and 5% (+)-verbenone. The Y96F/L244A/V247L mutant gave 55% (+)-cis-verbenol but interestingly also 32% (+)-verbenone, suggesting that it may be possible to engineer a P450(cam) mutant that could oxidize (+)-alpha-pinene directly to (+)-verbenone. Verbenol, verbenone, and myrtenol are naturally occurring plant fragrance and flavorings. The preparation of these compounds by selective enzymatic oxidation of (+)-alpha-pinene, which is readily available in large quantities, could have applications in synthesis. The results also show that the protein engineering of P450(cam) for high selectivity of substrate oxidation is more difficult than achieving high substrate turnover rates because of the subtle and dynamic nature of enzyme-substrate interactions.     

Biotransformation of the sesquiterpene (+)-valencene by cytochrome P450cam and P450BM-3

The sesquiterpenoids are a large class of naturally occurring compounds with biological functions and desirable properties. Oxidation of the sesquiterpene (+)-valencene by wild type and mutants of P450cam from Pseudomonas putida, and of P450BM-3 from Bacillus megaterium, have been investigated as a potential route to (+)-nootkatone, a fine fragrance. Wild type P450cam did not oxidise (+)-valencene but the mutants showed activities up to 9.8 nmol (nmol P450)(-1) min(-1), with (+)-trans-nootkatol and (+)-nootkatone constituting >85% of the products. Wild type P450BM-3 and mutants had higher activities (up to 43 min(-1)) than P450cam but were much less selective. Of the many products, cis- and trans-(+)-nootkatol, (+)-nootkatone, cis-(+)-valencene-1,10-epoxide, trans-(+)-nootkaton-9-ol, and (+)-nootkatone-13S,14-epoxide were isolated from whole-cell reactions and characterised. The selectivity patterns suggest that (+)-valencene has one binding orientation in P450cam but multiple orientations in P450BM-3.     

A scanning tunnelling study of immobilised cytochrome P450cam

A site-specifically engineered surface cysteine residue, located in a region where the haem moiety is closest to the surface, is used to anchor cytochrome P450cam enzyme molecules covalently to a gold electrode. More reproducibly ordered adsorption, at high coverage, occurs with this K344C mutant than with the wild-type enzyme. The subsequently formed close-packed monolayer arrays have been probed by scanning tunnelling microscopy under ambient conditions and under aqueous (buffered) solution at high resolution. Initial indications suggest that the immobilised enzyme is both electrochemically addressable and catalytically active.     

Predicting the product specificity and coupling of cytochrome P450cam

Summary We present an analysis of several molecular dynamics trajectories of substrate-bound cytochrome P450cam. Trajectories were calculated for the native substrate, camphor, as well as for the alternative substrates, norcamphor and thiocamphor. The system modeled consisted of the crystallographically resolved amino acids, the heme group with a single oxygen atom as the distal ligand, the bound substrate, and the crystallographic waters. These trajectories of the presumptive ferryl oxygen intermediate were used to predict regiospecificity of hydroxylation and coupling between NADH consumption and product formation. Simple geometric criteria in combination with electronic considerations were used to calculate the probability of hydroxylation at specific sites on the substrate. We found that for all the cases examined, the predicted product ratios were in good agreement with the experimentally observed values. We also determined that these simple geometric criteria can be used to predict the degree of coupling between NADH consumption and product formation for a given substrate, which was in good agreement with the experimental values.     

Mechanisms of reaction in cytochrome P450: Hydroxylation of camphor in P450cam

The fundamental nature of reactivity in cytochrome P450 enzymes is currently controversial. Modelling of bacterial P450cam has suggested an important role for the haem propionates in the catalysis, though this finding has been questioned. Understanding the mechanisms of this enzyme family is important both in terms of basic biochemistry and potentially in the prediction of drug metabolism. We have modelled the hydroxylation of camphor by P450cam, using combined quantum mechanics/molecular mechanics (QM/MM) methods. A set of reaction pathways in the enzyme was determined. We were able to pinpoint the source of the discrepancies in the previous results. We show that when a correct ionization state is assigned to Asp297, no spin density appears on the haem propionates and the protein structure in this region remains preserved. These results indicate that the haem propionates are not involved in catalysis.     

Active-site mobility inhibits reductive dehalogenation of 1,1,1-trichloroethane by cytochrome P450cam

Summary Recent studies by Wackett and co-workers have shown that cytochrome P450cam is capable of reductively dehalogenating hexachloroethane at a significant rate, but that no appreciable dehalogenation of 1,1,1-trichloroethane is observed. A growing body of evidence indicates that differences in intrinsic reactivity can not completely explain this observation. We therefore explored the possible role of differences in preferred binding orientation and in active-site mobility. A detailed analysis of molecular dynamics trajectories with each of these substrates bound at the active site of P450cam is presented. While the dynamics and overall time-average structure calculated for the protein are similar in the two trajectories, the two substrates behave quite differently. The smaller substrate, 1,1,1-trichloroethane, is significantly more mobile than hexachloroethane and has a preferred orientation in which the substituted carbon is generally far from the heme iron. In contrast, for hexachloroethane, one of the chlorine atoms is nearly always in van der Waals contact with the heme iron, which should favor the initial electron transfer step.     

Prediction of activation energies for aromatic oxidation by cytochrome P450

We have estimated the activation energy for aromatic oxidation by compound I in cytochrome P450 for a diverse set of 17 substrates using state-of-the-art density functional theory (B3LYP) with large basis sets. The activation energies vary from 60 to 87 kJ/mol. We then test if these results can be reproduced by computationally less demanding methods. The best methods (a B3LYP calculation of the activation energy of a methoxy-radical model or a partial least-squares model of the semiempirical AM1 bond dissociation energies and spin densities of the tetrahedral intermediate for both a hydroxyl-cation and a hydroxyl-radical model) give correlations with r(2) of 0.8 and mean absolute deviations of 3 kJ/mol. Finally, we apply these simpler methods on several sets of reactions for which experimental data are available and show that we can predict the reactive sites by combining calculations of the activation energies with the solvent-accessible surface area of each site.     

A role for Thr 252 in cytochrome P450cam oxygen activation

Molecular dynamics simulations of the ferrous dioxygen bound form of wild type cytochrome P450cam were performed and the results analyzed to reveal the time-dependent interactions of T252 with surrounding residues as well as with bound oxygen. The results indicate a time-dependent bimodal interaction of T252 with both G248 and the terminal oxygen of the bound dioxygen. The hydrogen bonding interaction of T252 with these two moieties is "anticorrelated" in the sense that the breaking of the T252-G248 hydrogen bond is concurrent with formation of the T252-dioxygen interaction. These simulations support the probability of a role of T252 in stabilization of the initial dioxygen bound complex and promotion of subsequent formation of compound I previously indicated by several experimental studies.     

Electron-transfer reactions and functionalization of cytochrome P450cam monooxygenase system in reverse micelles

Enzyme-based electron-transfer reactions involved in the cytochrome P450 monooxygenase system were investigated in nanostructural reverse micelles. A bacterial flavoprotein, putidaredoxin reductase (PdR), was activated and shown to be capable of catalyzing the electron transport from NADH to electron-carrier proteins such as cytochrome b5 (tCyt-b5) and putidaredoxin (Pdx) in reverse micelles. Ferric tCyt-b5 in reverse micelles was effectively converted to its ferrous form by the exogenous addition of separately prepared reverse micellar solution harboring PdR and NADH. The fact that direct interactions of macromolecular proteins should be possible in the reverse micellar system encouraged us to functionalize a multicomponent monooxygenase system composed of the bacterial cytochrome P450cam (P450cam), putidaredoxin (Pdx), and PdR in reverse micelles. The successful camphor hydroxylation reaction catalyzed by P450cam was significantly dependent on the coexistence of Pdx, PdR, and NADH but not H2O2, suggesting that the oxygen-transfer reactions proceeded via a "monooxygenation" mechanism. This is the first report of a multicomponent cytochrome P450 system exhibiting enzymatic activity in organic media.     

Enantioselective alpha-hydroxylation of 2-arylacetic acid derivatives and buspirone catalyzed by engineered cytochrome P450 BM-3

Here we report that an engineered microbial cytochrome P450 BM-3 (CYP102A subfamily) efficiently catalyzes the alpha-hydroxylation of phenylacetic acid esters. This P450 BM-3 variant also produces the authentic human metabolite of buspirone, R-6-hydroxybuspirone, with 99.5% ee.     

Hydrogen-bonding interactions in the active sites of cytochrome P450cam and its site-directed mutants.

Resonance Raman spectroscopy is applied to the cyanide adducts of cytochrome P450cam and its T252A and D251N site-directed mutants, both in their substrate-free and camphor-bound forms, to probe active-site heme structure and, in particular, interactions of the FeCN fragment with potential active-site H-bond donors. In contrast to the ferrous CO and ferric NO adducts, which form only essentially linear (slightly distorted) FeXY fragments, the spectra of the ferric CN(-) adducts provide clear evidence the for the existence of an additional, rather highly bent, conformer; that is, the cyanide complexes form both linear and bent conformers in both the substrate-free and substrate-bound forms. Formation of this bent conformer is most reasonably attributed to the presence of off-axis H-bond donors, which induce distortion on the FeCN fragment but not the FeCO and FeNO fragments, which are poorer H-bond acceptors. For all three proteins, the substrate-free form exhibits a complex spectral pattern which arises because one of the modes associated with the FeCN fragment is coupled with two heme macrocycle deformation modes. Significantly, no evidence for such coupling is observed in the spectra of the camphor-bound forms. While various unknown factors may possibly give rise to selective activation of such coupling in the substrate-free derivative, given the known facts about the active-site architecture of this enzyme, a plausible explanation is that the bent conformer is oriented toward the water-filled substrate-binding site in the substrate-free form, but oppositely, toward the proposed proton delivery shuttle, in the substrate-bound form. Sensitivity of the FeCN modes to H(2)O/D(2)O exchange in the two camphor-bound mutants, which is apparently absent for the camphor-bound native protein, is most reasonably attributed to the known presence of extra water in the active sites of these mutants.     

Comparison between electrochemistry/mass spectrometry and cytochrome P450 catalyzed oxidation reactions

The extent to which electrochemistry on-line with electrospray mass spectrometry can be used to mimic cytochrome P450 catalyzed oxidations has been investigated. Comparisons on the mechanistic level have been made for most reactions in an effort to explain why certain reactions can, and some cannot, be mimicked by electrochemical oxidations. The EC/MS/MS system used successfully mimics in cases where the P450 catalyzed reactions are supposed to proceed via a mechanism initiated by a one-electron oxidation, such as N-dealkylation, S-oxidation, P-oxidation, alcohol oxidation and dehydrogenation. The P450 catalyzed reactions initiated via direct hydrogen atom abstraction, such as O-dealkylation and hydroxylation of unsubstituted aromatic rings, generally had a too high oxidation potential to be electrochemically oxidized below the oxidation potential limit of water, and were not mimicked by the EC/MS/MS system. Even though the EC/MS/MS system is not able to mimic all oxidations performed by cytochrome P450, valuable information can be obtained concerning the sensitivity of the substrate towards oxidation and in which position of the molecule oxidations are likely to take place. For small-scale electrochemical synthesis of metabolites, starting from the drug, the EC/MS/MS system should be very useful for quick optimization of the electrochemical conditions. The simplicity of the system, and the ease and speed with which it can be applied to a large number of compounds, make it a useful tool in drug metabolism research.     

Electrochemical model of alkane oxidation by cytochrome P-450

Catalytic oxidation of alkanes to alcohols and ketones was shown to take place in an electrochemical cell with iron porphyrin deposited on a graphite cathode. The oxidation mechanism was assumed to be similar to that of cytochrome P-450 action.

---

, . , P-450.

     

Epoxidation of olefins by hydroperoxo-ferric cytochrome P450

The T252A mutant of cytochrome P450cam is unable to form the oxoferryl "active oxygen" intermediate, as judged by its inability to hydroxylate its normal substrate, camphor. In the present study, we demonstrate that T252A P450cam is nonetheless able to epoxidize olefins, due to the action of a second oxidant. However, as shown in earlier radiolytic studies and by the ability of T252A to reduce dioxygen to hydrogen peroxide, the mutant retains the ability to form the hydroperoxo-ferric reaction cycle intermediate. The present results provide strong evidence that hydroperoxo-ferric P450 can serve as a second electrophilic oxidant capable of olefin epoxidation.     

Nanosecond photoreduction of cytochrome p450cam by channel-specific Ru-diimine electron tunneling wires

We report the synthesis and characterization of Ru-diimine complexes designed to bind to cytochrome p450cam (CYP101). The sensitizer core has the structure [Ru(L(2))L'](2+), where L' is a perfluorinated biphenyl bridge (F(8)bp) connecting 4,4'-dimethylbipyridine to an enzyme substrate (adamantane, F(8)bp-Ad), a heme ligand (imidazole, F(8)bp-Im), or F (F(9)bp). The electron-transfer (ET) driving force (-deltaG degrees ) is varied by replacing the ancillary 2,2'-bipyridine ligands with 4,4',5,5'-tetramethylbipyridine (tmRu). The four complexes all bind p450cam tightly: Ru-F(8)bp-Ad (1, K(d) = 0.077 microM); Ru-F(8)bp-Im (2, K(d) = 3.7 microM); tmRu-F(9)bp (3, K(d) = 2.1 microM); and tmRu-F(8)bp-Im (4, K(d) = 0.48 microM). Binding is predominantly driven by hydrophobic interactions between the Ru-diimine wires and the substrate access channel. With Ru-F(8)bp wires, redox reactions can be triggered on the nanosecond time scale. Ru-wire 2, which ligates the heme iron, shows a small amount of transient heme photoreduction (ca. 30%), whereas the transient photoreduction yield for 4 is 76%. Forward ET with 4 occurs in roughly 40 ns (k(f) = 2.8 x 10(7) s(-)(1)), and back ET (Fe(II) --> Ru(III), k(b) approximately 1.7 x 10(8) s(-)(1)) is near the coupling-limited rate (k(max)). Direct photoreduction was not observed for 1 or 3. The large variation in ET rates among the Ru-diimine:p450 conjugates strongly supports a through-bond model of Ru-heme electronic coupling.     

Combined quantum mechanical/molecular mechanical study on the pentacoordinated ferric and ferrous cytochrome P450cam complexes

The pentacoordinated ferric and ferrous cytochrome P450(cam) complexes have been investigated by combined quantum mechanical/molecular mechanical (QM/MM) calculations in the presence of a protein/solvent environment and by QM calculations on the isolated QM regions with use of density functional theory. The B3LYP functional has been found more reliable than the BLYP and BHLYP functionals for estimating the relative state energies. The B3LYP/CHARMM calculations with an all-electron basis set for iron give high-spin ground states for the title complexes, in agreement with experiment. The comparison of the B3LYP/CHARMM results of the entire protein system with the B3LYP calculations on the naked QM regions shows that the amount of stabilization by the protein environment is largest for the intermediate-spin states, followed by the high-spin states of the complexes. The calculation of M?ssbauer parameters in the presence of the enzyme environment confirms the double occupation of the d(xz) orbital in the quintet spin state of the ferrous complex, consistent with the computed QM/MM energies in the enzyme environment, while the d(x)2(-)(y)2 orbital is doubly occupied in the gas-phase quintet state.