Common system setup for the entire catalytic cycle of cytochrome P450(cam) in quantum mechanical/molecular mechanical studies

We describe a system setup that is applicable to all species in the catalytic cycle of cytochrome P450(cam). The chosen procedure starts from the X-ray coordinates of the ferrous dioxygen complex and follows a protocol that includes the careful assignment of protonation states, comparison between different conceivable hydration schemes, and system preparation through a series of classical minimizations and molecular dynamics (MD) simulations. The resulting setup was validated by quantum mechanical/molecular mechanical (QM/MM) calculations on the resting state, the pentacoordinated ferric and ferrous complexes, Compound I, the transition state and hydroxo intermediate of the C--H hydroxylation reaction, and the product complex. The present QM/MM results are generally consistent with those obtained previously with individual setups. Concerning hydration, we find that saturating the protein interior with water is detrimental and leads to higher structural flexibility and catalytically inefficient active-site geometries. The MD simulations favor a low water density around Asp251 that facilitates side chain rotation of protonated Asp251 during the conversion of Compound 0 to Compound I. The QM/MM results for the two preferred hydration schemes (labeled SE-1 and SE-4) are similar, indicating that slight differences in the solvation close to the active site are not critical as long as camphor and the crystallographic water molecules preserve their positions in the experimental X-ray structures.     

Multireference ab initio quantum mechanics/molecular mechanics study on intermediates in the catalytic cycle of cytochrome P450(cam)

We have investigated the elusive reactive species of cytochrome P450(cam) (Compound I), the hydroxo complex formed during camphor hydroxylation, and the ferric hydroperoxo complex (Compound 0) by combined quantum mechanical/molecular mechanical (QM/MM) calculations, employing both density functional theory (DFT) and correlated ab initio methods. The first two intermediates appear multiconfigurational in character, especially in the doublet state and less so in the quartet state. DFT(B3LYP)/MM calculations reproduce the relative energies from correlated ab initio QM/MM treatments quite well, except for the splitting of the lowest A(1u)-A(2u) radical states. The inclusion of dynamic correlation is crucial for the proper ab initio treatment of these intermediates.     

Calculated electronic spectra of model ferric cytochrome P450 complexes

Restricted open-shell ground state properties and electronic spectra of two closely related low-spin, ferric, 6-coordinate, model cytochrome P450 complexes, one with a methyl mercaptide and the other a mercaptan as the second axial ligands, have been calculated with a newly modified, semiempirical INDO-SCF-CI method. The sensitivity of the calculated spectra to protonation of the sixth axial ligand, and the ability of the method to predict characteristic spectral features for the complexes investigated, are determined. Assignment of transitions, including xy- and z-polarized transitions, are made and compared with experimental observations where available. In particular, the origin of the anomalous split Soret spectrum observed in low-spin ferric complexes with mercaptide but not a mercaptan is investigated. Finally, a two part hypothesis is presented which provides a general explanation for the origin of both the observed split Soret and the red-shifted normal Soret in various ferrous and ferric P450 complexes in terms of the ground state orbital characters and simple symmetry considerations.     

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.     

Quantum and molecular mechanical study of the first proton transfer in the catalytic cycle of cytochrome P450cam and its mutant D251N

In the catalytic cycle of cytochrome P450cam, the hydroperoxo intermediate (Cpd 0) is formed by proton transfer from a reduced oxyheme complex (S5). This process is drastically slowed down when Asp251 is mutated to Asn (D251N). We report quantum mechanical/molecular mechanical (QM/MM) calculations that address this proton delivery in the doublet state through a hydrogen-bond network in the Asp251 channel, both for the wild-type enzyme and the D251N mutant, using four different active-site models. For the wild-type, we find a facile concerted mechanism for proton transfer from protonated Asp251 via Wat901 and Thr252 to the FeOO moiety, with a barrier of about 1 kcal/mol and a high exothermicity of more than 20 kcal/mol. In the D251N mutant with a neutral Asn251 residue, the proton transfer is almost thermoneutral or slightly exothermic in the three models considered. It is still very facile when the Asn251 residue adopts a conformation analogous to Asp251 in the wild-type enzyme, but the barrier increases significantly when the Asn251 side chain flips (as indicated by classical molecular dynamics simulations). This flip disrupts the hydrogen-bond network and hence the proton-transfer pathway, which causes a longer lifetime of S5 in the D251N mutant (consistent with experimental observations). The entry of an additional water molecule into the active site of D251N with flipped Asn251 regenerates the hydrogen-bond network and provides a viable mechanism for proton delivery in the mutant, with a moderate barrier of about 7 kcal/mol.     

Quantum mechanical/molecular mechanical investigation of the mechanism of C-H hydroxylation of camphor by cytochrome P450cam: theory supports a two-state rebound mechanism

