The catalytic cycle of catechol oxidase

Hybrid density functional theory with the B3LYP functional has been used to investigate the catalytic mechanism of catechol oxidase. Catechol oxidase belongs to a class of enzymes that has a copper dimer with histidine ligands at the active site. Another member of this class is tyrosinase, which has been studied by similar methods previously. An important advantage for the present study compared to the one for tyrosinase is that X-ray crystal structures exist for catechol oxidase. The most critical step in the mechanism for catechol oxidase is where the peroxide O–O bond is cleaved. In the suggested mechanism this cleavage occurs in concert with a proton transfer from the substrate. Shortly after the transition state is passed there is another proton transfer from the substrate, which completes the formation of a water molecule. An important feature of the mechanism, like the one for tyrosinase, is that no proton transfers to or from residues outside the metal complex are needed. The calculated energetics is in reasonable agreement with experiments. Comparisons are made to other similar enzymes studied previously.     

Mechanism for intein C-terminal cleavage: a proposal from quantum mechanical calculations

Inteins are autocatalytic protein cleavage and splicing elements. A cysteine to alanine mutation at the N-terminal of inteins inhibits splicing and isolates the C-terminal cleavage reaction. Experiments indicate an enhanced C-terminal cleavage reaction rate upon decreasing the solution pH for the cleavage mutant, which cannot be explained by the existing mechanistic framework. We use intein crystal structure data and the information about conserved amino acids to perform semiempirical PM3 calculations followed by high-level density functional theory calculations in both gas phase and implicit solvent environments. Based on these calculations, we propose a detailed “low pH” mechanism for intein C-terminal cleavage. Water plays an important role in the proposed reaction mechanism, acting as an acid as well as a base. The protonation of the scissile peptide bond nitrogen by a hydronium ion is an important first step in the reaction. That step is followed by the attack of the C-terminal asparagine side chain on its carbonyl carbon, causing succinimide formation and simultaneous peptide bond cleavage. The computed reaction energy barrier in the gas phase is ~33 kcal/mol and reduces to ~25 kcal/mol in solution, close to the 21 kcal/mol experimentally observed at pH 6.0. This mechanism is consistent with the observed increase in C-terminal cleavage activity at low pH for the cleavage mutant of the Mycobacterium tuberculosis RecA mini-intein.     

A density functional investigation of the extradiol cleavage mechanism in non-heme iron catechol dioxygenases

The mechanism for extradiol cleavage in non-heme iron catechol dioxygenase was modelled theoretically via density functional theory. Based on the Fe(II)-His,His,Glu motif observed in enzymes, an active site model complex, [Fe(acetate)(imidazole)(2)(catecholate)(O(2))](-), was optimized for states with six, four and two unpaired electrons (U6, U4 and U2, respectively). The transfer of the terminal atom of the coordinated dioxygen leading to "ferryl" Fe=O intermediates spontaneously generates an extradiol epoxide. The computed barriers range from 19 kcal mol(-1) on the U6 surface to approximately 25 kcal mol(-1) on the U4 surface, with overall reaction energies of +11.6, 6.3 and 7.1 kcal mol(-1) for U6, U4 and U2, respectively. The calculations for a protonated process reveal the terminal oxygen of O(2) to be the thermodynamically favoured site but subsequent oxygen transfer to the catechol has a barrier of approximately 30-40 kcal mol(-1), depending on the spin state. Instead, protonating the acetate group gives a slightly higher energy species but a subsequent barrier on the U4 surface of only 7 kcal mol(-1) relative to the hydroperoxide complex. The overall exoergicity increases to 13 kcal mol(-1). The favoured proton-assisted pathway does not involve significant radical character and has features reminiscent of a Criegee rearrangement which involves the participation of the aromatic ring pi-orbitals in the formation of the new carbon-oxygen bond. The subsequent collapse of the epoxide, attack by the coordinated hydroxide and final product formation proceeds with an overall exoergicity of approximately 75 kcal mol(-1) on the U4 surface.     

Base release in nucleosides induced by low-energy electrons: a DFT study

