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.     

Protein dynamics in cytochrome P450 molecular recognition and substrate specificity using 2D IR vibrational echo spectroscopy

Cytochrome (cyt) P450s hydroxylate a variety of substrates that can differ widely in their chemical structure. The importance of these enzymes in drug metabolism and other biological processes has motivated the study of the factors that enable their activity on diverse classes of molecules. Protein dynamics have been implicated in cyt P450 substrate specificity. Here, 2D IR vibrational echo spectroscopy is employed to measure the dynamics of cyt P450(cam) from Pseudomonas putida on fast time scales using CO bound at the active site as a vibrational probe. The substrate-free enzyme and the enzyme bound to both its natural substrate, camphor, and a series of related substrates are investigated to explicate the role of dynamics in molecular recognition in cyt P450(cam) and to delineate how the motions may contribute to hydroxylation specificity. In substrate-free cyt P450(cam), three conformational states are populated, and the structural fluctuations within a conformational state are relatively slow. Substrate binding selectively stabilizes one conformational state, and the dynamics become faster. Correlations in the observed dynamics with the specificity of hydroxylation of the substrates, the binding affinity, and the substrates' molecular volume suggest that motions on the hundreds of picosecond time scale contribute to the variation in activity of cyt P450(cam) toward different substrates.     

Surface-enhanced Raman scattering as a tool to probe cytochrome P450-catalysed substrate oxidation

Surface-enhanced Raman scattering was used as a spectroscopic tool to investigate the changes brought upon cytochrome P450BSß after fatty acid binding. Differences in the spectra of substrate-free and substrate-bound enzyme were observed indicating the potential for this method to be used in the screening of P450 substrates. In particular, the binding characteristics of myristic acid, an inherent substrate, and hydroxylauric acid, a product of fatty acid oxidation, towards P450BSß in the presence of H₂O₂ were investigated. Specific spectral changes could be assigned to changes in the heme environment only for myristic acid, indicating an occurrence of oxidation process characteristic for the enzymatic substrate.     

Structural and dynamical basis of broad substrate specificity, catalytic mechanism, and inhibition of cytochrome P450 3A4

Cytochrome P450 (CYP) 3A4 is responsible for the oxidative degradation of more than 50% of clinically used drugs. By means of molecular dynamics simulations with the newly developed force field parameters for the heme-thiolate group and its dioxygen adduct, we examine the differences in structural and dynamic properties between CYP3A4 in the resting form and its complexes with the substrate progesterone and the inhibitor metyrapone. The results indicate that the broad substrate specificity of CYP3A4 stems from the malleability of a loop (residues 211-218) that resides in the vicinity of the channel connecting the active site and bulk solvent. However, the high-amplitude motion of the flexible loop is found to be damped out upon binding of the inhibitor or the substrate in the active site. In the resting form of CYP3A4, a structural water molecule is bound to the sixth coordination position of the heme iron, stabilizing the octahedral coordination geometry. In addition to the direct coordination of metyrapone to the heme iron, the hydrogen bond interaction between the inhibitor carbonyl group and the side chain of Ser119 also contributes significantly to stabilizing the CYP3A4-metyrapone complex. On the other hand, progesterone is stabilized in the active site by the formation of two hydrogen bonds with Ser119 and Arg106, as well as by the van der Waals interactions with the heme and hydrophobic residues. The structural and dynamic features of the CYP3A4-progesterone complex indicate that the oxidative degradation of progesterone occurs through hydroxylation at the C16 position by the reactive oxygen coordinated to the heme iron.     

Exploring the substrate specificity of OxyB, a phenol coupling P450 enzyme involved in vancomycin biosynthesis

OxyB is a cytochrome P450 enzyme that catalyzes the first oxidative phenol coupling reaction during vancomycin biosynthesis. The preferred substrate is a linear peptide linked as a C-terminal thioester to a peptide carrier protein (PCP) domain of the glycopeptide antibiotic non-ribosomal peptide synthetase. Previous studies have shown that OxyB can efficiently oxidize a model hexapeptide-PCP conjugate (R-Leu(1)-R-Tyr(2)-S-Asn(3)-R-Hpg(4)-R-Hpg(5)-S-Tyr(6)-S-PCP) (Hpg = 4-hydroxyphenylglycine) into a macrocyclic product by phenolic coupling of the aromatic rings in residues-4 and -6. In this work, the substrate specificity of OxyB has been explored using a series of N-terminally truncated peptides related in sequence to this model hexapeptide-PCP conjugate. Deletion of one or three residues from the N-terminus afforded a penta- (Ac-Tyr-Asn-Hpg-Hpg-Tyr-S-PCP) and a tri- (Ac-Hpg-Hpg-Tyr-S-PCP) peptide that were also efficiently transformed into the corresponding macrocyclic cross-linked product by OxyB. The tripeptide, representing the core of the macrocycle in vancomycin created by OxyB, is thus sufficient, as a thioester with the PCP domain, for phenol coupling to occur. The related tetrapeptide-PCP thioester was not cyclized by OxyB, neither was a related model hexapeptide containing tryptophan in place of tyrosine-6, nor were tripeptides (related to the natural product K-13) with the sequence Ac-Tyr-Tyr-Tyr-S-PCP cross-linked by OxyB.     

