Mechanism of the highly effective peptide bond hydrolysis by MOF-808 catalyst under biologically relevant conditions

Efficient and selective hydrolysis of inert peptide bonds is of paramount importance. MOF-808, a metal-organic framework based on Zr₆ nodes, can hydrolyze peptide bonds efficiently under biologically relevant conditions. However, the details of the catalyst structure and of the underlying catalytic reaction mechanism are challenging to establish. By means of DFT calculations we first investigate the speciation of the Zr₆ nodes and identify the nature of ligands that bind to the Zr₆O₈H_4-x core in aqueous conditions. The core is predicted to strongly prefer a Zr₆O₈H₄ protonation state and to be predominantly decorated by bridging formate ligands, giving Zr₆(μ₃-O)₄(μ₃-OH)₄(BTC)₂(HCOO)₆ and Zr₆(μ₃-O)₄(μ₃-OH)₄(BTC)₂(HCOO)₅(OH)(H₂O) as the most favorable structures at physiological pH. The GlyGly peptide can bind MOF in several different ways, with the preferred structure involving coordination through the terminal carboxylate analogously to the binding mode of formate ligand. The pre-reactive binding mode in which the amide carbonyl oxygen coordinates the metal core lies 7 kcal higher in free energy. The preferred reaction pathway is predicted to have two close-lying transition states, either of which could be the rate-determining step: nucleophilic attack on the amide carbon atom and C-N bond breaking, with calculated relative free energies of 31 and 32 kcal mol^-1, respectively. Replacement of formate by water and hydroxide at the Zr₆ node is predicted to be possible, but does not appear to play a role in the hydrolysis mechanism.

Modern quantum chemistry with [Open]Molcas

MOLCAS/OpenMolcas is an ab initio electronic structure program providing a large set of computational methods from Hartree-Fock and density functional theory to various implementations of multiconfigurational theory. This article provides a comprehensive overview of the main features of the code, specifically reviewing the use of the code in previously reported chemical applications as well as more recent applications including the calculation of magnetic properties from optimized density matrix renormalization group wave functions.

Advances in the synthesis, characterisation, and mechanistic understanding of active sites in Fe-zeolites for redox catalysts

The recent research developments on the active sites in Fe-zeolites for redox catalysis are discussed. Building on the characterisation of the α-Fe/α-O active sites in the beta and chabazite zeolites, we demonstrate a bottom-up approach to successfully understand and develop Fe-zeolite catalysts. We use the room temperature benzene to phenol reaction as a relevant example. We then suggest how the spectroscopic identification of other monomeric and dimeric iron sites could be tackled. The challenges in the characterisation of active sites and intermediates in NOX selective catalytic reduction catalysts and further development of catalysts for mild partial methane oxidation are briefly discussed.

Electronic Structure of N-Bridged High-Valent Diiron-Oxo

Density functional theory (DFT) and an advanced ab initio technique based on density matrix renormalization group (DMRG-CASPT2) were employed to investigate a reactive N-bridged high-valent diiron-oxo species involved in H-abstraction reactions. We studied in detail two important doublet states, the ground state with two iron(IV) centers and a mixed valence Fe^V -Fe^IV excited state. We found that the latter state is low-lying. Furthermore, its electronic structure and spin density imply that it has significantly higher H-abstraction reactivity than the ground state. This low-lying excited state might be the reason behind the high oxidation reactivity of this diiron-oxo species towards methane.

Low-Lying Electromeric States in Chloro-Ligated Iron(IV)-Oxo Porphyrin as a Model for Compound I, Studied with Second-Order Perturbation Theory Based on Density Matrix Renormalization Group

Employing second-order perturbation theory based on the density matrix renormalization group (DMRG-CASPT2), this work aims at providing a quantitative description of the spin state energetics of a chloro-ligated iron(IV)-oxo porphyrin as a model for the cytochromes P450 active species, also known as compound I (Cpd I). We explored DMRG-CASPT2 to its full extent with an extensive active space (up to 31 active orbitals) as well as a large number of renormalized states m (up to 10000). Different flavors of DMRG-CASPT2, using either the costly exact 4-particle reduced density matrix (4-RDM) or the cheaper cumulant approximated 4-RDM (cu(4)), were analyzed. All flavors essentially converge to similar relative energies between different spin states. Including a correction for the protein environment, we found a quartet Fe^IVO ground state and, more importantly, a thermally accessible doublet Fe^VO excited state that might directly contribute to the reactivity of this iron-oxo species. Our results also showed that cheaper approaches, such as CASPT2 based on a smaller active space or the cumulant approximation DMRG-cu(4)-CASPT2, are capable of accurately describing the spin state energetics of this species.

