Metal-substituted protein MRI contrast agents engineered for enhanced relaxivity and ligand sensitivity

Engineered metalloproteins constitute a flexible new class of analyte-sensitive molecular imaging agents detectable by magnetic resonance imaging (MRI), but their contrast effects are generally weaker than synthetic agents. To augment the proton relaxivity of agents derived from the heme domain of cytochrome P450 BM3 (BM3h), we formed manganese(III)-containing proteins that have higher electron spin than their native ferric iron counterparts. Metal substitution was achieved by coexpressing BM3h variants with the bacterial heme transporter ChuA in Escherichia coli and supplementing the growth medium with Mn3+-protoporphyrin IX. Manganic BM3h variants exhibited up to 2.6-fold higher T1 relaxivities relative to native BM3h at 4.7 T. Application of ChuA-mediated porphyrin substitution to a collection of thermostable chimeric P450 domains resulted in a stable, high-relaxivity BM3h derivative displaying a 63% relaxivity change upon binding of arachidonic acid, a natural ligand for the P450 enzyme and an important component of biological signaling pathways. This work demonstrates that protein-based MRI sensors with robust ligand sensitivity may be created with ease by including metal substitution among the toolkit of methods available to the protein engineer.     

Identification of mutant Asp251Gly/Gln307His of cytochrome P450 BM3 for the generation of metabolites of diclofenac, ibuprofen and tolbutamide

The soluble, catalytically self-sufficient cytochrome P450 BM3 from Bacillus megaterium is a good candidate as biocatalyst for the synthesis of drug metabolites. To this end, error-prone polymerase chain reaction (PCR) was used to generate a library of P450 BM3 mutants with novel activities toward drugs. The double mutant Asp251Gly/Gln307His (A2) with activities towards diclofenac, ibuprofen and tolbutamide was identified by screening with the alkali method. This is based on the detection of NADPH oxidation during enzymatic turnover on whole Escherichia coli cells heterologously expressing the P450 BM3 mutants in the presence of the target substrates. The three drugs screened are marker substrates of human liver cytochromes P450 belonging to the 2C subfamily. Interestingly the mutations Asp251Gly/Gln307His are located on the protein surface and they are not directly involved in substrate binding and turnover. Dissociation constants and K(M) values of mutant A2 for diclofenac, ibuprofen and tolbutamide are in the micromolar range. Catalysis leads to hydroxylations in specific positions, producing 4'-hydroxydiclofenac, 2-hydroxyibuprofen and 4-hydroxytolbutamide, respectively.     

Light-initiated hydroxylation of lauric acid using hybrid P450 BM3 enzymes

We have developed hybrid P450 BM3 enzymes consisting of a Ru(II)-diimine photosensitizer covalently attached to non-native single cysteine residues of P450 BM3 heme domain mutants. These enzymes are capable, upon light activation, of selectively hydroxylating lauric acid with 40 times higher total turnover numbers compared to the peroxide shunt.     

Resonance surface plasmon spectroscopy: low molecular weight substrate binding to cytochrome p450

A new detection mechanism has been developed for low molecular weight substrate binding to heme proteins based on resonance localized surface plasmon spectroscopy. Cytochrome P450 has strong electronic transitions in the visible wavelength region. Upon binding of a substrate molecule (e.g., camphor), the absorption band of cytochrome P450 shifts to shorter wavelength. The event of camphor binding to a nanoparticle surface modified with cytochrome P450 protein receptors is monitored using UV-vis spectroscopy. It is observed for the first time that the binding of the substrate molecules to the protein receptor induces a blue-shift in the localized surface plasmon resonance (LSPR) of the nanosensors. The coupling between the molecular resonance of the substrate-free and substrate-bound cytochrome P450 proteins and the nanoparticles' LSPR leads to a highly wavelength-dependent LSPR response. When the LSPR of the nanoparticles is located at a wavelength distant from the cytochrome P450 resonance, an average of approximately 19 nm red-shift is observed upon cytochrome P450 binding to the nanoparticles and a approximately 6 nm blue-shift is observed upon camphor binding However, this response is significantly amplified approximately 3 to 5 times when the LSPR of the nanoparticles is located at a slightly longer wavelength than the cytochrome P450 resonance, that is, a 66.2 nm red-shift upon cytochrome P450 binding and a 34.7 nm blue-shift upon camphor binding. This is the first example of the detection of small molecules binding to a protein modified nanoparticle surface on the basis of LSPR.     

