Monooxygenation mechanism by cytochrome p-450

Regioselective oxidation of heterocyclic aromatic hydrocarbons catalyzed by cytochrome P450: A case study of carbazole

Recently there are increasing interests in accurately evaluating the health effects of heterocyclic PAHs. However, the activation mechanism and possible metabolites of heterocyclic PAHs catalyzed by human CYP1A1 is still elusive to a great extent. Here, leveraged to high level QM/MM calculations, the corresponding activation pathways of a representative heterocyclic PAHs, carbazole, were systematically explored. The first stage is electrophilic addition or hydrogen abstraction from N-H group. Electrophilic addition was evidenced to be more feasible and regioselectivity at C3 and C4 sites were identified. Correlations between energy barriers and key structural/electrostatic parameters reveal that O-Cα distance and Fe-O-Cα angle are the main origin for the catalytic regioselectivity. Electrophilic addition was determined as the rate-determining step and the subsequent possible reactions include epoxidation, NIH shift (the hydrogen migration from the site of hydroxylation to the adjacent carbon) and proton shuttle. The corresponding products are epoxides, ketones and hydroxylated carbazoles, respectively. The main metabolites (hydroxylated carbazoles) are estimated to be more toxic than carbazole. The regioselectivity of carbazole activated by CYP1A1 is different from the environmental processes (gas and aqueous phase). Collectively, these results will inform the in-depth understanding the metabolic processes of heterocyclic PAHs and aid the accurate evaluation of their health effects.

Synthesis and biological evaluation of coumarin derivatives as selective CYP2A6 inhibitors

Cytochrome P450 2A6 (CYP2A6) inhibitors are expected to be suitable as smoking cessation aids and for cancer prevention. Because the typical coumarin-based CYP2A6 inhibitor methoxsalen also inhibits CYP3A4, unintended drug-drug interactions are still a concern. Therefore, the development of selective CYP2A6 inhibitors is desirable. In this study, we synthesized coumarin-based molecules, determined the IC₅₀ values for CYP2A6 inhibition, verified the possibility of mechanism-based inhibition, and compared the selectivity for CYP2A6 versus CYP3A4. The results demonstrated that we developed CYP2A6 inhibitors that were more potent and selective than methoxsalen.

Water biocatalytic effect attenuates cytochrome P450-mediated carcinogenicity of diethylnitrosamine: A computational insight

The mechanism-based mutagenicity and carcinogenicity of diethylnitrosamine (DEN) are believed to act through interactions with cytochrome P450 (P450) enzymes. DFT calculations to explore the conceivable mechanisms underlying the reaction of P450 with DEN with and without water as a biocatalyst were performed. The results shed light on the biocatalytic role of water in lowering the H-abstraction energy barriers because of the electrostatic effect driven by hydrogen bonding. Our DFT analysis revealed how metabolites are formed in the dealkylation (toxification) and denitrosation (detoxification) pathways. Also, our findings uncovered the active position of DEN vulnerable to P450 interactions. Two factors control the toxification and detoxification rates: the stability of denitrosation products and the HS rebound barrier of the α-pathway. Thus, water biocatalytic attenuation of DEN carcinogenicity was attained by stabilizing denitrosation products and slowing the α-HS rebound process. Docking and MD simulations were performed to assess the binding modes of DEN to P450's active site and to inspect the denitrosation and dealkylation processes, respectively.

High warfarin sensitivity in carriers of CYP2C9*35 is determined by the impaired interaction with P450 oxidoreductase