Low-energy electrons are known to induce strand breaks and base damage in DNA and RNA through fragmentation of molecular bonding. Recently the glycosidic bond cleavage of nucleosides by low-energy electrons has been reported. These experimental results call for a theoretical investigation of the strength of the C(1)'-N link in nucleosides (dA, dC and dT) between the base and deoxyribose before and after electron attachment. Through density functional theory (DFT) calculations, we compare the C(1)'-N bond strength, i.e., the bond dissociation energy of the neutral and its anionic radical, and find that an excess electron effectively weakens the C(1)'- N bond strength in nucleosides by 61-75 kcal/mol in the gas phase and 76-83 kcal/mol in the solvated environment. As a result, electron-induced fragmentation of the C(1)'-N bond in the gas phase is exergonic for dA (DeltaG=-14 kcal/mol) and for dT (DeltaG=-6 kcal/mol) and is endergonic (DeltaG=+1 kcal/ mol) only for dC. In the gas phase all the anionic nucleosides are found to be in valence states. Solvation is found to increase the exergonic nature by an additional 20 kcal, making the fragmentation both exothermic and exergonic for all nucleoside anion radicals. Thus C(1)'-N bond breaking in nucleoside anion radicals is found to be thermodynamically favorable both in the gas phase and under solvation. The activation barrier for the C(1)'-N bond breaking process was found to be about 20 kcal/mol in every case examined, suggesting that a 1 eV electron would induce spontaneous cleavage of the bond and that stabilized anion radicals on the DNA strand would undergo base release at only a modest rate at room temperature. These results suggest that base release from nucleosides and DNA is an expected consequence of low-energy electron-induced damage but that the high barrier would inhibit this process in the stable anion radicals.     

Application of molecular dynamics and free energy perturbation methods to metalloporphyrin-ligand systems II: CO and dioxygen binding to myoglobin.

The protein contribution to the relative binding affinity of the ligands CO and O2 toward myoglobin (Mb) has been simulated using free energy perturbation calculations. The tautomers of the His E7 residue are different for the oxymyoglobin (MbO2) and carboxymyoglobin (MbCO) systems. This was modeled by performing two-step calculations that mutate the ligand and mutate the His E7 tautomers in separate steps. Differences in hydrogen bonding to the O2 and CO ligands were incorporated into the model. The O2 complex was calculated to be 2-3 kcal/mol more stable than the corresponding CO complex when compared to the same difference in an isolated heme control. This value agrees well with the experimental value of 2.0 kcal/mol. In qualitative agreement with experiments, the Fe-C-O bond is found to be bent (theta = 159.8 degrees) with a small tilt (theta = 6.2 degrees). The contributions made by each of the 29 residues--within the 9.0-A radius of the iron atom--to the free energy difference are separated into van der Waals and electrostatic contributions; the latter contributions are dominant. Aside from the proximal histidine and the heme group, the residues having the largest difference in free energy in mutating MbO2-->MbCO are His E7, Phe CD1, Phe CD4, Val E11, and Thr E10.     

Density functional models of the mechanism for decarboxylation in orotidine decarboxylase

The mechanism of orotidine 5-monophosphate decarboxylase (ODCase) has been modeled using hybrid Density Functional Theory (B3LYP functional). The main goal of the present study was to investigate if much larger quantum chemical models of the active site than previously used could shed new light on the mechanism. The models used include the five conserved amino acids expected to be the most important ones for catalysis. One result of this model is that a mechanism involving a direct cleavage of the C–C bond followed by a protonation of C6 by Lys93 appears unlikely, with a barrier for decarboxylation 20 kcal mol^–1 too high. Additional effects like electrostatic stress and ground-state destabilization have been estimated to have only a minor influence on the reaction barrier. The conclusion from the calculations is that the negative charge developing on the substrate during decarboxylation must be stabilized by a protonation of the carbonyl O2 of the substrate. For this mechanism, the addition of the catalytic amino acids decreases the reaction barrier by 25 kcal mol^–1, but full agreement with experimental results has still not been reached. Further modifications of this mechanism are discussed. Electronic supplementary material to this paper can be obtained by using the Springer LINK server located a http://dx.doi.org/10.1007/s00894-002-0080-2.Electronic Supplementary Material available.     

An integrated study of tyrosinase inhibition by rutin: progress using a computational simulation

