Theoretical Investigation of Steric and Electronic Effects in Coenzyme B12 Models

The inorganic chemistry of the cobalt corrinoids - an update

The inorganic chemistry of the cobalt corrinoids, derivatives of vitamin B₁₂, is reviewed, with particular emphasis on equilibrium constants for, and kinetics of, their axial ligand substitution reactions. The role the corrin ligand plays in controlling and modifying the properties of the metal ion is emphasised. Other aspects of the chemistry of these compounds, including their structure, corrinoid complexes with metals other than cobalt, the redox chemistry of the cobalt corrinoids and their chemical redox reactions, and their photochemistry are discussed. Their role as catalysts in non-biological reactions and aspects of their organometallic chemistry are briefly mentioned. Particular mention is made of the role that computational methods - and especially DFT calculations - have played in developing our understanding of the inorganic chemistry of these compounds. A brief overview of the biological chemistry of the B₁₂-dependent enzymes is also given for the reader's convenience.

Resonance Raman Optical Activity Shows Unusual Structural Sensitivity for Systems in Resonance with Multiple Excited States: Vitamin B12 Case

In this work, cobalamins with different upper axial substituents and a cobalamin derivative with a ring modification were studied using chiroptical spectroscopies, in particular resonance Raman optical activity (RROA), to shed light on the influence of structural modifications on RROA spectra in these strongly chiral systems in resonance with multiple excited states at 532 nm excitation. We have demonstrated that for these unique systems RROA possesses augmented structural specificity, surpassing resonance Raman spectroscopy and enabling at the same time measurement of cobalamins at fairy low concentrations of ∼10^-5 mol dm^-3. The enhanced structural specificity of RROA is a result of bisignate spectra due to resonance via more than one electronic state. The observation of increased structural capability of RROA for cobalamins opens a new perspective for studying chiral properties of other biological systems incorporating d-metal ions.

Theoretical analysis of C-F bond cleavage mediated by cob[I]alamin-based structures

In the present work, C-F bond cleavage mediated by the super-reduced form of cobalamin (i.e., Co^ICbl) was theoretically studied at the ONIOM(BP86/6-311++G(d,p):PM6) + SMD level of theory. Dispersion effects were introduced by employing Grimme's empirical dispersion at the ONIOM(BP86-D/6-311++G(d,p):PM6) + SMD level. In the first stage of the study, cobalamin was characterized in terms of the coordination number of the central cobalt atom. The ONIOM(BP86/6-311++G(d,p):PM6) results showed that the base-off form of the system is slightly more stable than its base-on counterpart (ΔE = E _base-off - E _base-on ~ -2 kcal/mol). The inclusion of dispersive forces in the description of the system stabilizes the base-on form, which becomes as stable as its base-off counterpart. Moreover, in the latter case, the energy barrier separating both structures was found to be negligible, with a computed value of 1.02 kcal/mol. In the second stage of the work, the reaction Co^ICbl + CH₃F → MeCbl + F^- was studied considering the base-off and the base-on forms of Co^ICbl. The reaction that occurs in the presence of the base-on form of Co^ICbl was found to be kinetically more favorable (ΔE ^≠ = 13.7 kcal/mol) than that occurring in the presence of the base-off form (ΔE ^≠ = 41.2 kcal/mol). Further reaction-force analyses of the processes showed that the energy barrier to C-F bond cleavage arises largely due to structural rearrangements when the reaction occurs on the base-on form of the Co^ICbl complex, but is mainly due to electronic rearrangements when the reaction takes place on the base-off form of the complex. The latter behavior emerges from differences in the synchronicity of the bond strengthening/weakening processes along the reaction path; the base-on mode of Co^ICbl is able to decrease the synchronicity of the chemical events. This work gives new molecular-level insights into the role of Cbl-based systems in the cleavage of C-F bonds. These insights have potential implications for research into processes for degrading fluorine-containing pollutants.

Co-C dissociation of adenosylcobalamin (coenzyme B12): role of dispersion, induction effects, solvent polarity, and relativistic and thermal corrections

Quantum-chemical cluster modeling is challenged in the limit of large, soft systems by the effects of dispersion and solvent, and well as other physical interactions. Adenosylcobalamin (AdoCbl, coenzyme B12), as one of the most complex cofactors in life, constitutes such a challenge. The cleavage of its unique organometallic Co-C bond has inspired multiple studies of this cofactor. This paper reports the fully relaxed potential energy surface of Co-C cleavage of AdoCbl, including for the first time all side-chain interactions with the dissociating Ado group. Various methods and corrections for dispersion, relativistic effects, solvent polarity, basis set superposition error, and thermal and vibrational effects were investigated, totaling more than 550 single-point energies for the large model. The results show immense variability depending on method, including solvation, functional type, and dispersion, challenging the conceived accuracy of methods used for such systems. In particular, B3LYP-D3 seems to severely underestimate the Co-C bond strength, consistent with previous results, and BP86 remains accurate for cobalamins when dispersion interactions are accounted for.