Cytochrome P450 2C9 (CYP2C9) metabolizes many clinically important drugs including warfarin and diclofenac. We have recently reported a new allelic variant, CYP2C9*35, found in a warfarin hypersensitive patient with Arg125Leu and Arg144Cys mutations. Here, we have investigated the molecular basis for the functional consequences of these polymorphic changes. CYP2C9.1 and CYP2C9-Arg144Cys expressed in human embryonic kidney 293 cells effectively metabolized both S-warfarin and diclofenac in NADPH-dependent reactions, whereas CYP2C9-Arg125Leu or CYP2C9.35 were catalytically silent. However, when NADPH was replaced by a direct electron donor to CYPs, cumene hydroperoxide, hereby bypassing the CYP oxidoreductase (POR), all variant enzymes were active, indicating unproductive interactions between CYP2C9.35 and POR. In silico analysis revealed a decrease of the electrostatic potential of CYP2C9-Arg125Leu-POR interacting surface and the loss of stabilizing salt bridges between these proteins. In conclusion, our data strongly suggest that the Arg125Leu mutation in CYP2C9.35 prevents CYP2C9-POR interactions resulting in the absence of NADPH-dependent CYP2C9-catalyzed activity in vivo, thus influencing the warfarin sensitivity in the carriers of this allele.

Ligand orientation of human neuroglobin obtained from solution NMR and molecular dynamics simulation as compared with X-ray crystallography

Neuroglobin, a new member of hemoprotein family, can reversibly bind oxygen and take part in many biological processes such as enzymatic reaction, signal transduction and the mitochondria function. Different from myoglobin and hemoglobin, it has a hexacoordinated heme environment, with histidyl imidazole of proximal His(96)(F8) and distal His(64)(E7) directly bound to the metal ion. In the present work, solution (1)H NMR spectroscopy was employed to investigate the electronic structure of heme center of wild-type met-human neuroglobin. The resonances of heme protons and key residues in the heme pocket were assigned. Two heme orientations resulting from a 180 degrees rotation about the alpha-gamma-meso axis with a population ratio about 2:1 were observed. Then the (1)H NMR chemical shifts of the ferriheme methyl groups were used to predict orientations of the axial ligand. The obtained axial ligand plane angle phi is consistent with that from the molecular dynamics simulation but not with those from the crystal data. Compared with mouse neuroglobin, the obtained average ligand orientation of human neuroglobin reflects the changeability of heme environment for the Ngb family.

S K-edge XAS and DFT calculations on cytochrome P450: covalent and ionic contributions to the cysteine-Fe bond and their contribution to reactivity

Experimental covalencies of the Fe-S bond for the resting low-spin and substrate-bound high-spin active site of cytochrome P450 are reported. DFT calculations on the active site indicate that one H-bonding interaction from the protein backbone is needed to reproduce the experimental values. The H-bonding to the thiolate from the backbone decreases the anisotropic pi covalency of the Fe-S bond lowering the barrier of free rotation of the exchangeable axial ligand, which is important for reactivity. The anionic axial thiolate ligand is calculated to lower the Fe(III/II) reduction potential of the active site by more than 1 V compared to a neutral imidazole ligand. About half of this derives from its covalent bonding and half from its electrostatic interaction with the oxidized Fe. This axial thiolate ligand increases the pK(a) of compound 0 (Fe(III)-hydroperoxo) favoring its protonation which promotes O-O bond heterolysis forming compound I. The reactivity of compound I is calculated to be relatively insensitive to the nature of the axial ligand due to opposing reduction potential and proton affinity contributions to the H-atom abstraction energy.

Theoretical calculation of triazolam hydroxylation and endogenous steroid inhibition in the active site of CYP3A4

CYP3A4 has unusual kinetic characteristics because it has a large active site. CYP3A4 produced more 4-hydroxytriazolam than alpha-hydroxytriazolam at concentrations of more than 60 muM triazolam, and different steroids had different inhibitory effects on the system. To clarify these interesting observations, the interactions between substrate and substrate/steroid were investigated by theoretical calculations. When two triazolam molecules were docked into the active site, the distance between the O-atom and the 4-hydroxylated site was less than the distance to the alpha-hydroxylated site because of interaction between the two triazolam molecules. Estradiol inhibited both alpha- and 4-hydroxytriazolam formation by 50%. Dehydroepiandrosterone (DHEA) inhibited alpha-hydroxylation more than 4-hydroxytriazolam formation, whereas aldosterone had no effect. When one triazolam molecule and one steroid molecule were simultaneously docked, estradiol increased the distance between the O-atom and the two hydroxylated sites, DHEA only increased the distance between the O-atom and alpha-hydroxylated site, and aldosterone did not change the distances. The relevant angles of Fe-O-C in the hydroxylated site of triazolam also widened, together with increased distance. These findings indicate that formation of a substrate and substrate/effector complex in the active site may be a factor for determining the enzyme kinetic parameters of CYP3A4.

