Corrected small basis set Hartree-Fock method for large systems

Chiral recognition of liquid phase dimers from gamma-valerolactone racemic mixture

Gamma-valerolactone (GVL) is a cyclic ester that can be considered a green alternative in chemical processes due to its environmentally friendly physical and chemical properties and low production cost from biomass. Although GVL is a chiral solvent, it is usually used as a racemic mixture, instead of its homochiral forms, which might improve its performance in enantioselective synthesis and chiral separation chromatographic techniques. This report presents the development and validation of an atomistic force field optimized to reproduce GVL liquid-phase properties via Monte Carlo (MC) and molecular dynamics (MD) simulation methods. The optimized force field improved the description of the interactions between pairs of molecules, which is a key aspect for a proper assessment of subtle interactions between the enantiomeric forms of GVL. Inspection of radial distribution functions (RDF) for correlations between RR, SS, and RS interactions found within GVL racemic mixture shows very subtle differences at the first solvation shell. Average interaction energies [Formula: see text], [Formula: see text] and [Formula: see text] for RR, SS, and RS dimer ensembles, respectively, were calculated with force field and also HF-3c and PBEh-3c quantum chemistry methods. For each methodology, resulting values obtained for [Formula: see text] and [Formula: see text] were almost the same and more negative than [Formula: see text]. Also, the average energy fluctuation obtained for RR and SS dimers were higher than the one obtained for RS.

Communication: Exact analytical derivatives for the domain-based local pair natural orbital MP2 method (DLPNO-MP2)

Electron correlation methods based on pair natural orbitals (PNOs) have gained an increasing degree of interest in recent years, as they permit energy calculations to be performed on systems containing up to many hundred atoms, while maintaining chemical accuracy for reaction energies. We present an approach for taking exact analytical first derivatives of the energy contributions in the simplest method of the family of Domain-based Local Pair Natural Orbital (DLPNO) methods, closed-shell DLPNO-MP2. The Lagrangian function contains constraints to account for the relaxation of PNOs. RI-MP2 reference geometries are reproduced accurately, as exemplified for four systems with a substantial degree of nonbonding interactions. By the example of electric field gradients, we demonstrate that omitting PNO-specific constraints can lead to dramatic errors for orbital-relaxed properties.

Intermolecular interactions in the condensed phase: Evaluation of semi-empirical quantum mechanical methods

To facilitate further development of approximate quantum mechanical methods for condensed phase applications, we present a new benchmark dataset of intermolecular interaction energies in the solution phase for a set of 15 dimers, each containing one charged monomer. The reference interaction energy in solution is computed via a thermodynamic cycle that integrates dimer binding energy in the gas phase at the coupled cluster level and solute-solvent interaction with density functional theory; the estimated uncertainty of such calculated interaction energy is ±1.5 kcal/mol. The dataset is used to benchmark the performance of a set of semi-empirical quantum mechanical (SQM) methods that include DFTB3-D3, DFTB3/CPE-D3, OM2-D3, PM6-D3, PM6-D3H+, and PM7 as well as the HF-3c method. We find that while all tested SQM methods tend to underestimate binding energies in the gas phase with a root-mean-squared error (RMSE) of 2-5 kcal/mol, they overestimate binding energies in the solution phase with an RMSE of 3-4 kcal/mol, with the exception of DFTB3/CPE-D3 and OM2-D3, for which the systematic deviation is less pronounced. In addition, we find that HF-3c systematically overestimates binding energies in both gas and solution phases. As most approximate QM methods are parametrized and evaluated using data measured or calculated in the gas phase, the dataset represents an important first step toward calibrating QM based methods for application in the condensed phase where polarization and exchange repulsion need to be treated in a balanced fashion.

A new near-linear scaling, efficient and accurate, open-shell domain-based local pair natural orbital coupled cluster singles and doubles theory

