Photodissociation of Co−C Bond in Methyl- and Ethylcobalamin: An Insight from TD-DFT Calculations

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.

Photoproduct formation in coenzyme B12-dependent CarH via a singlet pathway

The CarH photoreceptor exploits of the light-sensing ability of coenzyme B₁₂ ( adenosylcobalamin = AdoCbl) to perform its catalytic function, which includes large-scale structural changes to regulate transcription. In daylight, transcription is activated in CarH via the photo-cleavage of the Co-C_5' bond of coenzyme B₁₂. Subsequently, the photoproduct, 4',5'-anhydroadenosine (anhAdo) is formed inducing dissociation of the CarH tetramer from DNA. Several experimental studies have proposed that hydridocoblamin (HCbl) may be formed in process with anhAdo. The photolytic cleavage of the Co-C_5' bond of AdoCbl was previously investigated using photochemical techniques and the involvement of both singlet and triplet excited states were explored. Herein, QM/MM calculations were employed to probe (1) the photolytic processes which may involve singlet excited states, (2) the mechanism of anhAdo formation, and (3) whether HCbl is a viable intermediate in CarH. Time-dependent density functional theory (TD-DFT) calculations indicate that the mechanism of photodissociation of the Ado ligand involves the ligand field (LF) portion of the lowest singlet excited state (S₁) potential energy surface (PES). This is followed by deactivation to a point on the S₀ PES where the Co-C_5' bond remains broken. This species corresponds to a singlet diradical intermediate. From this point, the PES for anhAdo formation was explored, using the Co-C_5' and Co-C_4' bond distances as active coordinates, and a local minimum representing anhAdo and HCbl formation was found. The transition state (TS) for the formation of the Co-H bond of HCbl was located and its identity was confirmed by a single imaginary frequency of i1592 cm^-1. Comparisons to experimental studies and the potential role of rotation around the N-glycosidic bond of the Ado ligand were 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.

Visible-to-NIR-Light Activated Release: From Small Molecules to Nanomaterials

Photoactivatable (alternatively, photoremovable, photoreleasable, or photocleavable) protecting groups (PPGs), also known as caged or photocaged compounds, are used to enable non-invasive spatiotemporal photochemical control over the release of species of interest. Recent years have seen the development of PPGs activatable by biologically and chemically benign visible and near-infrared (NIR) light. These long-wavelength-absorbing moieties expand the applicability of this powerful method and its accessibility to non-specialist users. This review comprehensively covers organic and transition metal-containing photoactivatable compounds (complexes) that absorb in the visible- and NIR-range to release various leaving groups and gasotransmitters (carbon monoxide, nitric oxide, and hydrogen sulfide). The text also covers visible- and NIR-light-induced photosensitized release using molecular sensitizers, quantum dots, and upconversion and second-harmonic nanoparticles, as well as release via photodynamic (photooxygenation by singlet oxygen) and photothermal effects. Release from photoactivatable polymers, micelles, vesicles, and photoswitches, along with the related emerging field of photopharmacology, is discussed at the end of the review.

Ultrafast XANES Monitors Femtosecond Sequential Structural Evolution in Photoexcited Coenzyme B12

Polarized X-ray absorption near-edge structure (XANES) at the Co K-edge and broadband UV-vis transient absorption are used to monitor the sequential evolution of the excited-state structure of coenzyme B₁₂ (adenosylcobalamin) over the first picosecond following excitation. The initial state is characterized by sub-100 fs sequential changes around the central cobalt. These are polarized first in the y-direction orthogonal to the transition dipole and 50 fs later in the x-direction along the transition dipole. Expansion of the axial bonds follows on a ca. 200 fs time scale as the molecule moves out of the Franck-Condon active region of the potential energy surface. On the same 200 fs time scale there are electronic changes that result in the loss of stimulated emission and the appearance of a strong absorption at 340 nm. These measurements provide a cobalt-centered movie of the excited molecule as it evolves to the local excited-state minimum.

