Daniel Bím

Mono- and binuclear non-heme iron chemistry from a theoretical perspective.

In this minireview, we provide an account of the current state-of-the-art developments in the area of mono- and binuclear non-heme enzymes (NHFe and NHFe2) and the smaller NHFe(2) synthetic models, mostly from a theoretical and computational perspective. The sheer complexity, and at the same time the beauty, of the NHFe(2) world represents a challenge for experimental as well as theoretical methods. We emphasize that the concerted progress on both theoretical and experimental side is a conditio sine qua non for future understanding, exploration and utilization of the NHFe(2) systems. After briefly discussing the current challenges and advances in the computational methodology, we review the recent spectroscopic and computational studies of NHFe(2) enzymatic and inorganic systems and highlight the correlations between various experimental data (spectroscopic, kinetic, thermodynamic, electrochemical) and computations. Throughout, we attempt to keep in mind the most fascinating and attractive phenomenon in the NHFe(2) chemistry, which is the fact that despite the strong oxidative power of many reactive intermediates, the NHFe(2) enzymes perform catalysis with high selectivity. We conclude with our personal viewpoint and hope that further developments in quantum chemistry and especially in the field of multireference wave function methods are needed to have a solid theoretical basis for the NHFe(2) studies, mostly by providing benchmarking and calibration of the computationally efficient and easy-to-use DFT methods.

Accurate Prediction of One-Electron Reduction Potentials in Aqueous Solution by Variable-Temperature H-Atom Addition/Abstraction Methodology

A robust and efficient theoretical approach for calculation of the reduction potentials of charged species in aqueous solution is presented. Within this approach, the reduction potential of a charged complex (with a charge |n| ≥ 2) is probed by means of the reduction potential of its neutralized (protonated/deprotonated) cognate, employing one or several H-atom addition/abstraction thermodynamic cycles. This includes a separation of one-electron reduction from protonation/deprotonation through the temperature dependence. The accuracy of the method has been assessed for the set of 15 transition-metal complexes that are considered as highly challenging systems for computational electrochemistry. Unlike the standard computational protocol(s), the presented approach yields results that are in excellent agreement with experimental electrochemical data. Last but not least, the applicability and limitations of the approach are thoroughly discussed.

Cytosine radical cation: a gas‐phase study combining IRMPD spectroscopy, UV‐PD spectroscopy, ion‐molecule reactions, and theoretical calculations

The radical cation of cytosine (Cyt.+ ) is generated by dissociative oxidation from a ternary CuII complex in the gas phase. The radical cation is characterized by infrared multiple photon dissociation (IRMPD) spectroscopy in the fingerprint region, UV/Vis photodissociation (UVPD) spectroscopy, ion-molecule reactions, and theoretical calculations (density functional theory and ab initio). The experimental IRMPD spectrum features diagnostic bands for two enol-amino and two keto-amino tautomers of Cyt.+ that are calculated to be among the lowest energy isomers, in agreement with a previous study. Although the UVPD action spectrum can also be matched to a combination of the four lowest energy tautomers, the presence of a nonclassical distonic radical cation cannot be ruled out. Its formation is, however, unlikely due to the high energy of this isomer and the respective ternary CuII complex. Gas-phase ion-molecule reactions showed that Cyt.+ undergoes hydrogen-atom abstraction from 1-propanethiol, radical recombination reactions with nitric oxide, and electron transfer from dimethyl disulfide.

Copper(II) and Zinc(II) Complexes of Conformationally Constrained Polyazamacrocycles as Efficient Catalysts for RNA Model Substrate Cleavage in Aqueous Solution at Physiological pH

