Thermodynamic and kinetic characterization of Co-C bond homolysis catalyzed by coenzyme B(12)-dependent methylmalonyl-CoA mutase

Reactivity Tracking of an Enzyme Progress Coordinate

The reactivity of individual solvent-coupled protein configurations is used to track and resolve the progress coordinate for the core reaction sequence of substrate radical rearrangement and hydrogen atom transfer in the ethanolamine ammonia-lyase (EAL) enzyme from Salmonella enterica. The first-order decay of the substrate radical intermediate is the monitored reaction. Heterogeneous confinement from sucrose hydrates in the mesophase solvent surrounding the cryotrapped protein introduces distributed kinetics in the non-native decay of the substrate radical pair capture substate, which arise from an ensemble of configurational microstates. Reaction rates increase by >103-fold across the distribution to approach that for the native enabled substate for radical rearrangement, which reacts with monotonic kinetics. The native progress coordinate thus involves a collapse of the configuration space to generate optimized reactivity. Reactivity tracking reveals fundamental features of solvent-protein-reaction configurational coupling and leads to a model that refines the ensemble paradigm of enzyme catalysis for strongly adiabatic chemical steps.

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.

Electronic structure studies of free and enzyme-bound B12 species by magnetic circular dichroism and complementary spectroscopic techniques

Electronic absorption (Abs) and circular dichroism (CD) spectroscopic techniques have been used successfully for over half a century in studies of free and enzyme-bound B₁₂ species. More recently, magnetic circular dichroism (MCD) spectroscopy and other complementary techniques have provided an increasingly detailed understanding of the electronic structure of cobalamins. While CD spectroscopy measures the difference in the absorption of left- and right-circularly polarized light, MCD spectroscopy adds the application of a magnetic field parallel to the direction of light propagation. Transitions that are formally forbidden according to the Abs and CD selection rules, such as ligand field (or d→d) transitions, can gain MCD intensity through spin-orbit coupling. As such, MCD spectroscopy provides a uniquely sensitive probe of the different binding modes, Co oxidation states, and axial ligand environments of B₁₂ species in enzyme active sites, and thus the distinct reactivities displayed by these species. This chapter summarizes representative MCD studies of free and enzyme-bound B₁₂ species, including those present in adenosyltransferases, isomerases, and reductive dehalogenases. Complementary spectroscopic and computational data are also presented and discussed where appropriate.

Photolytic properties of B12-dependent enzymes: A theoretical perspective

The biologically active vitamin B₁₂ derivates, methylcobalamin (MeCbl) and adenosylcobalamin (AdoCbl), are ubiquitous organometallic cofactors. In addition to their key roles in enzymatic catalysis, B₁₂ cofactors have complex photolytic properties which have been the target of experimental and theoretical studies. With the recent discovery of B₁₂-dependent photoreceptors, there is an increased need to elucidate the underlying photochemical mechanisms of these systems. This book chapter summarizes the photolytic properties of MeCbl- and AdoCbl-dependent enzymes with particular emphasis on the effect of the environment of the cofactor on the excited state processes. These systems include isolated MeCbl and AdoCbl as well as the enzymes, ethanolamine ammonia-lyase (EAL), glutamate mutase (GLM), methionine synthase (MetH), and photoreceptor CarH. Central to determining the photodissociation mechanism of each system is the analysis of the lowest singlet excited state (S₁) potential energy surface (PES). Time-dependent density functional theory (TD-DFT), employing BP86/TZVPP, is widely used to construct such PESs. Regardless of the environment, the topology of the S₁ PES of AdoCbl or MeCbl is marked by characteristic features, namely the metal-to-ligand charge transfer (MLCT) and ligand field (LF) regions. Conversely, the relative energetics of these electronic states are affected by the environment. Applications and outlooks for Cbl photochemistry are also discussed.

Human B12-dependent enzymes: Methionine synthase and Methylmalonyl-CoA mutase