Probing the cytochrome P450-like reactivity of high-valent oxo iron intermediates in the gas phase

Electrospray ionization in combination with Fourier transform ion cyclotron resonance spectrometry is used to prepare and characterize at a molecular level high-valent oxoiron intermediates formed in the reaction of [(TPFPP)Fe(III)]Cl (TPFPP= meso-tetrakis(pentafluorophenyl)porphinato dianion) (1-Cl) with H(2)O(2) in methanol. The intrinsic reactivity in the gas phase of the iron(IV) oxo porphyrin cation radical complex, [(TPFPP)(.+)Fe(IV)=O](+), has been probed toward selected substrates (S), chosen among naturally occurring and model compounds. Whereas CO and cyclohexane proved to be unreactive, olefins, sulfides, amines, and phosphites all undergo oxygen atom transfer in the gas phase yielding the reduction product 1 and/or an adduct ion ([1-S](+)). The reaction efficiencies show a qualitative correlation with the oxophilic character of the active site of S. A notable exception is nitric oxide, which displays a remarkably high reactivity, in line with the important role of NO reactions with iron porphyrin complexes. Furthermore, subsidiary information on the neat association reaction of 1 with selected ligands (L) has been obtained by a kinetic study showing that both the efficiency and the extent of ligation toward the naked ion 1 depend on the electron-donating ability of L.     

Using a homology model of cytochrome P450 2D6 to predict substrate site of metabolism

CYP2D6 is an important enzyme that is involved in first pass metabolism and is responsible for metabolizing ~25% of currently marketed drugs. A homology model of CYP2D6 was built using X-ray structures of ligand-bound CYP2C5 complexes as templates. This homology model was used in docking studies to rationalize and predict the site of metabolism of known CYP2D6 substrates. While the homology model was generally found to be in good agreement with the recently solved apo (ligand-free) X-ray structure of CYP2D6, significant differences between the structures were observed in the B′ and F–G helical region. These structural differences are similar to those observed between ligand-free and ligand-bound structures of other CYPs and suggest that these conformational changes result from induced-fit adaptations upon ligand binding. By docking to the homology model using Glide, it was possible to identify the correct site of metabolism for a set of 16 CYP2D6 substrates 85% of the time when the 5 top scoring poses were examined. On the other hand, docking to the apo CYP2D6 X-ray structure led to a loss in accuracy in predicting the sites of metabolism for many of the CYP2D6 substrates considered in this study. These results demonstrate the importance of describing substrate-induced conformational changes that occur upon binding. The best results were obtained using Glide SP with van der Waals scaling set to 0.8 for both the receptor and ligand atoms. A discussion of putative binding modes that explain the distribution of metabolic sites for substrates, as well as a relationship between the number of metabolic sites and substrate size, are also presented. In addition, analysis of these binding modes enabled us to rationalize the typical hydroxylation and O-demethylation reactions catalyzed by CYP2D6 as well as the less common N-dealkylation.     

Structural fine-tuning of a multifunctional cytochrome P450 monooxygenase

AurH is a unique cytochrome P450 monooxygenase catalyzing the stepwise formation of a homochiral oxygen heterocycle, a key structural and pharmacophoric component of the antibiotic aureothin. The exceptional enzymatic reaction involves a tandem oxygenation process including a regio- and stereospecific hydroxylation, followed by heterocyclization. For the structural and biochemical basis of this unparalleled sequence, four crystal structures of AurH variants in different conformational states and in complex with the P450 inhibitor ancymidol were solved, which represent the first structures of the CYP151A group. Structural data in conjunction with computational docking, site-directed mutagenesis, and chemical analyses unveiled a switch function when recognizing the two substrates, deoxyaureothin and the hydroxylated intermediate, thus allowing the second oxygenation-heterocyclization step. Furthermore, we were able to modify the chemo- and regioselectivity of AurH, yielding mutants that catalyze the regioselective six-electron transfer of a nonactivated methyl group to a carboxylic acid via hydroxyl and aldehyde intermediates.     