Tyrosinase inhibition studies have recently gained the attention of researchers due to their potential application values. We simulated docking (binding energies for AutoDock Vina: -9.1 kcal/mol) and performed a molecular dynamics simulation to verify docking results between tyrosinase and rutin. The docking results suggest that rutin mostly interacts with histidine residues located in the active site. A 10 ns molecular dynamics simulation showed that one copper ion at the tyrosinase active site was responsible for the interaction with rutin. Kinetic analyses showed that rutin-mediated inactivation followed a first-order reaction and mono- and biphasic rate constants occurred with rutin. The inhibition was a typical competitive type with K(i) = 1.10±0.25 mM. Measurements of intrinsic and ANS-binding fluorescences showed that rutin showed a relatively strong binding affinity for tyrosinase and one possible binding site that could be a copper was detected accompanying with a hydrophobic exposure of tyrosinase. Cell viability testing with rutin in HaCaT keratinocytes showed that no toxic effects were produced. Taken together, rutin has the potential to be a potent anti-pigment agent. The strategy of predicting tyrosinase inhibition based on hydroxyl group number and computational simulation may prove useful for the screening of potential tyrosinase inhibitors.     

Mechanism of peptide hydrolysis by co-catalytic metal centers containing leucine aminopeptidase enzyme: a DFT approach

In this density functional theory study, reaction mechanisms of a co-catalytic binuclear metal center (Zn1–Zn2) containing enzyme leucine aminopeptidase for two different metal bridging nucleophiles (H₂O and –OH) have been investigated. In addition, the effects of the substrate (l-leucine-p-nitroanilide → l-leucyl-p-anisidine) and metal (Zn1 → Mg and Zn2 → Co, i.e., Mg1–Zn2 and Mg1–Co2 variants) substitutions on the energetics of the mechanism have been investigated. The general acid/base mechanism utilizing a bicarbonate ion followed by this enzyme is divided into two steps: (1) the formation of the gem-diolate intermediate, and (2) the cleavage of the peptide bond. With the computed barrier of 17.8 kcal/mol, the mechanism utilizing a hydroxyl nucleophile was found to be in excellent agreement with the experimentally measured barrier of 18.7 kcal/mol. The rate-limiting step for reaction with l-leucine-p-nitroanilide is the cleavage of the peptide bond with a barrier of 17.8 kcal/mol. However, for l-leucyl-p-anisidine all steps of the mechanism were found to occur with similar barriers (18.0–19.0 kcal/mol). For the metallovariants, cleavage of the peptide bond occurs in the rate-limiting step with barriers of 17.8, 18.0, and 24.2 kcal/mol for the Zn1–Zn2, Mg1–Zn2, and Mg1–Co2 enzymes, respectively. The nature of the metal ion was found to affect only the creation of the gem-diolate intermediate, and after that all three enzymes follow essentially the same energetics. The results reported in this study have elucidated specific roles of both metal centers, the nucleophile, indirect ligands, and substrates in the catalytic functioning of this important class of binuclear metallopeptidases.     

An Integrated Study of Tyrosinase Inhibition by Rutin: Progress using a Computational Simulation

Abstract

Tyrosinase inhibition studies have recently gained the attention of researchers due to their potential application values. We simulated docking (binding energies for AutoDock Vina: ?9.1 kcal/mol) and performed a molecular dynamics simulation to verify docking results between tyrosinase and rutin. The docking results suggest that rutin mostly interacts with histidine residues located in the active site. A 10 ns molecular dynamics simulation showed that one copper ion at the tyrosinase active site was responsible for the interaction with rutin. Kinetic analyses showed that rutin-mediated inactivation followed a first-order reaction and mono- and biphasic rate constants occurred with rutin. The inhibition was a typical competitive type with K_i = 1.10 ± 0.25 mM. Measurements of intrinsic and ANS-binding fluorescences showed that rutin showed a relatively strong binding affinity for tyrosinase and one possible binding site that could be a copper was detected accompanying with a hydrophobic exposure of tyrosinase. Cell viability testing with rutin in HaCaT keratinocytes showed that no toxic effects were produced. Taken together, rutin has the potential to be a potent antipigment agent. The strategy of predicting tyrosinase inhibition based on hydroxyl group number and computational simulation may prove useful for the screening of potential tyrosinase inhibitors.     

