The SQM/COSMO filter: reliable native pose identification based on the quantum-mechanical description of protein–ligand interactions and implicit COSMO solvation

Current virtual screening tools are fast, but reliable scoring is elusive. Here, we present the "SQM/COSMO filter", a novel scoring function featuring a quantitative semiempirical quantum mechanical (SQM) description of all types of noncovalent interactions coupled with implicit COSMO solvation. We show unequivocally that it outperforms eight widely used scoring functions. The accuracy and chemical generality of the SQM/COSMO filter make it a perfect tool for late stages of virtual screening

Efficient Neutron Production from a Novel Configuration of Deuterium Gas-Puff Z-Pinch

A novel configuration of a deuterium z pinch has been used to generate fusion neutrons. Injecting an outer hollow cylindrical plasma shell around an inner deuterium gas puff, neutron yields from DD reactions reached Y-n = (2.9 +/- 0.3) x 10(12) at 700 ns implosion time and 2.7 MA current. Such a neutron yield means a tenfold increase in comparison with previous deuterium gas puff experiments at the same current generator. The increase of beam-target yields was obtained by a larger amount of current assembled on the z-pinch axis, and subsequently by higher induced voltage and higher energies of deuterons. A stack of CR-39 track detectors on the z-pinch axis showed hydrogen ions up to 38 MeV. Maximum neutron energies of 15 and 22 MeV were observed by radial and axial time-of-flight detectors, respectively. The number of DD neutrons per one joule of stored plasma energy approached 5 x 10(7). This implies that deuterium gas puff z pinches belong to the most efficient plasma-based sources of DD neutrons

Noncovalent Interactions by QMC: Speedup by One-Particle Basis-Set Size Reduction

While it is empirically accepted that the fixed-node diffusion Monte-Carlo (FN-DMC) depends only weakly on the size of the one-particle basis sets used to expand its guiding functions, limits of this observation are not settled yet. Our recent work indicates that under the FN error cancellation conditions, augmented triple zeta basis sets are sufficient to achieve a benchmark level of 0.1 kcal/mol in a number of small noncovalent complexes. Here we report on a possibility of truncation of the one-particle basis sets used in FN-DMC guiding functions that has no visible effect on the accuracy of the production FN-DMC energy differences. The proposed scheme leads to no significant increase in the local energy variance, indicating that the total CPU cost of large-scale benchmark noncovalent interaction energy FN-DMC calculations may be reduced.Comment: ACS book chapter, accepte

S66: A Well-balanced Database of Benchmark Interaction Energies Relevant to Biomolecular Structures

With numerous new quantum chemistry methods being developed in recent years and the promise of even more new methods to be developed in the near future, it is clearly critical that highly accurate, well-balanced, reference data for many different atomic and molecular properties be available for the parametrization and validation of these methods. One area of research that is of particular importance in many areas of chemistry, biology, and material science is the study of noncovalent interactions. Because these interactions are often strongly influenced by correlation effects, it is necessary to use computationally expensive high-order wave function methods to describe them accurately. Here, we present a large new database of interaction energies calculated using an accurate CCSD(T)/CBS scheme. Data are presented for 66 molecular complexes, at their reference equilibrium geometries and at 8 points systematically exploring their dissociation curves; in total, the database contains 594 points: 66 at equilibrium geometries, and 528 in dissociation curves. The data set is designed to cover the most common types of noncovalent interactions in biomolecules, while keeping a balanced representation of dispersion and electrostatic contributions. The data set is therefore well suited for testing and development of methods applicable to bioorganic systems. In addition to the benchmark CCSD(T) results, we also provide decompositions of the interaction energies by means of DFT-SAPT calculations. The data set was used to test several correlated QM methods, including those parametrized specifically for noncovalent interactions. Among these, the SCS-MI-CCSD method outperforms all other tested methods, with a root-mean-square error of 0.08 kcal/mol for the S66 data set

Accurate Treatment of Large Supramolecular Complexes by Double-Hybrid Density Functionals Coupled with Nonlocal van der Waals Corrections

In this work, we present a thorough assessment of the performance of some representative double-hybrid density functionals (revPBE0-DH-NL and B2PLYP-NL) as well as their parent hybrid and GGA counterparts, in combination with the most modern version of the nonlocal (NL) van der Waals correction to describe very large weakly interacting molecular systems dominated by noncovalent interactions. Prior to the assessment, an accurate and homogeneous set of reference interaction energies was computed for the supramolecular complexes constituting the L7 and S12L data sets by using the novel, precise, and efficient DLPNO-CCSD(T) method at the complete basis set limit (CBS). The correction of the basis set superposition error and the inclusion of the deformation energies (for the S12L set) have been crucial for obtaining precise DLPNO-CCSD(T)/CBS interaction energies. Among the density functionals evaluated, the double-hybrid revPBE0-DH-NL and B2PLYP-NL with the three-body dispersion correction provide remarkably accurate association energies very close to the chemical accuracy. Overall, the NL van der Waals approach combined with proper density functionals can be seen as an accurate and affordable computational tool for the modeling of large weakly bonded supramolecular systems.Financial support by the “Ministerio de Economía y Competitividad” (MINECO) of Spain and European FEDER funds through projects CTQ2011-27253 and CTQ2012-31914 is acknowledged. The support of the Generalitat Valenciana (Prometeo/2012/053) is also acknowledged. J.A. thanks the EU for the FP7-PEOPLE-2012-IEF-329513 grant. J.C. acknowledges the “Ministerio de Educación, Cultura y Deporte” (MECD) of Spain for a predoctoral FPU grant

Unraveling the performance of dispersion-corrected functionals for the accurate description of weakly bound natural polyphenols

Long-range non-covalent interactions play a key role in the chemistry of natural polyphenols. We have previously proposed a description of supramolecular polyphenol complexes by the B3P86 density functional coupled with some corrections for dispersion. We couple here the B3P86 functional with the D3 correction for dispersion, assessing systematically the accuracy of the new B3P86-D3 model using for that the well-known S66, HB23, NCCE31, and S12L datasets for non-covalent interactions. Furthermore, the association energies of these complexes were carefully compared to those obtained by other dispersion-corrected functionals, such as B(3)LYP-D3, BP86-D3 or B3P86-NL. Finally, this set of models were also applied to a database composed of seven non-covalent polyphenol complexes of the most interest.FDM acknowledges financial support from the Swedish Research Council (Grant No. 621-2014-4646) and SNIC (Swedish National Infrastructure for Computing) for providing computer resources. The work in Limoges (IB and PT) is supported by the “Conseil Régional du Limousin”. PT gratefully acknowledges the support by the Operational Program Research and Development Fund (project CZ.1.05/2.1.00/03.0058 of the Ministry of Education, Youth and Sports of the Czech Republic). IB gratefully acknowledges financial support from “Association Djerbienne en France”

van der Waals Interactions in Material Modelling

Van der Waals (vdW) interactions stem from electronic zero-point fluctuations and are often critical for the correct description of structure, stability, and response properties of molecules and materials, including biomolecules, nanomaterials, and material interfaces. Here, we give a conceptual as well as mathematical overview of the current state of modeling vdW interactions,focusing in particular on the consequences of different approximations for practical applications. We present a systematic classification of approximate first-principles models based on the adiabatic-connection fluctuation-dissipation theorem, namely the nonlocal density functionals, interatomic methods, and methods based on the random-phase approximation. The applicability of these methods to different types of materials and material properties is discussed in connection with availability of theoretical and experimental benchmarks. We conclude with a roadmap of the open problems that remain to be solved to construct a universal, efficient, and accurate vdW model for realistic material modeling