CREST-A program for the exploration of low-energy molecular chemical space

Conformer-rotamer sampling tool (CREST) is an open-source program for the efficient and automated exploration of molecular chemical space. Originally developed in Pracht et al. [Phys. Chem. Chem. Phys. 22, 7169 (2020)] as an automated driver for calculations at the extended tight-binding level (xTB), it offers a variety of molecular- and metadynamics simulations, geometry optimization, and molecular structure analysis capabilities. Implemented algorithms include automated procedures for conformational sampling, explicit solvation studies, the calculation of absolute molecular entropy, and the identification of molecular protonation and deprotonation sites. Calculations are set up to run concurrently, providing efficient single-node parallelization. CREST is designed to require minimal user input and comes with an implementation of the GFNn-xTB Hamiltonians and the GFN-FF force-field. Furthermore, interfaces to any quantum chemistry and force-field software can easily be created. In this article, we present recent developments in the CREST code and show a selection of applications for the most important features of the program. An important novelty is the refactored calculation backend, which provides significant speed-up for sampling of small or medium-sized drug molecules and allows for more sophisticated setups, for example, quantum mechanics/molecular mechanics and minimum energy crossing point calculations.

Do Optimally Tuned Range-Separated Hybrid Functionals Require a Reparametrization of the Dispersion Correction? It Depends

For ground- and excited-state studies of large molecules, it is the state of the art to combine (time-dependent) DFT with dispersion-corrected range-separated hybrid functionals (RSHs), which ensures an asymptotically correct description of exchange effects and London dispersion. Specifically for studying excited states, it is common practice to tune the range-separation parameter ω (optimal tuning), which can further improve the accuracy. However, since optimal tuning essentially changes the functional, it is unclear if and how much the parameters used for the dispersion correction depend on the chosen ω value. To answer this question, we explore this interdependency by refitting the DFT-D4 dispersion model for six established RSHs over a wide range of ω values (0.05-0.45 a0-1) using a set of noncovalently bound molecular complexes. The results reveal some surprising differences among the investigated functionals: While PBE-based RSHs and ωB97M-D4 generally exhibit a weak interdependency and robust performance over a wide range of ω values, B88-based RSHs, specifically LC-BLYP, are strongly affected. For these, even a minor reduction of ω from the default value manifests in strong systematic overbinding and poor performance in the typical range of optimally tuned ω values. Finally, we discuss strategies to mitigate these issues and reflect the results in the context of the employed D4 parameter optimization algorithm and fit set, outlining strategies for future improvements.

High-throughput screening of spin states for transition metal complexes with spin-polarized extended tight-binding methods

The semiempirical GFNn-xTB ( n = 1 , 2 ) tight-binding methods are extended with a spin-dependent energy term (spin-polarization), enabling the fast and efficient screening of different spin states for transition metal complexes. While GFNn-xTB methods inherently can not differentiate properly between high-spin (HS) and low-spin (LS) states, this shortcoming is corrected with the presented methods termed spGFNn-xTB. The performance of spGFNn-xTB methods for spin state energy splittings is evaluated on a newly compiled benchmark set of 90 complexes (27 HS and 63 LS complexes) containing 3d, 4d, and 5d transition metals (termed TM90S) employing DFT references at the TPSSh-D4/def2-QZVPP level of theory. The challenging TM90S set contains complexes with charges between - 4 and +3, spin multiplicities between 1 and 6, and spin-splitting energies that range from - 47.8 to 146.6 kcal/mol with a mean average of 32.2 kcal/mol. On this set the (sp)GFNn-xTB methods, the PM6-D3H4 method, and the PM7 method are evaluated with spGFN1-xTB yielding the lowest MAD of 19.6 kcal/mol followed by spGFN2-xTB with 24.8 kcal/mol. While for the 4d and 5d subsets small or no improvements are observed with spin-polarization, large improvements are obtained for the 3d subset with spGFN1-xTB yielding the smallest MAD of 14.2 kcal/mol followed by spGFN2-xTB with 17.9 kcal/mol and PM6-D3H4 with 28.4 kcal/mol. The correct sign of the spin state splittings is obtained with spGFN2-xTB in 89% of all cases closely followed by spGFN1-xTB with 88%. On the full set, a pure semiempirical vertical spGFN2-xTB//GFN2-xTB-based workflow for screening purposes yields a slightly better MAD of 22.2 kcal/mol due to error compensation, while being qualitative correct for one additional case. In combination with their low computational cost (scanning spin states in seconds), the spGFNn-xTB methods represent robust tools for pre-screening steps of spin state calculations and high-throughput workflows.

Extended Conductor-like Polarizable Continuum Solvation Model (CPCM-X) for Semiempirical Methods