As part of a quest for efficient artificial catalysts of RNA phosphodiester bond cleavage, conformationally constrained mono- and bis-polyazamacrocycles in which tri- or tetraazaalkane chains link the ortho positions of a benzene ring were synthesized. The catalytic activities of mono- and dinuclear copper(II) and zinc(II) complexes of these polyazamacrocycles towards cleavage of the P-O bond in 2-hydroxypropyl-4-nitrophenylphosphate (HPNP) in aqueous solution at pH 7 have been determined. Only the complexes of the ligands incorporating three nitrogen atoms in a macrocycle proved to be capable of efficiently catalyzing HPNP transesterification. The dinuclear complexes were found to be approximately twice as efficient as their mononuclear counterparts, and exhibited Michaelis-Menten saturation kinetics with calculated rate constants of kcat ≈10(-4)  s(-1) . By means of quantum chemical calculations (DFT/COSMO-RS), several plausible reaction coordinates were described. By correlating the calculated barriers with the experimental kinetic data, two possible reaction scenarios were revealed, with activation free energies of 20-25 kcal mol(-1) .

Radical Reactions Affecting Polar Groups in Threonine Peptide Ions

Peptide cation-radicals containing the threonine residue undergo radical-induced dissociations upon collisional activation and photon absorption in the 210-400 nm range. Peptide cation-radicals containing a radical defect at the N-terminal residue, [•Ala-Thr-Ala-Arg+H]+, were generated by electron transfer dissociation (ETD) of peptide dications and characterized by UV-vis photodissociation action spectroscopy combined with time-dependent density functional theory (TD-DFT) calculations of absorption spectra, including thermal vibronic band broadening. The action spectrum of [•Ala-Thr-Ala-Arg+H]+ ions was indicative of the canonical structure of an N-terminally deaminated radical whereas isomeric structures differing in the position of the radical defect and amide bond geometry were excluded. This indicated that exothermic electron transfer to threonine peptide ions did not induce radical isomerizations in the fragment cation-radicals. Several isomeric structures, ion-molecule complexes, and transition states for isomerizations and dissociations were generated and analyzed by DFT and Møller-Plesset perturbational ab initio calculations to aid interpretation of the major dissociations by loss of water, hydroxyl radical, C3H6NO•, C3H7NO, and backbone cleavages. Born-Oppenheimer molecular dynamics (BOMD) in combination with DFT gradient geometry optimizations and intrinsic reaction coordinate analysis were used to search for low-energy cation-radical conformers and transition states. BOMD was also employed to analyze the reaction trajectory for loss of water from ion-molecule complexes.

Toward Accurate Conformational Energies of Smaller Peptides and Medium-Sized Macrocycles: MPCONF196 Benchmark Energy Data Set

A carefully selected set of acyclic and cyclic model peptides and several other macrocycles, comprising 13 compounds in total, has been used to calibrate the accuracy of the DFT(-D3) method for conformational energies, employing BP86, PBE0, PBE, B3LYP, BLYP, TPSS, TPSSh, M06-2X, B97-D, OLYP, revPBE, M06-L, SCAN, revTPSS, BH-LYP, and ωB97X-D3 functionals. Both high- and low-energy conformers, 15 or 16 for each compound adding to 196 in total, denoted as the MPCONF196 data set, were included, and the reference values were obtained by the composite protocol, yielding the CCSD(T)/CBS extrapolated energies or their DLPNO-CCSD(T)/CBS equivalents in the case of larger systems. The latter was shown to be in near-quantitative (∼0.10-0.15 kcal·mol-1) agreement with the canonical CCSD(T), provided the TightPNO setting is used, and, therefore, can be used as the reference for larger systems (likely up to 150-200 atoms) for the problem studied here. At the same time, it was found that many D3-corrected DFT functionals provide results of ∼1 kcal·mol-1 accuracy, which we consider as quite encouraging. This result implies that DFT-D3 methods can be, for example, safely used in efficient conformational sampling algorithms. Specifically, the DFT-D3/DZVP-DFT level of calculation seems to be the best trade-off between computational cost and accuracy. Based on the calculated data, we have not found any cheaper variant for the treatment of conformational energies, since the semiempirical methods (including DFTB) provide results of inferior accuracy (errors of 3-5 kcal·mol-1).

Beyond the classical thermodynamic contributions to hydrogen atom abstraction reactivity