A comparison of the reaction mechanisms of iron- and manganese-containing 2,3-HPCD: an important spin transition for manganese

Homoprotocatechuate (HPCA) dioxygenases are enzymes that take part in the catabolism of aromatic compounds in the environment. They use molecular oxygen to perform the ring cleavage of ortho-dihydroxylated aromatic compounds. A theoretical investigation of the catalytic cycle for HPCA 2,3-dioxygenase isolated from Brevibacterium fuscum (Bf 2,3-HPCD) was performed using hybrid DFT with the B3LYP functional, and a reaction mechanism is suggested. Models of different sizes were built from the crystal structure of the enzyme and were used in the search for intermediates and transition states. It was found that the enzyme follows a reaction pathway similar to that for other non-heme iron dioxygenases, and for the manganese-dependent analog MndD, although with different energetics. The computational results suggest that the rate-limiting step for the whole reaction of Bf 2,3-HPCD is the protonation of the activated oxygen, with an energy barrier of 17.4 kcal/mol, in good agreement with the experimental value of 16 kcal/mol obtained from the overall rate of the reaction. Surprisingly, a very low barrier was found for the O-O bond cleavage step, 11.3 kcal/mol, as compared to 21.8 kcal/mol for MndD (sextet spin state). This result motivated additional studies of the manganese-dependent enzyme. Different spin coupling between the unpaired electrons on the metal and on the evolving substrate radical was observed for the two enzymes, and therefore the quartet spin state potential energy surface of the MndD reaction was studied. The calculations show a crossing between the sextet and the quartet surfaces, and it was concluded that a spin transition occurs and determines a barrier of 14.4 kcal/mol for the O-O bond cleavage, which is found to be the rate-limiting step in MndD. Thus the two 83% identical enzymes, using different metal ions as co-factors, were found to have similar activation energies (in agreement with experiment), but different rate-limiting steps.     

O-O bond splitting mechanism in cytochrome oxidase

Hybrid density functional theory (DFT) calculations have been used to investigate different mechanisms for O–O bond splitting in cytochrome oxidase. It is shown that the requirement for a low activation barrier for the O–O bond splitting is that two protons, apart from the tyrosine hydroxyl proton, are available at the binuclear center. A mechanism is suggested for the transformation from a species with a molecularly coordinated O₂, to an O–O cleaved species with an oxo-ferryl group. The mechanism has a calculated activation barrier in reasonable agreement with experimental estimates, and the overall reaction is close to thermoneutral, in line with the requirement that the energy wasted as heat should be minimized. The rate limiting step in the mechanism occurs at the initial Fe–O₂ intermediate, consistent with experimental observations that the decay of the oxy intermediate parallels the increase of the oxo product. The formation of a radical at the cross-linked tyrosine–histidine structure is a possible source for one of the electrons required in the bond cleavage process. Possible sources for the two protons are discussed, including a suggested key role for the hydroxyl group on the farnesyl side chain of heme a₃.     

Theoretical study of the discrimination between O(2) and CO by myoglobin

Combined quantum chemical and molecular mechanics geometry optimisations have been performed on myoglobin without or with O(2) or CO bound to the haem group. The results show that the distal histidine residue is protonated on the N(epsilon 2) atom and forms a hydrogen bond to the haem ligand both in the O(2) and the CO complexes. We have also re-refined the crystal structure of CO[bond]myoglobin by a combined quantum chemical and crystallographic refinement. Thereby, we probably obtain the most accurate available structure of the active site of this complex, showing a Fe[bond]C[bond]O angle of 171 degrees, and Fe[bond]C and C[bond]O bond lengths of 170-171 and 116-117 pm. The resulting structures have been used to calculate the strength of the hydrogen bond between the distal histidine residue and O(2) or CO in the protein. This amounts to 31-33 kJ/mol for O(2) and 2-3 kJ/mol for CO. The difference in hydrogen-bond strength is 21-22 kJ/mol when corrected for entropy effects. This is slightly larger than the observed discrimination between O(2) or CO by myoglobin, 17 kJ/mol. We have also estimated the strain of the active site inside the protein. It is 2-4 kJ/mol larger for the O(2) complex than for the CO complex, independent of which crystal structure the calculations are based on. Together, these results clearly show that myoglobin discriminates between O(2) and CO mainly by electrostatic interactions, rather than by steric strain.     

