The enzymatic activation of coenzyme B12

The enzymatic "activation" of coenzyme B12 (5'-deoxyadenosylcobalamin, AdoCbl), in which homolysis of the carbon-cobalt bond of the coenzyme is catalyzed by some 10(9)- to 10(14)-fold, remains one of the outstanding problems in bioinorganic chemistry. Mechanisms which feature the enzymatic manipulation of the axial Co-N bond length have been investigated by theoretical and experimental methods. Classical mechanochemical triggering, in which steric compression of the long axial Co-N bond leads to increased upward folding of the corrin ring and stretching of the Co-C bond is found to be feasible by molecular modeling, but the strain induced in the Co-C bond seems to be too small to account for the observed catalytic power. The modeling study shows that the effect is a steric one which depends on the size of the axial nucleotide base, as substitution of imidazole (Im) for the normal 5,6-dimethylbenzimidazole (Bzm) axial base decreases the Co-C bond labilization considerably. An experimental test was thus devised using the coenzyme analog with Im in place of Bzm (Ado(Im)Cbl). Studies of the enzymatic activation of this analog by the B12-dependent ribonucleoside triphosphate reductase from Lactobacillus leichmannii coupled with studies of the non-enzymatic homolytic lability of the Co-C bond of Ado(Im)Cbl show that the enzyme is only slightly less efficient (3.8-fold, 0.8 kcal mol(-1)) at activating Ado(Im)Cbl than at activating AdoCbl itself. This suggests, in agreement with the modeling study, that mechanochemical triggering can make only a small contribution to the enzymatic activation of AdoCbl. Another possibility, electronic stabilization of the Co(II) homolysis product by compression of the axial Co-N bond, requires that enzymatic activation be sensitive to the basicity of the axial nucleotide. Preliminary studies of the enzymatic activation of a coenzyme analog with a 5-fluoroimidazole axial nucleotide suggest that the catalysis of Co-C bond homolysis may indeed be significantly slowed by the decrease in basicity.     

Accurate redetermination of the X-ray structure and electronic bonding in adenosylcobalamin

The electronic structure of adenosylcobalamin (B12 coenzyme, AdoCbl) has been calculated by a density functional method, using the orthogonalized linear combination of the atomic orbital method (OLCAO). Since a fixed accurately determined geometry was needed in such calculations, the crystal structure of adenosylcobalamin has been redone and refined to R = 0.065, using synchrotron diffraction data. Comparison with the recently reported electronic structures of cyano- (CNCbl) and methylcobalamin (MeCbl) shows that the net charges and bond orders vary only on the axial donors. The values in the three cobalamins suggest that the Co-C bond in MeCbl has a strength similar to that in AdoCbl, but it is significantly weaker that that in CNCbl. Present results are compared with those previously reported for the analogous corrin derivatives; i.e., simplified cobalamins with the side chains a-f replaced by H atoms. Despite a qualitative agreement, a discrepancy in the calculated HOMO-LUMO gap is found.     

DFT study of Co-C bond cleavage in the neutral and one-electron-reduced alkyl-cobalt(III) phthalocyanines

Density functional theory (DFT) has been applied to the analysis of the structural and electronic properties of the alkyl-cobalt(III) phthalocyanine complexes, [CoIIIPc]-R (Pc = phthalocyanine, R = Me or Et), and their pyridine adducts. The BP86/6-31G(d) level of theory shows good reliability for the optimized axial bond lengths and bond dissociation energies (BDEs). The mechanism of the reductive cleavage was probed for the [CoIIIPc]-Me complex which is known as a highly effective methyl group donor. In the present analysis, which follows a recent study on the reductive Co-C bond cleavage in methylcobalamin (J. Phys. Chem. B 2007, 111, 7638-7645), it is demonstrated that addition of an electron and formation of the pi-anion radical [CoIII(Pc*)]-Me- significantly lowers the energetic barrier required for homolytic Co-C bond dissociation. Such BDE lowering in [CoIII(Pc*)]-Me- arises from the involvement of two electronic states: upon electron addition, a quasi-degenerate pi*Pc state is initially formed, but when the cobalt-carbon bond is stretched, the unpaired electron moves to a sigma*Co-C state and the final cleavage involves the three-electron (sigma)2(sigma*)1 bond. As in corrin complexes, the pi*Pc-sigma*Co-C states crossing does not take place at the equilibrium geometry of [CoIII(Pc*)]-Me- but only when the Co-C bond is stretched to approximately 2.3 A. The DFT computed Co-C BDE of 23.3 kcal/mol in the one-electron-reduced phthalocyanine species, [CoIII(Pc*)]-Me-, is lowered by approximately 37% compared to the neutral Py-[CoIIIPc]-Me complex where BDE = 36.8 kcal/mol. A similar comparison for the corrin-containing complexes shows that a DFT computed BDE of 20.4 kcal/mol for [CoIII(corrin*)]-Me leads to approximately 45% bond strength reduction, in comparison to 37.0 kcal/mol for Im-[CoIII(corrin)]-Me+. These results suggest some preference by the alkylcorrinoids for the reductive cleavage mechanism.     

Solution structure and thermolysis of Cobeta-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).     

