Co-C dissociation of adenosylcobalamin (coenzyme B12): role of dispersion, induction effects, solvent polarity, and relativistic and thermal corrections

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.

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.

Electronic and photolytic properties of hydridocobalamin

Hydridocobalamin (HCbl), is a known member of the B₁₂ family of molecules (cobalamins, Cbls) yet unlike other well-studied Cbls, little is known of the electronic and photolytic properties of this species. Interest in HCbl has increased significantly in recent years when at least three experimentally proposed mechanisms implicate HCbl as an intermediary in the photoreaction of coenzyme B₁₂-dependent photoreceptor CarH. Specifically, cleavage of the Co-C_5' bond of coenzyme B₁₂ could lead to a β-hydride or β‑hydrogen elimination reaction to form HCbl. HCbl is known to be a transient species where the oxidation state of the Co is variable; Co(I)-H⁺ ↔ Co(II)-H ↔ Co(III)-H^-. Further, HCbl is a very unstable with a pKa of ~1. This complicates experimental studies and to the best of our knowledge there are no available crystal structures of HCbl - either for the isolated molecule or bound to an enzyme. In this study, the electronic structure, photolytic properties, and reactivity of HCbl were explored to determine the preferred oxidation state as well as its potential role in the formation of the photoproduct in CarH. Natural bond orbital (NBO) analysis was performed to determine the oxidation state of Co in isolated HCbl. Based on the NBO analysis of HCbl, Co clearly had excess negative charge, which is in stark contrast to other alkylCbls where the Co ion is marked by significant positive charge. In sum, NBO results indicate that the CoH bond is strongly polarized and almost ionic. It can be described as protonated Co(I). In addition, DFT was used to explore the bond dissociation energy of HCbl based on homolytic cleavage of the CoH bond. TD-DFT calculations were used to compare computed electronic transitions to the experimentally determined absorption spectrum. The photoreaction of CarH was explored using an isolated model system and a pathway for hydrogen transfer was found. Finally, quantum mechanics/molecular mechanics (QM/MM) calculations were employed to investigate the formation of HCbl in CarH.

Cobalt Catalyst Determines Regioselectivity in Ring Opening of Epoxides with Aryl Halides

Ring-opening of epoxides furnishing either linear or branched products belongs to the group of classic transformations in organic synthesis. However, the regioselective cross-electrophile coupling of aryl epoxides with aryl halides still represents a key challenge. Herein, we report that the vitamin B₁₂/Ni dual-catalytic system allows for the selective synthesis of linear products under blue-light irradiation, thus complementing methodologies that give access to branched alcohols. Experimental and theoretical studies corroborate the proposed mechanism involving alkylcobalamin as an intermediate in this reaction.

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.

Light Mediated Properties of a Thiolato-Derivative of Vitamin B12

Vitamin B₁₂ derivatives (Cbls = cobalamins) exhibit photolytic properties upon excitation with light. These properties can be modulated by several factors including the nature of the axial ligands. Upon excitation, homolytic cleavage of the organometallic bond to the upper axial ligand takes place in photolabile Cbls. The photosensitive nature of Cbls has made them potential candidates for light-activated drug delivery. The addition of a fluorophore to the nucleotide loop of thiolato Cbls has been shown to shift the region of photohomolysis to within the optical window of tissue (600-900 nm). With this possibility, there is a need to analyze photolytic properties of unique Cbls which contain a Co-S bond. Herein, the photodissociation of one such Cbl, namely, N-acetylcysteinylcobalamin (NACCbl), is analyzed based on density functional theory (DFT) and time-dependent DFT (TD-DFT) calculations. The S₀ and S₁ potential energy surfaces (PESs), as a function of axial bond lengths, were computed to determine the mechanism of photodissociation. Like other Cbls, the S₁ PES contains metal-to-ligand charge transfer (MLCT) and ligand field (LF) regions, but there are some unique differences. Interestingly, the S₁ PES of NACCbl contains three distinct minima regions opening several possibilities for the mechanism of radical pair (RP) formation. The mild photoresponsiveness, observed experimentally, can be attributed to the small gap in energy between the S₁ and S₀ PESs. Compared to other Cbls, the gap shown for NACCbl is neither exactly in line with the alkyl Cbls nor the nonalkyl Cbls.

