Insight into heme protein redox potential control and functional aspects of six-coordinate ligand-sensing heme proteins from studies of synthetic heme peptides

We describe detailed studies of peptide-sandwiched mesohemes PSMA and PSMW, which comprise two histidine (His)-containing peptides covalently attached to the propionate groups of iron mesoporphyrin II. Some of the energy produced by ligation of the His side chains to Fe in the PSMs is invested in inducing helical conformations in the peptides. Replacing an alanine residue in each peptide of PSMA with tryptophan (Trp) to give PSMW generates additional energy via Trp side chain-porphyrin interactions, which enhances the peptide helicity and stability of the His-ligated state. The structural change strengthened His-FeIII ligation to a greater extent than His-FeII ligation, leading to a 56-mV negative shift in the midpoint reduction potential at pH 8 (Em,8 value). This is intriguing because converting PSMA to PSMW decreased heme solvent exposure, which would normally be expected to stabilize FeII relative to FeIII. This and other results presented herein suggest that differences in stability may be at least as important as differences in porphyrin solvent exposure in governing redox potentials of heme protein variants having identical heme ligation motifs. Support for this possibility is provided by the results of studies from our laboratories comparing the microsomal and mitochondrial isoforms of mammalian cytochrome b5. Our studies of the PSMs also revealed that reduction of FeIII to FeII reversed the relative affinities of the first and second His ligands for Fe (K2III > K1III; K2II < K1II). We propose that this is a consequence of conformational mobility of the peptide components, coupled with the much greater ease with which FeII can be pulled from the mean plane of a porphyrin. An interesting consequence of this phenomenon, which we refer to as "dynamic strain", is that an exogenous ligand can compete with one of the His ligands in an FeII-PSM, a reaction accompanied by peptide helix unwinding. In this regard, the PSMs are better models of neuroglobin, CooA, and other six-coordinate ligand-sensing heme proteins than of stably bis(His)-ligated electron-transfer heme proteins such as cytochrome b5. Exclusive binding of exogenous ligands by the FeII form of PSMA led to positive shifts in its Em,8 value, which increases with increasing ligand strength. The possible relevance of this observation to the function of six-coordinate ligand-sensing heme proteins is discussed.     

Modulation of ligand-field parameters by heme ruffling in cytochromes c revealed by EPR spectroscopy

Electron paramagnetic resonance (EPR) spectra of variants of Hydrogenobacter thermophilus cytochrome c(552) (Ht c-552) and Pseudomonas aeruginosa cytochrome c(551) (Pa c-551) are analyzed to determine the effect of heme ruffling on ligand-field parameters. Mutations introduced at positions 13 and 22 in Ht c-552 were previously demonstrated to influence hydrogen bonding in the proximal heme pocket and to tune reduction potential (E(m)) over a range of 80 mV [Michel, L. V.; Ye, T.; Bowman, S. E. J.; Levin, B. D.; Hahn, M. A.; Russell, B. S.; Elliott, S. J.; Bren, K. L. Biochemistry 2007, 46, 11753-11760]. These mutations are shown here to also increase heme ruffling as E(m) decreases. The primary effect on electronic structure of increasing heme ruffling is found to be a decrease in the axial ligand-field term Δ/λ, which is proposed to arise from an increase in the energy of the d(xy) orbital. Mutations at position 7, previously demonstrated to influence heme ruffling in Pa c-551 and Ht c-552, are utilized to test this correlation between molecular and electronic structure. In conclusion, the structure of the proximal heme pocket of cytochromes c is shown to play a role in determining heme conformation and electronic structure.     

Solution 1h NMR characterization of equilibrium heme orientational disorder with functional consequences in mouse neuroglobin

The solution 1H NMR spectrum of oxidized (met) mouse neuroglobin, metNgb, demonstrates that it is low-spin and hexacoordinate with strong spectral similarities to ferricytochrome b5. The axial ligands are identified as His(F8) and His(E7), with the latter exhibiting an unstrained Fe-His bond. The presence of two sets of resonances is shown to arise from equilibrium heme orientational isomers ( approximately 2:1). The ligation of cyanide is shown to be extraordinarily slow with a factor approximately 2 difference in rate for the two heme orientations. Not only is Ngb the first mamalian globin with equilibrium heme disorder, but the disorder also has additional functional consequences.     