Electrochemistry of cytochrome P450 BM3 in sodium dodecyl sulfate films

Direct electrochemistry of the cytochrome P450 BM3 heme domain (BM3) was achieved by confining the protein within sodium dodecyl sulfate (SDS) films on the surface of basal-plane graphite (BPG) electrodes. Cyclic voltammetry revealed the heme FeIII/II redox couple at -330 mV (vs Ag/AgCl, pH 7.4). Up to 10 V/s, the peak current was linear with the scan rate, allowing us to treat the system as surface-confined within this regime. The standard heterogeneous rate constant determined at 10 V/s was estimated to be 10 s-1. Voltammograms obtained for the BM3-SDS-BPG system in the presence of dioxygen exhibited catalytic waves at the onset of FeIII reduction. The altered heme reduction potential of the BM3-SDS-graphite system indicates that SDS is likely bound in the enzyme active-site region. Compared to other P450-surfactant systems, we find redox potentials and electron-transfer rates that differ by approximately 100 mV and >10-fold, respectively, indicating that the nature of the surfactant environment has a significant effect on the observed heme redox properties.     

Design of a five-coordinate heme protein maquette: a spectroscopic model of deoxymyoglobin

The substitution of 1-methyl-l-histidine for the histidine heme ligands in a de novo designed four-alpha-helix bundle scaffold results in conversion of a six-coordinate cytochrome maquette into a self-assembled five-coordinate mono-(1-methyl-histidine)-ligated heme as an initial maquette for the dioxygen carrier protein myoglobin. UV-vis, magnetic circular dichroism, and resonance Raman spectroscopies demonstrate the presence of five-coordinate mono-(1-methyl-histidine) ligated ferrous heme spectroscopically similar to deoxymyoglobin. Thermodynamic analysis of the ferric and ferrous heme dissociation constants indicates greater destabilization of the ferric state than the ferrous state. The ferrous heme protein reacts with carbon monoxide to form a (1-methyl-histidine)-Fe(II)(heme)-CO complex; however, reaction with dioxygen leads to autoxidation and ferric heme dissociation. These results indicate that negative protein design can be used to generate a five-coordinate heme within a maquette scaffold.     

Laser photoexcitation of NAD(P)H induces reduction of P450 BM3 heme domain on the microsecond time scale

We demonstrate that photoexcitation of NAD(P)H at 355 nm using a Nd:YAG laser leads to rapid reduction of the heme domain of the Bacillus megaterium fatty acid hydroxylase flavocytochrome P450 BM3. An aqueous electron derived from photoexcited NAD(P)H is rapidly transferred to the heme domain, enabling the formation of a carbon monoxy complex of the ferrous P450 (FeII-CO) on the microsecond time scale. Using this approach we have determined the limiting rate constant (1770 s-1 for substrate-free heme domain) for formation of the FeII-CO complex. We find no dependence of the observed rate of FeII-CO complex formation on NAD(P)H concentration but demonstrate a hyperbolic dependence on carbon monoxide concentration. The apparent dissociation constant for the complex of carbon monoxide bound noncovalently to the ferric form of the BM3 heme domain (and with NADH as reductant) is 323 microM. Binding of a P450 substrate (N-palmitoylglycine) weakened the complex between carbon monoxide and the ferric BM3 heme domain (Kd increased to 1404 microM) but enhanced the rate of formation of the FeII-CO complex (3036 s-1 for substrate-free heme domain). This study demonstrates the applicability of NAD(P)H photoexcitation as a method for rapid electron delivery to P450 enzymes and provides a new route to probing the P450 catalytic cycle and its transient intermediates.     

Regioselective Hydroxylation of C12–C15 Fatty Acids with Fluorinated Substituents by Cytochrome P450 BM3

