A Density Functional Study on the Effect of the Trans Axial Ligand of Cobalamin on the Homolytic Cleavage of the Co−C Bond

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.

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.

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.

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.

Characterization of protein contributions to cobalt-carbon bond cleavage catalysis in adenosylcobalamin-dependent ethanolamine ammonia-lyase by using photolysis in the ternary complex

Protein contributions to the substrate-triggered cleavage of the cobalt-carbon (Co-C) bond and formation of the cob(II)alamin-5'-deoxyadenosyl radical pair in the adenosylcobalamin (AdoCbl)-dependent ethanolamine ammonia-lyase (EAL) from Salmonella typhimurium have been studied by using pulsed-laser photolysis of AdoCbl in the EAL-AdoCbl-substrate ternary complex, and time-resolved probing of the photoproduct dynamics by using ultraviolet-visible absorption spectroscopy on the 10(-7)-10(-1) s time scale. Experiments were performed in a fluid dimethylsulfoxide/water cryosolvent system at 240 K, under conditions of kinetic competence for thermal cleavage of the Co-C bond in the ternary complex. The static ultraviolet-visible absorption spectra of holo-EAL and ternary complex are comparable, indicating that the binding of substrate does not labilize the cofactor cobalt-carbon (Co-C) bond by significantly distorting the equilibrium AdoCbl structure. Photolysis of AdoCbl in EAL at 240 K leads to cob(II)alamin-5'-deoxyadenosyl radical pair quantum yields of <0.01 at 10(-6) s in both holo-EAL and ternary complex. Three photoproduct states are populated following a saturating laser pulse, and labeled, P(f), P(s), and P(c). The relative amplitudes and first-order recombination rate constants of P(f) (0.4-0.6; 40-50 s(-1)), P(s) (0.3-0.4; 4 s(-1)), and P(c) (0.1-0.2; 0) are comparable in holo-EAL and in the ternary complex. Time-resolved, full-spectrum electron paramagnetic resonance (EPR) spectroscopy shows that visible irradiation alters neither the kinetics of thermal cob(II)alamin-substrate radical pair formation, nor the equilibrium between ternary complex and cob(II)alamin-substrate radical pair, at 246 K. The results indicate that substrate binding to holo-EAL does not "switch" the protein to a new structural state, which promptly stabilizes the cob(II)alamin-5'-deoxyadenosyl radical pair photoproduct, either through an increased barrier to recombination, a decreased barrier to further radical pair separation, or lowering of the radical pair state free energy, or a combination of these effects. Therefore, we conclude that such a change in protein structure, which is independent of changes in the AdoCbl structure, and specifically the Co-C bond length, is not a basis of Co-C bond cleavage catalysis. The results suggest that, following the substrate trigger, the protein interacts with the cofactor to contiguously guide the cleavage of the Co-C bond, at every step along the cleavage coordinate, starting from the equilibrium configuration of the ternary complex. The cleavage is thus represented by a diagonal trajectory across a free energy surface, that is defined by chemical (Co-C separation) and protein configuration coordinates.

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 Mutation Analysis of Hydrogen Abstraction and Radical Rearrangement Steps in the Catalysis of Coenzyme B12-Dependent Diol Dehydratase

A mutation analysis of the catalytic functions of active-site residues of coenzyme B(12)-dependent diol dehydratase in the conversion of 1,2-propanediol to 1,1-propanediol has been carried out by using QM/MM computations. Mutants His143Ala, Glu170Gln, Glu170Ala, and Glu170Ala/Glu221Ala were considered to estimate the impact of the mutations of His143 and Glu170. In the His143Ala mutant the activation energy for OH migration increased to 16.4 from 11.5 kcal mol(-1) in the wild-type enzyme. The highest activation energy, 19.6 kcal mol(-1), was measured for hydrogen back-abstraction in this reaction. The transition state for OH migration is not sufficiently stabilized by the hydrogen-bonding interaction formed between the spectator OH group and Gln170 in the Glu170Gln mutant, which demonstrates that a strong proton acceptor is required to promote OH migration. In the Glu170Ala mutant, a new strong hydrogen bond is formed between the spectator OH group and Glu221. A computed activation energy of 13.6 kcal mol(-1) for OH migration in the Glu170Ala mutant is only 2.1 kcal mol(-1) higher than the corresponding barrier in the wild-type enzyme. Despite the low activation barrier, the Glu170Ala mutant is inactive because the subsequent hydrogen back-abstraction is energetically demanding in this mutant. OH migration is not feasible in the Glu170Ala/Glu221Ala mutant because the activation barrier for OH migration is greatly increased by the loss of COO(-) groups near the spectator OH group. This result indicates that the effect of partial deprotonation of the spectator OH group is the most important factor in reducing the activation barrier for OH migration in the conversion of 1,2-propanediol to 1,1-propanediol catalyzed by diol dehydratase.