Modeling the Reactions Catalyzed by Coenzyme B12 Dependent Enzymes: Accuracy and Cost-Quality Balance

The reactions catalyzed by coenzyme B₁₂ dependent enzymes are formally initiated by the homolytic cleavage of a carbon-cobalt bond and a subsequent or concerted H-atom-transfer reaction. A reasonable model chemistry for describing those reactions should, therefore, account for an accurate description of both reactions. The inherent limitation due to the necessary system size renders the coenzyme B₁₂ system a suitable candidate for DFT or hybrid QM/MM methods; however, the accurate description of both homolytic Co-C cleavage and H-atom-transfer reactions within this framework is challenging and can lead to controversial results with varying accuracy. We present an assessment study of 16 common density functionals applied to prototypical model systems for both reactions. H-abstraction reactions were modeled on the basis of four reference reactions designed to resemble a broad range of coenzyme B₁₂ reactions. The Co-C cleavage reaction is treated by an ONIOM(QM/MM) setup that is in excellent agreement with solution-phase experimental data and is as accurate as full DFT calculations on the complete model system. We find that the meta-GGAs TPSS-D3 and M06L-D3 and the meta-hybrid M06-D3 give the best overall performance with MUEs for both types of reactions below 10 kJ mol^-1. Our recommended model chemistry allows for a fast and accurate description of coenzyme B₁₂ chemistry that is readily applicable to study the reactions in an enzymatic framework.

Development of CHARMM-Compatible Force-Field Parameters for Cobalamin and Related Cofactors from Quantum Mechanical Calculations

Corrinoid cofactors such as cobalamin are used by many enzymes and are essential for most living organisms. Therefore, there is broad interest in investigating cobalamin-protein interactions with molecular dynamics simulations. Previously developed parameters for cobalamins are based mainly on crystal structure data. Here, we report CHARMM-compatible force field parameters for several corrinoids developed from quantum mechanical calculations. We provide parameters for corrinoids in three oxidation states, Co³⁺, Co²⁺, and Co¹⁺, and with various axial ligands. Lennard-Jones parameters for the cobalt center in the Co(II) and Co(I) states were optimized using a helium atom probe, and partial atomic charges were obtained with a combination of natural population analysis (NPA) and restrained electrostatic potential (RESP) fitting approaches. The Force Field Toolkit was used to optimize all bonded terms. The resulting parameters, determined solely from calculations of cobalamin alone or in water, were then validated by assessing their agreement with density functional theory geometries and by analyzing molecular dynamics simulation trajectories of several corrinoid proteins for which X-ray crystal structures are available. In each case, we obtained excellent agreement with the reference data. In comparison to previous CHARMM-compatible parameters for cobalamin, we observe a better agreement for the fold angle and lower RMSD in the cobalamin binding site. The approach described here is readily adaptable for developing CHARMM-compatible force-field parameters for other corrinoids or large biomolecules.

Energy vs. density on paths toward more exact density functionals

The density-energy balance of density functionals is quantified using trial densities. Density errors commonly amount to only a few kJ mol⁻¹.

Toward Quantitatively Accurate Calculation of the Redox-Associated Acid-Base and Ligand Binding Equilibria of Aquacobalamin