The effect of histidine residue modification on tyrosinase activity and conformation: inhibition kinetics and computational prediction

We found that the histidine chemical modification of tyrosinase conspicuously inactivated enzyme activity. The substrate reactions with diethylpyridinecarbamate showed slow-binding inhibition kinetics (K(I) = 0.24 +/- 0.03 mM). Bromoacetate, as another histidine modifier, was also applied in order to study inhibition kinetics. The bromoacetate directly induced the exposures of hydrophobic surfaces following by complete inactivation via ligand binding. For further insights, we predicted the 3D structure of tyrosinase and simulated the docking between tyrosinase and diethylpyridinecarbamate. The docking simulation was shown to the significant binding energy scores (-3.77 kcal/mol by AutoDock4 and -25.26 kcal/mol by Dock6). The computational prediction was informative to elucidate the role of free histidine residues at the active site, which are related to substrate accessibility during tyrosinase catalysis.     

Mixed-Type Inhibition of Tyrosinase from Agaricus bisporus by Terephthalic Acid: Computational Simulations and Kinetics

Tyrosinase inhibition studies are needed due to the agricultural and medicinal applications. For probing effective inhibitors of tyrosinase, a combination of computational prediction and enzymatic assay via kinetics were important. We predicted the 3D structure of tyrosinase from Agaricus bisporus, used a docking algorithm to simulate binding between tyrosinase and terephthalic acid (TPA) and studied the reversible inhibition of tyrosinase by TPA. Simulation was successful (binding energies for Autodock4 = -1.54 and Fred2.0 = -3.19 kcal/mol), suggesting that TPA interacts with histidine residues that are known to bind with copper ions at the active site. TPA inhibited tyrosinase in a mixed-type manner with a K ( i ) = 11.01 ± 2.12 mM. Measurements of intrinsic and ANS-binding fluorescences showed that TPA induced no changes in tertiary structure. The present study suggested that the strategy of predicting tyrosinase inhibition based on hydroxyl groups and orientation may prove useful for screening of potential tyrosinase inhibitors.     

Peroxynitrite reductase activity of selenoprotein glutathione peroxidase: a computational study

The peroxynitrite reductase activity of selenoprotein glutathione peroxidase (GPx) has been investigated using density functional theory calculations for peroxynitrite/peroxynitrous acid (ONOO-/ONOOH) substrates through two different "oxidation" and "nitration" pathways. In the oxidation pathway for ONOO-, the oxidation of GPx and the subsequent formation of the selenenic acid (E-Se-OH) occur through a concerted mechanism with an energy barrier of 4.7 (3.7) kcal/mol, which is in good agreement with the computed value of 7.1 kcal/mol for the drug ebselen and the experimentally measured barrier of 8.8 kcal/mol for both ebselen and GPx. For ONOOH, the formation of the E-Se-OH prefers a stepwise mechanism with an overall barrier of 6.9 (11.3) kcal/mol, which is 10.2 (11.2) kcal/mol lower than that for hydrogen peroxide (H2O2), indicating that ONOOH is a more efficient substrate for GPx oxidation. It has been demonstrated that the active site Gln83 residue plays a critical role during the oxidation process, which is consistent with the experimental suggestions. The nitration of GPx by ONOOH produces a nitro (E-Se-NO2) product via either of two different mechanisms, isomerization and direct, having almost the same barrier heights. A comparison between the rate-determining barriers of the oxidation and nitration pathways suggests that the oxidation of GPx by ONOOH is more preferable than its nitration. It was also shown that the rate-determining barriers remain the same, 21.5 (25.5) kcal/mol, in the peroxynitrite reductase and peroxidase activities of GPx.     

Effect of hesperetin on tyrosinase: inhibition kinetics integrated computational simulation study

