Quintet Spin States Research Articles

Carbon Dioxide: A Waste Product in the Catalytic Cycle of α-Ketoglutarate Dependent Halogenases Prevents the Formation of Hydroxylated By-Products

We present the first density functional theory study on alpha-ketoglutarate dependent halogenases and focus on the mechanism starting from the iron(IV)-oxo species. The studies show that the high-valent iron(IV)-oxo species reacts with substrates via an initial and rate determining hydrogen abstraction that is characterized by a large kinetic isotope effect (KIE) of 26.7 leading to a radical intermediate. This KIE value is in good agreement with experimental data. The reaction proceeds via two-state reactivity patterns on competing quintet and septet spin state surfaces with close lying hydrogen abstraction barriers. However, the septet spin radical intermediate gives very high barriers for hydroxylation and chlorination whereas the barriers on the quintet spin state surface are much lower. The calculations give extra information regarding the nature of the intermediates and a prediction of a new low-energy mechanism starting from the radical intermediate, whereby a waste product from an earlier step in the catalytic cycle (CO(2)) is recycled and takes the hydroxyl radical away to form bicarbonate via an OH trapping mechanism. As a consequence, this mechanism prevents the occurrence of hydroxylated byproduct and gives a rationale for the sole observance of halogenated products. By contrast, a direct halogenation reaction cannot compete with hydroxylation due to higher reaction barriers. Our findings support experimental work in the field and give a rationale for the lack of hydroxylation products in alpha-ketoglutarate dependent halogenases.

The Effect and Influence of cis-Ligands on the Electronic and Oxidizing Properties of Nonheme Oxoiron Biomimetics. A Density Functional Study

Density functional theory studies on the nature of the cis effect and cis influence of ligands on oxoiron nonheme complexes have been performed. A detailed analysis of the electronic and oxidizing properties of [Fe(IV)O(TPA)L](+) with L = F(-), Cl(-), and Br(-) and TPA = tris-(2-pyridylmethyl)amine are presented and compared with [Fe(IV)O(TPA)NCCH(3)](2+). The calculations show that the electronic cis effect is determined by favorable orbital overlap between first-row elements with the metal, which are missing between the metal and second- and third-row elements. As a consequence, the metal 3d block is split into a one-below-two set of orbitals with L = Cl(-) and Br(-), and the HOMO/LUMO energy gap is widened with respect to the system with L = F(-). However, this larger HOMO/LUMO gap does not lead to large differences in electron affinities of the complexes. Moreover, a quantum mechanical analysis of the binding of the ligand shows that it is built up from a large electric field effect of the ligand on the oxoiron species and a much smaller quantum mechanical effect due to orbital overlap. These contributions are of similar strength for the three tested halogen cis ligands and result in similar reactivity patterns with substrates. The calculations show that [Fe(IV)O(TPA)L](+) with L = F(-), Cl(-), and Br(-) have closely lying triplet and quintet spin states, but only the quintet spin state is reactive with substrates. Therefore, the efficiency of the oxidant will be determined by the triplet-quintet spin state crossing of the reaction. The reaction of styrene with a doubly charged reactant, that is, [Fe(V)O(TPA)L](2+) with L = F(-), Cl(-), and Br(-) or [Fe(V)O(TPA)NCCH(3)](3+), leads to an initial electron transfer from the substrate to the metal followed by a highly exothermic epoxidation mechanism. These reactivity differences are mainly determined by the overall charge of the system rather than the nature of the cis ligand.

High resolution electron paramagnetic resonance spectroscopy of quintet pyridyl-2,6-dinitrene in solid argon: Magnetic properties and molecular structure

