Molecular modeling of cytochrome P450 3A4

The three-dimensional structure of human cytochrome P450 3A4 was modeled based on crystallographic coordinates of four bacterial P450s: P450 BM-3, P450cam, P450terp, and P450eryF. The P450 3A4 sequence was aligned to those of the known proteins using a structure-based alignment of P450 BM-3, P450cam, P450terp, and P450eryF. The coordinates of the model were then calculated using a consensus strategy, and the final structure was optimized in the presence of water. The P450 3A4 model resembles P450 BM-3 the most, but the B helix is similar to that of P450eryF, which leads to an enlarged active site when compared with P450 BM-3, P450cam, and P450terp. The 3A4 residues equivalent to known substrate contact residues of the bacterial proteins and key residues of rat P450 2B1 are located in the active site or the substrate access channel. Docking of progesterone into the P450 3A4 model demonstrated that the substrate bound in a 6-orientation can interact with a number of active site residues, such as 114, 119, 301, 304, 305, 309, 370, 373, and 479, through hydrophobic interactions. The active site of the enzyme can also accommodate erythromycin, which, in addition to the residues listed for progesterone, also contacts residues 101, 104, 105, 214, 215, 217, 218, 374, and 478. The majority of 3A4 residues which interact with progesterone and/or erythromycin possess their equivalents in key residues of P450 2B enzymes, except for residues 297, 480 and 482, which do not contact either substrate in P450 3A4. The results from docking of progesterone and erythromycin into the enzyme model make it possible to pinpoint residues which may be important for 3A4 function and to target them for site-directed mutagenesis.     

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.     

Three-dimensional modelling of human cytochrome P450 1A2 and its interaction with caffeine and MeIQ

The three-dimensional modelling of proteins is a useful tool to fill the gap between the number of sequenced proteins and the number of experimentally known 3D structures. However, when the degree of homology between the protein and the available 3D templates is low, model building becomes a difficult task and the reliability of the results depends critically on the correctness of the sequence alignment. For this reason, we have undertaken the modelling of human cytochrome P450 1A2 starting by a careful analysis of several sequence alignment strategies (multiple sequence alignments and the TOPITS threading technique). The best results were obtained using TOPITS followed by a manual refinement to avoid unlikely gaps. Because TOPITS uses secondary structure predictions, several methods that are available for this purpose (Levin, Gibrat, DPM, NnPredict, PHD, SOPM and NNSP) have also been evaluated on cytochromes P450 with known 3D structures. More reliable predictions on -helices have been obtained with PHD, which is the method implemented in TOPITS. Thus, a 3D model for human cytochrome P450 1A2 has been built using the known crystal coordinates of P450 BM3 as the template. The model was refined using molecular mechanics computations. The model obtained shows a consistent location of the substrate recognition segments previously postulated for the CYP2 family members. The interaction of caffeine and a carcinogenic aromatic amine (MeIQ), which are characteristic P450 1A2 substrates, has been investigated. The substrates were solvated taking into account their molecular electrostatic potential distributions. The docking of the solvated substrates in the active site of the model was explored with the AUTODOCK programme, followed by molecular mechanics optimisation of the most interesting complexes. Stable complexes were obtained that could explain the oxidation of the considered substrates by cytochrome P450 1A2 and could offer an insight into the role played by water molecules.     

Substrate binding favors enhanced NO binding to P450cam

Ferric cytochrome P450cam from Pseudomonas putida (P450cam) in buffer solution at physiological pH 7.4 reversibly binds NO to yield the nitrosyl complex P450cam(NO). The presence of 1R-camphor affects the dynamics of NO binding to P450cam and enhances the association and dissociation rate constants significantly. In the case of the substrate-free form of P450cam, subconformers are evident and the NO binding kinetics are much slower than in the presence of the substrate. The association and dissociation processes were investigated by both laser flash photolysis and stopped-flow techniques at ambient and high pressure. Large and positive values of S and V observed for NO binding to and release from the substrate-free P450cam complex are consistent with the operation of a limiting dissociative ligand substitution mechanism, where the lability of coordinated water dominates the reactivity of the iron(III)-heme center with NO. In contrast, NO binding to P450cam in the presence of camphor displays negative activation entropy and activation volume values that support a mechanism dominated by a bond formation process. Volume profiles for the binding of NO appear to be a valuable approach to explain the differences observed for P450cam in the absence and presence of the substrate and enable the clarification of the underlying reaction mechanisms at a molecular level. Changes in spin state of the iron center during the binding/release of NO contribute significantly to the observed volume effects. The results are discussed in terms of relevance for the biological function of cytochrome P450 and in context to other investigations of the related reactions between NO and imidazole- and thiolate-ligated iron(III) hemoproteins.     

