1
|
Kříž K, Schmidt L, Andersson AT, Walz MM, van der Spoel D. An Imbalance in the Force: The Need for Standardized Benchmarks for Molecular Simulation. J Chem Inf Model 2023; 63:412-431. [PMID: 36630710 PMCID: PMC9875315 DOI: 10.1021/acs.jcim.2c01127] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2022] [Indexed: 01/12/2023]
Abstract
Force fields (FFs) for molecular simulation have been under development for more than half a century. As with any predictive model, rigorous testing and comparisons of models critically depends on the availability of standardized data sets and benchmarks. While such benchmarks are rather common in the fields of quantum chemistry, this is not the case for empirical FFs. That is, few benchmarks are reused to evaluate FFs, and development teams rather use their own training and test sets. Here we present an overview of currently available tests and benchmarks for computational chemistry, focusing on organic compounds, including halogens and common ions, as FFs for these are the most common ones. We argue that many of the benchmark data sets from quantum chemistry can in fact be reused for evaluating FFs, but new gas phase data is still needed for compounds containing phosphorus and sulfur in different valence states. In addition, more nonequilibrium interaction energies and forces, as well as molecular properties such as electrostatic potentials around compounds, would be beneficial. For the condensed phases there is a large body of experimental data available, and tools to utilize these data in an automated fashion are under development. If FF developers, as well as researchers in artificial intelligence, would adopt a number of these data sets, it would become easier to compare the relative strengths and weaknesses of different models and to, eventually, restore the balance in the force.
Collapse
Affiliation(s)
- Kristian Kříž
- Department
of Cell and Molecular Biology, Uppsala University, Box 596, SE-75124Uppsala, Sweden
| | - Lisa Schmidt
- Faculty
of Biosciences, University of Heidelberg, Heidelberg69117, Germany
| | - Alfred T. Andersson
- Department
of Cell and Molecular Biology, Uppsala University, Box 596, SE-75124Uppsala, Sweden
| | - Marie-Madeleine Walz
- Department
of Cell and Molecular Biology, Uppsala University, Box 596, SE-75124Uppsala, Sweden
| | - David van der Spoel
- Department
of Cell and Molecular Biology, Uppsala University, Box 596, SE-75124Uppsala, Sweden
| |
Collapse
|
2
|
Kojasoy V, Tantillo DJ. Impacts of noncovalent interactions involving sulfur atoms on protein stability, structure, folding, and bioactivity. Org Biomol Chem 2022; 21:11-23. [PMID: 36345987 DOI: 10.1039/d2ob01602h] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
This review discusses the various types of noncovalent interactions in which sulfur atoms participate and their effects on protein stability, structure, folding and bioactivity. Current approaches and recommendations for modelling these noncovalent interactions (in terms of both geometries and interaction energies) are highlighted.
Collapse
Affiliation(s)
- Volga Kojasoy
- Department of Chemistry, University of California, Davis, 1 Shields Avenue, Davis, CA, 95616, USA.
| | - Dean J Tantillo
- Department of Chemistry, University of California, Davis, 1 Shields Avenue, Davis, CA, 95616, USA.
| |
Collapse
|
3
|
Prasad VK, Otero-de-la-Roza A, DiLabio GA. Small-Basis Set Density-Functional Theory Methods Corrected with Atom-Centered Potentials. J Chem Theory Comput 2022; 18:2913-2930. [PMID: 35412817 DOI: 10.1021/acs.jctc.2c00036] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Density functional theory (DFT) is currently the most popular method for modeling noncovalent interactions and thermochemistry. The accurate calculation of noncovalent interaction energies, reaction energies, and barrier heights requires choosing an appropriate functional and, typically, a relatively large basis set. Deficiencies of the density-functional approximation and the use of a limited basis set are the leading sources of error in the calculation of noncovalent and thermochemical properties in molecular systems. In this article, we present three new DFT methods based on the BLYP, M06-2X, and CAM-B3LYP functionals in combination with the 6-31G* basis set and corrected with atom-centered potentials (ACPs). ACPs are one-electron potentials that have the same form as effective-core potentials, except they do not replace any electrons. The ACPs developed in this work are used to generate energy corrections to the underlying DFT/basis-set method such that the errors in predicted chemical properties are minimized while maintaining the low computational cost of the parent methods. ACPs were developed for the elements H, B, C, N, O, F, Si, P, S, and Cl. The ACP parameters were determined using an extensive training set of 118655 data points, mostly of complete basis set coupled-cluster level quality. The target molecular properties for the ACP-corrected methods include noncovalent interaction energies, molecular conformational energies, reaction energies, barrier heights, and bond separation energies. The ACPs were tested first on the training set and then on a validation set of 42567 additional data points. We show that the ACP-corrected methods can predict the target molecular properties with accuracy close to complete basis set wavefunction theory methods, but at a computational cost of double-ζ DFT methods. This makes the new BLYP/6-31G*-ACP, M06-2X/6-31G*-ACP, and CAM-B3LYP/6-31G*-ACP methods uniquely suited to the calculation of noncovalent, thermochemical, and kinetic properties in large molecular systems.
Collapse
Affiliation(s)
- Viki Kumar Prasad
- Department of Chemistry, University of British Columbia, Okanagan, 3247 University Way, Kelowna, British Columbia V1V 1V7, Canada
| | - Alberto Otero-de-la-Roza
- Departamento de Química Física y Analítica, Facultad de Química, Universidad de Oviedo, MALTA Consolider Team, Oviedo E-33006, Spain
| | - Gino A DiLabio
- Department of Chemistry, University of British Columbia, Okanagan, 3247 University Way, Kelowna, British Columbia V1V 1V7, Canada
| |
Collapse
|
4
|
Prasad VK, Otero-de-la-Roza A, DiLabio GA. Fast and Accurate Quantum Mechanical Modeling of Large Molecular Systems Using Small Basis Set Hartree-Fock Methods Corrected with Atom-Centered Potentials. J Chem Theory Comput 2022; 18:2208-2232. [PMID: 35313106 DOI: 10.1021/acs.jctc.1c01128] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
There has been significant interest in developing fast and accurate quantum mechanical methods for modeling large molecular systems. In this work, by utilizing a machine learning regression technique, we have developed new low-cost quantum mechanical approaches to model large molecular systems. The developed approaches rely on using one-electron Gaussian-type functions called atom-centered potentials (ACPs) to correct for the basis set incompleteness and the lack of correlation effects in the underlying minimal or small basis set Hartree-Fock (HF) methods. In particular, ACPs are proposed for ten elements common in organic and bioorganic chemistry (H, B, C, N, O, F, Si, P, S, and Cl) and four different base methods: two minimal basis sets (MINIs and MINIX) plus a double-ζ basis set (6-31G*) in combination with dispersion-corrected HF (HF-D3/MINIs, HF-D3/MINIX, HF-D3/6-31G*) and the HF-3c method. The new ACPs are trained on a very large set (73 832 data points) of noncovalent properties (interaction and conformational energies) and validated additionally on a set of 32 048 data points. All reference data are of complete basis set coupled-cluster quality, mostly CCSD(T)/CBS. The proposed ACP-corrected methods are shown to give errors in the tenths of a kcal/mol range for noncovalent interaction energies and up to 2 kcal/mol for molecular conformational energies. More importantly, the average errors are similar in the training and validation sets, confirming the robustness and applicability of these methods outside the boundaries of the training set. In addition, the performance of the new ACP-corrected methods is similar to complete basis set density functional theory (DFT) but at a cost that is orders of magnitude lower, and the proposed ACPs can be used in any computational chemistry program that supports effective-core potentials without modification. It is also shown that ACPs improve the description of covalent and noncovalent bond geometries of the underlying methods and that the improvement brought about by the application of the ACPs is directly related to the number of atoms to which they are applied, allowing the treatment of systems containing some atoms for which ACPs are not available. Overall, the ACP-corrected methods proposed in this work constitute an alternative accurate, economical, and reliable quantum mechanical approach to describe the geometries, interaction energies, and conformational energies of systems with hundreds to thousands of atoms.
Collapse
Affiliation(s)
- Viki Kumar Prasad
- Department of Chemistry, University of British Columbia, Okanagan, 3247 University Way, Kelowna, British Columbia, Canada V1V 1V7
| | - Alberto Otero-de-la-Roza
- MALTA Consolider Team, Departamento de Química Física y Analítica, Facultad de Química, Universidad de Oviedo, E-33006 Oviedo, Spain
| | - Gino A DiLabio
- Department of Chemistry, University of British Columbia, Okanagan, 3247 University Way, Kelowna, British Columbia, Canada V1V 1V7
| |
Collapse
|
5
|
Cuc NTT, An NT, Ngan VT, Chandra AK, Trung NT. Importance of water and intramolecular interaction governs substantial blue shift of Csp2–H stretching frequency in complexes between chalcogenoaldehydes and water. RSC Adv 2022; 12:1998-2008. [PMID: 35425273 PMCID: PMC8979115 DOI: 10.1039/d1ra07444j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2021] [Accepted: 01/04/2022] [Indexed: 11/25/2022] Open
Abstract
Geometrical structure, stability and cooperativity, and contribution of hydrogen bonds to the stability of complexes between chalcogenoaldehydes and water were thoroughly investigated using quantum chemical methods. The stability of the complexes increases significantly when one or more H2O molecules are added to the binary system, whereas it decreases sharply going from O to S, Se, or Te substitution. The O–H⋯O H-bond is twice as stable as Csp2–H⋯O and O–H⋯S/Se/Te H-bonds. It is found that a considerable blue-shift of Csp2–H stretching frequency in the Csp2–H⋯O H-bond is mainly determined by an addition of water into the complexes along with the low polarity of the Csp2–H covalent bond in formaldehyde and acetaldehyde. The Csp2–H stretching frequency shift as a function of net second hyperconjugative energy for the σ*(Csp2–H) antibonding orbital is observed. Remarkably, a considerable Csp2–H blue shift of 109 cm−1 has been reported for the first time. Upon the addition of H2O into the binary systems, halogenated complexes witness a decreasing magnitude of the Csp2–H stretching frequency blue-shift in the Csp2–H⋯O H-bond, whereas CH3-substituted complexes experience the opposite trend. The considerable blue shift of Csp2–H stretching frequency.![]()
Collapse
Affiliation(s)
- Nguyen Thi Thanh Cuc
- Laboratory of Computational Chemistry and Modelling (LCCM), Department of Chemistry, Faculty of Natural Sciences, Quy Nhon University, Vietnam
| | - Nguyen Truong An
- Laboratory of Computational Chemistry and Modelling (LCCM), Department of Chemistry, Faculty of Natural Sciences, Quy Nhon University, Vietnam
| | - Vu Thi Ngan
- Laboratory of Computational Chemistry and Modelling (LCCM), Department of Chemistry, Faculty of Natural Sciences, Quy Nhon University, Vietnam
| | - Asit. K. Chandra
- Department of Chemistry, North-Eastern Hill University, Shillong 793022, Meghalaya, India
| | - Nguyen Tien Trung
- Laboratory of Computational Chemistry and Modelling (LCCM), Department of Chemistry, Faculty of Natural Sciences, Quy Nhon University, Vietnam
| |
Collapse
|
6
|
Cuc NTT, Phan CTD, Nhung NTA, Nguyen MT, Trung NT, Ngan VT. Theoretical Aspects of Nonconventional Hydrogen Bonds in the Complexes of Aldehydes and Hydrogen Chalcogenides. J Phys Chem A 2021; 125:10291-10302. [PMID: 34818019 DOI: 10.1021/acs.jpca.1c06708] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Hydrogen bonds (H-bonds) in the complexes between aldehydes and hydrogen chalcogenides, XCHO...nH2Z with X = H, F, Cl, Br, and CH3, Z = O, S, Se, and Te, and n = 1,2, were investigated using high-level ab initio calculations. The Csp2-H...O H-bonds are found to be about twice as strong as the Csp2-H...S/Se/Te counterparts. Remarkably, the S/Se/Te-H...S/Se/Te H-bonds are 4.5 times as weak as the O-H...O ones. The addition of the second H2Z molecule into binary systems induces stronger complexes and causes a positive cooperative effect in ternary complexes. The blue shift of Csp2-H stretching frequency involving the Csp2-H...Z H-bond sharply increases when replacing one H atom in HCHO by a CH3 group. In contrast, when one H atom in HCHO is substituted with a halogen, the magnitude of blue-shifting of the Csp2-H...Z H-bond becomes smaller. The largest blue shift up to 92 cm-1 of Csp2-H stretching frequency in Csp2-H...O H-bond in CH3CHO...2H2O has rarely been observed and is much greater than that in the cases of the Csp2-H...S/Se/Te ones. The Csp2-H blue shift of Csp2-H...Z bonds in the halogenated aldehydes is converted into a red shift when H2O is replaced by a heavier analogue, such as H2S, H2Se, or H2Te. The stability and classification of nonconventional H-bonds including Csp2-H...Se/Te, Te-H...Te, and Se/Te-H...O have been established for the first time.
Collapse
Affiliation(s)
- Nguyen Thi Thanh Cuc
- Laboratory of Computational Chemistry and Modelling (LCCM), Faculty of Natural Sciences, Quy Nhon University, Quy Nhon 55100, Vietnam
| | - Cam-Tu Dang Phan
- Laboratory of Computational Chemistry and Modelling (LCCM), Faculty of Natural Sciences, Quy Nhon University, Quy Nhon 55100, Vietnam
| | - Nguyen Thi Ai Nhung
- Department of Chemistry, University of Sciences, Hue University, Hue 49000, Vietnam
| | | | - Nguyen Tien Trung
- Laboratory of Computational Chemistry and Modelling (LCCM), Faculty of Natural Sciences, Quy Nhon University, Quy Nhon 55100, Vietnam
| | - Vu Thi Ngan
- Laboratory of Computational Chemistry and Modelling (LCCM), Faculty of Natural Sciences, Quy Nhon University, Quy Nhon 55100, Vietnam
| |
Collapse
|
7
|
Ma X, Zhu Y, Yu J, Zhao G, Duanmu J, Yuan Y, Chang XP, Xu D, Zhou Q. Unprecedented observation and characterization of sulfur-centred bifurcated hydrogen bonds. Phys Chem Chem Phys 2021; 23:26519-26523. [PMID: 34807205 DOI: 10.1039/d1cp04601b] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Owing to the small electronegativity of the sulfur atom, it is commonly supposed that at most one weak H-bond can be formed between a sulfur atom and an H-bond donor. In this paper, an unprecedented 2 : 1 binding species generated from two molecules of phenol and a molecule of thioether was observed and characterized by various nuclear magnetic resonance (NMR) techniques, Fourier transform-infrared (FT-IR) techniques and density functional theory (DFT) calculations, revealing the formation of sulfur-centred O-H⋯S⋯H-O bifurcated H-bonds. This work may provide a simple and efficient method for the quantitative analysis of weak H-bonds between small organic molecules.
Collapse
Affiliation(s)
- Xiantao Ma
- College of Chemistry and Chemical Engineering, Green Catalysis & Synthesis Key Laboratory of Xinyang City, Xinyang Normal University, Xinyang, Henan 464000, China.
| | - Yingying Zhu
- College of Chemistry and Chemical Engineering, Green Catalysis & Synthesis Key Laboratory of Xinyang City, Xinyang Normal University, Xinyang, Henan 464000, China.
| | - Jing Yu
- College of Chemistry and Chemical Engineering, Green Catalysis & Synthesis Key Laboratory of Xinyang City, Xinyang Normal University, Xinyang, Henan 464000, China.
| | - Geng Zhao
- College of Chemistry and Chemical Engineering, Green Catalysis & Synthesis Key Laboratory of Xinyang City, Xinyang Normal University, Xinyang, Henan 464000, China.
| | - Jiaxin Duanmu
- College of Chemistry and Chemical Engineering, Green Catalysis & Synthesis Key Laboratory of Xinyang City, Xinyang Normal University, Xinyang, Henan 464000, China.
| | - Yiyun Yuan
- College of Chemistry and Chemical Engineering, Green Catalysis & Synthesis Key Laboratory of Xinyang City, Xinyang Normal University, Xinyang, Henan 464000, China.
| | - Xue-Ping Chang
- College of Chemistry and Chemical Engineering, Green Catalysis & Synthesis Key Laboratory of Xinyang City, Xinyang Normal University, Xinyang, Henan 464000, China.
| | - Dongli Xu
- College of Chemistry and Chemical Engineering, Green Catalysis & Synthesis Key Laboratory of Xinyang City, Xinyang Normal University, Xinyang, Henan 464000, China.
| | - Qiuju Zhou
- College of Chemistry and Chemical Engineering, Green Catalysis & Synthesis Key Laboratory of Xinyang City, Xinyang Normal University, Xinyang, Henan 464000, China.
| |
Collapse
|
8
|
Sparrow ZM, Ernst BG, Joo PT, Lao KU, DiStasio RA. NENCI-2021. I. A large benchmark database of non-equilibrium non-covalent interactions emphasizing close intermolecular contacts. J Chem Phys 2021; 155:184303. [PMID: 34773949 DOI: 10.1063/5.0068862] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this work, we present NENCI-2021, a benchmark database of ∼8000 Non-Equilibirum Non-Covalent Interaction energies for a large and diverse selection of intermolecular complexes of biological and chemical relevance. To meet the growing demand for large and high-quality quantum mechanical data in the chemical sciences, NENCI-2021 starts with the 101 molecular dimers in the widely used S66 and S101 databases and extends the scope of these works by (i) including 40 cation-π and anion-π complexes, a fundamentally important class of non-covalent interactions that are found throughout nature and pose a substantial challenge to theory, and (ii) systematically sampling all 141 intermolecular potential energy surfaces (PESs) by simultaneously varying the intermolecular distance and intermolecular angle in each dimer. Designed with an emphasis on close contacts, the complexes in NENCI-2021 were generated by sampling seven intermolecular distances along each PES (ranging from 0.7× to 1.1× the equilibrium separation) and nine intermolecular angles per distance (five for each ion-π complex), yielding an extensive database of 7763 benchmark intermolecular interaction energies (Eint) obtained at the coupled-cluster with singles, doubles, and perturbative triples/complete basis set [CCSD(T)/CBS] level of theory. The Eint values in NENCI-2021 span a total of 225.3 kcal/mol, ranging from -38.5 to +186.8 kcal/mol, with a mean (median) Eint value of -1.06 kcal/mol (-2.39 kcal/mol). In addition, a wide range of intermolecular atom-pair distances are also present in NENCI-2021, where close intermolecular contacts involving atoms that are located within the so-called van der Waals envelope are prevalent-these interactions, in particular, pose an enormous challenge for molecular modeling and are observed in many important chemical and biological systems. A detailed symmetry-adapted perturbation theory (SAPT)-based energy decomposition analysis also confirms the diverse and comprehensive nature of the intermolecular binding motifs present in NENCI-2021, which now includes a significant number of primarily induction-bound dimers (e.g., cation-π complexes). NENCI-2021 thus spans all regions of the SAPT ternary diagram, thereby warranting a new four-category classification scheme that includes complexes primarily bound by electrostatics (3499), induction (700), dispersion (1372), or mixtures thereof (2192). A critical error analysis performed on a representative set of intermolecular complexes in NENCI-2021 demonstrates that the Eint values provided herein have an average error of ±0.1 kcal/mol, even for complexes with strongly repulsive Eint values, and maximum errors of ±0.2-0.3 kcal/mol (i.e., ∼±1.0 kJ/mol) for the most challenging cases. For these reasons, we expect that NENCI-2021 will play an important role in the testing, training, and development of next-generation classical and polarizable force fields, density functional theory approximations, wavefunction theory methods, and machine learning based intra- and inter-molecular potentials.
Collapse
Affiliation(s)
- Zachary M Sparrow
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Brian G Ernst
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Paul T Joo
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Ka Un Lao
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Robert A DiStasio
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| |
Collapse
|
9
|
Loipersberger M, Bertels LW, Lee J, Head-Gordon M. Exploring the Limits of Second- and Third-Order Møller-Plesset Perturbation Theories for Noncovalent Interactions: Revisiting MP2.5 and Assessing the Importance of Regularization and Reference Orbitals. J Chem Theory Comput 2021; 17:5582-5599. [PMID: 34382394 PMCID: PMC9948597 DOI: 10.1021/acs.jctc.1c00469] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
This work systematically assesses the influence of reference orbitals, regularization, and scaling on the performance of second- and third-order Møller-Plesset perturbation theory wave function methods for noncovalent interactions (NCIs). Testing on 19 data sets (A24, DS14, HB15, HSG, S22, X40, HW30, NC15, S66, AlkBind12, CO2Nitrogen16, HB49, Ionic43, TA13, XB18, Bauza30, CT20, XB51, and Orel26rad) covers a wide range of different NCIs including hydrogen bonding, dispersion, and halogen bonding. Inclusion of potential energy surfaces from different hydrogen bonds and dispersion-bound complexes gauges accuracy for nonequilibrium geometries. Fifteen methods are tested. In notation where nonstandard choices of orbitals are denoted as methods:orbitals, these are MP2, κ-MP2, SCS-MP2, OOMP2, κ-OOMP2, MP3, MP2.5, MP3:OOMP2, MP2.5:OOMP2, MP3:κ-OOMP2, MP2.5:κ-OOMP2, κ-MP3:κ-OOMP2, κ-MP2.5:κ-OOMP2, MP3:ωB97X-V, and MP2.5:ωB97X-V. Furthermore, we compare these methods to the ωB97M-V and B3LYP-D3 density functionals, as well as CCSD. We find that the κ-regularization (κ = 1.45 au was used throughout) improves the energetics in almost all data sets for both MP2 (in 17 out of 19 data sets) and OOMP2 (16 out of 19). The improvement is significant (e.g., the root-mean-square deviation (RMSD) for the S66 data set is 0.29 kcal/mol for κ-OOMP2 versus 0.67 kcal/mol for MP2) and for interactions between stable closed-shell molecules, not strongly dependent on the reference orbitals. Scaled MP3 (with a factor of 0.5) using κ-OOMP2 reference orbitals (MP2.5:κ-OOMP2) provides significantly more accurate results for NCIs across all data sets with noniterative O(N6) scaling (S66 data set RMSD: 0.10 kcal/mol). Across the entire data set of 356 points, the improvement over standard MP2.5 is approximately a factor of 2: RMSD for MP3:κ-OOMP2 is 0.25 vs 0.50 kcal/mol for MP2.5. The use of high-quality density functional reference orbitals (ωB97X-V) also significantly improves the results of MP2.5 for NCI over a Hartree-Fock orbital reference. All our assessments and conclusions are based on the use of the medium-sized aug-cc-pVTZ basis to yield results that are directly compared against complete basis set limit reference values.
Collapse
Affiliation(s)
| | - Luke W. Bertels
- Department of Chemistry, University of California, Berkeley, California 94720, USA,Present Address: Department of Chemistry, Virginia Tech, Blacksburg, VA 24061, USA
| | - Joonho Lee
- Department of Chemistry, University of California, Berkeley, California 94720, USA,Present Address: Department of Chemistry, Columbia University, NY
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, USA,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
10
|
Abstract
The heavier chalcogen atoms S, Se, and Te can each participate in a range of different noncovalent interactions. They can serve as both proton donor and acceptor in H-bonds. Each atom can also act as electron acceptor in a chalcogen bond.
Collapse
Affiliation(s)
- Steve Scheiner
- Department of Chemistry and Biochemistry, Utah State University, Logan, Utah 84322-0300, USA
| |
Collapse
|
11
|
Otero-de-la-Roza A, DiLabio GA. Improved Basis-Set Incompleteness Potentials for Accurate Density-Functional Theory Calculations in Large Systems. J Chem Theory Comput 2020; 16:4176-4191. [PMID: 32470304 DOI: 10.1021/acs.jctc.0c00102] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The accurate calculation of chemical properties using density-functional theory (DFT) requires the use of a nearly complete basis set. In chemical systems involving hundreds to thousands of atoms, the cost of the calculations place practical limitations on the number of basis functions that can be used. Therefore, in most practical applications of DFT to large systems, there exists a basis-set incompleteness error (BSIE). In this article, we present the next iteration of the basis-set incompleteness potentials (BSIPs), one-electron potentials designed to correct for basis-set incompleteness error. The ultimate goal associated with the development of BSIPs is to allow the calculation of molecular properties using DFT with near-complete-basis-set results at a computational cost that is similar to a small basis set calculation. In this work, we develop BSIPs for 10 atoms in the first and second rows (H, B-F, Si-Cl) and 15 common basis sets of the Pople, Dunning, Karlsruhe, and Huzinaga types. Our new BSIPs are constructed to minimize BSIE in the calculation of reaction energies, barrier heights, noncovalent binding energies, and intermolecular distances. The BSIPs were obtained using a training set of 15 944 data points. The fitting approach employed a regularized linear least-squares method with variable selection (the LASSO method), which results in a much better fit to the training data than our previous BSIPs while, at the same time, reducing the computational cost of BSIP development. The proposed BSIPs are tested on various benchmark sets and demonstrate excellent performance in practice. Our new BSIPs are also transferable; i.e., they can be used to correct BSIE in calculations that employ density functionals other than the one used in the BSIP development (B3LYP). Finally, BSIPs can be used in any quantum chemistry program that have implemented effective-core potentials without changes to the software.
Collapse
Affiliation(s)
- A Otero-de-la-Roza
- Departamento de Quı́mica Fı́sica y Analítica and MALTA-Consolider Team, Facultad de Quı́mica, Universidad de Oviedo, 33006 Oviedo, Spain
| | - Gino A DiLabio
- Department of Chemistry, University of British Columbia, Okanagan, 3247 University Way, Kelowna, British Columbia V1V 1V7, Canada.,Faculty of Management, University of British Columbia, Okanagan, 1137 Alumni Avenue, Kelowna, British Columbia V1V 1V7, Canada
| |
Collapse
|
12
|
Perkins MA, Barlow KR, Dreux KM, Tschumper GS. Anchoring the hydrogen sulfide dimer potential energy surface to juxtapose (H2S)2 with (H2O)2. J Chem Phys 2020; 152:214306. [DOI: 10.1063/5.0008929] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Morgan A. Perkins
- Department of Chemistry and Biochemistry, University of Mississippi, University, Mississippi 38677-1848, USA
| | - Kayleigh R. Barlow
- Department of Chemistry and Biochemistry, University of Mississippi, University, Mississippi 38677-1848, USA
| | - Katelyn M. Dreux
- Department of Chemistry and Biochemistry, University of Mississippi, University, Mississippi 38677-1848, USA
| | - Gregory S. Tschumper
- Department of Chemistry and Biochemistry, University of Mississippi, University, Mississippi 38677-1848, USA
| |
Collapse
|
13
|
|
14
|
Gholivand K, Tizhoush SK, Kozakiewicz A, Eskandari K, Farshadfar K. Copper( i) complexes of functionalized sulfur-containing ligands: structural and theoretical insights into chalcogen bonding. CrystEngComm 2019. [DOI: 10.1039/c8ce02006j] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Four new copper(i) thiocyanate complexes were studied using geometrical parameters and the lump–hole approach for justification of the strength and nature of chalcogen bonding.
Collapse
Affiliation(s)
- Khodayar Gholivand
- Department of Chemistry
- Faculty of Science
- Tarbiat Modares University
- Tehran
- Iran
| | - Samaneh K. Tizhoush
- Department of Chemistry
- Faculty of Science
- Tarbiat Modares University
- Tehran
- Iran
| | - Anna Kozakiewicz
- Faculty of Chemistry
- Nicolaus Copernicus University in Toruń
- 87-100 Toruń
- Poland
| | - Kiamars Eskandari
- Department of Chemistry
- Isfahan University of Technology
- Isfahan 84156-83111
- Iran
| | - Kaveh Farshadfar
- Department of Chemistry
- Faculty of Science
- Tarbiat Modares University
- Tehran
- Iran
| |
Collapse
|
15
|
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
|
16
|
Dreux KM, Tschumper GS. Examination of the structures, energetics, and vibrational frequencies of small sulfur‐containing prototypical dimers, (H
2
S)
2
and H
2
O/H
2
S. J Comput Chem 2018; 40:229-236. [DOI: 10.1002/jcc.25578] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2018] [Revised: 08/13/2018] [Accepted: 08/13/2018] [Indexed: 11/08/2022]
Affiliation(s)
- Katelyn M. Dreux
- Department of Chemistry and Biochemistry University of Mississippi University Mississippi, 38677‐1848
| | - Gregory S. Tschumper
- Department of Chemistry and Biochemistry University of Mississippi University Mississippi, 38677‐1848
| |
Collapse
|
17
|
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
|
18
|
Orabi EA, English AM. Modeling Protein S-Aromatic Motifs Reveals Their Structural and Redox Flexibility. J Phys Chem B 2018. [PMID: 29533644 DOI: 10.1021/acs.jpcb.8b00089] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
S-aromatic motifs are important noncovalent forces for protein stability and function but remain poorly understood. Hence, we performed quantum calculations at the MP2(full)/6-311++G(d,p) level on complexes between Cys (H2S, MeSH) and Met (Me2S) models with models of Phe (benzene, toluene), Trp (indole, 3-methylindole), Tyr (phenol, 4-methylphenol), and His (imidazole, 4-methylimidazole). The most stable gas-phase conformers exhibit binding energies of -2 to -6 kcal/mol, and the S atom lies perpendicular to the ring plane. This reveals preferential interaction with the ring π-system, except in the imidazoles where S binds edge-on to an N atom. Complexation tunes the gas-phase vertical ionization potentials of the ligands over as much as 1 eV, and strong σ- or π-type H-bonding supports charge transfer to the H-bond donor, rendering it more oxidizable. When the S atom acts as an H-bond acceptor (N/O-Har···S), calibration of the CHARMM36 force field (by optimizing pair-specific Lennard-Jones parameters) is required. Implementing the optimized parameters in molecular dynamics simulations in bulk water, we find stable S-aromatic complexes with binding free energies of -0.6 to -1.1 kcal/mol at ligand separations up to 8 Å. The aqueous S-aromatics exhibit flexible binding conformations, but edge-on conformers are less stable in water. Reflecting this, only 0.3 to 10% of the S-indole, S-phenol, and S-imidazole structures are stabilized by N/O-Har···S or S-H···Oar/Nar σ-type H-bonding. The wide range of energies and geometries found for S-aromatic interactions and their tunable redox properties expose the versatility and variability of the S-aromatic motif in proteins and allow us to predict a number of their reported properties.
Collapse
Affiliation(s)
- Esam A Orabi
- Centre for Research in Molecular Modeling (CERMM) and PROTEO , Department of Chemistry and Biochemistry , Concordia University , 7141 Sherbrooke Street West , Montréal , Québec H4B 1R6 , Canada
| | - Ann M English
- Centre for Research in Molecular Modeling (CERMM) and PROTEO , Department of Chemistry and Biochemistry , Concordia University , 7141 Sherbrooke Street West , Montréal , Québec H4B 1R6 , Canada
| |
Collapse
|
19
|
Řezáč J. Empirical Self-Consistent Correction for the Description of Hydrogen Bonds in DFTB3. J Chem Theory Comput 2017; 13:4804-4817. [DOI: 10.1021/acs.jctc.7b00629] [Citation(s) in RCA: 41] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Jan Řezáč
- Institute of Organic Chemistry
and Biochemistry, Czech Academy of Sciences, 166 10 Prague, Czech Republic
| |
Collapse
|
20
|
Otero-de-la-Roza A, DiLabio GA. Transferable Atom-Centered Potentials for the Correction of Basis Set Incompleteness Errors in Density-Functional Theory. J Chem Theory Comput 2017. [DOI: 10.1021/acs.jctc.7b00300] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- A. Otero-de-la-Roza
- Department
of Chemistry and ‡Faculty of Management, University of British Columbia, Okanagan, 3247
University Way, Kelowna, British Columbia, Canada V1V 1V7
| | - Gino A. DiLabio
- Department
of Chemistry and ‡Faculty of Management, University of British Columbia, Okanagan, 3247
University Way, Kelowna, British Columbia, Canada V1V 1V7
| |
Collapse
|
21
|
Yan XC, Robertson MJ, Tirado-Rives J, Jorgensen WL. Improved Description of Sulfur Charge Anisotropy in OPLS Force Fields: Model Development and Parameterization. J Phys Chem B 2017. [PMID: 28627890 DOI: 10.1021/acs.jpcb.7b04233] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
The atomic point-charge model used in most molecular mechanics force fields does not represent well the electronic anisotropy that is featured in many directional noncovalent interactions. Sulfur participates in several types of such interactions with its lone pairs and σ-holes. The current study develops a new model, via the addition of off-atom charged sites, for a variety of sulfur compounds in the OPLS-AA and OPLS/CM5 force fields to address the lack of charge anisotropy. Parameter optimization is carried out to reproduce liquid-state properties, torsional and noncovalent energetics from reliable quantum mechanical calculations, and electrostatic potentials. Significant improvements are obtained for computed free energies of hydration, reducing the mean unsigned errors from ca. 1.4 to 0.4-0.7 kcal/mol. Enhanced accuracy in directionality and energetics is also obtained for molecular complexes with sulfur-containing hydrogen and halogen bonds. Moreover, the new model reproduces the unusual conformational preferences of sulfur-containing compounds with 1,4-intramolecular chalcogen bonds. Transferability of the new force field parameters to cysteine and methionine is verified via molecular dynamic simulations of blocked dipeptides. The study demonstrates the effectiveness of using off-atom charge sites to address electronic anisotropy, and provides a parametrization methodology that can be applied to other systems.
Collapse
Affiliation(s)
- Xin Cindy Yan
- Department of Chemistry, Yale University , New Haven, Connecticut 06520-8107, United States
| | - Michael J Robertson
- Department of Chemistry, Yale University , New Haven, Connecticut 06520-8107, United States
| | - Julian Tirado-Rives
- Department of Chemistry, Yale University , New Haven, Connecticut 06520-8107, United States
| | - William L Jorgensen
- Department of Chemistry, Yale University , New Haven, Connecticut 06520-8107, United States
| |
Collapse
|
22
|
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
|
23
|
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
|
24
|
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
|
25
|
Forbes CR, Sinha SK, Ganguly HK, Bai S, Yap GPA, Patel S, Zondlo NJ. Insights into Thiol-Aromatic Interactions: A Stereoelectronic Basis for S-H/π Interactions. J Am Chem Soc 2017; 139:1842-1855. [PMID: 28080040 PMCID: PMC5890429 DOI: 10.1021/jacs.6b08415] [Citation(s) in RCA: 67] [Impact Index Per Article: 9.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Thiols can engage favorably with aromatic rings in S-H/π interactions, within abiological systems and within proteins. However, the underlying bases for S-H/π interactions are not well understood. The crystal structure of Boc-l-4-thiolphenylalanine tert-butyl ester revealed crystal organization centered on the interaction of the thiol S-H with the aromatic ring of an adjacent molecule, with a through-space Hthiol···Caromatic distance of 2.71 Å, below the 2.90 Å sum of the van der Waals radii of H and C. The nature of this interaction was further examined by DFT calculations, IR spectroscopy, solid-state NMR spectroscopy, and analysis of the Cambridge Structural Database. The S-H/π interaction was found to be driven significantly by favorable molecular orbital interactions, between an aromatic π donor orbital and the S-H σ* acceptor orbital (a π → σ* interaction). For comparison, a structural analysis of O-H/π interactions and of cation/π interactions of alkali metal cations with aromatic rings was conducted. Na+ and K+ exhibit a significant preference for the centroid of the aromatic ring and distances near the sum of the van der Waals and ionic radii, as expected for predominantly electrostatic interactions. Li+ deviates substantially from Na+ and K+. The S-H/π interaction differs from classical cation/π interactions by the preferential alignment of the S-H σ* toward the ring carbons and an aromatic π orbital rather than toward the aromatic centroid. These results describe a potentially broadly applicable approach to understanding the interactions of weakly polar bonds with π systems.
Collapse
Affiliation(s)
- Christina R. Forbes
- Department of Chemistry and Biochemistry, University of Delaware, Newark DE 19716, United States
| | | | | | - Shi Bai
- Department of Chemistry and Biochemistry, University of Delaware, Newark DE 19716, United States
| | - Glenn P. A. Yap
- Department of Chemistry and Biochemistry, University of Delaware, Newark DE 19716, United States
| | - Sandeep Patel
- Department of Chemistry and Biochemistry, University of Delaware, Newark DE 19716, United States
| | - Neal J. Zondlo
- Department of Chemistry and Biochemistry, University of Delaware, Newark DE 19716, United States
| |
Collapse
|
26
|
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
|
27
|
Orabi EA, English AM. Sulfur-Aromatic Interactions: Modeling Cysteine and Methionine Binding to Tyrosinate and Histidinium Ions to Assess Their Influence on Protein Electron Transfer. Isr J Chem 2016. [DOI: 10.1002/ijch.201600047] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
- Esam A. Orabi
- Department of Chemistry and Biochemistry; Concordia University; 7141 Sherbrooke Street West Montréal Québec H4B 1R6 Canada
- Center for Research in Molecular Modeling (CERMM)Quebec; Network for Research on Protein Function, Engineering, and Applications (PROTEO)
- On leave from Department of Chemistry, Faculty of Science; Assiut University; Assiut 71516 Egypt
| | - Ann M. English
- Department of Chemistry and Biochemistry; Concordia University; 7141 Sherbrooke Street West Montréal Québec H4B 1R6 Canada
- Center for Research in Molecular Modeling (CERMM)Quebec; Network for Research on Protein Function, Engineering, and Applications (PROTEO)
| |
Collapse
|
28
|
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
|
29
|
Christensen A, Kubař T, Cui Q, Elstner M. Semiempirical Quantum Mechanical Methods for Noncovalent Interactions for Chemical and Biochemical Applications. Chem Rev 2016; 116:5301-37. [PMID: 27074247 PMCID: PMC4867870 DOI: 10.1021/acs.chemrev.5b00584] [Citation(s) in RCA: 246] [Impact Index Per Article: 30.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2015] [Indexed: 12/28/2022]
Abstract
Semiempirical (SE) methods can be derived from either Hartree-Fock or density functional theory by applying systematic approximations, leading to efficient computational schemes that are several orders of magnitude faster than ab initio calculations. Such numerical efficiency, in combination with modern computational facilities and linear scaling algorithms, allows application of SE methods to very large molecular systems with extensive conformational sampling. To reliably model the structure, dynamics, and reactivity of biological and other soft matter systems, however, good accuracy for the description of noncovalent interactions is required. In this review, we analyze popular SE approaches in terms of their ability to model noncovalent interactions, especially in the context of describing biomolecules, water solution, and organic materials. We discuss the most significant errors and proposed correction schemes, and we review their performance using standard test sets of molecular systems for quantum chemical methods and several recent applications. The general goal is to highlight both the value and limitations of SE methods and stimulate further developments that allow them to effectively complement ab initio methods in the analysis of complex molecular systems.
Collapse
Affiliation(s)
- Anders
S. Christensen
- Department
of Chemistry and Theoretical Chemistry Institute, University of Wisconsin—Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Tomáš Kubař
- Institute of Physical
Chemistry & Center for Functional Nanostructures and Institute of Physical
Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
| | - Qiang Cui
- Department
of Chemistry and Theoretical Chemistry Institute, University of Wisconsin—Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Marcus Elstner
- Institute of Physical
Chemistry & Center for Functional Nanostructures and Institute of Physical
Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
| |
Collapse
|
30
|
Řezáč J, Hobza P. Benchmark Calculations of Interaction Energies in Noncovalent Complexes and Their Applications. Chem Rev 2016; 116:5038-71. [DOI: 10.1021/acs.chemrev.5b00526] [Citation(s) in RCA: 281] [Impact Index Per Article: 35.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jan Řezáč
- Institute
of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, 166 10 Prague, Czech Republic
| | - Pavel Hobza
- Institute
of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, 166 10 Prague, Czech Republic
- Regional
Centre of Advanced Technologies and Materials, Department of Physical
Chemistry, Palacký University, 771 46 Olomouc, Czech Republic
| |
Collapse
|
31
|
Řezáč J. Cuby: An integrative framework for computational chemistry. J Comput Chem 2016; 37:1230-7. [DOI: 10.1002/jcc.24312] [Citation(s) in RCA: 111] [Impact Index Per Article: 13.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2015] [Revised: 01/05/2016] [Accepted: 01/07/2016] [Indexed: 12/12/2022]
Affiliation(s)
- Jan Řezáč
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic; 166 10 Prague Czech Republic
| |
Collapse
|
32
|
Peterson C, Penchoff D, Wilson A. Prediction of Thermochemical Properties Across the Periodic Table. ANNUAL REPORTS IN COMPUTATIONAL CHEMISTRY 2016. [DOI: 10.1016/bs.arcc.2016.04.001] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
33
|
Christensen AS, Elstner M, Cui Q. Improving intermolecular interactions in DFTB3 using extended polarization from chemical-potential equalization. J Chem Phys 2015; 143:084123. [PMID: 26328834 PMCID: PMC4552706 DOI: 10.1063/1.4929335] [Citation(s) in RCA: 43] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2015] [Accepted: 08/10/2015] [Indexed: 11/14/2022] Open
Abstract
Semi-empirical quantum mechanical methods traditionally expand the electron density in a minimal, valence-only electron basis set. The minimal-basis approximation causes molecular polarization to be underestimated, and hence intermolecular interaction energies are also underestimated, especially for intermolecular interactions involving charged species. In this work, the third-order self-consistent charge density functional tight-binding method (DFTB3) is augmented with an auxiliary response density using the chemical-potential equalization (CPE) method and an empirical dispersion correction (D3). The parameters in the CPE and D3 models are fitted to high-level CCSD(T) reference interaction energies for a broad range of chemical species, as well as dipole moments calculated at the DFT level; the impact of including polarizabilities of molecules in the parameterization is also considered. Parameters for the elements H, C, N, O, and S are presented. The Root Mean Square Deviation (RMSD) interaction energy is improved from 6.07 kcal/mol to 1.49 kcal/mol for interactions with one charged species, whereas the RMSD is improved from 5.60 kcal/mol to 1.73 for a set of 9 salt bridges, compared to uncorrected DFTB3. For large water clusters and complexes that are dominated by dispersion interactions, the already satisfactory performance of the DFTB3-D3 model is retained; polarizabilities of neutral molecules are also notably improved. Overall, the CPE extension of DFTB3-D3 provides a more balanced description of different types of non-covalent interactions than Neglect of Diatomic Differential Overlap type of semi-empirical methods (e.g., PM6-D3H4) and PBE-D3 with modest basis sets.
Collapse
Affiliation(s)
- Anders S Christensen
- Department of Chemistry, University of Wisconsin-Madison, 1101 University Ave., Madison, Wisconsin 53706, USA
| | - Marcus Elstner
- Theoretische Chemische Biologie, Universität Karlsruhe, Kaiserstr. 12, 76131 Karlsruhe, Germany
| | - Qiang Cui
- Department of Chemistry, University of Wisconsin-Madison, 1101 University Ave., Madison, Wisconsin 53706, USA
| |
Collapse
|
34
|
Manzer S, Horn PR, Mardirossian N, Head-Gordon M. Fast, accurate evaluation of exact exchange: The occ-RI-K algorithm. J Chem Phys 2015; 143:024113. [PMID: 26178096 PMCID: PMC4506295 DOI: 10.1063/1.4923369] [Citation(s) in RCA: 37] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2015] [Accepted: 06/22/2015] [Indexed: 11/15/2022] Open
Abstract
Construction of the exact exchange matrix, K, is typically the rate-determining step in hybrid density functional theory, and therefore, new approaches with increased efficiency are highly desirable. We present a framework with potential for greatly improved efficiency by computing a compressed exchange matrix that yields the exact exchange energy, gradient, and direct inversion of the iterative subspace (DIIS) error vector. The compressed exchange matrix is constructed with one index in the compact molecular orbital basis and the other index in the full atomic orbital basis. To illustrate the advantages, we present a practical algorithm that uses this framework in conjunction with the resolution of the identity (RI) approximation. We demonstrate that convergence using this method, referred to hereafter as occupied orbital RI-K (occ-RI-K), in combination with the DIIS algorithm is well-behaved, that the accuracy of computed energetics is excellent (identical to conventional RI-K), and that significant speedups can be obtained over existing integral-direct and RI-K methods. For a 4400 basis function C68H22 hydrogen-terminated graphene fragment, our algorithm yields a 14× speedup over the conventional algorithm and a speedup of 3.3× over RI-K.
Collapse
Affiliation(s)
- Samuel Manzer
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, and Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Paul R Horn
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, and Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Narbe Mardirossian
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, and Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, and Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
35
|
Mardirossian N, Head-Gordon M. Mapping the genome of meta-generalized gradient approximation density functionals: The search for B97M-V. J Chem Phys 2015; 142:074111. [DOI: 10.1063/1.4907719] [Citation(s) in RCA: 232] [Impact Index Per Article: 25.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] 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
|
36
|
Hussain HB, Wilson KA, Wetmore SD. Serine and Cysteine π-Interactions in Nature: A Comparison of the Frequency, Structure, and Stability of Contacts Involving Oxygen and Sulfur. Aust J Chem 2015. [DOI: 10.1071/ch14598] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Abstract
Despite many DNA–protein π-interactions in high-resolution crystal structures, only four X–H···π or X···π interactions were found between serine (Ser) or cysteine (Cys) and DNA nucleobase π-systems in over 100 DNA–protein complexes (where X = O for Ser and X = S for Cys). Nevertheless, 126 non-covalent contacts occur between Ser or Cys and the aromatic amino acids in many binding arrangements within proteins. Furthermore, Ser and Cys protein–protein π-interactions occur with similar frequencies and strengths. Most importantly, due to the great stability that can be provided to biological macromolecules (up to –20 kJ mol–1 for neutral π-systems or –40 kJ mol–1 for cationic π-systems), Ser and Cys π-interactions should be considered when analyzing protein stability and function.
Collapse
|
37
|
Alauddin M, Biswal HS, Gloaguen E, Mons M. Intra-residue interactions in proteins: interplay between serine or cysteine side chains and backbone conformations, revealed by laser spectroscopy of isolated model peptides. Phys Chem Chem Phys 2014; 17:2169-78. [PMID: 25482851 DOI: 10.1039/c4cp04449e] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Intra-residue interactions play an important role in proteins by influencing local folding of the backbone. Taking advantage of the capability of gas phase experiments to provide relevant information on the intrinsic H-bonding pattern of isolated peptide chains, the intra-residue interactions of serine and cysteine residues, i.e., OH/SH···OC(i) C6 and NH(i···)O/S C5 interactions in Ser/Cys residues, are probed by laser spectroscopy of isolated peptides. The strength of these local side chain-main chain interactions, elegantly documented from their IR spectral features for well-defined conformations of the main chain, demonstrates that a subtle competition exists between the two types of intra-residue bond: the C6 H-bond is the major interaction with Ser, in contrast to Cys where C5 interaction takes over. The restricted number of conformers observed in the gas phase experiment with Ser compared to Cys (where both extended and folded forms are observed) also suggests a significant mediation role of these intra-residue interactions on the competition between the several main chain folding patterns.
Collapse
Affiliation(s)
- Mohammad Alauddin
- CEA, IRAMIS, Laboratoire Interactions, Dynamique et Lasers, CEA Saclay, Bât 522, 91191 Gif-sur-Yvette, France.
| | | | | | | |
Collapse
|
38
|
Lian P, Guo HB, Riccardi D, Dong A, Parks JM, Xu Q, Pai EF, Miller SM, Wei DQ, Smith JC, Guo H. X-ray structure of a Hg2+ complex of mercuric reductase (MerA) and quantum mechanical/molecular mechanical study of Hg2+ transfer between the C-terminal and buried catalytic site cysteine pairs. Biochemistry 2014; 53:7211-22. [PMID: 25343681 PMCID: PMC4245977 DOI: 10.1021/bi500608u] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
![]()
Mercuric reductase, MerA, is a key
enzyme in bacterial mercury
resistance. This homodimeric enzyme captures and reduces toxic Hg2+ to Hg0, which is relatively unreactive and can
exit the cell passively. Prior to reduction, the Hg2+ is
transferred from a pair of cysteines (C558′ and C559′
using Tn501 numbering) at the C-terminus of one monomer
to another pair of cysteines (C136 and C141) in the catalytic site
of the other monomer. Here, we present the X-ray structure of the
C-terminal Hg2+ complex of the C136A/C141A double mutant
of the Tn501 MerA catalytic core and explore the
molecular mechanism of this Hg transfer with quantum mechanical/molecular
mechanical (QM/MM) calculations. The transfer is found to be nearly
thermoneutral and to pass through a stable tricoordinated intermediate
that is marginally less stable than the two end states. For the overall
process, Hg2+ is always paired with at least two thiolates
and thus is present at both the C-terminal and catalytic binding sites
as a neutral complex. Prior to Hg2+ transfer, C141 is negatively
charged. As Hg2+ is transferred into the catalytic site,
a proton is transferred from C136 to C559′ while C558′
becomes negatively charged, resulting in the net transfer of a negative
charge over a distance of ∼7.5 Å. Thus, the transport
of this soft divalent cation is made energetically feasible by pairing
a competition between multiple Cys thiols and/or thiolates for Hg2+ with a competition between the Hg2+ and protons
for the thiolates.
Collapse
Affiliation(s)
- Peng Lian
- The State Key Laboratory of Microbial Metabolism and College of Life Sciences and Biotechnology, Shanghai Jiao Tong University , Shanghai 200240, China
| | | | | | | | | | | | | | | | | | | | | |
Collapse
|
39
|
Shao Y, Gan Z, Epifanovsky E, Gilbert AT, Wormit M, Kussmann J, Lange AW, Behn A, Deng J, Feng X, Ghosh D, Goldey M, Horn PR, Jacobson LD, Kaliman I, Khaliullin RZ, Kuś T, Landau A, Liu J, Proynov EI, Rhee YM, Richard RM, Rohrdanz MA, Steele RP, Sundstrom EJ, Woodcock HL, Zimmerman PM, Zuev D, Albrecht B, Alguire E, Austin B, Beran GJO, Bernard YA, Berquist E, Brandhorst K, Bravaya KB, Brown ST, Casanova D, Chang CM, Chen Y, Chien SH, Closser KD, Crittenden DL, Diedenhofen M, DiStasio RA, Do H, Dutoi AD, Edgar RG, Fatehi S, Fusti-Molnar L, Ghysels A, Golubeva-Zadorozhnaya A, Gomes J, Hanson-Heine MW, Harbach PH, Hauser AW, Hohenstein EG, Holden ZC, Jagau TC, Ji H, Kaduk B, Khistyaev K, Kim J, Kim J, King RA, Klunzinger P, Kosenkov D, Kowalczyk T, Krauter CM, Lao KU, Laurent AD, Lawler KV, Levchenko SV, Lin CY, Liu F, Livshits E, Lochan RC, Luenser A, Manohar P, Manzer SF, Mao SP, Mardirossian N, Marenich AV, Maurer SA, Mayhall NJ, Neuscamman E, Oana CM, Olivares-Amaya R, O’Neill DP, Parkhill JA, Perrine TM, Peverati R, Prociuk A, Rehn DR, Rosta E, Russ NJ, Sharada SM, Sharma S, Small DW, Sodt A, Stein T, Stück D, Su YC, Thom AJ, Tsuchimochi T, Vanovschi V, Vogt L, Vydrov O, Wang T, Watson MA, Wenzel J, White A, Williams CF, Yang J, Yeganeh S, Yost SR, You ZQ, Zhang IY, Zhang X, Zhao Y, Brooks BR, Chan GK, Chipman DM, Cramer CJ, Goddard WA, Gordon MS, Hehre WJ, Klamt A, Schaefer HF, Schmidt MW, Sherrill CD, Truhlar DG, Warshel A, Xu X, Aspuru-Guzik A, Baer R, Bell AT, Besley NA, Chai JD, Dreuw A, Dunietz BD, Furlani TR, Gwaltney SR, Hsu CP, Jung Y, Kong J, Lambrecht DS, Liang W, Ochsenfeld C, Rassolov VA, Slipchenko LV, Subotnik JE, Van Voorhis T, Herbert JM, Krylov AI, Gill PM, Head-Gordon M. Advances in molecular quantum chemistry contained in the Q-Chem 4 program package. Mol Phys 2014. [DOI: 10.1080/00268976.2014.952696] [Citation(s) in RCA: 1769] [Impact Index Per Article: 176.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
|
40
|
Kaur D, Aulakh D, Khanna S, Singh H. Theoretical study on the nature of S···H and O ··· H hydrogen bonds. J Sulphur Chem 2014. [DOI: 10.1080/17415993.2013.874425] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
41
|
Mardirossian N, Head-Gordon M. ωB97X-V: A 10-parameter, range-separated hybrid, generalized gradient approximation density functional with nonlocal correlation, designed by a survival-of-the-fittest strategy. Phys Chem Chem Phys 2014; 16:9904-24. [DOI: 10.1039/c3cp54374a] [Citation(s) in RCA: 429] [Impact Index Per Article: 42.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
|
42
|
Jackson NE, Savoie BM, Kohlstedt KL, Olvera de la Cruz M, Schatz GC, Chen LX, Ratner MA. Controlling Conformations of Conjugated Polymers and Small Molecules: The Role of Nonbonding Interactions. J Am Chem Soc 2013; 135:10475-83. [DOI: 10.1021/ja403667s] [Citation(s) in RCA: 334] [Impact Index Per Article: 30.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Affiliation(s)
- Nicholas E. Jackson
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston,
Illinois 60208, United States
| | - Brett M. Savoie
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston,
Illinois 60208, United States
| | - Kevin L. Kohlstedt
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston,
Illinois 60208, United States
| | - Monica Olvera de la Cruz
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston,
Illinois 60208, United States
| | - George C. Schatz
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston,
Illinois 60208, United States
| | - Lin X. Chen
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston,
Illinois 60208, United States
- Chemical Sciences and Engineering
Division, Argonne National Laboratory,
9700 South Cass Avenue, Argonne, Illinois 60439, United States
| | - Mark A. Ratner
- Department of Chemistry, Northwestern University, 2145 Sheridan Road, Evanston,
Illinois 60208, United States
| |
Collapse
|