Protein induced singlet–triplet quasidegeneracy in the active site of [NiFe]-hydrogenase

Predicting Kinetics and Dynamics of Spin-Dependent Processes

ConspectusPredicting mechanisms and rates of nonadiabatic spin-dependent processes including photoinduced intersystem crossings, thermally activated spin-forbidden reactions, and spin crossovers in metal centers is a very active field of research. These processes play critical roles in transition-metal-based and metalloenzymatic catalysis, molecular magnets, light-harvesting materials, organic light-emitting diodes, photosensitizers for photodynamic therapy, and many other applications. Therefore, accurate modeling of spin-dependent processes in complex systems and on different time scales is important for many problems in chemistry, biochemistry, and materials sciences.Nonadiabatic statistical theory (NAST) and nonadiabatic molecular dynamics (NAMD) are two complementary approaches to modeling the kinetics and dynamics of spin-dependent processes. NAST predicts the probabilities and rate constants of nonradiative transitions between electronic states with different spin multiplicities using molecular properties at only few critical points on the potential energy surfaces (PESs), including the reactant minimum and the minimum energy crossing point (MECP) between two spin states. This makes it possible to obtain molecular properties for NAST calculations using accurate but often computationally expensive electronic structure methods, which is critical for predicting the rate constants of spin-dependent processes. Alternatively, NAST can be used to study spin-dependent processes in very large complex molecular systems using less computationally expensive electronic structure methods. The nuclear quantum effects, such as zero-point vibrational energy, tunneling, and interference between reaction paths can be easily incorporated. However, the statistical and local nature of NAST makes it more suitable for large systems and slow kinetics. In contrast, NAMD explores entire PESs of interacting electronic states, making it ideal for modeling fast barrierless spin-dependent processes. Because the knowledge of large portions of PESs is often needed, the simulations require a very large number of electronic structure calculations, which limits the NAMD applicability to relatively small molecular systems and ultrafast kinetics.In this Account, we discuss our contribution to the development of the NAST and NAMD approaches for predicting the rates and mechanism of spin-dependent processes. First, we briefly describe our NAST and NAMD implementations. The NAST implementation is an extension of the transition state theory to the processes involving two crossing potential energy surfaces of different spin multiplicities. The NAMD approach includes the trajectory surface hopping (TSH) and ab initio multiple spawning (AIMS) methods. Second, we discuss several applications of NAST and NAMD to model spin-dependent processes in different systems. The NAST applicability to large complex systems is demonstrated by the studies of the spin-forbidden isomerization of the active sites of metal-sulfur proteins. Our implementation of the MECP search algorithm within the fully ab initio fragment molecular orbital method allows applying NAST to systems with thousands of atoms, such as the solvated protein rubredoxin. Applications of NAMD to ultrafast spin-dependent processes are represented by the generalized AIMS simulations utilizing the fast GPU-based TeraChem electronic structure program to gain insight into the complex photoexcited state relaxation in 2-cyclopentenone.

Initial Steps and Thermochemistry of Unimolecular Decomposition of Oxadiazole Energetic Materials: Quantum Chemistry Modeling

In order to resolve the existing discrepancies in the mechanism and key intermediates of oxadiazole thermolysis, the initial decomposition pathways of oxadiazoles have been studied comprehensively using the M062X method for optimization and CBS-QB3 and DLPNO-CCSD(T) methods for energies. The transformation from the furoxan ring to nitro group was suggested as a potential decay channel of furoxan compounds. Results of thermochemistry calculations showed that the preferred decomposition reaction of oxadiazoles is the ring-opening through the cleavage of the O-C or O-N bond. The introduction of the nitro group has little effect on the preferential path of oxadiazole thermal decomposition, but a great impact on the energy barrier. The lowest energy barrier and bond dissociation energy of NO₂ loss of azoles were comprehensively studied based on the quantum chemistry calculations. The initial decay steps of 3,4-dinitrofurazanfuroxan and benzotrifuroxan were also studied to give insights into the mechanism of primary stages of thermal decomposition of oxadiazoles.

Quantum chemical approaches to [NiFe] hydrogenase

The mechanism by which [NiFe] hydrogenase catalyses the oxidation of molecular hydrogen is a significant yet challenging topic in bioinorganic chemistry. With far-reaching applications in renewable energy and carbon mitigation, significant effort has been invested in the study of these complexes. In particular, computational approaches offer a unique perspective on how this enzyme functions at an electronic and atomistic level. In this article, we discuss state-of-the art quantum chemical methods and how they have helped deepen our comprehension of [NiFe] hydrogenase. We outline the key strategies that can be used to compute the (i) geometry, (ii) electronic structure, (iii) thermodynamics and (iv) kinetic properties associated with the enzymatic activity of [NiFe] hydrogenase and other bioinorganic complexes.