Limits of Coupled-Cluster Calculations for Non-Heme Iron Complexes

In a large variety of studies, the coupled-cluster method with singles, doubles, and perturbative triples (CCSD(T)) is used as a reference for benchmarking the performance of density functional theory (DFT) functionals. In the case of open-shell species, this theory can be applied in different forms depending on the restricted or unrestricted treatment of spin. In this study, we show that these different approaches can produce results which deviate by ∼5 kcal/mol for different species on the potential energy surfaces. This was demonstrated for a simple model of the C-H activation carried out by non-heme iron enzymes. Assessing the limits of CCSD(T) prior to its use as a general benchmark tool is warranted. This was done using higher-order coupled-cluster calculations as well as multiconfigurational second-order perturbation theory (CASPT2), since iron-oxo species present some multireference character. Furthermore, we tested two different implementations of the local coupled-cluster method and compared them to the CCSD(T) results, showing that even though these novel approaches are promising, without further developments they appear not to be suitable for describing two-state reactivity of the system investigated in the current study. Additionally, we implemented and assessed the performance of the hotspot approach for the local unrestricted CCSD(T) scheme which aims at reducing the pair error for systems containing transition metals.

Reaction Mechanism of [NiFe] Hydrogenase Studied by Computational Methods

[NiFe] hydrogenases catalyze the reversible conversion of molecular hydrogen to protons and electrons. This seemingly simple reaction has attracted much attention because of the prospective use of H₂ as a clean fuel. In this paper, we have studied the full reaction mechanism of this enzyme with various computational methods. Geometries were obtained with combined quantum mechanical and molecular mechanics (QM/MM) calculations. To get more accurate energies and obtain a detailed account of the surroundings, we performed big-QM calculations with 819 atoms in the QM region. Moreover, QM/MM thermodynamic cycle perturbation calculations were performed to obtain free energies. Finally, density matrix renormalisation group complete active space self-consistent field calculations were carried out to study the electronic structures of the various states in the reaction mechanism. Our calculations indicate that the Ni-L state is not involved in the reaction mechanism. Instead, the Ni-C state is reduced by one electron and then the bridging hydride ion is transferred to the sulfur atom of Cys546 as a proton and the two electrons transfer to the Ni ion. This step turned out to be rate-determining with an energy barrier of 58 kJ/mol, which is consistent with the experimental rate of 750 ± 90 s^-1 (corresponding to ∼52 kJ/mol). The cleavage of the H-H bond is facile with an energy barrier of 33 kJ/mol, according to our calculations. We also find that the reaction energies are sensitive to the size of the QM system, the basis set, and the density functional theory method, in agreement with previous studies.

Spectroscopic Identification of the α-Fe/α-O Active Site in Fe-CHA Zeolite for the Low-Temperature Activation of the Methane C-H Bond

The formation of single-site α-Fe in the CHA zeolite topology is demonstrated. The site is shown to be active in oxygen atom abstraction from N₂O to form a highly reactive α-O, capable of methane activation at room temperature to form methanol. The methanol product can subsequently be desorbed by online steaming at 200 °C. For the intermediate steps of the reaction cycle, the evolution of the Fe active site is monitored by UV-vis-NIR and Mössbauer spectroscopy. A B3LYP-DFT model of the α-Fe site in CHA is constructed, and the ligand field transitions are calculated by CASPT2. The model is experimentally substantiated by the preferential formation of α-Fe over other Fe species, the requirement of paired framework aluminum and a MeOH/Fe ratio indicating a mononuclear active site. The simple CHA topology is shown to mitigate the heterogeneity of iron speciation found on other Fe-zeolites, with Fe₂O₃ being the only identifiable phase other than α-Fe formed in Fe-CHA. The α-Fe site is formed in the d6r composite building unit, which occurs frequently across synthetic and natural zeolites. Finally, through a comparison between α-Fe in Fe-CHA and Fe-*BEA, the topology's 6MR geometry is found to influence the structure, the ligand field, and consequently the spectroscopy of the α-Fe site in a predictable manner. Variations in zeolite topology can thus be used to rationally tune the active site properties.