光声量热法测定辅酶B12光解的焓变和体积变化

用时间分辨的光声量热法研究了浦防B12的光解反应，首次测定了脉冲激光诱导输酶B12光解反应的治变和反应体积变化，其值分别为131±12kJ·mol-1和6±1mL·mol-1，推测该反应体积变化可能与Co－C健断裂引起的柔性咕啉环构象变化有关．     

Co-C bond dissociation energy and reaction volume change of 2',5'-dideoxyadenosylcobalamin studied by laser-induced time-resolved photoacoustic calorimetry

Time resolved photoacoustic calorimetry (PAC) was applied to a study of the photolysis of a coenzyme B(12) analog 2',5'-dideoxyadenosylcobalamin, which lacks an -OH group at the 2' position of ribofuranose ring. In aqueous solution, we report for the first time the quantum yield Phi(d) (0.25+/-0.02), Co-C bond dissociation energy (BDE; 31.8+/-2.5 kcal mol(-1)) and reaction volume change deltaV(R) (6.5+/-0.5 ml mol(-1)) due to conformation changes of the corrin ring and its side chains accompanying the cleavage of the Co-C bond. These values for the analog are very similar to those for the natural cofactor. Based our results and previous studies, a possible explanation for the similarity in their structure and properties versus the large difference in their enzymatic activity is discussed.     

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.     

The solution structure of adenosylcobalamin and adenosylcobinamide determined by nOe-restrained molecular dynamics simulations

The solution structure of coenzyme B₁₂ (5′-deoxyadenosylcobalamin, AdoCbl) and the corresponding cobinamide, AdoCbi⁺, in which the axial nucleotide has been chemically removed, have been investigated using NMR-restrained molecular dynamics (MD) and simulated annealing (SA) calculations. The nOe cross peaks in the ROESY spectrum of both AdoCbl and AdoCbi⁺ are consistent with the presence of at least two principal conformations for each compound in solution. In the first, termed the southern conformation, the adenosyl (Ado) ligand is over the C ring of the molecule, the structure observed in the solid state. In the second, the Ado ligand has undergone an anticlockwise rotation and is over C10 in the eastern quadrant of the molecule. A two-state MD/SA simulation was used omitting nOes that arise only from the eastern conformation and that arise only from the southern conformation, respectively. Consensus structures were obtained by averaging the coordinates of 25 annealed structures of the southern and eastern conformations, respectively, of AdoCbl and AdoCbi⁺, followed by energy minimization. The consensus structure of the southern conformation of AdoCbl agrees well with the solid-state structure and has a very similar corrin fold angle. Several observations show that AdoCbl is considerably more rigid than AdoCbi⁺, and indeed is one of the most rigid cobalt corrinoids studied by these methods to date: the variability in the conformations of the corrin ring between the family of 25 annealed structure and the consensus structure is much smaller for AdoCbl than for AdoCbi⁺; during MD simulations, the previously demonstrated flexibility of the corrin ring as gauged by the corrin ruf angle (C5–Co–C15) is preserved for AdoCbi⁺ but is considerably diminished in AdoCbl because of a decrease in the maximum fold angle and an increase in the minimum fold angle achieved in the latter; the range of values of the Co–C bond length experienced in AdoCbi⁺ is substantially larger than in AdoCbl; the Ado ligand visits many more orientations relative to the corrin ring in AdoCbi⁺ than in AdoCbl; the pyrrole rings in AdoCbl undergo smaller deformations than in AdoCbi⁺; and the “breathing motion” of the corrin ring in which C5, C10 and C15 oscillate from above to below the mean corrin plane is much less pronounced in AdoCbl than in AdoCbi⁺. This rigidity is attributed to the presence of two bulky ligands in AdoCbl, the Ado ligand in the upper (β) axial position and the 5,6-dimethylbenzimidazole (bzm) ligand in the lower () axial ligand position, in contrast to the other structures which have only one or other of these two bulky ligands. The corrin fold angle in AdoCbl is significantly larger than that in AdoCbi⁺, a finding that is in agreement with a previous observation that CH₃Cbl has a larger fold angle than CH₃Cbi⁺; this implies that base-on corrins are under steric strain.     

