An efficient and near linear scaling pair natural orbital based local coupled cluster method

Going clean: structure and dynamics of peptides in the gas phase and paths to solvation

The gas phase is an artificial environment for biomolecules that has gained much attention both experimentally and theoretically due to its unique characteristic of providing a clean room environment for the comparison between theory and experiment. In this review we give an overview mainly on first-principles simulations of isolated peptides and the initial steps of their interactions with ions and solvent molecules: a bottom up approach to the complexity of biological environments. We focus on the accuracy of different methods to explore the conformational space, the connections between theory and experiment regarding collision cross section evaluations and (anharmonic) vibrational spectra, and the challenges faced in this field.

First UHF Implementation of the Incremental Scheme for Open-Shell Systems

The incremental scheme makes it possible to compute CCSD(T) correlation energies to high accuracy for large systems. We present the first extension of this fully automated black-box approach to open-shell systems using an Unrestricted Hartree-Fock (UHF) wave function, extending the efficient domain-specific basis set approach to handle open-shell references. We test our approach on a set of organic and metal organic structures and molecular clusters and demonstrate standard deviations from canonical CCSD(T) values of only 1.35 kJ/mol using a triple ζ basis set. We find that the incremental scheme is significantly more cost-effective than the canonical implementation even for relatively small systems and that the ease of parallelization makes it possible to perform high-level calculations on large systems in a few hours on inexpensive computers. We show that the approximations that make our approach widely applicable are significantly smaller than both the basis set incompleteness error and the intrinsic error of the CCSD(T) method, and we further demonstrate that incremental energies can be reliably used in extrapolation schemes to obtain near complete basis set limit CCSD(T) reaction energies for large systems.

Reactivity of Dimeric Tetrazirconium(IV) Wells-Dawson Polyoxometalate toward Dipeptide Hydrolysis Studied by a Combined Experimental and Density Functional Theory Approach

Detailed kinetic studies on the hydrolysis of glycylglycine (Gly-Gly) in the presence of the dimeric tetrazirconium(IV)-substituted Wells-Dawson-type polyoxometalate Na14[Zr4(P2W16O59)2(μ3-O)2(OH)2(H2O)4] · 57H2O (1) were performed by a combination of (1)H, (13)C, and (31)P NMR spectroscopies. The catalyst was shown to be stable under a broad range of reaction conditions. The effect of pD on the hydrolysis of Gly-Gly showed a bell-shaped profile with the fastest hydrolysis observed at pD 7.4. The observed rate constant for the hydrolysis of Gly-Gly at pD 7.4 and 60 °C was 4.67 × 10(-7) s(-1), representing a significant acceleration as compared to the uncatalyzed reaction. (13)C NMR data were indicative for coordination of Gly-Gly to 1 via its amide oxygen and amine nitrogen atoms, resulting in a hydrolytically active complex. Importantly, the effective hydrolysis of a series of Gly-X dipeptides with different X side chain amino acids in the presence of 1 was achieved, and the observed rate constant was shown to be dependent on the volume, chemical nature, and charge of the X amino acid side chain. To give a mechanistic explanation of the observed catalytic hydrolysis of Gly-Gly, a detailed quantum-chemical study was performed. The theoretical results confirmed the nature of the experimentally suggested binding mode in the hydrolytically active complex formed between Gly-Gly and 1. To elucidate the role of 1 in the hydrolytic process, both the uncatalyzed and the polyoxometalate-catalyzed reactions were examined. In the rate-determining step of the uncatalyzed Gly-Gly hydrolysis, a carboxylic oxygen atom abstracts a proton from a solvent water molecule and the nascent OH nucleophile attacks the peptide carbon atom. Analogous general-base activity of the free carboxylic group was found to take place also in the case of polyoxometalate-catalyzed hydrolysis as the main catalytic effect originates from the -C═O···Zr(IV) binding.