Time-resolved measurements of the photolysis and recombination of adenosylcobalamin bound to glutamate mutase

Femtosecond to nanosecond transient absorption spectroscopy is used to investigate the photolysis of 5'-deoxyadenosylcobalamin (coenzyme B12, AdoCbl) bound to glutamate mutase. The photochemistry of AdoCbl is found to be inherently dependent upon the environment of the cofactor. Excitation of AdoCbl bound to glutamate mutase results in formation of a metal-to-ligand charge transfer intermediate state which decays to form cob(II)alamin with a time constant of 105 ps. This observation is in contrast to earlier measurements in water where the photohomolysis proceeds through an intermediate state in which the axial dimethylbenzimidazole ligand appears to have dissociated, and measurements in ethylene glycol where prompt bond homolysis is observed (Yoder, L. M.; Cole, A. G.; Walker, L. A., II; Sension, R. J. J. Phys. Chem. B 2001, 105, 12180-12188). The quantum yield for formation of stable radical pairs in the enzyme is found to be phi = 0.05 +/- 0.03, and the resulting intrinsic rate constants for geminate recombination and "cage escape" are 1.0 +/- 0.1 and 0.05 +/- 0.03 ns(-1), respectively. The rate constant for geminate recombination is 30% less than that observed for AdoCbl in water or ethylene glycol. This reduction is insufficient to account for the 10(12)-fold increase in the homolysis rate observed when substrate is bound to the protein. Finally, the protein provides a cage to prevent diffusive loss of the adenosyl radical; however, the ultimate yield for long-lived radicals is determined by the evolution from a singlet to a triplet radical pair as proposed for AdoCbl in ethylene glycol.

Performance of DFT in modeling electronic and structural properties of cobalamins

Computational modeling of the enzymatic activity of B12-dependent enzymes requires a detailed understanding of the factors that influence the strength of the Co--C bond and the limits associated with a particular level of theory. To address this issue, a systematic analysis of the electronic and structural properties of coenzyme B12 models has been performed to establish the performance of three different functionals including B3LYP, BP86, and revPBE. In particular the cobalt-carbon bond dissociation energies, axial bond lengths, and selected stretching frequencies have been analyzed in detail. Current analysis shows that widely used B3LYP functional significantly underestimates the strength of the Co--C bond while the nonhybrid BP86 functional produces very consistent results in comparison to experimental data. To explain such different performance of these functionals molecular orbital analysis associated with axial bonds has been performed to show differences in axial bonding provided by hybrid and nonhybrid functionals.

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.

Influence of Environment on the Electronic Structure of Cob(III)alamins: Time-Resolved Absorption Studies of the S1 State Spectrum and Dynamics

Transient absorption spectroscopy has been used to elucidate the nature of the S1 intermediate state populated following excitation of cob(III)alamin (Cbl(III)) compounds. This state is sensitive both to axial ligation and to solvent polarity. The excited-state lifetime as a function of temperature and solvent environment is used to separate the dynamic and electrostatic influence of the solvent. Two distinct types of excited states are identified, both assigned to pi3d configurations. The spectra of both types of excited states are characterized by a red absorption band (ca. 600 nm) assigned to Co 3d --> 3d or Co 3d --> corrin pi* transitions and by visible absorption bands similar to the corrin pi-->pi* transitions observed for ground state Cbl(III) compounds. The excited state observed following excitation of nonalkyl Cbl(III) compounds has an excited-state spectrum characteristic of Cbl(III) molecules with a weakened bond to the axial ligand (Type I). A similar excited-state spectrum is observed for adenosylcobalamin (AdoCbl) in water and ethylene glycol. The excited-state spectrum of methyl, ethyl, and n-propylcobalamin is characteristic of a Cbl(III) species with a sigma-donating alkyl anion ligand (Type II). This Type II excited-state spectrum is also observed for AdoCbl bound to glutamate mutase. The results are discussed in the context of theoretical calculations of Cbl(III) species reported in the literature and highlight the need for additional calculations exploring the influence of the alkyl ligand on the electronic structure of cobalamins.

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.