We have developed a new method to accurately account for solvation effects in semiempirical quantum mechanics based on a polarizable continuum model (PCM). The extended conductor-like polarizable continuum model (CPCM-X) incorporates a computationally efficient domain decomposition conductor-like screening model (ddCOSMO) for extended tight binding (xTB) methods and uses a post-processing approach based on established solvation models, like the conductor-like screening model for real solvents (COSMO-RS) and the universal solvent model based on solute electron density (SMD). According to various benchmarks, the approach performs well across a broad range of systems and applications, including hydration free energies, non-aqueous solvation free energies, and large supramolecular association reactions of neutral and charged species. Our method for computing solvation free energies is much more accurate than the current methods in the xtb program package. It improves the accuracy of solvation free energies by up to 40% for larger supramolecular association reactions to match even the accuracy of higher-level DFT-based solvation models like COSMO-RS and SMD while being computationally more than 2 orders of magnitude faster. The proposed method and the underlying ddCOSMO model are readily available for a wide variety of solvents and are accessible in xtb for use in various computational applications.

ONIOM meets : efficient, accurate, and robust multi-layer simulations across the periodic table

The computational treatment of large molecular structures is of increasing interest in fields of modern chemistry. Accordingly, efficient quantum chemical approaches are needed to perform sophisticated investigations on such systems. This engaged the development of the well-established "Our own N-layered integrated molecular orbital and molecular mechanics" (ONIOM) multi-layer scheme [L. W. Chung et al., Chem. Rev., 2015, 115, 5678-5796]. In this work, we present the specific implementation of the ONIOM scheme into the xtb semi-empirical extended tight-binding program package and its application to challenging transition-metal complexes. The efficient and broadly applicable GFNn-xTB and -FF methods are applied in the ONIOM framework to elucidate reaction energies, geometry optimizations, and explicit solvation effects for metal-organic systems with up to several hundreds of atoms. It is shown that an ONIOM-based combination of density functional theory, semi-empirical, and force-field methods can be used to drastically reduce the computational costs and thus enable the investigation of huge systems at almost no significant loss in accuracy.

TBMaLT, a flexible toolkit for combining tight-binding and machine learning

Tight-binding approaches, especially the Density Functional Tight-Binding (DFTB) and the extended tight-binding schemes, allow for efficient quantum mechanical simulations of large systems and long-time scales. They are derived from ab initio density functional theory using pragmatic approximations and some empirical terms, ensuring a fine balance between speed and accuracy. Their accuracy can be improved by tuning the empirical parameters using machine learning techniques, especially when information about the local environment of the atoms is incorporated. As the significant quantum mechanical contributions are still provided by the tight-binding models, and only short-ranged corrections are fitted, the learning procedure is typically shorter and more transferable as it were with predicting the quantum mechanical properties directly with machine learning without an underlying physically motivated model. As a further advantage, derived quantum mechanical quantities can be calculated based on the tight-binding model without the need for additional learning. We have developed the open-source framework-Tight-Binding Machine Learning Toolkit-which allows the easy implementation of such combined approaches. The toolkit currently contains layers for the DFTB method and an interface to the GFN1-xTB Hamiltonian, but due to its modular structure and its well-defined interfaces, additional atom-based schemes can be implemented easily. We are discussing the general structure of the framework, some essential implementation details, and several proof-of-concept applications demonstrating the perspectives of the combined methods and the functionality of the toolkit.

Conformational Energy Benchmark for Longer n-Alkane Chains

We present the first benchmark set focusing on the conformational energies of highly flexible, long n-alkane chains, termed ACONFL. Unbranched alkanes are ubiquitous building blocks in nature, so the goal is to be able to calculate their properties most accurately to improve the modeling of, e.g., complex (biological) systems. Very accurate DLPNO-CCSD(T1)/CBS reference values are provided, which allow for a statistical meaningful evaluation of even the best available density functional methods. The performance of established and modern (dispersion corrected) density functionals is comprehensively assessed. The recently introduced r²SCAN-V functional shows excellent performance, similar to efficient composite DFT methods like B97-3c and r²SCAN-3c, which provide an even better cost-accuracy ratio, while almost reaching the accuracy of much more computationally demanding hybrid or double hybrid functionals with large QZ AO basis sets. In addition, we investigated the performance of common wave function methods, where MP2/CBS surprisingly performs worse compared to the simple D4 dispersion corrected Hartree-Fock. Furthermore, we investigate the performance of several semiempirical and force field methods, which are commonly used for the generation of conformational ensembles in multilevel workflows or in large scale molecular dynamics studies. Outstanding performance is obtained by the recently introduced general force field, GFN-FF, while other commonly applied methods like the universal force field yield large errors. We recommend the ACONFL as a helpful benchmark set for parametrization of new semiempirical or force field methods and machine learning potentials as well as a meaningful validation set for newly developed DFT or dispersion methods.