Mechanistic Aspects of Aryl-Halide Oxidative Addition, Coordination Chemistry, and Ring-Walking by Palladium

This contribution describes the reactivity of a zero-valent palladium phosphine complex with substrates that contain both an aryl halide moiety and an unsaturated carbon-carbon bond. Although η(2) -coordination of the metal center to a C=C or C≡C unit is kinetically favored, aryl halide bond activation is favored thermodynamically. These quantitative transformations proceed under mild reaction conditions in solution or in the solid state. Kinetic measurements indicate that formation of η(2) -coordination complexes are not nonproductive side-equilibria, but observable (and in several cases even isolated) intermediates en route to aryl halide bond cleavage. At the same time, DFT calculations show that the reaction with palladium may proceed through a dissociation-oxidative addition mechanism rather than through a haptotropic intramolecular process (i.e., ring walking). Furthermore, the transition state involves coordination of a third phosphine to the palladium center, which is lost during the oxidative addition as the C-halide bond is being broken. Interestingly, selective activation of aryl halides has been demonstrated by adding reactive aryl halides to the η(2) -coordination complexes. The product distribution can be controlled by the concentration of the reactants and/or the presence of excess phosphine.

Enhancing NMR Prediction for Organic Compounds Using Molecular Dynamics

NMR spectroscopy is a crucial tool in organic chemistry for the routine characterization of small molecules, structural elucidation of natural products, and study of reaction mechanisms. Although there is evidence that thermal motions strongly affect observed resonances, conventional predictions are performed only on stationary structures. Here we show that quasiclassical molecular dynamics provides a highly accurate and broadly applicable method for improving shielding predictions. Gas-phase values of the absolute shieldings of protons and carbons are predicted to nearly within experimental uncertainty, while the chemical shifts of large systems such as natural products are closely reproduced. Importantly, these results are obtained without the use of any empirical corrections. Our analysis suggests that the linear scaling factors currently employed are primarily a correction for vibrational effects. As a result, our method extends the reach of prediction methods to the study of molecules with unusual dynamics such as the iconic and controversial [18]annulene. Our predictions agree closely with experiment at both low and high temperatures and provide strong evidence that the equilibrium structure of [18]annulene is planar and aromatic.

Quantum chemical benchmark study on 46 RNA backbone families using a dinucleotide unit

We have created a benchmark set of quantum chemical structure-energy data denoted as UpU46, which consists of 46 uracil dinucleotides (UpU), representing all known 46 RNA backbone conformational families. Penalty-function-based restrained optimizations with COSMO TPSS-D3/def2-TZVP ensure a balance between keeping the target conformation and geometry relaxation. The backbone geometries are close to the clustering-means of their respective RNA bioinformatics family classification. High-level wave function methods (DLPNO-CCSD(T) as reference) and a wide-range of dispersion-corrected or inclusive DFT methods (DFT-D3, VV10, LC-BOP-LRD, M06-2X, M11, and more) are used to evaluate the conformational energies. The results are compared to the Amber RNA bsc0χOL3 force field. Most dispersion-corrected DFT methods surpass the Amber force field significantly in accuracy and yield mean absolute deviations (MADs) for relative conformational energies of ∼0.4-0.6 kcal/mol. Double-hybrid density functionals represent the most accurate class of density functionals. Low-cost quantum chemical methods such as PM6-D3H+, HF-3c, DFTB3-D3, as well as small basis set calculations corrected for basis set superposition errors (BSSEs) by the gCP procedure are also tested. Unfortunately, the presently available low-cost methods are struggling to describe the UpU conformational energies with satisfactory accuracy. The UpU46 benchmark is an ideal test for benchmarking and development of fast methods to describe nucleic acids, including force fields.

Control in the Rate-Determining Step Provides a Promising Strategy To Develop New Catalysts for CO2 Hydrogenation: A Local Pair Natural Orbital Coupled Cluster Theory Study