The origin of the low-spin character of the resting state of cytochrome P450cam investigated by means of active site analogues

Crown-capped iron(S-) porphyrins 1 x H2O and 2 x H2O and their corresponding Ba2+ complexes have been prepared as active site analogues of the resting state of cytochrome P450cam. cw-EPR studies and electronic structure calculations at the density functional theory (DFT) level of model systems suggest a functional role of the water cluster of P450cam.     

A virtual high throughput screen for high affinity cytochrome P450cam substrates. Implications for in silico prediction of drug metabolism

Structure-based virtual screening techniques require reliable scoring functions to discriminate potential substrates effectively. In this study we compared the performance of GOLD, PMF, DOCK and FlexX scoring functions in FlexX flexible docking to cytochrome P450cam binding site. Crystal structures of protein-substrate complexes were most effectively reproduced by the FlexX/PMF method. On the other hand, the FlexX/GOLD approach provided the best correlation between experimental binding constants and predicted scores. Binding modes selected by the FlexX/PMF approach were rescored by GOLD to obtain a reliable measure of binding energetics. The effectiveness of the FlexX/PMF/GOLD method was demonstrated by the correct classification of 32 out of the 33 experimentally studied compounds and also in a virtual HTS test on a library of 10,000 compounds. Although almost all the available functions were developed to be general, our study on cytochrome P450cam substrates suggests that careful selection or even tailoring the scoring function might increase the prediction power of virtual screens significantly. The FlexX/PMF/GOLD methodology was tested on cytochrome P450 3A4 substrates and inhibitors. This preliminary study revealed that the combined function was able to recognise 334 out of the 345 compounds bound to 3A4.     

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.     

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.     

Engineering cytochrome P-450s

Cytochrome P-450 isozymes represent a critical component of nature’s spectrum of detoxification catalysts that could be exploited for bioremediation. The ethanol-inducible human cytochrome P-450 2E1 serves as a model eukaryotic P-450 that complements the bacterial P-450 cam in dehalogenation and detoxification of environmental pollutants. We explored the construction of novel chimeric P-450s using cytochrome P-450 camC and 2E1 genes. For construction of chimera 1 (478 amino acids, 55.14 kDa), 145 amino acids from the N-terminus of P-450 2E1 protein (493 amino acids, 56.84 kDa) were replaced with 130 amino acids from the N-terminus of P-450 camC protein (415 amino acids, 46.66 kDa). In chimera 2 (525 amino acids, 60.24 kDa) the strategy involves replacement of 28 amino acids in the C-terminus of chimera 1 with 75 amino acids from the C-terminus of P-450 camC gene. Homology models of both the chimeric proteins were developed using SWISS-MODEL based on the known crystal structure of cytochrome P-450 camC, BM-3, 1DT6A, and 2C17A. The models indicated that the proposed heme-binding site was intact, which is inevitable for catalytic activity of cytochrome P-450s. The expression of chimera 1 and 2 genes in Escherichia coli DH5α was evident from light-pink cell pellets, protein band in sodium dodecyl sulfate polyacrylamide gel electrophoresis, and diagnostic carbon monoxide-difference spectra. Our studies show that strategies can be developed to exploit the natural diversity of the P-450 superfamily to generate chimeric biocatalysts that would provide new templates amenable to directed evolution.     

4D-QSAR analysis of a series of antifungal p450 inhibitors and 3D-pharmacophore comparisons as a function of alignment

