Proton-transfer dynamics in the activation of cytochrome P450eryF

Molecular dynamics simulations are combined with quantum chemistry calculations of instantaneous proton-transfer energy profiles to investigate proton-transfer events in the transient pathway of cytochrome P450eryF (6-deoxyerythronolide B hydroxylase; CYP107A1), from the oxyferrous species to the catalytically active ferryl oxygen species (compound I). This reaction is one of the most fundamental unresolved aspects in the mechanism of oxidation that is common to all cytochrome P450s. We find that this process involves an ultrafast proton transfer from the crystallographic water molecule W519 to the distal oxygen bound to the heme group, and a subsequent proton-transfer event from W564 to W519. Both proton-transfer events are found to be endothermic in the oxyferrous state, suggesting that the oxyferrous reduction is mechanistically linked to the proton-transfer dynamics. These findings indicate that the hydrogen bond network, proximate to the O(2)-binding cleft, plays a crucial functional role in the enzymatic activation of P450s. Our results are consistent with the effect of mutations on the enzymatic efficacy.     

Fluorescence detection of ligand binding to labeled cytochrome P450 BM3

The cytochrome P450 superfamily of monoxygenases are highly relevant for pharmaceutical, environmental and biocatalytical applications. The binding of a substrate to their catalytic site is usually detectable by UV-vis spectroscopy as a low-to-high spin state transition of the heme iron. However, the discovery of potential new substrates is limited by the fact that some compounds do not cause the typical spin-shift even if they are oxidised by P450 enzymes. Here we report a fluorescence-based method able to detect the binding of such substrates to the heme domain of cytochrome P450 BM3 from Bacillus megaterium. The protein was labeled with the fluorescent probe N,N'-dimethyl-N-(iodoacetyl)-N'-(7-nitrobenz-2-oxa-1,3-diazol-4-Yl)-ethylenediamine (IANBD). Arachidonic and lauric acids are substrates of P450 BM3 and were used to validate the method, as their binding can be detected both by a spin-shift of the Soret peak from 419 to 397 nm and by the fluorescence change of the labelled protein. The fluorescence emission of the probe linked to the protein increased by a value corresponding to 121 ± 9% and 52 ± 5% with respect to the initial one, upon titration with arachidonic or lauric acids respectively. The dissociation constants were calculated by both UV-vis and fluorescence spectroscopy. Three drugs, propranolol, chlorzoxazone and nifedipine, known to be oxidized by P450 BM3 and that bind without causing spin-shift, were also tested and the fluorescence emission of IANBD was found to decrease by 29 ± 5%, 21 ± 2% and 23 ± 3%, respectively, allowing the measurement of their dissociation constants.     

Observation of ligand binding to cytochrome P450 BM-3 by means of solid-state NMR spectroscopy

A broad understanding of the binding modes of ligands and inhibitors to cytochrome P450 is vital for the development of new drugs. We investigated ligand binding in a site-specific fashion on cytochrome P450 BM-3 from Bacillus megaterium, a 119 kDa paramagnetic enzyme, using solid-state magic angle spinning nuclear magnetic resonance methods. Selective labeling and longitudinal relaxation effects were utilized to identify the peaks in a site-specific fashion and to provide evidence for binding. Well-resolved one-dimensional and two-dimensional NMR spectra of cytochrome P450 BM-3 reveal shifts upon binding of its substrate, N-palmitoylglycine. These data are consistent with the crystallographic result that a biochemically important amino acid residue, Phe87, moves upon ligation. This experimental scheme provides a tool for probing ligand binding for complex systems.     

Proximal ligand electron donation and reactivity of the cytochrome P450 ferric-peroxo anion