Protein-Assisted Formation and Stabilization of Catalytically Active Polyoxometalate Species

The effect of the protein environment on the formation and stabilization of an elusive catalytically active polyoxometalate (POM) species, K₆ [Hf(α₂ -P₂ W₁₇ O₆₁ )] (1), is reported. In the co-crystal of hen egg-white lysozyme (HEWL) with 1, the catalytically active monomeric species is observed, originating from the dimeric 1:2 POM form, while it is intrinsically unstable under physiological pH conditions. The protein-assisted dissociation of the dimeric POM was rationalized by means of DFT calculations. The dissociation process is unfavorable in bulk water, but becomes favorable in the protein-POM complex due to the low dielectric response at the protein surface. The crystal structure shows that the monomeric form is stabilized by electrostatic and water-mediated hydrogen bonding interactions with the protein. It interacts at three distinct sites, close to the aspartate-containing hydrolysis sites, demonstrating high selectivity towards peptide bonds containing this residue.

Toward Highly Accurate Spin State Energetics in First-Row Transition Metal Complexes: A Combined CASPT2/CC Approach

In previous work on the performance of multiconfigurational second-order perturbation theory (CASPT2) in describing spin state energetics in first-row transition metal systems [ Pierloot et al. J. Chem. Theory Comput. 2017 , 13 , 537 - 553 ], we showed that standard CASPT2 works well for valence correlation but does not describe the metal semicore (3s3p) correlation effects accurately. This failure is partially responsible for the well-known bias toward high-spin states of CASPT2. In this paper, we expand our previous work and show that this bias could be partly removed with a combined CASPT2/CC approach: using high-quality CASPT2 with extensive correlation-consistent basis sets for valence correlation and low-cost CCSD(T) calculations with minimal basis sets for the metal semicore (3s3p) correlation effects. We demonstrate that this approach is efficient by studying the spin state energetics of a series of iron complexes modeling important intermediates in oxidative catalytic processes in chemistry and biochemistry. On the basis of a comparison with bare CCSD(T) results from this and previous work, the average error of the CASPT2/CC approach is estimated at around 2 kcal mol^-1 in favor of high spin states.

The dioxygen adducts of iron and manganese porphyrins: electronic structure and binding energy

The electronic structures of adducts of O₂ and metal porphyrins were thoroughly investigated by highly accurate DMRG-CASPT2.

H2binding to the active site of [NiFe] hydrogenase studied by multiconfigurational and coupled-cluster methods

CCSD(T) and DMRG-CASPT2 calculations show that H₂prefers to bind to Ni rather than to Fe in [NiFe] hydrogenase.

Phosphate Ester Bond Hydrolysis Promoted by Lanthanide-Substituted Keggin-type Polyoxometalates Studied by a Combined Experimental and Density Functional Theory Approach