Helices F-G are important for the substrate specificities of CYP3A7

CYP3A7 is a member of the human CYP3A family and a major form of P450 expressed in human fetal livers. Although CYP3A7 shares nearly 90% base sequence with CYP3A4, CYP3A7 shows striking functional differences in the catalytic preference for several substrates, such as dehydroepiandrosterone (DHEA) or dehydroepiandrosterone 3-sulfate (DHEA-3S). First, to clarify the reason for the differences between CYP3A7 and CYP3A4, a homology model of CYP3A7 was constructed using the CYP3A4 crystal structure. Because these two structures were similar, four kinds of chimeric enzymes were constructed to determine which sequences are important for exhibiting the characteristics of CYP3A7. The results of kinetic analysis of DHEA and DHEA-3S 16alpha-hydroxylations by CYP3A7, CYP3A4, and CYP3A chimeras suggested that the amino acid residues from Leu(210) to Glu(279) were important to express the specificity for substrates as CYP3A7. This region was on the F and G helices of the modeled CYP3A7. Furthermore, to assess which amino acid in this sequence is important for the substrate specificity of CYP3A7, a one-point mutation of CYP3A7 to CYP3A4 was made by site-directed mutagenesis. The mutants of K224T and K244E had lost DHEA and DHEA-3S 16alpha-hydroxylation activities. The mutants also greatly decreased the metabolism of testosterone, erythromycin, nevirapine, and triazolam relative to those activities of CYP3A7 wild-type enzyme. From these results, it is expected that CYP3A7 can recognize specific substrates using the lysines in F-G loops.

Theoretical investigations of the hydrolysis pathway of verdoheme to biliverdin

Conversion of iron(II) verdoheme to iron biliverdin in the presence of OH(-) was investigated using B3LYP method. Both 3-21G and 6-31G* basis sets were employed for geometry optimization calculation as well as energy stabilization estimation. Calculation at 6-31G* level was found necessary for a correct spin state estimation of the iron complexes. Two possible pathways for the conversion of iron verdoheme to iron biliverdin were considered. In one path the iron was six-coordinate while in the other it was considered to be five-coordinate. In the six-coordinated pathway, the ground state of bis imidazole iron verdoheme is singlet while that for open chain iron biliverdin it is triplet state with 4.86 kcal/mol more stable than the singlet state. The potential energy surface suggests that a spin inversion take place during the course of reaction after TS. The ring opening process in the six-coordinated pathway is in overall -2.26 kcal/mol exothermic with a kinetic barrier of 9.76 kcal/mol. In the five-coordinated pathway the reactant and product are in the ground triplet state. In this path, hydroxyl ion attacks the iron center to produce a complex, which is only 1.59 kcal/mol more stable than when OH(-) directly attacks the macrocycle. The activation barrier for the conversion of iron hydroxy species to the iron biliverdin complex by a rebound mechanism is estimated to be 32.68 kcal/mol. Large barrier for rebound mechanism, small barrier of 4.18 kcal/mol for ring opening process of the hydroxylated macrocycle, and relatively same stabilities for complexes resulted by the attack of nucleophile to the iron and macrocycle indicate that five-coordinated pathway with direct attack of nucleophile to the 5-oxo position of macrocycle might be the path for the conversion of verdoheme to biliverdin.

Mechanistic Study on the Oxidation of Anthocyanidin Synthase by Quantum Mechanical Calculation

