Reaction profile of the last step in cytochrome P-450 catalysis revealed by studies of model complexes

 A series of oxoiron(IV) porphyrin cation radical complexes was investigated as compound I analogs of cytochrome P-450. Both the spectroscopic features and the reactivities of the complexes in oxygen atom transfer to olefins were examined as a function of only one variable, the axial ligand trans to the oxoiron(IV) bond. The results disclosed two important kinetic steps – electron transfer from olefin to oxoiron(IV) and intramolecular electron transfer from metal to porphyrin radical – which are affected differently by the axial ligands. The large kinetic barrier of the latter step in the reaction of olefins with the perchlorato-bound oxoiron(IV) porphyrin cation radical complex enabled the trapping of a reaction intermediate in which the metal, but not the porphyrin radical, is reduced. The first electron transfer step is probably followed by σ-bond formation, which readily accounts for formation of isomerized organic products at low temperatures. It is finally postulated that part of the enhanced oxygenation activities of cytochrome P-450 monooxygenases and chloroperoxidases is due to a lowering of the energy barrier for the second electron transfer step via participation of their redox-active cysteinate ligand. Received: 16 January 1997 / Accepted: 24 May 1997     

Nonheme iron(II) complexes of macrocyclic ligands in the generation of oxoiron(IV) complexes and the catalytic epoxidation of olefins

Mononuclear nonheme oxoiron(IV) complexes bearing 15-membered macrocyclic ligands were generated from the reactions of their corresponding iron(II) complexes and iodosylbenzene (PhIO) in CH(3)CN. The oxoiron(IV) species were characterized with various spectroscopic techniques such as UV-vis spectrophotometer, electron paramagnetic resonance, electrospray ionization mass spectrometer, and resonance Raman spectroscopy. The oxoiron(IV) complexes were inactive in olefin epoxidation. In contrast, when iron(II) or oxoiron(IV) complexes were combined with PhIO in the presence of olefins, high yields of epoxide products were obtained. These results indicate that in addition to the oxoiron(IV) species, there must be at least one more active oxidant (e.g., Fe(IV)-OIPh adduct or oxoiron(V) species) that effects the olefin epoxidation. We have also demonstrated that the ligand environment of iron catalysts is an important factor in controlling the catalytic activity as well as the product selectivity in the epoxidation of olefins by PhIO.     

How does the push/pull effect of the axial ligand influence the catalytic properties of Compound I of catalase and cytochrome P450?

Density functional theory (DFT) calculations on the chemoselective epoxidation versus hydroxylation reactions of propene by oxoiron porphyrin models mimicking the active sites of catalase, cytochrome P450 (P450) and horseradish peroxidase Compound I (CpdI) are presented. The catalase reactions are concerted and proceed via two-state reactivity patterns on competing doublet and quartet spin state surfaces, but the lowest barrier is the one leading to epoxide products on the doublet spin surface. The results are compared with earlier DFT studies of models of cytochrome P450, horseradish peroxide (HRP), taurine/alpha-ketoglutarate dioxygenase and some synthetic oxoiron catalysts. The catalase barriers are midway in between those obtained for HRP and P450 models, so that tyrosinate ligated heme systems should be able to catalyze C-H hydroxylation and C=C epoxidation reactions. We show that for heme systems the barrier height of epoxidation linearly correlates with the electron affinity of Compound I as expected from the electron transfer mechanism of the rate determining step. Our studies show that the axial ligand does not influence the chemoselectivity of a reaction but that it does regulate the barrier heights and rate constants. Finally, we estimated the effect of the axial ligand on the oxoiron group and derived that it contributes from a field effect due to the charge of the ligand and a quantum mechanical effect as a result of orbital mixing. In catalase, the major component is the field effect, while the quantum mechanical effect is negligible. This is in contrast to P450 CpdI, where both effects are of similar order of magnitude.     

A Computer-Aided Radiopharmaceutical Drug Design Study Using Ab Initio and Molecular Mechanics Methods

