Combined QM/MM calculations of active-site vibrations in binding process of P450cam to putidaredoxin

Combined QM/MM calculations of the active-site of cytochrome P450cam have been performed before and after the binding of P450cam to putidaredoxin. The calculations were carried out for both a 5-coordinated and a 6-coordinated active-site of cytochrome P450cam, with either a water molecule or a carbon monoxide molecule as a 6th distal ligand. An experimentally observed increase in the Fe-S stretching frequency that occurs after cytochrome P450cam binds to putidaredoxin, has been reproduced in our study. Experimentally observed changes in the Fe-C and C-O vibration frequencies that occur after binding of both proteins, have also been reproduced in our study. The computed increase of the Fe-S and Fe-C stretching frequencies is correlated with a corresponding decrease of the Fe-S and Fe-C interatomic distances. According to our calculations, for the active-site with carbon monoxide in the triplet electronic state, the binding process increases the spin densities on the iron and sulfur atoms, which changes the Fe-C and C-O stretching frequencies in opposite directions, in agreement with experimental data.     

Linking of cytochrome P450cam and putidaredoxin by a co-ordination bridge

Cytochrome P450_cam (CYP101) catalyzes the oxidation of D(+)-camphor at the 5 position. The enzyme couples the reduction of dioxygen to the oxidation of the substrate. To transfer electrons from the reductant (NADH) to the cytochrome, two additional proteins are required. These are putidaredoxin (PdX) and putidaredoxin reductase (PdR). We have chemically linked a form of PdX with a histidine tag at the C-terminus to the P450. To accomplish this, we have modified cysteine 334 on P450 with a bipyridinyl group, and co-ordinated the C-terminal histidine tag of PdX by the addition of Ni²⁺ or Ru³⁺. The Ru³⁺ complex was the most stable. The non-linked system gave mostly 5-ketocamphor, a product of two consecutive hydroxylations, and H₂O₂, a product of 2-electron uncoupling. The Ni²⁺ complex gave both 5-exo-hydroxycamphor and 5-ketocamphor, but it also uncoupled. The Ru³⁺ complex gave a single product (5-exo-hydroxycamphor) and did not uncouple at the optimal PdR concentration. Our results are consistent with other studies of this system, in that strong binding of PdX to P450 is crucial for good coupling and for release of 5-exo-hydroxycamphor.     

Probing ligand binding modes of human cytochrome P450 2J2 by homology modeling, molecular dynamics simulation, and flexible molecular docking

Cytochrome P450 (P450) 2J2 catalyzes epoxidation of arachidonic acid to eicosatrienoic acids, which are related to a variety of diseases such as coronary artery disease, hypertension, and carcinogenesis. Recent experimental data also suggest that P450 2J2 could be a novel biomarker and a potential target for cancer therapy. However, the active site topology and substrate specificity of this enzyme remain unclear. In this study, a three-dimensional model of human P450 2J2 was first constructed on the basis of the crystal structure of human P450 2C9 in complex with a substrate using homology modeling method, and refined by molecular dynamics simulation. Flexible docking approaches were then employed to dock four ligands into the active site of P450 2J2 in order to probe the ligand-binding modes. By analyzing the results, active site architecture and certain key residues responsible for substrate specificity were identified on the enzyme, which might be very helpful for understanding the enzyme's biological role and providing insights for designing novel inhibitors of P450 2J2.     

Molecular basis of the recognition of arachidonic acid by cytochrome P450 2E1 along major access tunnel

Cytochrome P450 2E1 is widely known for its ability to oxidize both low molecular weight xenobiotics and endogenous fatty acids (e.g., arachidonic acid (AA)). In this study, we investigated the structural features of the AA‐bound CYP2E1 complex utilizing molecular dynamics (MD) and found that the distinct binding modes for both AA and fatty acid analog are conserved. Moreover, multiple random acceleration MD simulations and steered MD simulations uncovered the most possible tunnel for fatty acids. The main attractions are derived from three key residues, His107, Ala108, and His109, whose side chains reorient to keep ligands bound via hydrogen bonds during the initial unbinding process. More importantly, based on the calculated binding free energy results, we hypothesize that the hydrogen bonds between the receptor and the ligand are the most important contributors involved in the binding affinity. Thus, it is inferred that the hydrogen bonds between these three residues and the ligand may help offer insights into the structural basis of the different ligand egress mechanisms for fatty acids and small weight compounds. Our investigation provides detailed atomistic insights into the structural features of human CYP2E1–fatty acid complex structures. Furthermore, the ligand‐binding characteristics obtained in the present study are helpful for both experimental and computational studies of CYPs and may allow future researchers to achieve desirable changes in enzymatic activities. © 2014 Wiley Periodicals, Inc. Biopolymers 103: 53–66, 2015.     

