Photoproduct formation in coenzyme B12-dependent CarH photoreceptor via a triplet pathway

CarH is a cobalamin-based photoreceptor which has attracted significant interest due to its complex mechanism involving its organometallic coenzyme-B₁₂ chromophore. While several experimental and computational studies have sought to understand CarH's mechanism of action, there are still many aspects of the mechanism which remain unclear. While light is needed to activate the Co-C_5' bond, it is not entirely clear whether reaction pathway involves singlet or triplet diradical states. A recent experimental study implicated triplet pathway and importance of intersystem crossing (ISC) as a viable mechanistic route for photoproduct formation in CarH. Herein, a combined quantum mechanics/molecular mechanics approach (QM/MM) was used to explore the involvement of triplet states in CarH. Two possibilities were explored. The first possibility involved photo-induced homolytic cleavage of the Co-C_5' where the radical pair (RP) would deactivate to a triplet state (T₀) on the ground state potential energy surface (PES). However, a pathway for the formation of the photoproduct, 4',5'-anhydroadenosine (anhAdo), on the triplet ground state PES was not energetically feasible. The second possibility involved exploring a manifold of low-lying triplet excited states computed using TD-DFT within the QM/MM framework. Viable crossings of triplet excited states with singlet excited states were identified using semiclassical Landau-Zener theory and the effectiveness of spin-orbit coupling by El-Sayed rules. Several candidates along both the Co-N_Im potential energy curve (PEC) and Co-C_5'/Co-N_Im PES were identified, which appear to corroborate experimental findings and implicate the possible role of triplet states in CarH.

Theoretical modeling of the singlet-triplet spin transition in different Ni(II)-diketo-pyrphyrin-based metal-ligand octahedral complexes

The structural stability, charge transfer effects and strength of the spin-orbit couplings in different Ni(ii)-ligand complexes have been studied at the DFT (B3LYP and CAM-B3LYP) and coupled cluster (DLPNO-CCSD(T)) levels of theory. Accordingly, two different, porphyrin- and diketo-pyrphyrin-based four-coordination macrocycles as planar ligands as well as pyridine (or pyrrole) and mesylate anion molecular groups as vertical ligands were considered in order to build metal-organic complexes with octahedral coordination configurations. For each molecular system, the identification of equilibrium geometries and the intersystem crossing (the minimum energy crossing) points between the potential energy surfaces of the singlet and triplet spin states is followed by computing the spin-orbit couplings between the two spin states. Structures, based on the diketo-pyrphyrin macrocycle as the planar ligand, show stronger six-coordination metal-organic complexes due to the extra electrostatic interaction between the positively charged central metal cation and the negatively charged vertical ligands. The results also show that the magnitude of the spin-orbit coupling is influenced by the atomic positions of deprotonations of the ligands, and implicitly the direction of the charge transfer between the ligand and the central metal ion.

Effect of the Electron Density of the Heme Fe Atom on the Nature of Fe-O2 Bonding in Oxy Myoglobin

Mössbauer spectroscopy has been used to characterize oxygenated myoglobins (oxy Mbs) reconstituted with native and chemically modified ⁵⁷Fe-enriched heme cofactors with different electron densities of the heme Fe atom (ρ_Fe) and to elucidate the effect of a change in the ρ_Fe on the nature of the bond between heme Fe and oxygen (O₂), i.e., the Fe-O₂ bond, in the protein. Quadrupole splitting (ΔE_Q) was found to decrease with decreasing ρ_Fe, and the observed ρ_Fe-dependent ΔE_Q confirmed an increase in the contribution of the ferric-superoxide (Fe³⁺-O₂^-) form to the resonance hybrid of the Fe-O₂ fragment with decreasing ρ_Fe. These observations explicitly accounted for the lowering of O₂ affinity of the protein due to an increase in the O₂ dissociation rate and a decrease in the autoxidation reaction rate of oxy Mb through decreasing H⁺ affinity of the bound ligand with decreasing ρ_Fe. Therefore, the present study demonstrated the mechanism underlying the electronic control of O₂ affinity and the autoxidation of the protein through the heme electronic structure. Carbon monoxide (CO) adducts of reconstituted Mbs (CO-Mbs) were similarly characterized, and we found that the resonance between the two canonical forms of the Fe-CO fragment was also affected by a change in ρ_Fe. Thus, the nature of the Fe-ligand bond in the protein was found to be affected by the ρ_Fe.

Spin-inversion mechanisms in O2 binding to a model heme compound: A perspective from nonadiabatic wave packet calculations

