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

Key enzymes catalyzing glycerol to 1,3-propanediol

Biodiesel can replace petroleum diesel as it is produced from animal fats and vegetable oils, and it produces about 10 % (w/w) glycerol, which is a promising new industrial microbial carbon, as a major by-product. One of the most potential applications of glycerol is its biotransformation to high value chemicals such as 1,3-propanediol (1,3-PD), dihydroxyacetone (DHA), succinic acid, etc., through microbial fermentation. Glycerol dehydratase, 1,3-propanediol dehydrogenase (1,3-propanediol-oxydoreductase), and glycerol dehydrogenase, which were encoded, respectively, by dhaB, dhaT, and dhaD and with DHA kinase are encompassed by the dha regulon, are the three key enzymes in glycerol bioconversion into 1,3-PD and DHA, and these are discussed in this review article. The summary of the main research direction of these three key enzyme and methods of glycerol bioconversion into 1,3-PD and DHA indicates their potential application in future enzymatic research and industrial production, especially in biodiesel industry.

Glutamate 338 is an electrostatic facilitator of C-Co bond breakage in a dynamic/electrostatic model of catalysis by ornithine aminomutase

How cobalamin-dependent enzymes promote C–Co homolysis to initiate radical catalysis has been debated extensively. For the pyridoxal 5′-phosphate and cobalamin-dependent enzymes lysine 5,6-aminomutase and ornithine 4,5-aminomutase (OAM), large-scale re-orientation of the cobalamin-binding domain linked to C–Co bond breakage has been proposed. In these models, substrate binding triggers dynamic sampling of the B₁₂-binding Rossmann domain to achieve a catalytically competent ‘closed’ conformational state. In ‘closed’ conformations of OAM, Glu338 is thought to facilitate C–Co bond breakage by close association with the cobalamin adenosyl group. We investigated this using stopped-flow continuous-wave photolysis, viscosity dependence kinetic measurements, and electron paramagnetic resonance spectroscopy of a series of Glu338 variants. We found that substrate-induced C–Co bond homolysis is compromised in Glu388 variant forms of OAM, although photolysis of the C–Co bond is not affected by the identity of residue 338. Electrostatic interactions of Glu338 with the 5′-deoxyadenosyl group of B₁₂ potentiate C–Co bond homolysis in ‘closed’ conformations only; these conformations are unlocked by substrate binding. Our studies extend earlier models that identified a requirement for large-scale motion of the cobalamin domain. Our findings indicate that large-scale motion is required to pre-organize the active site by enabling transient formation of ‘closed’ conformations of OAM. In ‘closed’ conformations, Glu338 interacts with the 5′-deoxyadenosyl group of cobalamin. This interaction is required to potentiate C–Co homolysis, and is a crucial component of the approximately 10¹² rate enhancement achieved by cobalamin-dependent enzymes for C–Co bond homolysis.

Role of active site residues in promoting cobalt-carbon bond homolysis in adenosylcobalamin-dependent mutases revealed through experiment and computation

Adenosylcobalamin (AdoCbl) serves as a source of reactive free radicals that are generated by homolytic scission of the coenzyme's cobalt-carbon bond. AdoCbl-dependent enzymes accelerate AdoCbl homolysis by ∼10(12)-fold, but the mechanism by which this is accomplished remains unclear. We have combined experimental and computational approaches to gain molecular-level insight into this process for glutamate mutase. Two residues, glutamate 330 and lysine 326, form hydrogen bonds with the adenosyl group of the coenzyme. A series of mutations that impair the enzyme's ability to catalyze coenzyme homolysis and tritium exchange with the substrate by 2-4 orders of magnitude were introduced at these positions. These mutations, together with the wild-type enzyme, were also characterized in silico by molecular dynamics simulations of the enzyme-AdoCbl-substrate complex with AdoCbl modeled in the associated (Co-C bond formed) or dissociated [adenosyl radical with cob(II)alamin] state. The simulations reveal that the number of hydrogen bonds between the adenosyl group and the protein side chains increases in the homolytically dissociated state, with respect to the associated state, for both the wild-type and mutant enzymes. The mutations also cause a progressive increase in the mean distance between the 5'-carbon of the adenosyl radical and the abstractable hydrogen of the substrate. Interestingly, the distance between the 5'-carbon and substrate hydrogen, determined computationally, was found to inversely correlate with the log k for tritium exchange (r = 0.93) determined experimentally. Taken together, these results point to a dual role for these residues: they both stabilize the homolytic state through electrostatic interactions between the protein and the dissociated coenzyme and correctly position the adenosyl radical to facilitate the abstraction of hydrogen from the substrate.

