Extended X-ray Absorption Fine Structure Analysis of Coenzyme B12 Bound to Methylmalonyl-Coenzyme A Mutase Using Global Mapping Techniques

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.

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.

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.

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.

Spectroscopic and computational studies on the adenosylcobalamin-dependent methylmalonyl-CoA mutase: evaluation of enzymatic contributions to Co-C bond activation in the Co3+ ground state

Methylmalonyl-CoA mutase (MMCM) is an enzyme that utilizes the adenosylcobalamin (AdoCbl) cofactor to catalyze the rearrangement of methylmalonyl-CoA to succinyl-CoA. Despite many years of dedicated research, the mechanism by which MMCM and related AdoCbl-dependent enzymes accelerate the rate for homolytic cleavage of the cofactor's Co-C bond by approximately 12 orders of magnitude while avoiding potentially harmful side reactions remains one of the greatest subjects of debate among B(12) researchers. In this study, we have employed electronic absorption (Abs) and magnetic circular dichroism (MCD) spectroscopic techniques to probe cofactor/enzyme active site interactions in the Co(3+)Cbl "ground" state for MMCM reconstituted with both the native cofactor AdoCbl and its derivative methylcobalamin (MeCbl). In both cases, Abs and MCD spectra of the free and enzyme-bound cofactor are very similar, indicating that replacement of the intramolecular base 5,6-dimethylbenzimidazole (DMB) by a histidine residue from the enzyme active site has insignificant effects on the cofactor's electronic properties. Likewise, spectral perturbations associated with substrate (analogue) binding to holo-MMCM are minor, arguing against substrate-induced enzymatic Co-C bond activation. As compared to the AdoCbl data, however, Abs and MCD spectral changes for the sterically less constrained MeCbl cofactor upon binding to MMCM and treatment of holoenzyme with substrate (analogues) are much more substantial. Analysis of these changes within the framework of time-dependent density functional theory calculations provides uniquely detailed insight into the structural distortions imposed on the cofactor as the enzyme progresses through the reaction cycle. Together, our results indicate that, although the enzyme may serve to activate the cofactor in its Co(3+)Cbl ground state to a small degree, the dominant contribution to the enzymatic Co-C bond activation presumably comes through stabilization of the Co(2+)Cbl/Ado. post-homolysis products.

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.

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.

Coenzyme B(12) dependent glutamate mutase

The crystal structure of glutamate mutase with bound coenzyme B(12) suggests a radical shuttling mechanism within the active site of the enzyme. Quantum chemical calculations of the rearrangement in combination with kinetic and mutational studies suggest the catalytic mechanism of this enzyme to proceed via a fragmentation/recombination sequence with intermediates stabilized by partial protonation/deprotonation. Crucial residues in the active site have been identified. Solution structure studies indicate the mechanism of B(12) binding to the apoenzyme.

Adenosylcobalamin-dependent isomerases: new insights into structure and mechanism

Adenosylcobalamin-dependent isomerases catalyze a variety of chemically difficult 1,2-rearrangements that proceed through a mechanism involving free radical intermediates. These radicals are initially generated by homolysis of the cobalt-carbon bond of the coenzyme. Recently, the crystal structures of several of these enzymes have been solved, revealing two modes of coenzyme binding and highlighting the role of the protein in controlling the rearrangement of reactive substrate radical intermediates. Complementary data from kinetic, spectroscopic and theoretical studies have produced insights into the mechanism by which substrate radicals are generated at the active site, and the pathways by which they rearrange.

Assessment of the existence of hyper-long axial Co(II)-N bonds in cobinamide B(12) models by using electron paramagnetic resonance spectroscopy