Spin-inversion dynamics in O₂ binding to a model heme complex, which consisted of Fe(II)-porphyrin and imidazole, were studied using nonadiabatic wave packet dynamics calculations. We considered three active nuclear degrees of freedom in the dynamics, including the motions along the Fe-O distance, Fe-O-O angle, and Fe out-of-plane distance. Spin-free potential energy surfaces for the singlet, triplet, quintet, and septet states were developed using density functional theory calculations, and spin-orbit coupling elements were obtained from CASSCF-level electronic structure calculations. The spin-inversion mainly occurred between the singlet state and one of the triplet states due to large spin-orbit couplings and the contributions of other states were extremely small. The present quantum dynamics calculations suggested that the narrow crossing region model plays a dominant role in the O₂ binding dynamics. In addition, the one-dimensional Landau-Zener model underestimated the nonadiabatic transition probability.

Spin-inversion mechanisms in O2 binding to a model heme complex revisited by density function theory calculations

Spin-inversion mechanisms in O₂ binding to a model heme complex, consisting of Fe(II)-porphyrin and imidazole, were investigated using density-functional theory calculations. First, we applied the recently proposed mixed-spin Hamiltonian method to locate spin-inversion structures between different total spin multiplicities. Nine spin-inversion structures were successfully optimized for the singlet-triplet, singlet-quintet, triplet-quintet, and quintet-septet spin-inversion processes. We found that the singlet-triplet spin-inversion points are located around the potential energy surface region at short Fe-O distances, whereas the singlet-quintet and quintet-septet spin-inversion points are located at longer Fe-O distances. This suggests that both narrow and broad crossing models play roles in O₂ binding to the Fe-porphyrin complex. To further understand spin-inversion mechanisms, we performed on-the-fly Born-Oppenheimer molecular dynamics calculations. The reaction coordinates, which are correlated to the spin-inversion dynamics between different spin multiplicities, are also discussed.

First Principles Calculation of the Reaction Rates for Ligand Binding to Myoglobin: The Cases of NO and CO

Ligand binding by proteins is among the most fundamental processes in nature. Among these processes the binding of small gas molecules, such as O₂ , CO and NO to heme proteins has traditionally received vivid interest, which was further boosted by their recently recognized significant role in gas sensing in the body. At the heart of the binding of these ligands to the heme group is the spinforbidden reaction between high-spin iron(II) and the ligand yielding a low-spin adduct. We use computational means to address the complete mechanism of CO and NO binding by myoglobin. Considering that it involves several steps occurring on different time scales, molecular dynamics simulations were performed to address the diffusion of the ligand through the enzyme, and DFT calculations in combination with statistical rate calculation to investigate the spin-forbidden reaction. The calculations yielded rate constants in qualitative agreement with experiments and revealed that the bottleneck of NO and CO binding is different; for NO, diffusion was found to be rate-limiting, whereas for CO, the spin-forbidden step is the slowest.

Heme isomers substantially affect heme's electronic structure and function

Inspection of heme protein structures in the protein data bank reveals four isomers of heme characterized by different relative orientations of the vinyl side chains; remarkably, all these have been reported in multiple protein structures. Density functional theory computations explain this as due to similar energy of the isomers but with a sizable (25 kJ mol^-1) barrier to interconversion arising from restricted rotation around the conjugated bonds. The four isomers, EE, EZ, ZE, and ZZ, were then investigated as 4-coordinate hemes, as 5-coordinate deoxyhemes, in 6-coordinate O₂-adducts of globins and as compound I intermediates typical of heme peroxidases. Substantial differences were observed in electronic properties relevant to heme function: notably, the spin state energy gap of O₂-heme adducts, important for fast reversible binding of O₂, depends on the isomer state, and O₂-binding enthalpies change by up to 16 kJ mol^-1; redox potentials change by up to 0.2 V depending on the isomer, and the doublet-quartet energy splitting of compound I, central to "two-state" reactivity, is affected by up to ∼15 kJ mol^-1. These effects are consistently seen with three distinct density functionals, i.e. the effects are not method-dependent. Thus, the nature of the isomer state is an important but overlooked feature of heme chemistry and function, and previous and future studies of hemes may be reconsidered in this new context.

Spin crossover dynamics studies on the thermally activated molecular oxygen binding mechanism on a model copper complex

The theoretical description of the primary dioxygen (O₂) binding and activation step in many copper or iron enzymes, suffers from the instrinsically electronic non-adiabaticity of the spin flip events of the triplet dioxygen molecule (³O₂), mediated by spin–orbit couplings.

Empirical Force Fields for Mechanistic Studies of Chemical Reactions in Proteins

Following chemical reactions in atomistic detail is one of the most challenging aspects of current computational approaches to chemistry. In this chapter the application of adiabatic reactive MD (ARMD) and its multistate version (MS-ARMD) are discussed. Both methods allow to study bond-breaking and bond-forming processes in chemical and biological processes. Particular emphasis is put on practical aspects for applying the methods to investigate the dynamics of chemical reactions. The chapter closes with an outlook of possible generalizations of the methods discussed.

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.

A quantitative model for oxygen uptake and release in a family of hemeproteins

MOTIVATION

Hemeproteins have many diverse functions that largely depend on the rate at which they uptake or release small ligands, like oxygen. These proteins have been extensively studied using either simulations or experiments, albeit only qualitatively and one or two proteins at a time.