Photochemical Spin Dynamics of the Vitamin B12 Derivative, Methylcobalamin

Derivatives of vitamin B₁₂ are six-coordinate cobalt corrinoids found in humans, other animals, and microorganisms. By acting as enzymatic cofactors and photoreceptor chromophores, they serve vital metabolic and photoprotective functions. Depending on the context, the chemical mechanisms of the biologically active derivatives of B₁₂-methylcobalamin (MeCbl) and 5'-deoxyadenosylcobalamin (AdoCbl)-can be very different from one another. The extent to which this chemistry is tuned by the upper axial ligand, however, is not yet clear. Here, we have used a combination of time-resolved Fourier transform-electron paramagnetic resonance (FT-EPR), magnetic field effect experiments, and spin dynamic simulations to reveal that the upper axial ligand alone only results in relatively minor changes to the photochemical spin dynamics of B₁₂. By studying the photolysis of MeCbl, we find that, similar to AdoCbl, the initial (or "geminate") radical pairs (RPs) are born predominantly in the singlet spin state and thus originate from singlet excited-state precursors. This is in contrast to the triplet RPs and precursors proposed previously. Unlike AdoCbl, the extent of geminate recombination is limited following MeCbl photolysis, resulting in significant distortions to the FT-EPR signal caused by polarization from spin-correlated methyl-methyl radical "f-pairs" formed following rapid diffusion. Despite the photophysical mechanism that precedes photolysis of MeCbl showing wavelength dependence, the subsequent spin dynamics appear to be largely independent of excitation wavelength, again similar to AdoCbl. Our data finally provide clarity to what in the literature to date has been a confused and contradictory picture. We conclude that, although the upper axial position of MeCbl and AdoCbl does impact their reactivity to some extent, the remarkable biochemical diversity of these fascinating molecules is most likely a result of tuning by their protein environment.

How does the mutation in the cap domain of methylcobalamin-dependent methionine synthase influence the photoactivation of the Co–C bond?

The topology of the S₁ PES is modulated by introducing a mutation at the F708 position. The mutation influences the photoactivation of the Co–C bond by decreasing the rate of geminate recombination and altering the rate of radical pair formation.

The photochemistry and photobiology of vitamin B12

This Perspective provides the first detailed overview of the photoresponse of vitamin B₁₂ and its derivatives, from the early, photophysical events to the burgeoning area of B₁₂-dependent photobiology.

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.

Fluorophore Assisted Photolysis of Thiolato-Cob(III)alamins

Cobalamins are known to react with thiols to yield stable β-axial Co(III)-S bonded thiolato-cobalamin complexes. However, in stark contrast to the Co-C bond in alkylcobalamins, the photolability of the Co-S bond in thiolato-cobalamins remains undetermined. We have investigated the photolysis of N-acetylcysteinyl cob(III)alamin at several wavelengths within the ultraviolet and visible spectrum. To aid in photolysis, we show that attaching fluorophore "antennae" to the cobalamin scaffold can improve photolytic efficiency by up to an order of magnitude. Additionally, electron paramagnetic resonance confirms previous conjectures that the photolysis of thiolato-cobalamins at wavelengths as long as 546 nm produces thiyl radicals.

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.

Theoretical spectroscopy and photodynamics of a ruthenium nitrosyl complex

