Stabilization of the adenosyl radical in coenzyme B12 – a theoretical study

Computational investigations of B12-dependent enzymatic reactions

Nature employs two biologically active forms of vitamin B₁₂, adenosylcobalamin (or coenzyme B₁₂) and methylcobalamin, as cofactors in molecular transformations both in bacteria and mammals. Computational chemistry, guided by experimental data, has been used to explore fundamental characteristics of these enzymatic reactions. In particular, the quantum mechanics/molecular mechanics (QM/MM) method has proven to be a powerful tool in elucidating important characteristics of B₁₂-dependent enzymatic reactions. Herein, we will present a brief tutorial in conducting QM/MM calculations for B₁₂ enzymatic reactions. We will summarize recent contributions that target the use of QM/MM calculations in both photochemical and enzymatic reactions including AdoCbl-dependent ethanolamine ammonia lyase, glutamate mutase, and photoreceptor CarH.

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.

Binding interaction of a ring-hydroxylating dioxygenase with fluoranthene in Pseudomonas aeruginosa DN1

Pseudomonas aeruginosa DN1 can efficiently utilize fluoranthene as its sole carbon source, and the initial reaction in the biodegradation process is catalyzed by a ring-hydroxylating dioxygenase (RHD). To clarify the binding interaction of RHD with fluoranthene in the strain DN1, the genes encoding alpha subunit (RS30940) and beta subunit (RS05115) of RHD were functionally characterized through multi-technique combination such as gene knockout and homology modeling as well as molecular docking analysis. The results showed that the mutants lacking the characteristic alpha subunit and/or beta subunit failed to degrade fluoranthene effectively. Based on the translated protein sequence and Ramachandran plot, 96.5% of the primary amino-acid sequences of the alpha subunit in the modeled structure of the RHD were in the permitted region, 2.3% in the allowed region, but 1.2% in the disallowed area. The catalytic mechanism mediated by key residues was proposed by the simulations of molecular docking, wherein the active site of alpha subunit constituted a triangle structure of the mononuclear iron atom and the two oxygen atoms coupled with the predicted catalytic ternary of His₂₁₇-His₂₂₂-Asp₃₇₂ for the dihydroxylation reaction with fluoranthene. Those amino acid residues adjacent to fluoranthene were nonpolar groups, and the C₇-C₈ positions on the fluoranthene ring were estimated to be the best oxidation sites. The distance of C₇-O and C₈-O was 3.77 Å and 3.04 Å respectively, and both of them were parallel. The results of synchronous fluorescence and site-directed mutagenesis confirmed the roles of the predicted residues during catalysis. This binding interaction could enhance our understanding of the catalytic mechanism of RHDs and provide a solid foundation for further enzymatic modification.

Surprising Homolytic Gas Phase Co-C Bond Dissociation Energies of Organometallic Aryl-Cobinamides Reveal Notable Non-Bonded Intramolecular Interactions

Aryl-cobalamins are a new class of organometallic structural mimics of vitamin B12 designed as potential 'antivitamins B12 '. Here, the first cationic aryl-cobinamides are described, which were synthesized using the newly developed diaryl-iodonium method. The aryl-cobinamides were obtained as pairs of organometallic coordination isomers, the stereo-structure of which was unambiguously assigned based on homo- and heteronuclear NMR spectra. The availability of isomers with axial attachment of the aryl group, either at the 'beta' or at the 'alpha' face of the cobalt-center allowed for an unprecedented comparison of the organometallic reactivity of such pairs. The homolytic gas-phase bond dissociation energies (BDEs) of the coordination-isomeric phenyl- and 4-ethylphenyl-cobinamides were determined by ESI-MS threshold CID experiments, furnishing (Co-Csp 2 )-BDEs of 38.4 and 40.6 kcal mol-1 , respectively, for the two β-isomers, and the larger BDEs of 46.6 and 43.8 kcal mol-1 for the corresponding α-isomers. Surprisingly, the observed (Co-Csp 2 )-BDEs of the Coβ -aryl-cobinamides were smaller than the (Co-Csp 3 )-BDE of Coβ -methyl-cobinamide. DFT studies and the magnitudes of the experimental (Co-Csp 2 )-BDEs revealed relevant contributions of non-bonded interactions in aryl-cobinamides, notably steric strain between the aryl and the cobalt-corrin moieties and non-bonded interactions with and among the peripheral sidechains.

Three Dimensional Triply Resonant Sum Frequency Spectroscopy Revealing Vibronic Coupling in Cobalamins: Toward a Probe of Reaction Coordinates

Triply resonant sum frequency (TRSF) spectroscopy is a fully coherent mixed vibrational-electronic spectroscopic technique that is ideally suited for probing the vibrational-electronic couplings that become important in driving reactions. We have used cyanocobalamin (CNCbl) and deuterated aquacobalamin (D₂OCbl⁺) as model systems for demonstrating the feasibility of using the selectivity of coherent multidimensional spectroscopy to resolve electronic states within the broad absorption spectra of transition metal complexes and identify the nature of the vibrational and electronic state couplings. We resolve three short and long axis vibrational modes in the vibrationally congested 1400-1750 cm^-1 region that are individually coupled to different electronic states in the 18 000-21 000 cm^-1 region but have minimal coupling to each other. Double resonance with the individual vibrational fundamentals and their overtones selectively enhances the corresponding electronic resonances and resolves features within the broad absorption spectrum. This work demonstrates the feasibility of identifying coupling between different pairs of vibrational states with different electronic states that together form the reaction coordinate surface of transition metal enzymes.

Introducing Stable Radicals into Molecular Machines