Redox processes in complex transition metal-containing species are often intimately associated with changes in ligand protonation states and metal coordination number. A major challenge is therefore to develop consistent computational approaches for computing pH-dependent redox and ligand dissociation properties of organometallic species. Reduction of the Co center in the vitamin B12 derivative aquacobalamin can be accompanied by ligand dissociation, protonation, or both, making these properties difficult to compute accurately. We examine this challenge here by using density functional theory and continuum solvation to compute Co-ligand binding equilibrium constants (Kon/off), pKas, and reduction potentials for models of aquacobalamin in aqueous solution. We consider two models for cobalamin ligand coordination: the first follows the hexa, penta, tetra coordination scheme for Co(III), Co(II), and Co(I) species, respectively, and the second model features saturation of each vacant axial coordination site on Co(II) and Co(I) species with a single, explicit water molecule to maintain six directly interacting ligands or water molecules in each oxidation state. Comparing these two coordination schemes in combination with five dispersion-corrected density functionals, we find that the accuracy of the computed properties is largely independent of the scheme used, but including only a continuum representation of the solvent yields marginally better results than saturating the first solvation shell around Co throughout. PBE performs best, displaying balanced accuracy and superior performance overall, with RMS errors of 80 mV for seven reduction potentials, 2.0 log units for five pKas and 2.3 log units for two log Kon/off values for the aquacobalamin system. Furthermore, we find that the BP86 functional commonly used in corrinoid studies suffers from erratic behavior and inaccurate descriptions of Co-axial ligand binding, leading to substantial errors in predicted pKas and Kon/off values. These findings demonstrate the effectiveness of the present approach for computing electrochemical and thermodynamic properties of a complex transition metal-containing cofactor.

Mechanism of Co-C photodissociation in adenosylcobalamin

A mechanism of Co-C bond photodissociation in the base-on form of adenosylcobalamin (AdoCbl) was investigated by time-dependent density functional theory (TD-DFT). The key mechanistic step involves singlet radical pair (RP) generation from the first electronically excited state (S1). To connect TD-DFT calculations with ultra-fast excited state dynamics, the potential energy surface (PES) of the S1 state was constructed using Co-C and Co-NIm axial coordinates. The S1 PES can be characterized by two minima separated by a seam resulting from the crossing of two surfaces, of metal-to-ligand charge transfer (MLCT) character near the minimum, and a shallow ligand field (LF) surface at elongated axial bond distances. Only one possible pathway for photolysis (path A) was identified based on energetic grounds. This pathway is characterized by the first elongation of the Co-C bond, followed by photolysis from an LF state where the axial base is partially detached. A new perspective on the photolysis of AdoCbl is then gained by connecting TD-DFT results with available experimental observations.

Theoretical Study of Spin Crossover in 30 Iron Complexes

Iron complexes are important spin crossover (SCO) systems with vital roles in oxidative metabolism and promising technological potential. The SCO tendency depends on the free energy balance of high- and low-spin states, which again depends on physical effects such as dispersion, relativistic effects, and vibrational entropy. This work studied 30 different iron SCO systems with experimentally known thermochemical data, using 12 different density functionals. Remarkably general entropy-enthalpy compensation across SCO systems was identified (R = 0.82, p = 0.002) that should be considered in rational SCO design. Iron(II) complexes displayed higher ΔH and ΔS values than iron(III) complexes and also less steep compensation effects. First-coordination sphere ΔS values computed from numerical frequencies reproduce most of the experimental entropy and should thus be included when modeling spin-state changes in inorganic chemistry (R = 0.52, p = 3.4 × 10(-3); standard error in TΔS ≈ 4.4 kJ/mol at 298 K vs 16 kJ/mol of total TΔS on average). Zero-point energies favored high-spin states by 9 kJ/mol on average. Interestingly, dispersion effects are surprisingly large for the SCO process (average: 9 kJ/mol, but up to 33 kJ/mol) and favor the more compact low-spin state. Relativistic effects favor low-spin by ∼9 kJ/mol on average, but up to 24 kJ/mol. B3LYP*, TPSSh, B2PLYP, and PW6B95 performed best for the typical calculation scheme that includes ZPE. However, if relativistic and dispersion effects are included, only B3LYP* remained accurate. On average, high-spin was favored by LYP by 11-15 kJ/mol relative to other correlation functionals, and by 4.2 kJ/mol per 1% HF exchange in hybrids. 13% HF exchange was optimal without dispersion, and 15% was optimal with all effects included for these systems.

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.