Structural evidence for direct hydride transfer from NADH to cytochrome P450nor

Nitric oxide reductase cytochrome P450nor catalyzes an unusual reaction, direct electron transfer from NAD(P)H to bound heme. Here, we succeeded in determining the crystal structure of P450nor in a complex with an NADH analogue, nicotinic acid adenine dinucleotide, which provides conclusive evidence for the mechanism of the unprecedented electron transfer. Comparison of the structure with those of dinucleotide-free forms revealed a global conformational change accompanied by intriguing local movements caused by the binding of the pyridine nucleotide. Arg64 and Arg174 fix the pyrophosphate moiety upon the dinucleotide binding. Stereo-selective hydride transfer from NADH to NO-bound heme was suggested from the structure, the nicotinic acid ring being fixed near the heme by the conserved Thr residue in the I-helix and the upward-shifted propionate side-chain of the heme. A proton channel near the NADH channel is formed upon the dinucleotide binding, which should direct continuous transfer of the hydride and proton. A salt-bridge network (Glu71-Arg64-Asp88) was shown to be crucial for a high catalytic turnover.     

Structural analysis of P450 AmphL from Streptomyces nodosus provides insights into substrate selectivity of polyene macrolide antibiotic biosynthetic P450s

AmphL is a cytochrome P450 enzyme that catalyzes the C8 oxidation of 8-deoxyamphotericin B to the polyene macrolide antibiotic, amphotericin B. To understand this substrate selectivity, we solved the crystal structure of AmphL to a resolution of 2.0 Å in complex with amphotericin B and performed molecular dynamics (MD) simulations. A detailed comparison with the closely related P450, PimD, which catalyzes the epoxidation of 4,5-desepoxypimaricin to the macrolide antibiotic, pimaricin, reveals key catalytic structural features responsible for stereo- and regio-selective oxidation. Both P450s have a similar access channel that runs parallel to the active site I helix over the surface of the heme. Molecular dynamics simulations of substrate binding reveal PimD can “pull” substrates further into the P450 access channel owing to additional electrostatic interactions between the protein and the carboxyl group attached to the hemiketal ring of 4,5-desepoxypimaricin. This substrate interaction is absent in AmphL although the additional substrate -OH groups in 8-deoxyamphotericin B help to correctly position the substrate for C8 oxidation. Simulations of the oxy-complex indicates that these -OH groups may also participate in a proton relay network required for O₂ activation as has been suggested for two other macrolide P450s, PimD and P450eryF. These findings provide experimentally testable models that can potentially contribute to a new generation of novel macrolide antibiotics with enhanced antifungal and/or antiprotozoal efficacy.     

Productive and non-productive complexes in cytochrome P450-containing systems

The equilibrium dissociation constants K_D, the complex association / dissociation rate constants (k _on /k _off) and lifetimes of the complexes of redox partners were measured for three cytochrome P450-containing monooxygenase systems (P450cam, P450scc, and P450 2B4) under hydroxylation conditions. The Q parameter representing the ratio of protein-protein complex lifetime (τ _lT ) to time required for a single hydroxylation cycle (τ_turnover) was introduced for estimation of productivity of complexes formed within the systems studied. The Q parameter was insignificantly changed upon transition from the oxidation to hydroxylation conditions. Lifetimes (τ _lT ) for the binary complexes formed within the P450cam and the P450scc systems obligatory requiring an intermediate electron transfer protein between the reductase and cytochrome P450 could not realize hydroxylation reactions for substrates with known τ_turnover and so they were non-productive while the binary complexes formed within the P450 2B4 system, not requiring such intermediate electron-transfer protein, appeared to be productive. Formation of ternary complexes was demonstrated under hydroxylation conditions in all three systems. Analysis of Q values led to the conclusion that the ternary complexes formed within the P450cam and the P450scc systems were productive. In the case of the P450 2B4 system, more than half (about 60%) ternary complexes were also found to be productive.     