Humans have only two known cobalamin or B₁₂-dependent enzymes: cytoplasmic methionine synthase and mitochondrial methylmalonyl-CoA mutase. A complex intracellular B₁₂ trafficking pathway, comprising a multitude of chaperones, process and deliver cobalamin to the two target enzymes. Methionine synthase catalyzes the transfer of a methyl group from N⁵-methytetrahydrofolate to homocysteine, generating tetrahydrofolate and methionine. Cobalamin serves as an intermediate methyl group carrier and cycles between methylcobalamin and cob(I)alamin. Methylmalonyl-CoA mutase uses the 5'-deoxyadenosylcobalamin form of the cofactor and catalyzes the 1,2 rearrangement of methylmalonyl-CoA to succinyl-CoA. Two chaperones, CblA (or MMAA) and CblB (or MMAB, also known as adenosyltransferase), serve the mutase and ensure that the fidelity of the cofactor loading and unloading processes is maintained. This chapter focuses on assays for purifying and measuring the activities of methionine synthase and methylmalonyl-CoA mutase.

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.

Resonance Raman spectroscopic study of the interaction between Co(II)rrinoids and the ATP:corrinoid adenosyltransferase PduO from Lactobacillus reuteri

The human-type ATP:corrinoid adenosyltransferase PduO from Lactobacillus reuteri (LrPduO) catalyzes the adenosylation of Co(II)rrinoids to generate adenosylcobalamin (AdoCbl) or adenosylcobinamide (AdoCbi(+)). This process requires the formation of "supernucleophilic" Co(I)rrinoid intermediates in the enzyme active site which are properly positioned to abstract the adeonsyl moiety from co-substrate ATP. Previous magnetic circular dichroism (MCD) spectroscopic and X-ray crystallographic analyses revealed that LrPduO achieves the thermodynamically challenging reduction of Co(II)rrinoids by displacing the axial ligand with a non-coordinating phenylalanine residue to produce a four-coordinate species. However, relatively little is currently known about the interaction between the tetradentate equatorial ligand of Co(II)rrinoids (the corrin ring) and the enzyme active site. To address this issue, we have collected resonance Raman (rR) data of Co(II)rrinoids free in solution and bound to the LrPduO active site. The relevant resonance-enhanced vibrational features of the free Co(II)rrinoids are assigned on the basis of rR intensity calculations using density functional theory to establish a suitable framework for interpreting rR spectral changes that occur upon Co(II)rrinoid binding to the LrPduO/ATP complex in terms of structural perturbations of the corrin ring. To complement our rR data, we have also obtained MCD spectra of Co(II)rrinoids bound to LrPduO complexed with the ATP analogue UTP. Collectively, our results provide compelling evidence that in the LrPduO active site, the corrin ring of Co(II)rrinoids is firmly locked in place by several amino acid side chains so as to facilitate the dissociation of the axial ligand.

Photolytic properties of cobalamins: a theoretical perspective

This Perspective Article highlights recent theoretical developments, and summarizes the current understanding of the photolytic properties of cobalamins from a computational point of view. The primary focus is on two alkyl cobalamins, methylcobalamin (MeCbl) and adenosylcobalamin (AdoCbl), as well as two non-alkyl cobalamins, cyanocobalamin (CNCbl) and hydroxocobalamin (HOCbl). Photolysis of alkyl cobalamins involves low-lying singlet excited states where photodissociation of the Co-C bond leads to formation of singlet-born alkyl/cob(ii)alamin radical pairs (RPs). Potential energy surfaces (PESs) associated with cobalamin low-lying excited states as functions of both axial bonds, provide the most reliable tool for initial analysis of their photochemical and photophysical properties. Due to the complexity, and size limitations associated with the cobalamins, the primary method for calculating ground state properties is density functional theory (DFT), while time-dependent DFT (TD-DFT) is used for electronically excited states. For alkyl cobalamins, energy pathways on the lowest singlet surface, connecting metal-to-ligand charge transfer (MLCT) and ligand field (LF) minima, can be associated with photo-homolysis of the Co-C bond observed experimentally. Additionally, energy pathways between minima and seams associated with crossing of S1/S0 surfaces, are the most efficient for internal conversion (IC) to the ground state. Depending on the specific cobalamin, such IC may involve simultaneous elongation of both axial bonds (CNCbl), or detachment of axial base followed by corrin ring distortion (MeCbl). The possibility of intersystem crossing, and the formation of triplet RPs is also discussed based on Landau-Zener theory.

