51
|
Matczak P, Mlostoń G, Hamera-Fałdyga R, Görls H, Weigand W. Structure of Diferrocenyl Thioketone: From Molecule to Crystal. MOLECULES (BASEL, SWITZERLAND) 2019; 24:molecules24213950. [PMID: 31683693 PMCID: PMC6864675 DOI: 10.3390/molecules24213950] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/12/2019] [Revised: 10/28/2019] [Accepted: 10/29/2019] [Indexed: 01/07/2023]
Abstract
Ferrocenyl-functionalized thioketones have recently been recognized as useful building blocks for sulfur-containing compounds with potential applications in materials chemistry. This work is devoted to a single representative of such thioketones, namely diferrocenyl thioketone (Fc2CS), whose structure has been determined here for the first time. Both X-ray crystallography and a wide variety of quantum-chemical methods were used to explore the structure of Fc2CS. In addition to the X-ray structure determination, intermolecular interactions occurring in the crystal structure of Fc2CS were examined in detail by quantum-chemical methods. These methods were also an invaluable tool in studying the molecular structure of Fc2CS, from the gas phase to solutions and to its crystal. Intramolecular interactions governing the conformational behavior of an isolated Fc2CS molecule were deduced from quantum-chemical analyses carried out in orbital space and real space. Our experimental and theoretical results indicate that the main structural features of an isolated Fc2CS molecule in its lowest-energy geometry are retained both upon solvation and in the crystal. The tilt of ferrocenyl groups is only slightly affected by crystal packing forces that are dominated by dispersion. Nonetheless, a network of intermolecular interactions, such as H···H, C···H and S···H, was detected in the Fc2CS crystal but each of them is fairly weak.
Collapse
Affiliation(s)
- Piotr Matczak
- Department of Physical Chemistry, Faculty of Chemistry, University of Lodz, Pomorska 163/165, 90236 Lodz, Poland.
| | - Grzegorz Mlostoń
- Department of Organic and Applied Chemistry, Faculty of Chemistry, University of Lodz, Tamka 12, 91403 Lodz, Poland.
| | - Róża Hamera-Fałdyga
- Department of Organic and Applied Chemistry, Faculty of Chemistry, University of Lodz, Tamka 12, 91403 Lodz, Poland.
| | - Helmar Görls
- Institute of Inorganic and Analytical Chemistry, Friedrich-Schiller-University Jena, Humboldtstrasse 8, 07743 Jena, Germany.
| | - Wolfgang Weigand
- Institute of Inorganic and Analytical Chemistry, Friedrich-Schiller-University Jena, Humboldtstrasse 8, 07743 Jena, Germany.
| |
Collapse
|
52
|
Wappett DA, Goerigk L. Toward a Quantum-Chemical Benchmark Set for Enzymatically Catalyzed Reactions: Important Steps and Insights. J Phys Chem A 2019; 123:7057-7074. [DOI: 10.1021/acs.jpca.9b05088] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
| | - Lars Goerigk
- School of Chemistry, The University of Melbourne, Victoria 3010, Australia
| |
Collapse
|
53
|
Grauso L, Teta R, Esposito G, Menna M, Mangoni A. Computational prediction of chiroptical properties in structure elucidation of natural products. Nat Prod Rep 2019; 36:1005-1030. [PMID: 31166350 DOI: 10.1039/c9np00018f] [Citation(s) in RCA: 57] [Impact Index Per Article: 11.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
Covering: up to 2019This review covers the current status of the quantum mechanical prediction of chiroptical properties, such as electronic CD and optical rotation, as needed for stereochemical assignments in new natural products. The reliability of the prediction of chiroptical properties is steadily increasing, with a parallel decrease in the required computational resources. Now, quantum mechanical calculations for a medium-sized natural product can be reliably performed by natural product chemists on a mainstream PC. This review is aimed to guide natural product chemists through the numerous steps involved in such calculations. Through a concise, but comprehensive, discussion of the current computational practice, enriched by a few illustrative examples, this review provides readers with the theoretical background and practical knowledge needed to select the most appropriate parameters for performing the calculations, to anticipate possible problems, and to critically evaluate the reliability of their computational results. Common reasons for mistakes are also discussed; in particular, the importance of the correct evaluation of conformational ensembles of flexible molecules (an aspect often overlooked in current research) is stressed.
Collapse
Affiliation(s)
- Laura Grauso
- Department of Agricultural Sciences, University of Naples Federico II, Via Università 100, 80055 Portici (NA), Italy
| | - Roberta Teta
- Department of Pharmacy, University of Naples Federico II, Via Domenico Montesano 49, 80131 Napoli, Italy.
| | - Germana Esposito
- Department of Pharmacy, University of Naples Federico II, Via Domenico Montesano 49, 80131 Napoli, Italy.
| | - Marialuisa Menna
- Department of Pharmacy, University of Naples Federico II, Via Domenico Montesano 49, 80131 Napoli, Italy.
| | - Alfonso Mangoni
- Department of Pharmacy, University of Naples Federico II, Via Domenico Montesano 49, 80131 Napoli, Italy.
| |
Collapse
|
54
|
Santra G, Sylvetsky N, Martin JML. Minimally Empirical Double-Hybrid Functionals Trained against the GMTKN55 Database: revDSD-PBEP86-D4, revDOD-PBE-D4, and DOD-SCAN-D4. J Phys Chem A 2019; 123:5129-5143. [PMID: 31136709 PMCID: PMC9479158 DOI: 10.1021/acs.jpca.9b03157] [Citation(s) in RCA: 231] [Impact Index Per Article: 46.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
![]()
We
present a family of minimally empirical double-hybrid DFT functionals
parametrized against the very large and diverse GMTKN55 benchmark.
The very recently proposed ωB97M(2) empirical double hybrid
(with 16 adjustable parameters) has the lowest WTMAD2 (weighted mean
absolute deviation over GMTKN55) ever reported at 2.19 kcal/mol. However,
refits of the DSD-BLYP and DSD-PBEP86 spin-component-scaled, dispersion-corrected
double hybrids can achieve WTMAD2 values as low as 2.33 with the very
recent D4 dispersion correction (2.42 kcal/mol with the D3(BJ) dispersion
term) using just a handful of adjustable parameters. If we use full
DFT correlation in the initial orbital evaluation, the xrevDSD-PBEP86-D4
functional reaches WTMAD2 = 2.23 kcal/mol, statistically indistinguishable
from ωB97M(2) but using just four nonarbitrary adjustable parameters
(and three semiarbitrary ones). The changes from the original DSD
parametrizations are primarily due to noncovalent interaction energies
for large systems, which were undersampled in the original parametrization
set. With the new parametrization, same-spin correlation can be eliminated
at minimal cost in performance, which permits revDOD-PBEP86-D4 and
revDOD-PBE-D4 functionals that scale as N4 or even N3 with the size of the system.
Dependence of WTMAD2 for DSD functionals on the percentage of HF exchange
is roughly quadratic; it is sufficiently weak that any reasonable
value in the 64% to 72% range can be chosen semiarbitrarily. DSD-SCAN
and DOD-SCAN double hybrids involving the SCAN nonempirical meta-GGA
as the semilocal component have also been considered and offer a good
alternative if one wishes to eliminate either the empirical dispersion
correction or the same-spin correlation component. noDispSD-SCAN66
achieves WTMAD2 = 3.0 kcal/mol, compared to 2.7 kcal/mol for DOD-SCAN66-D4.
However, the best performance without dispersion corrections (WTMAD2
= 2.8 kcal/mol) is reached by revωB97X-2, a slight reparametrization
of the Chai–Head-Gordon range-separated double hybrid. Finally,
in the context of double-hybrid functionals, the very recent D4 dispersion
correction is clearly superior over D3(BJ).
Collapse
Affiliation(s)
- Golokesh Santra
- Department of Organic Chemistry, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| | - Nitai Sylvetsky
- Department of Organic Chemistry, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| | - Jan M. L. Martin
- Department of Organic Chemistry, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| |
Collapse
|
55
|
Bannwarth C, Ehlert S, Grimme S. GFN2-xTB-An Accurate and Broadly Parametrized Self-Consistent Tight-Binding Quantum Chemical Method with Multipole Electrostatics and Density-Dependent Dispersion Contributions. J Chem Theory Comput 2019; 15:1652-1671. [PMID: 30741547 DOI: 10.1021/acs.jctc.8b01176] [Citation(s) in RCA: 1433] [Impact Index Per Article: 286.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
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.
Collapse
Affiliation(s)
- Christoph Bannwarth
- Mulliken Center for Theoretical Chemistry , Universität Bonn , Beringstr. 4 , 53115 Bonn , Germany.,Department of Chemistry , Stanford University , Stanford , California 94305 , United States
| | - Sebastian Ehlert
- Mulliken Center for Theoretical Chemistry , Universität Bonn , Beringstr. 4 , 53115 Bonn , Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry , Universität Bonn , Beringstr. 4 , 53115 Bonn , Germany
| |
Collapse
|
56
|
PEPCONF, a diverse data set of peptide conformational energies. Sci Data 2019; 6:180310. [PMID: 30667382 PMCID: PMC6343515 DOI: 10.1038/sdata.2018.310] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2018] [Accepted: 11/30/2018] [Indexed: 12/13/2022] Open
Abstract
We present an extensive and diverse database of peptide conformational energies. Our database contains five different classes of model geometries: dipeptides, tripeptides, and disulfide-bridged, bioactive, and cyclic peptides. In total, the database consists of 3775 conformational energy data points and 4530 conformer geometries. All the reference energies have been calculated at the LC-ωPBE-XDM/aug-cc-pVTZ level of theory, which is shown to yield conformational energies with an accuracy in the order of tenths of a kcal/mol when compared to complete-basis-set coupled-cluster reference data. The peptide conformational data set (PEPCONF) is presented as a high-quality reference set for the development and benchmarking of molecular-mechanics and semi-empirical electronic structure methods, which are the most commonly used techniques in the modeling of medium to large proteins.
Collapse
|
57
|
Perera AS, Cheramy J, Poopari MR, Xu Y. Aggregation of lactic acid in cold rare-gas matrices and the link to solution: a matrix isolation-vibrational circular dichroism study. Phys Chem Chem Phys 2019; 21:3574-3584. [DOI: 10.1039/c8cp04748k] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Crucial insight into lactic acid self-aggregation in solution is obtained by following its unique VCD spectral features in cold matrices.
Collapse
Affiliation(s)
| | - Joseph Cheramy
- Chemistry Department
- The University of Alberta
- Edmonton
- Canada
| | | | - Yunjie Xu
- Chemistry Department
- The University of Alberta
- Edmonton
- Canada
| |
Collapse
|
58
|
Morgante P, Peverati R. Statistically representative databases for density functional theory via data science. Phys Chem Chem Phys 2019; 21:19092-19103. [DOI: 10.1039/c9cp03211h] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Cluster analysis applied to quantum chemistry: a new broad database of chemical properties with a reasonable computational cost.
Collapse
|
59
|
Karton A. Thermochemistry of Guanine Tautomers Re-Examined by Means of High-Level CCSD(T) Composite Ab Initio Methods. Aust J Chem 2019. [DOI: 10.1071/ch19276] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
Abstract
We obtained accurate gas-phase tautomerization energies for a set of 14 guanine tautomers by means of high-level thermochemical procedures approximating the CCSD(T) energy at the complete basis set (CBS) limit. For the five low-lying tautomers, we use the computationally demanding W1-F12 composite method for obtaining the tautomerization energies. The relative W1-F12 tautomerization enthalpies at 298K are: 0.00 (1), 2.37 (2), 2.63 (3), 4.03 (3′), and 14.31 (4) kJmol−1. Thus, as many as four tautomers are found within a small energy window of less than 1.0kcalmol−1 (1kcalmol−1=4.184kJmol−1). We use these highly accurate W1-F12 tautomerization energies to evaluate the performance of a wide range of lower-level composite ab initio procedures. The Gn composite procedures (G4, G4(MP2), G4(MP2)-6X, G3, G3B3, G3(MP2), and G3(MP2)B3) predict that the enol tautomer (3) is more stable than the keto tautomer (2) by amounts ranging from 0.36 (G4) to 1.28 (G3(MP2)) kJmol−1. We also find that an approximated CCSD(T)/CBS energy calculated as HF/jul-cc-pV{D,T}Z+CCSD/jul-cc-pVTZ+(T)/jul-cc-pVDZ results in a root-mean-square deviation (RMSD) of merely 0.11kJmol−1 relative to the W1-F12 reference values. We use this approximated CCSD(T)/CBS method to obtain the tautomerization energies of 14 guanine tautomers. The relative tautomerization enthalpies at 298K are: 0.00 (1), 2.20 (2), 2.51 (3), 4.06 (3′), 14.30 (4), 25.65 (5), 43.78 (4′), 53.50 (6′), 61.58 (6), 77.37 (7), 82.52 (8′), 86.02 (9), 100.70 (10), and 121.01 (8) kJmol−1. Using these tautomerization enthalpies, we evaluate the performance of standard and composite methods for the entire set of 14 guanine tautomers. The best-performing procedures emerge as (RMSDs are given in parentheses): G4(MP2)-6X (0.51), CCSD(T)+ΔMP2/CBS (0.52), and G4(MP2) (0.64kJmol−1). The worst performers are CCSD(T)/AVDZ (1.05), CBS-QB3 (1.24), and CBS-APNO (1.38kJmol−1).
Collapse
|
60
|
Morgante P, Peverati R. ACCDB: A collection of chemistry databases for broad computational purposes. J Comput Chem 2018; 40:839-848. [DOI: 10.1002/jcc.25761] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/05/2018] [Revised: 10/09/2018] [Accepted: 11/11/2018] [Indexed: 01/04/2023]
Affiliation(s)
- Pierpaolo Morgante
- Chemistry Program; Florida Institute of Technology, 150 W. University Blvd.; Melbourne Florida, 32901
| | - Roberto Peverati
- Chemistry Program; Florida Institute of Technology, 150 W. University Blvd.; Melbourne Florida, 32901
| |
Collapse
|
61
|
Kesharwani MK, Sylvetsky N, Köhn A, Tew DP, Martin JML. Do CCSD and approximate CCSD-F12 variants converge to the same basis set limits? The case of atomization energies. J Chem Phys 2018; 149:154109. [PMID: 30342453 DOI: 10.1063/1.5048665] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
While the title question is a clear "yes" from purely theoretical arguments, the case is less clear for practical calculations with finite (one-particle) basis sets. To shed further light on this issue, the convergence to the basis set limit of CCSD (coupled cluster theory with all single and double excitations) and of different approximate implementations of CCSD-F12 (explicitly correlated CCSD) has been investigated in detail for the W4-17 thermochemical benchmark. Near the CBS ([1-particle] complete basis set) limit, CCSD and CCSD(F12*) agree to within their respective uncertainties (about ±0.04 kcal/mol) due to residual basis set incompleteness error, but a nontrivial difference remains between CCSD-F12b and CCSD(F12*), which is roughly proportional to the degree of static correlation. The observed basis set convergence behavior results from the superposition of a rapidly converging, attractive, CCSD[F12]-CCSD-F12b difference (consisting mostly of third-order terms) and a more slowly converging, repulsive, fourth-order difference between CCSD(F12*) and CCSD[F12]. For accurate thermochemistry, we recommend CCSD(F12*) over CCSD-F12b if at all possible. There are some indications that the nZaPa family of basis sets exhibits somewhat smoother convergence than the correlation consistent family.
Collapse
Affiliation(s)
- Manoj K Kesharwani
- Department of Organic Chemistry, Weizmann Institute of Science, 76100 Reḥovot, Israel
| | - Nitai Sylvetsky
- Department of Organic Chemistry, Weizmann Institute of Science, 76100 Reḥovot, Israel
| | - Andreas Köhn
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - David P Tew
- Max Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Jan M L Martin
- Department of Organic Chemistry, Weizmann Institute of Science, 76100 Reḥovot, Israel
| |
Collapse
|
62
|
Thomas J, Carrillo MJ, Serrato A, Xie F, Jäger W, Xu Y, Lin W. Microwave spectrum of the complex of 3,3,3-trifluoro-2-(trifluoromethyl)propanoic acid and formic acid. Mol Phys 2018. [DOI: 10.1080/00268976.2018.1532539] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Affiliation(s)
- Javix Thomas
- Department of Chemistry, University of Alberta, Edmonton, Alberta, Canada
| | - Michael J. Carrillo
- Department of Chemistry, The University of Texas Rio Grande Valley, Brownsville, US
| | - Agapito Serrato
- Department of Chemistry, The University of Texas Rio Grande Valley, Brownsville, US
| | - Fan Xie
- Department of Chemistry, University of Alberta, Edmonton, Alberta, Canada
| | - Wolfgang Jäger
- Department of Chemistry, University of Alberta, Edmonton, Alberta, Canada
| | - Yunjie Xu
- Department of Chemistry, University of Alberta, Edmonton, Alberta, Canada
| | - Wei Lin
- Department of Chemistry, The University of Texas Rio Grande Valley, Brownsville, US
| |
Collapse
|
63
|
Najibi A, Goerigk L. The Nonlocal Kernel in van der Waals Density Functionals as an Additive Correction: An Extensive Analysis with Special Emphasis on the B97M-V and ωB97M-V Approaches. J Chem Theory Comput 2018; 14:5725-5738. [DOI: 10.1021/acs.jctc.8b00842] [Citation(s) in RCA: 108] [Impact Index Per Article: 18.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Asim Najibi
- School of Chemistry, The University of Melbourne, Parkville, Victoria 3010, Australia
| | - Lars Goerigk
- School of Chemistry, The University of Melbourne, Parkville, Victoria 3010, Australia
| |
Collapse
|
64
|
Loboda OA, Dolgonos GA, Boese AD. Towards hybrid density functional calculations of molecular crystals via fragment-based methods. J Chem Phys 2018; 149:124104. [DOI: 10.1063/1.5046908] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Oleksandr A. Loboda
- Institute of Chemistry, University of Graz, Heinrichstrasse 28/IV, A-8010 Graz, Austria
| | - Grygoriy A. Dolgonos
- Institute of Chemistry, University of Graz, Heinrichstrasse 28/IV, A-8010 Graz, Austria
| | - A. Daniel Boese
- Institute of Chemistry, University of Graz, Heinrichstrasse 28/IV, A-8010 Graz, Austria
| |
Collapse
|
65
|
Mehta N, Casanova-Páez M, Goerigk L. Semi-empirical or non-empirical double-hybrid density functionals: which are more robust? Phys Chem Chem Phys 2018; 20:23175-23194. [PMID: 30062343 DOI: 10.1039/c8cp03852j] [Citation(s) in RCA: 89] [Impact Index Per Article: 14.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The development of non-empirical double-hybrid density functionals (DHDFs) is a very active research area with the number of approaches in this field having increased rapidly. At the same time, there is a lack of published work that provides a fair assessment and comparison between non-empirical and semi-empirical DHDFs on an equal footing. Herein, we close this gap and present a thorough analysis of both classes of DHDFs on the large GMTKN55 benchmark database for general main-group thermochemistry, kinetics, and noncovalent interactions [Goerigk et al., Phys. Chem. Chem. Phys., 2017, 19, 32184-32215]. In total, 115 variations of dispersion-corrected and -uncorrected DHDFs are tested, which will be condensed to an in-depth assessment of 31 methods: 19 non-empirical and 12 semi-empirical DHDFs. As such, our study represents the largest DHDF study ever conducted and can serve as an important benchmark informing method developers and users alike. Our results show that semi-empirical DHDFs are the most robust density functional approximations and more reliable and accurate than non-empirical ones. In fact, some non-empirical approaches are even outperformed by hybrid approaches or even dispersion-corrected and -uncorrected MP2 and SCS-MP2. SOS0-PBE0-2-D3(BJ) is the only exception and the only non-empirical DHDF that we can safely recommend for general applicability. However, it is still outperformed by six semi-empirical DHDFs, of which we would like to particularly recommend the following five: ωB97X-2-D3(BJ), DSD-BLYP-D3(BJ), DSD-PBEP86-D3(BJ), B2NC-PLYP-D3(BJ), and B2GPPLYP-D3(BJ). Our findings seriously question current trends in the field and they highlight that novel strategies have to be found in order to outperform the currently best density functional theory methods on the market. We hope that our study can function as an important cornerstone inspiring such a change of direction in the field.
Collapse
Affiliation(s)
- Nisha Mehta
- School of Chemistry, The University of Melbourne, Parkville, Australia.
| | | | | |
Collapse
|
66
|
Mardirossian N, Head-Gordon M. Survival of the most transferable at the top of Jacob's ladder: Defining and testing the ωB97M(2) double hybrid density functional. J Chem Phys 2018; 148:241736. [PMID: 29960332 PMCID: PMC5991970 DOI: 10.1063/1.5025226] [Citation(s) in RCA: 117] [Impact Index Per Article: 19.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/08/2018] [Accepted: 05/09/2018] [Indexed: 01/19/2023] Open
Abstract
A meta-generalized gradient approximation, range-separated double hybrid (DH) density functional with VV10 non-local correlation is presented. The final 14-parameter functional form is determined by screening trillions of candidate fits through a combination of best subset selection, forward stepwise selection, and random sample consensus (RANSAC) outlier detection. The MGCDB84 database of 4986 data points is employed in this work, containing a training set of 870 data points, a validation set of 2964 data points, and a test set of 1152 data points. Following an xDH approach, orbitals from the ωB97M-V density functional are used to compute the second-order perturbation theory correction. The resulting functional, ωB97M(2), is benchmarked against a variety of leading double hybrid density functionals, including B2PLYP-D3(BJ), B2GPPLYP-D3(BJ), ωB97X-2(TQZ), XYG3, PTPSS-D3(0), XYGJ-OS, DSD-PBEP86-D3(BJ), and DSD-PBEPBE-D3(BJ). Encouragingly, the overall performance of ωB97M(2) on nearly 5000 data points clearly surpasses that of all of the tested density functionals. As a Rung 5 density functional, ωB97M(2) completes our family of combinatorially optimized functionals, complementing B97M-V on Rung 3, and ωB97X-V and ωB97M-V on Rung 4. The results suggest that ωB97M(2) has the potential to serve as a powerful predictive tool for accurate and efficient electronic structure calculations of main-group chemistry.
Collapse
Affiliation(s)
- Narbe Mardirossian
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| |
Collapse
|
67
|
Exploring conformational preferences of alanine tetrapeptide by CCSD(T), MP2, and dispersion-corrected DFT methods. Chem Phys Lett 2018. [DOI: 10.1016/j.cplett.2018.05.006] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
68
|
Kraus P, Frank I. Density Functional Theory for Microwave Spectroscopy of Noncovalent Complexes: A Benchmark Study. J Phys Chem A 2018; 122:4894-4901. [PMID: 29750513 DOI: 10.1021/acs.jpca.8b03345] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
In this work, we compare the results obtained with 89 computational methods for predicting noncovalent bond lengths in weakly bound complexes. Evaluations for the performance in noncovalent interaction energies and covalent bond lengths obtained from five other data sets are included. The overall best performing density functional is the ωB97M-V method, achieving balanced results across all three categories. For noncovalent geometries, the best methods include B97M-V, B3LYP-D3(BJ) and DSD-PBEPBE-D3(BJ). The effects of systematic improvement of the density functional approximation and of dispersion corrections are also discussed.
Collapse
Affiliation(s)
- P Kraus
- Institut für Physikalische Chemie und Elektrochemie , Leibniz Universität Hannover , Callinstraße 3A , 30167 Hannover , Germany
| | - I Frank
- Institut für Physikalische Chemie und Elektrochemie , Leibniz Universität Hannover , Callinstraße 3A , 30167 Hannover , Germany
| |
Collapse
|
69
|
Cavasin AT, Hillisch A, Uellendahl F, Schneckener S, Göller AH. Reliable and Performant Identification of Low-Energy Conformers in the Gas Phase and Water. J Chem Inf Model 2018; 58:1005-1020. [PMID: 29717870 DOI: 10.1021/acs.jcim.8b00151] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
Prediction of compound properties from structure via quantitative structure-activity relationship and machine-learning approaches is an important computational chemistry task in small-molecule drug research. Though many such properties are dependent on three-dimensional structures or even conformer ensembles, the majority of models are based on descriptors derived from two-dimensional structures. Here we present results from a thorough benchmark study of force field, semiempirical, and density functional methods for the calculation of conformer energies in the gas phase and water solvation as a foundation for the correct identification of relevant low-energy conformers. We find that the tight-binding ansatz GFN-xTB shows the lowest error metrics and highest correlation to the benchmark PBE0-D3(BJ)/def2-TZVP in the gas phase for the computationally fast methods and that in solvent OPLS3 becomes comparable in performance. MMFF94, AM1, and DFTB+ perform worse, whereas the performance-optimized but far more expensive functional PBEh-3c yields energies almost perfectly correlated to the benchmark and should be used whenever affordable. On the basis of our findings, we have implemented a reliable and fast protocol for the identification of low-energy conformers of drug-like molecules in water that can be used for the quantification of strain energy and entropy contributions to target binding as well as for the derivation of conformer-ensemble-dependent molecular descriptors.
Collapse
Affiliation(s)
| | - Alexander Hillisch
- Bayer AG , Drug Discovery, Chemical Research , 42096 Wuppertal , Germany
| | - Felix Uellendahl
- Bayer AG , Drug Discovery, Chemical Research , 42096 Wuppertal , Germany
| | - Sebastian Schneckener
- Bayer AG , Engineering & Technology, Applied Mathematics , 51368 Leverkusen , Germany
| | - Andreas H Göller
- Bayer AG , Drug Discovery, Chemical Research , 42096 Wuppertal , Germany
| |
Collapse
|
70
|
Řezáč J, Bím D, Gutten O, Rulíšek L. Toward Accurate Conformational Energies of Smaller Peptides and Medium-Sized Macrocycles: MPCONF196 Benchmark Energy Data Set. J Chem Theory Comput 2018; 14:1254-1266. [DOI: 10.1021/acs.jctc.7b01074] [Citation(s) in RCA: 55] [Impact Index Per Article: 9.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Affiliation(s)
- Jan Řezáč
- Gilead Sciences Research Center and The Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic
| | - Daniel Bím
- Gilead Sciences Research Center and The Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic
| | - Ondrej Gutten
- Gilead Sciences Research Center and The Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic
| | - Lubomír Rulíšek
- Gilead Sciences Research Center and The Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic
| |
Collapse
|
71
|
Abdo YA, Weeks JW, Layfield W, Tremlett WM, Graham JW, Tabor ME, Causey SE, Carr JM, Tschumper GS. Intramolecular Hydrogen Bonding in α-Epoxy Alcohols: A Conformational Analysis of 1,2-Dialkyl-2,3-epoxycyclopentanol Diastereomers. CHEM LETT 2018. [DOI: 10.1246/cl.170932] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Affiliation(s)
- Yasmeen A. Abdo
- Department of Chemistry and Biochemistry, The University of Mississippi, 322 Coulter Hall, University, Mississippi, 38677-1848, USA
| | - Justin W. Weeks
- Chemistry Department, Central Alabama Community College, 1675 Cherokee Road, Alexander, Alabama 35010, USA
| | - William Layfield
- Chemistry Department, Central Alabama Community College, 1675 Cherokee Road, Alexander, Alabama 35010, USA
| | - William M. Tremlett
- Chemistry Department, Central Alabama Community College, 1675 Cherokee Road, Alexander, Alabama 35010, USA
| | - Joshua W. Graham
- Chemistry Department, Central Alabama Community College, 1675 Cherokee Road, Alexander, Alabama 35010, USA
| | - Mikayla E. Tabor
- Chemistry Department, Central Alabama Community College, 1675 Cherokee Road, Alexander, Alabama 35010, USA
| | - Sarah E. Causey
- Chemistry Department, Central Alabama Community College, 1675 Cherokee Road, Alexander, Alabama 35010, USA
| | - Jeremy M. Carr
- Chemistry Department, Central Alabama Community College, 1675 Cherokee Road, Alexander, Alabama 35010, USA
| | - Gregory S. Tschumper
- Department of Chemistry and Biochemistry, The University of Mississippi, 322 Coulter Hall, University, Mississippi, 38677-1848, USA
| |
Collapse
|
72
|
Kesharwani MK, Karton A, Sylvetsky N, Martin JML. The S66 Non-Covalent Interactions Benchmark Reconsidered Using Explicitly Correlated Methods Near the Basis Set Limit. Aust J Chem 2018. [DOI: 10.1071/ch17588] [Citation(s) in RCA: 33] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
The S66 benchmark for non-covalent interactions has been re-evaluated using explicitly correlated methods with basis sets near the one-particle basis set limit. It is found that post-MP2 ‘high-level corrections’ are treated adequately well using a combination of CCSD(F12*) with (aug-)cc-pVTZ-F12 basis sets on the one hand, and (T) extrapolated from conventional CCSD(T)/heavy-aug-cc-pV{D,T}Z on the other hand. Implications for earlier benchmarks on the larger S66×8 problem set in particular, and for accurate calculations on non-covalent interactions in general, are discussed. At a slight cost in accuracy, (T) can be considerably accelerated by using sano-V{D,T}Z+ basis sets, whereas half-counterpoise CCSD(F12*)(T)/cc-pVDZ-F12 offers the best compromise between accuracy and computational cost.
Collapse
|
73
|
Chan B, Radom L. Modelling the Effect of Conformation on Hydrogen-Atom Abstraction from Peptides. Aust J Chem 2018. [DOI: 10.1071/ch17621] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Abstract
Computational quantum chemistry is used to examine the effect of conformation on the kinetics of hydrogen-atom abstraction by HO• from amides of glycine and proline as peptide models. In accord with previous findings, it is found that there are substantial variations possible in the conformations and the corresponding energies, with the captodative effect, hydrogen bonding, and solvation being some of the major features that contribute to the variations. The ‘minimum-energy-structure-pathway’ strategy that is often employed in theoretical studies of peptide chemistry with small models certainly provides valuable fundamental information. However, one may anticipate different reaction outcomes in structurally constrained systems due to modified reaction thermodynamics and kinetics, as demonstrated explicitly in the present study. Thus, using a ‘consistent-conformation-pathway’ approach may indeed be more informative in such circumstances, and in this regard theory provides information that would be difficult to obtain from experimental studies alone.
Collapse
|
74
|
Song C, Martínez TJ. Analytical gradients for tensor hyper-contracted MP2 and SOS-MP2 on graphical processing units. J Chem Phys 2017; 147:161723. [DOI: 10.1063/1.4997997] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Chenchen Song
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, USA and SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
| | - Todd J. Martínez
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, USA and SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
| |
Collapse
|
75
|
Vamhindi BSDR, Karton A. Can DFT and ab initio methods adequately describe binding energies in strongly interacting C6X6⋯C2X π–π complexes? Chem Phys 2017. [DOI: 10.1016/j.chemphys.2017.05.020] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
76
|
Witte J, Neaton JB, Head-Gordon M. Effective empirical corrections for basis set superposition error in the def2-SVPD basis: gCP and DFT-C. J Chem Phys 2017. [DOI: 10.1063/1.4986962] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023] Open
|
77
|
Mardirossian N, Head-Gordon M. Thirty years of density functional theory in computational chemistry: an overview and extensive assessment of 200 density functionals. Mol Phys 2017. [DOI: 10.1080/00268976.2017.1333644] [Citation(s) in RCA: 709] [Impact Index Per Article: 101.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Affiliation(s)
- Narbe Mardirossian
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA, USA
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA, USA
| |
Collapse
|
78
|
Chehelamirani M, Salahub DR. Effect of dispersion corrections on covalent and non-covalent interactions in DFTB calculations. Struct Chem 2017. [DOI: 10.1007/s11224-017-0976-1] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
79
|
Manna D, Kesharwani MK, Sylvetsky N, Martin JML. Conventional and Explicitly Correlated ab Initio Benchmark Study on Water Clusters: Revision of the BEGDB and WATER27 Data Sets. J Chem Theory Comput 2017; 13:3136-3152. [DOI: 10.1021/acs.jctc.6b01046] [Citation(s) in RCA: 69] [Impact Index Per Article: 9.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Affiliation(s)
- Debashree Manna
- Department of Organic Chemistry, Weizmann Institute of Science, 76100 Rechovot, Israel
| | - Manoj K. Kesharwani
- Department of Organic Chemistry, Weizmann Institute of Science, 76100 Rechovot, Israel
| | - Nitai Sylvetsky
- Department of Organic Chemistry, Weizmann Institute of Science, 76100 Rechovot, Israel
| | - Jan M. L. Martin
- Department of Organic Chemistry, Weizmann Institute of Science, 76100 Rechovot, Israel
| |
Collapse
|
80
|
Witte J, Mardirossian N, Neaton JB, Head-Gordon M. Assessing DFT-D3 Damping Functions Across Widely Used Density Functionals: Can We Do Better? J Chem Theory Comput 2017; 13:2043-2052. [DOI: 10.1021/acs.jctc.7b00176] [Citation(s) in RCA: 56] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
| | | | - Jeffrey B. Neaton
- Kavli Energy
Nanosciences
Institute at Berkeley, Berkeley, California 94720, United States
| | | |
Collapse
|
81
|
Song C, Martínez TJ. Atomic orbital-based SOS-MP2 with tensor hypercontraction. II. Local tensor hypercontraction. J Chem Phys 2017; 146:034104. [DOI: 10.1063/1.4973840] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Chenchen Song
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, USA and SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
| | - Todd J. Martínez
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, USA and SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
| |
Collapse
|
82
|
Sancho-García JC, Brémond É, Savarese M, Pérez-Jiménez AJ, Adamo C. Partnering dispersion corrections with modern parameter-free double-hybrid density functionals. Phys Chem Chem Phys 2017; 19:13481-13487. [DOI: 10.1039/c7cp00709d] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
The PBE-QIDH and SOS1-PBE-QIDH double-hybrid density functionals are merged with a pair of dispersion corrections, namely the pairwise additive D3(BJ) and the non-local correlation functional VV10, leading to the corresponding dispersion-corrected models.
Collapse
Affiliation(s)
| | - É. Brémond
- CompuNet
- Istituto Italiano di Tecnologia
- I-16163 Genoa
- Italy
| | - M. Savarese
- CompuNet
- Istituto Italiano di Tecnologia
- I-16163 Genoa
- Italy
| | | | - C. Adamo
- CompuNet
- Istituto Italiano di Tecnologia
- I-16163 Genoa
- Italy
- Institut de Recherche de Chimie Paris
| |
Collapse
|
83
|
Goerigk L, Hansen A, Bauer C, Ehrlich S, Najibi A, Grimme S. A look at the density functional theory zoo with the advanced GMTKN55 database for general main group thermochemistry, kinetics and noncovalent interactions. Phys Chem Chem Phys 2017; 19:32184-32215. [DOI: 10.1039/c7cp04913g] [Citation(s) in RCA: 854] [Impact Index Per Article: 122.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
We present the updated and extended GMTKN55 benchmark database for more accurate and extensive energetic evaluation of density functionals and other electronic structure methods with detailed guidelines for method users.
Collapse
Affiliation(s)
- Lars Goerigk
- School of Chemistry
- The University of Melbourne
- Parkville
- Australia
| | - Andreas Hansen
- Universität Bonn
- Mulliken Center for Theoretical Chemistry
- Bonn
- Germany
| | - Christoph Bauer
- Universität Bonn
- Mulliken Center for Theoretical Chemistry
- Bonn
- Germany
| | - Stephan Ehrlich
- Universität Bonn
- Mulliken Center for Theoretical Chemistry
- Bonn
- Germany
| | - Asim Najibi
- School of Chemistry
- The University of Melbourne
- Parkville
- Australia
| | - Stefan Grimme
- Universität Bonn
- Mulliken Center for Theoretical Chemistry
- Bonn
- Germany
| |
Collapse
|
84
|
Mardirossian N, Head-Gordon M. How Accurate Are the Minnesota Density Functionals for Noncovalent Interactions, Isomerization Energies, Thermochemistry, and Barrier Heights Involving Molecules Composed of Main-Group Elements? J Chem Theory Comput 2016; 12:4303-25. [PMID: 27537680 DOI: 10.1021/acs.jctc.6b00637] [Citation(s) in RCA: 286] [Impact Index Per Article: 35.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Abstract
The 14 Minnesota density functionals published between the years 2005 and early 2016 are benchmarked on a comprehensive database of 4986 data points (84 data sets) involving molecules composed of main-group elements. The database includes noncovalent interactions, isomerization energies, thermochemistry, and barrier heights, as well as equilibrium bond lengths and equilibrium binding energies of noncovalent dimers. Additionally, the sensitivity of the Minnesota density functionals to the choice of basis set and integration grid is explored for both noncovalent interactions and thermochemistry. Overall, the main strength of the hybrid Minnesota density functionals is that the best ones provide very good performance for thermochemistry (e.g., M06-2X), barrier heights (e.g., M08-HX, M08-SO, MN15), and systems heavily characterized by self-interaction error (e.g., M06-2X, M08-HX, M08-SO, MN15), while the main weakness is that none of them are state-of-the-art for the full spectrum of noncovalent interactions and isomerization energies (although M06-2X is recommended from the 10 hybrid Minnesota functionals). Similarly, the main strength of the local Minnesota density functionals is that the best ones provide very good performance for thermochemistry (e.g., MN15-L), barrier heights (e.g., MN12-L), and systems heavily characterized by self-interaction error (e.g., MN12-L and MN15-L), while the main weakness is that none of them are state-of-the-art for the full spectrum of noncovalent interactions and isomerization energies (although M06-L is clearly the best from the four local Minnesota functionals). As an overall guide, M06-2X and MN15 are perhaps the most broadly useful hybrid Minnesota functionals, while M06-L and MN15-L are perhaps the most broadly useful local Minnesota functionals, although each has different strengths and weaknesses.
Collapse
Affiliation(s)
- Narbe Mardirossian
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California , Berkeley, California 94720, United States.,Chemical Sciences Division, Lawrence Berkeley National Laboratory , Berkeley, California 94720, United States
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California , Berkeley, California 94720, United States.,Chemical Sciences Division, Lawrence Berkeley National Laboratory , Berkeley, California 94720, United States
| |
Collapse
|
85
|
Mardirossian N, Head-Gordon M. ωB97M-V: A combinatorially optimized, range-separated hybrid, meta-GGA density functional with VV10 nonlocal correlation. J Chem Phys 2016; 144:214110. [DOI: 10.1063/1.4952647] [Citation(s) in RCA: 370] [Impact Index Per Article: 46.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023] Open
Affiliation(s)
- Narbe Mardirossian
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
86
|
First-principles data set of 45,892 isolated and cation-coordinated conformers of 20 proteinogenic amino acids. Sci Data 2016; 3:160009. [PMID: 26881946 PMCID: PMC4755128 DOI: 10.1038/sdata.2016.9] [Citation(s) in RCA: 41] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2015] [Accepted: 01/15/2016] [Indexed: 11/24/2022] Open
Abstract
We present a structural data set of the 20 proteinogenic amino acids and their amino-methylated and acetylated (capped) dipeptides. Different protonation states of the backbone (uncharged and zwitterionic) were considered for the amino acids as well as varied side chain protonation states. Furthermore, we studied amino acids and dipeptides in complex with divalent cations (Ca2+, Ba2+, Sr2+, Cd2+, Pb2+, and Hg2+). The database covers the conformational hierarchies of 280 systems in a wide relative energy range of up to 4 eV (390 kJ/mol), summing up to a total of 45,892 stationary points on the respective potential-energy surfaces. All systems were calculated on equal first-principles footing, applying density-functional theory in the generalized gradient approximation corrected for long-range van der Waals interactions. We show good agreement to available experimental data for gas-phase ion affinities. Our curated data can be utilized, for example, for a wide comparison across chemical space of the building blocks of life, for the parametrization of protein force fields, and for the calculation of reference spectra for biophysical applications.
Collapse
|
87
|
Brauer B, Kesharwani MK, Kozuch S, Martin JML. The S66x8 benchmark for noncovalent interactions revisited: explicitly correlated ab initio methods and density functional theory. Phys Chem Chem Phys 2016; 18:20905-25. [DOI: 10.1039/c6cp00688d] [Citation(s) in RCA: 137] [Impact Index Per Article: 17.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
The S66x8 dataset for noncovalent interactions of biochemical relevance has been re-examined by means of CCSD(F12*)(T), DFT, and SAPT methods.
Collapse
Affiliation(s)
- Brina Brauer
- Department of Organic Chemistry
- Weizmann Institute of Science
- 76100 Rehovot
- Israel
| | - Manoj K. Kesharwani
- Department of Organic Chemistry
- Weizmann Institute of Science
- 76100 Rehovot
- Israel
| | - Sebastian Kozuch
- Department of Chemistry
- Ben-Gurion University of the Negev
- 84105 Beer-Sheva
- Israel
| | - Jan M. L. Martin
- Department of Organic Chemistry
- Weizmann Institute of Science
- 76100 Rehovot
- Israel
| |
Collapse
|
88
|
Brandenburg JG, Caldeweyher E, Grimme S. Screened exchange hybrid density functional for accurate and efficient structures and interaction energies. Phys Chem Chem Phys 2016; 18:15519-23. [DOI: 10.1039/c6cp01697a] [Citation(s) in RCA: 42] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
HSE-3c: a computationally efficient and numerically robust screened hybrid functional that can be applied to periodic small gap systems.
Collapse
Affiliation(s)
| | - Eike Caldeweyher
- Mulliken Center for Theoretical Chemistry
- University of Bonn
- 53115 Bonn
- Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry
- University of Bonn
- 53115 Bonn
- Germany
| |
Collapse
|