Coherence spectroscopy investigations of the low-frequency vibrations of heme: effects of protein-specific perturbations

Femtosecond coherence spectroscopy is used to probe the low-frequency (20-200 cm(-1)) vibrational modes of heme proteins in solution. Horseradish peroxidase (HRP), myoglobin (Mb), and Campylobacter jejuni globin (Cgb) are compared and significant differences in the coherence spectra are revealed. It is concluded that hydrogen bonding and ligand charge do not strongly affect the low-frequency coherence spectra and that protein-specific deformations of the heme group lower its symmetry and control the relative spectral intensities. Such deformations potentially provide a means for proteins to tune heme reaction coordinates, so that they can perform a broad array of specific functions. Native HRP displays complex spectral behavior above approximately 50 cm(-1) and very weak activity below approximately 50 cm(-1). Binding of the substrate analog, benzhydroxamic acid, leads to distinct changes in the coherence and Raman spectra of HRP that are consistent with the stabilization of a heme water ligand. The CN derivatives of the three proteins are studied to make comparisons under conditions of uniform heme coordination and spin-state. MbCN is dominated by a doming mode near 40 cm(-1), while HRPCN displays a strong oscillation at higher frequency (96 cm(-1)) that can be correlated with the saddling distortion observed in the X-ray structure. In contrast, CgbCN displays low-frequency coherence spectra that contain strong modes near 30 and 80 cm(-1), probably associated with a combination of heme doming and ruffling. HRPNO displays a strong doming mode near 40 cm(-1) that is activated by photolysis. The damping of the coherent motions is significantly reduced when the heme is shielded from solvent fluctuations by the protein material and reduced still further when T approximately < 50 K, as pure dephasing processes due to the protein-solvent phonon bath are frozen out.     

Quantum mechanical/molecular mechanical study of mechanisms of heme degradation by the enzyme heme oxygenase: the strategic function of the water cluster

Heme degradation by heme oxygenase (HO) enzymes is important in maintaining iron homeostasis and prevention of oxidative stress, etc. In response to mechanistic uncertainties, we performed quantum mechanical/molecular mechanical investigations of the heme hydroxylation by HO, in the native route and with the oxygen surrogate donor H2O2. It is demonstrated that H2O2 cannot be deprotonated to yield Fe(III)OOH, and hence the surrogate reaction starts from the FeHOOH complex. The calculations show that, when starting from either Fe(III)OOH or Fe(III)HOOH, the fully concerted mechanism involving O-O bond breakage and O-C(meso) bond formation is highly disfavored. The low-energy mechanism involves a nonsynchronous, effectively concerted pathway, in which the active species undergoes first O-O bond homolysis followed by a barrier-free (small with Fe(III)HOOH) hydroxyl radical attack on the meso position of the porphyrin. During the reaction of Fe(III)HOOH, formation of the Por+*FeIV=O species, compound I, competes with heme hydroxylation, thereby reducing the efficiency of the surrogate route. All these conclusions are in accord with experimental findings (Chu, G. C.; Katakura, K.; Zhang, X.; Yoshida, T.; Ikeda-Saito, M. J. Biol. Chem. 1999, 274, 21319). The study highlights the role of the water cluster in the distal pocket in creating "function" for the enzyme; this cluster affects the O-O cleavage and the O-Cmeso formation, but more so it is responsible for the orientation of the hydroxyl radical and for the observed alpha-meso regioselectivity of hydroxylation (Ortiz de Montellano, P. R. Acc. Chem. Res. 1998, 31, 543). Differences/similarities with P450 and HRP are discussed.     

Compound I of heme oxygenase cannot hydroxylate its heme meso-carbon