Hydrogenase Enzymes and Their Synthetic Models: The Role of Metal Hydrides

Hydrogenase enzymes efficiently process H2 and protons at organometallic FeFe, NiFe, or Fe active sites. Synthetic modeling of the many H2ase states has provided insight into H2ase structure and mechanism, as well as afforded catalysts for the H2 energy vector. Particularly important are hydride-bearing states, with synthetic hydride analogues now known for each hydrogenase class. These hydrides are typically prepared by protonation of low-valent cores. Examples of FeFe and NiFe hydrides derived from H2 have also been prepared. Such chemistry is more developed than mimicry of the redox-inactive monoFe enzyme, although functional models of the latter are now emerging. Advances in physical and theoretical characterization of H2ase enzymes and synthetic models have proven key to the study of hydrides in particular, and will guide modeling efforts toward more robust and active species optimized for practical applications.

Hydride bridge in [NiFe]-hydrogenase observed by nuclear resonance vibrational spectroscopy

The metabolism of many anaerobes relies on [NiFe]-hydrogenases, whose characterization when bound to substrates has proven non-trivial. Presented here is direct evidence for a hydride bridge in the active site of the ⁵⁷Fe-labelled fully reduced Ni-R form of Desulfovibrio vulgaris Miyazaki F [NiFe]-hydrogenase. A unique ‘wagging' mode involving H⁻ motion perpendicular to the Ni(μ-H)⁵⁷Fe plane was studied using ⁵⁷Fe-specific nuclear resonance vibrational spectroscopy and density functional theory (DFT) calculations. On Ni(μ-D)⁵⁷Fe deuteride substitution, this wagging causes a characteristic perturbation of Fe–CO/CN bands. Spectra have been interpreted by comparison with Ni(μ-H/D)⁵⁷Fe enzyme mimics [(dppe)Ni(μ-pdt)(μ-H/D)⁵⁷Fe(CO)₃]⁺ and DFT calculations, which collectively indicate a low-spin Ni(II)(μ-H)Fe(II) core for Ni-R, with H⁻ binding Ni more tightly than Fe. The present methodology is also relevant to characterizing Fe–H moieties in other important natural and synthetic catalysts.

Understanding the catalytic mechanism of redox-active hydrogenases is a key to efficient hydrogen production and consumption. Here, the authors use nuclear resonance vibrational spectroscopy to study [NiFe]-hydrogenase, and observe a bridging hydride structure in an EPR silent intermediate.

Modulation of active site electronic structure by the protein matrix to control [NiFe] hydrogenase reactivity

Control of the reactivity of the nickel center of the [NiFe] hydrogenase and other metalloproteins commonly involves outer coordination sphere ligands that act to modify the geometry and physical properties of the active site metal centers. We carried out a combined set of classical molecular dynamics and quantum/classical mechanics calculations to provide quantitative estimates of how dynamic fluctuations of the active site within the protein matrix modulate the electronic structure at the catalytic center. Specifically we focused on the dynamics of the inner and outer coordination spheres of the cysteinate-bound Ni-Fe cluster in the catalytically active Ni-C state. There are correlated movements of the cysteinate ligands and the surrounding hydrogen-bonding network, which modulate the electron affinity at the active site and the proton affinity of a terminal cysteinate. On the basis of these findings, we hypothesize a coupling between protein dynamics and electron and proton transfer reactions critical to dihydrogen production.

Atomic partitioning of M-H2 bonds in [NiFe] hydrogenase--a test case of concurrent binding

The possibility of simultaneous addition of η(2)-H2 to both the metals (Ni and Fe) in the active site of the as isolated state of the enzyme (Ni-SI) is examined here by an atom-by-atom electronic energy partitioning based on the QTAIM method. Results show that the 4LS state prefers H2 removal than addition. Destabilization of the atomic basins of the thiolate bridges and decrease of the electrophilicity of the Fe and Ni, resulting in poor back donation to the CO ligand, are the bottlenecks that hamper dihydrogen activation simultaneously. The study helps to understand why such states are seldom accessed in the activation of dihydrogen. Moreover, Ni has been found to be the natural choice for the dihydrogen binding.

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.