CYP125 from Mycobacterium tuberculosis catalyzes sequential oxidation of the cholesterol side-chain terminal methyl group to the alcohol, aldehyde, and finally acid. Here, we demonstrate that CYP125 simultaneously catalyzes the formation of five other products, all of which result from deformylation of the sterol side chain. The aldehyde intermediate is shown to be the precursor of both the conventional acid metabolite and the five deformylation products. The acid arises by protonation of the ferric-peroxo anion species and formation of the ferryl-oxene species, also known as Compound I, followed by hydrogen abstraction and oxygen transfer. The deformylation products arise by addition of the same ferric-peroxo anion to the aldehyde intermediate to give a peroxyhemiacetal that leads to C-C bond cleavage. This bifurcation of the catalytic sequence has allowed us to examine the effect of electron donation by the proximal ligand on the properties of the ferric-peroxo anion. Replacement of the cysteine thiolate iron ligand by a selenocysteine results in UV-vis, EPR, and resonance Raman spectral changes indicative of an increased electron donation from the proximal selenolate ligand to the iron. Analysis of the product distribution in the reaction of the selenocysteine substituted enzyme reveals a gain in the formation of the acid (Compound I pathway) at the expense of deformylation products. These observations are consistent with an increase in the pK(a) of the ferric-peroxo anion, which favors its protonation and, therefore, Compound I formation.     

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.     

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

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

Analysis of the cryptophycin P450 epoxidase reveals substrate tolerance and cooperativity

Cryptophycins are potent anticancer agents isolated from Nostoc sp. ATCC 53789 and Nostoc sp. GSV 224. The most potent natural cryptophycin analogues retain a beta-epoxide at the C2'-C3' position of the molecule. A P450 epoxidase encoded by c rpE recently identified from the cryptophycin gene cluster was shown to install this key functional group into cryptophycin-4 (Cr-4) to produce cryptophycin-2 (Cr-2) in a regio- and stereospecific manner. Here we report a detailed characterization of the CrpE epoxidase using an engineered maltose binding protein (MBP)-CrpE fusion. The substrate tolerance of the CrpE polypeptide was investigated with a series of structurally related cryptophycin analogues generated by chemoenzymatic synthesis. The enzyme specifically installed a beta-epoxide between C2' and C3' of cyclic cryptophycin analogues. The kcat/Km values of the enzyme were determined to provide further insights into the P450 epoxidase catalytic efficiency affected by substrate structural variation. Finally, binding analysis revealed cooperativity of MBP-CrpE toward natural and unnatural desepoxy cryptophycin substrates.     

Configurational entropy and cooperativity between ligand binding and dimerization in glycopeptide antibiotics

Oligomerization and ligand binding are thermodynamically cooperative processes in many biochemical systems, and the mechanisms giving rise to cooperative behavior are generally attributed to changes in structure. In glycopeptide antibiotics, however, these cooperative processes are not accompanied by significant structural changes. To investigate the mechanism by which cooperativity arises in these compounds, fully solvated molecular dynamics simulations and quasiharmonic normal-mode analysis were performed on chloroeremomycin, vancomycin, and dechlorovancomycin. Configurational entropies were derived from the vibrational modes recovered from ligand-free and ligand-bound forms of the monomeric and dimeric species. Results indicate that both ligand binding and dimerization incur an entropic cost as vibrational activity in the central core of the antibiotic is shifted to higher frequencies with lower amplitudes. Nevertheless, ligand binding and dimerization are cooperative because the entropic cost of both processes occurring together is less than the cost of these processes occurring separately. These reductions in configurational entropy are more than sufficient in magnitude to account for the experimentally observed cooperativity between dimerization and ligand binding. We conclude that biochemical cooperativity can be mediated through changes in vibrational activity, irrespective of the presence or absence of concomitant structural change. This may represent a general mechanism of allostery underlying cooperative phenomena in diverse macromolecular systems.     

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.     

Investigation of the proton-assisted pathway to formation of the catalytically active, ferryl species of P450s by molecular dynamics studies of P450eryF