Synthesis, spectroscopic characterization, axial base coordination equilibrium and photolytic kinetics studies of a new coenzyme B12 analogue-3'-deoxy-2',3'-anhydrothymidylcobalamin

A new coenzyme B12 (AdoCbl) analogue, 3'-deoxy-2',3'-didehydrothymidylcobalamin (2',3'-anThyCbl) was prepared by the reaction of 5'-iodo-3'-deoxy-2',3'-dihydrothmidine with reduced B12a, and characterized by UV-Vis, CD, ESI-MS and NMR spectroscopies. Its axial base (dbzm) coordination equilibria with pH's and temperatures were investigated and showed similar features to those of coenzyme B12. Photolytic dynamics studies under homolytic and heterolytic conditions demonstrated that the Co-C bond of the analogue is slightly more photolabile relative to coenzyme B12.     

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.     

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

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

Influence of environment on the electronic structure of Cob(III)alamins: time-resolved absorption studies of the S(1) state spectrum and dynamics

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

A localized orbital analysis of the thermochemical errors in hybrid density functional theory: achieving chemical accuracy via a simple empirical correction scheme

This paper describes an empirical localized orbital correction model which improves the accuracy of density functional theory (DFT) methods for the prediction of thermochemical properties for molecules of first and second row elements. The B3LYP localized orbital correction version of the model improves B3LYP DFT atomization energy calculations on the G3 data set of 222 molecules from a mean absolute deviation (MAD) from experiment of 4.8 to 0.8 kcal/mol. The almost complete elimination of large outliers and the substantial reduction in MAD yield overall results comparable to the G3 wave-function-based method; furthermore, the new model has zero additional computational cost beyond standard DFT calculations. The following four classes of correction parameters are applied to a molecule based on standard valence bond assignments: corrections to atoms, corrections to individual bonds, corrections for neighboring bonds of a given bond, and radical environmental corrections. Although the model is heuristic and is based on a 22 parameter multiple linear regression to experimental errors, each of the parameters is justified on physical grounds, and each provides insight into the fundamental limitations of DFT, most importantly the failure of current DFT methods to accurately account for nondynamical electron correlation.     

Structure and vibrational spectra of the enol form of hexafluoro-acetylacetone. A density functional theoretical study

Molecular and vibrational structure of 1,1,1,6,6,6-hexafluoropentane-2,4-dione (hexafluoro-acetylacetone) have been investigated by means of density functional theory (DFT) calculations and have been compared with those of acetylacetone, the parent molecule. According to the theoretical calculations HFAA has an asymmetric structure with hydrogen bond strength of about 12 kcal mol(-1), about 6 kcal mol(-1) less than that of acetylacetone. This weakening of hydrogen bond is consistent with frequency shifts for OH/OD stretching, OH/OD out of plane bending and O...O stretching modes upon substitution of methyl hydrogen atoms with fluorine atoms. The symmetric structure based on electron diffraction data is interpreted as superposition of two asymmetric structures.     

ER-BOC calculations of crystal bulk properties from smallest-cluster models