The development of efficient catalysts with base metals for CO2 hydrogenation has always been a major thrust of interest. A series of experimental and theoretical work has revealed that the catalytic cycle typically involves two key steps, namely, base-promoted heterolytic H2 splitting and hydride transfer to CO2, either of which can be the rate-determining step (RDS) of the entire reaction. To explore the determining factor for the nature of RDS, we present herein a comparative mechanistic investigation on CO2 hydrogenation mediated by [M(H)(η(2)-H2)(PP3(Ph))](n+) (M = Fe(II), Ru(II), and Co(III); PP3(Ph) = tris(2-(diphenylphosphino)phenyl)phosphine) type complexes. In order to construct reliable free energy profiles, we used highly correlated wave function based ab initio methods of the coupled cluster type alongside the standard density functional theory. Our calculations demonstrate that the hydricity of the metal-hydride intermediate generated by H2 splitting dictates the nature of the RDS for the Fe(II) and Co(III) systems, while the RDS for the Ru(II) catalyst appears to be ambiguous. CO2 hydrogenation catalyzed by the Fe(II) complex that possesses moderate hydricity traverses an H2-splitting RDS, whereas the RDS for the high-hydricity Co(III) species is found to be the hydride transfer. Thus, our findings suggest that hydricity can be used as a practical guide in future catalyst design. Enhancing the electron-accepting ability of low-hydricity catalysts is likely to improve their catalytic performance, while increasing the electron-donating ability of high-hydricity complexes may speed up CO2 conversion. Moreover, we also established the active roles of base NEt3 in directing the heterolytic H2 splitting and assisting product release through the formation of an acid-base complex.

Computational study of a model system of enzyme-mediated [4+2] cycloaddition reaction

A possible mechanistic pathway related to an enzyme-catalyzed [4+2] cycloaddition reac-tion was studied by theoretical calculations at density functional (B3LYP, O3LYP, M062X) and semiempirical levels (PM6-DH2, PM6) performed on a model system. The calculations were carried out for the key [4+2] cycloaddition step considering enzyme-catalyzed biosynthesis of Spinosyn A in a model reaction, where a reliable example of a biological Diels-Alder reaction was reported experimentally. In the present study it was demonstrated that the [4+2] cycloaddition reaction may benefit from moving along the energetically balanced reaction coordinate, which enabled the catalytic rate enhancement of the [4+2] cycloaddition pathway involving a single transition state. Modeling of such a system with coordination of three amino acids indicated a reliable decrease of activation energy by ~18.0 kcal/mol as compared to a non-catalytic transformation.

Noncovalent interactions determine the conformation of aurophilic complexes with 2-mercapto-4-methyl-5-thiazoleacetic acid ligands

The synthesis of three new gold(i) salts is reported. The intermolecular interactions stabilizing the aurophilic dimers are examined by several first principles interpretative techniques.

Incremental evaluation of coupled cluster dipole polarizabilities

In this work we present the first implementation of the incremental scheme for coupled cluster linear-response frequency-dependent dipole polarizabilities.

Towards biochemically relevant QM computations on nucleic acids: controlled electronic structure geometry optimization of nucleic acid structural motifs using penalty restraint functions

Recent developments in dispersion-corrected density functional theory methods allow for the first time the description of large fragments of nucleic acids (hundreds of atoms) with an accuracy clearly surpassing the accuracy of common biomolecular force fields. Such calculations can significantly improve the description of the potential energy surface of nucleic acid molecules, which may be useful for studies of molecular interactions and conformational preferences of nucleic acids, as well as verification and parameterization of other methods. The first of such studies, however, demonstrated that successful applications of accurate QM calculations to larger nucleic acid building blocks are hampered by difficulties in obtaining geometries that are biochemically relevant and are not biased by non-native structural features. We present an approach that can greatly facilitate large-scale QM studies on nucleic acids, namely electronic structure geometry optimization of nucleic acid fragments utilizing a penalty function to restrain key internal coordinates with a specific focus on the torsional backbone angles. This work explores the viability of these restraint optimizations for DFT-D3, PM6-D3H and HF-3c optimizations on a set of examples (a UpA dinucleotide, a DNA G-quadruplex and a B-DNA fragment). Evaluation of different penalty function strengths reveals only a minor system-dependency and reasonable restraint values range from 0.01 to 0.05 Eh rad(-2) for the backbone torsions. Restraints are crucial to perform the QM calculations on biochemically relevant conformations in implicit solvation and gas phase geometry optimizations. The reasons for using restrained instead of constrained or unconstrained optimizations are explained and an open-source external optimizer is provided.