Phenylvinylcobalamin: an alkenylcobalamin featuring a ligand with a large trans influence

Cob(I)alamin reacts with phenylacetylene to produce two diastereomers in which the organic ligand is coordinated to the upper (β) and lower (α) face of the corrin ring, respectively. The isomers were separated chromatographically and characterised by ESI-MS and, in the case of the β isomer, by (1)H and (13)C NMR. Only the β isomer crystallised and its molecular structure, determined by X-ray diffraction, shows that the organic ligand coordinates Co(III) through the β carbon of the phenylvinyl ligand. The Co-C bond length is 2.004(8) Å while the Co-N bond length to the trans 5,6-dimethylbenzimidazole (dmbzm) base is 2.217(8) Å, one of the longest Co-Ndmbzm bond lengths known in an organocobalamin. Unlike benzylcobalamin (BzCbl), phenylvinylcobalamin (PhVnCbl) is stable towards homolysis. DFT calculations (BP86/TZVP) on model compounds of BzCbl and PhVnCbl show that the Co-C bond dissociation energy for homolysis to Co(II) and an organic radical in the former is 8 kcal mol(-1) lower than in the latter. An analysis of the electron density at the Co-C bond critical point using Bader's QTAIM approach shows that the Co-C bond in PhVnCbl is shorter, stronger and somewhat more covalent than that in BzCbl, and has some multiple bond character. Together with calculations that show that the benzyl radical is more stable than the phenylvinyl radical, this rationalises the stability of PhVnCbl compared to BzCbl. The phenylvinyl ligand has a large trans influence. The pKa for deprotonation of dmbzm and its coordination by the metal in β-PhVnCbl is 4.60 ± 0.01, one of the highest values reported to date in cobalamin chemistry. The displacement of dmbzm ligand by CN(-) in β-PhVnCbl occurs with log K = 0.7 ± 0.1; the trans influence order of C-donor ligands is therefore CN(-) < CCH < CHCH2 = PhVn < Me < Et.

Combined spectroscopic/computational studies of vitamin B12 precursors: geometric and electronic structures of cobinamides

Vitamin B(12) (cyanocobalamin) and its biologically active derivatives, methylcobalamin and adenosylcobalamin, are members of the family of corrinoids, which also includes cobinamides. As biological precursors to cobalamins, cobinamides possess the same structural core, consisting of a low-spin Co(3+) ion that is ligated equatorially by the four nitrogens of a highly substituted tetrapyrrole macrocycle (the corrin ring), but differ with respect to the lower axial ligation. Specifically, cobinamides possess a water molecule instead of the nucleotide loop that coordinates axially to Co(3+)cobalamins via its dimethylbenzimidazole (DMB) base. Compared to the cobalamin species, cobinamides have proven much more difficult to study experimentally, thus far eluding characterization by X-ray crystallography. In this study, we have utilized combined quantum mechanics/molecular mechanics (QM/MM) computations to generate complete structural models of a representative set of cobinamide species with varying upper axial ligands. To validate the use of this approach, analogous QM/MM geometry optimizations were carried out on entire models of the cobalamin counterparts for which high-resolution X-ray structural data are available. The accuracy of the cobinamide structures was assessed further by comparing electronic absorption spectra computed using time-dependent density functional theory to those obtained experimentally. Collectively, the results obtained in this study indicate that the DMB → H(2)O lower axial ligand switch primarily affects the energies of the Co 3d(z(2))-based molecular orbital (MO) and, to a lesser extent, the other Co 3d-based MOs as well as the corrin π-based highest energy MO. Thus, while the energy of the lowest-energy electronic transition of cobalamins changes considerably as a function of the upper axial ligand, it is nearly invariant for the cobinamides.