Anthocyanidin synthase (ANS), a member of the 2-oxoglutarate-dependent dioxygenase family in flavonoid biosynthesis, catalyzes the conversion of leucoanthocyanidins (e.g. 2R,3S,4S-cis-leucocyanidin, LCD) to flav-2-en-3,4-diols, a direct precursor of colored anthocyanidins via flavan-3,3,4-triols. The detailed oxygenation mechanism of 2R,3S,4S-cis-LCD to flav-2-en-3,4-diols was investigated using the density functional theory method. An initial model for the calculation was constructed from a structure obtained by a 100-ps molecular dynamics simulation of Arabidopsis ANS under physiological conditions. This model consisted of an LCD molecule as the substrate together with an iron atom, two histidine residues, an aspartic acid residue, a succinate, and an oxygen atom as ligands of the iron atom. The results of the calculation indicated that both the C-3 and C-4 positions of LCD can be oxidized, although C-4 oxidation is preferable. The C-3 oxidation required several steps to form flavan-3,3,4-triol: 1) formation of Fe(III)-OH and a substrate C-3 radical via hydrogen atom abstraction by Fe(IV)=O, 2) formation of a C-3 ketone and a water molecule, 3) addition of OH(-) into the C-3 position of the ketone, and 4) addition of H(+) to form flavan-3,3,4-triol. On the other hand, C-4 oxidation of 2R,3S,4S-cis-LCD resulted in the direct formation of 2R,3R-trans-dihydroquercetin. These results suggest that the oxidation at C-3 of LCD, a key reaction for coloring in anthocyanin biosynthesis, can be regarded as a "side reaction" from the viewpoint of quantum mechanics of enzymatic reactions. Molecular evolutional implications of ANS and related proteins are discussed in terms of reaction dynamics.

Effect of the axial cysteine ligand on the electronic structure and reactivity of high-valent iron(IV) oxo-porphyrins (Compound I): A theoretical study

The effect of axial ligands on the reactivity of high-valent iron(IV) oxo-porphyrins (Compound I) was investigated using the B3LYP hybrid density functional method. We studied alkane hydroxylation using four models: Compound I with thiolate, imidazole, phenolate, and chloride anions as axial ligands. The first three ligands were employed as models for cysteinate, histidine, and tyrosinate, respectively. Our calculations show that anionic ligands and neutral ligands favor different electronic states for stationary points in the reaction coordinate, and the calculated energy barrier and energy of several reaction intermediates show similar values. A remarkable effect of axial ligands was found in the final product release step. Our calculations show that the thiolate ligand weakens a bond between heme and an alcohol. In contrast, the imidazole ligand significantly increases the interaction between heme and an alcohol, which causes the catalytic cycle to be less efficient.

Electronic structure of compound I in human isoforms of cytochrome P450 from QM/MM modeling

Human cytochromes P450 play a vital role in drug metabolism. The key step in substrate oxidation involves hydrogen atom abstraction or C=C bond addition by the oxygen atom of the Compound I intermediate. The latter has three unpaired electrons, two on the Fe-O center and one shared between the porphyrin ring and the proximal cysteinyl sulfur atom. Changes in its electronic structure have been suggested to affect reactivity. The electronic and geometric structure of Compound I in three important human subfamilies of cytochrome P450 (P450, 2C, 2B, and 3A) that are major contributors to drug metabolism is characterized here using combined quantum mechanical/molecular mechanical (QM/MM) calculations at the B3LYP:CHARMM27 level. Compound I is remarkably similar in all isoforms, with the third unpaired electron located mainly on the porphyrin ring, and this prediction is not very sensitive to details of the QM/MM methodology, such as the DFT functional, the basis set, or the size of the QM region. The presence of substrate also has no effect. The main source of variability in spin density on the cysteinyl sulfur (from 26 to 50%) is the details of the system setup, such as the starting protein geometry used for QM/MM minimization. This conformational effect is larger than the differences between human isoforms, which are therefore not distinguishable on electronic grounds, so it is unlikely that observed large differences in substrate selectivity can be explained to a large extent in these terms.