DFT analysis of co-alkyl and co-adenosyl vibrational modes in B12-cofactors

Density functional theory (DFT)-based normal mode calculations have been carried out on models for B12-cofactors to assign reported isotope-edited resonance Raman spectra, which isolate vibrations of the organo-Co group. Interpretation is straightforward for alkyl-Co derivatives, which display prominent Co-C stretching vibrational bands. DFT correctly reproduces Co-C distances and frequencies for the methyl and ethyl derivatives. However, spectra are complex for adenosyl derivatives, due to mixing of Co-C stretching with a ribose deformation coordinate and to activation of modes involving Co-C-C bending and Co-adenosyl torsion. Despite this complexity, the computed spectra provide a satisfactory re-assignment of the experimental data. Reported trends in adenosyl-cobalamin spectra upon binding to the methylmalonyl CoA mutase enzyme, as well as on subsequent binding of substrates and inhibitors, provide support for an activation mechanism involving substrate-induced deformation of the adenosyl ligand.

First Principles Study of Coenzyme B12. Crystal Packing Forces Effect on Axial Bond Lengths

In this work we analyze the structure of coenzyme B12 (AdoCbl) by means of periodic density functional theory (DFT) in order to elucidate the influence of the corrin side chains and the crystalline environment on the properties of axial bonds. The Co-Nax axial bond is very weak and its strength of less than 8 kcal/mol is four times smaller than Co-C which in solution is approximately 31 kcal/mol. The proper description of the Co-Nax distance has been problematic in previous DFT calculations and the source of disagreement between experiment and theory remained unexplained. To resolve this discrepancy, periodic DFT calculations within the Car-Parrinello molecular dynamics (CPMD) framework were carried out on three different structural models of increased complexity. The simplest model (DBI-Ado+) contains the naked corrin ring with a total of 96 atoms. The second model is the full coenzyme B12 (AdoCbl) with 209 atoms which has been taken from crystallographic analysis. To understand the extent to which the crystal packing forces influence the structural properties of AdoCbl the complete crystal consisting of four AdoCbl molecules plus 48 water molecules periodically repeated in space was analyzed (1008 atoms). The results show that the properties associated with the Co-C bond can be well reproduced using truncated models. This does not apply to the Co-Nax axial bond and the presence of the local environment appears to be essential for the correct prediction of its bond length. The most interesting outcome of the present analysis is the finding that the actual length of the Co-Nax bond (2.262 A) is largely influenced by crystal packing forces.

Characterization of Co−C Bonding in Dichlorovinylcobaloxime Complexes

This study combines theory and experiment in an examination of Co-C bonding and reductive Co-C cleavage in cobalt dichlorovinyl complexes. It is motivated by the role of dichlorovinyl complexes as intermediates in the dechlorination of trichloroethylene by cobalamin and cobalamin model complexes. A series of seven cis-1,2-dichlorovinyl(L)cobaloxime complexes were prepared (L = m- and p-substituted pyridines; cobaloxime = bis(dimethylglyoximato)cobalt). The complexes were characterized using 1H NMR, 13C NMR, cyclic voltammetry, and X-ray crystallography. Examination of the metrical parameters of the Co-C=C unit across the series shows very little change in the C=C bond length and a slight increase in the Co-C bond length with increasing electron-donating ability of the pyridine ligand. These structural changes along with electronic structure calculations indicate that Co-C pi-bonding is not important in these complexes. The stronger Co-C bonds of vinylcobaloximes compared to those of alkylcobaloximes are best explained by the higher s character at C. Changes in the reduction potential across the series indicate that the pyridine-bound form is the primary electrochemically active species. Theoretical examination of the Co-C cleavage following reduction supports the direct formation of the cis-1,2-dichlorovinyl anion and not the cis-1,2-dichlorovinyl radical.

Carbon−Cobalt Bond Distance and Bond Cleavage in One-Electron Reduced Methylcobalamin: A Failure of the Conventional DFT Method