The Coupled-Cluster expansion, truncated after single and double excitations (CCSD), provides accurate and reliable molecular electronic wave functions and energies for many molecular systems around their equilibrium geometries. However, the high computational cost, which is well-known to scale as O(N⁶) with system size N, has limited its practical application to small systems consisting of not more than approximately 20-30 atoms. To overcome these limitations, low-order scaling approximations to CCSD have been intensively investigated over the past few years. In our previous work, we have shown that by combining the pair natural orbital (PNO) approach and the concept of orbital domains it is possible to achieve fully linear scaling CC implementations (DLPNO-CCSD and DLPNO-CCSD(T)) that recover around 99.9% of the total correlation energy [C. Riplinger et al., J. Chem. Phys. 144, 024109 (2016)]. The production level implementations of the DLPNO-CCSD and DLPNO-CCSD(T) methods were shown to be applicable to realistic systems composed of a few hundred atoms in a routine, black-box fashion on relatively modest hardware. In 2011, a reduced-scaling CCSD approach for high-spin open-shell unrestricted Hartree-Fock reference wave functions was proposed (UHF-LPNO-CCSD) [A. Hansen et al., J. Chem. Phys. 135, 214102 (2011)]. After a few years of experience with this method, a few shortcomings of UHF-LPNO-CCSD were noticed that required a redesign of the method, which is the subject of this paper. To this end, we employ the high-spin open-shell variant of the N-electron valence perturbation theory formalism to define the initial guess wave function, and consequently also the open-shell PNOs. The new PNO ansatz properly converges to the closed-shell limit since all truncations and approximations have been made in strict analogy to the closed-shell case. Furthermore, given the fact that the formalism uses a single set of orbitals, only a single PNO integral transformation is necessary, which offers large computational savings. We show that, with the default PNO truncation parameters, approximately 99.9% of the total CCSD correlation energy is recovered for open-shell species, which is comparable to the performance of the method for closed-shells. UHF-DLPNO-CCSD shows a linear scaling behavior for closed-shell systems, while linear to quadratic scaling is obtained for open-shell systems. The largest systems we have considered contain more than 500 atoms and feature more than 10 000 basis functions with a triple-ζ quality basis set.

Solid state frustrated Lewis pair chemistry

In solution the PCy3/B(C6F5)3 pair is rapidly deactivated by nucleophilic aromatic substitution. In the solid state deactivation is effectively suppressed and the active frustrated phosphane/borane Lewis pair splits dihydrogen or adds to sulfur dioxide. A variety of phosphane/B(C6F5)3 pairs have been used to carry out active FLP reactions in the solid state. The reactions were analyzed by DFT calculations and by solid state NMR spectroscopy. The solid state dihydrogen splitting reaction was also carried out under near to ambient conditions with suspensions of the non-quenched phosphane/borane mixtures in the fluorous liquid perfluoromethylcyclohexane.

Simplified DFT methods for consistent structures and energies of large systems

Kohn-Sham density functional theory (DFT) is routinely used for the fast electronic structure computation of large systems and will most likely continue to be the method of choice for the generation of reliable geometries in the foreseeable future. Here, we present a hierarchy of simplified DFT methods designed for consistent structures and non-covalent interactions of large systems with particular focus on molecular crystals. The covered methods are a minimal basis set Hartree-Fock (HF-3c), a small basis set screened exchange hybrid functional (HSE-3c), and a generalized gradient approximated functional evaluated in a medium-sized basis set (B97-3c), all augmented with semi-classical correction potentials. We give an overview on the methods design, a comprehensive evaluation on established benchmark sets for geometries and lattice energies of molecular crystals, and highlight some realistic applications on large organic crystals with several hundreds of atoms in the primitive unit cell.

Quantum Chemical Approaches in Structure-Based Virtual Screening and Lead Optimization