Dispersion corrected r2SCAN based global hybrid functionals: r2SCANh, r2SCAN0, and r2SCAN50

The regularized and restored semilocal meta-generalized gradient approximation (meta-GGA) exchange-correlation functional r²SCAN [Furness et al., J. Phys. Chem. Lett. 11, 8208-8215 (2020)] is used to create three global hybrid functionals with varying admixtures of Hartree-Fock "exact" exchange (HFX). The resulting functionals r²SCANh (10% HFX), r²SCAN0 (25% HFX), and r²SCAN50 (50% HFX) are combined with the semi-classical D4 London dispersion correction. The new functionals are assessed for the calculation of molecular geometries, main-group, and metalorganic thermochemistry at 26 comprehensive benchmark sets. These include the extensive GMTKN55 database, ROST61, and IONPI19 sets. It is shown that a moderate admixture of HFX leads to relative improvements of the mean absolute deviations for thermochemistry of 11% (r²SCANh-D4), 16% (r²SCAN0-D4), and 1% (r²SCAN50-D4) compared to the parental semi-local meta-GGA. For organometallic reaction energies and barriers, r²SCAN0-D4 yields an even larger mean improvement of 35%. The computation of structural parameters (geometry optimization) does not systematically profit from the HFX admixture. Overall, the best variant r²SCAN0-D4 performs well for both main-group and organometallic thermochemistry and is better or on par with well-established global hybrid functionals, such as PW6B95-D4 or PBE0-D4. Regarding systems prone to self-interaction errors (SIE4x4), r²SCAN0-D4 shows reasonable performance, reaching the quality of the range-separated ωB97X-V functional. Accordingly, r²SCAN0-D4 in combination with a sufficiently converged basis set [def2-QZVP(P)] represents a robust and reliable choice for general use in the calculation of thermochemical properties of both main-group and organometallic chemistry.

Simultaneous on-line vacuum single- and multi-photon ionization on an orthogonal acceleration time-of-flight mass spectrometer platform

New instrumental development for robust process monitoring with two soft ionization methods working in parallel.

Quantum Chemistry Common Driver and Databases (QCDB) and Quantum Chemistry Engine (QCEngine): Automation and interoperability among computational chemistry programs

Community efforts in the computational molecular sciences (CMS) are evolving toward modular, open, and interoperable interfaces that work with existing community codes to provide more functionality and composability than could be achieved with a single program. The Quantum Chemistry Common Driver and Databases (QCDB) project provides such capability through an application programming interface (API) that facilitates interoperability across multiple quantum chemistry software packages. In tandem with the Molecular Sciences Software Institute and their Quantum Chemistry Archive ecosystem, the unique functionalities of several CMS programs are integrated, including CFOUR, GAMESS, NWChem, OpenMM, Psi4, Qcore, TeraChem, and Turbomole, to provide common computational functions, i.e., energy, gradient, and Hessian computations as well as molecular properties such as atomic charges and vibrational frequency analysis. Both standard users and power users benefit from adopting these APIs as they lower the language barrier of input styles and enable a standard layout of variables and data. These designs allow end-to-end interoperable programming of complex computations and provide best practices options by default.

Calculation of improved enthalpy and entropy of vaporization by a modified partition function in quantum cluster equilibrium theory

In this work, we present an altered partition function that leads to an improved calculation of the enthalpy and entropy of vaporization in the framework of quantum cluster equilibrium theory. The changes are based on a previously suggested modification [S. Grimme, Chem. Eur. J. 18, 9955-9964 (2012)] of the molecular entropy calculation in the gas phase. Here, the low energy vibrational frequencies in the vibrational partition function are treated as hindered rotations instead of vibrations. The new scheme is tested on a set of nine organic solvents for the calculation of the enthalpy and entropy of vaporization. The enthalpies and entropies of vaporization show improvements from 6.5 error to 3.3 kJ mol^-1 deviation to experiment and from 28.4 error to 13.5 J mol^-1 K^-1 deviation to experiment, respectively. The effect of the corrected partition function is visible in the different populations of clusters, which become physically more meaningful in that larger clusters are higher populated in the liquid phase and the gas phase is mainly populated by the monomers. Furthermore, the corrected partition function also overcomes technical difficulties and leads to an increased stability of the calculations in regard to the size of the cluster set.

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

Frustrated Lewis-Pair Neighbors at the Xanthene Framework: Epimerization at Phosphorus and Cooperative Formation of Macrocyclic Adduct Structures