A training set of 55 antifungal p450 analogue inhibitors was used to construct receptor-independent four-dimensional quantitative structure-activity relationship (RI 4D-QSAR) models. Ten different alignments were used to build the models, and one alignment yields a significantly better model than the other alignments. Two different methodologies were used to measure the similarity of the best 4D-QSAR models of each alignment. One method compares the residual of fit between pairs of models using the cross-correlation coefficient of their residuals of fit as a similarity measure. The other method compares the spatial distributions of the IPE types (3D-pharmacophores) of pairs of 4D-QSAR models from different alignments. Optimum models from several different alignments have nearly the same correlation coefficients, r(2), and cross-validation correlation coefficients, xv-r(2), yet the 3D-pharmacophores of these models are very different from one another. The highest 3D-pharmacophore similarity correlation coefficient between any pair of 4D-QSAR models from the 10 alignments considered is only 0.216. However, the best 4D-QSAR models of each alignment do contain some proximate common pharmacorphore sites. A test set of 10 compounds was used to validate the predictivity of the best 4D-QSAR models of each alignment. The "best" model from the 10 alignments has the highest predictivity. The inferred active sites mapped out by the 4D-QSAR models suggest that hydrogen bond interactions are not prevalent when this class of P450 analogue inhibitors binds to the receptor active site. This feature of the 4D-QSAR models is in agreement with the crystal structure results that indicate no ligand-receptor hydrogen bonds are formed.     

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.     

A Compound I Mimic Reveals the Transient Active Species of a Cytochrome P450 Enzyme: Insight into the Stereoselectivity of P450-Catalysed Oxidations

Catching the structure of cytochrome P450 enzymes in flagrante is crucial for the development of P450 biocatalysts, as most structures collected are found trapped in a precatalytic conformation. At the heart of P450 catalysis lies Cpd I, a short-lived, highly reactive intermediate, whose recalcitrant nature has thwarted most attempts at capturing catalytically relevant poses of P450s. We report the crystal structure of P450BM3 mimicking the state in the precise moment preceding epoxidation, which is in perfect agreement with the experimentally observed stereoselectivity. This structure was attained by incorporation of the stable Cpd I mimic oxomolybdenum mesoporphyrin IX into P450BM3 in the presence of styrene. The orientation of styrene to the Mo-oxo species in the crystal structures sheds light onto the dynamics involved in the rotation of styrene to present its vinyl group to Cpd I. This method serves as a powerful tool for predicting and modelling the stereoselectivity of P450 reactions.     

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.     

The elusive oxidant species of cytochrome P450 enzymes: characterization by combined quantum mechanical/molecular mechanical (QM/MM) calculations

The primary oxidant of cytochrome P450 enzymes, Compound I, is hard to detect experimentally; in the case of cytochrome P450(cam), this intermediate does not accumulate in solution during the catalytic cycle even at temperatures as low as 200 K (ref 4). Theory can play an important role in characterizing such elusive species. We present here combined quantum mechanical/molecular mechanical (QM/MM) calculations of Compound I of cytochrome P450(cam) in the full enzyme environment as well as density functional studies of the isolated QM region. The calculations assign the ground state of the species, quantify the effect of polarization and hydrogen bonding on its properties, and show that the protein environment and its specific hydrogen bonding to the cysteinate ligand are crucial for sustaining the Fe-S bond and for preventing the full oxidation of the sulfur.     

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.     

Fluorescent probes for cytochrome p450 structural characterization and inhibitor screening

We have synthesized two luminescent probes (D-4-Ad and D-8-Ad) that target cytochrome P450cam. D-4-Ad luminescence is quenched by F?rster energy transfer upon binding (Kd = 0.83 muM) but is restored when the probe is displaced from the active site by camphor. In contrast, D-8-Ad (Kd approximately 0.02 muM) is not displaced from the enzyme, even in the presence of a large excess of camphor. The 2.2 A resolution crystal structure of the D-8-Ad:P450cam complex reveals extensive hydrophobic contacts between the probe and the enzyme, which result from the conformational flexibility of the B', F, and G helices. Probes with properties similar to those of D-4-Ad potentially could be useful for screening P450 inhibitors.     