The stereospecific cytochrome P450-catalyzed hydroxylation of the C(5)-H((5-exo)) bond in camphor has been studied theoretically by a combined quantum mechanical/molecular mechanical (QM/MM) approach. Density functional theory is employed to treat the electronic structure of the active site (40-100 atoms), while the protein and solvent environment (ca. 24,000 atoms) is described by the CHARMM force field. The calculated energy profile of the hydrogen-abstraction oxygen-rebound mechanism indicates that the reaction takes place in two spin states (doublet and quartet), as has been suggested earlier on the basis of calculations on simpler models ("two-state reactivity"). While the reaction on the doublet potential energy surface is nonsynchronous, yet effectively concerted, the quartet pathway is truly stepwise, including formation of a distinct intermediate substrate radical and a hydroxo-iron complex. Comparative calculations in the gas phase demonstrate the effect of the protein environment on the geometry and relative stability of intermediates (in terms of spin states and redox electromers) through steric constraints and electronic polarization.     

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.     

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.     

Butane and propane oxidation by engineered cytochrome P450cam

The haem monooxygenase cytochrome P450cam has been engineered to oxidise the gaseous alkanes butane and propane to butan-2-ol and propan-2-ol, respectively, by the use of bulky amino acid substitutions to reduce the volume of the substrate pocket and thus improve the enzyme-substrate fit: the F87W/Y96F/T101L/V247L mutant oxidizes butane with a turnover rate of 750 min-1 and 95% yield based on NADH consumed while the wild-type enzyme has an activity of 0.4 min-1 with 4% yield.     

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.     

Nitric oxide binding to ferric cytochrome P450: a computational study

The interaction between nitric oxide (NO) and the active site of ferric cytochrome P450 was studied by means of density functional theory (DFT), at the generalized gradient approximation level, and of the SAM1 semiempirical method. The electrostatic effects of the protein environment were included in our DFT scheme by using a hybrid quantum classical approach. The active-site model consisted of an iron(III) porphyrin, the adjacent cysteine residue, and one coordinated water molecule. For this system, spin populations and relative energies for selected spin states were computed. Interestingly, the unpaired electron density, the HOMO, and the LUMO were found to be highly localized on the iron and in an appreciable extent on the sulfur coordinated to the metal. This provides central information about the reactivity of nitric oxide with the active site. Since the substitution of a molecule of H2O by NO has been proposed as being responsible for the inhibition of the cytochrome in the presence of nitric oxide, we have analyzed the thermodynamic feasibility of the ligand exchange process. The structure of the nitrosylated active site was partially optimized using SAM1. A low-spin ground state was obtained for the nitrosyl complex, with a linear Fe-N-O angle. The trends found in Fe-N-O angles and Fe-N lengths of the higher energy spin states provided a notable insight into the electronic configuration of the complex within the framework of the Enemark and Feltham formalism. In relation to the protein environment, it was assessed that the electrostatic field has significant effects on several computed properties. However, in both vacuum and protein environments, the ligand exchange reaction turned out to be exergonic and the relative orders of spin states of the relevant species were the same.     

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.     

Cytochrome P450CAM enzymatic catalysis cycle: a quantum mechanics/molecular mechanics study

The catalytic pathway of cytochrome P450cam is studied by means of a hybrid quantum mechanics/molecular mechanics method. Our results reveal an active role of the enzyme in the different catalytic steps. The protein initially controls the energy gap between the high- and low-spin states in the substrate binding process, allowing thermodynamic reduction by putidaredoxin reductase and molecular oxygen addition. A second electron reduction activates the delivery of protons to the active site through a selective interaction of Thr252 and the distal oxygen causing the O--O cleavage. Finally, the protein environment catalyzes the substrate hydrogen atom abstraction step with a remarkably low free energy barrier ( approximately 8 kcal/mol). Our results are consistent with the effect of mutations on the enzymatic efficacy and provide a satisfactory explanation for the experimental failure to trap the proposed catalytically competent species, a ferryl Fe(IV) heme.     

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.     

Multi-state epoxidation of ethene by cytochrome P450: a quantum chemical study.

The epoxidation of ethene by a model for Compound I of cytochrome P450, studied by the use of density functional B3LYP calculations, involves two-state reactivity (TSR) with multiple electromer species, hence "multi-state epoxidation". The reaction is found to proceed in stepwise and effectively concerted manners. Several reactive states are involved; the reactant is an (oxo)iron(IV) porphyrin cation radical complex with two closely lying spin states (quartet and doublet), both of which react with ethene to form intermediate complexes with a covalent C-O bond and a carbon-centered radical (radical intermediates). The radical intermediates exist in two electromers that differ in the oxidation state of iron; Por(+)(*)Fe(III)OCH(2)CH(2)(*) and PorFe(IV)OCH(2)CH(2)(*) (Por = porphyrin). These radical intermediates exist in both the doublet- and quartet spin states. The quartet spin intermediates have substantial barriers for transformation to the quartet spin PorFe(III)-epoxide complex (2.3 kcal mol(-)(1) for PorFe(IV)OCH(2)CH(2)(*) and 7.2 kcal mol(-)(1) for Por(+)(*)Fe(III)OCH(2)CH(2)(*)). In contrast, the doublet spin radicals collapse to the corresponding PorFe(III)-epoxide complex with virtually no barriers. Consequently, the lifetimes of the radical intermediates are much longer on the quartet- than on the doublet spin surface. The loss of isomeric identity in the epoxide and rearrangements to other products arise therefore mostly, if not only, from the quartet process, while the doublet state epoxidation is effectively concerted (Scheme 7). Experimental trends are discussed in the light of the computed mechanistic scheme, and a comparison is made with closely related mechanistic schemes deduced from experiment.     

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.