A Theoretical Study on the Substrate Deacylation Mechanism of Class C β-Lactamase

The whole reaction of the deacylation of class C beta-lactamase was investigated by performing quantum chemical calculations under physiological conditions. In this study, the X-ray crystallographic structure of the inhibitor moxalactam-bound class C beta-lactamase (Patera et al. J. Am. Chem. Soc. 2000, 122, 10504-10512.) was utilized and moxalactam was changed into the substrate cefaclor. A model for quantum chemical calculations was constructed using an energy-minimized structure of the substrate-bound enzyme obtained by molecular mechanics calculation, in which the enzyme was soaked in thousands of TIP3P water molecules. It was found that the deacylation reaction consisted of two elementary processes. The first process was formation of a tetrahedral intermediate, which was initiated by the activation of catalytic water by Tyr150, and the second process was detachment of the hydroxylated substrate from the enzyme, which associated with proton transfer from the side chain of Lys67 to Ser64O(gamma). The first process is a rate-determining process, and the activation energy was estimated to be 30.47 kcal/mol from density functional theory calculations considering electron correlation (B3LYP/6-31G**). The side chain of Tyr150 was initially in a deprotonated state and was stably present in the active site of the acyl-enzyme complex, being held by Lys67 and Lys315 cooperatively.

Metabolism of 5'alpha,8'-cycloabscisic acid, a highly potent and long-lasting abscisic acid analogue, in radish seedlings

We synthesized 5'alpha,8'-cycloabscisic acid (CycloABA), a highly potent and long-lasting abscisic acid (ABA) analogue, by a different method from that reported before. CycloABA fed to radish seedlings had more metabolic tolerance than ABA. The major metabolite of CycloABA was the glucose conjugate, which was the minor metabolite of ABA. The 8'-hydroxylated metabolite and its cyclized isomer, which were major metabolites of ABA, were not found as metabolites of CycloABA. The present results suggest that the highly potent and long-lasting activity of CycloABA is caused by resistance to ABA 8'-hydroxylase, and that CycloABA is partially metabolized to the glucose conjugate by ABA glucosyltransferase.

Synthesis of chlorophyll b: localization of chlorophyllide a oxygenase and discovery of a stable radical in the catalytic subunit

BACKGROUND

Assembly of stable light-harvesting complexes (LHCs) in the chloroplast of green algae and plants requires synthesis of chlorophyll (Chl) b, a reaction that involves oxygenation of the 7-methyl group of Chl a to a formyl group. This reaction uses molecular oxygen and is catalyzed by chlorophyllide a oxygenase (CAO). The amino acid sequence of CAO predicts mononuclear iron and Rieske iron-sulfur centers in the protein. The mechanism of synthesis of Chl b and localization of this reaction in the chloroplast are essential steps toward understanding LHC assembly.

RESULTS

Fluorescence of a CAO-GFP fusion protein, transiently expressed in young pea leaves, was found at the periphery of mature chloroplasts and on thylakoid membranes by confocal fluorescence microscopy. However, when membranes from partially degreened cells of Chlamydomonas reinhardtii cw15 were resolved on sucrose gradients, full-length CAO was detected by immunoblot analysis only on the chloroplast envelope inner membrane. The electron paramagnetic resonance spectrum of CAO included a resonance at g = 4.3, assigned to the predicted mononuclear iron center. Instead of a spectrum of the predicted Rieske iron-sulfur center, a nearly symmetrical, approximately 100 Gauss peak-to-trough signal was observed at g = 2.057, with a sensitivity to temperature characteristic of an iron-sulfur center. A remarkably stable radical in the protein was revealed by an isotropic, 9 Gauss peak-to-trough signal at g = 2.0042. Fragmentation of the protein after incorporation of 125I- identified a conserved tyrosine residue (Tyr-422 in Chlamydomonas and Tyr-518 in Arabidopsis) as the radical species. The radical was quenched by chlorophyll a, an indication that it may be involved in the enzymatic reaction.