Heme oxygenase (HO) catalyzes heme catabolism through three successive oxygenation steps where the substrate heme itself activates O2. It has been thought that the reactive species responsible for the first heme oxygenation, meso-hydroxylation, is the hydroperoxy-ferric heme intermediate (Fe-OOH) rather than an oxo ferryl porphyrin cation radical, so-called compound I. A recent theoretical study (Kamachi, T.; Yoshizawa, K. J. Am. Chem. Soc. 2005, 127, 10686), however, proposed that compound I can oxidize its meso-carbon atom with the assistance of a bridging water molecule. In this communication, we report the first direct observation of compound I of a heme-HO-1 complex, generated by reaction of ferric-HO-1 with m-chloroperbenzoic acid. HO compound I slowly decays to compound II without producing any meso-hydroxylated products. It does react with guaiacol and thioanisole, however. Our findings unambiguously rule out involvement of compound I in the HO catalysis.     

Ligand-field dependence of the excited state dynamics of Hangman bisporphyrin dyad complexes

A new Hangman porphyrin architecture has been developed to interrogate the ligand-field dependence of photoinduced PCET versus excitation energy transfer and intersystem crossing in PZn(II)-PFe(III)-OH dyads (P = porphyrin). In this design, a hanging carboxylic acid group establishes a hydrogen-bonding network to anchor the weak-field OH- ligand in the distal site of the PFe(III)-OH acceptor, whereas the proximal site is left available to accept strong-field imidazole ligands. Thus, controlling the tertiary coordination environment gives access to the first synthetic example of a porphyrin dyad with a biologically relevant weak-field/strong-field configuration of axial ligands at the heme. Transient absorption spectroscopy has been employed to probe the fate of the initial PZn(II)-based S1 excited state, revealing rapid S1 quenching for all dyads in the presence and absence of strong-field imidazole ligands (tau = 6-50 ps). The absence of a (P*+)Zn(II) signal that would complement photoinduced PCET at the PFe(III)-OH subunit (i.e., PFe(III)-OH --> PFe(II)-OH2) shows that excitation energy transfer and intersystem crossing channels dominate the quenching, regardless of whether proximal strong field ligands are present. Moreover, this photophysical assignment is independent of the solvent dielectric constant and whether a phenylene or biphenylene spacer is used to span the two porphyrin subunits. Electronic structure calculations suggest that the structural reorganization attendant to reductive PCET at the high-spin Fe(III)-OH center imposes a severe kinetic cost that can only be alleviated by inducing a low-spin electronic configuration with two strong-field axial ligands.     

Small molecule directed aggregation of a heme peptide on gold: an STM study

Microperoxidase 11 was adsorbed on Au(111) from basic aqueous solutions containing pure heme peptide and co-added stoichiometric amounts of exogenous neutral and ionic ligands. The addition of small molecules to MP11 produced different aggregate structures that were easily differentiated by STM. In the absence of a complexing agent, the MP-11 formed large clusters of metallopeptide molecules on the gold surface. With neutral imidazole in solution the MP11 aggregated into regular elongated structures (nano-épis) on the substrate. When S(2-) is used as coupling agent, single heme peptide molecules are isolated with identifiable porphyrin ring and substructure in the peptide chain.     

Kinetics of histidine deligation from the heme in GuHCl-unfolded Fe(III) cytochrome C studied by a laser-induced pH-jump technique.

We have developed an instrumental setup that uses transient absorption to monitor protein folding/unfolding processes following a laser-induced, ultrafast release of protons from o-nitrobenzaldehyde. The resulting increase in [H(+)], which can be more than 100 microM, is complete within a few nanoseconds. The increase in [H(+)] lowers the pH of the solution from neutrality to approximately 4 at the highest laser pulse energy used. Protein structural rearrangements can be followed by transient absorption, with kinetic monitoring over a broad time range (approximately 10 ns to 500 ms). Using this pH-jump/transient absorption technique, we have examined the dissociation kinetics of non-native axial heme ligands (either histidine His26 or His33) in GuHCl-unfolded Fe(III) cytochrome c (cyt c). Deligation of the non-native ligands following the acidic pH-jump occurs as a biexponential process with different pre-exponential factors. The pre-exponential factors markedly depend on the extent of the pH-jump, as expected from differences in the pK(a) values of His26 and His33. The two lifetimes were found to depend on temperature but were not functions of either the magnitude of the pH-jump or the pre-pulse pH of the solution. The activation energies of the deligation processes support the suggestion that GuHCl-unfolded cyt c structures with non-native histidine axial ligands represent kinetic traps in unfolding.     