Protein control of cobalt-axial nitrogen ligand bond length has been proposed to modulate the reactivity of vitamin B(12) coenzyme during the catalytic cycle of B(12)-dependent enzymes. In particular, hyper-long Co-N bonds may favor homolytic cleavage of the trans-cobalt-carbon bond in the coenzyme. X-ray crystallographic studies point to hyper-long bonds in two B(12) holoenzymes; however, mixed redox and ligand states in the crystals thwart clear conclusions. Since EPR theory predicts an increase in Co(II) hyperfine splitting as donation from the axial N-donor ligand decreases, EPR spectroscopy could clarify the X-ray results. However, the theory is apparently undermined by the similar splitting reported for the 2-picoline (2-pic) and pyridine (py) adducts of Co(II) cobinamide (Co(II)Cbi(+)), adducts thought to have long and normal Co-N axial bond lengths, respectively. Cobinamides, with the B(12) 5,6-dimethylbenzimidazole loop removed, are excellent B(12) models. We studied Co(II)Cbi(+) adducts of unhindered 4-substituted pyridines (4-X-py's) in ethylene glycol to separate orbital size effects from Co-N axial distance effects on these splittings. The linear increase in splitting with the decrease in 4-X-py basicity found is consistent with the theoretically predicted increase in unpaired electron spin density as axial N lone pair donation to Co(II) decreases. No adduct (and hence no hyper-long Co(II)-N axial bond) was formed even by 8 M 2-pic, if the 2-pic was purified by a novel Co(III)-affinity distillation procedure designed to remove trace nitrogenous ligand impurities present in 2-pic distilled in the regular manner. Adducts formed by impurities in 2-pic and other hindered pyridines misled previous investigators into attributing results to adducts with long Co-N bonds. We find that many 2-substituted py's known to form adducts with simple synthetic Co models do not bind Co(II)Cbi(+). Thus, the equatorial corrin ring sterically impedes binding, making Co(II)Cbi(+) a highly selective binding agent for unhindered sp(2) N-donor ligands. Our results resolve the apparent conflict between EPR experiment and theory. The reported Co(II) hyperfine splitting of the enzyme-bound cofactor in five B(12) enzymes is similar to that of the relevant free cofactor. The most reasonable interpretation of this similarity is that the Co-N axial bond of the bound cofactor is not hyper-long in any of the five cases.

Glutamate mutase from Clostridium cochlearium: the structure of a coenzyme B12-dependent enzyme provides new mechanistic insights

BACKGROUND

Glutamate mutase (Glm) equilibrates (S)-glutamate with (2S,3S)-3-methylaspartate. Catalysis proceeds with the homolytic cleavage of the organometallic bond of the cofactor to yield a 5'-desoxyadenosyl radical. This radical then abstracts a hydrogen atom from the protein-bound substrate to initiate the rearrangement reaction. Glm from Clostridium cochlearium is a heterotetrameric molecule consisting of two sigma and two epsilon polypeptide chains.

RESULTS

We have determined the crystal structures of inactive recombinant Glm reconstituted with either cyanocobalamin or methylcobalamin. The molecule shows close similarity to the structure of methylmalonyl CoA mutase (MCM), despite poor sequence similarity between its catalytic epsilon subunit and the corresponding TIM-barrel domain of MCM. Each of the two independent B12 cofactor molecules is associated with a substrate-binding site, which was found to be occupied by a (2S,3S)-tartrate ion. A 1:1 mixture of cofactors with cobalt in oxidation states II and III was observed in both crystal structures of inactive Glm.

CONCLUSIONS

The long axial cobalt-nitrogen bond first observed in the structure of MCM appears to result from a contribution of the species without upper ligand. The tight binding of the tartrate ion conforms to the requirements of tight control of the reactive intermediates and suggests how the enzyme might use the substrate-binding energy to initiate cleavage of the cobalt-carbon bond. The cofactor does not appear to have a participating role during the radical rearrangement reaction.

Coenzyme B(12) Axial-Base Chemical Precedent Studies. Adenosylcobinamide Plus Sterically Hindered Axial-Base Co-C Bond Cleavage Product and Kinetic Studies: Evidence for the Dominance of Axial-Base Transition-State Effects and for Co-N(Axial-Base) Distance-Dependent, Competing sigma and pi Effects

Adenosylcobinamide (AdoCbi(+)) plus the sterically hindered bases 1,2-dimethylimidazole, 2-methylpyridine, and 2,6-dimethylpyridine, as well as control experiments with imidazolate and 4-methylimidazolate, have been investigated to provide chemical precedent for the benzimidazole base-off, protein histidine imidazole base-on form of adenosylcobalamin (AdoCbl, also coenzyme B(12)). This imidazole base-on form of AdoCbl was observed in the recent X-ray crystallographic structural study of methylmalonyl-CoA (MMCoA) mutase; of interest to the present work is the fact that MMCoA mutase contains a long, ca. 2.5 Å, Co-N(imidazole) axial bond, at least in the enzyme's crystallographically characterized Co(II)/Co(III) state and conformation. In the present studies, upper limits for the axial-base binding K(assoc) parameters to form [AdoCbi.bulky base](+) BF(4)(-) have been obtained; these thermodynamic studies reveal that sterically hindered bases do not bind detectably to AdoCbi(+) in the ground state, which results in negligible ground-state free-energy stabilization via the formation of [AdoCbi.bulky base](+). The sterically hindered bases do, however, bind to Co(II)Cbi(+), a good energetic model of the [Ado. - - -.CoCbi](+) homolysis transition state. Kinetic studies demonstrate that the sterically hindered bases are involved in the rate-determining step of Co-C bond homolysis, accelerating it by 200-fold; hence, Co-C cleavage does occur via the low-level and otherwise nondetectable amount of [AdoCbi.bulky base](+) formed in solution. Product studies reveal (i) that both Co-C heterolysis and homolysis occur, and (ii) that there is no simple correlation between the ratio of Co-C heterolysis to homolysis and the Co-N(axial-base) bond length. Overall, the results provide strong evidence for the dominance of axial-base transition-state effects on Co-C bond cleavage, and reveal a subtle interplay of sigma and pi effects as a function of the Co-N(axial-base) bond length.

Binding of Cob(II)alamin to the adenosylcobalamin-dependent ribonucleotide reductase from Lactobacillus leichmannii. Identification of dimethylbenzimidazole as the axial ligand

The ribonucleoside triphosphate reductase (RTPR) from Lactobacillus leichmannii catalyzes the reduction of nucleoside 5'-triphosphates to 2'-deoxynucleoside 5'-triphosphates and uses coenzyme B12, adenosylcobalamin (AdoCbl), as a cofactor. Use of a mechanism-based inhibitor, 2'-deoxy-2'-methylenecytidine 5'-triphosphate, and isotopically labeled RTPR and AdoCbl in conjunction with EPR spectroscopy has allowed identification of the lower axial ligand of cob(II)alamin when bound to RTPR. In common with the AdoCbl-dependent enzymes catalyzing irreversible heteroatom migrations and in contrast to the enzymes catalyzing reversible carbon skeleton rearrangements, the dimethylbenzimidazole moiety of the cofactor is not displaced by a protein histidine upon binding to RTPR.