Spectroelectrochemistry of cytochrome P450cam

The spectroelectrochemistry of camphor-bound cytochrome P450cam (P450cam) using gold electrodes is described. The electrodes were modified with either 4,4(')-dithiodipyridin or sodium dithionite. Electrolysis of P450cam was carried out when the enzyme was in solution, while at the same time UV-visible absorption spectra were recorded. Reversible oxidation and reduction could be observed with both 4,4(')-dithiodipyridin and dithionite modified electrodes. A formal potential (E(0')) of -373mV vs Ag/AgCl 1M KCl was determined. The spectra of P450cam complexed with either carbon monoxide or metyrapone, both being inhibitors of P450 catalysis, clearly indicated that the protein retained its native state in the electrochemical cell during electrolysis.     

P450cam biocatalysis in surfactant-stabilized two-phase emulsions

Cytochrome P450 monooxygenases (P450s) are powerful biocatalysts that have the ability to oxidize a broad range of substrates, often at non-reactive carbon centers. However, incorporation of P450s into synthetic schemes has so far been limited to a few whole-cell transformations. P450 substrates are often hydrophobic and have low water solubility, limiting the amount of product that can be produced. To help overcome this limitation, we have examined P450cam activity in two-phase hexane/water emulsions with and without the anionic surfactant, bis(2-ethylhexyl) sulfosuccinate sodium salt (AOT). Hydroxylation of camphor to hydroxycamphor by the three- component P450cam system was chosen as the model reaction, and regeneration of NADH was accomplished with yeast alcohol dehydrogenase (YADH). P450cam was activated in the surfactant-free emulsions, and addition of AOT improved the activity even further, at least over the range of camphor concentrations for which initial rates were readily measurable in all media. The largest observed rate enhancement was 4.5-fold. Nearly 50-times more product was formed in the surfactant-stabilized emulsions than was achieved in aqueous buffer, with total turnover numbers reaching 28,900 for P450cam and 11,800 for YADH. In the absence of surfactant, the two-phase reaction appeared to be mass-transfer limited, while inclusion of AOT alleviated transport limitations and/or afforded a larger interfacial area for P450 activation. The oxidation of hydroxycamphor to 2,5-diketocamphane was also observed, owing to the large concentration of hydroxycamphor relative to camphor in the aqueous phase of the two-phase emulsion. This competing reaction was accompanied by the uncoupled oxidation of NADH (i.e., NADH oxidation without formation of 2,5-diketocamphane), which reduced the availability of NADH for camphor oxidation and further limited the yield of hydroxycamphor in the two-phase emulsions. These results indicate that a surfactant-stabilized two-phase emulsion is a promising reaction medium for practical P450 biocatalysis, although its effectiveness for a given P450/substrate combination can depend on several factors, including competitive or sequential reactions, product inhibition, and NAD(P)H uncoupling.     

The solution structure of a gallium-substituted putidaredoxin mutant: GaPdx C85S

The Fe₂S₂ cluster of the ferredoxin putidaredoxin (Pdx) can be replaced by a single gallium ion, giving rise to a colorless, diamagnetic protein in which, apart from the metal binding site, the major structural features of the native ferredoxin are conserved. The solution structure of the C85S variant of gallium putidaredoxin (C85S GaPdx), in which a non-ligand cysteine is replaced by a serine, has been determined via multidimensional NMR methods using uniformly ¹⁵N,¹³C labeled samples of C85S GaPdx. Stereospecific assignments of leucine and valine methyl resonances were made using ¹³C,¹H HSQC spectra obtained with fractionally ¹³C-labeled samples, and backbone dihedral angle restraints were obtained using a combination of two-dimensional J-modulated ¹⁵N,¹H HSQC and three-dimensional (HN)CO(CO)NH experiments. A total of 1117 NOE-derived distance restraints were used in the calculations, including 454 short range ($i - j 3$), 456 long range (i - j 4) interresidue restraints and 207 non-trivial intraresidue restraints. 97 and 55 ₁ angular restraints were also included in the calculation of a family of 20 structures using a combined distance geometry-simulated annealing protocol. Most regions of the protein are well defined in the calculations, with an RMSD of 0.525 Å for backbone atoms excluding the metal binding loop (residues 34–48) and the last three C-terminal residues (residues 103–106). Where comparison is possible, these regions show an increase in dynamic behavior over the native protein, as does the loop containing residues 74–76. Structural and dynamic differences between native Pdx and GaPdx are discussed in relation to charge and packing of the metal binding site.     

Crystal structure of putidaredoxin reductase from Pseudomonas putida, the final structural component of the cytochrome P450cam monooxygenase

The crystal structure of recombinant putidaredoxin reductase (Pdr), an FAD-containing NADH-dependent flavoprotein component of the cytochrome P450cam monooxygenase from Pseudomonas putida, has been determined to 1.90 A resolution. The protein has a fold similar to that of disulfide reductases and consists of the FAD-binding, NAD-binding, and C-terminal domains. Compared to homologous flavoenzymes, the reductase component of biphenyl dioxygenase (BphA4) and apoptosis-inducing factor, Pdr lacks one of the arginine residues that compensates partially for the negative charge on the pyrophosphate of FAD. This uncompensated negative charge is likely to decrease the electron-accepting ability of the flavin. The aromatic side-chain of the "gatekeeper" Tyr159 is in the "out" conformation and leaves the nicotinamide-binding site of Pdr completely open. The presence of electron density in the NAD-binding channel indicates that NAD originating from Escherichia coli is partially bound to Pdr. A structural comparison of Pdr with homologous flavoproteins indicates that an open and accessible nicotinamide-binding site, the presence of an acidic residue in the middle part of the NAD-binding channel that binds the nicotinamide ribose, and multiple positively charged arginine residues surrounding the entrance of the NAD-binding channel are the special structural elements that assist tighter and more specific binding of the oxidized pyridine nucleotide by the BphA4-like flavoproteins. The crystallographic model of Pdr explains differences in the electron transfer mechanism in the Pdr-putidaredoxin redox couple and their mammalian counterparts, adrenodoxin reductase and adrenodoxin.     

Structural evidence for a functionally relevant second camphor binding site in P450cam: model for substrate entry into a P450 active site

P450cam has long served as a prototype for the cytochrome P450 (CYP) gene family. But, little is known about how substrate enters its active site pocket, and how access is achieved in a way that minimizes exposure of the reactive heme. We hypothesize that P450cam may first bind substrate transiently near the mobile F-G helix that covers the active site pocket. Such a two-step binding process is kinetically required if P450cam rarely populates an open conformation-as suggested by previous literature and the inability to obtain a crystal structure of P450cam in an open conformation. Such a mechanism would minimize exposure of the heme by allowing P450cam to stay in a closed conformation as long as possible, since only brief flexing into an open conformation would be required to allow substrate entry. To test this model, we have attempted to dock a second camphor molecule into the crystal structure of camphor-bound P450cam. The docking identified only one potential entry site pocket, a well-defined cavity on the F-helix side of the F-G flap, 16 A from the heme iron. Location of this entry site pocket is consistent with our NMR T1 relaxation-based measurements of distances for a camphor that binds in fast exchange (active site camphor is known to bind in slow exchange). Presence of a second camphor binding site is also confirmed with [(1)H-(13)C] HSQC titrations of (13)CH3-threonine labeled P450cam. To confirm that camphor can bind outside of the active site pocket, (13)CH3-S-pyridine was bound to the heme iron to physically block the active site, and to serve as an NMR chemical shift probe. Titration of this P450cam-pyridine complex confirms that camphor can bind to a site outside the active site pocket, with an estimated Kd of 43 microM. The two-site binding model that is proposed based on these data is analogous to that recently proposed for CYP3A4, and is consistent with recent crystal structures of P450cam bound to tethered-substrates, which force a partially opened conformation.     

Molecular dynamics simulations of phenylimidazole inhibitor complexes of cytochrome P450 cam

Molecular dynamics simulations have been performed on three phenylimidazole inhibitor complexes ofP450 _cam, utilizing the X-ray structures and the AMBER suite of programs. Compared to their corresponding optimized X-ray structures, very similar features were observed for the 1-phenylimidazole (1-PI) and 2-phenylimidazole (2-PI) complexes during a 100 ps MD simulation. The 1-PI inhibitor binds as a Type II complex with the imidazole nitrogen as a ligand of the heme iron. Analysis of the inhibitor-enzyme interctions during the MD simulations reveals that electrostatic interactions of the imidazole with the heme and van der Waals interactions of the phenyl ring with nearby hydrophobic residues are dominant. By contrast, 2-PI binds as a Type I inhibitor in the substrate binding pocket, but not as a ligand of the iron. The interactions of this inhibitor are qualitatively different from that of the Type II 1-PI, being mainly electrostatic/H-bonding interactions with a bound water and polar residues. Although the third compound, 4-PI, in common with 1-PI, also binds as a Type II inhibitor, with one nitrogen of the imidazole as a ligand to the iron, the MD average binding orientation deviates significantly from the X-ray structure. The most important changes observed include: (1) the rotation of the imidazole ring of this inhibitor by about 90° to enhance electrostatic interactions of the imidazole NH group with the carbonyl group of LEU244, and (2) the rotation of the carbonyl group of ASP251 to form a H-bond with VAL254. An analysis of the H-bonding network surrounding this substrate in the optimized crystal structure revealed that there is no H-bonding partner either for the free polar NH group in the imidazole ring of 4-phenylimidazole or for the polar carbonyl group of the nearby ASP251 residue. The deviation of the dynamically averaged inhibitor-enzyme structure of the 4-PI complex from the optimized crystal structure can therefore be rationalized as a consequence of the optimization of the electrostatic interactions among the polar groups.     

Molecular dynamics of the P450cam–Pdx complex reveals complex stability and novel interface contacts

Cytochrome P450cam catalyzes the stereo and regiospecific hydroxylation of camphor to 5‐exo‐hydroxylcamphor. The two electrons for the oxidation of camphor are provided by putidaredoxin (Pdx), a Fe₂S₂ containing protein. Two recent crystal structures of the P450cam–Pdx complex, one solved with the aid of covalent cross‐linking and one without, have provided a structural picture of the redox partner interaction. To study the stability of the complex structure and the minor differences between the recent crystal structures, a 100 nanosecond molecular dynamics (MD) simulation of the cross‐linked structure, mutated in silico to wild type and the linker molecule removed, was performed. The complex was stable over the course of the simulation though conformational changes including the movement of the C helix of P450cam further toward Pdx allowed for the formation of a number of new contacts at the complex interface that remained stable throughout the simulation. While several minor crystal contacts were lost in the simulation, all major contacts that had been experimentally studied previously were maintained. The equilibrated MD structure contained a mixture of contacts resembling both the cross‐linked and noncovalent structures and the newly identified interactions. Finally, the reformation of the P450cam Asp251–Arg186 ion pair in the MD simulation mirrors the ion pair observed in the more promiscuous CYP101D1 and suggests that the Asp251–Arg186 ion pair may be important.     

Exploring the structure characteristics and major channels of cytochrome P450 2A6, 2A13, and 2E1 with pilocarpine

The majority of cytochromes P450 play a critical role in metabolism of endogenous and exogenous substrates, some of its products are carcinogens. Therefore, inhibition of P450 enzymes activity can promote the detoxification and elimination of chemical carcinogens. In this study, molecular dynamics (MD) simulations and adaptive steered molecular dynamics (ASMD) simulations were performed to explore the structure features and channel dynamics of three P450 isoforms 2A6, 2A13, and 2E1 bound with the common inhibitor pilocarpine. The binding free energy results combined with the PMF calculations give a reasonable ranking of binding affinity, which are consistent with the experimental data. Our results uncover how a sequence divergence of different CYP2 enzymes causes individual variations in major channel selections. On the basis of channel bottleneck and energy decomposition analysis, we propose a gating mechanism of their respective major channels in three enzymes, which may be attributed to a reversal of Phe209 in CYP2A6/2A13, as well as the rotation of Phe116 and Phe298 in CYP2E1. The hydrophobic residues not only make strong hydrophobic interactions with inhibitor, but also act as gatekeeper to regulate the opening of channel. The present study provides important insights into the structure–function relationships of three cytochrome P450s and the molecular basis for development of potent inhibitors.     

Exploring coumarin egress channels in human cytochrome P450 2A6 by random acceleration and steered molecular dynamics simulations

The kinetic analysis of coumarin oxidation by CYP2A6 suggested that substrate binding and release occurred in the multiple steps and such events proceeded rapidly. However, the crystal structure of the CYP2A6-coumarin complex reveals that no obvious channel is open enough to allow coumarin to pass through. Thus, an intriguing and important question arises: how coumarin enters and exits the active site, which is deeply buried at the center of CYP2A6 fold. In this study, geometric analysis of the potential openings was first performed on all the available crystal structures of CYP2A6. And then, random acceleration molecular dynamics simulations were used to explore the possible substrate egress channels in CYP2A6. Two channels were most frequently observed. Afterwards, steered molecular dynamics simulations were performed and potentials of mean force were constructed to compare the preference of the two channels serving as the substrate egress channel. The results showed that channel 2c, which is located between helices I and G and the helix B'-C region, was the most likely channel for coumarin egress. The opening of channel 2c was characterized by a rotation of Phe111 together with a bending of helix B'. Our findings will not only be helpful for understanding the unbinding mechanism of coumarin and for identifying structural determinants related to the biological function of CYP2A6, but also provide further insight into the channel selectivity of P450s.     

Continuous B-epitope maps of cytochrome P450cam (CYP101) obtained by peptide scanning: correlation to spatial structure

Protein continuous B-epitopes can be revealed using short synthetic peptides that overlap a known protein sequence. Since the whole protein surface is considered to possess antigenic properties, a question that arises is whether a set of linear B-epitopes determined by peptide scanning correlates with a protein spatial structure. We have chosen cytochrome P450cam (CYP101) of Pseudomonas putida, with known 3D structure, as a template. Sera of two rabbits and antibody egg yolk preparations from three chickens were produced against the P450cam molecule. These polyclonals were analyzed separately in ELISA with 409 overlapping P450cam hexapeptides. The whole set of continuous antigenic sites of P450cam covered about 45% of the P450cam sequence. However, immunodominant sites (those revealed with more than 50% antibody preparations), the so-called "antigenic core," represent only 9% of the protein sequence. While the amount of water-accessible residues in the total antigenic map (42%) was close to that in the whole native P450cam molecule (39%), the amount of water-accessible residues in the antigenic core was significantly higher (64%). These results led to the conclusion that antigenic core epitopes can be associated to the molecular surface, whereas epitopes with low detection frequency may partly correspond to unfolded regions of the protein molecule.     

Crystal structure of putidaredoxin, the [2Fe-2S] component of the P450cam monooxygenase system from Pseudomonas putida

Stability of the [2Fe-2S]-containing putidaredoxin (Pdx), the electron donor to cytochrome P450cam in Pseudomonas putida, was improved by mutating non-ligating cysteine residues, Cys73 and Cys85, to serine singly and in combination. The increasing order of stability is Cys73Ser/Cys85Ser>Cys73Ser>Cys85Ser>WT Pdx. Crystal structures of Cys73Ser/Cys85Ser and Cys73Ser mutants of Pdx, solved by single-wavelength anomalous dispersion phasing using the [2Fe-2S] iron atoms to 1.47 A and 1.65 A resolution, respectively, are nearly identical and very similar to those of bovine adrenodoxin (Adx) and Escherichia coli ferredoxin. However, unlike the Adx structure, no motion between the core and interaction domains of Pdx is observed. This higher conformational stability of Pdx might be due to the presence of a more extensive hydrogen bonding network at the interface between the two structural domains around the conserved His49. In particular, formation of a hydrogen bond between the side-chain of Tyr51 and the carbonyl oxygen atom of Glu77 and the presence of two well-ordered water molecules linking the interaction domain and the C-terminal peptide to the core of the molecule are unique to Pdx. The folding topology of the NMR model is similar to that of the X-ray structure of Pdx. The overall rmsd of Calpha positions between the two models is 1.59 A. The largest positional differences are observed for residues 18-21 and 33-37 in the loop regions and the C terminus. The latter two peptides display conformational heterogeneity in the crystal structures. Owing to flexibility, the aromatic ring of the C-terminal Trp106 can closely approach the side-chains of Asp38 and Thr47 (3.2-3.9 A) or move away and leave the active site solvent-exposed. Therefore, Trp106, previously shown to be important in the Pdr-to-Pdx and Pdx-to-P450cam electron transfer reactions is in a position to regulate and/or mediate electron transfer to or from the [2Fe-2S] center of Pdx.     

Solution NMR structure of putidaredoxin-cytochrome P450cam complex via a combined residual dipolar coupling-spin labeling approach suggests a role for Trp106 of putidaredoxin in complex formation

The 58-kDa complex formed between the [2Fe-2S] ferredoxin, putidaredoxin (Pdx), and cytochrome P450cam (CYP101) from the bacterium Pseudomonas putida has been investigated by high-resolution solution NMR spectroscopy. Pdx serves as both the physiological reductant and effector for CYP101 in the enzymatic reaction involving conversion of substrate camphor to 5-exo-hydroxycamphor. In order to obtain an experimental structure for the oxidized Pdx-CYP101 complex, a combined approach using orientational data on the two proteins derived from residual dipolar couplings and distance restraints from site-specific spin labeling of Pdx has been applied. Spectral changes for residues in and near the paramagnetic metal cluster region of Pdx in complex with CYP101 have also been mapped for the first time using ¹⁵N and ¹³C NMR spectroscopy, leading to direct identification of the residues strongly affected by CYP101 binding. The new NMR structure of the Pdx-CYP101 complex agrees well with results from previous mutagenesis and biophysical studies involving residues at the binding interface such as formation of a salt bridge between Asp38 of Pdx and Arg112 of CYP101, while at the same time identifying key features different from those of earlier modeling studies. Analysis of the binding interface of the complex reveals that the side chain of Trp106, the C-terminal residue of Pdx and critical for binding to CYP101, is located across from the heme-binding loop of CYP101 and forms non-polar contacts with several residues in the vicinity of the heme group on CYP101, pointing to a potentially important role in complex formation.     

Establishment of a yeast system that stably expresses human cytochrome P450 reductase: application for the study of drug metabolism of cytochrome P450s in vitro

Cytochrome P450s (CYPs) hold a balance in studying pharmacokinetics, toxico-kinetics, drug metabolism, and drug-drug interactions, which require association with cytochrome P450 reductase (CPR) to achieve optimal activity. A novel system of Saccharomyces cerevisiae useful for expression studies of mammalian microsomal CYPs was established. Human CPR (hCPR) was co-expressed with human CYP3A4 (hCYP3A4) in this system, and two expression plasmids pTpLC and pYeplac195-3A4 containing the cDNA of hCPR and hCYP3A4 were constructed, respectively. The two plasmids were applied first and controlled by phosphoglycerate kinase (PGK) promoter. S. cerevisiae BWG1-7alpha transformed with the expression plasmids produced the respective proteins in the expected molecular sizes reactive with both anti-hCYP3A4 immunoglobulin (Ig) and anti-hCPR Ig. The activity of hCPR in yeast BWG-CPR was 443.2 nmol reduced cytochrome c/min/mg, which was about three times the CPR activity of the microsome prepared from the parental yeast. The protein amount of hCYP3A4 in BWG-CPR/3A4 was 35.53 pmol/mg, and the 6beta-hydroxylation testosterone formation activity of hCYP3A4 expressed was 7.5 nmol/min/nmol CYP, 30 times higher than the activity of hCYP3A4 expressed in the parental yeast, and almost two times the activity of hCYP3A4 from homologous human liver microsome. Meanwhile, BWG-CPR/3A4 retained 100 generations under nonselective culture conditions, indicating this yeast was a mitotically stable transformant. BWG-CPR was further tested daily by the PCR amplification of hCPR of yeast genome, Western blot analysis, and the activity assay of hCPR of yeast microsome. This special expression host for CYPs was validated to be stable and efficient for the expression of CYPs, applying as an effective selection model for the drug metabolism in vitro.