This is an investigation of technetium ligands and their complexes with [TcO]³⁺ using ab initio population analysis and molecular mechanics conformational searching methods. Calculated atomic electronic populations on the technetium atom in complexes with a number of ligands gauge the degree of covalent bonding between technetium and these ligands. Here a reduction in the positive charge on the [TcO]³⁺ moiety by complexation with a given ligand is correlated with covalent bonding. Our ab initio results suggest that ligands with more sulphur atoms have better covalent bonding to technetium than do other ligands. A conformational analysis of the uncomplexed ligands indicates that conformational reorganization before complexation correlates inversely with stable complex formation. This conformational analysis shows that ligands with ethylene carbonyl bridges have low energy conformations closer to the final complexation geometries than do ligands with ethylene, propylene or propylene carbonyl bridges. The presence of these low energy conformations facilitates a faster complexation of the ethylene carbonyl [TcO]³⁺ moiety. This result produces a kinetic explaination why ethylene carbonyl bridged ligands form stable complexes while many other ligands do not [1]. The conclusion is that kinetic and thermodynamic considerations play a role in stable complex formation between these ligands and technetium.     

Oxoiron(IV) porphyrin π-cation radical complexes with a chameleon behavior in cytochrome P450 model reactions

There is an intriguing, current controversy on the involvement of multiple oxidizing species in oxygen transfer reactions by cytochromes P450 and iron porphyrin complexes. The primary evidence for the multiple oxidants theory was that products and/or product distributions obtained in the catalytic oxygenations were different depending on reaction conditions such as catalysts, oxidants, and solvents. In the present work, we carried out detailed mechanistic studies on competitive olefin epoxidation, alkane hydroxylation, and C=C epoxidation versus allylic C–H hydroxylation in olefin oxygenation with in situ generated oxoiron(IV) porphyrin -cation radicals (1) under various reaction conditions. We found that the products and product distributions were markedly different depending on the reaction conditions. For example, 1 bearing different axial ligands showed different product selectivities in competitive epoxidations of cis-olefins and trans-olefins and of styrene and para-substituted styrenes. The hydroxylation of ethylbenzene by 1 afforded different products, such as 1-phenylethanol and ethylbenzoquinone, depending on the axial ligands of 1 and substrates. Moreover, the regioselectivity of C=C epoxidation versus C–H hydroxylation in the oxygenation of cyclohexene by 1 changed dramatically depending on the reaction temperatures, the electronic nature of the iron porphyrins, and substrates. These results demonstrate that 1 can exhibit diverse reactivity patterns under different reaction conditions, leading us to propose that the different products and/or product distributions observed in the catalytic oxygenation reactions by iron porphyrin models might not arise from the involvement of multiple oxidizing species but from 1 under different circumstances. This study provides strong evidence that 1 can behave like a chameleon oxidant that changes its reactivity and selectivity under the influence of environmental changes.Electronic Supplementary Material Supplementary material is available for this article at .     

Comparing the spectroscopic and electrochemical properties of ruthenium and osmium complexes of the tridentate polyazine ligands 2,2′:6′,2″-terpyridine and 2,3,5,6-tetrakis(2-pyridyl)pyrazine

A number of osmium and ruthenium complexes of the tridentate ligands 2,2′:6′,2″-terpyridine (tpy) and 2,3,5,6-tetrakis(2-pyridyl)pyrazine (tpp) have been prepared and characterized by our laboratory. All these complexes possess metal based oxidations and ligand based reductions localized on each polyazine ligand. Polymetallic complexes bridged by the tpp ligand exhibit two sequential tpp based reductions prior to the reduction of other polyazine ligands in these complexes. The spectroscopy of these complexes is dominated by ligand based π-π^* transitions in the ultraviolet and MLCT (metal-to-ligand charge transfer) bands terminating on each polyzine ligand in the visible. For the complexes reported herein the lowest lying optical transitionis a M → BL CT band. For most of the complexes reported, occupation of this excited state gives rise to an observable emission at room temperature. For ruthenium complexes of these tridentate ligands, the presence of a low-lying LF state shortens the excited state lifetimes of these chromophores. This gives rise to ruthenium complexes that display shorter excited state lifetimes than the analogous osmium based systems. Incorporation of tpp based chromophores into polymetallic frameworks leads to the production of bimetallic species with long-lived excited states, 100 ns at room temperature. This makes these chromophores good candidates for the development of stereochemically defined supramolecular complexes. It is possible to measure an electrochemical HOMO-LUMO energy gap and a correlation between this electrochemically measured energy gap and the spectroscopic energy associated with this HOMO→LUMO transition are reported herein (HOMO== highest occupied molecular orbital, LUMO = lowest unoccupied molecular orbital).     