Enhanced Reactivities of Iron(IV)‐Oxo Porphyrin π‐Cation Radicals in Oxygenation Reactions by Electron‐Donating Axial Ligands

The proximal axial ligand in heme iron enzymes plays an important role in tuning the reactivities of iron(IV)‐oxo porphyrin π‐cation radicals in oxidation reactions. The present study reports the effects of axial ligands in olefin epoxidation, aromatic hydroxylation, alcohol oxidation, and alkane hydroxylation, by [(tmp)^+. Fe^IV(O)(p‐Y‐PyO)]⁺ ( 1 ‐Y) (tmp=meso‐tetramesitylporphyrin, p‐Y‐PyO=para‐substituted pyridine N‐oxides, and Y=OCH₃, CH₃, H, Cl). In all of the oxidation reactions, the reactivities of 1 ‐Y are found to follow the order 1 ‐OCH₃ > 1 ‐CH₃ > 1 ‐H > 1 ‐Cl; negative Hammett ρ values of ?1.4 to ?2.7 were obtained by plotting the reaction rates against the σ_p values of the substituents of p‐Y‐PyO. These results, as well as previous ones on the effect of anionic nucleophiles, show that iron(IV)‐oxo porphyrin π‐cation radicals bearing electron‐donating axial ligands are more reactive in oxo‐transfer and hydrogen‐atom abstraction reactions. These results are counterintuitive since iron(IV)‐oxo porphyrin π‐cation radicals are electrophilic species. Theoretical calculations of anionic and neutral ligands reproduced the counterintuitive experimental findings and elucidated the root cause of the axial ligand effects. Thus, in the case of anionic ligands, as the ligand becomes a better electron donor, it strengthens the FeO? H bond and thereby enhances its H‐abstraction activity. In addition, it weakens the Fe?O bond and encourages oxo‐transfer reactivity. Both are Bell–Evans–Polanyi effects, however, in a series of neutral ligands like p‐Y‐PyO, there is a relatively weak trend that appears to originate in two‐state reactivity (TSR). This combination of experiment and theory enabled us to elucidate the factors that control the reactivity patterns of iron(IV)‐oxo porphyrin π‐cation radicals in oxidation reactions and to resolve an enigmatic and fundamental problem.     

Low-frequency mode activity of heme: femtosecond coherence spectroscopy of iron porphine halides and nitrophorin

The low-frequency mode activity of metalloporphyrins has been studied for iron porphine-halides (Fe(P)(X), X = Cl, Br) and nitrophorin 4 (NP4) using femtosecond coherence spectroscopy (FCS) in combination with polarized resonance Raman spectroscopy and density functional theory (DFT). It is confirmed that the mode symmetry selection rules for FCS are the same as for Raman scattering and that both Franck-Condon and Jahn-Teller mode activities are observed for Fe(P)(X) under Soret resonance conditions. The DFT-calculated low-frequency (20-400 cm (-1)) modes, and their frequency shifts upon halide substitution, are in good agreement with experimental Raman and coherence data, so that mode assignments can be made. The doming mode is located at approximately 80 cm (-1) for Fe(P)(Cl) and at approximately 60 cm (-1) for Fe(P)(Br). NP4 is also studied with coherence techniques, and the NO-bound species of ferric and ferrous NP4 display a mode at approximately 30-40 cm (-1) that is associated with transient heme doming motion following NO photolysis. The coherence spectra of three ferric derivatives of NP4 with different degrees of heme ruffling distortion are also investigated. We find a mode at approximately 60 cm (-1) whose relative intensity in the coherence spectra depends quadratically on the magnitude of the ruffling distortion. To quantitatively account for this correlation, a new "distortion-induced" Raman enhancement mechanism is presented. This mechanism is unique to low-frequency "soft modes" of the molecular framework that can be distorted by environmental forces. These results demonstrate the potential of FCS as a sensitive probe of dynamic and functionally important nonplanar heme vibrational excitations that are induced by the protein environmental forces or by the chemical reactions in the aqueous phase.     

Ibuprofen induces an allosteric conformational transition in the heme complex of human serum albumin with significant effects on heme ligation