Geometry optimizations at the HF, B3LYP, and CASSCF levels of electronic structure theory have been performed for methylcobalamin (MeCbl) model compounds in both the Co(III) (MeCbl(III)) and Co(II) (MeCbl(II)) formal oxidation states. Since the HOMO-LUMO and C-Co sigma-sigma MO gaps are significantly smaller in the MeCbl(II) compounds compared with MeCbl(III), a pseudo-Jahn Teller effect is possible. CASSCF calculations show that there is strong coupling between C-Co sigma-sigma MOs for the MeCbl(II) models leading to strong state mixing with significant total charge density transfer (approximately 0.4 e-), mainly from the C-Co sigma MO to C-Co sigma MO (approximately 0.3 e-). CASSCF(9:7) calculations show that the strong state mixing leads to an increase in the C-Co bond length for MeCbl(II) model compounds from 1.969 A (DFT and HF calculations) to 2.164 A in the base-on MeCbl(II) model and from 1.938 A to 2.144 A in the base-off MeCbl(II) model. Concomitantly, the Co-N axial bond length increases from 2.121 A (DFT) to 2.344 A in the CASSCF calculation. This coupling interaction between states can be used to explain the much lower Co-C bond dissociation enthalpy and much faster bond cleavage rate for the one-electron reduced methylcobalamin radical anion compared to MeCbl(III). It may also be important for axial bond distances in other Co(II) compounds.

Reaction mechanism of porphyrin metallation studied by theoretical methods

We have studied the reaction mechanism for the insertion of Mg2+ and Fe2+ into a porphyrin ring with density functional calculations with large basis set and including solvation, zero-point and thermal effects. We have followed the reaction from the outer-sphere complex, in which the metal is coordinated with six water molecules and the porphyrin is doubly protonated, until the metal ion is inserted into the deprotonated porphyrin ring with only one water ligand remaining. This reaction involves the stepwise displacement of five water molecules and the removal of two protons from the porphyrin ring. In addition, a step seems to be necessary in which a porphyrin pyrrolenine nitrogen atom changes its interaction from a hydrogen bond to a metal-bound solvent molecule to a direct coordination to the metal ion. If the protons are taken up by a neutral imidazole molecule, the deprotonation reactions are exothermic with minimal barriers. However, with a water molecule as an acceptor, they are endothermic. The ligand exchange reactions were approximately thermoneutral (+/-20 kJ mol(-1), with one exception) with barriers of up to 72 kJ mol(-1) for Mg and 51 kJ mol(-1) for Fe. For Mg, the highest barrier was found for the formation of the first bond to the porphyrin ring. For Fe, a higher barrier was found for the formation of the second bond to the porphyrin ring, but this barrier is probably lower in solution. No evidence was found for an initial pre-equilibrium between a planar and a distorted porphyrin ring. Instead, the porphyrin becomes more and more distorted as the number of metal-porphyrin bonds increase (by up to 191 kJ mol(-1)). This strain is released when the porphyrin becomes deprotonated and the metal moves into the ring plane. Implications of these findings for the chelatase enzymes are discussed.

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.

Electroreduction of a Series of Alkylcobalamins: Mechanism of Stepwise Reductive Cleavage of the Co−C Bond

The electrochemical (EC) reduction mechanism of methylcobalamin (Me-Cbl) in a mixed DMF/MeOH solvent in 0.2 M tetrabutylammonium fluoroborate electrolyte was studied as a function of temperature and solvent ratio vs a nonaqueous Ag/AgCl/Cl(-) reference electrode. Double-potential-step chronoamperometry allowed the rate constant of the subsequent homogeneous reaction to be measured over the temperature range from 0 to -80 degrees C in 40:60 and 50:50 DMF:MeOH ratios. Activation enthalpies are 5.8 +/- 0.5 and 7.6 +/- 0.3 kcal/mol in the 40:60 and 50:50 mixtures of DMF/MeOH, respectively. Digital simulation and curve-fitting for an EC mechanism using a predetermined homogeneous rate constant of 5.5 x 10(3) s(-1) give E degrees' = -1.466 V, k degrees = 0.016 cm/s, and alpha = 0.77 at 20 degrees C for a quasi-reversible electrode process. Digital simulation of the results of Lexa and Savéant (J. Am. Chem. Soc. 1978, 100, 3220-3222) shows that the mechanism is a series of stepwise homogeneous equilibrium processes with an irreversible step following the initial electron transfer (ET) and allows estimation of the equilibrium and rate constants of these reactions. An electron coupling matrix element of H(kA) = (4.7 +/- 1.1) x 10(-4) eV ( approximately 46 J/mol) is calculated for the nonadiabatic ET step for reduction to the radical anion. A reversible bond dissociation enthalpy for homolytic cleavage of Me-Cbl is calculated as 31 +/- 2 kcal/mol. The voltammetry of the ethyl-, n-propyl-, n-butyl-, isobutyl-, and adenosyl-substituted cobalamin was studied, and estimated reversible redox potentials were correlated with Co-C bond distances as determined by DFT (B3LYP/ LANL2DZ) calculations.