Photoactive transition-metal nitrosyl complexes are particularly interesting as potential drugs that deliver nitric oxide (NO) upon UV-light irradiation to be used, e.g., in photodynamic therapy. It is well-recognized that quantum-chemical calculations can guide the rational design and synthesis of molecules with specific functions. In this contribution, it is shown how electronic structure calculations and dynamical simulations can provide a unique insight into the photodissociation mechanism of NO. Exemplarily, [Ru(PaPy3)(NO)](2+) is investigated in detail, as a prototype of a particularly promising class of photoactive metal nitrosyl complexes. The ability of time-dependent density functional theory (TD-DFT) to obtain reliable excited-state energies compared with more sophisticated multiconfigurational spin-corrected calculations is evaluated. Moreover, a TD-DFT-based trajectory surface-hopping molecular dynamics study is employed to reveal the details of the radiationless decay of the molecule via internal conversion and intersystem crossing. Calculations show that the ground state of [Ru(PaPy3)(NO)](2+) includes a significant admixture of the Ru(III)(NO)(0) electronic configuration, in contrast to the previously postulated Ru(II)(NO)(+) structure of similar metal nitrosyls. Moreover, the lowest singlet and triplet excited states populate the antibonding metal d → πNO* orbitals, favoring NO dissociation. Molecular dynamics show that intersystem crossing is ultrafast (<10 fs) and dissociation is initiated in less than 50 fs. The competing relaxation to the lowest S1 singlet state takes place in less than 100 fs and thus competes with NO dissociation, which mostly takes place in the higher-lying excited triplet states. All of these processes are accompanied by bending of the NO ligand, which is not confined to any particular state.

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.

Quantum mechanical methods for the investigation of metalloproteins and related bioinorganic compounds

It is well known that transition metal ions are often bound to proteins, conveying very specific functional properties. In fact, metalloproteins play crucial biological roles in the transport and activation of small molecules such as H2, O2, and N2, as well as in several other biochemical processes. However, even if the presence of transition metals in the active site of proteins allows a very rich biochemistry, the experimental disclosure of structure-activity relationships in metalloproteins is generally difficult exactly because of the presence of transition metals, which are intrinsically characterized by a very versatile and often elusive chemistry. For this reason, computational methods are becoming very popular tools in the characterization of metalloproteins. In particular, since computing power is becoming less and less expensive, due to the continuous technological development of CPUs, the computational tools suited to investigate metalloproteins are becoming more accessible and therefore more commonly used also in molecular biology and biochemistry laboratories. Here, we present the main procedures and computational methods based on quantum mechanics, which are commonly used to study the structural, electronic, and reactivity properties of metalloproteins and related bioinspired compounds, with a specific focus on the practical and technical aspects that must be generally tackled to properly study such biomolecular systems.

Theoretical study on the photodegradation mechanism of nona-BDEs in methanol

Polybrominated diphenyl ethers (PBDEs) have received special environmental concern due to their potential toxicity to humans and wildlife worldwide, however, it is difficult to reveal their dominant photochemical degradation pathways by experiment. We explored the reaction mechanisms of photochemical degradation-debromination of three nona-BDEs in methanol using theoretical calculations, in which time-dependent density functional theory (TDDFT) combined with the polarizable continuum (PCM) model is applied. The selectivity of debromination was studied, and the major octa-BDE products photochemically debrominated from nona-BDEs were identified. We find that the debromination reaction results from the electronic transitions from π to σ* orbitals when nona-BDEs are exposed to UV-light in the sunlight region, at which point the two low-lying excited states for each nona-BDE are πσ*(5Br) and πσ*(4Br), which correlate to the σ* orbitals located on the penta-Br and tetra-Br substituted phenyls, respectively. Our calculations indicate that each nona-BDE may degrade to form three kinds of octa-BDE products via the πσ*(5Br) state, whereas only one kind of octa-BDEs can be formed via the πσ*(4Br) state. Our calculations can interpret the recent experiments successfully.

Theoretical Investigation of the Electronic Structure of Fe(II) Complexes at Spin-State Transitions