Enhanced ‘in situ’ catalysis via microwave selective heating: catalytic chain transfer polymerisation

The formation of a CCTP active catalyst, which is promoted by adopting microwave heating and an “in situ” synthesis methodology.

Mercury methylation by HgcA: theory supports carbanion transfer to Hg(II)

Many proteins use corrinoid cofactors to facilitate methyl transfer reactions. Recently, a corrinoid protein, HgcA, has been shown to be required for the production of the neurotoxin methylmercury by anaerobic bacteria. A strictly conserved Cys residue in HgcA was predicted to be a lower-axial ligand to Co(III), which has never been observed in a corrinoid protein. Here, we use density functional theory to study homolytic and heterolytic Co-C bond dissociation and methyl transfer to Hg(II) substrates with model methylcobalamin complexes containing a lower-axial Cys or His ligand to cobalt, the latter of which is commonly found in other corrinoid proteins. We find that Cys thiolate coordination to Co facilitates both methyl radical and methyl carbanion transfer to Hg(II) substrates, but carbanion transfer is more favorable overall in the condensed phase. Thus, our findings are consistent with HgcA representing a new class of corrinoid protein capable of transferring methyl groups to electrophilic substrates.

Spectroscopic characterization of active-site variants of the PduO-type ATP:corrinoid adenosyltransferase from Lactobacillus reuteri: insights into the mechanism of four-coordinate Co(II)corrinoid formation

The PduO-type adenosine 5'-triphosphate (ATP):corrinoid adenosyltransferase from Lactobacillus reuteri (LrPduO) catalyzes the transfer of the adenosyl-group of ATP to Co(1+)cobalamin (Cbl) and Co(1+)cobinamide (Cbi) substrates to synthesize adenosylcobalamin (AdoCbl) and adenosylcobinamide (AdoCbi(+)), respectively. Previous studies revealed that to overcome the thermodynamically challenging Co(2+) → Co(1+) reduction, the enzyme drastically weakens the axial ligand-Co(2+) bond so as to generate effectively four-coordinate (4c) Co(2+)corrinoid species. To explore how LrPduO generates these unusual 4c species, we have used magnetic circular dichroism (MCD) and electron paramagnetic resonance (EPR) spectroscopic techniques. The effects of active-site amino acid substitutions on the relative yield of formation of 4c Co(2+)corrinoid species were examined by performing eight single-amino acid substitutions at seven residues that are involved in ATP-binding, an intersubunit salt bridge, and the hydrophobic region surrounding the bound corrin ring. A quantitative analysis of our MCD and EPR spectra indicates that the entire hydrophobic pocket below the corrin ring, and not just residue F112, is critical for the removal of the axial ligand from the cobalt center of the Co(2+)corrinoids. Our data also show that a higher level of coordination among several LrPduO amino acid residues is required to exclude the dimethylbenzimidazole moiety of Co(II)Cbl from the active site than to remove the water molecule from Co(II)Cbi(+). Thus, the hydrophilic interactions around and above the corrin ring are more critical to form 4c Co(II)Cbl than 4c Co(II)Cbi(+). Finally, when ATP analogues were used as cosubstrate, only "unactivated" five-coordinate (5c) Co(II)Cbl was observed, disclosing an unexpectedly large role of the ATP-induced active-site conformational changes with respect to the formation of 4c Co(II)Cbl. Collectively, our results indicate that the level of control exerted by LrPduO over the timing for the formation of the 4c Co(2+)corrinoid intermediates is even more exquisite than previously anticipated.

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.

Comparison of chemical properties of iron, cobalt, and nickel porphyrins, corrins, and hydrocorphins

Density functional calculations have been used to compare the geometric, electronic, and functional properties of the three important tetrapyrrole systems in biology, heme, coenzyme B 12, and coenzyme F430, formed from iron porphyrin ( Por ), cobalt corrin ( Cor ), and nickel hydrocorphin ( Hcor ). The results show that the flexibility of the ring systems follows the trend Hcor > Cor > Por and that the size of the central cavity follows the trend Cor < Por < Hcor . Therefore, low-spin Co I, Co II, and Co III fit well into the Cor ring, whereas Por seems to be more ideal for the higher spin states of iron, and the cavity in Hcor is tailored for the larger Ni ion, especially in the high-spin Ni II state. This is confirmed by the thermodynamic stabilities of the various combinations of metals and ring systems. Reduction potentials indicate that the +I and +III states are less stable for Ni than for the other metal ions. Moreover, Ni – C bonds are appreciably less stable than Co - C bonds. However, it is still possible that a Ni – CH 3 bond is formed in F 430 by a heterolytic methyl transfer reaction, provided that the donor is appropriate, e.g. if coenzyme M is protonated. This can be facilitated by the adjacent SO 3− group in this coenzyme and by the axial glutamine ligand, which stabilizes the Ni III state. Our results also show that a Ni III– CH 3 complex is readily hydrolysed to form a methane molecule and that the Ni III hydrolysis product can oxidize coenzyme B and M to a heterodisulphide in the reaction mechanism of methyl coenzyme M reductase.

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.