We demonstrate herein that wild‐type cytochrome P450 BM3 can recognize non‐natural substrates, such as fluorinated C₁₂–C₁₅ chain‐length fatty acids, and show better catalysis for their efficient conversion. Although the binding affinities for fluorinated substrates in the P450 BM3 pocket are marginally lower than those for non‐fluorinated substrates, spin‐shift measurements suggest that fluoro substituents at the ω‐position can facilitate rearrangement of the dynamic structure of the bulk‐water network within the hydrophobic pocket through a micro desolvation process to expel the water ligand of the heme iron that is present in the resting state. A lowering of the Michaelis–Menten constant (K_m), however, indicates that fluorinated fatty acids are indeed better substrates compared with their non‐fluorinated counterparts. An enhancement of the turnover frequencies (k_cat) for electron transfer from NADPH to the heme iron and for C? H bond oxidation by compound I (Cpd I) to yield the product suggests that the activation energies associated with going from the enzyme–substrate (ES state) to the corresponding transition state (ES^≠ state) are significantly lowered for both steps in the case of the fluorinated substrates. Delicate control of the regioselectivity by the fluorinated terminal methyl groups of the C₁₂–C₁₅ fatty acids has been noted. Despite the fact that residues Arg47/Tyr51/Ser72 exert significant control over the hydroxylation of the subterminal carbon atoms toward the hydrocarbon tail, the fluorine substituent(s) at the ω‐position affects the regioselective hydroxylation. For substrate hydroxylation, we have found that fluorinated lauric acids probably give a better structural fit for the heme pocket than fluorinated pentadecanoic acid, even though pentadecanoic acid is by far the best substrate among the reported fatty acids. Interestingly, 12‐fluorododecanoic acid, with only one fluorine atom at the terminal methyl group, exhibits a comparable turnover frequency to that of pentadecanoic acid. Thus, fluorination of the terminal methyl group introduces additional interactions of the substrate within the hydrophobic pocket, which influence the electron transfers for both dioxygen activation and the controlled oxidation of aliphatics mediated by high‐valent oxoferryl species.     

Are branched chain fatty acids the natural substrates for P450(BM3)?

Branched chain fatty acids are substrates for cytochrome P450(BM3) (CYP102) from Bacillus megaterium; oxidation of C15 and C17 iso and anteiso fatty acids by P450(BM3) leads to the formation of hydroxylated products that possess high levels of regiochemical and stereochemical purity.     

Spectra and kinetic studies of the compound I derivative of cytochrome P450 119

The Compound I derivative of cytochrome P450 119 (CYP119) was produced by laser flash photolysis of the corresponding Compound II derivative, which was first prepared by reaction of the resting enzyme with peroxynitrite. The UV-vis spectrum of the Compound I species contained an asymmetric Soret band that could be resolved into overlapping transitions centered at approximately 367 and approximately 416 nm and a Q band with lambda(max) approximately 650 nm. Reactions of the Compound I derivative with organic substrates gave epoxidized (alkene oxidation) and hydroxylated (C-H oxidation) products, as demonstrated by product studies and oxygen-18 labeling studies. The kinetics of oxidations by CYP119 Compound I were measured directly; the reactions included hydroxylations of benzyl alcohol, ethylbenzene, Tris buffer, lauric acid, and methyl laurate and epoxidations of styrene and 10-undecenoic acid. Apparent second-order rate constants, equal to the product of the equilibrium binding constant (K(bind)) and the first-order oxidation rate constant (k(ox)), were obtained for all of the substrates. The oxidations of lauric acid and methyl laurate displayed saturation kinetic behavior, which permitted the determination of both K(bind) and k(ox) for these substrates. The unactivated C-H positions of lauric acid reacted with a rate constant of k(ox) = 0.8 s(-1) at room temperature. The CYP119 Compound I derivative is more reactive than model Compound I species [iron(IV)-oxo porphyrin radical cations] and similar in reactivity to the Compound I derivative of the heme-thiolate enzyme chloroperoxidase. Kinetic isotope effects (kH/kD) for oxidations of benzyl alcohol and ethylbenzene were small, reflecting the increased reactivity of the Compound I derivative in comparison to models. Nonetheless, CYP119 Compound I apparently is much less reactive than the oxidizing species formed in the P450 cam reaction cycle. Studies of competition kinetics employing CYP119 activated by hydrogen peroxide indicated that the same oxidizing transient is formed in the photochemical reaction and in the hydrogen peroxide shunt reaction.     

Enhanced electron transfer and lauric acid hydroxylation by site-directed mutagenesis of CYP119