Tyrosinase inhibitors have potential applications in medicine, cosmetics and agriculture to prevent hyperpigmentation or browning effects. Some of the flavonoids mostly found in herbal plants and fruits are revealed as tyrosinase inhibitors. We studied the inhibitory effects of one such flavonoid, hesperetin, on mushroom tyrosinase using inhibition kinetics and computational simulation. Hesperetin reversibly inhibited tyrosinase in a competitive manner with K_i = 4.03 ± 0.26 mM. Measurements of ANS-binding fluorescence showed that hesperetin induced the hydrophobic disruption of tyrosinase. For further insight, we used the docking algorithms to simulate binding between tyrosinase and hesperetin. Simulation was successful (binding energies for Dock6.3: −34.41 kcal/mol and for AutoDock4.2: −5.67 kcal/mol) and showed that a copper ion coordinating with 3 histidine residues (HIS61, HIS85, and HIS259) within the active site pocket was chelated via hesperetin binding. Our study provides insight into the inhibition of tyrosinase in response to flavonoids. A combination of inhibition kinetics and computational prediction may facilitate the identification of potential natural tyrosinase inhibitors such as flavonoids and the prediction of their inhibitory mechanisms.     

Inhibition of tyrosinase by gastrodin: An integrated kinetic-computational simulation analysis

We studied the inhibitory effect of gastrodin on tyrosinase using inhibition kinetics and computational simulation. Gastrodin reversibly inhibited tyrosinase in a mixed-type manner with K_i = 123.8 ± 20.2 mM. Time-interval kinetics revealed the inhibition to be a first-order process with mono- and bi-phasic components. Using AutoDock Vina, we calculated a binding energy of ?6.3 kcal/mol for gastrodin and tyrosinase, and we performed a molecular dynamics simulation of the tyrosinase–gastrodin interaction. The simulation results suggested that gastrodin interacts primarily with histidine residues in the active site. A 10-ns molecular dynamics simulation showed that one copper ion in the tyrosinase active site was responsible for the interaction with gastrodin. Our study provides insight into the inhibition of tyrosinase by the hydroxyl groups of gastrodin. A combination of inhibition kinetics and computational calculations may help to confirm the inhibitory action of gastrodin on tyrosinase and define the mechanisms of inhibition.     

Refined crystal structure of ascorbate oxidase at 1.9 A resolution.

The crystal structure of the fully oxidized form of ascorbate oxidase (EC 1.10.3.3) from Zucchini has been refined at 1.90 A (1 A = 0.1 nm) resolution, using an energy-restrained least-squares refinement procedure. The refined model, which includes 8764 protein atoms, 9 copper atoms and 970 solvent molecules, has a crystallographic R-factor of 20.3% for 85,252 reflections between 8 and 1.90 A resolution. The root-mean-square deviation in bond lengths and bond angles from ideal values is 0.011 A and 2.99 degrees, respectively. The subunits of 552 residues (70,000 Mr) are arranged as tetramers with D2 symmetry. One of the dyads is realized by the crystallographic axis parallel to the c-axis giving one dimer in the asymmetric unit. The dimer related about this crystallographic axis is suggested as the dimer present in solution. Asn92 is the attachment site for one of the two N-linked sugar moieties, which has defined electron density for the N-linked N-acetyl-glucosamine ring. Each subunit is built up by three domains arranged sequentially on the polypeptide chain and tightly associated in space. The folding of all three domains is of a similar beta-barrel type and related to plastocyanin and azurin. An analysis of intra- and intertetramer hydrogen bond and van der Waals interactions is presented. Each subunit has four copper atoms bound as mononuclear and trinuclear species. The mononuclear copper has two histidine, a cysteine and a methionine ligand and represents the type-1 copper. It is located in domain 3. The bond lengths of the type-1 copper centre are comparable to the values for oxidized plastocyanin. The trinuclear cluster has eight histidine ligands symmetrically supplied from domain 1 and 3. It may be subdivided into a pair of copper atoms with histidine ligands whose ligating N-atoms (5 NE2 atoms and one ND1 atom) are arranged trigonal prismatic. The pair is the putative type-3 copper. The remaining copper has two histidine ligands and is the putative spectroscopic type-2 copper. Two oxygen atoms are bound to the trinuclear species as OH- or O2- and bridging the putative type-3 copper pair and as OH- or H2O bound to the putative type-2 copper trans to the copper pair. The bond lengths within the trinuclear copper site are similar to comparable binuclear model compounds. The putative binding site for the reducing substrate is close to the type-1 copper.(ABSTRACT TRUNCATED AT 400 WORDS)     