The recently determined crystal structure of cytochrome P450eryF (6-deoxyerythronolide B hydroxylase; CYP107A1) in its ferric heme substrate-bound form has been used to address one of the most fundamental unresolved aspects of the mechanism of oxidation common to this ubiquitous family of metabolizing heme proteins, the pathway from the twice reduced dioxygen species to the putative catalytically active ferryl oxygen species. Both of these species are too transient to have been characterized experimentally, and the transformation from one to the other has been only partially characterized. The observed requirement of two protons and the formation of water in this transformation suggests a proton-assisted dioxygen bond cleavage as a plausible pathway. However, this pathway is difficult to establish by experiment alone, and the source of the protons in the largely hydrophobic binding pocket of the P450s remains unclear. In this work we have performed molecular dynamics simulations of the twice reduced dioxygen substrate-bound form of this isozyme in order to (i) determine the plausibility of the proposed pathway to compound I formation, a proton-assisted cleavage of the dioxygen bond, and (ii) investigate the possible source of these protons. The analysis of the molecular dynamics trajectories of this species does indeed provide further evidence for this pathway and points to a source of protons. Specifically, two dynamically stable hydrogen bonds to the distal oxygen atom of the dioxygen ligand, one by the substrate and the other by a bound water, are found, consistent with the proposed proton-assisted cleavage of the bond and formation of water. In addition, an extensive dynamically stable hydrogen bond network is formed that connects the distal oxygen to Glu 360, a well-conserved residue in a channel accessible to solvent that could be the ultimate source of protons. The simulations were done for both a protonated and unprotonated Glu and led to a proposed mechanism of proton transfer by it to the distal oxygen atom. In order to validate the procedures used for the simulation of this transient twice-reduced species, we have used these same procedures to perform molecular dynamics simulations of two other forms of P450eryF, the ferric and ferryl substrate-bound species, and compared the results with experiment. The results for the ferric substrate-bound species were assessed by comparisons to the experimentally determined X-ray structure and fluctuations, and good agreement was found. The simulations performed for the ferryl substrate-bound species led to the correct prediction of the observed regio- and stereospecific hydroxylation of its natural substrate, 6-deoxyerythronolide B (6-DEB) at the 6S position. The results of these two additional studies lend credibility to the important mechanistic inferences from the simulations of the transient twice reduced dioxygen species: further evidence for a proton-assisted pathway from it to the catalytically active ferryl species and a possible source of the protons.     

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.     

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.     

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.     

Engineering and assaying of cytochrome P450 biocatalysts

Cytochrome P450s constitute a highly fascinating superfamily of enzymes which catalyze a broad range of reactions. They are essential for drug metabolism and promise industrial applications in biotechnology and biosensing. The constant search for cytochrome P450 enzymes with enhanced catalytic performances has generated a large body of research. This review will concentrate on two key aspects related to the identification and improvement of cytochrome P450 biocatalysts, namely the engineering and assaying of these enzymes. To this end, recent advances in cytochrome P450 development are reported and commonly used screening methods are surveyed.     

Structural dynamics of the cooperative binding of organic molecules in the human cytochrome P450 3A4

Cytochrome P450 3A4 (CYP3A4) is a key enzyme responsible for the metabolism of 50% of all orally administered drugs which exhibit an intriguing kinetic behavior typified by a sigmoidal dependence of the reaction velocity on the substrate concentration. There is evidence for the binding of two substrates in the active site of the enzyme, but the mechanism of this cooperative binding is unclear. Diazepam is such a drug that undergoes metabolism by CYP3A4 with sigmoidal dependence. Metabolism is initiated by hydrogen atom abstraction from the drug. To understand the factors that determine the cooperative binding and the juxtaposition of the C-H bond undergoing abstraction, we carried out molecular dynamics simulations for two enzymatic conformers and examined the differences between the substrate-free and the bound enzymes, with one and two diazepam molecules. Our results indicate that the effector substrate interacts both with the active substrate and with the enzyme, and that this interaction results in side chain reorientation with relatively minor long-range effects. In accord with experiment, we find that F304, in the interface between the active and effector binding sites, is a key residue in the mechanism of cooperative binding. The addition of the effector substrate stabilizes F304 and its environment, especially F213, and induces a favorable orientation of the active substrate, leading to a short distance between the targeted hydrogen for abstraction and the active species of the enzyme. In addition, in one conformer of the enzyme, residue R212 may strongly interact with F304 and counteract the effector's impact on the enzyme.     

The role of thermodynamics and kinetics in ligand binding to G-quadruplex DNA

Molecular dynamics simulations were used to investigate the binding of four different 2,4,6-triarylpyridines to G-quadruplex DNA. Both the binding free energies, and the kinetics of binding are required to explain the measured degree of ligand induced stabilisation of the compounds, with bulky substituents having the potential to prevent the ligand from reaching the lowest energy binding site.     

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.