CYP119, a cytochrome P450 from a thermophilic organism for which a crystal structure is available, is shown here to hydroxylate lauric acid in a reaction supported by putidaredoxin and putidaredoxin reductase. This fatty acid hydroxylation activity is increased 15-fold by T214V and D77R mutations. The T214V mutation increases the rate by facilitating substrate binding and enhancing the associated spin state change, whereas the D77R mutation improves binding of the heterologous redox partner putidaredoxin to CYP119 and the rate of electron transfer from it to the heme group. A sequence alignment with P450(cam) can, therefore, be used to identify a part of the binding site for putidaredoxin on an unrelated P450 enzyme. This information can be used to engineer by mutagenesis an improved complementarity of the protein-protein interface that results in improved electron transfer from putidaredoxin to the P450 enzyme. As a result, the catalytic activity of the thermo- and barostable CYP119 has been incorporated into a catalytic system that hydroxylates fatty acids.     

Oxygen activation and electron transfer in flavocytochrome P450 BM3

In flavocytochrome P450 BM3, there is a conserved phenylalanine residue at position 393 (Phe393), close to Cys400, the thiolate ligand to the heme. Substitution of Phe393 by Ala, His, Tyr, and Trp has allowed us to modulate the reduction potential of the heme, while retaining the structural integrity of the enzyme's active site. Substrate binding triggers electron transfer in P450 BM3 by inducing a shift from a low- to high-spin ferric heme and a 140 mV increase in the heme reduction potential. Kinetic analysis of the mutants indicated that the spin-state shift alone accelerates the rate of heme reduction (the rate determining step for overall catalysis) by 200-fold and that the concomitant shift in reduction potential is only responsible for a modest 2-fold rate enhancement. The second step in the P450 catalytic cycle involves binding of dioxygen to the ferrous heme. The stabilities of the oxy-ferrous complexes in the mutant enzymes were also analyzed using stopped-flow kinetics. These were found to be surprisingly stable, decaying to superoxide and ferric heme at rates of 0.01-0.5 s(-)(1). The stability of the oxy-ferrous complexes was greater for mutants with higher reduction potentials, which had lower catalytic turnover rates but faster heme reduction rates. The catalytic rate-determining step of these enzymes can no longer be the initial heme reduction event but is likely to be either reduction of the stabilized oxy-ferrous complex, i.e., the second flavin to heme electron transfer or a subsequent protonation event. Modulating the reduction potential of P450 BM3 appears to tune the two steps in opposite directions; the potential of the wild-type enzyme appears to be optimized to maximize the overall rate of turnover. The dependence of the visible absorption spectrum of the oxy-ferrous complex on the heme reduction potential is also discussed.     

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.     

Reaction of cytochrome P450BM3 and peroxynitrite yields nitrosyl complex

Peroxynitrite has come into the spotlight in recent years. Its effects on proteins have been implicated in several diseases such as acute lung injury, rheumatoid arthritis, implant rejection, artherosclerosis, Parkinson's disease, and Alzheimer's disease. Peroxynitrite is thought to inactivate a variety of proteins including thiolate-ligated heme proteins such as cytochrome P450 2B1 and PGI2 synthase, through the nitration of tyrosine residues. In previous studies it was reported that thiolate-ligated heme enzymes react with peroxynitrite to form a ferryl intermediate. In an effort to spectroscopically characterize this species in P450BM3, we discovered that the peroxynitrite-generated intermediate is not an FeIVoxo, but rather an iron-nitrosyl [FeNO]6 complex. We present density functional calculations as well as M?ssbauer and stopped-flow spectroscopic characterizations of the peroxynitrite-generated intermediate in P450BM3.     

Whole‐Cell Biotransformation of Benzene to Phenol Catalysed by Intracellular Cytochrome P450BM3 Activated by External Additives

An Escherichia coli whole‐cell biocatalyst for the direct hydroxylation of benzene to phenol has been developed. By adding amino acid derivatives as decoy molecules to the culture medium, wild‐type cytochrome P450BM3 (P450BM3) expressed in E.coli can be activated and non‐native substrates hydroxylated, without supplementing with NADPH. The yield of phenol reached 59 % when N‐heptyl‐l ‐prolyl‐l ‐phenylalanine (C7‐Pro‐Phe) was employed as the decoy molecule. It was shown that decoy molecules, especially those lacking fluorination, reached the cytosol of E. coli, thus imparting in vivo catalytic activity for the oxyfunctionalisation of non‐native substrates to intracellular P450BM3.     

