Axial bonding in alkylcobalamins: DFT analysis of the inverse versus normal trans influence

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 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.

Trans Influence and Substituent Effects on the HOMO-LUMO Energy Gap and Stokes Shift in Ru Mono-Diimine Derivatives

The ground (S0) and excited triplet (T1) electronic states and corresponding optical spectra of a series of cationic complexes [RuH(CO)L(PPh3)2]+ (L=2,2´-bipyridyl) (Rubpy), 4,4´-dicarboxylic-2,2´-bipyridyl (Rudcbpy), bis-4,4'-(N-methylamide)-2,2´-bipyridyl (Rudamidebpy), bis-4,4'-(methyl)-2,2´-bipyridyl (RudMebpy), [Ru(CO)2dcbpy(PPh3)2]2+ (Ru(2CO)dcbpy), and [Ru(H)2dcbpy(PPh3)2] (Ru(2H)dcbpy) have been studied by combined Density Functional/Time-Dependent Density Functional (DFT/TDDFT) techniques using different combinations of DFT exchange-correlation functionals and basis sets. PBE0/LANL2DZ provided more accurate geometries to describe S0 whereas B3LYP/LANL2DZ predicted spectral energies that correlated better with the available experiment data. The Ru (II) complexes with different substituents emit photons ranging from 560-610 nm in the series RudMebpy, Rubpy, Rudamidebpy, Rudcbpy. The calculations predicted a maximum emission at about 540 nm for the complex constructed from two carbonyl π-acceptors ligands trans to the dcbpy, while an emission in the far infrared region is calculated when two H σ-donor ligands trans to the dcbpy. Our calculation results show correlations between HOMO-LUMO energy gap, Stokes shift, and T1 distortion, which reflect the different effects of electron-withdrawing and donating groups. We proposed that these correlations can be used to predict the photophysical properties for new complexes.

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.

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.

TD-DFT insight into photodissociation of the Co-C bond in coenzyme B12

Coenzyme B₁₂ (AdoCbl) is one of the most biologically active forms of vitamin B₁₂, and continues to be a topic of active research interest. The mechanism of Co-C bond cleavage in AdoCbl, and the corresponding enzymatic reactions are however, not well understood at the molecular level. In this work, time-dependent density functional theory (TD-DFT) has been applied to investigate the photodissociation of coenzyme B₁₂. To reduce computational cost, while retaining the major spectroscopic features of AdoCbl, a truncated model based on ribosylcobalamin (RibCbl) was used to simulate Co-C photodissociation. Equilibrium geometries of RibCbl were obtained by optimization at the DFT/BP86/TZVP level of theory, and low-lying excited states were calculated by TD-DFT using the same functional and basis set. The calculated singlet states, and absorption spectra were simulated in both the gas phase, and water, using the polarizable continuum model (PCM). Both spectra were in reasonable agreement with experimental data, and potential energy curves based on vertical excitations were plotted to explore the nature of Co-C bond dissociation. It was found that a repulsive ³(σ_Co−C → σ^*_Co−C) triplet state became dissociative at large Co-C bond distance, similar to a previous observation for methylcobalamin (MeCbl). Furthermore, potential energy surfaces (PESs) obtained as a function of both Co-C_Rib and Co-N_Im distances, identify the S₁ state as a key intermediate generated during photoexcitation of RibCbl, attributed to a mixture of a metal-to-ligand charge transfer (MLCT) and a σ bonding-ligand charge transfer (SBLCT) states.

Phenylvinylcobalamin: an alkenylcobalamin featuring a ligand with a large trans influence