Oxoiron(IV) complexes of the tris(2-pyridylmethyl)amine ligand family: effect of pyridine α-substituents

The oxoiron(IV) complexes of two 6-substituted tris(2-pyridylmethyl)amine ligand derivatives have been generated and characterized with respect to their spectroscopic and reactivity properties. The introduction of an α-substituent maintains the low-spin nature of the oxoiron(IV) unit but weakens the ligand field, as evidenced by red shifts in its characteristic near-IR chromophore. While its hydrogen-atom abstraction ability is only slightly affected, the oxo-transfer reactivity of the oxoiron(IV) center is significantly enhanced relative to that of the parent complex. These results demonstrate that the ligand environment plays a key role in modulating the reactivity of this important biological oxidant. Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users. An erratum to this article can be found at     

Inhibition of the catecholase activity of biomimetic dinuclear copper complexes by kojic acid

The inhibition of the catechol oxidase activity exhibited by three dinuclear copper(II) complexes, derived from different diaminotetrabenzimidazole ligands, by kojic acid [5-hydroxy-2-(hydroxymethyl)-γ-pyrone] has been studied. The catalytic mechanism of the catecholase reaction proceeds in two steps and for both of these inhibition by kojic acid is of competitive type. The inhibitor binds strongly to the dicopper(II) complex in the first step and to the dicopper-dioxygen adduct in the second step, preventing in both cases the binding of the catechol substrate. Binding studies of kojic acid to the dinuclear copper(II) complexes and a series of mononuclear analogs, carried out spectrophotometrically and by NMR, enable us to propose that the inhibitor acts as a bridging ligand between the metal centers in the dicopper(II) catalysts. Received: 23 August 1999 / Accepted: 20 January 2000     

Parallel mechanistic studies on the counterion effect of manganese salen and porphyrin complexes on olefin epoxidation by iodosylarenes

Counterions of manganese(III) porphyrin complexes influence diastereoselectivity in cis-stilbene epoxidation and product distribution in cyclohexene epoxidation markedly. In the epoxidation of cis-stilbene by iodosylbenzene carried out in a solvent mixture of CH(3)CN and CH(2)Cl(2), trans-stilbene oxide is the major product in the reaction of manganese complexes bearing a ligating anion (i.e., Cl(-)), whereas cis-stilbene oxide is the dominant product in the reactions of manganese complexes bearing a poorly-ligating anion (i.e., CF(3)SO(4)(-)). In cyclohexene epoxidation, the yields of allylic oxidation products such as cyclohexenol and cyclohexenone are higher when the counterion of the manganese catalysts is Cl(-) than when the counterion is CF(3)SO(4)(-). The product selectivities are also dependent on the nature of iodosylarenes and the axial and porphyrin ligands of the manganese porphyrin catalysts. The observation that product selectivities are different depending on the iodosylarenes may indicate the involvement of multiple oxidants in oxygen atom transfer reactions. These results are compared with those observed in manganese salen-catalyzed epoxidation of olefins by iodosylarenes.     

A new hybridocytochemical method based on mercurated nucleic acid probes and sulfhydryl-hapten ligands

Summary In the preceding paper, a method to detect specific DNA sequences with mercurated nucleic acid probes and sulfhydryl-hapten ligands has been described. Due to the instability of the bond between mercury and a negatively charged sulfhydryl-hapten ligand (trinitrophenyl-glutathione), the in situ formed hybrid could not be detected. On basis of model system experiments it was suggested that this mercury-sulfhydryl bond could be stabilized by an extra polar interaction between ligand and nucleic acid. This was achreved by reversing the net charge of the ligand. Such ligands were synthesized by reacting aliphatic diamines to the carboxyl groups of Tnp-glutathione using a water soluble carbodiimide. Gel chromatographic analysis of mercurated polynucleotide-ligand complexes showed that the stability of the mercury-sulfhydryl bond is increased by the reversal of the net charge of the ligand.In situ hybridized mercurated mouse satellite DNA to mouse liver nuclei and mercurated kinetoplast cRNA hybridized to Crithidia fasciculata were immunocytochemically detected after the introduction of these positively charged ligands.The described method is applicable for RNA and DNA probes. It has a sensitivity comparable to other non-autoradiographic methods, is relatively simple to perform and can be carried out with ordinary laboratory chemicals.This investigation was supported by the Netherlands Foundation for Medical Research Fungo (grant nr 13-54-21)     