The electronic structure relevant to low spin (LS)↔high spin (HS) transitions in Fe(II) coordination compounds with a FeN₆ core are studied. The selected [Fe(tz)₆]²⁺ (1) (tz = 1H-tetrazole), [Fe(bipy)₃]²⁺ (2) (bipy = 2,2′-bipyridine), and [Fe(terpy)₂]²⁺ (3) (terpy = 2,2′:6′,2″-terpyridine) complexes have been actively studied experimentally, and with their respective mono-, bi-, and tridentate ligands, they constitute a comprehensive set for theoretical case studies. The methods in this work include density functional theory (DFT), time-dependent DFT (TD-DFT), and multiconfigurational second order perturbation theory (CASPT2). We determine the structural parameters as well as the energy splitting of the LS–HS states (ΔE_HL) applying the above methods and comparing their performance. We also determine the potential energy curves representing the ground and low-energy excited singlet, triplet, and quintet d⁶ states along the mode(s) that connect the LS and HS states. The results indicate that while DFT is well suited for the prediction of structural parameters, an accurate multiconfigurational approach is essential for the quantitative determination of ΔE_HL. In addition, a good qualitative agreement is found between the TD-DFT and CASPT2 potential energy curves. Although the TD-DFT results might differ in some respect (in our case, we found a discrepancy at the triplet states), our results suggest that this approach, with due care, is very promising as an alternative for the very expensive CASPT2 method. Finally, the two-dimensional (2D) potential energy surfaces above the plane spanned by the two relevant configuration coordinates in [Fe(terpy)₂]²⁺ were computed at both the DFT and CASPT2 levels. These 2D surfaces indicate that the singlet–triplet and triplet–quintet states are separated along different coordinates, i.e., different vibration modes. Our results confirm that in contrast to the case of complexes with mono- and bidentate ligands, the singlet–quintet transitions in [Fe(terpy)₂]²⁺ cannot be described using a single configuration coordinate.

Photolability of NO in designed metal nitrosyls with carboxamido-N donors: a theoretical attempt to unravel the mechanism

During the past few years, photoactive metal nitrosyls (NO complexes of metals) have drawn attention as potential drugs for delivery of nitric oxide (NO) to biological targets under the control of light. Major success in this area has been achieved with designed metal nitrosyls derived from ligands that contain carboxamide group(s). A number of iron, manganese and ruthenium {MNO}(6) nitrosyls of such kind exhibit excellent NO photolability under low-power visible and near-IR light. The results of theoretical studies on these NO-donors have provided insight into (a) the electronic transitions that lead to photorelease of NO and (b) the structural features of the ligands that dictate the sensitivity of the nitrosyls to light of specific wavelengths. In addition, the results have afforded clear understanding of the electronic configurations of the various nitrosyls. This article highlights these results in a coherent manner. Good matches between the predicted and observed spectral features and NO photolability strongly suggest that theoretical studies should be an integral part of the smart design of such NO-donors in the future research.

Continuous wave photolysis magnetic field effect investigations with free and protein-bound alkylcobalamins

The activation of the Co-C bond in adenosylcobalamin-dependent enzymes generates a singlet-born Co(II)-adenosyl radical pair. Two of the salient questions regarding this process are: (1) What is the origin of the considerable homolysis rate enhancement achieved by this class of enzyme? (2) Are the reaction dynamics of the resultant radical pair sensitive to the application of external magnetic fields? Here, we present continuous wave photolysis magnetic field effect (MFE) data that reveal the ethanolamine ammonia lyase (EAL) active site to be an ideal microreactor in which to observe enhanced magnetic field sensitivity in the adenosylcobalamin radical pair. The observed field dependence is in excellent agreement with that calculated from published hyperfine couplings for the constituent radicals, and the magnitude of the MFE (<18%) is almost identical to that observed in a solvent containing 67% glycerol. Similar augmentation is not observed, however, in the equivalent experiments with EAL-bound methylcobalamin, where all field sensitivity observed in the free cofactor is washed out completely. Parallels are drawn between the latter case and the loss of field sensitivity in the EAL holoenzyme upon substrate binding (Jones et al. J. Am. Chem. Soc. 2007, 129, 15718-15727). Both are attributed to the rapid removal of the alkyl radical immediately after homolysis, such that there is inadequate radical pair recombination for the observation of field effects. Taken together, these results support the notion that rapid radical quenching, through the coupling of homolysis and hydrogen abstraction steps, and subsequent radical pair stabilization make a contribution to the observed rate acceleration of Co-C bond homolysis in adenosylcobalamin-dependent enzymes.