Cob(I)alamin reacts with phenylacetylene to produce two diastereomers in which the organic ligand is coordinated to the upper (β) and lower (α) face of the corrin ring, respectively. The isomers were separated chromatographically and characterised by ESI-MS and, in the case of the β isomer, by (1)H and (13)C NMR. Only the β isomer crystallised and its molecular structure, determined by X-ray diffraction, shows that the organic ligand coordinates Co(III) through the β carbon of the phenylvinyl ligand. The Co-C bond length is 2.004(8) Å while the Co-N bond length to the trans 5,6-dimethylbenzimidazole (dmbzm) base is 2.217(8) Å, one of the longest Co-Ndmbzm bond lengths known in an organocobalamin. Unlike benzylcobalamin (BzCbl), phenylvinylcobalamin (PhVnCbl) is stable towards homolysis. DFT calculations (BP86/TZVP) on model compounds of BzCbl and PhVnCbl show that the Co-C bond dissociation energy for homolysis to Co(II) and an organic radical in the former is 8 kcal mol(-1) lower than in the latter. An analysis of the electron density at the Co-C bond critical point using Bader's QTAIM approach shows that the Co-C bond in PhVnCbl is shorter, stronger and somewhat more covalent than that in BzCbl, and has some multiple bond character. Together with calculations that show that the benzyl radical is more stable than the phenylvinyl radical, this rationalises the stability of PhVnCbl compared to BzCbl. The phenylvinyl ligand has a large trans influence. The pKa for deprotonation of dmbzm and its coordination by the metal in β-PhVnCbl is 4.60 ± 0.01, one of the highest values reported to date in cobalamin chemistry. The displacement of dmbzm ligand by CN(-) in β-PhVnCbl occurs with log K = 0.7 ± 0.1; the trans influence order of C-donor ligands is therefore CN(-) < CCH < CHCH2 = PhVn < Me < Et.

Observation of the inverse trans influence (ITI) in a uranium(V) imide coordination complex: an experimental study and theoretical evaluation

An inverse trans influence has been observed in a high-valent U(V) imide complex, [(((Ad)ArO)(3)N)U(NMes)]. A thorough theoretical evaluation has been employed in order to corroborate the ITI in this unusual complex. Computations on the target complex, [(((Ad)ArO)(3)N)U(NMes)], and the model complexes [(((Me)ArO)(3)N)U(NMes)] and [(NMe(3))(OMe(2))(OMe)(3)U(NPh)] are discussed along with synthetic details and supporting spectroscopic data. Additionally, the syntheses and full characterization data of the related U(V) trimethylsilylimide complex [(((Ad)ArO)(3)N)U(NTMS)] and U(IV) azide complex [(((Ad)ArO)(3)N)U(N(3))] are also presented for comparison.

Vitamin B12: unique metalorganic compounds and the most complex vitamins

The chemistry and biochemistry of the vitamin B(12) compounds (cobalamins, XCbl) are described, with particular emphasis on their structural aspects and their relationships with properties and function. A brief history of B(12), reveals how much the effort of chemists, biochemists and crystallographers have contributed in the past to understand the basic properties of this very complex vitamin. The properties of the two cobalamins, the two important B(12) cofactors Ado- and MeCbl are described, with particular emphasis on how the Co-C bond cleavage is involved in the enzymatic mechanisms. The main structural features of cobalamins are described, with particular reference to the axial fragment. The structure/property relationships in cobalamins are summarized. The recent studies on base-off/base-on equilibrium are emphasized for their relevance to the mode of binding of the cofactor to the protein scaffold. The absorption, transport and cellular uptake of cobalamins and the structure of the B(12) transport proteins, IF and TC, in mammals are reviewed. The B(12) transport in bacteria and the structure of the so far determined proteins are briefly described. The currently accepted mechanisms for the catalytic cycles of the AdoCbl and MeCbl enzymes are reported. The structure and function of B(12) enzymes, particularly the important mammalian enzymes methyltransferase (MetH) and methyl-malonyl-coenzyme A mutase (MMCM), are described and briefly discussed. Since fast proliferating cells require higher amount of vitamin B(12) than that required by normal cells, the study of B(12 )conjugates as targeting agents has recently gained importance. Bioconjugates have been studied as potential agents for delivering radioisotopes and NMR probes or as various cytotoxic agents towards cancer cells in humans and the most recent studies are described. Specifically, functionalized bioconjugates are used as "Trojan horses" to carry into the cell the appropriate antitumour or diagnostic label. Possible future developments of B(12) work are summarized.