Resolution and Characterization of Chemical Steps in Enzyme Catalytic Sequences by Using Low-Temperature and Time-Resolved, Full-Spectrum EPR Spectroscopy in Fluid Cryosolvent and Frozen Solution Systems

Approaches to the resolution and characterization of individual chemical steps in enzyme catalytic sequences, by using temperatures in the cryogenic range of 190-250 K, and kinetics measured by time-resolved, full-spectrum electron paramagnetic resonance spectroscopy in fluid cryosolvent and frozen solution systems, are described. The preparation and performance of the adenosylcobalamin-dependent ethanolamine ammonia-lyase enzyme from Salmonella typhimurium in the two systems exemplifies the biochemical and spectroscopic methods. General advantages of low-temperature studies are (1) slowing of reaction steps, so that measurements can be made by using straightforward T-step kinetic methods and commercial instrumentation, (2) resolution of individual reaction steps, so that first-order kinetic analysis can be applied, and (3) accumulation of intermediates that are not detectable at room temperatures. The broad temperature range from room temperature to 190 K encompasses three regimes: (1) temperature-independent mean free energy surface (corresponding to native behavior); (2) the narrow temperature region of a glass-like transition in the protein, over which the free energy surface changes, revealing dependence of the native reaction on collective protein/solvent motions; and (3) the temperature range below the glass transition region, for which persistent reaction corresponds to nonnative, alternative reaction pathways, in the vicinity of the native configurational envelope. Representative outcomes of low-temperature kinetics studies are portrayed on Eyring and free energy surface (landscape) plots, and guidelines for interpretations are presented.

Unprecedented Mechanism Employed by the Salmonella enterica EutT ATP:Co(I)rrinoid Adenosyltransferase Precludes Adenosylation of Incomplete Co(II)rrinoids

Three distinct families of ATP:corrinoid adenosyltransferases (ACATs) exist that are capable of converting vitamin B12 derivatives into coenzyme B12 by catalyzing the thermodynamically challenging reduction of Co(II) rrinoids to form "supernucleophilic" Co(I) intermediates. While the structures and mechanisms of two of the ACAT families have been studied extensively, little is known about the EutT enzymes beyond the fact that they exhibit a unique requirement for a divalent metal cofactor for enzymatic activity. In this study we have obtained compelling evidence that EutT converts cob(II)alamin into an effectively four-coordinate Co(II) species so as to facilitate Co(II)→Co(I) reduction. Intriguingly, EutT fails to promote axial ligand dissociation from the substrate analogue cob(II)inamide, a natural precursor of cob(II)alamin. This unique substrate specificity of EutT has important physiological implications.

The entropic contributions in vitamin B12 enzymes still reflect the electrostatic paradigm

The catalytic power of enzymes containing coenzyme B12 has been, in some respects, the "last bastion" for the strain hypothesis. Our previous study of this system established by a careful sampling that the major part of the catalytic effect is due to the electrostatic interaction between the ribose of the ado group and the protein and that the strain contribution is very small. This finding has not been sufficiently appreciated due to misunderstandings of the power of the empirical valence bond (EVB) calculations and the need of sufficient sampling. Furthermore, some interesting new experiments point toward entropic effects as the source of the catalytic power, casting doubt on the validity of the electrostatic idea, at least, in the case of B12 enzymes. Here, we focus on the observation of the entropic effects and on analyzing their origin. We clarify that our EVB approach evaluates free energies rather than enthalpies and demonstrate by using the restraint release (RR) approach that the observed entropic contribution to the activation barrier is of electrostatic origin. Our study illustrates the power of the RR approach by evaluating the entropic contributions to catalysis and provides further support to our paradigm for the origin of the catalytic power of B12 enzymes. Overall, our study provides major support to our electrostatic preorganization idea and also highlights the basic requirements from ab initio quantum mechanics/molecular mechanics calculations of activation free energies of enzymatic reactions.

Methylmalonic acid quantified in dried blood spots provides a precise, valid, and stable measure of functional vitamin B-12 status in healthy women