How are hydrogen bonds modified by metal binding?

We have used density functional theory calculations to investigate how the hydrogen-bond strength is modified when a ligand is bound to a metal using over 60 model systems involving six metals and eight ligands frequently encountered in metalloproteins. We study how the hydrogen-bond geometry and energy vary with the nature of metal, the oxidation state, the coordination number, the ligand involved in the hydrogen bond, other first-sphere ligands, and different hydrogen-bond probe molecules. The results show that, in general, the hydrogen-bond strength is increased for neutral ligands and decreased for negatively charged ligands. The size of the effect is mainly determined by the net charge of the metal complex, and all effects are typically decreased when the model is solvated. In water solution, the hydrogen-bond strength can increase by up to 37 kJ/mol for neutral ligands, and that of negatively charged ligands can increase (for complexes with a negative net charge) or decrease (for positively charged complexes). If the net charge of the complex does not change, there is normally little difference between different metals or different types of complexes. The only exception is observed for sulphur-containing ligands (Met and Cys) and if the ligand is redox-active (e.g. high-valence Fe–O complexes).     

Bonding and electronic structures in W@Au12AE complexes (AE= NO+, CO, BF, CN–, or BO–): analogies among ligands isoelectronic to carbon monoxide

A theoretical study on the geometries and electronic structures of W@Au₁₂AE (AE=NO⁺, BF, CN^–, or BO^–) was carried out to gain insight into interactions between W@Au₁₂ and ligands isoelectronic with CO. The best configuration for the adsorption site is on-top type for all five complexes. After complexing with boron ligands (BF or BO^–), the axial Au–W bond distance in W@Au₁₂ is lengthened notably, but NO⁺ has the opposite effect on the axial Au–W bond. A charge transfer and energy decomposition analysis shows that the metal–ligand bonds have enhanced σ-donation strength from NO⁺ to BO^–. Furthermore, the A–E bond strength in the complexes becomes weaker with stronger π-back-donation interactions. Finally, W@Au₁₂CO has the largest HOMO–LUMO gap, making it the most stable in terms of kinetic stability.     

Synthetic analogs of the active site of cytochrome P-450cam characterization of thiolpeptide fragment-hemin complexes by optical and ESR spectrometries

Thiol-containing penta (Leu-Ala-Cys-Ser-Leu) and nona (Leu-Ala-Cys-Ser-Leu-Ile-Glu-Ser-Leu) peptides corresponding to residues 132-136 and 132-140, respectively, of apo-P450 from Psuedomonas putida were synthesized to examine heme-binding by the enzymes in the oxidized (ferric) form. The peptide-hemin complexes prepared in solution were characterized by their optical and ESR spectra. In these complexes without nitrogenous ligands, no ferric complexes in the low-spin state were observed, suggesting that simultaneous axial coordination of Cys-134 (thiolate) and Ser-139 (hydroxyl) of apo-P450cam to the heme is unlikely to occur. In the presence of nitrogenous ligands, such as py, Im and MeIm, the resulting complexes were good models of apo-P450cam-nitrogenous ligand adducts in the low-spin ferric state retaining a thiolate-Fe(III)-nitrogen axial coordination mode, as judged by their spectral pattern as well as by crystal field analysis of ESR g-values.     

Enantioselective epoxidation of non-functionalized alkenes using carbohydrate based salen-Mn(III) complexes

Three new salen ligands with carbohydrate moieties were prepared from a salicylaldehyde derivative obtained by reaction of 1,2:5,6-di-O-isopropylidene-alpha-D-glucofuranose with 3-tert-butyl-5-(chloro-methyl)-2-hydroxybenzaldehyde. These ligands were coordinated with Mn(III) to give three chiral salen-Mn(III) complexes. The complexes were characterized and employed in the asymmetric epoxidation of unfunctionalized alkenes. Catalytic results showed that although there are no chiral groups on the diimine bridge, these complexes had some enantioselectivity, which indicates the carbohydrate moiety has an asymmetric inducing effect in the epoxidation reaction.     