Attachment of a pair of P-stereogenic mesityl(alkynyl)phosphanyl groups at the 4- and 5-positions of a 9,9-dimethylxanthene framework gave mixtures of the respective rac- and meso-bisphosphanyl diastereoisomers. They slowly epimerized in a thermally induced reaction with Gibbs activation barriers of about 25 kcal mol^-1 at room temperature (measured and DFT calculated). The reaction of the meso-mesityl(tert-butylethynyl)phosphanyl derivative with two molar equivalents of Piers' borane [HB(C₆ F₅ )₂ ] led to the formation of the alkylidene-bridged geminal bisphosphane/borane-frustrated Lewis pair system. The compound was obtained enriched (>85 %) in the rac diastereoisomer. With a variety of bifunctional donor substrates, the rac-bis-P/B FLP formed macrocyclic compounds. They were all formally derived from meso-configurated diastereoisomers of the bisphosphanylxanthene backbone.

Robust and Efficient Implicit Solvation Model for Fast Semiempirical Methods

We present a robust and efficient method to implicitly account for solvation effects in modern semiempirical quantum mechanics and force fields. A computationally efficient yet accurate solvation model based on the analytical linearized Poisson-Boltzmann (ALPB) model is parameterized for the extended tight binding (xTB) and density functional tight binding (DFTB) methods as well as for the recently proposed GFN-FF general force field. The proposed methods perform well over a broad range of systems and applications, from conformational energies over transition-metal complexes to large supramolecular association reactions of charged species. For hydration free energies of small molecules, GFN1-xTB(ALPB) is reaching the accuracy of sophisticated explicitly solvated approaches, with a mean absolute deviation of only 1.4 kcal/mol compared to the experiment. Logarithmic octanol-water partition coefficients (log K_ow) are computed with a mean absolute deviation of about 0.65 using GFN2-xTB(ALPB) compared to experimental values indicating a consistent description of differential solvent effects. Overall, more than twenty solvents for each of the six semiempirical methods are parameterized and tested. They are readily available in the xtb and dftb+ programs for diverse computational applications.

Puff-resolved analysis and selected quantification of chemicals in the gas phase of E-Cigarettes, Heat-not-Burn devices and conventional cigarettes using single photon ionization time-of-flight mass spectrometry (SPI-TOFMS): A comparative study

INTRODUCTION

A wide array of alternative nicotine delivery devices (ANDD) has been developed and they are often described as less harmful than combustible cigarettes. This work compares the chemical emissions of three ANDD in comparison to cigarette smoke. All the tested ANDD are characterized by not involving combustion of tobacco.

METHOD

Single photon ionization time-of-flight mass spectrometry (SPI-TOFMS) is coupled to a linear smoking machine, which allows a comprehensive, online analysis of the gaseous phase of the ANDD aerosol and the conventional cigarette smoke (CC). The following devices were investigated in this study: a tobacco cigarette with a glowing piece of coal as a heating source, an electric device for heating tobacco and a first-generation electronic cigarette. Data obtained from a standard 2R4F research cigarette are taken as a reference.

RESULTS

The puff-by-puff profile of all products was recorded. The ANDD show a substantial reduction or complete absence of known harmful and potentially harmful substances compared to the CC. In addition, tar substances (i.e. semivolatile and low volatile aromatic and phenolic compounds) are formed to a much lower extent. Nicotine, however, is supplied in comparable amounts except for the investigated electronic cigarette.

CONCLUSIONS

The data shows that consumers switching from CC to ANDD are exposed to lower concentrations of harmful and potentially harmful substances. However, toxicological and epidemiological studies must deliver conclusive results if these reduced exposures are beneficial for users.

IMPLICATION

The comparison of puff-resolved profiles of emissions from different tobacco products, traditional and alternative, may help users switch to lower emission products. Puff-resolved comparison overcomes technical changes, use modes between products and may help in their regulation.

r2SCAN-3c: A "Swiss army knife" composite electronic-structure method

The recently proposed r²SCAN meta-generalized-gradient approximation (mGGA) of Furness and co-workers is used to construct an efficient composite electronic-structure method termed r²SCAN-3c. To this end, the unaltered r²SCAN functional is combined with a tailor-made triple-ζ Gaussian atomic orbital basis set as well as with refitted D4 and geometrical counter-poise corrections for London-dispersion and basis set superposition error. The performance of the new method is evaluated for the GMTKN55 database covering large parts of chemical space with about 1500 data points, as well as additional benchmarks for non-covalent interactions, organometallic reactions, and lattice energies of organic molecules and ices, as well as for the adsorption on polar salt and non-polar coinage-metal surfaces. These comprehensive tests reveal a spectacular performance and robustness of r²SCAN-3c: It by far surpasses its predecessor B97-3c at only twice the cost and provides one of the best results of all semi-local density-functional theory (DFT)/QZ methods ever tested for the GMTKN55 database at one-tenth of the cost. Specifically, for reaction and conformational energies as well as non-covalent interactions, it outperforms prominent hybrid-DFT/QZ approaches at two to three orders of magnitude lower cost. Perhaps, the most relevant remaining issue of r²SCAN-3c is self-interaction error (SIE), owing to its mGGA nature. However, SIE is slightly reduced compared to other (m)GGAs, as is demonstrated in two examples. After all, this remarkably efficient and robust method is chosen as our new group default, replacing previous composite DFT and partially even expensive high-level methods in most standard applications for systems with up to several hundreds of atoms.