Today computational chemistry is a consolidated tool in drug lead discovery endeavors. Due to methodological developments and to the enormous advance in computer hardware, methods based on quantum mechanics (QM) have gained great attention in the last 10 years, and calculations on biomacromolecules are becoming increasingly explored, aiming to provide better accuracy in the description of protein-ligand interactions and the prediction of binding affinities. In principle, the QM formulation includes all contributions to the energy, accounting for terms usually missing in molecular mechanics force-fields, such as electronic polarization effects, metal coordination, and covalent binding; moreover, QM methods are systematically improvable, and provide a greater degree of transferability. In this mini-review we present recent applications of explicit QM-based methods in small-molecule docking and scoring, and in the calculation of binding free-energy in protein-ligand systems. Although the routine use of QM-based approaches in an industrial drug lead discovery setting remains a formidable challenging task, it is likely they will increasingly become active players within the drug discovery pipeline.

Random versus Systematic Errors in Reaction Enthalpies Computed Using Semiempirical and Minimal Basis Set Methods

The connectivity-based hierarchy (CBH) protocol for computing accurate reaction enthalpies developed by Sengupta and Raghavachari is tested for fast ab initio methods (PBEh-3c, HF-3c, and HF/STO-3G), tight-binding density functional theory (DFT) methods (GFN-xTB, DFTB, and DFTB-D3), and neglect-of-diatomic-differential-overlap (NDDO)-based semiempirical methods (AM1, PM3, PM6, PM6-DH+, PM6-D2, PM6-D3H+, PM6-D3H4X, PM7, and OM2) using the same set of 25 reactions as in the original study. For the CBH-2 scheme, which reflects the change in the immediate chemical environment of all of the heavy atoms, the respective mean unsigned error relative to G4 for PBEh-3c, HF-3c, HF/STO-3G, GFN-xTB, DFTB-D3, DFTB, PM3, AM1, PM6, PM6-DH+, PM6-D3, PM6-D3H+, PM6-D3H4X, PM7, and OM2 are 1.9, 2.4, 3.0, 3.9, 3.7, 4.5, 4.8, 5.5, 5.4, 5.3, 5,4, 6.5, 5.3, 5.2, and 5.9 kcal/mol, with a single outlier removed for HF-3c, PM6, PM6-DH+, PM6-D3, PM6-D3H4X, and PM7. The increase in accuracy for the NDDO-based methods is relatively modest due to the random errors in predicted heats for formation.

Influence of Solvent Representation on Nuclear Shielding Calculations of Protonation States of Small Biological Molecules

In this study, we assess the influence of solvation on the accuracy and reliability of isotropic nuclear magnetic shielding calculations for amino acids in comparison to experimental data. We focus particularly on the performance of solvation methods for different protonation states, as biological molecules occur almost exclusively in aqueous solution and are subject to protonation with pH. We identify significant shortcomings of current implicit solvent models and present a hybrid solvation approach that improves agreement with experimental data by taking into account the presence of direct interactions between amino acid protonation state and water molecules.

NAMD goes quantum: an integrative suite for hybrid simulations