Kinetic isotope effects implicate the iron-oxene as the sole oxidant in P450-catalyzed N-dealkylation

Multiple oxidants have been implicated as playing a role in cytochrome P450-mediated oxidations. Herein, we report results on N-dealkylation, one of the most facile reactions mediated by P450 enzymes. We have employed the N-oxides of a series of para-substituted 13C2H2-labeled N,N-dimethylanilines to function as both substrates and surrogate oxygen atom donors for P450cam and P4502E1. Kinetic isotope effect profiles obtained using the N-oxide system were found to closely match the profiles produced using the complete NAD(P)H/NAD(P)-P450 reductase/O2 system. The results are consistent with oxidation occurring solely through an iron-oxene species.     

Evidence for basic ferryls in cytochromes P450

Using a combination of M?ssbauer spectroscopy and density functional calculations, we have determined that the ferryl forms of P450(BM3) and P450cam are protonated at physiological pH. Density functional calculations were performed on large active-site models of these enzymes to determine the theoretical M?ssbauer parameters for the ferryl and protonated ferryl (Fe(IV)OH) species. These calculations revealed a significant enlargement of the quadrupole splitting parameter upon protonation of the ferryl unit. The calculated quadrupole splittings for the protonated and unprotonated ferryl forms of P450(BM3) are DeltaE(Q) = 2.17 mm/s and DeltaE(Q) = 1.05 mm/s, respectively. For P450cam, they are DeltaE(Q) = 1.84 mm/s and DeltaE(Q) = 0.66 mm/s, respectively. The experimentally determined quadrupole splittings (P450(BM3), DeltaE(Q) = 2.16 mm/s; P450cam, DeltaE(Q) = 2.06 mm/s) are in good agreement with the values calculated for the protonated forms of the enzymes. Our results suggest that basic ferryls are a natural consequence of thiolate-ligated hemes.     

Cytochrome P450cam-monoterpene interactions

The use of the molecular mechanics AMBER force field (FF) to predict product profiles for the hydroxylation of the monoterpenes 1R-camphor, 1S-camphor, 1R-norcamphor, 1S-norcamphor and camphane by the enzyme cytochrome P450cam from the soil bacterium Pseudomonas putida was investigated. Predictions were carried out by applying multiple substrate (starting) orientations in the enzyme pocket in two procedures: a procedure based on molecular dynamics (MD) and a procedure based on short MD simulations followed by geometry optimisations. The latter (GO) procedure is faster and enabled the use of more monoterpene starting orientations. Monoterpene orientations were transformed into product profiles by applying both energetic and geometrical criteria appropriate for the (monoterpene) hydrogen abstraction reaction. Good predictions compared to experimental data were obtained for most compounds in both the MD and GO procedures. Prior to the product profile calculations, the FF was calibrated by reproducing the experimental data for the binding energy of 1R-camphor and 1S-camphor to P450cam and the energy of vaporisation of water. Focus of the calibration was on the value for the scaling factor for the electrostatic interactions.     

Protein side-chain motion and hydration in proton-transfer pathways. Results for cytochrome p450cam

Proton-transfer reactions form an integral part of bioenergetics and enzymatic catalysis. The identification of proton-conducting pathways inside a protein is a key to understanding the mechanisms of biomolecular proton transfer. Proton pathways are modeled here as hydrogen bonded networks of proton-conducting groups, including proton-exchanging groups of amino acid side chains and bound water molecules. We focus on the identification of potential proton-conducting pathways inside a protein of known structure. However, consideration of the static structure alone is often not sufficient to detect suitable proton-transfer paths, leading, for example, from the protein surface to the active site buried inside the protein. We include dynamic fluctuations of amino acid side chains and water molecules into our analysis. To illustrate the method, proton transfer into the active site of cytochrome P450cam is studied. The cooperative rotation of amino acids and motion of water molecules are found to connect the protein surface to the molecular oxygen. Our observations emphasize the intrinsic dynamical nature of proton pathways where critical connections in the network may be transiently provided by mobile groups.