DFTB+, a software package for efficient approximate density functional theory based atomistic simulations

DFTB+ is a versatile community developed open source software package offering fast and efficient methods for carrying out atomistic quantum mechanical simulations. By implementing various methods approximating density functional theory (DFT), such as the density functional based tight binding (DFTB) and the extended tight binding method, it enables simulations of large systems and long timescales with reasonable accuracy while being considerably faster for typical simulations than the respective ab initio methods. Based on the DFTB framework, it additionally offers approximated versions of various DFT extensions including hybrid functionals, time dependent formalism for treating excited systems, electron transport using non-equilibrium Green's functions, and many more. DFTB+ can be used as a user-friendly standalone application in addition to being embedded into other software packages as a library or acting as a calculation-server accessed by socket communication. We give an overview of the recently developed capabilities of the DFTB+ code, demonstrating with a few use case examples, discuss the strengths and weaknesses of the various features, and also discuss on-going developments and possible future perspectives.

Extension and evaluation of the D4 London-dispersion model for periodic systems

We present an extension of the DFT-D4 model [J. Chem. Phys., 2019, 150, 154122] for periodic systems. The main new ingredients are additional reference polarizabilities for highly-coordinated group 1-5 elements derived from pseudo-periodic electrostatically-embedded cluster calculations. To illustrate the performance of the updated method, several test cases are considered, for which we compare D4 to its predecessor D3(BJ), as well as to a comprehensive set of other dispersion-corrected methods. The largest improvements are observed for solid-state polarizabilities of 16 inorganic salts, where the D4 model achieves an unprecedented accuracy, surpassing its predecessor as well as other, computationally much more demanding approaches. For cell volumes and lattice energies of two sets of chemically diverse molecular crystals, the accuracy gain is less pronounced compared to the already excellently performing D3(BJ) method. For the challenging adsorption energies of small organic molecules on metallic as well as on ionic surfaces, DFT-D4 provides values in good agreement with experimental and/or high-level references. These results suggest the application of the proposed D4 model as a physically improved yet computationally efficient dispersion correction for standard DFT calculations as well as low-cost approaches like semi-empirical or even force-field models.

A generally applicable atomic-charge dependent London dispersion correction

The so-called D4 model is presented for the accurate computation of London dispersion interactions in density functional theory approximations (DFT-D4) and generally for atomistic modeling methods. In this successor to the DFT-D3 model, the atomic coordination-dependent dipole polarizabilities are scaled based on atomic partial charges which can be taken from various sources. For this purpose, a new charge-dependent parameter-economic scaling function is designed. Classical charges are obtained from an atomic electronegativity equilibration procedure for which efficient analytical derivatives with respect to nuclear positions are developed. A numerical Casimir-Polder integration of the atom-in-molecule dynamic polarizabilities then yields charge- and geometry-dependent dipole-dipole dispersion coefficients. Similar to the D3 model, the dynamic polarizabilities are precomputed by time-dependent DFT and all elements up to radon (Z = 86) are covered. The two-body dispersion energy expression has the usual sum-over-atom-pairs form and includes dipole-dipole as well as dipole-quadrupole interactions. For a benchmark set of 1225 molecular dipole-dipole dispersion coefficients, the D4 model achieves an unprecedented accuracy with a mean relative deviation of 3.8% compared to 4.7% for D3. In addition to the two-body part, three-body effects are described by an Axilrod-Teller-Muto term. A common many-body dispersion expansion was extensively tested, and an energy correction based on D4 polarizabilities is found to be advantageous for larger systems. Becke-Johnson-type damping parameters for DFT-D4 are determined for more than 60 common density functionals. For various standard energy benchmark sets, DFT-D4 slightly but consistently outperforms DFT-D3. Especially for metal containing systems, the introduced charge dependence of the dispersion coefficients improves thermochemical properties. We suggest (DFT-)D4 as a physically improved and more sophisticated dispersion model in place of DFT-D3 for DFT calculations as well as other low-cost approaches like semi-empirical models.

GFN2-xTB-An Accurate and Broadly Parametrized Self-Consistent Tight-Binding Quantum Chemical Method with Multipole Electrostatics and Density-Dependent Dispersion Contributions

