Time-dependent density functional theory study of cobalt corrinoids: Electronically excited states of methylcobalamin

Vibronic Coupling in Vitamin B12: A Combined Spectroscopic and Computational Study

Understanding the diverse reactivities of vitamin B12 and its derivatives, collectively called cobalamins, requires detailed knowledge of their geometric and electronic structures. Electronic absorption (Abs) and resonance Raman (rR) spectroscopies have proven invaluable in this area, particularly when used in concert with computational techniques such as density functional theory (DFT). There remain, however, lingering uncertainties in the computational description of electronic excited states of cobalamins, particularly surrounding the vibronic coupling that impacts the Abs bandshapes and gives rise to rR enhancement of vibrational modes. Past computational analyses of the vibrational spectra of cobalamins have either neglected rR enhancement or calculated rR enhancement for only a small number of modes. In the present study, we used the recently developed ORCA_ASA computational tool in conjunction with the popular B3LYP and BP86 functionals to predict Abs bandshapes and rR spectra for vitamin B12. The ORCA_ASA/B3LYP-computed Abs envelope in the visible spectral region and rR spectra of vitamin B12 agree remarkably well with our experimental data, while BP86 fails to reproduce both. This finding represents a significant advance in our understanding of how these two commonly used density functionals differently model the electronic properties of cobalamins. Guided by the computed frequencies for the Co-C stretching and Co-C-N bending modes, we identified, for the first time, isotope-sensitive features in our rR spectra of 12CNCbl and 13CNCbl that can be assigned to these modes. A normal coordinate analysis of the experimentally determined Co-C stretching and Co-C-N bending frequencies indicates that the Co-C force constant for vitamin B12 is 2.67 mdyn/Å, considerably larger than the Co-C force constants reported for alkylcobalamins.

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.

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.

H2 Evolution from Electrocatalysts with Redox-Active Ligands: Mechanistic Insights from Theory and Experiment vis-à-vis Co-Mabiq

Electrocatalytic hydrogen production via transition metal complexes offers a promising approach for chemical energy storage. Optimal platforms to effectively control the proton and electron transfer steps en route to H₂ evolution still need to be established, and redox-active ligands could play an important role in this context. In this study, we explore the role of the redox-active Mabiq (Mabiq = 2-4:6-8-bis(3,3,4,4-tetramethlyldihydropyrrolo)-10-15-(2,2-biquinazolino)-[15]-1,3,5,8,10,14-hexaene1,3,7,9,11,14-N₆) ligand in the hydrogen evolution reaction (HER). Using spectro-electrochemical studies in conjunction with quantum chemical calculations, we identified two precatalytic intermediates formed upon the addition of two electrons and one proton to [Co^II(Mabiq)(THF)](PF₆) (Co_Mbq). We further examined the acid strength effect on the generation of the intermediates. The generation of the first intermediate, Co_Mbq-H¹, involves proton addition to the bridging imine-nitrogen atom of the ligand and requires strong proton activity. The second intermediate, Co_Mbq-H², acquires a proton at the diketiminate carbon for which a weaker proton activity is sufficient. We propose two decoupled H₂ evolution pathways based on these two intermediates, which operate at different overpotentials. Our results show how the various protonation sites of the redox-active Mabiq ligand affect the energies and activities of HER intermediates.

Resonance Raman Optical Activity Spectroscopy in Probing Structural Changes Invisible to Circular Dichroism Spectroscopy: A Study on Truncated Vitamin B12 Derivatives

This work demonstrates resonance Raman optical activity (RROA) spectra of three truncated vitamin B₁₂ derivatives modified within the nucleotide loop. Since truncated cobalamins possess sufficiently high rotational strength in the range of ROA excitation (532 nm), it was possible to record their spectra in the resonance condition. They showed several distinct spectral features allowing for the distinguishing of studied compounds, in contrast to other methods, i.e., UV-Vis absorption, electronic circular dichroism, and resonance Raman spectroscopy. The improved capacity of the RROA method is based here on the excitation of molecules via more than two electronic states, giving rise to the bisignate RROA spectrum, significantly distinct from a parent Raman spectrum. This observation is an important step in the dissemination of using RROA spectroscopy in studying the complex structure of corrinoids which may prove crucial for a better understanding of their biological role.

In vitro drug interaction of levocetirizine and diclofenac: Theoretical and spectroscopic studies

Levocetirizine dihydrochloride is known to interact with some anti-inflammatory drugs. We report here a comprehensive integrated theoretical and experimental study for the in vitro drug interaction between levocetirizine dihydrochloride (LEV) and diclofenac sodium (DIC). The interaction of the two drugs was confirmed by the molecular ion peak obtained from the mass spectrum of the product. Moreover, FTIR and ¹HNMR spectra of the individual drugs and their interaction product were inspected to allocate the possible sites of interaction. In addition, quantum mechanical DFT calculations were performed to search for the interaction sites and to verify the types of interactions deduced from the spectroscopic studies such as charge-transfer and non-bonding π-π interactions. It was found that the studied drugs interact with each other in aqueous solution via four types of interactions, namely, ion-pair formation, three weak hydrogen bonds, non-bonding π-π interactions and charge-transfer from DIC to LEV.

Enhancement of biocompatibility and photoacoustic contrast activity of metal clusters

Organometallic carbonyl clusters (OMCC) of group VIII elements are water soluble, bio-compatible and stable high-contrast photoacoustic agents for live cell imaging. But, they have limited application due to weak absorption within 700-1000nm wavelength which is known as the biological window of absorption. In this article, we report that hexa-nuclear iron (Fe₆) carbonyl cluster derivatized with sodium thio-propanoate has very good absorption within 700-1600nm wave length. This modeled compound is water soluble and bio-compatible. The bio-compatibility of this compound is tested through cytotoxicity, LogP and metabolic probability at CYP450-2D6 enzyme.