Preparation, characterization, and substrate metabolism of gold-immobilized cytochrome P450 2C9

The cytochrome P450 enzymes represent an important class of heme-containing enzymes. There is considerable interest in immobilizing these enzymes on a surface so that interactions between a single enzyme and other species can be studied with respect to electron transfer, homodimer or heterodimer interactions, or for construction of biological-based chips for standardizing cytochrome P450 metabolism or for high-throughput screening of pharmaceutical agents. Previous studies have generally immobilized P450 enzymes in a matrix or on a surface. Here, we have attached CYP2C9 to gold substrates such that the resulting construct maintains the ability to bind and metabolize substrates in the presence of NADPH and cytochrome P450 reductase. The activity of these chips is directly dependent upon the linkers used to attach CYP2C9 and to the presence of key molecules in the active site during enzyme attachment. A novel method to detect substrate-enzyme binding, namely, superconducting quantum interference device (SQUID) magnetometry, was used to monitor the binding of substrates. Most significantly, conditions that allow measurable CYP2C9 metabolism to occur have been developed.     

Identification of broad specificity P450CAM variants by primary screening against indole as substrate

High-throughput screening of cytochrome P450CAM libraries, for their ability to oxidise indole to indigo and indirubin, has resulted in the identification of variants with activity towards the structurally unrelated substrate diphenylmethane.     

Direct electrochemistry of enzymes from the cytochrome P450 2C family

The human cytochromes P450 are responsible for the clearance of ∼90% of xenobiotics yet comparatively little is known about their electrochemistry. Here we report the first direct electrochemistry of P450s from the 2C subfamily; one of the major groups of enzymes from this family. Specifically, the proteins that we have examined are recombinant human P450s 2C9, 2C18 and 2C19 and reversible Fe^III/II couples are seen in the absence of dioxygen. Even in the presence of trace amounts of dioxygen, a pronounced cathodic response is seen which is assigned to catalytic reduction of the bound dioxygen ligand by the ferrous P450.     

细胞色素P450的电化学研究进展

细胞色素P450的电化学研究从一个侧面反映了为使细胞色素P450达到工业催化剂的最终目的人们所作的不懈努力。本文从细胞色素P450在电极上的电子转移研究,隧道扫描显微镜的微观成像研究和使用电极作为细胞色素P450的电子给体从而实现细胞色素P450底物转化三方面,评述了近年来细胞色素P450的电化学研究进展。     

Synthesis of highly enantioenriched hydroxy- and dihydroxy-fatty esters: substrate precursors for cytochrome P450BioI

A series of highly enantioenriched hydroxy- and dihydroxy-fatty esters were required as part of our ongoing investigation into cytochrome P450_BioI. This mediates the biosynthesis of pimelic acid via C–C bond cleavage of long chain fatty acids within Bacillus subtilis. Herein we report the synthesis of various stereoisomers of methyl 7-hydroxytetradecanoate, methyl 8-hydroxytetradecanoate, and methyl 7,8-dihydroxytetradecanoate in highly enantioenriched form, using a combination of asymmetric synthesis and a preparative enantioselective HPLC is reported.     

Oxidizing species in the mechanism of cytochrome P450

This review discusses the mechanisms of oxygen activation by cytochrome P450 enzymes, the possible catalytic roles of the various iron--oxygen species formed in the catalytic cycle, and progress in understanding the mechanisms of hydrocarbon hydroxylation, heteroatom oxidation, and olefin epoxidation. The focus of the review is on recent results, but earlier work is discussed as appropriate. The literature through to February 2002 is surveyed, and 175 referenced are cited.     

Stochastic simulations of the cytochrome P450 catalytic cycle

Cytochrome P450 enzymes involve a complex reaction cycle which has been described, for the first time, by a stochastic simulation of the system in the present work. A series of models are developed for a basic catalytic cycle, employing a set of microscopic rate constants for the oxidation of p-alkoxyacylanilides catalyzed by the cytochrome P450 1A2. By analyzing the effects of low concentrations of enzyme and substrate on the system, and the dependence of the system on several rate constants, it is discovered that the system evolves along relatively stable patterns from its initial state, as indicated from different runs of simulations. Strong fluctuations appear at the entrance and exit of the pathway, with very weak fluctuations in the middle sections of the cycle. Although noises are apparent when the reactant populations are very low, basically, the fundamental feature of the P450 cycle based on a microscopic view is that it is deterministic in nature. Meanwhile, the mathematical models we developed are qualitatively validated by a comparison with those experimental results of the P450 cycle. The findings of this work will be helpful for a further deeper understanding of the catalytic mechanism of cytochrome P450 enzymes.     

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.