An extended semiempirical tight-binding model is presented, which is primarily designed for the fast calculation of structures and noncovalent interaction energies for molecular systems with roughly 1000 atoms. The essential novelty in this so-called GFN2-xTB method is the inclusion of anisotropic second order density fluctuation effects via short-range damped interactions of cumulative atomic multipole moments. Without noticeable increase in the computational demands, this results in a less empirical and overall more physically sound method, which does not require any classical halogen or hydrogen bonding corrections and which relies solely on global and element-specific parameters (available up to radon, Z = 86). Moreover, the atomic partial charge dependent D4 London dispersion model is incorporated self-consistently, which can be naturally obtained in a tight-binding picture from second order density fluctuations. Fully analytical and numerically precise gradients (nuclear forces) are implemented. The accuracy of the method is benchmarked for a wide variety of systems and compared with other semiempirical methods. Along with excellent performance for the "target" properties, we also find lower errors for "off-target" properties such as barrier heights and molecular dipole moments. High computational efficiency along with the improved physics compared to its precursor GFN-xTB makes this method well-suited to explore the conformational space of molecular systems. Significant improvements are furthermore observed for various benchmark sets, which are prototypical for biomolecular systems in aqueous solution.

Understanding and Quantifying London Dispersion Effects in Organometallic Complexes

Quantum chemical methods are nowadays able to determine properties of larger chemical systems with high accuracy and Kohn-Sham density functional theory (DFT) in particular has proven to be robust and suitable for everyday applications of electronic structure theory. A clear disadvantage of many established standard density functional approximations like B3LYP is their inability to describe long-range electron correlation effects. The inclusion of such effects, also termed London dispersion, into DFT has been extensively researched in recent years, resulting in some efficient and routinely used correction schemes. The well-established D3 method has demonstrated its efficiency and accuracy in numerous applications since 2010. Recently, it was improved by developing the successor (termed D4) which additionally includes atomic partial charge information for the generation of pairwise dispersion coefficients. These coefficients determine the leading-order (two-body) and higher-order (three- or many-body) terms of the D4 dispersion energy which is simply added to a standard DFT energy. With its excellent accuracy-to-cost ratio, the DFT-D4 method is well suited for the determination of structures and chemical properties for molecules of most kinds. While dispersion effects in organic molecules are nowadays well studied, much less is known for organometallic complexes. For such systems, there has been a growing interest in designing dispersion-controlled reactions especially in the field of homogeneous catalysis. Here, efficient electronic structure methods are necessary for screening of promising model complexes and quantifying dispersion effects. In this Account, we describe the quality of calculated structural and thermodynamic properties in gas-phase obtained with DFT-D4 corrected methods, specifically for organometallic complexes. The physical effects leading to London dispersion interactions are briefly discussed in the picture of second-order perturbation theory. Subsequently, basic theoretical aspects of the D4 method are introduced followed by selected case studies. Several chemical examples are presented starting with the analysis of transition metal thermochemistry and noncovalent interactions for small, heavy element containing main group compounds. Computed reaction energies can only match highly accurate reference values when all energy contributions are included in the DFT treatment, thus highlighting the major role of dispersion interactions for the accurate description of thermochemistry in gas-phase. Furthermore, the correlation between structural and catalytic properties is emphasized where the accessibility of high quality structures is essential for reaction planning and catalyst design. We present calculations for aggregates of organometallic systems with intrinsically large repulsive electrostatic interactions which can be stabilized by London dispersion effects. The newly introduced inclusion of atomic charge information in the DFT-D4 model robustly leads to quantitatively improved dispersion energies in particular for metallic systems. By construction it yields results which are easily understandable due to a clear separation into hybridization and charge (oxidation) state and two- and many-body effects, respectively. Due to its high computational efficiency, the D4 dispersion model is even applicable to low-cost classical and semiempirical theoretical methods.

Exploring the chemical nature of super-heavy main-group elements by means of efficient plane-wave density-functional theory

Presenting an accurate yet efficient plane-wave DFT approach for the computational exploration of the bulk properties of the super-heavy main-group elements including copernicium (Cn–Og, Z = 112–118).

Aggregation Behavior of a Six-Membered Cyclic Frustrated Phosphane/Borane Lewis Pair: Formation of a Supramolecular Cyclooctameric Macrocyclic Ring System

A new six-membered cyclic frustrated phosphane/borane Lewis pair was liberated from its HB(C₆ F₅ )₂ adduct by treatment with vinylcyclohexane. The system is an active frustrated Lewis pair that undergoes cycloaddition reactions with suitable π reagents and it splits dihydrogen. At room temperature in solution the new compound is a monomer, however, in the crystal and in solution at low temperature it aggregates to a thermodynamically favoured supramolecular macrocyclic cyclooctamer.

High-resolution time and spatial imaging of tobacco and its pyrolysis products during a cigarette puff by microprobe sampling photoionisation mass spectrometry