Axial bonding in alkylcobalamins: DFT analysis of the inverse versus normal trans influence

Density functional theory has been applied to study the origin of the inverse and normal trans influence in alkylcobalamins. In order to cover the X-ray structural data available for alkylcobalamins with a variety of axial substituents, geometries of 28 related corrin-containing models have been optimized and analyzed. The BP86/6-31G(d) level of theory was applied which showed good reliability in reproducing the axial bond lengths. Comparison of experimental and calculated data allowed to conclude that the inverse trans influence is not a general feature of cobalamins, as it appeared from the experimental data analysis alone. Inverse trans influence is observed for the series of R groups with increasing bulk and electron donating ability. For the series of R groups having similar medium bulk, but differing significantly in the electron donating ability, normal trans influence was found. Finally, it was determined, that the axial bond lengths correlate well but differently in the two series of R groups with the orbital energies of the six molecular orbitals essential in axial interligand bonding.

DFT and ONIOM(DFT:MM) studies on Co-C bond cleavage and hydrogen transfer in B12-dependent methylmalonyl-CoA mutase. Stepwise or concerted mechanism?

The considerable protein effect on the homolytic Co-C bond cleavage to form the 5'-deoxyadenosyl (Ado) radical and cob(II)alamin and the subsequent hydrogen transfer from the methylmalonyl-CoA substrate to the Ado radical in the methylmalonyl-CoA mutase (MMCM) have been extensively studied by DFT and ONIOM(DFT/MM) methods. Several quantum models have been used to systematically study the protein effect. The calculations have shown that the Co-C bond dissociation energy is very much reduced in the protein, compared to that in the gas phase. The large protein effect can be decomposed into the cage effect, the effect of coenzyme geometrical distortion, and the protein MM effect. The largest contributor is the MM effect, which mainly consists of the interaction of the QM part of the coenzyme with the MM part of the coenzyme and the surrounding residues. In particular, Glu370 plays an important role in the Co-C bond cleavage process. These effects tremendously enhance the stability of the Co-C bond cleavage state in the protein. The initial Co-C bond cleavage and the subsequent hydrogen transfer were found to occur in a stepwise manner in the protein, although the concerted pathway for the Co-C bond cleavage coupled with the hydrogen transfer is more favored in the gas phase. The assumed concerted transition state in the protein has more deformation of the coenzyme and the substrate and has less interaction with the protein than the stepwise route. Key factors and residues in promoting the enzymatic reaction rate have been discussed in detail.

X-ray structural characterization of imidazolylcobalamin and histidinylcobalamin: cobalamin models for aquacobalamin bound to the B12 transporter protein transcobalamin

The X-ray structures of imidazolylcobalamin (ImCbl) and histidinylcobalamin (HisCbl) are reported. These structures are of interest given that the recent structures of human and bovine transcobalamin prepared in their holo forms from aquacobalamin show a histidine residue of the metalloprotein bound at the beta-axial site of the cobalamin (Wuerges, J. et al. Proc. Natl. Acad. Sci. U.S.A. 2006, 103, 4386-4391). The beta-axial Co-N bond distances for ImCbl and HisCbl are 1.94(1) and 1.951(7) A, respectively. The alpha-axial Co-N bond distances to the 5,6-dimethylbenzimidazole are 2.01(1) and 1.979(8) A for ImCbl and HisCbl, respectively, and are typical for cobalamins with weak sigma-donor ligands at the beta-axial site. The corrin fold angles of 11.8(3) degrees (ImCbl) and 12.0(3) degrees (HisCbl) are smaller than those typically observed for cobalamins.

NMR observations of 13C-enriched coenzyme B12 bound to the ribonucleotide reductase from Lactobacillus leichmannii

The 13C NMR resonance and one-bond 1H-13C coupling constants of coenzyme B12 enriched in 13C in the cobalt-bound carbon have been observed in the complex of the coenzyme with the B12-dependent ribonucleotide reductase from Lactobacillus leichmannii. Neither the 13C NMR chemical shift nor the 1H-13C coupling constants are significantly altered by binding of the coenzyme to the enzyme. The results suggest that ground-state Co-C bond distortion is not utilized by this enzyme to activate coenzyme B12 for C-Co bond homolysis.