Quantum chemical approach to estimating the thermodynamics of metabolic reactions

Thermodynamics plays an increasingly important role in modeling and engineering metabolism. We present the first nonempirical computational method for estimating standard Gibbs reaction energies of metabolic reactions based on quantum chemistry, which can help fill in the gaps in the existing thermodynamic data. When applied to a test set of reactions from core metabolism, the quantum chemical approach is comparable in accuracy to group contribution methods for isomerization and group transfer reactions and for reactions not including multiply charged anions. The errors in standard Gibbs reaction energy estimates are correlated with the charges of the participating molecules. The quantum chemical approach is amenable to systematic improvements and holds potential for providing thermodynamic data for all of metabolism.

The thermochemistry of london dispersion-driven transition metal reactions: getting the 'right answer for the right reason'

Reliable thermochemical measurements and theoretical predictions for reactions involving large transition metal complexes in which long-range intramolecular London dispersion interactions contribute significantly to their stabilization are still a challenge, particularly for reactions in solution. As an illustrative and chemically important example, two reactions are investigated where a large dipalladium complex is quenched by bulky phosphane ligands (triphenylphosphane and tricyclohexylphosphane). Reaction enthalpies and Gibbs free energies were measured by isotherm titration calorimetry (ITC) and theoretically 'back-corrected' to yield 0 K gas-phase reaction energies (ΔE). It is shown that the Gibbs free solvation energy calculated with continuum models represents the largest source of error in theoretical thermochemistry protocols. The ('back-corrected') experimental reaction energies were used to benchmark (dispersion-corrected) density functional and wave function theory methods. Particularly, we investigated whether the atom-pairwise D3 dispersion correction is also accurate for transition metal chemistry, and how accurately recently developed local coupled-cluster methods describe the important long-range electron correlation contributions. Both, modern dispersion-corrected density functions (e.g., PW6B95-D3(BJ) or B3LYP-NL), as well as the now possible DLPNO-CCSD(T) calculations, are within the 'experimental' gas phase reference value. The remaining uncertainties of 2-3 kcal mol(-1) can be essentially attributed to the solvation models. Hence, the future for accurate theoretical thermochemistry of large transition metal reactions in solution is very promising.

Formation and structure of two luminescent salts of [Au(SCSN3)2]- obtained through the [2 + 3] cyclization of carbon disulfide and azide ion

The colorless, two-coordinate gold(i) complex, [(H2O)3Na][Au(SCSN3)2], has been synthesized through the [2 + 3] cyclic reaction of carbon disulfide and sodium azide in the presence of the labile complex (tht)AuCl. Metathesis of [(H2O)3Na][Au(SCSN3)2], with tetra(phenyl)arsonium chloride produced colorless needles of (Ph4As)[Au(SCSN3)2]. The structure of [(H2O)3Na][Au(SCSN3)2] involves linear gold coordination by two exocyclic sulfur atoms of the 1,2,3,4-thiatriazole-5-thiolate anions. These two-coordinate anions self-associate to form extended, zig-zag chains that are connected by aurophilic bonding with Au∙∙∙Au distances of 3.2653(3) Å and 3.3090(3) Å. Remarkably, the individual S-Au-S units that are connected though aurophilic interactions are eclipsed. The structure of (Ph4As)[Au(SCSN3)2] also contains linear, two-coordinate gold ions with bonding to the 1,2,3,4-thiatriazole-5-thiolate anionic ligands through the exocyclic sulfur atoms. However, in this salt, the anions self-associate through AuAu bonds (Au∙∙∙Au distance of 3.2007(3) Å) to form simple dimers, which also have an eclipsed arrangement of the ligands. Electronic structure calculations strongly suggest that the staggered geometry for the [(Au(SCSN3))2](+) dimer is energetically favored relative to the eclipsed geometry. However, attractive π-stacking interactions appear to promote the observed eclipsed arrangement of the ligands.