The time- and space-resolved chemical signatures of gases and vapours formed in solid-state combustion processes are difficult to examine using recent analytical techniques. A machine-smoked cigarette represents a very reproducible model system for dynamic solid-state combustion. By using a special sampling system (microprobe unit) that extracts the formed gases from inside of the burning cigarette, which is coupled to a photoionisation mass spectrometer, it was possible to study the evolution of organic gases during a 2-s cigarette puff. The concentrations of various pyrolysis and combustion products such as 1,3-butadiene, toluene, acetaldehyde and phenol were monitored on-line at different sampling points within cigarettes. A near-microscopic-scale spatial resolution and a 200-ms time resolution were achieved. Finally, the recorded information was combined to generate time-resolved concentration maps, showing the formation and destruction zones of the investigated compounds in the burning cigarette. The combustion zone at the tip of cigarette, where e.g. 1,3-butadiene is predominately formed, was clearly separable from the pyrolysis zones. Depending on the stability of the precursor (e.g. lignin or cellulose), the position of pyrolytic formation varies. In conclusion, it was demonstrated that soft photoionisation mass spectrometry in conjunction with a microprobe sampling device can be used for time- and space-resolved analysis of combustion and pyrolysis reactions. In addition to studies on the model cigarette, further model systems may be studied with this approach. This may include further studies on the combustion of biomass or coal chunks, on heterogeneously catalysed reactions or on spray, dust and gas combustion processes.

The three-dimensional structure of cystathionine beta-lyase from Arabidopsis and its substrate specificity

The pyridoxal 5'-phosphate-dependent enzyme cystathionine beta-lyase (CBL) catalyzes the penultimate step in the de novo biosynthesis of Met in microbes and plants. Absence of CBL in higher organisms makes it an important target for the development of antibiotics and herbicides. The three-dimensional structure of cystathionine beta-lyase from Arabidopsis was determined by Patterson search techniques, using the structure of tobacco (Nicotiana tabacum) cystathionine gamma-synthase as starting point. At a resolution of 2.3 A, the model was refined to a final crystallographic R-factor of 24.9%. The overall structure is very similar to other pyridoxal 5'-phosphate-dependent enzymes of the gamma-family. Exchange of a few critical residues within the active site causes the different substrate preferences between Escherichia coli and Arabidopsis CBL. Loss of interactions at the alpha-carboxyl site is the reason for the poorer substrate binding of Arabidopsis CBL. In addition, the binding pocket of Arabidopsis CBL is larger than that of E. coli CBL, explaining the similar binding of L-cystathionine and L-djenkolate in Arabidopsis CBL in contrast to E. coli CBL, where the substrate binding site is optimized for the natural substrate cystathionine.

Structural investigations of the highly flexible recombinant ribosomal protein L12 from Thermotoga maritima

Ribosomal protein L7/L12, the only multicopy component of the ribosome, is involved in translation factor binding and in the ribosomal GTPase center. The gene for L7/L12 from Thermotoga maritima was cloned and the protein expressed at high levels in Escherichia coli. Purification of L7/L12 was achieved under non-denaturing conditions via heat treatment and two chromatographic steps. Circular dichroism melting profiles were monitored at 222 nm, showing the melting temperature of the protein at pH 7.5 around 110 degrees C, compared to approximately 60 degrees C for the highly homologous Escherichia coli protein. The unfolding was reversible and renaturation closely followed the path of the thermal melting. Dynamic light scattering, gel filtration chromatography, and crosslinking experiments suggested that under physiological buffer conditions Thermotoga maritima L7/L12 exists as a tetramer. The protein was crystallized under two conditions, yielding an orthorhombic (C222(1)) and a cubic (12(1)3) space group with an estimated two and three to four L7/L12 molecules per asymmetric unit, respectively. The crystals contained the full-length protein, and cryogenic buffers were developed which improved the mosaic spreads and the resolution limits. For the structure solution isoleucine was mutated to methionine at two separate positions, the mutant forms expressed as selenomethionine variants and crystallized.

Characterization of recombinant Arabidopsis thaliana threonine synthase