Hydrolytic cleavage of 4-nitrophenyl phosphate (NPP), a commonly used DNA model substrate, was examined in the presence of series of lanthanide-substituted Keggin-type polyoxometalates (POMs) [Me₂NH₂]₁₁[Ce^III(PW₁₁O₃₉)₂], [Me₂NH₂]₁₀[Ce^IV(PW₁₁O₃₉)₂] (abbreviated as (Ce^IV(PW₁₁)₂), and K₄[EuPW₁₁O₃₉] by means of NMR and luminescence spectroscopies and density functional theory (DFT) calculations. Among the examined complexes, the Ce(IV)-substituted Keggin POM (Ce^IV(PW₁₁)₂) showed the highest reactivity, and its aqueous speciation was fully determined under different conditions of pD, temperature, concentration, and ionic strength by means of ³¹P and ³¹P diffusion-ordered NMR spectroscopy. The cleavage of the phosphoester bond of NPP in the presence of (Ce^IV(PW₁₁)₂) proceeded with an observed rate constant k_obs = (5.31 ± 0.06) × 10^-6 s^-1 at pD 6.4 and 50 °C. The pD dependence of NPP hydrolysis exhibits a bell-shaped profile, with the fastest rate observed at pD 6.4. The formation constant (K_f = 127 M^-1) and catalytic rate constant (k_c = 19.41 × 10^-5 s^-1) for the NPP-Ce(IV)-Keggin POM complex were calculated, and binding between Ce^IV(PW₁₁)₂ and the phosphate group of NPP was also evidenced by the change of the chemical shift of the ³¹P nucleus in NPP upon addition of the POM complex. DFT calculations revealed that binding of NPP to the parent catalyst Ce^IV(PW₁₁)₂ is thermodynamically unlikely. On the contrary, formation of complexes with the monomeric 1:1 species, Ce^IVPW₁₁, is considered to be more favorable, and the most stable complex, [Ce^IVPW₁₁(H₂O)₂(NPP-κO)₂]^7-, was found to involve two NPP ligands coordinated to the Ce^IVcenter of Ce^IVPW₁₁ in the monodentate fashion. The formation of such species is considered to be responsible for the hydrolytic activity of Ce^IV(PW₁₁)₂ toward phosphomonoesters. On the basis of these findings a principle mechanism for the hydrolysis of NPP by the POM is proposed.

The active site of low-temperature methane hydroxylation in iron-containing zeolites

An efficient catalytic process for converting methane into methanol could have far-reaching economic implications. Iron-containing zeolites (microporous aluminosilicate minerals) are noteworthy in this regard, having an outstanding ability to hydroxylate methane rapidly at room temperature to form methanol. Reactivity occurs at an extra-lattice active site called α-Fe(ii), which is activated by nitrous oxide to form the reactive intermediate α-O; however, despite nearly three decades of research, the nature of the active site and the factors determining its exceptional reactivity are unclear. The main difficulty is that the reactive species-α-Fe(ii) and α-O-are challenging to probe spectroscopically: data from bulk techniques such as X-ray absorption spectroscopy and magnetic susceptibility are complicated by contributions from inactive 'spectator' iron. Here we show that a site-selective spectroscopic method regularly used in bioinorganic chemistry can overcome this problem. Magnetic circular dichroism reveals α-Fe(ii) to be a mononuclear, high-spin, square planar Fe(ii) site, while the reactive intermediate, α-O, is a mononuclear, high-spin Fe(iv)=O species, whose exceptional reactivity derives from a constrained coordination geometry enforced by the zeolite lattice. These findings illustrate the value of our approach to exploring active sites in heterogeneous systems. The results also suggest that using matrix constraints to activate metal sites for function-producing what is known in the context of metalloenzymes as an 'entatic' state-might be a useful way to tune the activity of heterogeneous catalysts.

Theoretical Study of the Dissociation Energy of First-Row Metallocenium Ions

The bond dissociation energy of a series of metallocenium ions, i.e., the energy difference of the reaction MCp2(+) → MCp(+) + Cp· (with M = Ti, V, Cr, Mn, Fe, Co, and Ni), was studied by means of multiconfigurational perturbation theory (CASPT2, RASPT2, NEVPT2) and restricted coupled cluster theory (CCSD(T)). From a comparison between the results obtained from these different methods, and a detailed analysis of their treatment of electron correlation effects, a set of MCp(+)-Cp binding energies are proposed with an accuracy of 5 kcal/mol. The computed results are in good agreement with the experimental data measured by threshold photoelectron photoion coincidence (TPEPICO) spectroscopy but disagree with the more recent threshold collision-induced dissociation (TCID) experiments.

Reactivity of Dimeric Tetrazirconium(IV) Wells-Dawson Polyoxometalate toward Dipeptide Hydrolysis Studied by a Combined Experimental and Density Functional Theory Approach

Detailed kinetic studies on the hydrolysis of glycylglycine (Gly-Gly) in the presence of the dimeric tetrazirconium(IV)-substituted Wells-Dawson-type polyoxometalate Na14[Zr4(P2W16O59)2(μ3-O)2(OH)2(H2O)4] · 57H2O (1) were performed by a combination of (1)H, (13)C, and (31)P NMR spectroscopies. The catalyst was shown to be stable under a broad range of reaction conditions. The effect of pD on the hydrolysis of Gly-Gly showed a bell-shaped profile with the fastest hydrolysis observed at pD 7.4. The observed rate constant for the hydrolysis of Gly-Gly at pD 7.4 and 60 °C was 4.67 × 10(-7) s(-1), representing a significant acceleration as compared to the uncatalyzed reaction. (13)C NMR data were indicative for coordination of Gly-Gly to 1 via its amide oxygen and amine nitrogen atoms, resulting in a hydrolytically active complex. Importantly, the effective hydrolysis of a series of Gly-X dipeptides with different X side chain amino acids in the presence of 1 was achieved, and the observed rate constant was shown to be dependent on the volume, chemical nature, and charge of the X amino acid side chain. To give a mechanistic explanation of the observed catalytic hydrolysis of Gly-Gly, a detailed quantum-chemical study was performed. The theoretical results confirmed the nature of the experimentally suggested binding mode in the hydrolytically active complex formed between Gly-Gly and 1. To elucidate the role of 1 in the hydrolytic process, both the uncatalyzed and the polyoxometalate-catalyzed reactions were examined. In the rate-determining step of the uncatalyzed Gly-Gly hydrolysis, a carboxylic oxygen atom abstracts a proton from a solvent water molecule and the nascent OH nucleophile attacks the peptide carbon atom. Analogous general-base activity of the free carboxylic group was found to take place also in the case of polyoxometalate-catalyzed hydrolysis as the main catalytic effect originates from the -C═O···Zr(IV) binding.

Multinuclear diffusion NMR spectroscopy and DFT modeling: a powerful combination for unraveling the mechanism of phosphoester bond hydrolysis catalyzed by metal-substituted polyoxometalates

A detailed reaction mechanism is proposed for the hydrolysis of the phosphoester bonds in the DNA model substrate bis(4-nitrophenyl) phosphate (BNPP) in the presence of the Zr(IV)-substituted Keggin type polyoxometalate (Et2NH2)8[{α-PW11O39Zr(μ-OH)(H2O)}2]⋅7 H2O (ZrK 2:2) at pD 6.4. Low-temperature (31)P DOSY spectra at pD 6.4 gave the first experimental evidence for the presence of ZrK 1:1 in fast equilibrium with ZrK 2:2 in purely aqueous solution. Moreover, theoretical calculations identified the ZrK 1:1 form as the potentially active species in solution. The reaction intermediates involved in the hydrolysis were identified by means of (1)H/(31)P NMR studies, including EXSY and DOSY NMR spectroscopy, which were supported by DFT calculations. This experimental/theoretical approach enabled the determination of the structures of four intermediate species in which the starting compound BNPP, nitrophenyl phosphate (NPP), or the end product phosphate (P) is coordinated to ZrK 1:1. In the proposed reaction mechanism, BNPP initially coordinates to ZrK 1:1 in a monodentate fashion, which results in hydrolysis of the first phosphoester bond in BNPP and formation of NPP. EXSY NMR studies showed that the bidentate complex between NPP and ZrK 1:1 is in equilibrium with monobound and free NPP. Subsequently, hydrolysis of NPP results in P, which is in equilibrium with its monobound form.

A computational study of the glycylserine hydrolysis at physiological pH: a zwitterionic versus anionic mechanism

The hydrolysis of GlySer at physiological pH was investigated by modeling the most feasible reaction mechanisms in aqueous phase at the MP2/6-311+(2df,2p)//SMD-M06/6-311+(2df,2p) level of the theory. To refine the energies of the most relevant transition states along the reaction paths the cluster-continuum concept was adopted. The hydrolytic process could proceed through two competitive mechanisms involving either the zwitterionic or the anionic form of GlySer. The calculations suggest that at physiological pH the actual mechanism is most probably mixed, anionic-zwitterionic. In this reaction scheme the first stage of N→O acyl transfer involves the anionic form whereas the second stage, during which the resultant ester is hydrolyzed, most likely involves the zwitterionic ester form of GlySer. The energy requirement for the first reaction stage is estimated to be slightly lower than for the second one. The calculated activation parameters (e.g. ΔG(#) = 27.8 kcal mol(-1)) for the nucleophilic addition of a water molecule to the ester carbonyl group of the zwitterionic ester are in good agreement with the experimentally determined values at pD 7.4 (ΔG(#) = 28.7 kcal mol(-1)).

Accurate calculations of geometries and singlet–triplet energy differences for active-site models of [NiFe] hydrogenase

The relative stability of singlet and triplet state models of [NiFe] hydrogenase have been studied by advanced quantum-mechanical methods.

Theoretical modelling of photoswitching of hyperpolarisabilities in ruthenium complexes

Static excited-state polarisabilities and hyperpolarisabilities of three Ru^II ammine complexes are computed at the density functional theory (DFT) and several correlated ab initio levels. Most accurate modelling of the low energy electronic absorption spectrum is obtained with the hybrid functionals B3LYP, B3P86 or M06 for the complex [Ru^II(NH₃)₅(MeQ⁺)]³⁺ (MeQ⁺=N-methyl-4,4′-bipyridinium, 3) in acetonitrile. The match with experimental data is less good for [Ru^II(NH₃)₅L]³⁺ (L=N-methylpyrazinium, 2; N-methyl-4-{E,E-4-(4-pyridyl)buta-1,3-dienyl}pyridinium, 4). These calculations confirm that the first dipole- allowed excited state (FDAES) has metal-to-ligand charge-transfer (MLCT) character. Both the solution and gas-phase results obtained for 3 by using B3LYP, B3P86 or M06 are very similar to those from restricted active-space SCF second-order perturbation theory (RASPT2) with a very large basis set and large active space. However, the time-dependent DFT λ_max predictions from the long-range corrected functionals CAM-B3LYP, LC-ωPBE and wB97XB and also the fully ab initio resolution of identity approximate coupled-cluster method (gas-phase only) are less accurate for all three complexes. The ground state (GS) two-state approximation first hyperpolarisability β_2SA for 3 from RASPT2 is very close to that derived experimentally via hyper-Rayleigh scattering, whereas the corresponding DFT-based values are considerably larger. The β responses calculated by using B3LYP, B3P86 or M06 increase markedly as the π-conjugation extends on moving along the series 2→4, for both the GS and FDAES species. All three functionals predict substantial FDAES β enhancements for each complex, increasing with the π-conjugation, up to about sevenfold for 4. Also, the computed second hyperpolarisabilities γ generally increase in the FDAES, but the results vary between the different functionals.

Fourteen-electron ring model and the anomalous magnetic circular dichroism of meso-triarylsubporphyrins

The MCD spectra of meso-triarylsubporphyrins show a sign anomaly which is correlated with the acceptor properties of the aryl substituent. From the spectra, magnetic moments of the excited states are determined. In the context of a simplified orbital model, the sign change is attributed to the quenching of the magnetic moment of the LUMO by acceptor orbitals of the substituent. The actual calculation of this moment presents a major challenge to computational methods. It is shown that wave function techniques based on CASSCF underestimate the covalency effects that are responsible for the quenching. In contrast, a CI method based on DFT orbitals yields excellent results, which fully support the orbital model.

Multiconfigurational Second-Order Perturbation Theory Restricted Active Space (RASPT2) Method for Electronic Excited States: A Benchmark Study

The recently developed second-order perturbation theory restricted active space (RASPT2) method has been benchmarked versus the well-established complete active space (CASPT2) approach. Vertical excitation energies for valence and Rydberg excited states of different groups of organic (polyenes, acenes, heterocycles, azabenzenes, nucleobases, and free base porphin) and inorganic (nickel atom and copper tetrachloride dianion) molecules have been computed at the RASPT2 and multistate (MS) RASPT2 levels using different reference spaces and compared with CASPT2, CCSD, and experimental data in order to set the accuracy of the approach, which extends the applicability of multiconfigurational perturbation theory to much larger and complex systems than previously. Relevant aspects in multiconfigurational excited state quantum chemistry such as the valence-Rydberg mixing problem in organic molecules or the double d-shell effect for first-row transition metals have also been addressed.

Multireference ab initio calculations of g tensors for trinuclear copper clusters in multicopper oxidases

EPR spectroscopy has proven to be an indispensable tool in elucidating the structure of metal sites in proteins. In recent years, experimental EPR data have been complemented by theoretical calculations, which have become a standard tool of many quantum chemical packages. However, there have only been a few attempts to calculate EPR g tensors for exchange-coupled systems with more than two spins. In this work, we present a quantum chemical study of structural, electronic, and magnetic properties of intermediates in the reaction cycle of multicopper oxidases and of their inorganic models. All these systems contain three copper(II) ions bridged by hydroxide or O(2-) anions and their ground states are antiferromagnetically coupled doublets. We demonstrate that only multireference methods, such as CASSCF/CASPT2 or MRCI can yield qualitatively correct results (compared to the experimental values) and consider the accuracy of the calculated EPR g tensors as the current benchmark of quantum chemical methods. By decomposing the calculated g tensors into terms arising from interactions of the ground state with the various excited states, the origin of the zero-field splitting is explained. The results of the study demonstrate that a truly quantitative prediction of the g tensors of exchange-coupled systems is a great challenge to contemporary theory. The predictions strongly depend on small energy differences that are difficult to predict with sufficient accuracy by any quantum chemical method that is applicable to systems of the size of our target systems.