The high resolution X-band electron para magnetic resonance (EPR) spectrum of quintet pyridyl-2,6-dinitrene was recorded after the photolysis of 4-amino-2,6-diazido-3,5-dichloropyridine in solid argon matrix at 15 K. This spectrum represents a new type of powder EPR spectra that are characteristic for quintet spin states with zero-field splitting parameters |E(q)/D(q)| approximately 1/4. All EPR lines of the quintet dinitrene were unambiguously assigned based on the eigenfield calculations of the Zeeman energy levels and angular dependencies of resonance magnetic fields. Owing to the high resolution of the experimental EPR spectrum, zero-field splitting parameters of the quintet dinitrene were determined with a high accuracy: D(q)=0.2100+/-0.0005 cm(-1) and E(q)=-0.0560+/-0.0002 cm(-1). These parameters provide correct information regarding the molecular angle Theta and distance r between two triplet sites in the molecule of quintet dinitrene. The measured molecular angle Theta=114.2 degrees+/-0.2 degrees is in excellent agreement with results of the density functional theory calculations. The analysis of the magnetic parameters shows that the spin population on the nitrene units in the quintet dinitrene is greater than that on the nitrene unit in the triplet nitrene.

Photochemical generation of triplet—triplet nitrene pairs in aromatic diazide crystals

The molecular and crystal structures of 4-amino-2,6-diazido-3,5-dichloropyridine and 6-amino-2,4-diazido-1,3,5-triazine, as well as the paramagnetic photolysis products of their crystals at 77 K, were studied using X-ray diffraction analysis and ESR spectroscopy. Triplet nitrenes generated during the photolysis of diazidopyridine form triplet—triplet nitrene pairs, whose ESR spectrum corresponds to the quintet spin state. The high-spin state (S = 2) results from the exchange interaction between two triplet molecules with the zero-field splitting parameters |D| = 1.0280 cm−1 and |E| = 0.0038 cm−1 and the γ angle between two C—N nitrene bonds equal to 133°. This angle is close to an angle of 136.2° between the C-N bonds of two adjacent molecules in the crystal structure. No formation of the triplet—triplet nitrene pairs is observed during the photolysis of crystalline diazidotriazine, whose molecules lie in the parallel planes.

An N-substituted aza[14]metacyclophane tetracation: a spin-quintet tetraradical with four para-phenylenediamine-based semi-quinone moieties

A fully dianisylaminophenylated aza[1(4)]metacyclophane has been synthesized that exhibits four reversible two-electron oxidation processes, and its dicationic and tetracationic species have been found to be in spin-triplet and spin-quintet states.

EPR spectroscopy of quintet 4-amino-3,5-dichloropyridine-2,6-diyldinitrene isolated in solid argon

An EPR spectrum of solid Ar isolated quintet 4-amino-3,5-dichloropyridin-2,6-diyldinitrene that formed by the photolysis of 4-amino-2,6-diazido-3,5-dichloropyridine at 15 K was recorded. Using computer simulation based on numerical diagonalization of the quintet spin Hamiltonian matrices, it was established that this EPR spectrum corresponds to a quintet spin state with the magnetic parameters g = 2.0023, |D q| = 0.2100 cm−1, and |E q| = 0.0560 cm−1. Owing to high resolution of the experimental spectrum, the zero-field splitting parameters of the quintet intermediate were determined to an accuracy of at least 5·10−4 cm−1. Calculations of the fine-structure energy levels in external magnetic field and the dependences of the EPR signal positions and intensities of the quintet dinitrene on the direction of external magnetic field were performed for the first time. This allowed unambiguous assignment of all EPR lines of quintet molecules having both in-principal-axis and off-principal-axis orientations.

The Mechanism of Cysteine Oxygenation by Cysteine Dioxygenase Enzymes

We present here the first density functional theoretic study into the mechanism of cysteine dioxygenation by a model of cysteine dioxygenase enzymes. A large active site model containing the ligands bound to iron plus amino acid residues that are involved in hydrogen bonding interactions with the substrate is used. The reaction takes place via multi-state reactivity patterns on competing singlet, triplet, and quintet spin states, whereby the latter is the ground state in most complexes. Several new intermediates have been predicted, which have not been anticipated before. The dioxygen-bound complex is in a singlet spin ground state, and a state crossing to the quintet spin state leads to an FeOOS ring structure that splits into a cysteinyloxide radical that reorients and abstracts an electron from the iron center. In the final step, the oxoiron donates the oxygen atom to the substrate to produce cysteine sulfinic acid in a highly exothermic process. The rate-determining step is the initial step in the reaction mechanism on the quintet spin state surface.

Electronic Structure of Bispidine Iron(IV) Oxo Complexes