RESULTS

We present a physical-chemical model, which uses data obtained exclusively from computer simulations, to describe the uptake and release of oxygen in a family of hemeproteins, called truncated hemoglobins (trHbs). Through a rigorous statistical analysis we demonstrate that our model successfully recaptures all the reported experimental oxygen association and dissociation kinetic rate constants, thus allowing us to establish the key factors that determine the rates at which these hemeproteins uptake and release oxygen. We found that internal tunnels as well as the distal site water molecules control ligand uptake, whereas oxygen stabilization by distal site residues controls ligand release. Because these rates largely determine the functions of these hemeproteins, these approaches will also be important tools in characterizing the trHbs members with unknown functions.

CONTACT

lboechi@ic.fcen.uba.ar

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A DFT and multi-configurational perturbation theory study on O2 binding to a model heme compound via the spin-change barrier

Potential energy surface of the O₂ binding to a model heme compound via minimum energy intersystem crossing point (MEISCP) was investigated with DFT and MS-CASPT2. The porphyrin's symmetric shrinking vibration mode contributes to reach the MESICP from the T₁ minimum.

Abrupt versus Gradual Spin-Crossover in Fe(II)(phen)2(NCS)2 and Fe(III)(dedtc)3 Compared by X-ray Absorption and Emission Spectroscopy and Quantum-Chemical Calculations

Molecular spin-crossover (SCO) compounds are attractive for information storage and photovoltaic technologies. We compared two prototypic SCO compounds with Fe(II)N6 (1, [Fe(phen)2(NCS)2], with phen = 1,10-phenanthroline) or Fe(III)S6 (2, [Fe(dedtc)3], with dedtc = N,N'-diethyldithiocarbamate) centers, which show abrupt (1) or gradual (2) thermally induced SCO, using K-edge X-ray absorption and Kβ emission spectroscopy (XAS/XES) in a 8-315 K temperature range, single-crystal X-ray diffraction (XRD), and density functional theory (DFT). Core-to-valence and valence-to-core electronic transitions in the XAS/XES spectra and bond lengths change from XRD provided benchmark data, verifying the adequacy of the TPSSh/TZVP DFT approach for the description of low-spin (LS) and high-spin (HS) species. Determination of the spin densities, charge distributions, bonding descriptors, and valence-level configurations, as well as similar experimental and calculated enthalpy changes (ΔH), suggested that the varying metal-ligand bonding properties and deviating electronic structures converge to similar enthalpic contributions to the free-energy change (ΔG) and thus presumably are not decisive for the differing SCO behavior of 1 and 2. Rather, SCO seems to be governed by vibrational contributions to the entropy changes (ΔS) in both complexes. Intra- and intermolecular interactions in crystals of 1 and 2 were identified by atoms-in-molecules analysis. Thermal excitation of individual dedtc ligand vibrations accompanies the gradual SCO in 2. In contrast, extensive inter- and intramolecular phen/NCS vibrational mode coupling may be an important factor in the cooperative SCO behavior of 1.

A two-in-one pincer ligand and its diiron(II) complex showing spin state switching in solution through reversible ligand exchange

A novel pyrazolate-bridged ligand providing two {PNN} pincer-type compartments has been synthesized. Its diiron(II) complex LFe2(OTf)3(CH3CN) (1; Tf = triflate) features, in solid state, two bridging triflate ligands, with a terminal triflate and a MeCN ligand completing the octahedral coordination spheres of the two high-spin metal ions. In MeCN solution, 1 is shown to undergo a sequential, reversible, and complete spin transition to the low-spin state upon cooling. Detailed UV/Vis and (19)F NMR spectroscopic studies as well as magnetic measurements have unraveled that spin state switching correlates with a rapid multistep triflate/MeCN ligand exchange equilibrium. The spin transition temperature can be continuously tuned by varying the triflate concentration in solution.

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

Quantum-chemical cluster modeling is challenged in the limit of large, soft systems by the effects of dispersion and solvent, and well as other physical interactions. Adenosylcobalamin (AdoCbl, coenzyme B12), as one of the most complex cofactors in life, constitutes such a challenge. The cleavage of its unique organometallic Co-C bond has inspired multiple studies of this cofactor. This paper reports the fully relaxed potential energy surface of Co-C cleavage of AdoCbl, including for the first time all side-chain interactions with the dissociating Ado group. Various methods and corrections for dispersion, relativistic effects, solvent polarity, basis set superposition error, and thermal and vibrational effects were investigated, totaling more than 550 single-point energies for the large model. The results show immense variability depending on method, including solvation, functional type, and dispersion, challenging the conceived accuracy of methods used for such systems. In particular, B3LYP-D3 seems to severely underestimate the Co-C bond strength, consistent with previous results, and BP86 remains accurate for cobalamins when dispersion interactions are accounted for.