Electronic reorganization: Origin of sigma trans promotion effect

Binding of two ligands trans to each other by some transition metal complexes may be cooperative [Khoroshun et al., Mol Phys 2002, 100, 523]. Several interesting consequent effects include (i) inverse relationship between bond strength and binding affinity; (ii) smaller coordination barriers to formation of weaker bonds; (iii) enhancement of Lewis acidity with increased number of ligands. We describe a simple model, sigma trans promotion effect (TPE), which considers electronic reorganization between two Lewis structures, and predicts the above-mentioned effects. The applied result of present study is the unified perspective on several facts of heme chemistry. Particularly, we reiterate an important but often overlooked notion, developed previously within the spin pairing model [Drago and Corden, Acc Chem Res 1980, 13, 353], that, in hemoproteins, the proximal histidine and the distal ligand such as O2 or CO cooperate in promoting electronic reorganization. As a result, depopulation of dz2 orbital upon ligand binding contributes to the phenomenon of hemoglobin cooperativity. The presented density functional (B3LYP) calculations on realistic models, the processes of carbon monoxide binding by Fe(II) porphyrins and dinitrogen binding by triamido/triamidoamine Mo(III) complexes, particularly the evaluation of the coordination barriers due to spin-state change by location of the minima on seams of crossing, support the TPE model predictions. From a broader theoretical perspective, the present study would hopefully stimulate the development of much needed frameworks and tools for facile comparisons of wave functions and their properties between different geometries, species, and electronic states. Advancement of practical wave function comparisons may yield fresh qualitative perspectives on chemical reactivity, and promote better understanding of related concepts such as electronic reorganization.

How the Co-C bond is cleaved in coenzyme B12 enzymes: a theoretical study

The homolytic cleavage of the organometallic Co-C bond in vitamin B12-dependent enzymes is accelerated by a factor of approximately 10(12) in the protein compared to that of the isolated cofactor in aqueous solution. To understand this much debated effect, we have studied the Co-C bond cleavage in the enzyme glutamate mutase with combined quantum and molecular mechanics methods. We show that the calculated bond dissociation energy (BDE) of the Co-C bond in adenosyl cobalamin is reduced by 135 kJ/mol in the enzyme. This catalytic effect can be divided into four terms. First, the adenosine radical is kept within 4.2 angstroms of the Co ion in the enzyme, which decreases the BDE by 20 kJ/mol. Second, the surrounding enzyme stabilizes the dissociated state by 42 kJ/mol using electrostatic and van der Waals interactions. Third, the protein itself is stabilized by 11 kJ/mol in the dissociated state. Finally, the coenzyme is geometrically distorted by the protein, and this distortion is 61 kJ/mol larger in the Co(III) state. This deformation of the coenzyme is caused mainly by steric interactions, and it is especially the ribose moiety and the Co-C5'-C4' angle that are distorted. Without the polar ribose group, the catalytic effect is much smaller, e.g. only 42 kJ/mol for methyl cobalamin. The deformation of the coenzyme is caused mainly by the substrate, a side chain of the coenzyme itself, and a few residues around the adenosine part of the coenzyme.

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.

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.

The structure of sitting-atop complexes of metalloporphyrins studied by theoretical methods

The metallation of tetrapyrroles is believed to proceed via a sitting-atop (SAT) complex, in which some of the pyrrole nitrogen atoms are still protonated and the metal ion resides above the ring plane. No crystal structure of such a complex has been presented, but NMR and extended X-ray absorption fine structure (EXAFS) data has been reported for Cu(2+) in acetonitrile. We have used density functional calculations to obtain reasonable models for SAT complexes of porphyrins with Mg(2+), Fe(2+), and Cu(2+). The results show that there are many possible SAT complexes with 1-5 solvent molecules, one or two metal ions, and cis or trans protonation of the porphyrin ring. Many of these have similar energies and their relative stabilities vary with the metal ion. A complex with two cis pyrrolenine nitrogens atoms and 2-4 solvent molecules coordinated to Cu(2+) fits the NMR and EXAFS data best. However, we cannot fully exclude the possibility that what is observed is rather a mixture of a doubly protonated porphyrin and the copper porphyrin. Mg(2+) has a lower affinity for porphyrin and stronger affinity for water, so a complex with five water molecules and only one bond to porphyrin seems to be most stable. For Fe(2+), a cis structure with two first-sphere water molecules and four interactions to the porphyrin seems to be most likely.