CONCLUSION

CAO was found on the chloroplast envelope and thylakoid membranes in mature chloroplasts but only on the envelope inner membrane in dark-grown C. reinhardtii cells. Such localization provides further support for the envelope membranes as the initial site of Chl b synthesis and assembly of LHCs during chloroplast development. Identification of a tyrosine radical in the protein provides insight into the mechanism of Chl b synthesis.

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.

Mechanism and structure–reactivity relationships for aromatic hydroxylation by cytochrome P450

Cytochrome P450 enzymes play a central role in drug metabolism, and models of their mechanism could contribute significantly to pharmaceutical research and development of new drugs. The mechanism of cytochrome P450 mediated hydroxylation of aromatics and the effects of substituents on reactivity have been investigated using B3LYP density functional theory computations in a realistic porphyrin model system. Two different orientations of substrate approach for addition of Compound I to benzene, and also possible subsequent rearrangement pathways have been explored. The rate-limiting Compound I addition to an aromatic carbon atom proceeds on the doublet potential energy surface via a transition state with mixed radical and cationic character. Subsequent formation of epoxide, ketone and phenol products is shown to occur with low barriers, especially starting from a cation-like rather than a radical-like tetrahedral adduct of Compound I with benzene. Effects of ring substituents were explored by calculating the activation barriers for Compound I addition in the meta and para-position for a range of monosubstituted benzenes and for more complex polysubstituted benzenes. Two structure-reactivity relationships including 8 and 10 different substituted benzenes have been determined using (i) experimentally derived Hammett sigma-constants and (ii) a theoretical scale based on bond dissociation energies of hydroxyl adducts of the substrates, respectively. In both cases a dual-parameter approach that employs a combination of radical and cationic electronic descriptors gave good relationships with correlation coefficients R2 of 0.96 and 0.82, respectively. These relationships can be extended to predict the reactivity of other substituted aromatics, and thus can potentially be used in predictive drug metabolism models.

NMR study on the structural changes of cytochrome P450cam upon the complex formation with putidaredoxin. Functional significance of the putidaredoxin-induced structural changes

We investigated putidaredoxin-induced structural changes in carbonmonoxy P450cam by using NMR spectroscopy. The resonance from the beta-proton of the axial cysteine was upfield shifted by 0.12 ppm upon the putidaredoxin binding, indicating that the axial cysteine approaches to the heme-iron by about 0.1 A. The approach of the axial cysteine to the heme-iron would enhance the electronic donation from the axial thiolate to the heme-iron, resulting in the enhanced heterolysis of the dioxygen bond. In addition to the structural perturbation on the axial ligand, the structural changes in the substrate and ligand binding site were observed. The resonances from the 5-exo- and 9-methyl-protons of d-camphor, which were newly identified in this study, were upfield shifted by 1.28 and 0.20 ppm, respectively, implying that d-camphor moves to the heme-iron by 0.15-0.7 A. Based on the radical rebound mechanism, the approach of d-camphor to the heme-iron could promote the oxygen transfer reaction. On the other hand, the downfield shift of the resonance from the gamma-methyl group of Thr-252 reflects the movement of the side chain away from the heme-iron by approximately 0.25 A. Because Thr-252 regulates the heterolysis of the dioxygen bond, the positional rearrangement of Thr-252 might assist the scission of the dioxygen bond. We, therefore, conclude that putidaredoxin induces the specific heme environmental changes of P450cam, which would facilitate the oxygen activation and the oxygen transfer reaction.

A theoretical study on the mechanism of camphor hydroxylation by compound I of cytochrome p450