Methylmalonic acid (MMA) is a sensitive and specific functional biomarker of vitamin B-12 status, commonly assessed in plasma or serum. Dried blood spots (DBSs) allow simpler and more cost-efficient blood sampling than plasma. To facilitate convenient testing for vitamin B-12 deficiency in large-scale surveys and in population groups from remote areas, we developed a method for MMA quantification in DBSs and tested its applicability as well as the long-term stability of MMA in DBSs at various temperatures. MMA was extracted from an 8-mm DBS punch with water:methanol (95:5, v:v) and methyl-d3-malonic acid as the internal standard. After sample cleanup by ultrafiltration and hexane extraction, MMA was quantified by using reversed-phase LC-tandem mass spectrometry. Extraction conditions were optimized to maximize the detection signal and achieve DBS extract concentrations above the lowest limit of quantification (signal-to-noise ratio ≥ 10) of 10 nmol/L. Recovery was between 93% and 96%. Intra- and interassay variation (CV%) for DBS MMA was 0.49% and 2.3%, respectively. Calibrators showed linearity (R(2) = 0.998) between 10 and 10,000 nmol/L. In 94 healthy women, MMA concentrations in DBS extract (min-max: 10.2-80.5 nmol/L) and plasma (min-max: 68-950 nmol/L) were correlated (ρ = 0.90) (P < 0.001). MMA concentrations in DBSs were stable at room temperature for 1 wk, in the refrigerator for 8 wk, and at -80°C for at least 1 y. This simple and robust method allows quantification of MMA in DBSs of healthy individuals. The linear relation between plasma and DBS MMA suggests that DBS MMA could predict plasma MMA, the current reference indicator for functional vitamin B-12 deficiency. With the advantages of minimally invasive specimen collection and no need for laborious blood processing steps, this method has the potential to be a reliable, convenient, and field-applicable alternative for assessment of vitamin B-12 status.

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.

A mechanochemical switch to control radical intermediates

B₁₂-dependent enzymes employ radical species with exceptional prowess to catalyze some of the most chemically challenging, thermodynamically unfavorable reactions. However, dealing with highly reactive intermediates is an extremely demanding task, requiring sophisticated control strategies to prevent unwanted side reactions. Using hybrid quantum mechanical/molecular mechanical simulations, we follow the full catalytic cycle of an AdoB₁₂-dependent enzyme and present the details of a mechanism that utilizes a highly effective mechanochemical switch. When the switch is "off", the 5'-deoxyadenosyl radical moiety is stabilized by releasing the internal strain of an enzyme-imposed conformation. Turning the switch "on," the enzyme environment becomes the driving force to impose a distinct conformation of the 5'-deoxyadenosyl radical to avoid deleterious radical transfer. This mechanochemical switch illustrates the elaborate way in which enzymes attain selectivity of extremely chemically challenging reactions.

Entropic origin of cobalt-carbon bond cleavage catalysis in adenosylcobalamin-dependent ethanolamine ammonia-lyase

Adenosylcobalamin-dependent enzymes accelerate the cleavage of the cobalt-carbon (Co-C) bond of the bound coenzyme by >10(10)-fold. The cleavage-generated 5'-deoxyadenosyl radical initiates the catalytic cycle by abstracting a hydrogen atom from substrate. Kinetic coupling of the Co-C bond cleavage and hydrogen-atom-transfer steps at ambient temperatures has interfered with past experimental attempts to directly address the factors that govern Co-C bond cleavage catalysis. Here, we use time-resolved, full-spectrum electron paramagnetic resonance spectroscopy, with temperature-step reaction initiation, starting from the enzyme-coenzyme-substrate ternary complex and (2)H-labeled substrate, to study radical pair generation in ethanolamine ammonia-lyase from Salmonella typhimurium at 234-248 K in a dimethylsulfoxide/water cryosolvent system. The monoexponential kinetics of formation of the (2)H- and (1)H-substituted substrate radicals are the same, indicating that Co-C bond cleavage rate-limits radical pair formation. Analysis of the kinetics by using a linear, three-state model allows extraction of the microscopic rate constant for Co-C bond cleavage. Eyring analysis reveals that the activation enthalpy for Co-C bond cleavage is 32 ± 1 kcal/mol, which is the same as for the cleavage reaction in solution. The origin of Co-C bond cleavage catalysis in the enzyme is, therefore, the large, favorable activation entropy of 61 ± 6 cal/(mol·K) (relative to 7 ± 1 cal/(mol·K) in solution). This represents a paradigm shift from traditional, enthalpy-based mechanisms that have been proposed for Co-C bond-breaking in B12 enzymes. The catalysis is proposed to arise from an increase in protein configurational entropy along the reaction coordinate.