The effect of axial ligand plane orientation on the contact and pseudocontact shifts of low-spin ferriheme proteins

 The effect of axial ligand nodal plane orientation on the contact and pseudocontact shifts of a symmetrical low-spin octamethylferriheme center has been calculated as a function of the angle of the axial ligand. Simple Hückel techniques have been used to estimate the contact contribution, and values obtained from model hemes, together with counter-rotation of the g-tensor, have been used to estimate the pseudocontact contribution, for the eight β-pyrrole methyl and four meso-H positions. It is found that the maximum and minimum contact shifts occur when the axial ligand is aligned at an angle of ±15° to the meso-H axes of the heme, rather than when the axial ligand plane lies along the porphyrin nitrogens, as assumed previously by some investigators. For systems having one planar axial ligand or two ligands in parallel planes, the contact and pseudocontact contributions at the meso-H positions are comparable in size (at least on the basis of simple Hückel estimates), while the contact contribution clearly dominates the isotropic shifts of the heme methyls. Allowing for the substituent effect of the 2,4-vinyls of protohemin, or the 2,4-thioethers of hemin c, as well as the average diamagnetic shifts of the heme methyls and meso-H, plots of the predicted shifts as a function of axial ligand nodal plane orientation have been constructed for hemin b- and c-containing proteins. Excellent agreement in the order of shifts, and reasonable agreement in the sizes of the observed shifts, is observed in the majority of the ferriheme proteins for which the methyl and meso-H resonances have been assigned and proton shifts reported. Plots have also been constructed for hemin c-containing proteins having the two axial ligand nodal planes oriented at relative angles of 40°, 70°, and 80°. Excellent agreement in the order of shifts, and reasonable agreement in the magnitudes of the observed shifts, is observed in all cases of bacterial cytochromes which do not fit the plots that assume the ligands are in parallel planes, except one – the cytochrome c-552 of Nitrosomonas europae. Except for this case, where the order of the predicted methyl shifts at any angle of the axial ligands disagrees with the observed, the reasons can usually be attributed to a large dihedral angle between two axial ligand nodal planes, to strong H-bonding interactions involving His and/or CN^– ligands, or to off-axis binding of one (or both) axial ligand(s). Ruffling of the porphyrin ring may also contribute to the contact shift in as yet undefined ways. Hence, despite the simplicity of the calculations, the agreement with observed data is highly satisfying and the concept of the importance of axial ligand plane orientation on the observed proton shifts of heme proteins is fully confirmed. Received: 15 June 1998 / Accepted: 6 August 1998     

Mimicking the electron donor side of Photosystem II in artificial photosynthesis

This review focuses on our recent efforts in synthetic ruthenium–tyrosine–manganese chemistry mimicking the donor side reactions of Photosystem II. Tyrosine and tryptophan residues were linked to ruthenium photosensitizers, which resulted in model complexes for proton-coupled electron transfer from amino acids. A new mechanistic model was proposed and used to design complexes in which the mechanism could be switched between concerted and step-wise proton-coupled electron transfer. Moreover, a manganese dimer linked to a ruthenium complex could be oxidized in three successive steps, from Mn₂^II,II to Mn₂^III,IV by the photo-oxidized ruthenium sensitizer. This was possible thanks to a charge compensating ligand exchange in the manganese complex. Detailed studies of the ligand exchange suggested that at high water concentrations, each oxidation step is coupled to a proton-release of water-derived ligands, analogous to the oxidation steps of the manganese cluster of Photosystem II.     

Ligand-supported homology modelling of protein binding-sites using knowledge-based potentials