Computational Studies of Bioorganometallic Enzymes and Cofactors

Because of their complex geometric and electronic structures, the active sites and cofactors of bioorganometallic enzymes, which are characterized by their metal–carbon bonds, pose a major challenge for computational chemists. However, recent progress in computer technology and theoretical chemistry, along with insights gained from mechanistic, spectroscopic, and X-ray crystallographic studies, have established an excellent foundation for the successful completion of computational studies aimed at elucidating the electronic structures and catalytic cycles of these species. This chapter briefly reviews the most popular computational approaches employed in theoretical studies of bioorganometallic species and summarizes important information obtained from computational studies of (i) the enzymatic formation and cleavage of the Co–C bond of coenzyme B12; (ii) the catalytic cycle of methyl-coenzyme M reductase and its nickel-containing cofactor F430; (iii) the polynuclear active-site clusters of the bifunctional enzyme carbon monoxide dehydrogenase/acetyl-coenzyme A synthase; and (iv) the magnetic properties of the active-site cluster of Fe-only hydrogenases.

Electronic structure of Cob(I)alamin: the story of an unusual nucleophile

The electronic structure and the ligand-field spectrum of cobalt(I) corrin is reported using complete active space multiconfigurational perturbation theory (CASPT2) to address some inconsistencies and the nature of the cobalt(I) "supernucleophile", cob(I)alamin. An assignment of six of the seven intense lines in the experimental spectrum is obtained at a root-mean-square accuracy of 0.14 eV and largest error of 0.21 eV. Agreement is significantly better for CASPT2 than density functional theory (DFT), but DFT does surprisingly well. The correlated wave function implies that the ground state of Co(I) corrin is severely multiconfigurational, with only 67% of the d(8) reference configuration and prominent contributions of 20% from open-shell metal-to-ligand charge-transfer configurations. The ground state exhibits a fascinating degree of covalency between cobalt and the nitrogen orbitals, described by the bonding and antibonding orbital pair of a cobalt d-orbital and a delta-orbital linearly combined from nitrogen orbitals. Thus, the standard description of the d(8) supernucleophile is not completely valid. From a biological perspective, the mixing in of Co(II) configurations in cob(I)alamin may be an important reason for the redox accessibility of the formal Co(I) state of the cofactor, which again provides the catalytic power for one half-reaction of enzymes such as cobalamin-dependent methionine synthase.

Computational insights into the mechanism of radical generation in B12-dependent methylmalonyl-CoA mutase

ONIOM calculations have provided novel insights into the mechanism of homolytic Co-C5' bond cleavage in the 5'-deoxyadenosylcobalamin cofactor catalyzed by methylmalonyl-CoA mutase. We have shown that it is a stepwise process in which conformational changes in the 5'-deoxyadenosine moiety precede the actual homolysis step. In the transition state structure for homolysis, the Co-C5' bond elongates by approximately 0.5 Angstroms from the value found in the substrate-bound reactant complex. The overall barrier to homolysis is approximately 10 kcal/mol, and the radical products are approximately 2.5 kcal/mol less stable than the initial ternary complex of enzyme, substrate, and cofactor. The movement of the deoxyadenosine moiety during the homolysis step positions the resulting 5'-deoxyadenosyl radical for the subsequent hydrogen atom transfer from the substrate, methylmalonyl-CoA.

Computational study on the difference between the Co–C bond dissociation energy in methylcobalamin and adenosylcobalamin

The bond dissociation energies of the Co-C bonds in the cobalamin cofactors methylcobalamin and adenosylcobalamin were calculated using the hybrid quantum mechanics/molecular mechanics method IMOMM (integrated molecular orbital and molecular mechanics). Calculations were performed on models of differing complexities as well as on the full systems. We investigated the origin of the different experimental values for the Co-C bond dissociation energies in methylcobalamin and adenosylcobalamin, and have provided an explanation for the difficulties encountered when we attempt to reproduce this difference in quantum chemistry. Additional calculations have been performed using the Miertus-Scrocco-Tomasi method in order to estimate the influence of solvent effects on the homolytic Co-C bond cleavage. Introduction of these solvation effects is shown to be necessary for the correct reproduction of experimental trends in bond dissociation energies in solution, which consequently have no direct correlation with dissociation processes in the enzyme.

Electronic structure and bonding in hydroxocobalamin