Reduction of dioxygen catalyzed by pyrene-wired heme domain cytochrome P450 BM3 electrodes

We have electronically wired the cytochrome P450 BM3 heme domain to a graphite electrode with the use of a pyrene-terminated tether. AFM images clearly reveal that pyrene-wired enzyme molecules are adsorbed to the electrode surface. The enzyme-electrode system undergoes rapid and reversible electron transfer, displaying a standard rate constant higher than that of any other P450-electrode system. We also show that the graphite-pyrene-BM3 system catalyzes the four-electron reduction of dioxygen to water.     

Free Energy Calculations Give Insight into the Stereoselective Hydroxylation of α-Ionones by Engineered Cytochrome P450 BM3 Mutants

Previously, stereoselective hydroxylation of α-ionone by Cytochrome P450 BM3 mutants M01 A82W and M11 L437N was observed. While both mutants hydroxylate α-ionone in a regioselective manner at the C3 position, M01 A82W catalyzes formation of trans-3-OH-α-ionone products whereas M11 L437N exhibits opposite stereoselectivity, producing trans-(3S,6S)-OH-α-ionone and cis-(3S,6R)-OH-α-ionone. Here, we explore the stereoselective C3 hydroxylation of α-ionone by Cytochrome P450 BM3 mutants M01 A82W and M11 L437N using molecular dynamics-based free energy calculations to study the interaction between the enzyme and both the substrates and the products. The one-step perturbation approach is applied using an optimized reference state for substrates and products. While the free energy differences between the substrates free in solution amount to ～0 kJ mol(-1), the differences in mutant M01 A82W agree with the experimentally obtained dissociation constants K(d). Moreover, a correlation with experimentally observed trends in product formation is found in both mutants. The trans isomers show the most favorable relative binding free energy in the range of all four possible hydroxylated diastereomers for mutant M01 A82W, while the trans product from (6S)-α-ionone and the cis product from (6R)-α-ionone show highest affinity for mutant M11 L437N. Marcus theory is subsequently used to relate the thermodynamic stability to transition state energies and rates of formation.     

Oxygen activation by cytochrome p450: a thermodynamic analysis

Electrode potentials for every intermediate in the cytochrome P450 cycle were estimated and evaluated by means of an oxidation state diagram. By this approach, and within the uncertainties of the approximations, the superoxide complex of cytochrome P450 at pH 7 is oxidizing: E degrees ' (P450FeO(2)2+, H+/P450FeOOH2+) = +0.93 V, and the Gibbs energy for the reaction of the hydroperoxo complex of cytochrome P450 to form compound I and water, P450FeOOH2+ + H+ = P450FeO2+ por(*+) + H2O, is 0 kJ/mol. Although cytochrome P450FeOOH2+ and cytochrome P450FeO2+ por(*+) are approximately isoenergetic, they are likely to react at different rates with substrates and may yield different products. Homolysis of the hydroperoxo complex of cytochrome P450 to compound II and the hydroxyl radical, P450FeOOH2+ = P450FeO2+ + HO(*), is unfavorable (DeltaG degrees ' = +92 kJ/mol), as is the dissociation into HOO- and cytochrome P450Fe3+ (+73 kJ/mol). It is shown that the sum of the Gibbs energy of association for cytochrome P450Fe3+ with the hydroperoxo anion and the Gibbs energy for the one-electron reduction of cytochrome P450FeOOH2+, relative to NHE, is constant (-203 kJ/mol). While the estimated E degrees ' (P450FeO(2)2+, H+/P450FeOOH2+) = +0.93 V at pH 7 is larger than necessary to effect reduction of cytochrome P450FeO(2)2+, the magnitude of this electrode potential implies that the binding constant for cytochrome P450Fe3+ with hydrogen peroxide is ca. 3 x 106 M(-1) at pH 7. An association constant of this magnitude ensures that a fraction of cytochrome P450FeOOH2+ is available to form compound I or to react with substrates directly, while a larger one would imply that compound I is too weak an oxidant. In general, the energetics of the reduction of dioxygen to water determines the energetics of catalysis of hydroxylations by cytochrome P450. These results enable calibration of energy levels obtained for intermediates in the cytochrome P450 reaction cycle obtained by ab initio calculations and provide insights into the catalytic efficiency of cytochrome P450 and guidelines for the development of competent hydroxylation catalysts.     

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.     

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.