Threonine synthase (TS) catalyses the last step in the biosynthesis of threonine, the pyridoxal 5'-phosphate dependent conversion of L-homoserine phosphate (HSerP) into L-threonine and inorganic phosphate. Recombinant Arabidopsis thaliana TS (aTS) was characterized to compare a higher plant TS with its counterparts from Escherichia coli and yeast. This comparison revealed several unique properties of aTS: (a) aTS is a regulatory enzyme whose activity was increased up to 85-fold by S-adenosyl-L-methionine (SAM) and specifically inhibited by AMP; (b) HSerP analogues shown previously to be potent inhibitors of E. coli TS failed to inhibit aTS; and (c) aTS was a dimer, while the E. coli and yeast enzymes are monomers. The N-terminal region of aTS is essential for its regulatory properties and protects against inhibition by HSerP analogues, as an aTS devoid of 77 N-terminal residues was neither activated by SAM nor inhibited by AMP, but was inhibited by HSerP analogues. The C-terminal region of aTS seems to be involved in dimer formation, as the N-terminally truncated aTS was also found to be a dimer. These conclusions are supported by a multiple amino-acid sequence alignment, which revealed the existence of two TS subfamilies. aTS was classified as a member of subfamily 1 and its N-terminus is at least 35 residues longer than those of any nonplant TS. Monomeric E. coli and yeast TS are members of subfamily 2, characterized by C-termini extending about 50 residues over those of subfamily 1 members. As a first step towards a better understanding of the properties of aTS, the enzyme was crystallized by the sitting drop vapour diffusion method. The crystals diffracted to beyond 0.28 nm resolution and belonged to the space group P222 (unit cell parameters: a = 6.16 nm, b = 10.54 nm, c = 14.63 nm, alpha = beta = gamma = 90 degrees).

Isolation, cloning, sequence analysis and X-ray structure of dimethyl sulfoxide/trimethylamine N-oxide reductase from Rhodobacter capsulatus

The periplasmic enzyme dimethyl sulfoxide/trimethylamine N-oxide reductase (DMSOR/TMAOR) from the photosynthetic purple bacterium Rhodobacter capsulatus functions as the terminal electron acceptor in its respiratory chain. The enzyme catalyzes the reduction of highly oxidized substrates like dimethyl sulfoxide (DMSO) or trimethylamine N-oxide (TMAO). At a molybdenum redox centre, two single electrons are transferred from cytochrome c556 to the substrate, e.g. DMSO, generating dimethyl sulfide (DMS) and water. The operon encoding this enzyme was isolated, cloned and sequenced, and its chromosomal location determined. It was shown by analytical and crystallographic data that DMSOR and TMAOR are identical enzymes. Degenerate primers were derived from short peptide sequences and a 700 bp fragment was amplified by nested PCR, subsequently cloned and radioactively labeled to screen a prepared lambda DASH library. Positive lambda clones were subcloned into pBluescript and subsequently transformed into Escherichia coli to sequence the DMSOR/TMAOR operon. By an optimized protein purification high yields (5 mg protein/l culture) with a specific activity of 30 U/mg were obtained. The molecular mass was experimentally determined by electrospray mass spectroscopy (MS) to be 85034 Da and from the deduced amino acid sequence of the apoenzyme to be 85033 Da. The enzyme was crystallized in space group P4(1)2(1)2 with unit cell dimensions of a = b = 80.7 A and c = 229.2 A diffracting beyond 1.8 A. The three-dimensional structure was solved by a combination of multiple isomorphous replacement (MIR) and molecular replacement techniques. The atomic model was refined to an R-factor of 0.169 for 57394 independent reflections. The spherical protein consists of four domains with a funnel-like cavity that leads to the freely accessible metal-ion redox center. The sole bis(molybdopterin guanine dinucleotide)molybdenum cofactor (1541 Da) of the single chain protein has the molybdenum ion bound to the cis-dithiolene group of only one molybdopterin guanine dinucleotide (MGD) molecule. In addition, two oxo ligands and the oxygen of a serine side chain are bound to the molybdenum ion.

Isolation, cloning, sequence analysis and localization of the operon encoding dimethyl sulfoxide/trimethylamine N-oxide reductase from Rhodobacter capsulatus

The operon encoding the periplasmic enzymes dimethyl sulfoxide reductase (DMSOR) and trimethylamine N-oxide reductase (TMAOR) from the purple, non-sulfur, photosynthetic bacterium Rhodobacter capsulatus was isolated, cloned and sequenced, and its chromosomal location determined. It was shown by analytical and crystallographic data that DMSOR and TMAOR are identical enzymes. Degenerate primers were derived from short peptide sequences generated by automated Edman degradation and a 700 bp fragment was amplified by nested PCR, subsequently cloned and radioactively labeled to screen a prepared lambda DASH library. Positive lambda clones were subcloned into pBluescript and subsequently transformed into Escherichia coli to sequence the DMSOR/ TMAOR operon. The promoter consisted of an A + T-rich region followed by a -35 region, a putative ribosome binding site, and a leader sequence of 13 amino acid residues. The transcription terminator was a G + C-rich dyad sequence capable of forming a hairpin structure, which may act rho-independently. An optimized protein purification of the wild-type enzyme is also described, giving high yields (5 mg protein per liter of culture) and a specific activity of 30 units/mg. The molecular mass was determined by electrospray mass spectrometry to be 85,034 Da; from the deduced amino acid sequence the molecular mass of the apoenzyme was 85,033 Da.