Combined spectroscopic and computational analysis of the vibrational properties of vitamin B12 in its Co3+, Co2+, and Co1+ oxidation states

While the geometric and electronic structures of vitamin B12 (cyanocobalamin, CNCbl) and its reduced derivatives Co(2+)cobalamin (Co(2+)Cbl) and Co(1+)cobalamin (Co(1+)Cbl(-)) are now reasonably well established, their vibrational properties, in particular their resonance Raman (rR) spectra, have remained quite poorly understood. The goal of this study was to establish definitive assignments of the corrin-based vibrational modes that dominate the rR spectra of vitamin B12 in its Co(3+), Co(2+), and Co(1+) oxidation states. rR spectra were collected for all three species with laser excitation in resonance with the most intense corrin-based π → π* transitions. These experimental data were used to validate the computed vibrational frequencies, eigenvector compositions, and relative rR intensities of the normal modes of interest as obtained by density functional theory (DFT) calculations. Importantly, the computational methodology employed in this study successfully reproduces the experimental observation that the frequencies and rR excitation profiles of the corrin-based vibrational modes vary significantly as a function of the cobalt oxidation state. Our DFT results suggest that this variation reflects large differences in the degree of mixing between the occupied Co 3d orbitals and empty corrin π* orbitals in CNCbl, Co(2+)Cbl, and Co(1+)Cbl(-). As a result, vibrations mainly involving stretching of conjugated C-C and C-N bonds oriented along one axis of the corrin ring may, in fact, couple to a perpendicularly polarized electronic transition. This unusual coupling between electronic transitions and vibrational motions of corrinoids greatly complicates an assignment of the corrin-based normal modes of vibrations on the basis of their rR excitation profiles.

Charge Separation Propensity of the Coenzyme B12-Tyrosine Complex in Adenosylcobalamin-Dependent Methylmalonyl-CoA Mutase Enzyme

We report the electrophilic Fukui function analysis based on density functional reactivity theory (DFRT) to demonstrate the feasibility of the proton-coupled electron transfer (PCET) mechanism. To characterize the charge propensity of an electron-transfer site other than the proton-acceptor site of the coenzyme B12-tyrosine complex, several structural models (ranging from minimal to actual enzyme scaffolds) have been employed at DFT and QM/MM computations. It is shown, based on the methylmalonyl-CoA mutase (MCM) enzyme that substrate binding plays a significant role in displacing the phenoxyl proton of the tyrosine (Y89), which initiates the electron transfer from Y89 to coenzyme B12. PCET-based enzymatic reaction implies that one electron-reduced form of the AdoCbl cofactor induces the cleavage of the Co-C bond, as an alternative to its neutral analogue, which can assist in understanding the origin of the observed trillion-fold rate enhancement in MCM enzyme.

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.

Biochemistry of B12-cofactors in human metabolism

Vitamin B12, the "antipernicious anaemia factor", is a crystallisable cobalt-complex, which belongs to a group of unique "complete" corrinoids, named cobalamins (Cbl). In humans, instead of the "vitamin", two organometallic B12-forms are coenzymes in two metabolically important enzymes: Methyl-cobalamin, the cofactor of methionine synthase, and coenzyme B12 (adenosyl-cobalamin), the cofactor of methylmalonyl-CoA mutase. The cytoplasmatic methionine synthase catalyzes the transfer of a methyl group from N-methyl-tetrahydrofolate to homocysteine to yield methionine and to liberate tetrahydrofolate. In the mitochondrial methylmalonyl-CoA mutase a radical process transforms methylmalonyl-CoA (a remains e.g. from uneven numbered fatty acids) into succinyl-CoA, for further metabolic use. In addition, in the human mitochondria an adenosyl-transferase incorporates the organometallic group of coenzyme B12. In all these enzymes, the bound B12-derivatives engage (or are formed) in exceptional organometallic enzymatic reactions. This chapter recapitulates the physiological chemistry of vitamin B12, relevant in the context of the metabolic transformation of B12-derivatives into the relevant coenzyme forms and their use in B12-dependent enzymes.