Human serum albumin (HSA), the most prominent protein in blood plasma, is able to bind a wide range of endogenous and exogenous compounds. Among the endogenous ligands, HSA is a significant transporter of heme, the heme-HSA complex being present in blood plasma. Drug binding to heme-HSA affects allosterically the heme affinity for HSA and vice versa. Heme-HSA, heme, and their complexes with ibuprofen have been characterized by electronic absorption, resonance Raman, and electron paramagnetic resonance (EPR) spectroscopy. Comparison of the results for the heme and heme-HSA systems has provided insight into the structural consequences on the heme pocket of ibuprofen binding. The pentacoordinate tyrosine-bound heme coordination of heme-HSA, observed in the absence of ibuprofen, becomes hexacoordinate low spin upon ibuprofen binding, and heme dissociates at increasing drug levels. The electronic absorption spectrum and nu(Fe-CO)/nu(CO) vibrational frequencies of the CO-heme-HSA-ibuprofen complex, together with the observation of a Fe-His Raman mode at 218 cm(-1) upon photolysis of the CO complex and the low spin EPR g values indicate that a His residue is one of the low spin axial ligands, the sixth ligand probably being Tyr161. The only His residue in the vicinity of the heme Fe atom is His146, 9 A distant in the absence of the drug. This indicates that drug binding to heme-HSA results in a significant rearrangement of the heme pocket, implying that the conformational adaptability of HSA involves more than the immediate vicinity of the drug binding site. As a whole, the present spectroscopic investigation supports the notion that HSA could be considered as the prototype of monomeric allosteric proteins.     

Low-energy electron attachment to 5'-thymidine monophosphate: modeling single strand breaks through dissociative electron attachment

Mechanisms of low-energy electron (LEE) attachment and subsequent single-strand break (SSB) formation are investigated by density functional theory treatment of a simple model for DNA, i.e., the nucleotide, 5'-thymidine monophosphate (5'-dTMPH). In the present study, the C5'-O5' bond dissociation due to LEE attachment has been followed along the adiabatic as well as on the vertical (electron attached to the optimized geometry of the neutral molecule) anionic surfaces using B3LYP functional and 6-31G* and 6-31++G** basis sets. Surprisingly, it is found that the PES of C5'-O5' bond dissociation in the anion radicals have approximately the same barrier for both adiabatic and vertical pathways. These results provide support for the hypothesis that transiently bound electrons (shape resonances) to the virtual molecular orbitals of the neutral molecule likely play a key role in the cleavage of the sugar-phosphate C5'-O5' bond in DNA resulting in the direct formation of single strand breaks without significant molecular relaxation. To take into account the solvation effects, we considered the neutral and anion radical of 5'-dTMP surrounded by 5 or 11 water molecules with Na+ as a counterion. These structures were optimized using the B3LYP/6-31G** level of theory. We find the barrier height for adiabatic C5'-O5' bond dissociation of 5'-dTMP anion radical in aqueous environment is so substantially higher than in the gas phase that the adiabatic route will not contribute to DNA strand cleavage in aqueous systems. This result is in agreement with experiment.     

Loss of charged versus neutral heme from gaseous holomyoglobin ions.

The dissociation of holomyoglobin ions ranging in charge state from +10 to +2 has been studied using collisional activation in a quadrupole ion trap. Collisional activation times and amplitudes were varied to investigate the effects of these variables on dissociation of the heme group from the holoprotein. The onset of neutral heme loss occurs at a lower activation amplitude than loss of charged heme. For solutions of ferri-myoglobin, charged heme loss was prominent for +10 to +4 holomyoglobin ions, while neutral heme loss product was found to be dominant for charge states +3 and +2. For any given charge state, activation of holomyoglobin ions from a solution containing primarily ferro-myoglobin yielded significantly more abundant neutral heme loss products than was observed for activation of ions from solutions containing primarily ferri-myoglobin. The relative concentrations of the two oxidation states were shown to be affected by redox chemistry within the nano-electrospray emitter used in this work. Results from a double activation experiment revealed that the precursor ions of a given charge state contained a mixture of two populations, with ferro-myoglobin giving rise to neutral heme loss upon dissociation and ferri-myoglobin yielding charged heme. No evidence for electron transfer upon collisional activation of ferri-myoglobin ions was observed. Furthermore, little or no evidence for electron transfer associated with ion/ion reactions with anions derived from perfluoro-1,3-dimethylcyclohexane was observed. Definitive results could not be drawn for the lowest precursor ion charge states (+3 and +2) due to low dissociation efficiencies.     