Mechanistic and energetic aspects for the conversion of camphor to 5-exo-hydroxycamphor by the compound I iron-oxo species of cytochrome P450 are discussed from B3LYP DFT calculations. This reaction occurs in a two-step manner along the lines that the oxygen rebound mechanism suggests. The activation energy for the first transition state of the H atom abstraction at the C5 atom of camphor is computed to be more than 20 kcal/mol. This H atom abstraction is the rate-determining step in this hydroxylation reaction, leading to a reaction intermediate that involves a carbon radical species and the iron-hydroxo species. The second transition state of the rebound step that connects the reaction intermediate and the product alcohol complex lies a few kcal/mol below that for the H atom abstraction on the doublet and quartet potential energy surfaces. This energetic feature allows the virtually barrierless recombination in both spin states, being consistent with experimentally observed high stereoselectivity and brief lifetimes of the reaction intermediate. The overall energetic profile of the catalytic mechanism of camphor hydroxylation particularly with respect to why the high activation energy for the H atom abstraction is accessible under physiological conditions is also considered and calculated. According to a proton source model involving Thr252, Asp251, and two solvent water molecules (Biochemistry 1998, 37, 9211), the energetics for the conversion of the iron-peroxo species to compound I is studied. A significant energy over 50 kcal/mol is released in the course of this dioxygen activation process. The energy released in this chemical process is an important driving force in alkane hydroxylation by cytochrome P450. This energy is used for the access to the high activation energy for the H atom abstraction.

Reactivities of organic phase biosensors: 6. Square-wave and differential pulse studies of genetically engineered cytochrome P450(cam) (CYP101) bioelectrodes in selected solvents

Cytochrome P450(cam) (CYP101) bioelectrodes suitable for application in organic phases were prepared from genetically engineered CYP101 and vesicular dispersions of didodecyldimethylammonium bromide. The amperometric biosensor system was characterised under anaerobic conditions by cyclic and square-wave voltammetric methods. Cyclic- and square-wave-voltammetry studies showed that the biosensors exhibited direct reversible electron transfer between the haem iron atom and the glassy carbon electrode surface. The formal redox potential estimated for the electrode in acetonitrile was -380 mV/Ag-AgCl. The formal potential shifted anodically as the organic phase biosensor responded irreversibly to substrate (camphor) under anaerobic and aerobic conditions in acetonitrile. Differential pulse analysis of the reactivities of the CYP101 enzyme electrode confirmed the square-wave voltammetry result, which showed that the binding of substrate decreased the redox potential necessary for initiating the monooxygenation reaction of cytochrome P450(cam).

Crystallographic and spectroscopic characterization of a nonheme Fe(IV)-O complex

Following the heme paradigm, it is often proposed that dioxygen activation by nonheme monoiron enzymes involves an iron(IV)=oxo intermediate that is responsible for the substrate oxidation step. Such a transient species has now been obtained from a synthetic complex with a nonheme macrocyclic ligand and characterized spectroscopically. Its high-resolution crystal structure reveals an iron-oxygen bond length of 1.646(3) angstroms, demonstrating that a terminal iron(IV)=oxo unit can exist in a nonporphyrin ligand environment and lending credence to proposed mechanisms of nonheme iron catalysis.

Revisiting the mechanism of P450 enzymes with the radical clocks norcarane and spiro[2,5]octane

Norcarane (1) and spiro[2.5]octane (2) yield different product distributions depending on whether they are oxidized via concerted, radical, or cationic mechanisms. For this reason, these two probes were used to investigate the mechanisms of hydrocarbon hydroxylation by two mammalian and two bacterial cytochrome P450 enzymes. Products indicative of a radical intermediate with a lifetime ranging from 16 to 52 ps were detected during the oxidation of norcarane by P450(cam) (CYP101), P450(BM3) (CYP102), CYP2B1, and CYP2E1. Trace amounts of the cation rearrangement product were observed with norcarane for all but CYP2E1, while no cation or radical rearrangement products were observed for spiro[2.5]octane. The results for the oxidation of norcarane with a radical rearrangement rate of 2 x 10(8) s(-1) are consistent with the involvement of a two-state radical rebound mechanism, while for the slower (5 x 10(7) s(-1)) spiro[2,5]oct-4-yl radical rearrangement products were beyond detection. Taken together with earlier data for the hydroxylation of bicyclo[2.1.0]pentane, which also suggested a 50 ps radical lifetime, these three structurally similar and functionally simple substrates show a consistent pattern of rearrangement that supports a radical rebound mechanism for this set of cytochrome P450 enzymes.