The fragmentation-recombination mechanism of the enzyme glutamate mutase studied by QM/MM simulations

The radical mechanism of the conversion of glutamate to methylaspartate catalyzed by glutamate mutase is studied with quantum mechanical/molecular mechanical (QM/MM) simulations based on density functional theory (DFT/MM). The hydrogen transfer between the substrate and the cofactor is found to be rate limiting with a barrier of 101.1 kJ mol(-1). A careful comparison to the uncatalyzed reaction in water is performed. The protein influences the reaction predominantly electrostatically and to a lesser degree sterically. Our calculations shed light on the atomistic details of the reaction mechanism. The well-known arginine claw and Glu 171 ( Clostridium cochlearium notation) are found to have the strongest influence on the reaction. However, a catalytic role of Glu 214, Lys 322, Gln 147, Glu 330, Lys 326, and Met 294 is found as well. The arginine claw keeps the intermediates in place and is probably responsible for the enantioselectivity. Glu 171 temporarily accepts a proton from the glutamyl radical intermediate and donates it back at the end of the reaction. We relate our results to experimental data when available. Our simulations lead to further understanding of how glutamate mutase catalyzes the carbon skeleton rearrangement of glutamate.

On the importance of ribose orientation in the substrate activation of the coenzyme B12-dependent mutases

The degree to which the corrin ring portion of coenzyme B(12) can facilitate the H-atom-abstraction step in the glutamate mutase (GM)-catalyzed reaction of (S)-glutamate has been investigated with density functional theory. The crystal structure of GM identifies two possible orientations of the ribose portion of coenzyme B(12). In one orientation (A), the OH groups of the ribose extend away from the corrin ring, whereas in the other orientation (B) the OH groups, especially that involving O3', are instead directed towards the corrin ring. Our calculations identify a sizable stabilization amounting to about 30 kJ mol(-1) in the transition structure (TS) complex corresponding to orientation B (TS(B)CorIm). In the TS complex where the ribose instead is positioned in orientation A, no such effect is manifested. The observed stabilization in TS(B)CorIm appears to be the result of favorable interactions involving O3' and the corrin ring, including a C-HO hydrogen bond. We find that the degree of stabilization is not particularly sensitive to the Co-C distance. Our calculations show that any potential stabilization afforded to the H-atom-abstraction step by coenzyme B(12) is sensitive to the orientation of the ribose moiety.

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.

Cobalamin- and corrinoid-dependent enzymes

This chapter reviews the literature on cobalamin- and corrinoid-containing enzymes. These enzymes fall into two broad classes, those using methylcobalamin or related methylcorrinoids as prosthetic groups and catalyzing methyl transfer reactions, and those using adenosylcobalamin as the prosthetic group and catalyzing the generation of substrate radicals that in turn undergo rearrangements and/or eliminations.

Kinetic and spectroscopic studies of the ATP:corrinoid adenosyltransferase PduO from Lactobacillus reuteri: substrate specificity and insights into the mechanism of Co(II)corrinoid reduction