1H NMR study of the magnetic properties and electronic structure of the hydroxide complex of substrate-bound heme oxygenase from Neisseria meningitidis: influence of the axial water deprotonation on the distal H-bond network

The substrate and active site residues of the low-spin hydroxide complex of the protohemin complex of Neisseria meningitidis heme oxygenase (NmHO) have been assigned by saturation transfer between the hydroxide and previously characterized aquo complex. The available dipolar shifts allowed the quantitation of both the orientation and anisotropy of the paramagnetic susceptibility tensor. The resulting positive sign, and reduced magnitude of the axial anisotropy relative to the cyanide complex, dictate that the orbital ground state is the conventional "d(pi)" (d(2)(xy)(d(xz), d(yz))(3)); and not the unusual "d(xy)" (d(2)(xz)d(2)(yz)d(xy)) orbital ground state reported for the hydroxide complex of the homologous heme oxygenase (HO) from Pseudomonas aeruginosa (Caignan, G.; Deshmukh, R.; Zeng, Y.; Wilks, A.; Bunce, R. A.; Rivera, M. J. Am. Chem. Soc. 2003, 125, 11842-11852) and proposed as a signature of the HO distal cavity. The conservation of slow labile proton exchange with solvent from pH 7.0 to 10.8 confirms the extraordinary dynamic stability of NmHO complexes. Comparison of the diamagnetic contribution to the labile proton chemical shifts in the aquo and hydroxide complexes reveals strongly conserved bond strengths in the distal H-bond network, with the exception of the distal His53 N(epsilon)(1)H. The iron-ligated water is linked to His53 primarily by a pair of nonligated, ordered water molecules that transmit the conversion of the ligated H-bond donor (H(2)O) to a H-bond acceptor (OH(-)), thereby increasing the H-bond donor strength of the His53 side chain.     

Solution 1H NMR characterization of axial interactions of the proximal and distal His in the cyanomet complexes of the isolated chains and the 65 kDa intact tetramer of human hemoglobin A.

Solution 1H NMR has been used to investigate the axial bonding of the proximal His and the hydrogen-bonding of the distal His to the bound ligand in the isolated chains as well as the subunits of intact, tetrameric, cyanomet human hemoglobin A. The complete proximal His, including all ring protons necessary to monitor bonding in each subunit, could be definitively assigned by 1D/2D methods despite the large size (approximately 65 kDa) and severe relaxation (to T(1) approximately 3 ms, line width approximately 1.5 kHz) of two of the protons. The complete distal His E7 ring was assigned in the alpha-chain and alpha-subunit of HbA, and the dipolar shifts and relaxation were analyzed to reveal a disposition intermediate between the positions adopted in HbCO and HbO2 that is optimal for forming a hydrogen bond with bound cyanide. The lability of the alpha-subunit His E7 NepsilonH is found to be similar to that in sperm whale cyanomet myoglobin. The orientation of the distal His E7 in the beta-subunit is found to be consistent with that seen in either HbCO or HbO2. While the His E7 labile NepsilonH proton signal could not be detected in either the beta-chain or subunit, it is concluded that this more likely reflects increased lability over that of the alpha-subunit, and not the absence of a hydrogen bond to the bound ligand. Analysis of the heme mean methyl hyperfine shift, which has been shown to be very sensitive to the presence of distal hydrogen bonds to bound cyanide (Nguyen, B. D.; Xia, Z.; Cutruzzolá, F.; Travaglini Allocatelli, C.; Brancaccio, A.; Brunori, M.; La Mar, G. N. J. Biol. Chem. 2000, 275, 742-751), directly supports the presence of a distal His E7 hydrogen bond to cyanide in the beta-chain and beta-subunit which is weaker than the same hydrogen bond in the alpha-subunit. The potential for the proximal His hyperfine shifts in serving as indicators of axial strain in the allosteric transition of HbA is discussed.     

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.     

Ab initio molecular dynamics of heme in cytochrome c