Hybrid methods that combine quantum mechanics (QM) and molecular mechanics (MM) can be applied to studies of reaction mechanisms in locations ranging from active sites of small enzymes to multiple sites in large bioenergetic complexes. By combining the widely used molecular dynamics and visualization programs NAMD and VMD with the quantum chemistry packages ORCA and MOPAC, we created an integrated, comprehensive, customizable, and easy-to-use suite (http://www.ks.uiuc.edu/Research/qmmm). Through the QwikMD interface, setup, execution, visualization, and analysis are streamlined for all levels of expertise.

Composite and Low-Cost Approaches for Molecular Crystal Structure Prediction

Molecular crystal structure prediction (CSP) requires evaluating differences in lattice energy between candidate crystal structures accurately and efficiently. In this work, we explore and compare several low-cost alternatives to dispersion-corrected density-functional theory (DFT) in the plane-waves/pseudopotential approximation, the most accurate and general approach used for CSP at present. Three types of low-cost methods are considered: DFT with a small basis set of finite-support numerical orbitals (the SIESTA method), dispersion-corrected Gaussian small or minimal-basis-set Hartree-Fock and DFT with additional empirical corrections (HF-3c and PBEh-3c), and self-consistent-charge dispersion-corrected density-functional tight binding (SCC-DFTB3-D3). In addition, we study the performance of composite methods that comprise a geometry optimization using a low-cost approach followed by a single-point calculation using the accurate but comparatively expensive B86bPBE-XDM functional. All methods were tested for their abilities to produce absolute lattice energies, relative lattice energies, and crystal geometries. We show that assessing various methods by their ability to produce absolute lattice energies can be misleading and that relative lattice energies are a much better indicator of performance in CSP. The EE14 set of relative solubilities of homochiral and heterochiral chiral crystals is proposed for relative lattice-energy benchmarking. Our results show that PBE-D2 plus a DZP basis set of numerical orbitals coupled with a final B86bPBE-XDM single-point calculation gives excellent performance at a fraction of the cost of a full B86bPBE-XDM calculation, although the results are sensitive to the particular details of the computational protocol. The B86bPBE-XDM//PBE-D2/DZP method was subsequently tested in a practical CSP application from our recent work on the crystal structure of the enantiopure and racemate forms of 1-aza[6]helicene, a chiral organic semiconductor. Our results show that this multilevel method is able to correctly reproduce the energy ranking of both crystal forms.

Catalytic enantioselective acyl transfer: the case for 4-PPY with a C-3 carboxamide peptide auxiliary based on synthesis and modelling studies

A series of 4-pyrrolidinopyridine (4-PPY) C-3 carboxamides containing peptide-based side chains have been synthesised and evaluated in the kinetic resolution of a small library of chiral benzylic secondary alcohols. A key design element was the incorporation of a tryptophan residue in the peptide side chain for promoting π-stacking between peptide side chain and the pyridinium ring of the N-acyl intermediate, in which modelling was used as a structure-based guiding tool. Together, a catalyst containing a LeuTrp-N-Boc side chain (catalyst 8) was identified that achieved s-values up to and in slight excess of 10. A transition-state model based on the modelling is proposed to explain the origin of enantioselectivity. This study establishes the usefulness of modelling as a structure-based guiding tool for enantioselectivity optimization as well as the potential for developing scalable peptide-based DMAP-type catalysts for large-scale resolution work.

Quantum-Chemical Insights into the Self-Assembly of Carbon-Based Supramolecular Complexes

Understanding how molecular systems self-assemble to form well-organized superstructures governed by noncovalent interactions is essential in the field of supramolecular chemistry. In the nanoscience context, the self-assembly of different carbon-based nanoforms (fullerenes, carbon nanotubes and graphene) with, in general, electron-donor molecular systems, has received increasing attention as a means of generating potential candidates for technological applications. In these carbon-based systems, a deep characterization of the supramolecular organization is crucial to establish an intimate relation between supramolecular structure and functionality. Detailed structural information on the self-assembly of these carbon-based nanoforms is however not always accessible from experimental techniques. In this regard, quantum chemistry has demonstrated to be key to gain a deep insight into the supramolecular organization of molecular systems of high interest. In this review, we intend to highlight the fundamental role that quantum-chemical calculations can play to understand the supramolecular self-assembly of carbon-based nanoforms through a limited selection of supramolecular assemblies involving fullerene, fullerene fragments, nanotubes and graphene with several electron-rich π-conjugated systems.

Polydopamine and eumelanin models in various oxidation states

The most stable molecular structures of PDA and eumelanin are proposed.

First-principles stability ranking of molecular crystal polymorphs with the DFT+MBD approach

The ability to accurately calculate the relative stabilities of numerous polymorphs of a given molecular crystal is crucial for the success of any molecular crystal structure prediction (CSP) approach. We have recently presented a hierarchical CSP procedure based on van-der-Waals-inclusive density functional theory [Hoja et al., 2018, arXiv:1803.07503], which yields excellent stability rankings for molecular crystals involving rigid molecules, salts, co-crystals, and highly polymorphic drug-like molecules. This approach includes many-body dispersion effects, exact exchange, as well as vibrational free energies. Here, we discuss in detail the impact of these effects on the obtained stability rankings. In addition, we assess the impact of the approximations used in our hierarchical procedure. We show that our procedure is generally robust to 1-2 kJ mol-1 for the systems in the latest CSP blind test but vibrational free energies for crystals involving flexible molecules would benefit from directly including many-body dispersion interactions. In addition, we also discuss the effect of temperature on the structure of molecular crystals and a simple but effective method for estimating anharmonic effects.

Anharmonic vibrational spectra from double incremental potential energy and dipole surfaces

Using incremental approaches, size limitations for property surface generations are pushed significantly, enabling accurate large molecule anharmonic vibrational spectra calculations.

Computational Chemistry: The Fate of Current Methods and Future Challenges

"Where do we go from here?" is the underlying question regarding the future (perhaps foreseeable) developments in computational chemistry. Although this young discipline has already permeated practically all of chemistry, it is likely to become even more powerful with the rapid development of computational hard- and software.

Intramolecular London Dispersion Interaction Effects on Gas-Phase and Solid-State Structures of Diamondoid Dimers

The covalent diamantyl (C₂₈H₃₈) and oxadiamantyl (C₂₆H₃₄O₂) dimers are stabilized by London dispersion attractions between the dimer moieties. Their solid-state and gas-phase structures were studied using a multitechnique approach, including single-crystal X-ray diffraction (XRD), gas-phase electron diffraction (GED), a combined GED/microwave (MW) spectroscopy study, and quantum chemical calculations. The inclusion of medium-range electron correlation as well as the London dispersion energy in density functional theory is essential to reproduce the experimental geometries. The conformational dynamics computed for C₂₆H₃₄O₂ agree well with solution NMR data and help in the assignment of the gas-phase MW data to individual diastereomers. Both in the solid state and the gas phase the central C-C bond is of similar length for the diamantyl [XRD, 1.642(2) Å; GED, 1.630(5) Å] and the oxadiamantyl dimers [XRD, 1.643(1) Å; GED, 1.632(9) Å; GED+MW, 1.632(5) Å], despite the presence of two oxygen atoms. Out of a larger series of quantum chemical computations, the best match with the experimental reference data is achieved with the PBEh-3c, PBE0-D3, PBE0, B3PW91-D3, and M06-2X approaches. This is the first gas-phase confirmation that the markedly elongated C-C bond is an intrinsic feature of the molecule and that crystal packing effects have only a minor influence.

Efficient and Accurate Born-Oppenheimer Molecular Dynamics for Large Molecular Systems

An efficient scheme for the calculation of Born-Oppenheimer molecular dynamics (BOMD) simulations is introduced. It combines the corrected small basis set Hartree-Fock (HF-3c) method by Sure and Grimme [J. Comput. Chem. 2013, 43, 1672], extended Lagrangian BOMD (XL-BOMD) by Niklasson et al. [J. Chem. Phys. 2009, 130, 214109], and the calculation of the two electron integrals on graphics processing units (GPUs) [J. Chem. Phys. 2013, 138, 134114; J. Chem. Theory Comput. 2015, 11, 918]. To explore the parallel performance of our strong scaling implementation of the method, we present timings and extract, as its validation and first illustrative application, high-quality vibrational spectra from simulated trajectories of β-carotene, paclitaxel, and liquid water (up to 500 atoms). We conclude that the presented BOMD scheme may be used as a cost-efficient and reliable tool for computing vibrational spectra and thermodynamics of large molecular systems including explicit solvent molecules containing 500 atoms and more. Simulating 50 ps of maitotoxin (nearly 500 atoms) employing time steps of 0.5 fs requires ∼3 weeks on 12 CPUs (Intel Xeon E5 2620 v3) with 24 GPUs (AMD FirePro 3D W8100).

Thermal Expansion of Carbamazepine: Systematic Crystallographic Measurements Challenge Quantum Chemical Calculations

We report systematic temperature-dependent X-ray measurements on the most stable carbamazepine polymorph. This active pharmaceutical ingredient is used to demonstrate how the thermal expansion can probe certain intermolecular interactions resulting in anisotropic expansion behavior. We show that most structural features can be captured by electronic structure calculations at the quasi-harmonic approximation (QHA) provided a dispersion-corrected density functional based method is employed. The impact of thermal expansion on the phonon modes and hence free energy contributions is large enough to impact the relative stability of different polymorphs.

Synthesis and DFT Calculations of Novel Vanillin-Chalcones and Their 3-Aryl-5-(4-(2-(dimethylamino)-ethoxy)-3-methoxyphenyl)-4,5-dihydro-1H-pyrazole-1-carbaldehyde Derivatives as Antifungal Agents

Novel (E)-1-(aryl)-3-(4-(2-(dimethylamino)ethoxy)-3-methoxyphenyl) prop-2-en-1-ones 4 were synthesized by a Claisen-Schmidt reaction of 4-(2-(dimethylamino)ethoxy)-3-methoxy-benzaldehyde (2) with several acetophenone derivatives 3. Subsequently, cyclocondensation reactions of chalcones 4 with hydrazine hydrate afforded the new racemic 3-aryl-5-(4-(2-(dimethylamino)ethoxy)-3-methoxyphenyl)-4,5-dihydro-1H-pyrazole-1-carbaldehydes 5 when the reaction was carried out in formic acid. The antifungal activity of both series of compounds against eight fungal species was determined. In general, chalcone derivatives 4 showed better activities than pyrazolines 5 against all tested fungi. None of the compounds 4a–g and 5a–g showed activity against the three Aspergillus spp. In contrast, most of the compounds 4 showed moderate to high activities against three dermatophytes (MICs 31.25–62.5 µg/mL), being 4a followed by 4c the most active structures. Interestingly, 4a and 4c possess fungicidal rather than fungistatic activities, with MFC values between 31.25 and 62.5 μg/mL. The comparison of the percentages of inhibition of C. neoformans by the most active compounds 4, allowed us to know the role played by the different substituents of the chalcones’ A-ring. Also the most anti-cryptococcal compounds 4a–c and 4g, were tested in a second panel of five clinical C. neoformans strains in order to have an overview of their inhibition capacity not only of standardized but also of clinical C. neoformans strains. DFT calculations showed that the electrophilicity is the main electronic property to explain the differences in antifungal activities for the synthesized chalcones and pyrazolines compounds. Furthermore, a quantitative reactivity analysis showed that electron-withdrawing substituted chalcones presented the higher electrophilic character and hence, the greater antifungal activities among compounds of series 4.

On the hydrogen activation by frustrated Lewis pairs in the solid state: benchmark studies and theoretical insights

Recently, the concept of small molecule activation by frustrated Lewis pairs (FLPs) has been expanded to the solid state showing a variety of interesting reactivities. Therefore, there is a need to establish a computational protocol to investigate such systems theoretically. In the present study, we selected several FLPs and applied multiple levels of theory, ranging from a semi-empirical tight-binding Hamiltonian to dispersion corrected hybrid density functionals. Their performance is benchmarked for the computation of crystal geometries, thermostatistical contributions, and reaction energies. We show that the computationally efficient HF-3c method gives accurate crystal structures and is numerically stable and sufficiently fast for routine applications. This method also gives reliable values for the thermostatistical contributions to Gibbs free energies. The meta-generalized gradient approximated TPSS-D3 evaluated in a projector augmented plane wave basis set is able to produce sufficiently accurate reaction electronic energies. The established protocol is intended to support experimental studies and to predict new reactions in the emerging field of solid-state FLPs.This article is part of the themed issue 'Frustrated Lewis pair chemistry'.

Accuracy of Frequencies Obtained with the Aid of Explicitly Correlated Wave Function Based Methods

We asses the basis set convergence of harmonic frequencies using different explicitly correlated wave function based methods. All commonly available CCSD(T) variants as well as MP2-F12 and MP4(F12*) are considered, and a hierarchy of the different approaches is established. As for reaction and atomization energies, CCSD(F12*)(T*) is a close approximation to CCSD(F12)(T*) and clearly superior to the other tested approximations. The used scaling for the triples correction enhances the accuracy relative to CCSD(F12*)(T) especially for small basis sets and is very attractive since no additional computational costs are added. However, this scaling slightly breaks size consistency, and therefore we additionally study the accuracy of CCSD(F12*)(T*) and CCSD(F12*)(T) in the context of calculating anharmonic frequencies to check if this causes problems in the generation of the potential energy surface (PES). We find a fast basis set convergence for harmonic and anharmonic frequencies. Already in the cc-pVDZ-F12 basis, the RMSD to the CBS limit is only around 4-5 cm^-1.

Psi4 1.1: An Open-Source Electronic Structure Program Emphasizing Automation, Advanced Libraries, and Interoperability

Psi4 is an ab initio electronic structure program providing methods such as Hartree-Fock, density functional theory, configuration interaction, and coupled-cluster theory. The 1.1 release represents a major update meant to automate complex tasks, such as geometry optimization using complete-basis-set extrapolation or focal-point methods. Conversion of the top-level code to a Python module means that Psi4 can now be used in complex workflows alongside other Python tools. Several new features have been added with the aid of libraries providing easy access to techniques such as density fitting, Cholesky decomposition, and Laplace denominators. The build system has been completely rewritten to simplify interoperability with independent, reusable software components for quantum chemistry. Finally, a wide range of new theoretical methods and analyses have been added to the code base, including functional-group and open-shell symmetry adapted perturbation theory, density-fitted coupled cluster with frozen natural orbitals, orbital-optimized perturbation and coupled-cluster methods (e.g., OO-MP2 and OO-LCCD), density-fitted multiconfigurational self-consistent field, density cumulant functional theory, algebraic-diagrammatic construction excited states, improvements to the geometry optimizer, and the "X2C" approach to relativistic corrections, among many other improvements.

Tannins and Tannin-Related Derivatives Enhance the (Pseudo-)Halogenating Activity of Lactoperoxidase

Several hydrolyzable tannins, proanthocyanidins, tannin derivatives, and a tannin-rich plant extract of tormentil rhizome were tested for their potential to regenerate the (pseudo-)halogenating activity, i.e., the oxidation of SCN^- to hypothiocyanite ^-OSCN, of lactoperoxidase (LPO) after hydrogen peroxide-mediated enzyme inactivation. Measurements were performed using 5-thio-2-nitrobenzoic acid in the presence of tannins and related substances in order to determine kinetic parameters and to trace the LPO-mediated ^-OSCN formation. The results were combined with docking studies and molecular orbital analysis. The ^-OSCN-regenerating effect of tannin derivatives relates well with their binding properties toward LPO as well as their occupied molecular orbitals. Especially simple compounds like ellagic acid or methyl gallate and the complex plant extract were found as potent enzyme-regenerating compounds. As the (pseudo-)halogenating activity of LPO contributes to the maintenance of oral bacterial homeostasis, the results provide new insights into the antibacterial mode of action of tannins and related compounds. Furthermore, chemical properties of the tested compounds that are important for efficient enzyme-substrate interaction and regeneration of the ^-OSCN formation by LPO were identified.