The PduO-type ATP:corrinoid adenosyltransferase from Lactobacillus reuteri ( LrPduO) catalyzes the formation of the essential Co-C bond of adenosylcobalamin (coenzyme B 12) by transferring the adenosyl group from cosubstrate ATP to a transient Co (1+)corrinoid species generated in the enzyme active site. While PduO-type enzymes have previously been believed to be capable of adenosylating only Co (1+)cobalamin (Co (1+)Cbl (-)), our kinetic data obtained in this study provide in vitro evidence that LrPduO can in fact also utilize the incomplete corrinoid Co (1+)cobinamide (Co (1+)Cbi) as an alternative substrate. To explore the mechanism by which LrPduO overcomes the thermodynamically challenging reduction of its Co (2+)corrinoid substrates, we have examined how the enzyme active site alters the geometric and electronic properties of Co (2+)Cbl and Co (2+)Cbi (+) by using electronic absorption, magnetic circular dichroism, and electron paramagnetic resonance spectroscopic techniques. Our data reveal that upon binding to LrPduO that was preincubated with ATP, both Co (2+)corrinoids undergo a partial ( approximately 40-50%) conversion to distinct paramagnetic Co (2+) species. The spectroscopic signatures of these species are consistent with essentially four-coordinate, square-planar Co (2+) complexes, based on a comparison with the results obtained in our previous studies of related enzymes. Consequently, it appears that the general strategy employed by adenosyltransferases for effecting Co (2+) --> Co (1+) reduction involves the formation of an "activated" Co (2+)corrinoid intermediate that lacks any significant axial bonding interactions, to stabilize the redox-active, Co 3d z (2) -based molecular orbital.

Coupling of hydrogenic tunneling to active-site motion in the hydrogen radical transfer catalyzed by a coenzyme B12-dependent mutase

Hydrogen transfer reactions catalyzed by coenzyme B(12)-dependent methylmalonyl-CoA mutase have very large kinetic isotope effects, indicating that they proceed by a highly quantal tunneling mechanism. We explain the kinetic isotope effect by using a combined quantum mechanical/molecular mechanical potential and semiclassical quantum dynamics calculations. Multidimensional tunneling increases the magnitude of the calculated intrinsic hydrogen kinetic isotope effect by a factor of 3.6 from 14 to 51, in excellent agreement with experimental results. These calculations confirm that tunneling contributions can be large enough to explain even a kinetic isotope effect >50, not because the barrier is unusually thin but because corner-cutting tunneling decreases the distance over which the system tunnels without a comparable increase in either the effective potential barrier or the effective mass for tunneling.

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.

Quantum catalysis in B12-dependent methylmalonyl-CoA mutase: experimental and computational insights

B12-dependent methylmalonyl-CoA mutase catalyses the interchange of a hydrogen atom and the carbonyl-CoA group on adjacent carbons of methylmalonyl-CoA to give the rearranged product, succinyl-CoA. The first step in this reaction involves the transient generation of cofactor radicals by homolytic rupture of the cobalt-carbon bond to generate the deoxyadenosyl radical and cob(II)alamin. This step exhibits a curious sensitivity to isotopic substitution in the substrate, methylmalonyl-CoA, which has been interpreted as evidence for kinetic coupling. The magnitude of the isotopic discrimination is large and a deuterium isotope effect ranging from 35.6 at 20 degrees C to 49.9 at 5 degrees C has been recorded. Arrhenius analysis of the temperature dependence of this isotope effect provides evidence for quantum tunnelling in this hydrogen transfer step. The mechanistic complexity of the observed rate constant for cobalt-carbon bond homolysis together with the spectroscopically silent nature of many of the component steps limits the insights that can be derived by experimental approaches alone. Computational studies using a newly developed geometry optimization scheme that allows determination of the transition state in the full quantum mechanical/molecular mechanical coordinate space have yielded novel insights into the strategy deployed for labilizing the cobalt-carbon bond and poising the resulting deoxyadenosyl radical for subsequent hydrogen atom abstraction.

Alternative pathways for radical dissipation in an active site mutant of B12-dependent methylmalonyl-CoA mutase