Ever since their discovery, stable organic radicals have received considerable attention from chemists because of their unique optical, electronic, and magnetic properties. Currently, one of the most appealing challenges for the chemical community is to develop sophisticated artificial molecular machines that can do work by consuming external energy, after the manner of motor proteins. In this context, radical-pairing interactions are important in addressing the challenge: they not only provide supramolecular assistance in the synthesis of molecular machines but also open the door to developing multifunctional systems relying on the various properties of the radical species. In this Outlook, by taking the radical cationic state of 1,1'-dialkyl-4,4'-bipyridinium (BIPY^•+) as an example, we highlight our research on the art and science of introducing radical-pairing interactions into functional systems, from prototypical molecular switches to complex molecular machines, followed by a discussion of the (i) limitations of the current systems and (ii) future research directions for designing BIPY^•+-based molecular machines with useful functions.

Cloning, Expression, Mutagenesis Library Construction of Glycerol Dehydratase, and Binding Mode Simulation of Its Reactivase with Ligands

The production of 1, 3-propanediol (1, 3-PD) and 3-hydroxypropionaldehyde (3-HPA) by enzyme reaction has been a hot field, and glycerol dehydratase (GDHt) is the key and rate-limiting enzyme involved in their biosynthesis. The gldABC gene encoding GDHt was cloned from Klebsiella pneumoniae, and the activity of the corresponding proteins expressed extracellularly and intracellularly was 6.8 and 3.2 U/mg, respectively, about six and three times higher than that of the wild strain. The change of amino acids for the β subunit can adjust the length of the Co-N bond and affect the homolysis rate of the Co-C bond to change GDHt activity. The expression plasmid, pET-32a-gldAC (containing no gldB which encodes the β subunit of GDHt), was constructed to build the mutagenesis library to improve the GDHt activity. The binding models of glycerol dehydratase reactivation factor (GDHtR) with ATP, CTP, or GTP were simulated by semi-flexible docking, respectively, and there was almost no difference between them. This research provided the basis for studying the quantitative structure-activity relationships between GDHtR and its ligands, as well as searching inexpensive ligands to replace ATP. These results and methods are of great use in economical and highly efficient production of 3-HPA and 1, 3-PD by the enzyme method.

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.

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.

Investigating the intermediates in the reaction of ribonucleoside triphosphate reductase from Lactobacillus leichmannii: An application of HF EPR-RFQ technology

In this investigation high-frequency electron paramagnetic resonance spectroscopy (HFEPR) in conjunction with innovative rapid freeze-quench (RFQ) technology is employed to study the exchange-coupled thiyl radical-cob(II)alamin system in ribonucleotide reductase from a prokaryote Lactobacillus leichmannii. The size of the exchange coupling (Jex) and the values of the thiyl radical g tensor are refined, while confirming the previously determined (Gerfen et al. (1996) [20]) distance between the paramagnets. Conclusions relevant to ribonucleotide reductase catalysis and the architecture of the active site are presented. A key part of this work has been the development of a unique RFQ apparatus for the preparation of millisecond quench time RFQ samples which can be packed into small (0.5 mm ID) sample tubes used for CW and pulsed HFEPR--lack of this ability has heretofore precluded such studies. The technology is compatible with a broad range of spectroscopic techniques and can be readily adopted by other laboratories.

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.

Structural basis for adenosylcobalamin activation in AdoCbl-dependent ribonucleotide reductases

Class II ribonucleotide reductases (RNR) catalyze the formation of an essential thiyl radical by homolytic cleavage of the Co-C bond in their adenosylcobalamin (AdoCbl) cofactor. Several mechanisms for the dramatic acceleration of Co-C bond cleavage in AdoCbl-dependent enzymes have been advanced, but no consensus yet exists. We present the structure of the class II RNR from Thermotoga maritima in three complexes: (i) with allosteric effector dTTP, substrate GDP, and AdoCbl; (ii) with dTTP and AdoCbl; (iii) with dTTP, GDP, and adenosine. Comparison of these structures gives the deepest structural insights so far into the mechanism of radical generation and transfer for AdoCbl-dependent RNR. AdoCbl binds to the active site pocket, shielding the substrate, transient 5'-deoxyadenosyl radical and nascent thiyl radical from solution. The e-propionamide side chain of AdoCbl forms hydrogen bonds directly to the α-phosphate group of the substrate. This interaction appears to cause a "locking-in" of the cofactor, and it is the first observation of a direct cofactor-substrate interaction in an AdoCbl-dependent enzyme. The structures support an ordered sequential reaction mechanism with release or relaxation of AdoCbl on each catalytic cycle. A conformational change of the AdoCbl adenosyl ribose is required to allow hydrogen transfer to the catalytic thiol group. Previously proposed mechanisms for radical transfer in B12-dependent enzymes cannot fully explain the transfer in class II RNR, suggesting that it may form a separate class that differs from the well-characterized eliminases and mutases.

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.

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.

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.

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.

A new paradigm for electrostatic catalysis of radical reactions in vitamin B12 enzymes

The catalytic power of enzymes containing coenzyme B(12) cofactor has been, in some respects, the "last bastion" for the strain hypothesis. The present work explores the origin of this effect by using simulation methods that overcome the sampling difficulties of previous energy minimization studies. It is found that the major part of the catalytic effect is due to the electrostatic interaction between the ribose and the protein, and that the strain contribution is very small. Remarkably, enzymes can use electrostatic effects even in a radical process, when the charge distribution of the reacting fragments does not change significantly during the reaction. Electrostatic catalysis can, in such cases, be obtained by attaching a polar group to the leaving fragment and designing an active site that interacts more strongly with this group in the product state than in the reactant state. The finding that evolution had to use this trick provides further evidence to the observation that it is extremely hard to catalyze enzymatic reactions by nonelectrostatic factors. The trick used by B(12) enzymes may, in fact, be a very powerful new strategy in enzyme design.

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.

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.

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.