De novo structural prediction of transition metal complexes: application to technetium

De novo structural prediction of transition metal complexes is investigated. Technetium complexes are chosen given their importance in medical imaging and nuclear waste remediation and for the chemical diversity they display. A new conformational searching algorithm (LIGB) for transition metals is described that allows one to search for different conformational and geometric isomers within a single simulation. In the preponderance of cases, both conformational searching techniques (LIGB and high-temperature molecular dynamics/simulated annealing) provide comparable results, while LIGB is superior for macrocyclic complexes. A genetic algorithm-optimized PM3(tm) parametrization for Tc is compared with the standard implementation and found to yield a significant improvement in predictive ability for the most prevalent Tc structural motifs. The utility of a coupled molecular mechanics-semiempirical quantum mechanics protocol is demonstrated for very rapid, efficient, and effective de novo prediction of transition metal complex geometries.

Synthetic analogue of the [Fe(2)(mu-OH)(2)(mu-O(2)CR)](3+) core of soluble methane monooxygenase hydroxylase via synthesis and dioxygen reactivity of carboxylate-bridged diiron(II) complexes

We describe the synthesis and dioxygen reactivity of diiron(II) tetracarboxylate complexes [Fe(2)(mu-O(2)CAr(Tol))(2)(O(2)CAr(Tol))(2)(N,N-Me(2)en)(2)] (2) and [Fe(2)(mu-O(2)CAr(Tol))(2)(O(2)CAr(Tol))(2)(N,N-Bn(2)en)(2)] (6), where Ar(Tol)CO(2)(-) = 2,6-di(p-tolyl)benzoate. These complexes were prepared as models for the diiron(II) center in the hydroxylase component of soluble methane monooxygenase (MMOH). Compound 6 reacts with dioxygen to afford PhCHO in approximately 60(5)% yield, following oxidative N-dealkylation of the pendant benzyl group on the diamine ligand. The diiron(III) complex [Fe(2)(mu-OH)(2)(mu-O(2)CAr(Tol))(O(2)CAr(Tol))(3)(N-Bnen)(N,N-Bn(2)en)] (8) was isolated from the reaction mixture. The 4.2 K Mössbauer spectrum of 8 displays a single quadrupole doublet with parameters delta = 0.48(2) mm s(-1) and Delta E(Q) = 0.61(2) mm s(-1). The [Fe(2)(mu-OH)(2)(mu-O(2)CR)](3+) core structure in 8 matches that of the fully oxidized form of MMOH. The conversion of 6 to 8 closely parallels the chemistry of MMOH in which an O(2)-derived oxygen atom is inserted into the C-H bond of methane. Several reaction pathways are considered to account for this novel chemical transformation, and these are compared with mechanistic frameworks previously developed for related cytochrome P450 and copper(I) dioxygen chemistry.

[Molecular science of the living organism: the case of G-proteins]

The concept of Molecular Science of the Living Organism was described, where the living state is explained as the purposive flows of the quantum mechanically controlled chemical reaction systems which support the homeostasis of the living organism. In the 21st century, the post genomic sequence era, the concept may be a self-evident truth. Molecular Science of the Living Organism was presented in the case of G-proteins: i.e., the atomically controlled mechanism of 1. the carcinogenesis which originates from the point mutation of ras p21, 2. the activation of a receptor protein at the cell membrane, especially in the case of bacteriorhodopsin, 3. the activation of an inactive G-protein by the activated receptor protein.