Calculations of Mode-Specific Tunneling of Double-Hydrogen Transfer in Porphycene Agree with and Illuminate Experiment

We report a theoretical study of mode-specific tunneling splittings in double-hydrogen transfer in trans-porphycene. We use a novel, mode-specific "Qim path method", in which the reaction coordinate is the imaginary-frequency normal mode of the saddle point separating the equivalent minima. The model considers all 108 normal modes and uses no adjustable parameters. The method gives the ground vibrational-state tunneling splitting, as well the increase in the splitting upon excitation of certain modes, in good agreement with experiment. Interpretation of these results is also transparent with this method. In addition, predictions are made for mode excitations not investigated experimentally. Results for d1 and d2 isotopolgues are also in agreement with experiment.

Improved cluster-in-molecule local correlation approach for electron correlation calculation of large systems

An improved cluster-in-molecule (CIM) local correlation approach is developed to allow electron correlation calculations of large systems more accurate and faster. We have proposed a refined strategy of constructing virtual LMOs of various clusters, which is suitable for basis sets of various types. To recover medium-range electron correlation, which is important for quantitative descriptions of large systems, we find that a larger distance threshold (ξ) is necessary for highly accurate results. Our illustrative calculations show that the present CIM-MP2 (second-order Møller-Plesser perturbation theory, MP2) or CIM-CCSD (coupled cluster singles and doubles, CCSD) scheme with a suitable ξ value is capable of recovering more than 99.8% correlation energies for a wide range of systems at different basis sets. Furthermore, the present CIM-MP2 scheme can provide reliable relative energy differences as the conventional MP2 method for secondary structures of polypeptides.

Tractability gains in symmetry-adapted perturbation theory including coupled double excitations: CCD+ST(CCD) dispersion with natural orbital truncations

This work focuses on efficient and accurate treatment of the intermolecular dispersion interaction using the CCD+ST(CCD) dispersion approach formulated by Williams et al. [J. Chem. Phys. 103, 4586 (1995)]. We apply natural orbital truncation techniques to the solution of the monomer coupled-cluster double (CCD) equations, yielding substantial accelerations in this computationally demanding portion of the SAPT2+(CCD), SAPT2+(3)(CCD), and SAPT2+3(CCD) analyses. It is shown that the wholly rate-limiting dimer-basis particle-particle ladder term can be computed in a reduced natural virtual space which is essentially the same size as the monomer-basis virtual space, with an error on the order of a few thousandths of 1 kcal mol(-1). Coupled with our existing natural orbital techniques for the perturbative triple excitation contributions [E. G. Hohenstein and C. D. Sherrill, J. Chem. Phys. 133, 104107 (2010)], this technique provides speedups of greater than an order of magnitude for the evaluation of the complete SAPT2+3(CCD) decomposition, with a total error of a few hundredths of 1 kcal mol(-1). The combined approach yields tractability gains of almost 2× in the system size, allowing for SAPT2+3(CCD)/aug-cc-pVTZ analysis to be performed for systems such as adenine-thymine for the first time. Natural orbital based SAPT2+3(CCD)/aug-cc-pVTZ results are presented for stacked and hydrogen-bonded configurations of uracil dimer and the adenine-thymine dimer.