A new approach, MOBILE, is presented that models protein binding-sites including bound ligand molecules as restraints. Initially generated, homology models of the target protein are refined iteratively by including information about bioactive ligands as spatial restraints and optimising the mutual interactions between the ligands and the binding-sites. Thus optimised models can be used for structure-based drug design and virtual screening. In a first step, ligands are docked into an averaged ensemble of crude homology models of the target protein. In the next step, improved homology models are generated, considering explicitly the previously placed ligands by defining restraints between protein and ligand atoms. These restraints are expressed in terms of knowledge-based distance-dependent pair potentials, which were compiled from crystallographically determined protein-ligand complexes. Subsequently, the most favourable models are selected by ranking the interactions between the ligands and the generated pockets using these potentials. Final models are obtained by selecting the best-ranked side-chain conformers from various models, followed by an energy optimisation of the entire complex using a common force-field. Application of the knowledge-based pair potentials proved efficient to restrain the homology modelling process and to score and optimise the modelled protein-ligand complexes. For a test set of 46 protein-ligand complexes, taken from the Protein Data Bank (PDB), the success rate of producing near-native binding-site geometries (rmsd<2.0A) with MODELLER is 70% when the ligand restrains the homology modelling process in its native orientation. Scoring these complexes with the knowledge-based potentials, in 66% of the cases a pose with rmsd <2.0A is found on rank 1. Finally, MOBILE has been applied to two case studies modelling factor Xa based on trypsin and aldose reductase based on aldehyde reductase.     

Influence of solvent, anion and presence of nitrogen in the ring structure on the mechanism of complexation of alkali metal cations with crown ethers

The mechanism of complexation of alkali metal cations with macrocyclic ligands such as the simple crown ethers and the role of desolvation vs. ligand rearrangement are discussed. The unique role of water solvent in the rate-determining step of complexations in aqueous solutions is brought into focus. The competitive role of the anion, which becomes of paramount importance in solvents of low permittivity, is reiterated. Monoazo crown ethers are shown to possess isomeric equilibria in methanol solvent. The rate-determining process for the first step of complexation of these macrocycles with Na+ in methanol appears to be the rearrangement of the ligand through inversion to an exo position of the nitrogen lone electron pair. The rate-determining step of the overall complexation is the entrance of the Na+ into the ring with (possibly) concomitant rotation of the lone electron of the nitrogen to an endo configuration.     

Electrochemical properties of copper complexes with unsubstituted and substituted 1,10-o-phenanthrolines in N,N-dimethylformamide solvent

The electrochemical reduction of a series of copper(II) complexes with 1,10-o-phenanthrolines, namely the 1:1 and 1:2 metal:ligand complexes with 2,9-dimethylphenanthroline, 4,7-dimethylphenanthroline and unsubstituted phenanthroline, respectively, has been studied in N,N-dimethylformamide using platinum electrodes. As to the 1:2 complexes, the effect of the presence of substituents with different electronic and steric effects on the phenanthroline ligands has been studied with the aim of rationalizing the different values of the standard potentials which have been measured. Furthermore, the possibility of electrogenerating neutral species, with a formally zerovalent copper centre, exhibiting different stability depending on the nature of the ligands, has been ascertained. In out solvent medium, 1:1 complexes have been found to be in equilibrium with the corresponding 1:2 complexes. A scheme for the reduction of solutions of these compounds, including the different equilibria associated to the electrode charge transfers, has been outlined.     

The axial ligand effect of oxo-iron porphyrin catalysts. How does chloride compare to thiolate?

We have performed density functional theory calculations on an oxo-iron porphyrin catalyst with chloride as an axial ligand and tested its reactivity toward propene. The reactions proceed via multistate reactivity on competing doublet and quartet spin surfaces. Close-lying epoxidation and hydroxylation mechanisms are identified, whereby in the gas phase the epoxidation reaction is dominant, while in environments with a large dielectric constant the hydroxylation pathways become competitive. By contrast to reactions with thiolate as an axial ligand all low-lying pathways have small ring-closure and rebound barriers, so it is expected that side products and rearrangements will be unlikely with Fe=O(porphyrin)Cl, whereas with Fe=O(porphyrin)SH some side products were predicted. The major differences in the electronic configurations of Fe=O(porphyrin)Cl and Fe=O(porphyrin)SH are due to strong mixing of thiolate orbitals with iron 3d orbitals, a mixing which is much less with chloride as an axial ligand. Predictions of the reactivity of ethylbenzene-h ₁₂ versus ethylbenzene-d ₁₂ are made. Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.