Molecular modeling of the mechanochemical triggering mechanism for catalysis of carbon-cobalt bond homolysis in coenzyme B12

The possible contributions of the mechanochemical triggering effect to the enzymatic activation of the carbon-cobalt bond of coenzyme B12 (5'-deoxyadenosylcobalamin, AdoCbl) for homolytic cleavage have been studied by molecular modeling and semiempirical molecular orbital calculations. Classically, this effect has envisioned enzymatic compression of the axial Co-N bond in the ground state to cause upward folding of the corrin ring and subsequent sterically induced distortion of the Co-C bond leading to its destabilization. The models of this process show that in both methylcobalamin (CH3Cbl) and AdoCbl, compression of the axial Co-N bond does engender upward folding of the corrin ring, and that the extent of such upward folding is smaller in an analog in which the normal 5,6-dimethylbenzimidazole axial ligand is replaced by the sterically smaller ligand, imidazole (CH3(lm)Cbl and Ado(lm)Cbl). Furthermore, in AdoCbl, this upward folding of the corrin is accompanied by increases in the carbon-cobalt bond length and in the Co-C-C bond angle (which are also less pronounced in Ado(Im)Cbl), and which indicate that the Co-C bond is indeed destabilized by this mechanism. However, these effects on the Co-C bond are small, and destabilization of this bond by this mechanism is unlikely to contribute more than ca. 3 kcal mol(-1) towards the enzymatic catalysis of Co-C bond homolysis, far short of the observed ca. 14 kcal mol(-1). A second version of mechanochemical triggering, in which compression of the axial Co-N bond in the transition state for Co-C bond homolysis stabilizes the transition state by increased Co-N orbital overlap, has also been investigated. Stretching the Co-C bond to simulate the approach to the transition state was found to result in an upward folding of the corrin ring, a slight decrease in the axial Co-N bond length, a slight displacement of the metal atom from the plane of the equatorial nitrogens towards the "lower" axial ligand, and a decrease in strain energy amounting to about 8 kcal mol(-1) for both AdoCbl and Ado(Im)Cbl. In such modeled transition states, compression of the axial Co-N bond to just below 2.0 A (the distance subsequently found to provide maximal stabilization of the transition state by increased orbital overlap) required about 4 kcal mol(-1) for AdoCbl, and about 2.5 kcal mol(-1) for Ado(Im)Cbl. ZINDO/1 calculations on slightly simplified structures showed that maximal electronic stabilization of the transition state by about 10 kcal mol(-1) occurred at an axial Co-N bond distance of 1.96 A for both AdoCbl and Ado(Im)Cbl. The net result is that this type of transition state mechanochemical triggering can provide 14 kcal mol(-1) of transition state stabilization for AdoCbl, and about 15.5 kcal mol(-1) for the Ado(Im)Cbl, enough to completely explain the observed enzymatic catalysis. These results are discussed in the light of current knowledge about class I AdoCbl-dependent enzymes, in which the coenzyme is bound in its "base-off" conformation, with the lower axial ligand position occupied by the imidazole moiety of an active site histidine residue, and the class II enzymes, in which AdoCbl binds to the enzyme in its "base-on" conformation, and the pendent 5,6-dimethylbenzimidazole base remains coordinated to the metal during Co-C bond activation.     

O-O bond cleavage and alkane hydroxylation in methane monooxygenase

Several new aspects of the O-O bond cleavage and alkane hydroxylation mechanisms have been studied by hybrid density functional theory in this reinvestigation of methane monooxygenase. As concerning key intermediates in these reactions, a new important low-lying state is found, described either as Fe2(III,V) or as Fe2(III,IV)O. A fully optimized transition state for O-O bond cleavage has been determined. It is suggested that the large difference in optimal size (as determined in gas phase) of the complex, before and after the O-O bond cleavage, leads to an additional driving force for the reaction, not considered previously. The strain of the enzyme is estimated to lead to a driving force in the forward direction of about 5 kcal/mol, which could explain some of the pH dependence found in recent experiments. For the hydroxylation reaction, a clean hydrogen abstraction transition state leading to a substrate radical is again found, in contrast to interpretations of radical clock experiments. An explanation, based on new results, is suggested that could account for both the experimental and theoretical results.