Significance Hydrogen atom abstraction reactivity is a key property of many important biological and synthetic compounds that depends on their proton-coupled reduction potentials. These potentials quantify the ability of species to acquire an electron and proton pair. Intuitively, a species with a higher proton-coupled reduction potential abstracts hydrogen atoms more easily, which translates into a lower reaction barrier. Beyond this classical thermodynamic effect on reactivity, we discovered a significant contribution arising from a factor reflecting propensity for (a)synchronicity in concerted H+/e− transfers, which stems directly from the reduction potentials and acidity constants of reactants and products. We show that the most synchronous hydrogen atom abstractions tend to pass over the highest barriers, as exemplified by computations on FeIVO oxidants. Hydrogen atom abstraction (HAA) reactions are cornerstones of chemistry. Various (metallo)enzymes performing the HAA catalysis evolved in nature and inspired the rational development of multiple synthetic catalysts. Still, the factors determining their catalytic efficiency are not fully understood. Herein, we define the simple thermodynamic factor η by employing two thermodynamic cycles: one for an oxidant (catalyst), along with its reduced, protonated, and hydrogenated form; and one for the substrate, along with its oxidized, deprotonated, and dehydrogenated form. It is demonstrated that η reflects the propensity of the substrate and catalyst for (a)synchronicity in concerted H+/e− transfers. As such, it significantly contributes to the activation energies of the HAA reactions, in addition to a classical thermodynamic (Bell–Evans–Polanyi) effect. In an attempt to understand the physicochemical interpretation of η, we discovered an elegant link between η and reorganization energy λ from Marcus theory. We discovered computationally that for a homologous set of HAA reactions, λ reaches its maximum for the lowest |η|, which then corresponds to the most synchronous HAA mechanism. This immediately implies that among HAA processes with the same reaction free energy, ΔG0, the highest barrier (≡ΔG≠) is expected for the most synchronous proton-coupled electron (i.e., hydrogen) transfer. As proof of concept, redox and acidobasic properties of nonheme FeIVO complexes are correlated with activation free energies for HAA from C−H and O−H bonds. We believe that the reported findings may represent a powerful concept in designing new HAA catalysts.

Macrocycle Conformational Sampling by DFT-D3/COSMO-RS Methodology

To find and calibrate a robust and reliable computational protocol for mapping conformational space of medium-sized molecules, exhaustive conformational sampling has been carried out for a series of seven macrocyclic compounds of varying ring size and one acyclic analogue. While five of them were taken from the MD/LLMOD/force field study by Shelley and co-workers ( Watts , K. S. ; Dalal , P. ; Tebben , A. J. ; Cheney , D. L. ; Shelley , J. C. Macrocycle Conformational Sampling with MacroModel . J. Chem. Inf. MODEL 2014 , 54 , 2680 - 2696 ), three represent potential macrocyclic inhibitors of human cyclophilin A. The free energy values (GDFT/COSMO-RS) for all of the conformers of each compound were obtained by a composite protocol based on in vacuo quantum mechanics (DFT-D3 method in a large basis set), standard gas-phase thermodynamics, and the COSMO-RS solvation model. The GDFT/COSMO-RS values were used as the reference for evaluating the performance of conformational sampling algorithms: standard and extended MD/LLMOD search (simulated-annealing molecular dynamics with low-lying eigenvector following algorithms, employing the OPLS2005 force field plus GBSA solvation) available in MacroModel and plain molecular dynamics (MD) sampling at high temperature (1000 K) using the semiempirical quantum mechanics (SQM) potential SQM(PM6-D3H4/COSMO) followed by energy minimization of the snapshots. It has been shown that the former protocol (MD/LLMOD) may provide a more complete set of initial structures that ultimately leads to the identification of a greater number of low-energy conformers (as assessed by GDFT/COSMO-RS) than the latter (i.e., plain SQM MD). The CPU time needed to fully evaluate one medium-sized compound (∼100 atoms, typically resulting in several hundred or a few thousand conformers generated and treated quantum-mechanically) is approximately 1,000-100,000 CPU hours on todays computers, which transforms to 1-7 days on a small-sized computer cluster with a few hundred CPUs. Finally, our data sets based on the rigorous quantum-chemical approach allow us to formulate a composite conformational sampling protocol with multiple checkpoints and truncation of redundant structural data that offers superior insights at affordable computational cost.