The electronic structure, based on DFT calculations, of a range of FeIV=O complexes with two tetra- (L1 and L2) and two isomeric pentadentate bispidine ligands (L3 and L4) is discussed with special emphasis on the relative stability of the two possible spin states (S = 1, triplet, intermediate-spin, and S = 2, quintet, high-spin; bispidines are very rigid diazaadamantane-derived 3,7-diazabicyclo[3.3.1]nonane ligands with two tertiary amine and two or three pyridine donors, leading to cis-octahedral [(X)(L)FeIV=O]2+ complexes, where X = NCCH3, OH2, OH-, and pyridine, and where X = pyridine is tethered to the bispidine backbone in L3, L4). The two main structural effects are a strong trans influence, exerted by the oxo group in both the triplet and the quintet spin states, and a Jahn-Teller-type distortion in the plane perpendicular to the oxo group in the quintet state. Due to the ligand architecture the two sites for substrate coordination in complexes with the tetradentate ligands L1 and L2 are electronically very different, and with the pentadentate ligands L3 and L4, a single isomer is enforced in each case. Because of the rigidity of the bispidine ligands and the orientation of the "Jahn-Teller axis", which is controlled by the sixth donor X, the Jahn-Teller-type distortion in the high-spin state of the two isomers is quite different. It is shown how this can be used as a design principle to tune the relative stability of the two spin states.

Combined Experimental and Theoretical Study on Aromatic Hydroxylation by Mononuclear Nonheme Iron(IV)−Oxo Complexes

The hydroxylation of aromatic compounds by mononuclear nonheme iron(IV)-oxo complexes, [FeIV(Bn-tpen)(O)]2+ (Bn-tpen=N-benzyl-N,N',N'-tris(2-pyridylmethyl)ethane-1,2-diamine) and [FeIV(N4Py)(O)]2+ (N4Py=N,N-bis(2-pyridylmethyl)-N-bis(2-pyridyl)methylamine), has been investigated by a combined experimental and theoretical approach. In the experimental work, we have performed kinetic studies of the oxidation of anthracene with nonheme iron(IV)-oxo complexes generated in situ, thereby determining kinetic and thermodynamic parameters, a Hammett rho value, and a kinetic isotope effect (KIE) value. A large negative Hammett rho value of -3.9 and an inverse KIE value of 0.9 indicate that the iron-oxo group attacks the aromatic ring via an electrophilic pathway. By carrying out isotope labeling experiments, the oxygen in oxygenated products was found to derive from the nonheme iron(IV)-oxo species. In the theoretical work, we have conducted density functional theory (DFT) calculations on the hydroxylation of benzene by [FeIV(N4Py)(O)]2+. The calculations show that the reaction proceeds via two-state reactivity patterns on competing triplet and quintet spin states via an initial rate determining electrophilic substitution step. In analogy to heme iron(IV)-oxo catalysts, the ligand is noninnocent and actively participates in the reaction mechanism by reshuttling a proton from the ipso position to the oxo group. Calculated kinetic isotope effects of C6H6 versus C6D6 confirm an inverse isotope effect for the electrophilic substitution pathway. Based on the experimental and theoretical results, we have concluded that the aromatic ring oxidation by mononuclear nonheme iron(IV)-oxo complexes does not occur via a hydrogen atom abstraction mechanism but involves an initial electrophilic attack on the pi-system of the aromatic ring to produce a tetrahedral radical or cationic sigma-complex.

What Factors Influence the Ratio of CH Hydroxylation versus CC Epoxidation by a Nonheme Cytochrome P450 Biomimetic?