Structure-energy relations in methylcobalamin with and without bound axial base

The properties of the Co-C bond in methylcobalamin (MeCbl) are analyzed by means of first-principles molecular dynamics. The optimized structure is in very good agreement with experiments, reproducing the bent-up deformation of the corrin ring as well as the metal-ligand bond distances. The analysis of the binding energies, bond orders, and vibrational stretching frequencies shows that the axial base slightly weakens the Co-C bond (by 4%), while the alkyl ligand substantially reinforces the Co-axial base bond (by 90%). These findings support several experiments and provide insight into the conversion between the base-on and base-off forms of the MeCbl cofactor.

Theoretical Investigations into the Intermediacy of Chlorinated Vinylcobalamins in the Reductive Dehalogenation of Chlorinated Ethylenes

The reductive dehalogenation of perchloroethylene and trichloroethylene by vitamin B(12) produces approximately 95% (Z)-dichloroethylene (DCE) and small amounts of (E)-DCE and 1,1-DCE, which are further reduced to ethylene and ethane. Chloroacetylene and acetylene have been detected as intermediates, but not dichloroacetylene. Organocobalamins (RCbls) have been proposed to be intermediates in this process. Density functional theory based approaches were employed to investigate the properties of chlorinated vinylcobalamins and chlorinated vinyl radicals. They reveal that all vinyl radicals studied have reduction potentials more positive (E degrees >or= -0.49) than that of the Co(II)/Co(I) couple of B(12) (E degrees = -0.61 V), indicating that any (chlorinated) vinyl radicals formed in the reductive dehalogenation process should be reduced to the corresponding anions by cob(I)alamin in competition with their combination with Co(II) to yield the corresponding vinylcobalamins. The computed Co-C homolytic bond dissociation enthalpies (BDEs) of the latter complexes range from 33.4 to 45.8 kcal/mol. The substituent effects on the BDEs are affected by the stabilities of the vinyl radicals as well as steric interactions between (Z)-chloro substituents and the corrin ring. The calculated E degrees values of the cobalamin models were within approximately 200 mV of one another since electron attachment is to a corrin ring pi-orbital, whose energy is relatively unaffected by chloride substitution of the vinyl ligand, and all were >500 mV more negative than that of the Co(II)/Co(I) couple of B(12). Reduction of the base-off forms of vinyl- and chlorovinylcobalamin models also involves the corrin pi* orbital, but reduction of the base-off dichlorovinyl- and trichlorovinylcobalamin models occurs with electron attachment to the sigma(Co)(-)(C*) orbital, yielding calculated E degrees values more positive than that of the calculated Co(II)/Co(I) couple of B(12). Thus, cob(I)alamin is expected to reduce these base-off vinyl-Cbls. Heterolytic cleavage of the Co-C bonds is much more favorable than homolysis (>21 kcal/mol) and is significantly more exergonic when coupled to chloride elimination.

Catalytic Roles of Active-Site Amino Acid Residues of Coenzyme B12-Dependent Diol Dehydratase: Protonation State of Histidine and Pull Effect of Glutamate

The hydrogen abstraction and the OH migration processes catalyzed by diol dehydratase are discussed by means of a quantum mechanical/molecular mechanical method. To evaluate the push effect of His143 and the pull effect of Glu170, we considered three kinds of whole-enzyme model, the protonated and two unprotonated His143 models. A calculated activation energy for the hydrogen abstraction by the adenosyl radical is 15.6 (13.6) kcal/mol in the protonated (unprotonated) His143 model. QM/MM calculational results show that the mechanism of the OH migration is significantly changed by the protonation of His143. In the protonated His143 model, the OH group migration triggered by the full proton donation from the imidazolium to the migrating OH group occurs by a stepwise OH abstraction/re-addition process in which the water production reduces the barrier for the C-O bond cleavage. On the other hand, the OH migration in the unprotonated His143 model proceeds in a concerted manner, as we previously proposed using a simple model including only K+ ion and substrate. The latter mechanism seems to be kinetically more favorable from the calculated energy profiles and is consistent with experimental results. The activation barrier of the OH group migration step is only 1.6 kcal/mol reduced by the hydrogen-bonding interaction between the O2 of the substrate and unprotonated His143. Thus, it is predicted that His143 is not protonated, and therefore the main active-site amino acid residue that lowers the energy of the transition state for the OH group migration is determined to be Glu170.