Ab initio calculation of bulk properties of crystals with a high accuracy, which is a long-time goal of solid chemistry and physics, is still difficult and expensive because a large cluster is required as a crystal structure model. This article proposes a model based on density functional theory (DFT) quantum chemistry calculations and the assumption that the bond order of a given atom with its nearest atoms in a compound is conserved over the entire range from its diatomic molecules to clusters and further to crystals. This entire range bond order conservation (ER-BOC) provides an effective way to correlate bulk properties of crystals with those of the corresponding molecules and small clusters. By combining this ER-BOC principle with hybrid DFT quantum chemistry calculations, accurate predictions of the bulk bond lengths of a crystal can be made using calculations on small clusters.     

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.     

Interaction of metal porphyrins with fullerene C60: a new insight

The electronic structure and bonding in the noncovalent, supramolecular complexes of fullerene C60 with a series of first-row transition metal porphines MP (M=Fe, Co, Ni, Cu, Zn) have been re-examined with DFT methods. A dispersion correction was made for the C60-MP binding energy through an empirical method (J. Comput. Chem. 2004, 25, 1463). Several density functionals and two types of basis sets were employed in the calculations. Our calculated results are rather different from those obtained in a recent paper (J. Phys. Chem. A 2005, 109, 3704). The ground state of C60.FeP is predicted to be high spin (S=2); the low-spin (S=0), closed-shell state is even higher in energy than the intermediate-spin (S=1) state. With only one electron in the Co-dz2 orbital, the calculated Co-C60 distance is in fact rather short, about 0.1 A longer than the Fe-C60 distance in high-spin C60.FeP. Double occupation of an M-dz2 orbital in MP prevents close association of any axial ligand, and so the Ni-C60, Cu-C60, and Zn-C60 distances are much longer than the Co-C60 one. The evaluated MP-C60 binding energies (Ebind) are 0.8 eV (18.5 kcal/mol) for M=Fe/Co and 0.5 eV (11.5 kcal/mol) for M=Ni/Cu/Zn (Ebind is about 0.2 eV larger in the case of C60-MTPP). They are believed to be reliable and accurate based on our dispersion-corrected DFT calculations that included the counterpoise (CP) correction. The effects of the C60 contact on the redox properties of MP were also examined.     

Conformational analysis: crystallographic,NMR, and molecular mechanics studies of flexible sulfonic esters

Two novel X-ray structures of the sulfonic ester derivatives 2-(6-iodo-1,3-benzodioxol-5-yl)ethyl 4-nitrobenzenesulfonate, 3, and 2-(6-iodo-1,3-benzodioxol-5-yl)ethyl 4-methylbenzenesulfonate, 4, have been obtained in a study aimed at analyzing the structures and conformations of sulfonic ester derivatives that are routinely used in alkaloid syntheses. The crystal structure of 4 is highly unusual, containing four independent molecules that belong to two distinct conformational types: (1) a hairpin conformation (stabilized mainly by intramolecular pi-stacking) and (2) a stepped conformation (stabilized mainly by intermolecular pi-stacking). Compound 3, on the other hand, crystallizes exclusively as the hairpin conformer. New MM+ force field parameters for sulfonic esters have been developed using the X-ray data, empirical rules, and DFT calculations to estimate the bond dipole parameters. Grid searches of conformational space for 3 and 4 using MM methods show that there are several gas-phase conformations within 5 kcal/mol of the global minimum and that the lowest energy conformations (by approximately 4.6 kcal/mol) are of the hairpin type. Analysis of the MM conformational energies suggests that the dominant intramolecular interaction stabilizing the hairpin conformations of 3 and 4 is van der Waals attraction. Moreover, the lattice energies for packing the hairpin conformations of 3 and 4 are approximately 4 kcal/mol more favorable than for the stepped conformations. Various intermolecular interactions contribute to the complexity of the observed crystal structures of 3 and 4, including electrostatic attraction between O and I atoms in neighboring molecules. Langevin dynamics (LD) simulations at several temperatures (6.0 ns, friction coefficient = 2.5 ps(-1)) indicate that the conformational exchange rates are approximately 10(10)-10(11) s(-1) over the temperature range 213-400 K, accounting for the temperature-independent (1)H NMR spectra of 3 and 4.