The enzymatic activation of coenzyme B12

The enzymatic "activation" of coenzyme B12 (5'-deoxyadenosylcobalamin, AdoCbl), in which homolysis of the carbon-cobalt bond of the coenzyme is catalyzed by some 10(9)- to 10(14)-fold, remains one of the outstanding problems in bioinorganic chemistry. Mechanisms which feature the enzymatic manipulation of the axial Co-N bond length have been investigated by theoretical and experimental methods. Classical mechanochemical triggering, in which steric compression of the long axial Co-N bond leads to increased upward folding of the corrin ring and stretching of the Co-C bond is found to be feasible by molecular modeling, but the strain induced in the Co-C bond seems to be too small to account for the observed catalytic power. The modeling study shows that the effect is a steric one which depends on the size of the axial nucleotide base, as substitution of imidazole (Im) for the normal 5,6-dimethylbenzimidazole (Bzm) axial base decreases the Co-C bond labilization considerably. An experimental test was thus devised using the coenzyme analog with Im in place of Bzm (Ado(Im)Cbl). Studies of the enzymatic activation of this analog by the B12-dependent ribonucleoside triphosphate reductase from Lactobacillus leichmannii coupled with studies of the non-enzymatic homolytic lability of the Co-C bond of Ado(Im)Cbl show that the enzyme is only slightly less efficient (3.8-fold, 0.8 kcal mol(-1)) at activating Ado(Im)Cbl than at activating AdoCbl itself. This suggests, in agreement with the modeling study, that mechanochemical triggering can make only a small contribution to the enzymatic activation of AdoCbl. Another possibility, electronic stabilization of the Co(II) homolysis product by compression of the axial Co-N bond, requires that enzymatic activation be sensitive to the basicity of the axial nucleotide. Preliminary studies of the enzymatic activation of a coenzyme analog with a 5-fluoroimidazole axial nucleotide suggest that the catalysis of Co-C bond homolysis may indeed be significantly slowed by the decrease in basicity.

Validation of semiempirical methods for modeling of corrinoid systems

Several semiempirical methods (MNDO-d, PM3tm, PM3-d, PM5, PM6, and AM1-d) have been tested against experimental data and density functional theory (DFT) results in search for the best methods that can be used for quantum-mechanical-molecular mechanics (QM/MM) modeling of corrinoid systems of vitamin B(12) co-factor. It has been found that the PM6 parametrization in its present form gives results closest to hybrid DFT calculations that are most widely used thus far. In comparison with pure DFT and experimental data the best agreement is obtained for PM3tm parametrization, while PM6 yields slightly worse results. AM1-d yields bad geometry of the corrin moiety. The worst performance was observed for MNDO-d, which has severe problem with position and orientation of the alpha-ligands.

Solution Structure and Thermolysis of Coβ-5‘-Deoxyadenosylimidazolylcobamide, a Coenzyme B12 Analogue with an Imidazole Axial Nucleoside

The solution structure of Cobeta-5'-deoxyadenosylimidazolylcobamide, Ado(Im)Cbl, the coenzyme B(12) analogue in which the axial 5,6-dimethylbenzimidazole (Bzm) ligand is replaced by imidazole, has been determined by NMR-restrained molecular modeling. A two-state model, in which a conformation with the adenosyl moiety over the southern quadrant of the corrin and a conformation with the adenosyl ligand over the eastern quadrant of the corrin are both populated at room temperature, was required by the nOe data. A rotation profile and molecular dynamics simulations suggest that the eastern conformation is the more stable, in contrast to AdoCbl itself in which the southern conformation is preferred. Consensus structures of the two conformers show that the axial Co-N bond is slightly shorter and the corrin ring is less folded in Ado(Im)Cbl than in AdoCbl. A study of the thermolysis of Ado(Im)Cbl in aqueous solution (50-125 degrees C) revealed competing homolytic and heterolytic pathways as for AdoCbl but with heterolysis being 9-fold faster and homolysis being 3-fold slower at 100 degrees C than for AdoCbl. Determination of the pK(a)'s for the Ado(Im)Cbl base-on/base-off reaction and for the detached imidazole ribonucleoside as a function of temperature permitted correction of the homolysis and heterolysis rate constants for the temperature-dependent presence of the base-off species of Ado(Im)Cbl. Activation analysis of the resulting rate constants for the base-on species show that the entropy of activation for Ado(Im)Cbl homolysis (13.7 +/- 0.9 cal mol(-1) K(-1)) is identical with that of AdoCbl (13.5 +/- 0.7 cal mol(-1) K(-1)) but that the enthalpy of activation (34.8 kcal mol(-1)) is 1.0 +/- 0.4 kcal mol(-1) larger. The opposite effect is seen for heterolysis, where the enthalpies of activation are identical but the entropy of activation is 5 +/- 1 cal mol(-1) K(-1) less negative for Ado(Im)Cbl. Extrapolation to 37 degrees C provides a rate constant for Ado(Im)Cbl homolysis of 2.1 x 10(-9) s(-1), 4.3-fold smaller than for AdoCbl. Combined with earlier results for the enzyme-induced homolysis of Ado(Im)Cbl by the ribonucleoside triphosphate reductase from Lactobacillus leichmannii, the catalytic efficiency of the enzyme for homolysis of Ado(Im)Cbl at 37 degrees C can be calculated to be 4.0 x 10(8), 3.8-fold, or 0.8 kcal mol(-1), smaller than for AdoCbl. Thus, the bulky Bzm ligand makes at best a <1 kcal mol(-1) contribution to the enzymatic activation of coenzyme B(12).

Photolysis and recombination of adenosylcobalamin bound to glutamate mutase

Ultrafast spectroscopic measurements are used to determine the kinetics of homolysis and recombination for adenosylcobalamin bound in the active site of glutamate mutase. These are the first such measurements on an adenosylcobalamin-dependent enzyme. A short-lived intermediate is formed prior to formation of the cob(II)alamin radical. This intermediate was not observed upon photolysis of adenosylcobalamin in free solution. The intrinsic rate constant for geminate recombination for adenosylcobalamin bound to glutamate mutase is 1.08 +/- 0.10 ns-1, only 16% smaller than the rate constant measured in free solution, 1.39 +/- 0.06 ns-1, suggesting the protein does not greatly perturb the stability of the cobalt-carbon bond upon binding the coenzyme.

Adenosylcobinamide Plus Exogenous, Sterically Hindered, Putative Axial Bases: A Reinvestigation into the Cause of Record Levels of Co−C Heterolysis

A reinvestigation of an earlier Ph.D. thesis (Sirovatka, J. M. Ph.D. Thesis, Colorado State University, Fort Collins, CO, 1999) is reported herein. That thesis examined the thermolysis reaction of AdoCbi(+)BF(4)(-) in ethylene glycol solution with exogenous bases, N-methylimidazole (N-Me-Im) and the sterically hindered 1,2-dimethylimidazole, (1,2-Me(2)-Im), 2-methylpyridine (2-Me-py), and 2,6-dimethylpyridine (2,6-Me(2)-py). In the present work, multiple purities of each base have been utilized as a check to see if impurities in the nitrogenous bases are causing the observed homolysis and heterolysis product distributions as others have implied (Trommel, J. S.; Warncke, K.; Marzilli, L. G. J. Am. Chem. Soc. 2001, 123, 3358). The "impurity hypothesis" is disproven by a series of results, including the following: N-Me-Im displays an invariant 52 +/- 10% heterolysis and the 1,2-Me(2)-Im system displays an invariant 83 +/- 7% heterolysis as a function of different base purification methods. Moreover, 2-Me-py and 2,6-Me(2)-py also display an invariant approximately 16 +/- 5% heterolysis as a function of different purification methods. What is responsible for the high levels of Co-C heterolysis in the AdoCbi(+) plus sterically bulky base thermolyses was uncovered via a revisitation of our four, earlier alternative hypotheses for the enhanced Co-C heterolysis (Sirovatka, J. M.; Finke, R. G. Inorg. Chem. 1999, 38, 1697). Our prior number one alternative hypothesis is shown to be correct: the added bases simply deprotonate the ethylene glycol solvent, forming ethylene glycolate anion and base-H(+)() as the key agents behind the previously ill-understood Co-C heterolyses. Also reported are Co(II)Cbi(+) titrations with five bases (1,2-Me(2)-Im, N-Me-Im, pyridine, 2-MePy, and 2,6-Me(2)-py). These experiments confirm Marzilli and co-workers' (op. cit.) results by showing that sterically hindered bases do not bind to Co(II)Cbi(+); therefore, Co(II)Cbi(+) EPR literature showing binding of bulky pyridines is erroneous as is the previously reported binding of bulky pyridine bases to Co(II)Cbi(+) by UV-vis spectroscopy (Sirovatka, J. Ph.D. Thesis, op. cit.). Also reported is our current best synthesis and purification of AdoCbi(+)BF(4)(-), work that builds off our 1987 synthesis of AdoCbi(+)BF(4)(-) (Hay, B. P.; Finke, R. G. J. Am. Chem. Soc. 1987, 109, 8012). Finally, the multiple, compounding errors which have caused problems in this project are listed, notably errors in the protein X-ray crystallography literature, the EXAFS literature, the Co(II)Cbi(+) plus bulky-bases EPR literature, the misleading B(12)-model literature, the erroneous experimental work (Sirovatka, op. cit.) and thus incorrect conclusions in one of our prior papers, as well as the erroneous implications in parts of the Marzilli and co-workers paper (op. cit.). It is hoped that a forthright reporting of these errors will help others avoid similar mistakes in the future when studying complex, bioinorganic systems.

A Combined Density Functional Theory and Molecular Mechanics Study of the Relationship between the Structure of Coenzyme B12 and Its Binding to Methylmalonyl-CoA Mutase

A combined density functional theory (DFT) and molecular mechanics (MM) approach was applied to investigate the relationship between the structure of a free coenzyme B12, and bound to methylmalonyl-CoA mutase. It was found that, upon coenzyme binding to apoenzyme, the Co-C bond remains intact, while the C-Naxial bond becomes slightly elongated and labilized. The labilization of the Co-Naxial bond that takes place in coenzyme B12-dependent enzymes is most likely necessary for fine-tuning of the cobalt-nitrogen (axial base) distance. The controlling of this distance is important to inhibit abiological site reaction involving heterolysis of the Co-C bond but is not important for biologically relevant Co-C bond homolysis.

Spectroscopic and computational studies of Co3+-corrinoids: spectral and electronic properties of the B12 cofactors and biologically relevant precursors

The B(12) cofactors methylcobalamin (MeCbl) and 5'-deoxyadenosylcobalamin (AdoCbl) have long fascinated chemists because of their complex structures and unusual reactivities in biological systems; however, their electronic absorption (Abs) spectra have remained largely unassigned. In this study, we have used Abs, circular dichroism (CD), magnetic CD (MCD), and resonance Raman spectroscopic techniques to probe the electronic excited states of Co(3+)Cbl species that differ with respect to their upper axial ligand, including MeCbl, AdoCbl, aquacobalamin (H(2)OCbl(+)), and vitamin B(12) (cyanocobalamin, CNCbl). Also included to probe the effect of the lower axial ligand on the electronic properties of Cbls is Ado-cobinamide (AdoCbi(+)), an AdoCbl derivative that lacks the tethered base 5,6-dimethylbenzimidazole (DMB) and instead binds a water molecule in the lower axial position. Spectroscopic data for each species are analyzed within the framework of time-dependent density functional theory (TD-DFT) to assign the major spectral features (the so-called alpha/beta, D/E, and gamma bands) and to generate experimentally validated electronic-structure descriptions. These studies reveal that the "unique" Abs spectra of MeCbl and AdoCbl, which differ considerably from the "typical" Abs spectra of H(2)OCbl(+) and CNCbl, reflect the high degree of sigma-donation from the alkyl ligand to the Co center and the consequent destabilization of all Co 3d orbitals. They reveal further that with increasing sigma-donor strength of the upper axial ligand, the contribution from the formally unoccupied Co 3d(z(2)) orbital to the HOMO increases, which induces a strong Co[bond]N(DMB) sigma-antibonding interaction, consistent with the experimentally observed lengthening of this bond from H(2)OCbl(+) to CNCbl and MeCbl. Alternatively, our spectroscopic and computational data for MeCbl and MeCbi(+) reveal that substitution of the DMB by a water molecule in the lower axial position has negligible effects on the Co[bond]C. A simple model is presented that explains why the identity of the upper axial ligand has a major effect on the Co[bond]N(ax) strength, whereas the lower axial ligand does not appreciably modulate the nature of the Co[bond]C. Implications of these results with respect to enzymatic Co[bond]C activation are discussed.

Quantum chemical modeling of Co--C bond activation in B(12)-dependent enzymes

Recent progress in computational modeling of the catalytic activation of cobalt-carbon bond cleavage shows that quantum chemical calculations could be an important part of coenzyme B(12) research. Particular emphasis has been placed on density functional theory, which is now emerging as a powerful tool to elucidate the electronic structure and spectroscopic properties of the active sites of metalloenzymes.

Enzymatic activity of coenzyme B(12) derivatives with altered axial nucleotides: probing the mechanochemical triggering hypothesis in ribonucleotide reductase

Theoretical studies (J. Inorg. Biochem. 2001, 83, 121) of the involvement of the bulky 5,6-dimethylbenzimidazole (Dmbz) ligand of coenzyme B(12) (5'-deoxyadenosylcobalamin, AdoCbl) in the mechanism of activation of the carbon-cobalt bond of the coenzyme for homolytic cleavage by AdoCbl-dependent enzymes (the "mechanochemical triggering" mechanisms) have shown that a purely steric, ground-state mechanism can supply only a few kilocalories per mole (of the observed 13-16 kcal mol(-1)) of activation, but that an electronic mechanism, operating to stabilize the transition state, can explain all of the observed catalytic effect. To address these mechanisms experimentally, analogues of AdoCbl in which the Dmbz ligand is replaced by benzimidazole (Ado(Bzim)Cbl) or by imidazole (Ado(Im)Cbl) have been prepared and characterized. Both of these analogues support turnover in the AdoCbl-dependent ribonucleoside triphosphate reductase (RTPR) from Lactobacillus leichmannii at 100% of the activity of AdoCbl itself, but the Ado(Im)Cbl analogue has a significantly higher K(m). 5'-Deoxyadenosylcobinamide, the analogue in which the axial nucleotide has been chemically removed, in contrast, is inactive in the spectrophotometric assay, which indicates that it has at most 1% of the activity of AdoCbl. Stopped-flow spectrophotometric measurements of the formation of cob(II)alamin at the enzyme active site show that RTPR binds Ado(Bzim)Cbl slightly more weakly than it does AdoCbl, but binds Ado(Im)Cbl 8-fold more weakly. While the equilibrium constant for cob(II)alamin formation is nearly the same for Ado(Bzim)Cbl and AdoCbl, it is 5-fold smaller for Ado(Im)Cbl. Finally, the forward rate constant for enzyme-induced Co-C bond homolysis was about the same for Ado(Bzim)Cbl and for AdoCbl but was 17-fold smaller for Ado(Im)Cbl. These results are consistent with a small contribution from ground-state mechanochemical triggering, but they do not in themselves rule out transition-state mechanical triggering.

Theoretical evaluation of the hydrogen kinetic isotope effect on the first step of the methylmalonyl-CoA mutase reaction

We have calculated hydrogen kinetic isotope effects (KIEs) for the first step of the methylmalonyl-CoA mutase reaction, including multidimensional tunneling correction at the zero curvature (ZCT) level, and compared them with the experimental values. Both alternative mechanisms of this step, concerted and stepwise, can be accommodated. It turned out to be essential to include Arg207 hydrogen-bonded to the reactant in the mechanism predicting simultaneous breaking of the Co-C bond of AdoCbl and hydrogen atom transfer. The consequence of the stepwise mechanism is a much larger facilitation of the homolytic dissociation of the carbon-cobalt bond by the enzyme than currently appreciated; our results suggest lowering of the activation energy by about 23 kcal mol(-1). We have also shown that large hydrogen KIEs of tunneling origin do not necessarily break the Swain-Schaad equation. Furthermore, when this equation does not hold, the exponent may be smaller in the presence of tunneling than it is at the semi-classical limit, indicating that nonclassical behavior may be a more common phenomenon than expected.