Density functional calculations on a nonheme biomimetic (Fe=O(TMCS+) have been performed and its catalytic properties versus propene investigated. Our studies show that this catalyst is able to chemoselectively hydroxylate C=H bonds even in the presence of C=C double bonds. This phenomenon has been analyzed and found to occur due to Pauli repusions between protons on the TMCS ligand with protons attached to the approaching substrate. The geometries of the rate determining transition states indicate that the steric hindrance is larger in the epoxidation transition states than in the hydroxylation ones with much shorter distances; hence the hydroxylation pathway is favored over the epoxidation. Although, the reactant experiences close lying triplet and quintet spin states, the dominant reaction mechanism takes place on the quintet spin state surface; i.e., Fe=O(TMCS)+ reacts via single-state reactivity. Our calculations show that this spin state selectivity is the result of geometric orientation of the transition state structures, whereby the triplet ones are destabilized by electrostatic repulsions between the substrate and the ligand while the quintet spin transition states are aligned along the ideal axis. The reactivity patterns and geometries are compared with oxoiron species of dioxygenase and monoxygenase enzymes. Thus, Fe=O(TMCS)+ shows some similarities with P450 enzyme reactivity: it chemoselectively hydroxylates C=H bonds even in the presence of a C=C double bond and therefore is an acceptable P450 biomimetic. However, the absolute barriers of substrate oxidation by Fe=O(TMCS)+ are higher than the ones obtained with heme enzymes, but the chemoselectivity is lesser affected by external perturbations such as hydrogen bonding of a methanol molecule toward the thiolate sulfur or a dielectric constant. This is the first oxoiron complex whereby we calculated a chemoselective hydroxylation over epoxidation in the gas phase.

Propene Activation by the Oxo-Iron Active Species of Taurine/α-Ketoglutarate Dioxygenase (TauD) Enzyme. How Does the Catalysis Compare to Heme-Enzymes?

Density functional calculations on the oxygenation reaction of propene by a model for taurine/alpha-ketoglutarate dioxygenase (TauD) enzyme are presented. The oxo-iron active species of TauD is shown to be a powerful and aggressive oxidant, which is able to hydroxylate C-H bonds and epoxidize C=C bonds with low barriers. In the case of propene oxygenation, the hydroxylation and epoxidation mechanisms are competitive on a dominant quintet spin state surface. We have compared the mechanism and thermodynamics of TauD with oxo-iron heme catalysts, such as the cytochromes P450, and found some critical differences. The TauD model is found to be much more reactive toward oxygenation of substrates than oxo-iron complexes in a heme environment with much lower reaction barriers. We have analyzed this and assigned this to the strength of the O-H bond formed after hydrogen abstraction from a substrate, which is at least 10 kcal mol(-)(1) stronger in five-coordinated oxo-iron nonheme complexes than in six-coordinated oxo-iron heme complexes. Since, the metal in TauD enzymes is five-coordinated, whereas in heme-enzymes it is six-coordinated there are some critical differences in the valence molecular orbitals. Thus, in oxo-iron heme catalysts one of the antibonding pi orbitals is replaced by a low-lying nonbonding delta orbital resulting in a lower overall spin state. Moreover, heme-enzymes have an extra oxidation equivalent located on the heme, which is missing in non-heme oxo-iron catalysts. As a result, the oxo-iron species of TauD reacts via single-state reactivity on a dominant quintet spin state surface, whereas oxo-iron heme catalysts react via two-state reactivity on competing doublet and quartet spin states.

First observation of the hyperfine structure of an excited quintet state in liquid solution

After pulsed photoexcitation of a new fullerene-linked bisnitroxide, a well resolved transient EPR spectrum is detected which is assigned to an excited quintet spin state generated by spin coupling of the nitroxides and the fullerene excited triplet.

Combined Quantum Mechanical/Molecular Mechanical Study on the Pentacoordinated Ferric and Ferrous Cytochrome P450cam Complexes

The pentacoordinated ferric and ferrous cytochrome P450(cam) complexes have been investigated by combined quantum mechanical/molecular mechanical (QM/MM) calculations in the presence of a protein/solvent environment and by QM calculations on the isolated QM regions with use of density functional theory. The B3LYP functional has been found more reliable than the BLYP and BHLYP functionals for estimating the relative state energies. The B3LYP/CHARMM calculations with an all-electron basis set for iron give high-spin ground states for the title complexes, in agreement with experiment. The comparison of the B3LYP/CHARMM results of the entire protein system with the B3LYP calculations on the naked QM regions shows that the amount of stabilization by the protein environment is largest for the intermediate-spin states, followed by the high-spin states of the complexes. The calculation of Mössbauer parameters in the presence of the enzyme environment confirms the double occupation of the d(xz) orbital in the quintet spin state of the ferrous complex, consistent with the computed QM/MM energies in the enzyme environment, while the d(x)2(-)(y)2 orbital is doubly occupied in the gas-phase quintet state.

Thermochemistry of spin-crossover Fe(II)complexes calculated with density functional methods

The optimized geometries, vibrational frequencies and thermodynamic functions of the singlet and quintet spin states of three ferrous complexes: Fe(1,10-phenanthroline) 2(NCS) 2, [Fe(2-picolylamine) 3] 2+ and [Fe(bis(1,4,7-triazacyclononane) 2] 2+ were calculated by the BP86/6-31G* density functional method. The characteristics of the spin-crossover transitions in solution were simulated by the mean-field approximation. The model uses the ideal gas enthalpy and entropy functions based on DFT results and two parameters: Λ, which scales the electronic energy difference Δ E elec and mimics solute–solvent interactions and Γ, which describes interactions between solute molecules. Since the fitted quantities are unique functions of the two parameters their experimental values could be exactly recovered.

Possible Photoinduced Spin Transitions in Bis(phenylmethylenyl)[2.2]paracyclophanes. A Spin−Orbit Coupling Study

A possible mechanism for the spin transitions in various stacking conformations of bis(phenylmethylenyl)[2.2]paracyclophanes, which have close-lying lowest singlet, triplet, and quintet spin states, is theoretically investigated by using diphenylcarbene dimers as models. Spin−orbit coupling (SOC) matrix elements, which play an essential role in the spin transition phenomena, are calculated with the effective one-electron spin−orbit Hamiltonian. The SOC between the first excited singlet state and the first excited triplet state and that between the first excited triplet state and the lowest quintet state are strong. The SOC between the first excited quintet state and the first excited triplet state and that between the first excited triplet state and the lowest singlet state are also strong. These results demonstrate that the spin conversion between the low-spin singlet state and the high-spin quintet state can occur via the first excited intermediate-spin triplet state. We propose that possible photoinduced spin-crossover phenomena can be observed in these organic molecular systems.

Cis-Stilbene and (1 alpha,2 beta,3 alpha)-(2-ethenyl-3-methoxycyclopropyl)benzene as mechanistic probes in the Mn(III)(salen)-catalyzed epoxidation: influence of the oxygen source and the counterion on the diastereoselectivity of the competitive concerted and radical-type oxygen transfer.

cis-Stilbene (1) has been epoxidized by a set of diverse oxygen donors [OxD], catalyzed by the Mn(III)(salen)X complexes 3 (X = Cl, PF(6)), to afford a mixture of cis- and trans-epoxides 2. The cis/trans ratios range from 29:71 (extensive isomerization) to 92:8, which depends both on the oxygen source [OxD] and on the counterion X of the catalyst. When (1 alpha,2 beta,3 alpha)-(2-ethenyl-3-methoxycyclopropyl)-benzene (4) is used as substrate, a mechanistic probe which differentiates between radical and cationic intermediates, no cationic ring-opening products are found in this epoxidation reaction; thus, isomerized epoxide product arises from intermediary radicals. The dependence of the diastereoselectivity on the oxygen source is rationalized in terms of a bifurcation step in the catalytic cycle, in which concerted Lewis-acid-activated oxygen transfer competes with stepwise epoxidation by the established Mn(V)(oxo) species. The experimental counterion effect is attributed to the computationally assessed ligand-dependent reaction profiles and stereoselectivities of the singlet, triplet, and quintet spin states available to the manganese species.

π-Topology and Spin Alignment Utilizing the Excited Molecular Field: Observation of the Excited High-Spin Quartet (S = 3/2) and Quintet (S = 2) States on Purely Organic π-Conjugated Spin Systems

As a model system for the photoinduced/photoswitched spin alignment in a purely organic pi-conjugated spin system, 9-[4-(4,4,5,5-tetramethyl-1-yloxyimidazolin-2-yl)phenyl]anthracene (1a), 9-[3-(4,4,5,5-tetramethyl-1-yloxyimidazolin-2-yl)phenyl]anthracene (1b), 9,10-bis[4-(4,4,5,5-tetramethyl-1-yloxyimidazolin-2-yl)phenyl]anthracene (2a), and 9,10-bis[3-(4,4,5,5-tetramethyl-1-yloxyimidazolin-2-yl)phenyl]anthracene (2b) were designed and synthesized. In these spin systems, 9-phenylanthracene and 9,10-diphenylanthracene were chosen as photo spin couplers and iminonitroxide was chosen as a dangling stable radical. Time-resolved electron spin resonance (TRESR) spectra of the first excited states with resolved fine-structure splittings were observed for 1a and 2a in an EPA or a 2-MTHF rigid glass matrix. Using the spectral simulation based on the eigenfield method, the observed TRESR spectra for 1a and 2a were unambiguously assigned as an excited quartet (S = 3/2) spin state (Q) and an excited quintet (S = 2) spin state (Qu), respectively. The g value and fine-structure splitting for the quartet state of 1a were determined to be g(Q) = 2.0043, D(Q) = 0.0235 cm(-1), and E(Q) = 0.0 cm(-1). The relative populations (polarization) of each M(S)() sublevel in Q were determined to be P(+1/2') = P(-1/2') = 0.5 and P(+3/2') = P(-3/2') = 0.0 with an increasing order of energy in zero magnetic field. The spin Hamiltonian parameters for Qu are g = 2.0043, D = 0.0130 cm(-1), and E = 0.0 cm(-1), and the relative populations in Qu were determined to be P(0') = 0.30, P(-1') = P(+1') = 0.35 and P(-2') = P(+2') = 0.0. These are the first observations of a photoexcited quartet and a quintet high-spin state in pi-conjugated triplet-radical pair systems. In contrast high-spin excited states were not observed for 1b and 2b, the pi-topological isomers of 1a and 2a, showing the role of pi-topology in the spin alignment of the excited states. Since a weak antiferromagnetic exchange interaction was observed in the ground state of 2a, the clear detection of the excited quintet high-spin state shows that the effective exchange coupling between the two dangling radicals through the diphenylanthracene spin coupler has been changed from antiferromagnetic to ferromagnetic upon photoexcitation. Thus, a photoinduced spin alignment utilizing the excited triplet molecular field was realized for the first time in the purely organic pi-conjugated spin system. Furthermore, the mechanism for the generation of dynamic electron spin polarization was investigated for the observed quartet and quintet states, and a plausible mechanism of the enhanced selective intersystem crossing was proposed. Ab initio molecular orbital calculations based on density functional theory were carried out to determine the electronic structures of the excited high-spin states and to understand the mechanism of the spin alignment utilizing the excited molecular field. The role of the spin delocalization and the spin polarization mechanisms were revealed on the photoexcited state.

EPR spectral determination of electronic substituent effects on the D values of hydrocarbon polyradicals (Quintet and septet spin states) composed of localized 1,3-cyclopentanediyl spin-carrying units linked by 1,3-Di- and 1,3,5-trimethylenebenzene ferromagnetic couplers

The parent and p-nitrophenyl-substituted diradicals D-3a,b (triplets), tetraradicals T-3a,b (quintets), and hexaradicals H-3a,b (septets) were photochemically generated in matrix-isolated form (toluene, 77 K) by successive denitrogenation of the trisazoalkanes 3a,b and EPR spectrally characterized. In these high-spin polyradicals the spin-spin interaction within the localized spin-carrying 1,3-cyclopentanediyl diradical unit is much stronger than within the cross-conjugated ferromagnetic coupling unit. Accordingly, a change of the electronic properties in the cyclopentanediyl unit affects decisively the D value of the whole polyradical. Therefore, the spin-accepting p-nitro group reduces the D value of the tetra- and hexaradical in the same amount as that of the diradical. Thus, irrespective of the spin multiplicity, the substituent stabilizes electronically the triplet (D-3a,b), quintet (T-3a,b), and septet (H-3a,b) species with equal efficacy.

Metal-Metal Bonding in d(1)d(1) and d(2)d(2) Bioctahedral Dimer Systems: A Density Functional Study of Face-Shared M(2)X(9)(3)(-) (M = Ti, Zr, Hf, V, Nb, Ta) Complexes.

Density functional theory is used to investigate the electronic and geometric structures and periodic trends in metal-metal bonding of d(1)d(1) and d(2)d(2) face-shared M(2)X(9)(3)(-) dimers of Ti, Zr, Hf (d(1)d(1)) and V, Nb, Ta (d(2)d(2)). For these systems three distinct coupling modes can be recognized, depending on the occupation of the trigonal t(2g)(a(1) + e) single-ion orbitals, which determine the ground-state geometry and extent of metal-metal bonding. For Ti(2)Cl(9)(3)(-), the [a(1) x a(1)] broken-symmetry optimized structure, corresponding to significant delocalization of the metal-based sigma electrons, nicely rationalizes the strong antiferromagnetic coupling reported for Cs(3)Ti(2)Cl(9). The ground-state geometries for Zr(2)Cl(9)(3)(-) and Hf(2)Cl(9)(3)(-) correspond to complete delocalization of the metal-based electrons in a metal-metal sigma bond. For V(2)Cl(9)(3)(-), the global minimum is found to be the ferromagnetic [a(1)e x e(2)] spin-quintet state giving rise to a long V-V separation, consistent with the known structure and reported weak ferromagnetic behavior of Cs(3)V(2)Cl(9). For Nb(2)X(9)(3)(-) (X = Cl, Br, I) and Ta(2)Cl(9)(3)(-), the [a(1)e x a(1)e] spin-triplet state, where complete delocalization of the sigma and delta(pi) electrons occur in a metal-metal double bond, is found to be the global minimum and consequently relatively short internuclear distances result, again, in good agreement with experiment. The periodic trends in metal-metal bonding in these and the isovalent d(3)d(3) complexes can be rationalized in terms of the energetic contributions of orbital overlap (DeltaE(ovlp)) and spin polarization (DeltaE(spe)), the difference DeltaE(spe) - DeltaE(ovlp) determining the tendency of the metal-based electrons to delocalize in the dimer. For d(1)d(1) systems, DeltaE(ovlp) is always greater than DeltaE(spe) and therefore delocalized ground states result for all complexes of the titanium triad. Across the first transition series, the dramatic increase in DeltaE(spe) dominates DeltaE(ovlp) and therefore V(2)Cl(9)(3)(-) and Cr(2)Cl(9)(3)(-) have localized ground states. For the second and third transition series, the much larger DeltaE(ovlp) term ensures that all these complexes remain delocalized.

Ab initio calculations on iron-porphyrin model systems for intermediates in the oxidative cycle of cytochrome P450s.

Geometry optimizations for several spin states of the iron(III)-S-methyl- porphyrin complex, the iron (III)-oxo-S-methyl-porphyrin complex and the respective anions were performed in order to examine models for intermediates in the oxidative cycle of cytochrome P450. The aim of this study was to obtain insights into the ground states of the intermediates of this catalytic cycle and to use the ab initio calculated geometries and charge distributions to suggest better and more realistic parameters for forcefields which are generally used for modeling P450s. The results indicate that the ground states of both the iron(III)-S-methyl-porphyrin complex and the iron(III)-oxo-S-methyl-porphyrin complex are sextet spin states (high spin). The ground states of the anions of both complexes are probably quintet spin states. The fact that experimentally a shift from low spin to high spin is observed upon binding of the substrate suggests that the ab initio calculations for the iron(III)-S-methyl-porphyrin complex in vacuum give a correct representation of the (hydrophobic) substrate-bound state of the active site of P450. The ab initio geometries of the iron-porphyrin complexes are very similar to the experimentally observed geometries, except for the longer iron-sulfur bond in ab initio calculations, which is probably caused by the omission of polarization functions on the sulfur atom during the geometry optimization. The charge distribution in all ab initio calculated complexes can be described by a series of concentric rings of alternating charge, thus allowing a relatively large positive charge on the iron atom. The commonly used forcefields generally underestimate the charge differences between the iron atom and the different parts of the porphyrin moiety or ignore the charges completely. Although forcefield calculations can reproduce the experimental geometry of iron-porphyrin moieties, extension of the forcefields with charges obtained from ab initio calculations should give a better description of the heme moiety in protein modeling and docking experiments.