The electronic structure of hydroxocobalamin (OHCbl) has been calculated by a density functional method, using the orthogonalized linear combination of the atomic orbitals method (OLCAO). The X-ray crystal structure has been determined from synchrotron X-ray diffraction data and the geometry determined was used in the calculations. Comparison with the recently reported electronic structures of cyanocobalamin (CNCbl), methylcobalamin (MeCbl) and adenosylcobalamin (AdoCbl) shows that Mulliken charges (Q*) and bond orders (BO) vary only on the axial fragment.

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.

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.

Conversion of Homocysteine to Methionine by Methionine Synthase: A Density Functional Study

This communication reports a theoretical study of the conversion of homocysteine to methionine by methionine synthase. The reaction pathway is based on density functional calculations with large basis sets, including thermodynamic, relativistic, and solvent effects. We find that the suggested SN2 mechanism explains well the experimentally observed reaction rate. The results show that the reaction is highly polar, as reflected in the change of charge density along the reaction coordinate. It is enhanced in the protein by two effects: deprotonation of the bound substrate and desolvation of substrate and cofactor in the rate-determining step.

Comparison of the chemical properties of iron and cobalt porphyrins and corrins

Density functional calculations have been used to compare various geometric, electronic and functional properties of iron and cobalt porphyrin (Por) and corrin (Cor) species. The investigation is focussed on octahedral M(II/III) complexes (where M is the metal) with two axial imidazole ligands (as a model of b and c type cytochromes) or with one imidazole and one methyl ligand (as a model of methylcobalamin). However, we have also studied some five-coordinate M(II) complexes with an imidazole ligand and four-coordinate M(I/II) complexes without any axial ligands as models of other intermediates in the reaction cycle of coenzyme B12. The central cavity of the corrin ring is smaller than that of porphine. We show that the cavity of corrin is close to ideal for low-spin Co(III), Co(II), and Co(I) with the axial ligands encountered in biology, whereas the cavity in porphine is better suited for intermediate-spin states. Therefore, the low-spin state of Co is strongly favoured in complexes with corrins, whereas there is a small energy difference between the various spin states in iron porphyrin species. There are no clear differences for the reduction potentials of the octahedral complexes, but [Co(I)Cor] is more easily formed (by at least 40 kJ mole(-1)) than [Fe(I)Por]. Cobalt and corrin form a strong Cobond;C bond that is more stable against hydrolysis than iron and porphine. Finally, Fe(II/III) gives a much lower reorganisation energy than Co(II/III); this is owing to the occupied d(z2) orbital in Co(II). Altogether, these results give some clues about how nature has chosen the tetrapyrrole rings and their central metal ion.

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.

Equilibrium and kinetic studies on the reactions of alkylcobalamins with cyanide

Ligand substitution equilibria of different alkylcobalamins (RCbl, R = Me, CH(2)Br, CH(2)CF(3), CHF(2), CF(3)) with cyanide have been studied. It was found that CN(-) first substitutes the 5,6-dimethylbenzimidazole (Bzm) moiety in the alpha-position, followed by substitution of the alkyl group in the beta-position trans to Bzm. The formation constants K(CN) for the 1:1 cyanide adducts (R(CN)Cbl) were found to be 0.38 +/- 0.03, 0.43 +/- 0.03, and 123 +/- 9 M(-1) for R = Me, CH(2)Br, and CF(3), respectively. In the case of R = CH(2)CF(3), the 1:1 adduct decomposes in the dark with CN(-) to give (CN)(2)Cbl. The unfavorable formation constants for R = Me and CH(2)Br indicate the requirement of very high cyanide concentrations to produce the 1:1 complex, which cause the kinetics of the displacement of Bzm to be too fast to follow kinetically. The kinetics of the displacement of Bzm by CN(-) could be followed for R = CH(2)CF(3) and CF(3) to form CF(3)CH(2)(CN)Cbl and CF(3)(CN)Cbl, respectively, in the rate-determining step. Both reactions show saturation kinetics at high cyanide concentration, and the limiting rate constants are characterized by the activation parameters: R = CH(2)CF(3), DeltaH = 71 +/- 1 kJ mol(-1), DeltaS = -25 +/- 4 J K(-1) mol(-1), and DeltaV = +8.9 +/- 1.0 cm(3) mol(-1); R = CF(3), DeltaH = 77 +/- 3 kJ mol(-1), DeltaS = +44 +/- 11 J K(-1) mol(-1), and DeltaV = +14.8 +/- 0.8 cm(3) mol(-1), respectively. These parameters are interpreted in terms of an I(d) and D mechanism for R = CH(2)CF(3) and CF(3), respectively. The results of the study enable the formulation of a general mechanism that can account for the substitution behavior of all investigated alkylcobalamins including coenzyme B(12).

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.