Ab initio molecular dynamics (AIMD) calculations, based on the Car-Parrinello method, have been carried out for three models of heme c that is present in cytochrome c. Both the reduced (Fe(II)) and oxidized (Fe(III)) forms have been analyzed. The simplest models (1R and 1O, respectively) consist of a unsubstituted porphyrin (with no side chains) and two axially coordinated imidazole and ethylmethylthioether ligands. Density functional theory optimizations of these models confirm the basic electronic features and are the starting point for building more complex derivatives. AIMD simulations were performed after reaching the thermal stability at T = 300 K. The evolution of the Fe-L(ax) bond strengths is examined together with the relative rotations of the imidazole and methionine about the axial vector, which appear rather independent from each other. The next models (2R and 2O) contain side chains at the heme to better simulate the actual active site. It is observed that two adjacent propionate groups induce some important effects. The axial Fe-Sdelta bond is only weakened in 2R but is definitely cleaved in the oxidized species 2O. Also the mobility of the Im ligand seems to be reduced by the formation of a strong hydrogen bond that involves the Im Ndelta1-Hdelta1 bond and one carboxylate group. In 2O the interaction becomes so strong that a proton transfer occurs and the propionic acid is formed. Finally, the models 3 include a free N-methyl-acetamide molecule to mimic a portion of the protein backbone. This influences the orientation of carboxylate groups and limits the amount of their hydrogen bonding with the Im ligand. Residual electrostatic interactions are maintained, which are still able to modulate the dissociation of the methionine from the heme.     

A proton-shuttle mechanism mediated by the porphyrin in benzene hydroxylation by cytochrome p450 enzymes

Benzene hydroxylation is a fundamental process in chemical catalysis. In nature, this reaction is catalyzed by the enzyme cytochrome P450 via oxygen transfer in a still debated mechanism of considerable complexity. The paper uses hybrid density functional calculations to elucidate the mechanisms by which benzene is converted to phenol, benzene oxide, and ketone, by the active species of the enzyme, the high-valent iron-oxo porphyrin species. The effects of the protein polarity and hydrogen-bonding donation to the active species are mimicked, as before (Ogliaro, F.; Cohen, S.; de Visser, S. P.; Shaik, S. J. Am. Chem. Soc. 2000, 122, 12892-12893). It is verified that the reaction does not proceed either by hydrogen abstraction or by initial electron transfer (Ortiz de Montellano, P. R. In Cytochrome P450: Structure, Mechanism and Biochemistry, 2nd ed.; Ortiz de Montellano, P. R., Ed.; Plenum Press: New York, 1995; Chapter 8, pp 245-303). In accord with the latest experimental conclusions, the theoretical calculations show that the reactivity is an interplay of electrophilic and radicalar pathways, which involve an initial attack on the pi-system of the benzene to produce sigma-complexes (Korzekwa, K. R.; Swinney, D. C.; Trager, W. T. Biochemistry 1989, 28, 9019-9027). The dominant reaction channel is electrophilic and proceeds via the cationic sigma-complex,( 2)3, that involves an internal ion pair made from a cationic benzene moiety and an anionic iron porphyrin. The minor channel proceeds by intermediacy of the radical sigma-complex, (2)2, in which the benzene moiety is radicalar and the iron-porphyrin moiety is neutral. Ring closure in these intermediates produces the benzene oxide product ((2)4), which does not rearrange to phenol ((2)7) or cyclohexenone ((2)6). While such a rearrangement can occur post-enzymatically under physiological conditions by acid catalysis, the computations reveal a novel mechanism whereby the active species of the enzyme catalyzes directly the production of phenol and cyclohexenone. This enzymatic mechanism involves proton shuttles mediated by the porphyrin ring through the N-protonated intermediate, (2)5, which relays the proton either to the oxygen atom to form phenol ((2)7) or to the ortho-carbon atom to produce cyclohexenone product ((2)6). The formation of the phenol via this proton-shuttle mechanism will be competitive with the nonenzymatic conversion of benzene oxide to phenol by external acid catalysis. With the assumption that (2)5 is not fully thermalized, this novel mechanism would account also for the observation that there is a partial skeletal retention of the original hydrogen of the activated C-H bond, due to migration of the hydrogen from the site of hydroxylation to the adjacent carbon (so-called "NIH shift" (Jerina, D. M.; Daly, J. W. Science 1974, 185, 573-582)). Thus, in general, the computationally discovered mechanism of a porphyrin proton shuttle suggests thatthere is an enzymatic pathway that converts benzene directly to a phenol and ketone, in addition to nonenzymatic production of these species by conversion of arene oxide to phenol and ketone. The potential generality of protonated porphyrin intermediates in P450 chemistry is discussed in the light of the H/D exchange observed during some olefin epoxidation reactions (Groves, J. T.; Avaria-Neisser, G. E.; Fish, K. M.; Imachi, M.; Kuczkowski, R. J. Am. Chem. Soc. 1986, 108, 3837-3838) and the general observation of heme alkylation products (Kunze, K. L.; Mangold, B. L. K.; Wheeler, C.; Beilan, H. S.; Ortiz de Montellano, P. R. J. Biol. Chem. 1983, 258, 4202-4207). The competition, similarities, and differences between benzene oxidation viz. olefin epoxidation and alkanyl C-H hydroxylation are discussed, and comparison is made with relevant experimental and computational data. The dominance of low-spin reactivity in benzene hydroxylation viz. two-state reactivity (Shaik, S.; de Visser, S. P.; Ogliaro, F.; Schwarz, H.; Schr?der, D. Curr. Opin. Chem. Biol. 2002, 6, 556-567) in olefin epoxidation and alkane hydroxylation is traced to the loss of benzene resonance energy during the bond activation step.     