Energy cascades, excited state dynamics, and photochemistry in cob(III)alamins and ferric porphyrins

Porphyrins and the related chlorins and corrins contain a cyclic tetrapyrrole with the ability to coordinate an active metal center and to perform a variety of functions exploiting the oxidation state, reactivity, and axial ligation of the metal center. These compounds are used in optically activated applications ranging from light harvesting and energy conversion to medical therapeutics and photodynamic therapy to molecular electronics, spintronics, optoelectronic thin films, and optomagnetics. Cobalt containing corrin rings extend the range of applications through photolytic cleavage of a unique axial carbon-cobalt bond, permitting spatiotemporal control of drug delivery. The photochemistry and photophysics of cyclic tetrapyrroles are controlled by electronic relaxation dynamics including internal conversion and intersystem crossing. Typically the electronic excitation cascades through ring centered ππ* states, ligand to metal charge transfer (LMCT) states, metal to ligand charge transfer (MLCT) states, and metal centered states. Ultrafast transient absorption spectroscopy provides a powerful tool for the investigation of the electronic state dynamics in metal containing tetrapyrroles. The UV-visible spectrum is sensitive to the oxidation state, electronic configuration, spin state, and axial ligation of the central metal atom. Ultrashort broadband white light probes spanning the range from 270 to 800 nm, combined with tunable excitation pulses, permit the detailed unravelling of the time scales involved in the electronic energy cascade. State-of-the-art theoretical calculations provide additional insight required for precise assignment of the states. In this Account, we focus on recent ultrafast transient absorption studies of ferric porphyrins and corrin containing cob(III)alamins elucidating the electronic states responsible for ultrafast energy cascades, excited state dynamics, and the resulting photoreactivity or photostability of these compounds. Iron tetraphenyl porphyrin chloride (Fe((III))TPPCl) exhibits picosecond decay to a metal centered d → d* (4)T state. This state decays on a ca. 16 ps time scale in room temperature solution but persists for much longer in a cryogenic glass. The photoreactivity of the (4)T state may lead to novel future applications for these compounds. In contrast, the nonplanar cob(III)alamins contain two axial ligands to the central cobalt atom. The upper axial ligand can be an alkyl group as in the two biologically active coenzymes or a nonalkyl ligand such as -CN in cyanocobalamin (vitamin B12) or -OH in hydroxocobalamin. The electronic structure, energy cascade, and bond cleavage of these compounds is sensitive to the details of the axial ligand. Nonalkylcobalamins exhibit ultrafast internal conversion to a low-lying state of metal to ligand or ligand to metal charge transfer character. The compounds are generally photostable with ground state recovery complete on a time scale of 2-7 ps in room temperature aqueous solution. Alkylcobalamins exhibit ultrafast internal conversion to an S1 state of d/π → π* character. Most compounds undergo bond cleavage from this state with near unit quantum yield within ∼100 ps. Recent theoretical calculations provide a potential energy surface accounting for these observations. Conformation dependent mixing of the corrin π and cobalt d orbitals plays a significant role in the observed photochemistry and photophysics.

Computational investigation of the ligand field effect to improve the photoacoustic properties of organometallic carbonyl clusters

Organometallic nitrosyl carbonyl clusters are stable and better high-contrast photoacoustic contrast agents (PACAs) than organometallic carbonyl clusters.

Theoretical modeling of low-energy electronic absorption bands in reduced cobaloximes

The reduced Co(I) states of cobaloximes are powerful nucleophiles that play an important role in the hydrogen-evolving catalytic activity of these species. In this work we analyze the low-energy electronic absorption bands of two cobaloxime systems experimentally and use a variety of density functional theory and molecular orbital ab initio quantum chemical approaches. Overall we find a reasonable qualitative understanding of the electronic excitation spectra of these compounds but show that obtaining quantitative results remains a challenging task.

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.

Cob(I)alamin: insight into the nature of electronically excited states elucidated via quantum chemical computations and analysis of absorption, CD and MCD data

The nature of electronically excited states of the super-reduced form of vitamin B(12) (i.e., cob(I)alamin or B(12s)), a ubiquitous B(12) intermediate, was investigated by performing quantum-chemical calculations within the time-dependent density functional theory (TD-DFT) framework and by establishing their correspondence to experimental data. Using response theory, the electronic absorption (Abs), circular dichroism (CD) and magnetic CD (MCD) spectra of cob(I)alamin were simulated and directly compared with experiment. Several issues have been taken into considerations while performing the TD-DFT calculations, such as strong dependence on the applied exchange-correlation (XC) functional or structural simplification imposed on the cob(I)alamin. In addition, the low-lying transitions were also validated by performing CASSCF/MC-XQDPT2 calculations. By comparing computational results with existing experimental data a new level of understanding of electronic excitations has been established at the molecular level. The present study extends and confirms conclusions reached for other cobalamins. In particular, the better performance of the BP86 functional, rather than hybrid-type, was observed in terms of the excitations associated with both Co d and corrin π localized transitions. In addition, the lowest energy band was associated with multiple metal-to-ligand charge transfer excitations as opposed to the commonly assumed view of a single π → π* transition followed by vibrational progression. Finally, the use of the full cob(I)alamin structure, instead of simplified molecular models, shed new light on the spectral analyses of cobalamin systems and revealed new challenges of this approach related to long-range charge transfer excitations involving side chains.