Methylmalonyl-CoA mutase catalyzes the adenosylcobalamin-dependent rearrangement of (2R)-methylmalonyl-CoA to succinyl-CoA. The crystal structure of the enzyme reveals that Y243 is in van der Waals contact with the methyl group of the substrate and suggests a possible role for it in the stereochemical control of the reaction. This hypothesis was tested by designing a molecular hole by replacing the phenolic side chain of Y243 with the methyl group of alanine. The Y243A mutation lowered the catalytic efficiency >(4 x 10(4))-fold compared to wild-type enzyme, the K(M)app for the cofactor approximately 4-fold, and the cob(II)alamin concentration under steady-state turnover conditions approximately 2-fold. However, the mutation did not appear to lead to loss of the stereochemical preference for the substrate. The Y243A mutation is expected to create a cavity and should, in principle, allow accommodation of bulkier substrates. To test this, we used ethylmalonyl-CoA and allylmalonyl-CoA as alternate substrates. Surprisingly, both analogues resulted in suicidal inactivation, albeit in an O(2)-dependent and O(2)-independent fashion, respectively. The inactivation by allylmalonyl-CoA was further investigated, and revealed formation of cob(II)alamin at an approximately 1.5-fold higher rate than with wild-type mutase under single-turnover conditions. Product analysis revealed a stoichiometric mixture of 5'-deoxyadenosine, aquocobalamin, and allylmalonyl-CoA. Taken together, these results are consistent with an internal electron transfer from cob(II)alamin to the substrate analogue radical. These studies serve to emphasize the fine control exerted by Y243 in the vicinity of the substrate to minimize radical extinction in side reactions.

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.

Co-C bond activation in methylmalonyl-CoA mutase by stabilization of the post-homolysis product Co2+ cobalamin

Despite decades of research, the mechanism by which coenzyme B12 (adenosylcobalamin, AdoCbl)-dependent enzymes promote homolytic cleavage of the cofactor's Co-C bond to initiate catalysis has continued to elude researchers. In this work, we utilized magnetic circular dichroism spectroscopy to explore how the electronic structure of the reduced B12 cofactor (i.e., the post-homolysis product Co2+ Cbl) is modulated by the enzyme methylmalonyl-CoA mutase. Our data reveal a fairly uniform stabilization of the Co 3d orbitals relative to the corrin pi/pi*-based molecular orbitals when Co2+ Cbl is bound to the enzyme active site, particularly in the presence of substrate. Contrastingly, our previous studies (Brooks, A. J.; Vlasie, M.; Banerjee, R.; Brunold, T. C. J. Am. Chem. Soc. 2004, 126, 8167-8180.) showed that when AdoCbl is bound to the MMCM active site, no enzymatic perturbation of the Co3+ Cbl electronic structure occurs, even in the presence of substrate (analogues). Collectively, these observations provide direct evidence that enzymatic Co-C bond activation involves stabilization of the post-homolysis product, Co2+ Cbl, rather than destabilization of the Co3+ Cbl "ground" state.

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.

Spectroscopic and Computational Studies of the ATP:Corrinoid Adenosyltransferase (CobA) fromSalmonella enterica: Insights into the Mechanism of Adenosylcobalamin Biosynthesis

CobA from Salmonella enterica is a member of an enzymatic system responsible for the de novo biosynthesis of adenosylcobalamin (AdoCbl), catalyzing the formation of the essential Co-C bond by transferring the adenosyl group from a molecule of ATP to a transient Co(1+)corrinoid species generated in the enzyme active site. A particularly fascinating aspect of this reaction is that the flavodoxin in vivo reducing agent that serves as the electron donor to CobA possesses a reduction potential that is considerably more positive than that of the Co(2+/1+) couple of the corrinoid substrate. To explore how CobA may overcome this challenge, we have employed electronic absorption, magnetic circular dichroism, and electron paramagnetic resonance (EPR) spectroscopies to probe the interaction between Co(3+)- and Co(2+)corrinoids and the enzyme active site. Our data reveal that while Co(3+)corrinoids interact only weakly with CobA, Co(2+)corrinoids undergo partial conversion to a new paramagnetic species that can be obtained in nearly quantitative yield when CobA is preincubated with the co-substrate ATP. This "activated" species is characterized by a distinct set of ligand field transitions in the near-IR spectral region and EPR parameters that are unprecedented for Co(2+)corrinoids. Analysis of these data on the basis of qualitative spectral correlations and density functional theory computations reveals that this unique Co(2+)corrinoid species possesses an essentially square-planar Co(2+) center that lacks any significant axial bonding interactions. Possible implications of these findings for the mechanism of Co(2+) --> Co(1+) reduction employed by CobA and Co-C bond-forming enzymes in general are explored.

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

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