Mutual control of axial and equatorial ligands: model studies with [Ni]-bacteriochlorophyll-a

Modification of the metal's electronic environment by ligand association and dissociation in metalloenzymes is considered cardinal to their catalytic activity. We have recently presented a novel system that utilizes the bacteriochlorophyll (BChl) macrocycle as a ligand and reporter. This system allows for charge mobilization in the equatorial plane and experimental estimate of changes in the electronic charge density around the metal with no modification of the metal's chemical environment. The unique spectroscopy, electrochemistry and coordination chemistry of [Ni]-bacteriochlorophyll ([Ni]-BChl) enable us to follow directly fine details and steps involved in the function of the metal redox center. This approach is utilized here whereby electro-chemical reduction of [Ni]-BChl to the monoanion [Ni]-BChl(-) results in reversible dissociation of biologically relevant axial ligands. Similar ligand dissociation was previously detected upon photoexcitation of [Ni]-BChl (Musewald, C.; Hartwich, G.; Lossau, H.; Gilch, P.; Pollinger-Dammer, F.; Scheer, H.; Michel-Beyerle, M. E. J. Phys. Chem. B 1999, 103, 7055-7060 and Noy, D.; Yerushalmi, R.; Brumfeld, V.; Ashur, I.; Baldridge, K. K.; Scheer, H.; Scherz, A. J. Am. Chem. Soc. 2000, 122, 3937-3944). The electrochemical measurements and quantum mechanical (QM) calculations performed here for the neutral, singly reduced, monoligated, and singly reduced, monoligated [Ni]-BChl suggest the following: (a) Electroreduction, although resulting in a pi anion [Ni]-BChl(-) radical, causes electron density migration to the [Ni]-BChl core. (b) Reduction of nonligated [Ni]-BChl does not change the macrocycle conformation, whereas axial ligation results in a dramatic expansion of the metal core and a flattening of the highly ruffled macrocycle conformation. (c) In both the monoanion and singly excited [Ni]-BChl ([Ni]-BChl*), the frontier singly occupied molecular orbital (SOMO) has a small but nonnegligible metal character. Finally, (d) computationally, we found that a reduction of [Ni]-BChl*imidazole results in a weaker metal-axial ligand bond. Yet, it remains weakly bound in the gas phase. The experimentally observed ligand dissociation is accounted for computationally when solvation is considered. On the basis of the experimental observations and QM calculations, we propose a mechanism whereby alterations in the equatorial pi system and modulation of sigma bonding between the axial ligands and the metal core are mutually correlated. Such a mechanism highlights the dynamic role of axial ligands in regulating the activity of metal centers such as factor F430 (F430), a nickel-based coenzyme that is essential in methanogenic archea.