1
|
Metcalf DP, Glick ZL, Bortolato A, Jiang A, Cheney DL, Sherrill CD. Directional Δ G Neural Network (DrΔ G-Net): A Modular Neural Network Approach to Binding Free Energy Prediction. J Chem Inf Model 2024; 64:1907-1918. [PMID: 38470995 DOI: 10.1021/acs.jcim.3c02054] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/14/2024]
Abstract
The protein-ligand binding free energy is a central quantity in structure-based computational drug discovery efforts. Although popular alchemical methods provide sound statistical means of computing the binding free energy of a large breadth of systems, they are generally too costly to be applied at the same frequency as end point or ligand-based methods. By contrast, these data-driven approaches are typically fast enough to address thousands of systems but with reduced transferability to unseen systems. We introduce DrΔG-Net (or simply Dragnet), an equivariant graph neural network that can blend ligand-based and protein-ligand data-driven approaches. It is based on a 3D fingerprint representation of the ligand alone and in complex with the protein target. Dragnet is a global scoring function to predict the binding affinity of arbitrary protein-ligand complexes, but can be easily tuned via transfer learning to specific systems or end points, performing similarly to common 2D ligand-based approaches in these tasks. Dragnet is evaluated on a total of 28 validation proteins with a set of congeneric ligands derived from the Binding DB and one custom set extracted from the ChEMBL Database. In general, a handful of experimental binding affinities are sufficient to optimize the scoring function for a particular protein and ligand scaffold. When not available, predictions from physics-based methods such as absolute free energy perturbation can be used for the transfer learning tuning of Dragnet. Furthermore, we use our data to illustrate the present limitations of data-driven modeling of binding free energy predictions.
Collapse
Affiliation(s)
- Derek P Metcalf
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| | - Zachary L Glick
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| | - Andrea Bortolato
- Molecular Structure and Design, Bristol-Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, United States
| | - Andy Jiang
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| | - Daniel L Cheney
- Molecular Structure and Design, Bristol-Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, United States
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| |
Collapse
|
2
|
Spronk SA, Glick ZL, Metcalf DP, Sherrill CD, Cheney DL. A quantum chemical interaction energy dataset for accurately modeling protein-ligand interactions. Sci Data 2023; 10:619. [PMID: 37699937 PMCID: PMC10497680 DOI: 10.1038/s41597-023-02443-1] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2023] [Accepted: 08/03/2023] [Indexed: 09/14/2023] Open
Abstract
Fast and accurate calculation of intermolecular interaction energies is desirable for understanding many chemical and biological processes, including the binding of small molecules to proteins. The Splinter ["Symmetry-adapted perturbation theory (SAPT0) protein-ligand interaction"] dataset has been created to facilitate the development and improvement of methods for performing such calculations. Molecular fragments representing commonly found substructures in proteins and small-molecule ligands were paired into >9000 unique dimers, assembled into numerous configurations using an approach designed to adequately cover the breadth of the dimers' potential energy surfaces while enhancing sampling in favorable regions. ~1.5 million configurations of these dimers were randomly generated, and a structurally diverse subset of these were minimized to obtain an additional ~80 thousand local and global minima. For all >1.6 million configurations, SAPT0 calculations were performed with two basis sets to complete the dataset. It is expected that Splinter will be a useful benchmark dataset for training and testing various methods for the calculation of intermolecular interaction energies.
Collapse
Affiliation(s)
- Steven A Spronk
- Molecular Structure and Design, Bristol Myers Squibb Company, P. O. Box 5400, Princeton, NJ, 08543, USA.
| | - Zachary L Glick
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, GA, 30332-0400, USA
| | - Derek P Metcalf
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, GA, 30332-0400, USA
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, GA, 30332-0400, USA.
| | - Daniel L Cheney
- Molecular Structure and Design, Bristol Myers Squibb Company, P. O. Box 5400, Princeton, NJ, 08543, USA
| |
Collapse
|
3
|
Nelson PM, Glick ZL, Sherrill CD. Approximating large-basis coupled-cluster theory vibrational frequencies using focal-point approximations. J Chem Phys 2023; 159:094104. [PMID: 37655773 DOI: 10.1063/5.0168608] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2023] [Accepted: 08/09/2023] [Indexed: 09/02/2023] Open
Abstract
The focal-point approximation can be used to estimate a high-accuracy, slow quantum chemistry computation by combining several lower-accuracy, faster computations. We examine the performance of focal-point methods by combining second-order Møller-Plesset perturbation theory (MP2) with coupled-cluster singles, doubles, and perturbative triples [CCSD(T)] for the calculation of harmonic frequencies and that of fundamental frequencies using second-order vibrational perturbation theory (VPT2). In contrast to standard CCSD(T), the focal-point CCSD(T) method approaches the complete basis set (CBS) limit with only triple-ζ basis sets for the coupled-cluster portion of the computation. The predicted harmonic and fundamental frequencies were compared with the experimental values for a set of 20 molecules containing up to six atoms. The focal-point method combining CCSD(T)/aug-cc-pV(T + d)Z with CBS-extrapolated MP2 has mean absolute errors vs experiment of only 7.3 cm-1 for the fundamental frequencies, which are essentially the same as the mean absolute error for CCSD(T) extrapolated to the CBS limit using the aug-cc-pV(Q + d)Z and aug-cc-pV(5 + d)Z basis sets. However, for H2O, the focal-point procedure requires only 3% of the computation time as the extrapolated CCSD(T) result, and the cost savings will grow for larger molecules.
Collapse
Affiliation(s)
- Philip M Nelson
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Zachary L Glick
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| |
Collapse
|
4
|
Sargent CT, Kasera R, Glick ZL, Sherrill CD, Cheney DL. A quantitative assessment of deformation energy in intermolecular interactions: How important is it? J Chem Phys 2023; 158:244106. [PMID: 37352421 DOI: 10.1063/5.0155895] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/24/2023] [Accepted: 06/01/2023] [Indexed: 06/25/2023] Open
Abstract
Dimer interaction energies have been well studied in computational chemistry, but they can offer an incomplete understanding of molecular binding depending on the system. In the current study, we present a dataset of focal-point coupled-cluster interaction and deformation energies (summing to binding energies, De) of 28 organic molecular dimers. We use these highly accurate energies to evaluate ten density functional approximations for their accuracy. The best performing method (with a double-ζ basis set), B97M-D3BJ, is then used to calculate the binding energies of 104 organic dimers, and we analyze the influence of the nature and strength of interaction on deformation energies. Deformation energies can be as large as 50% of the dimer interaction energy, especially when hydrogen bonding is present. In most cases, two or more hydrogen bonds present in a dimer correspond to an interaction energy of -10 to -25 kcal mol-1, allowing a deformation energy above 1 kcal mol-1 (and up to 9.5 kcal mol-1). A lack of hydrogen bonding usually restricts the deformation energy to below 1 kcal mol-1 due to the weaker interaction energy.
Collapse
Affiliation(s)
- Caroline T Sargent
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Raina Kasera
- Molecular Structure and Design, Bristol Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - Zachary L Glick
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Daniel L Cheney
- Molecular Structure and Design, Bristol Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| |
Collapse
|
5
|
Borca CH, Glick ZL, Metcalf DP, Burns LA, Sherrill CD. Benchmark coupled-cluster lattice energy of crystalline benzene and assessment of multi-level approximations in the many-body expansion. J Chem Phys 2023; 158:234102. [PMID: 37318167 DOI: 10.1063/5.0159410] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2023] [Accepted: 05/25/2023] [Indexed: 06/16/2023] Open
Abstract
The many-body expansion (MBE) is promising for the efficient, parallel computation of lattice energies in organic crystals. Very high accuracy should be achievable by employing coupled-cluster singles, doubles, and perturbative triples at the complete basis set limit [CCSD(T)/CBS] for the dimers, trimers, and potentially tetramers resulting from the MBE, but such a brute-force approach seems impractical for crystals of all but the smallest molecules. Here, we investigate hybrid or multi-level approaches that employ CCSD(T)/CBS only for the closest dimers and trimers and utilize much faster methods like Møller-Plesset perturbation theory (MP2) for more distant dimers and trimers. For trimers, MP2 is supplemented with the Axilrod-Teller-Muto (ATM) model of three-body dispersion. MP2(+ATM) is shown to be a very effective replacement for CCSD(T)/CBS for all but the closest dimers and trimers. A limited investigation of tetramers using CCSD(T)/CBS suggests that the four-body contribution is entirely negligible. The large set of CCSD(T)/CBS dimer and trimer data should be valuable in benchmarking approximate methods for molecular crystals and allows us to see that a literature estimate of the core-valence contribution of the closest dimers to the lattice energy using just MP2 was overbinding by 0.5 kJ mol-1, and an estimate of the three-body contribution from the closest trimers using the T0 approximation in local CCSD(T) was underbinding by 0.7 kJ mol-1. Our CCSD(T)/CBS best estimate of the 0 K lattice energy is -54.01 kJ mol-1, compared to an estimated experimental value of -55.3 ± 2.2 kJ mol-1.
Collapse
Affiliation(s)
- Carlos H Borca
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
- Department of Chemical and Biological Engineering, School of Engineering and Applied Science, Princeton University, 41 Olden Street, Princeton, New Jersey 08544, USA
| | - Zachary L Glick
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - Derek P Metcalf
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - Lori A Burns
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
- School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| |
Collapse
|
6
|
Kumawat RL, Sherrill CD. High-Order Quantum-Mechanical Analysis of Hydrogen Bonding in Hachimoji and Natural DNA Base Pairs. J Chem Inf Model 2023; 63:3150-3157. [PMID: 37125692 DOI: 10.1021/acs.jcim.3c00428] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/02/2023]
Abstract
High-order quantum chemistry is applied to hydrogen-bonded natural DNA nucleobase pairs [adenine:thymine (A:T) and guanine:cytosine (G:C)] and non-natural Hachimoji nucleobase pairs [isoguanine:1-methylcytosine (B:S) and 2-aminoimidazo[1,2a][1,3,5]triazin-4(1H)-one:6-amino-5-nitropyridin-2-one (P:Z)] to see how the intermolecular interaction energies and their energetic components (electrostatics, exchange-repulsion, induction/polarization, and London dispersion interactions) vary among the base pairs. We examined the Hoogsteen (HG) geometries in addition to the traditional Watson-Crick (WC) geometries. Coupled-cluster theory through perturbative triples [CCSD(T)] extrapolated to the complete basis set (CBS) limit and high-order symmetry-adapted perturbation theory (SAPT) at the SAPT2+(3)(CCD)δMP2/aug-cc-pVTZ level are used to estimate highly accurate noncovalent interaction energies. Electrostatic interactions are the most attractive component of the interaction energies, but the sum of induction/polarization and London dispersion is nearly as large, for all base pairs and geometries considered. Interestingly, the non-natural Hachimoji base pairs interact more strongly than the corresponding natural base pairs, by -21.8 (B:S) and -0.3 (P:Z) kcal mol-1 in the WC geometries, according to CCSD(T)/CBS. This is consistent with the H-bond distances being generally shorter in the non-natural base pairs. The natural base pairs are energetically more stabilized in their Hoogsteen geometries than in their WC geometries. The Hoogsteen geometry makes the A:T base pair slightly more stable, by -0.8 kcal mol-1, and it greatly stabilizes the G:C+ base pair, by -15.3 kcal mol-1. The G:C+ stabilization is mainly due to the fact that C has typically added a proton when found in Hoogsteen geometries. By contrast, Hoogsteen geometries are substantially less favorable than WC geometries for non-natural Hachimoji base pairs, by 17.3 (B:S) and 13.8 (P:Z) kcal mol-1.
Collapse
Affiliation(s)
- Rameshwar L Kumawat
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| |
Collapse
|
7
|
Xie Y, Glick ZL, Sherrill CD. Assessment of three-body dispersion models against coupled-cluster benchmarks for crystalline benzene, carbon dioxide, and triazine. J Chem Phys 2023; 158:094110. [PMID: 36889937 DOI: 10.1063/5.0143712] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/09/2023] Open
Abstract
To study the contribution of three-body dispersion to crystal lattice energies, we compute the three-body contributions to the lattice energies for crystalline benzene, carbon dioxide, and triazine using various computational methods. We show that these contributions converge quickly as the intermolecular distances between the monomers grow. In particular, the smallest value among the three pairwise intermonomer closest-contact distances, Rmin, shows a strong correlation with the three-body contribution to the lattice energy, and, here, the largest of the closest-contact distances, Rmax, serves as a cutoff criterion to limit the number of trimers to be considered. We considered all trimers up to Rmax=15Å. The trimers with Rmin<4Å contribute 90.4%, 90.6%, and 93.9% of the total three-body contributions for crystalline benzene, carbon dioxide, and triazine, respectively, for the coupled-cluster singles, doubles, and perturbative triples [CCSD(T)] method. For trimers with Rmin>4Å, the second-order Møller-Plesset perturbation theory (MP2) supplemented with the Axilrod-Teller-Muto (ATM) three-body dispersion correction reproduces the CCSD(T) values for the cumulative three-body contributions with errors of less than 0.1 kJ mol-1. Moreover, three-body contributions are converged within 0.15 kJ mol-1 by Rmax=10Å. From these results, it appears that in molecular crystals where dispersion dominates the three-body contribution to the lattice energy, the trimers with Rmin>4Å can be computed with the MP2+ATM method to reduce the computational cost, and those with Rmax>10Å appear to be basically negligible.
Collapse
Affiliation(s)
- Yi Xie
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Zachary L Glick
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| |
Collapse
|
8
|
Ceriotti M, Jensen L, Manolopoulos DE, Martinez T, Reichman DR, Sciortino F, Sherrill CD, Shi Q, Vega C, Wang LS, Weiss EA, Zhu X, Stein J, Lian T. 2021 JCP Emerging Investigator Special Collection. J Chem Phys 2023; 158:060401. [PMID: 36792492 DOI: 10.1063/5.0143234] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/15/2023] Open
|
9
|
Sargent CT, Metcalf DP, Glick ZL, Borca CH, Sherrill CD. Benchmarking two-body contributions to crystal lattice energies and a range-dependent assessment of approximate methods. J Chem Phys 2023; 158:054112. [PMID: 36754814 DOI: 10.1063/5.0141872] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023] Open
Abstract
Using the many-body expansion to predict crystal lattice energies (CLEs), a pleasantly parallel process, allows for flexibility in the choice of theoretical methods. Benchmark-level two-body contributions to CLEs of 23 molecular crystals have been computed using interaction energies of dimers with minimum inter-monomer separations (i.e., closest contact distances) up to 30 Å. In a search for ways to reduce the computational expense of calculating accurate CLEs, we have computed these two-body contributions with 15 different quantum chemical levels of theory and compared these energies to those computed with coupled-cluster in the complete basis set (CBS) limit. Interaction energies of the more distant dimers are easier to compute accurately and several of the methods tested are suitable as replacements for coupled-cluster through perturbative triples for all but the closest dimers. For our dataset, sub-kJ mol-1 accuracy can be obtained when calculating two-body interaction energies of dimers with separations shorter than 4 Å with coupled-cluster with single, double, and perturbative triple excitations/CBS and dimers with separations longer than 4 Å with MP2.5/aug-cc-pVDZ, among other schemes, reducing the number of dimers to be computed with coupled-cluster by as much as 98%.
Collapse
Affiliation(s)
- Caroline T Sargent
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Derek P Metcalf
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Zachary L Glick
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Carlos H Borca
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| |
Collapse
|
10
|
Wang Z, Zou X, Xie Y, Zhang H, Hu L, Chan CCS, Zhang R, Guo J, Kwok RTK, Lam JWY, Williams ID, Zeng Z, Wong KS, Sherrill CD, Ye R, Tang BZ. A nonconjugated radical polymer with stable red luminescence in the solid state. Mater Horiz 2022; 9:2564-2571. [PMID: 35880529 DOI: 10.1039/d2mh00808d] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Organic radicals are unstable and stable radicals usually display non-luminescent properties. Luminescent radicals possess the all-in-one properties of optoelectronics, electronics, and magnetics. To date, the reported structures of luminescent radicals are limited to triphenylmethyl radical derivatives and their analogues, which are stabilized with extended π-conjugation. Here, we demonstrate the first example of a nonconjugated luminescent radical. In spite of the lack of delocalized π-stabilization, the radical polymer readily emits red luminescence in the solid state. A traditional luminescent quencher, 2,2,6,6-tetramethylpiperidin-1-yl turned into a red chromophore when grafted onto a polymer backbone. Experimental data confirm that the emission is associated with the nitroxide radicals and is also affected by the packing of the polymer. This work discloses a novel class of luminescent radicals and a distinctive pathway for luminescence from open-shell materials.
Collapse
Affiliation(s)
- Zhaoyu Wang
- Department of Chemistry, Hong Kong Branch of Chinese National Engineering Research Center for Tissue Restoration and Reconstruction, Institute for Advanced Study, Guangdong-Hong Kong-Macao Joint Laboratory of Optoelectronic and Magnetic Functional Materials, Department of Chemical and Biological Engineering, and Department of Physics, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China
| | - Xinhui Zou
- Department of Chemistry, Hong Kong Branch of Chinese National Engineering Research Center for Tissue Restoration and Reconstruction, Institute for Advanced Study, Guangdong-Hong Kong-Macao Joint Laboratory of Optoelectronic and Magnetic Functional Materials, Department of Chemical and Biological Engineering, and Department of Physics, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China
| | - Yi Xie
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia, 30332-0400, USA
| | - Haoke Zhang
- Department of Chemistry, Hong Kong Branch of Chinese National Engineering Research Center for Tissue Restoration and Reconstruction, Institute for Advanced Study, Guangdong-Hong Kong-Macao Joint Laboratory of Optoelectronic and Magnetic Functional Materials, Department of Chemical and Biological Engineering, and Department of Physics, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China
| | - Lianrui Hu
- Department of Chemistry, Hong Kong Branch of Chinese National Engineering Research Center for Tissue Restoration and Reconstruction, Institute for Advanced Study, Guangdong-Hong Kong-Macao Joint Laboratory of Optoelectronic and Magnetic Functional Materials, Department of Chemical and Biological Engineering, and Department of Physics, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China
| | - Christopher C S Chan
- Department of Chemistry, Hong Kong Branch of Chinese National Engineering Research Center for Tissue Restoration and Reconstruction, Institute for Advanced Study, Guangdong-Hong Kong-Macao Joint Laboratory of Optoelectronic and Magnetic Functional Materials, Department of Chemical and Biological Engineering, and Department of Physics, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China
| | - Ruoyao Zhang
- Department of Chemistry, Hong Kong Branch of Chinese National Engineering Research Center for Tissue Restoration and Reconstruction, Institute for Advanced Study, Guangdong-Hong Kong-Macao Joint Laboratory of Optoelectronic and Magnetic Functional Materials, Department of Chemical and Biological Engineering, and Department of Physics, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China
| | - Jing Guo
- State Key Laboratory of Chemo/Biosensing and Chemometrics, College of Chemistry and Chemical Engineering, Hunan University, Changsha, 410082, P. R. China
| | - Ryan T K Kwok
- Department of Chemistry, Hong Kong Branch of Chinese National Engineering Research Center for Tissue Restoration and Reconstruction, Institute for Advanced Study, Guangdong-Hong Kong-Macao Joint Laboratory of Optoelectronic and Magnetic Functional Materials, Department of Chemical and Biological Engineering, and Department of Physics, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China
| | - Jacky W Y Lam
- Department of Chemistry, Hong Kong Branch of Chinese National Engineering Research Center for Tissue Restoration and Reconstruction, Institute for Advanced Study, Guangdong-Hong Kong-Macao Joint Laboratory of Optoelectronic and Magnetic Functional Materials, Department of Chemical and Biological Engineering, and Department of Physics, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China
| | - Ian D Williams
- Department of Chemistry, Hong Kong Branch of Chinese National Engineering Research Center for Tissue Restoration and Reconstruction, Institute for Advanced Study, Guangdong-Hong Kong-Macao Joint Laboratory of Optoelectronic and Magnetic Functional Materials, Department of Chemical and Biological Engineering, and Department of Physics, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China
| | - Zebing Zeng
- State Key Laboratory of Chemo/Biosensing and Chemometrics, College of Chemistry and Chemical Engineering, Hunan University, Changsha, 410082, P. R. China
| | - Kam Sing Wong
- Department of Chemistry, Hong Kong Branch of Chinese National Engineering Research Center for Tissue Restoration and Reconstruction, Institute for Advanced Study, Guangdong-Hong Kong-Macao Joint Laboratory of Optoelectronic and Magnetic Functional Materials, Department of Chemical and Biological Engineering, and Department of Physics, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia, 30332-0400, USA
| | - Ruquan Ye
- Department of Chemistry, State Key Laboratory of Marine Pollution, City University of Hong Kong, Hong Kong, China.
| | - Ben Zhong Tang
- School of Science and Engineering, Shenzhen Institute of Aggregate Science and Technology, The Chinese University of Hong Kong, Shenzhen, Guangdong, 518172, China.
| |
Collapse
|
11
|
Metcalf DP, Smith AJ, Glick ZL, Sherrill CD. Range-dependence of two-body intermolecular interactions and their energy components in molecular crystals. J Chem Phys 2022; 157:084503. [DOI: 10.1063/5.0103644] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Routinely assessing the stability of molecular crystals with high accuracy remains an open challenge in the computational sciences. The many-body expansion decomposes computation of the crystal lattice energy into an embarrassingly parallel collection of computations over molecular dimers, trimers, and so forth, making quantum chemistry techniques tractable for many crystals of small organic molecules. By examining the range-dependence of different types of energetic contributions to the crystal lattice energy, we can glean qualitative understanding of solid-state intermolecular interactions as well as practical, exploitable reductions in the number of computations required for accurate energies. Here, we assess the range-dependent character of two-body interactions of 24 small organic molecular crystals using the physically interpretable components from symmetry-adapted perturbation theory (electrostatics, exchange repulsion, induction/polarization, and London dispersion). We also examine correlations between the convergence rates of electrostatics and London dispersion terms with molecular dipole moments and polarizabilities, to provide guidance for estimating convergence rates in other molecular crystals.
Collapse
Affiliation(s)
- Derek P Metcalf
- Chemistry & Biochemistry, Georgia Institute of Technology, United States of America
| | | | - Zachary Lee Glick
- Chemistry and Biochemistry, Georgia Institute of Technology College of Sciences, United States of America
| | - C. David Sherrill
- School of Chemistry and Biochemistry, Georgia Institute of Technology College of Sciences, United States of America
| |
Collapse
|
12
|
Xie Y, Smith DGA, Sherrill CD. Implementation of Symmetry-Adapted Perturbation Theory based on density functional theory and using hybrid exchange-correlation kernels for dispersion terms. J Chem Phys 2022; 157:024801. [DOI: 10.1063/5.0090688] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We report the implementation of a symmetry-adapted perturbation theory algorithm based on a density functional theory description of the monomers [SAPT(DFT)]. The implementation adopts a density-fitting treatment of hybrid exchange-correlation kernels to enable the description of monomers with hybrid functionals, as in the algorithm by Bukowski, Podeszwa, and Szalewicz [Chem. Phys. Lett. 414, 111(2005)]. We have improved the algorithm by increasing numerical stability with QR factorization, and optimized the computation of the exchange-correlation kernel with its 2-index density-fitted representation. The algorithm scales as O(N5) formally and is usable for systems with up to ∼3000 basis functions, as demonstrated forthe C60-buckycatcher complex with the aug-cc-pVDZ basis set. The hybrid-kernel-based SAPT(DFT) algorithm is shown to be as accurate as SAPT(DFT) implementations based on local effective exact exchange potentials obtained from the local Hartree-Fock (LHF) method, while avoiding the lower-scaling [O(N4)] but iterative and sometimes hard-to-converge LHF process. The hybrid-kernel algorithm outperforms Hartree-Fock-based SAPT (SAPT0) for the S66 test set, and its accuracy is comparable to the many-body perturbation theory based SAPT2+ approach, which scales as O(N7), although SAPT2+ exhibits a more narrow distribution of errors.
Collapse
Affiliation(s)
- Yi Xie
- Georgia Institute of Technology College of Sciences, United States of America
| | | | - C. David Sherrill
- School of Chemistry and Biochemistry, Georgia Institute of Technology College of Sciences, United States of America
| |
Collapse
|
13
|
Sirianni DA, Zhu X, Sitkoff DF, Cheney DL, Sherrill CD. The influence of a solvent environment on direct non-covalent interactions between two molecules: A symmetry-adapted perturbation theory study of polarization tuning of π-π interactions by water. J Chem Phys 2022; 156:194306. [PMID: 35597646 DOI: 10.1063/5.0087302] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
High-level quantum chemical computations have provided significant insight into the fundamental physical nature of non-covalent interactions. These studies have focused primarily on gas-phase computations of small van der Waals dimers; however, these interactions frequently take place in complex chemical environments, such as proteins, solutions, or solids. To better understand how the chemical environment affects non-covalent interactions, we have undertaken a quantum chemical study of π-π interactions in an aqueous solution, as exemplified by T-shaped benzene dimers surrounded by 28 or 50 explicit water molecules. We report interaction energies (IEs) using second-order Møller-Plesset perturbation theory, and we apply the intramolecular and functional-group partitioning extensions of symmetry-adapted perturbation theory (ISAPT and F-SAPT, respectively) to analyze how the solvent molecules tune the π-π interactions of the solute. For complexes containing neutral monomers, even 50 explicit waters (constituting a first and partial second solvation shell) change total SAPT IEs between the two solute molecules by only tenths of a kcal mol-1, while significant changes of up to 3 kcal mol-1 of the electrostatic component are seen for the cationic pyridinium-benzene dimer. This difference between charged and neutral solutes is attributed to large non-additive three-body interactions within solvated ion-containing complexes. Overall, except for charged solutes, our quantum computations indicate that nearby solvent molecules cause very little "tuning" of the direct solute-solute interactions. This indicates that differences in binding energies between the gas phase and solution phase are primarily indirect effects of the competition between solute-solute and solute-solvent interactions.
Collapse
Affiliation(s)
- Dominic A Sirianni
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Xiao Zhu
- Molecular Structure and Design, Bristol Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - Doree F Sitkoff
- Molecular Structure and Design, Bristol Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - Daniel L Cheney
- Molecular Structure and Design, Bristol Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| |
Collapse
|
14
|
Lachance-Brais C, Hennecker CD, Alenaizan A, Luo X, Toader V, Taing M, Sherrill CD, Mittermaier AK, Sleiman HF. Tuning DNA Supramolecular Polymers by the Addition of Small, Functionalized Nucleobase Mimics. J Am Chem Soc 2021; 143:19824-19833. [PMID: 34783562 DOI: 10.1021/jacs.1c08972] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Nucleobase mimicking small molecules able to reconfigure DNA are a recently discovered strategy that promises to extend the structural and functional diversity of nucleic acids. However, only simple, unfunctionalized molecules such as cyanuric acid and melamine have so far been used in this approach. In this work, we show that the addition of substituted cyanuric acid molecules can successfully program polyadenine strands to assemble into supramolecular fibers. Unlike conventional DNA nanostructure functionalization, which typically end-labels DNA strands, our approach incorporates functional groups into DNA with high density using small molecules and results in new DNA triple helices coated with alkylamine or alcohol units that grow into micrometer-long fibers. We find that small changes in the small molecule functional group can result in large structural and energetic variation in the overall assembly. A combination of circular dichroism, atomic force microscopy, molecular dynamics simulations, and a new thermodynamic method, transient equilibrium mapping, elucidated the molecular factors behind these large changes. In particular, we identify substantial DNA sugar and phosphate group deformations to accommodate a hydrogen bond between the phosphate and the small-molecule functional groups, as well as a critical chain length of the functional group which switches this interaction from intra- to interfiber. These parameters allow the controlled formation of hierarchical, hybrid DNA assemblies simply through the addition and variation of small, functionalized molecules.
Collapse
Affiliation(s)
| | - Christopher D Hennecker
- Department of Chemistry, McGill University, 801 Sherbrooke St. W., Montreal, QC H3A0B8, Canada
| | - Asem Alenaizan
- School of Chemistry and Biochemistry, Georgia Institute of Technology, 901 Atlantic Drive, Atlanta, Georgia 30332-0400, United States
| | - Xin Luo
- Department of Chemistry, McGill University, 801 Sherbrooke St. W., Montreal, QC H3A0B8, Canada
| | - Violeta Toader
- Department of Chemistry, McGill University, 801 Sherbrooke St. W., Montreal, QC H3A0B8, Canada
| | - Monica Taing
- Department of Chemistry, McGill University, 801 Sherbrooke St. W., Montreal, QC H3A0B8, Canada
| | - C David Sherrill
- School of Chemistry and Biochemistry, Georgia Institute of Technology, 901 Atlantic Drive, Atlanta, Georgia 30332-0400, United States
| | - Anthony K Mittermaier
- Department of Chemistry, McGill University, 801 Sherbrooke St. W., Montreal, QC H3A0B8, Canada
| | - Hanadi F Sleiman
- Department of Chemistry, McGill University, 801 Sherbrooke St. W., Montreal, QC H3A0B8, Canada
| |
Collapse
|
15
|
Smith DGA, Lolinco AT, Glick ZL, Lee J, Alenaizan A, Barnes TA, Borca CH, Di Remigio R, Dotson DL, Ehlert S, Heide AG, Herbst MF, Hermann J, Hicks CB, Horton JT, Hurtado AG, Kraus P, Kruse H, Lee SJR, Misiewicz JP, Naden LN, Ramezanghorbani F, Scheurer M, Schriber JB, Simmonett AC, Steinmetzer J, Wagner JR, Ward L, Welborn M, Altarawy D, Anwar J, Chodera JD, Dreuw A, Kulik HJ, Liu F, Martínez TJ, Matthews DA, Schaefer HF, Šponer J, Turney JM, Wang LP, De Silva N, King RA, Stanton JF, Gordon MS, Windus TL, Sherrill CD, Burns LA. Quantum Chemistry Common Driver and Databases (QCDB) and Quantum Chemistry Engine (QCEngine): Automation and interoperability among computational chemistry programs. J Chem Phys 2021; 155:204801. [PMID: 34852489 DOI: 10.1063/5.0059356] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
Community efforts in the computational molecular sciences (CMS) are evolving toward modular, open, and interoperable interfaces that work with existing community codes to provide more functionality and composability than could be achieved with a single program. The Quantum Chemistry Common Driver and Databases (QCDB) project provides such capability through an application programming interface (API) that facilitates interoperability across multiple quantum chemistry software packages. In tandem with the Molecular Sciences Software Institute and their Quantum Chemistry Archive ecosystem, the unique functionalities of several CMS programs are integrated, including CFOUR, GAMESS, NWChem, OpenMM, Psi4, Qcore, TeraChem, and Turbomole, to provide common computational functions, i.e., energy, gradient, and Hessian computations as well as molecular properties such as atomic charges and vibrational frequency analysis. Both standard users and power users benefit from adopting these APIs as they lower the language barrier of input styles and enable a standard layout of variables and data. These designs allow end-to-end interoperable programming of complex computations and provide best practices options by default.
Collapse
Affiliation(s)
- Daniel G A Smith
- Molecular Sciences Software Institute, Blacksburg, Virginia 24060, USA
| | | | - Zachary L Glick
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - Jiyoung Lee
- Department of Chemistry, Iowa State University, Ames, Iowa 50011, USA
| | - Asem Alenaizan
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - Taylor A Barnes
- Molecular Sciences Software Institute, Blacksburg, Virginia 24060, USA
| | - Carlos H Borca
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - Roberto Di Remigio
- Department of Chemistry, Centre for Theoretical and Computational Chemistry, UiT, The Arctic University of Norway, N-9037 Tromsø, Norway
| | - David L Dotson
- Open Force Field Initiative, University of Colorado Boulder, Boulder, Colorado 80309, USA
| | - Sebastian Ehlert
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Universität Bonn, Beringstraße 4, D-53115 Bonn, Germany
| | - Alexander G Heide
- Center for Computational Quantum Chemistry, University of Georgia, Athens, Georgia 30602, USA
| | - Michael F Herbst
- Applied and Computational Mathematics, RWTH Aachen University, Schinkelstr. 2, 52062 Aachen, Germany
| | - Jan Hermann
- FU Berlin, Department of Mathematics and Computer Science, 14195 Berlin, Germany
| | - Colton B Hicks
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Joshua T Horton
- Department of Chemistry, Lancaster University, Lancaster LA1 4YW, United Kingdom
| | - Adrian G Hurtado
- Institute for Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794-5250, USA
| | - Peter Kraus
- School of Molecular and Life Sciences, Curtin University, GPO Box U1987, Perth 6845, WA, Australia
| | - Holger Kruse
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | | | - Jonathon P Misiewicz
- Center for Computational Quantum Chemistry, University of Georgia, Athens, Georgia 30602, USA
| | - Levi N Naden
- Molecular Sciences Software Institute, Blacksburg, Virginia 24060, USA
| | | | - Maximilian Scheurer
- Interdisciplinary Center for Scientific Computing, Heidelberg University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Jeffrey B Schriber
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - Andrew C Simmonett
- Laboratory of Computational Biology, National Institutes of Health-National Heart, Lung and Blood Institute, Bethesda, Maryland 20892, USA
| | - Johannes Steinmetzer
- Institute of Physical Chemistry, Friedrich Schiller University Jena, Jena, Germany
| | - Jeffrey R Wagner
- Open Force Field Initiative, University of Colorado Boulder, Boulder, Colorado 80309, USA
| | - Logan Ward
- Data Science and Learning Division, Argonne National Laboratory, Lemont, Illinois 60439, USA
| | - Matthew Welborn
- Molecular Sciences Software Institute, Blacksburg, Virginia 24060, USA
| | - Doaa Altarawy
- Molecular Sciences Software Institute, Blacksburg, Virginia 24060, USA
| | - Jamshed Anwar
- Department of Chemistry, Lancaster University, Lancaster LA1 4YW, United Kingdom
| | - John D Chodera
- Computational and Systems Biology Program, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, New York 10065, USA
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Heidelberg University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Fang Liu
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Todd J Martínez
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Devin A Matthews
- The Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, Texas 78712, USA
| | - Henry F Schaefer
- Center for Computational Quantum Chemistry, University of Georgia, Athens, Georgia 30602, USA
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Justin M Turney
- Center for Computational Quantum Chemistry, University of Georgia, Athens, Georgia 30602, USA
| | - Lee-Ping Wang
- Department of Chemistry, University of California Davis, Davis, California 95616, USA
| | - Nuwan De Silva
- Department of Chemistry, Iowa State University, Ames, Iowa 50011, USA
| | - Rollin A King
- Department of Chemistry, Bethel University, St. Paul, Minnesota 55112, USA
| | - John F Stanton
- Quantum Theory Project, The University of Florida, 2328 New Physics Building, Gainesville, Florida 32611-8435, USA
| | - Mark S Gordon
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Theresa L Windus
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - Lori A Burns
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| |
Collapse
|
16
|
Affiliation(s)
- C David Sherrill
- School of Chemistry and Biochemistry and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - David E Manolopoulos
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, Oxford University, South Parks Road, Oxford OX1 3QZ, United Kingdom
| | - Todd J Martínez
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, USA
| | - Michele Ceriotti
- Laboratory of Computational Science and Modeling, Institute of Materials, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Angelos Michaelides
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| |
Collapse
|
17
|
Schriber JB, Sirianni DA, Smith DGA, Burns LA, Sitkoff D, Cheney DL, Sherrill CD. Optimized damping parameters for empirical dispersion corrections to symmetry-adapted perturbation theory. J Chem Phys 2021; 154:234107. [PMID: 34241276 DOI: 10.1063/5.0049745] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Symmetry-adapted perturbation theory (SAPT) has become an invaluable tool for studying the fundamental nature of non-covalent interactions by directly computing the electrostatics, exchange (steric) repulsion, induction (polarization), and London dispersion contributions to the interaction energy using quantum mechanics. Further application of SAPT is primarily limited by its computational expense, where even its most affordable variant (SAPT0) scales as the fifth power of system size [O(N5)] due to the dispersion terms. The algorithmic scaling of SAPT0 is reduced from O(N5)→O(N4) by replacing these terms with the empirical D3 dispersion correction of Grimme and co-workers, forming a method that may be termed SAPT0-D3. Here, we optimize the damping parameters for the -D3 terms in SAPT0-D3 using a much larger training set than has previously been considered, namely, 8299 interaction energies computed at the complete-basis-set limit of coupled cluster through perturbative triples [CCSD(T)/CBS]. Perhaps surprisingly, with only three fitted parameters, SAPT0-D3 improves on the accuracy of SAPT0, reducing mean absolute errors from 0.61 to 0.49 kcal mol-1 over the full set of complexes. Additionally, SAPT0-D3 exhibits a nearly 2.5× speedup over conventional SAPT0 for systems with ∼300 atoms and is applied here to systems with up to 459 atoms. Finally, we have also implemented a functional group partitioning of the approach (F-SAPT0-D3) and applied it to determine important contacts in the binding of salbutamol to G-protein coupled β1-adrenergic receptor in both active and inactive forms. SAPT0-D3 capabilities have been added to the open-source Psi4 software.
Collapse
Affiliation(s)
- Jeffrey B Schriber
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Dominic A Sirianni
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Daniel G A Smith
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Lori A Burns
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Doree Sitkoff
- Molecular Structure and Design, Bristol-Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - Daniel L Cheney
- Molecular Structure and Design, Bristol-Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| |
Collapse
|
18
|
Glick ZL, Koutsoukas A, Cheney DL, Sherrill CD. Cartesian message passing neural networks for directional properties: Fast and transferable atomic multipoles. J Chem Phys 2021; 154:224103. [PMID: 34241239 DOI: 10.1063/5.0050444] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
The message passing neural network (MPNN) framework is a promising tool for modeling atomic properties but is, until recently, incompatible with directional properties, such as Cartesian tensors. We propose a modified Cartesian MPNN (CMPNN) suitable for predicting atom-centered multipoles, an essential component of ab initio force fields. The efficacy of this model is demonstrated on a newly developed dataset consisting of 46 623 chemical structures and corresponding high-quality atomic multipoles, which was deposited into the publicly available Molecular Sciences Software Institute QCArchive server. We show that the CMPNN accurately predicts atom-centered charges, dipoles, and quadrupoles and that errors in the predicted atomic multipoles have a negligible effect on multipole-multipole electrostatic energies. The CMPNN is accurate enough to model conformational dependencies of a molecule's electronic structure. This opens up the possibility of recomputing atomic multipoles on the fly throughout a simulation in which they might exhibit strong conformational dependence.
Collapse
Affiliation(s)
- Zachary L Glick
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Alexios Koutsoukas
- Molecular Structure and Design, Bristol Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - Daniel L Cheney
- Molecular Structure and Design, Bristol Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| |
Collapse
|
19
|
Schriber JB, Nascimento DR, Koutsoukas A, Spronk SA, Cheney DL, Sherrill CD. CLIFF: A component-based, machine-learned, intermolecular force field. J Chem Phys 2021; 154:184110. [PMID: 34241025 DOI: 10.1063/5.0042989] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
Computation of intermolecular interactions is a challenge in drug discovery because accurate ab initio techniques are too computationally expensive to be routinely applied to drug-protein models. Classical force fields are more computationally feasible, and force fields designed to match symmetry adapted perturbation theory (SAPT) interaction energies can remain accurate in this context. Unfortunately, the application of such force fields is complicated by the laborious parameterization required for computations on new molecules. Here, we introduce the component-based machine-learned intermolecular force field (CLIFF), which combines accurate, physics-based equations for intermolecular interaction energies with machine-learning models to enable automatic parameterization. The CLIFF uses functional forms corresponding to electrostatic, exchange-repulsion, induction/polarization, and London dispersion components in SAPT. Molecule-independent parameters are fit with respect to SAPT2+(3)δMP2/aug-cc-pVTZ, and molecule-dependent atomic parameters (atomic widths, atomic multipoles, and Hirshfeld ratios) are obtained from machine learning models developed for C, N, O, H, S, F, Cl, and Br. The CLIFF achieves mean absolute errors (MAEs) no worse than 0.70 kcal mol-1 in both total and component energies across a diverse dimer test set. For the side chain-side chain interaction database derived from protein fragments, the CLIFF produces total interaction energies with an MAE of 0.27 kcal mol-1 with respect to reference data, outperforming similar and even more expensive methods. In applications to a set of model drug-protein interactions, the CLIFF is able to accurately rank-order ligand binding strengths and achieves less than 10% error with respect to SAPT reference values for most complexes.
Collapse
Affiliation(s)
- Jeffrey B Schriber
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30318, USA
| | - Daniel R Nascimento
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30318, USA
| | - Alexios Koutsoukas
- Molecular Structure and Design, Bristol Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - Steven A Spronk
- Molecular Structure and Design, Bristol Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - Daniel L Cheney
- Molecular Structure and Design, Bristol Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30318, USA
| |
Collapse
|
20
|
Alenaizan A, Borca CH, Karunakaran SC, Kendall AK, Stubbs G, Schuster GB, Sherrill CD, Hud NV. X-ray Fiber Diffraction and Computational Analyses of Stacked Hexads in Supramolecular Polymers: Insight into Self-Assembly in Water by Prospective Prebiotic Nucleobases. J Am Chem Soc 2021; 143:6079-6094. [DOI: 10.1021/jacs.0c12010] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Asem Alenaizan
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
- NSF-NASA Center for Chemical Evolution, Atlanta, Georgia 30332-0400, United States
- Center for Computational Molecular Science and Technology, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| | - Carlos H. Borca
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
- Center for Computational Molecular Science and Technology, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| | - Suneesh C. Karunakaran
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
- NSF-NASA Center for Chemical Evolution, Atlanta, Georgia 30332-0400, United States
| | - Amy K. Kendall
- Department of Biological Sciences and Center for Structural Biology, Vanderbilt University, Nashville, Tennessee 37235, United States
| | - Gerald Stubbs
- Department of Biological Sciences and Center for Structural Biology, Vanderbilt University, Nashville, Tennessee 37235, United States
| | - Gary B. Schuster
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| | - C. David Sherrill
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
- NSF-NASA Center for Chemical Evolution, Atlanta, Georgia 30332-0400, United States
- Center for Computational Molecular Science and Technology, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
- School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0765, United States
| | - Nicholas V. Hud
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
- NSF-NASA Center for Chemical Evolution, Atlanta, Georgia 30332-0400, United States
| |
Collapse
|
21
|
Alenaizan A, Fauché K, Krishnamurthy R, Sherrill CD. Noncovalent Helicene Structure between Nucleic Acids and Cyanuric Acid. Chemistry 2021; 27:4043-4052. [DOI: 10.1002/chem.202004390] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2020] [Revised: 11/08/2020] [Indexed: 12/16/2022]
Affiliation(s)
- Asem Alenaizan
- School of Chemistry and Biochemistry Georgia Institute of Technology Atlanta GA 30332-0400 USA
- Center for Computational Molecular Science and Technology Georgia Institute of Technology Atlanta GA 30332-0400 USA
- NSF-NASA Center for Chemical Evolution Atlanta GA 30332 USA
| | - Kévin Fauché
- NSF-NASA Center for Chemical Evolution Atlanta GA 30332 USA
- Department of Chemistry The Scripps Research Institute La Jolla CA 92037 USA
| | - Ramanarayanan Krishnamurthy
- NSF-NASA Center for Chemical Evolution Atlanta GA 30332 USA
- Department of Chemistry The Scripps Research Institute La Jolla CA 92037 USA
| | - C. David Sherrill
- School of Chemistry and Biochemistry Georgia Institute of Technology Atlanta GA 30332-0400 USA
- Center for Computational Molecular Science and Technology Georgia Institute of Technology Atlanta GA 30332-0400 USA
- NSF-NASA Center for Chemical Evolution Atlanta GA 30332 USA
- School of Computational Science and Engineering Georgia Institute of Technology Atlanta GA 30332-0765 USA
| |
Collapse
|
22
|
Alenaizan A, Barnett JL, Hud NV, Sherrill CD, Petrov AS. The proto-Nucleic Acid Builder: a software tool for constructing nucleic acid analogs. Nucleic Acids Res 2021; 49:79-89. [PMID: 33300028 PMCID: PMC7797056 DOI: 10.1093/nar/gkaa1159] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2020] [Revised: 11/09/2020] [Accepted: 11/13/2020] [Indexed: 11/13/2022] Open
Abstract
The helical structures of DNA and RNA were originally revealed by experimental data. Likewise, the development of programs for modeling these natural polymers was guided by known structures. These nucleic acid polymers represent only two members of a potentially vast class of polymers with similar structural features, but that differ from DNA and RNA in the backbone or nucleobases. Xeno nucleic acids (XNAs) incorporate alternative backbones that affect the conformational, chemical, and thermodynamic properties of XNAs. Given the vast chemical space of possible XNAs, computational modeling of alternative nucleic acids can accelerate the search for plausible nucleic acid analogs and guide their rational design. Additionally, a tool for the modeling of nucleic acids could help reveal what nucleic acid polymers may have existed before RNA in the early evolution of life. To aid the development of novel XNA polymers and the search for possible pre-RNA candidates, this article presents the proto-Nucleic Acid Builder (https://github.com/GT-NucleicAcids/pnab), an open-source program for modeling nucleic acid analogs with alternative backbones and nucleobases. The torsion-driven conformation search procedure implemented here predicts structures with good accuracy compared to experimental structures, and correctly demonstrates the correlation between the helical structure and the backbone conformation in DNA and RNA.
Collapse
Affiliation(s)
- Asem Alenaizan
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, GA 30332-0400, USA.,Center for Computational Molecular Science and Technology, Georgia Institute of Technology, Atlanta, GA 30332-0400, USA
| | - Joshua L Barnett
- School of Physics, Georgia Institute of Technology, Atlanta, GA 30332-0430, USA
| | - Nicholas V Hud
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, GA 30332-0400, USA
| | - C David Sherrill
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, GA 30332-0400, USA.,Center for Computational Molecular Science and Technology, Georgia Institute of Technology, Atlanta, GA 30332-0400, USA.,School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, GA 30332-0765, USA
| | - Anton S Petrov
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, GA 30332-0400, USA
| |
Collapse
|
23
|
Metcalf DP, Jiang A, Spronk SA, Cheney DL, Sherrill CD. Electron-Passing Neural Networks for Atomic Charge Prediction in Systems with Arbitrary Molecular Charge. J Chem Inf Model 2020; 61:115-122. [PMID: 33326247 DOI: 10.1021/acs.jcim.0c01071] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Abstract
Atomic charges are critical quantities in molecular mechanics and molecular dynamics, but obtaining these quantities requires heuristic choices based on atom typing or relatively expensive quantum mechanical computations to generate a density to be partitioned. Most machine learning efforts in this domain ignore total molecular charges, relying on overfitting and arbitrary rescaling in order to match the total system charge. Here, we introduce the electron-passing neural network (EPNN), a fast, accurate neural network atomic charge partitioning model that conserves total molecular charge by construction. EPNNs predict atomic charges very similar to those obtained by partitioning quantum mechanical densities but at such a small fraction of the cost that they can be easily computed for large biomolecules. Charges from this method may be used directly for molecular mechanics, as features for cheminformatics, or as input to any neural network potential.
Collapse
Affiliation(s)
- Derek P Metcalf
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| | - Andy Jiang
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| | - Steven A Spronk
- Molecular Structure and Design, Bristol Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, United States
| | - Daniel L Cheney
- Molecular Structure and Design, Bristol Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, United States
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| |
Collapse
|
24
|
Ediger MD, Jensen L, Manolopoulos DE, Martinez TJ, Michaelides A, Reichman DR, Sherrill CD, Shi Q, Straub JE, Vega C, Wang LS, Brigham EC, Lian T. JCP Emerging Investigator Special Collection 2019. J Chem Phys 2020; 153:110402. [PMID: 32962387 DOI: 10.1063/5.0021946] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Affiliation(s)
- Mark D Ediger
- Department of Chemistry, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| | - Lasse Jensen
- Department of Chemistry, Pennsylvania State University, University Park, Pennsylvania 16802, USA
| | - David E Manolopoulos
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford OX1 3QZ, United Kingdom
| | - Todd J Martinez
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Angelos Michaelides
- Thomas Young Centre and London Centre for Nanotechnology, 17-19 Gordon Street, London WC1H 0AH, United Kingdom and Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - David R Reichman
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Qiang Shi
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, Institute of Chemistry, Chinese Academy of Sciences, Zhongguancun, Beijing 100190, China; University of Chinese Academy of Sciences, Beijing 100049, China; and Physical Science Laboratory, Huairou National Comprehensive Science Center, Beijing 101407, China
| | - John E Straub
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, USA
| | - Carlos Vega
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Lai-Sheng Wang
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, USA
| | | | - Tianquan Lian
- Department of Chemistry, Emory University, Atlanta, Georgia 30322, USA
| |
Collapse
|
25
|
Affiliation(s)
- C. David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia
30332-0400, USA
| | - David E. Manolopoulos
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, Oxford University, South Parks Road, Oxford OX1
3QZ, United Kingdom
| | - Todd J. Martínez
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305,
USA
| | - Angelos Michaelides
- Thomas Young Centre, London Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| |
Collapse
|
26
|
Glick ZL, Metcalf DP, Koutsoukas A, Spronk SA, Cheney DL, Sherrill CD. AP-Net: An atomic-pairwise neural network for smooth and transferable interaction potentials. J Chem Phys 2020; 153:044112. [DOI: 10.1063/5.0011521] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Affiliation(s)
- Zachary L. Glick
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Derek P. Metcalf
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Alexios Koutsoukas
- Molecular Structure and Design, Bristol Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - Steven A. Spronk
- Molecular Structure and Design, Bristol Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - Daniel L. Cheney
- Molecular Structure and Design, Bristol Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - C. David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| |
Collapse
|
27
|
Smith DGA, Burns LA, Simmonett AC, Parrish RM, Schieber MC, Galvelis R, Kraus P, Kruse H, Di Remigio R, Alenaizan A, James AM, Lehtola S, Misiewicz JP, Scheurer M, Shaw RA, Schriber JB, Xie Y, Glick ZL, Sirianni DA, O’Brien JS, Waldrop JM, Kumar A, Hohenstein EG, Pritchard BP, Brooks BR, Schaefer HF, Sokolov AY, Patkowski K, DePrince AE, Bozkaya U, King RA, Evangelista FA, Turney JM, Crawford TD, Sherrill CD. Psi4 1.4: Open-source software for high-throughput quantum chemistry. J Chem Phys 2020; 152:184108. [PMID: 32414239 PMCID: PMC7228781 DOI: 10.1063/5.0006002] [Citation(s) in RCA: 301] [Impact Index Per Article: 75.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2020] [Accepted: 04/12/2020] [Indexed: 12/13/2022] Open
Abstract
PSI4 is a free and open-source ab initio electronic structure program providing implementations of Hartree-Fock, density functional theory, many-body perturbation theory, configuration interaction, density cumulant theory, symmetry-adapted perturbation theory, and coupled-cluster theory. Most of the methods are quite efficient, thanks to density fitting and multi-core parallelism. The program is a hybrid of C++ and Python, and calculations may be run with very simple text files or using the Python API, facilitating post-processing and complex workflows; method developers also have access to most of PSI4's core functionalities via Python. Job specification may be passed using The Molecular Sciences Software Institute (MolSSI) QCSCHEMA data format, facilitating interoperability. A rewrite of our top-level computation driver, and concomitant adoption of the MolSSI QCARCHIVE INFRASTRUCTURE project, makes the latest version of PSI4 well suited to distributed computation of large numbers of independent tasks. The project has fostered the development of independent software components that may be reused in other quantum chemistry programs.
Collapse
Affiliation(s)
| | - Lori A. Burns
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| | - Andrew C. Simmonett
- National Institutes of Health – National Heart,
Lung and Blood Institute, Laboratory of Computational Biology, Bethesda,
Maryland 20892, USA
| | - Robert M. Parrish
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| | - Matthew C. Schieber
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| | | | - Peter Kraus
- School of Molecular and Life Sciences, Curtin
University, Kent St., Bentley, Perth, Western Australia 6102,
Australia
| | - Holger Kruse
- Institute of Biophysics of the Czech Academy of
Sciences, Královopolská 135, 612 65 Brno, Czech
Republic
| | - Roberto Di Remigio
- Department of Chemistry, Centre for Theoretical
and Computational Chemistry, UiT, The Arctic University of Norway, N-9037
Tromsø, Norway
| | - Asem Alenaizan
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| | - Andrew M. James
- Department of Chemistry, Virginia
Tech, Blacksburg, Virginia 24061, USA
| | - Susi Lehtola
- Department of Chemistry, University of
Helsinki, P.O. Box 55 (A. I. Virtasen aukio 1), FI-00014 Helsinki,
Finland
| | - Jonathon P. Misiewicz
- Center for Computational Quantum Chemistry,
University of Georgia, Athens, Georgia 30602, USA
| | - Maximilian Scheurer
- Interdisciplinary Center for Scientific
Computing, Heidelberg University, D-69120 Heidelberg,
Germany
| | - Robert A. Shaw
- ARC Centre of Excellence in Exciton Science,
School of Science, RMIT University, Melbourne, VIC 3000,
Australia
| | - Jeffrey B. Schriber
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| | - Yi Xie
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| | - Zachary L. Glick
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| | - Dominic A. Sirianni
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| | - Joseph Senan O’Brien
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| | - Jonathan M. Waldrop
- Department of Chemistry and Biochemistry, Auburn
University, Auburn, Alabama 36849, USA
| | - Ashutosh Kumar
- Department of Chemistry, Virginia
Tech, Blacksburg, Virginia 24061, USA
| | - Edward G. Hohenstein
- SLAC National Accelerator Laboratory, Stanford
PULSE Institute, Menlo Park, California 94025,
USA
| | | | - Bernard R. Brooks
- National Institutes of Health – National Heart,
Lung and Blood Institute, Laboratory of Computational Biology, Bethesda,
Maryland 20892, USA
| | - Henry F. Schaefer
- Center for Computational Quantum Chemistry,
University of Georgia, Athens, Georgia 30602, USA
| | - Alexander Yu. Sokolov
- Department of Chemistry and Biochemistry, The
Ohio State University, Columbus, Ohio 43210, USA
| | - Konrad Patkowski
- Department of Chemistry and Biochemistry, Auburn
University, Auburn, Alabama 36849, USA
| | - A. Eugene DePrince
- Department of Chemistry and Biochemistry,
Florida State University, Tallahassee, Florida 32306-4390,
USA
| | - Uğur Bozkaya
- Department of Chemistry, Hacettepe
University, Ankara 06800, Turkey
| | - Rollin A. King
- Department of Chemistry, Bethel
University, St. Paul, Minnesota 55112, USA
| | | | - Justin M. Turney
- Center for Computational Quantum Chemistry,
University of Georgia, Athens, Georgia 30602, USA
| | | | - C. David Sherrill
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| |
Collapse
|
28
|
Warden CE, Smith DGA, Burns LA, Bozkaya U, Sherrill CD. Efficient and automated computation of accurate molecular geometries using focal-point approximations to large-basis coupled-cluster theory. J Chem Phys 2020; 152:124109. [DOI: 10.1063/5.0004863] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Constance E. Warden
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Daniel G. A. Smith
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Lori A. Burns
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Uğur Bozkaya
- Department of Chemistry, Hacettepe University, Ankara 06800, Turkey
| | - C. David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| |
Collapse
|
29
|
Metcalf DP, Koutsoukas A, Spronk SA, Claus BL, Loughney DA, Johnson SR, Cheney DL, Sherrill CD. Approaches for machine learning intermolecular interaction energies and application to energy components from symmetry adapted perturbation theory. J Chem Phys 2020; 152:074103. [DOI: 10.1063/1.5142636] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Affiliation(s)
- Derek P. Metcalf
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Alexios Koutsoukas
- Molecular Structure and Design, Bristol-Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - Steven A. Spronk
- Molecular Structure and Design, Bristol-Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - Brian L. Claus
- Molecular Structure and Design, Bristol-Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - Deborah A. Loughney
- Molecular Structure and Design, Bristol-Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - Stephen R. Johnson
- Molecular Structure and Design, Bristol-Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - Daniel L. Cheney
- Molecular Structure and Design, Bristol-Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - C. David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| |
Collapse
|
30
|
Abstract
This paper presents techniques for Fock matrix construction that are designed for high performance on shared and distributed memory parallel computers when using Gaussian basis sets. Four main techniques are considered. (1) To calculate electron repulsion integrals, we demonstrate batching together the calculation of multiple shell quartets of the same angular momentum class so that the calculation of large sets of primitive integrals can be efficiently vectorized. (2) For multithreaded summation of entries into the Fock matrix, we investigate using a combination of atomic operations and thread-local copies of the Fock matrix. (3) For distributed memory parallel computers, we present a globally accessible matrix class for accessing distributed Fock and density matrices. The new matrix class introduces a batched mode for remote memory access that can reduce the synchronization cost. (4) For density fitting, we exploit both symmetry (of the Coulomb and exchange matrices) and sparsity (of 3-index tensors) and give a performance comparison of density fitting and the conventional direct calculation approach. The techniques are implemented in an open-source software library called GTFock.
Collapse
Affiliation(s)
- Hua Huang
- School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-4017, USA
| | - C David Sherrill
- School of Chemistry and Biochemistry and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Edmond Chow
- School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-4017, USA
| |
Collapse
|
31
|
Borca CH, Bakr BW, Burns LA, Sherrill CD. CrystaLattE: Automated computation of lattice energies of organic crystals exploiting the many-body expansion to achieve dual-level parallelism. J Chem Phys 2019; 151:144103. [PMID: 31615262 DOI: 10.1063/1.5120520] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023] Open
Abstract
We present an algorithm to compute the lattice energies of molecular crystals based on the many-body cluster expansion. The required computations on dimers, trimers, etc., within the crystal are independent of each other, leading to a naturally parallel approach. The algorithm exploits the long-range three-dimensional periodic order of crystals to automatically detect and avoid redundant or unnecessary computations. For this purpose, Coulomb-matrix descriptors from machine learning applications are found to be efficient in determining whether two N-mers are identical. The algorithm is implemented as an open-source Python program, CrystaLattE, that uses some of the features of the Quantum Chemistry Common Driver and Databases library. CrystaLattE is initially interfaced with the quantum chemistry package Psi4. With CrystaLattE, we have applied the fast, dispersion-corrected Hartree-Fock method HF-3c to the lattice energy of crystalline benzene. Including all 73 symmetry-unique dimers and 7130 symmetry-unique trimers that can be formed from molecules within a 15 Å cutoff from a central reference monomer, HF-3c plus an Axilrod-Teller-Muto estimate of three-body dispersion exhibits an error of only -1.0 kJ mol-1 vs the estimated 0 K experimental lattice energy of -55.3 ± 2.2 kJ mol-1. The convergence of the HF-3c two- and three-body contributions to the lattice energy as a function of intermonomer distance is examined.
Collapse
Affiliation(s)
- Carlos H Borca
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - Brandon W Bakr
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - Lori A Burns
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - C David Sherrill
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| |
Collapse
|
32
|
Chen XK, Bakr BW, Auffray M, Tsuchiya Y, Sherrill CD, Adachi C, Bredas JL. Intramolecular Noncovalent Interactions Facilitate Thermally Activated Delayed Fluorescence (TADF). J Phys Chem Lett 2019; 10:3260-3268. [PMID: 31141375 DOI: 10.1021/acs.jpclett.9b01220] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/19/2023]
Abstract
In the conventional molecular design of thermally activated delayed fluorescence (TADF) organic emitters, simultaneously achieving a fast rate of reverse intersystem crossing (RISC) from the triplet to the singlet manifold and a fast rate of radiative decay is a challenging task. A number of recent experimental data, however, point to TADF emitters with intramolecular π-π interactions as a potential pathway to overcome the issue. Here, we report a comprehensive investigation of TADF emitters with intramolecular π···π or lone-pair···π noncovalent interactions. We uncover the impact of those intramolecular noncovalent interactions on the TADF properties. In particular, we find that folded geometries in TADF molecules can trigger lone-pair···π interactions, introduce a n → π* character of the relevant transitions, enhance the singlet-triplet spin-orbit coupling, and ultimately greatly facilitate the RISC process. This work provides a robust foundation for the molecular design of a novel class of highly efficient TADF emitters in which intramolecular noncovalent interactions play a critical function.
Collapse
Affiliation(s)
- Xian-Kai Chen
- School of Chemistry and Biochemistry and Center for Organic Photonics and Electronics , Georgia Institute of Technology , Atlanta , Georgia 30332-0400 , United States
| | - Brandon W Bakr
- School of Chemistry and Biochemistry and Center for Organic Photonics and Electronics , Georgia Institute of Technology , Atlanta , Georgia 30332-0400 , United States
| | - Morgan Auffray
- Center for Organic Photonics and Electronics Research (OPERA) , Kyushu University , 744 Motooka, Nishi , Fukuoka 819-0395 , Japan
| | - Youichi Tsuchiya
- Center for Organic Photonics and Electronics Research (OPERA) , Kyushu University , 744 Motooka, Nishi , Fukuoka 819-0395 , Japan
| | - C David Sherrill
- School of Chemistry and Biochemistry and Center for Organic Photonics and Electronics , Georgia Institute of Technology , Atlanta , Georgia 30332-0400 , United States
| | - Chihaya Adachi
- Center for Organic Photonics and Electronics Research (OPERA) , Kyushu University , 744 Motooka, Nishi , Fukuoka 819-0395 , Japan
- International Institute for Carbon Neutral Energy Research (WPI-I2CNER) , Kyushu University , 744 Motooka, Nishi-ku , Fukuoka 819-0395 , Japan
| | - Jean-Luc Bredas
- School of Chemistry and Biochemistry and Center for Organic Photonics and Electronics , Georgia Institute of Technology , Atlanta , Georgia 30332-0400 , United States
| |
Collapse
|
33
|
Brahmachari U, Gonthier JF, Sherrill CD, Barry BA. Water Bridges Conduct Sequential Proton Transfer in Photosynthetic Oxygen Evolution. J Phys Chem B 2019; 123:4487-4496. [DOI: 10.1021/acs.jpcb.9b01523] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
|
34
|
Hwang J, Li P, Smith MD, Warden CE, Sirianni DA, Vik EC, Maier JM, Yehl CJ, Sherrill CD, Shimizu KD. Tipping the Balance between S-π and O-π Interactions. J Am Chem Soc 2018; 140:13301-13307. [PMID: 30251855 DOI: 10.1021/jacs.8b07617] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
A comprehensive experimental survey consisting of 36 molecular balances was conducted to compare 18 pairs of S-π versus O-π interactions over a wide range of structural, geometric, and solvent parameters. A strong linear correlation was observed between the folding energies of the sulfur and oxygen balances across the entire library of balance pairs. The more stable interaction systematically switched from the O-π to S-π interaction. Computational studies of bimolecular PhSCH3-arene and PhOCH3-arene complexes were able to replicate the experimental trends in the molecular balances. The change in preference for the O-π to S-π interaction was due to the interplay of stabilizing (dispersion and solvophobic) and destabilizing (exchange-repulsion) terms arising from the differences in size and polarizability of the oxygen and sulfur atoms.
Collapse
Affiliation(s)
- Jungwun Hwang
- Department of Chemistry and Biochemistry , University of South Carolina , Columbia , South Carolina 29208 , United States
| | - Ping Li
- Department of Chemistry and Biochemistry , University of South Carolina , Columbia , South Carolina 29208 , United States
| | - Mark D Smith
- Department of Chemistry and Biochemistry , University of South Carolina , Columbia , South Carolina 29208 , United States
| | | | | | - Erik C Vik
- Department of Chemistry and Biochemistry , University of South Carolina , Columbia , South Carolina 29208 , United States
| | - Josef M Maier
- Department of Chemistry and Biochemistry , University of South Carolina , Columbia , South Carolina 29208 , United States
| | - Christopher J Yehl
- Department of Chemistry and Biochemistry , University of South Carolina , Columbia , South Carolina 29208 , United States
| | | | - Ken D Shimizu
- Department of Chemistry and Biochemistry , University of South Carolina , Columbia , South Carolina 29208 , United States
| |
Collapse
|
35
|
Bakr BW, Sherrill CD. Analysis of transition state stabilization by non-covalent interactions in organocatalysis: application of atomic and functional-group partitioned symmetry-adapted perturbation theory to the addition of organoboron reagents to fluoroketones. Phys Chem Chem Phys 2018; 20:18241-18251. [PMID: 29947381 DOI: 10.1039/c8cp02029a] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
This work seeks to apply symmetry-adapted perturbation theory (SAPT) to the recent study of Hoveyda and co-workers [K. A. Lee et al., Nat. Chem. 2016, 8, 768] where an allyl addition to a ketone became enantioselective when the ketone was fluorinated. Through the application of atomic SAPT (A-SAPT) and functional-group SAPT (F-SAPT), the non-covalent interactions between specific atoms and functional groups in the transition states associated with the fluoroketone reactions can be quantified. Our A-SAPT analysis confirms that a HF contact thought to enhance stereoselectivity shows a strong preference for one of the transition states leading to the experimentally observed product enantiomer. Other key atom-atom contacts invoked to rationalize relative transition state energies are also found to behave as expected based on chemical intuition and contact distances. On the other hand, hypothesized steric clashes between substrate phenyl or ortho-methyl phenyl groups and the catalyst are not supported by F-SAPT computations, and indeed, these are actually favorable π-π interactions.
Collapse
Affiliation(s)
- Brandon W Bakr
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, GA 30332-0400, USA.
| | | |
Collapse
|
36
|
Burns LA, Faver JC, Zheng Z, Marshall MS, Smith DGA, Vanommeslaeghe K, MacKerell AD, Merz KM, Sherrill CD. The BioFragment Database (BFDb): An open-data platform for computational chemistry analysis of noncovalent interactions. J Chem Phys 2018; 147:161727. [PMID: 29096505 DOI: 10.1063/1.5001028] [Citation(s) in RCA: 74] [Impact Index Per Article: 12.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Accurate potential energy models are necessary for reliable atomistic simulations of chemical phenomena. In the realm of biomolecular modeling, large systems like proteins comprise very many noncovalent interactions (NCIs) that can contribute to the protein's stability and structure. This work presents two high-quality chemical databases of common fragment interactions in biomolecular systems as extracted from high-resolution Protein DataBank crystal structures: 3380 sidechain-sidechain interactions and 100 backbone-backbone interactions that inaugurate the BioFragment Database (BFDb). Absolute interaction energies are generated with a computationally tractable explicitly correlated coupled cluster with perturbative triples [CCSD(T)-F12] "silver standard" (0.05 kcal/mol average error) for NCI that demands only a fraction of the cost of the conventional "gold standard," CCSD(T) at the complete basis set limit. By sampling extensively from biological environments, BFDb spans the natural diversity of protein NCI motifs and orientations. In addition to supplying a thorough assessment for lower scaling force-field (2), semi-empirical (3), density functional (244), and wavefunction (45) methods (comprising >1M interaction energies), BFDb provides interactive tools for running and manipulating the resulting large datasets and offers a valuable resource for potential energy model development and validation.
Collapse
Affiliation(s)
- Lori A Burns
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - John C Faver
- Quantum Theory Project, The University of Florida, 2328 New Physics Building, Gainesville, Florida 32611-8435, USA
| | - Zheng Zheng
- Quantum Theory Project, The University of Florida, 2328 New Physics Building, Gainesville, Florida 32611-8435, USA
| | - Michael S Marshall
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Daniel G A Smith
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Kenno Vanommeslaeghe
- Department of Analytical Chemistry and Pharmaceutical Technology (FABI), Center for Pharmaceutical Research (CePhaR), Vrije Universiteit Brussel (VUB), Laarbeeklaan 103, B-1090 Brussels, Belgium
| | - Alexander D MacKerell
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, Maryland 21201, USA
| | - Kenneth M Merz
- Quantum Theory Project, The University of Florida, 2328 New Physics Building, Gainesville, Florida 32611-8435, USA
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| |
Collapse
|
37
|
Bozkaya U, Sherrill CD. Analytic energy gradients for the coupled-cluster singles and doubles with perturbative triples method with the density-fitting approximation. J Chem Phys 2018; 147:044104. [PMID: 28764345 DOI: 10.1063/1.4994918] [Citation(s) in RCA: 30] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
An efficient implementation of analytic gradients for the coupled-cluster singles and doubles with perturbative triples [CCSD(T)] method with the density-fitting (DF) approximation, denoted as DF-CCSD(T), is reported. For the molecules considered, the DF approach substantially accelerates conventional CCSD(T) analytic gradients due to the reduced input/output time and the acceleration of the so-called "gradient terms": formation of particle density matrices (PDMs), computation of the generalized Fock-matrix (GFM), solution of the Z-vector equation, formation of the effective PDMs and GFM, back-transformation of the PDMs and GFM, from the molecular orbital to the atomic orbital (AO) basis, and computation of gradients in the AO basis. For the largest member of the molecular test set considered (C6H14), the computational times for analytic gradients (with the correlation-consistent polarized valence triple-ζ basis set in serial) are 106.2 [CCSD(T)] and 49.8 [DF-CCSD(T)] h, a speedup of more than 2-fold. In the evaluation of gradient terms, the DF approach completely avoids the use of four-index two-electron integrals. Similar to our previous studies on DF-second-order Møller-Plesset perturbation theory and DF-CCSD gradients, our formalism employs 2- and 3-index two-particle density matrices (TPDMs) instead of 4-index TPDMs. Errors introduced by the DF approximation are negligible for equilibrium geometries and harmonic vibrational frequencies.
Collapse
Affiliation(s)
- Uğur Bozkaya
- Department of Chemistry, Hacettepe University, Ankara 06800, Turkey
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| |
Collapse
|
38
|
Sirianni DA, Alenaizan A, Cheney DL, Sherrill CD. Assessment of Density Functional Methods for Geometry Optimization of Bimolecular van der Waals Complexes. J Chem Theory Comput 2018; 14:3004-3013. [PMID: 29763302 DOI: 10.1021/acs.jctc.8b00114] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We explore the suitability of three popular density functionals (B97-D3, B3LYP-D3, M05-2X) for producing accurate equilibrium geometries of van der Waals (vdW) complexes with diverse binding motifs. For these functionals, optimizations using Dunning's aug-cc-pVDZ basis set best combine accuracy and a reasonable computational expense. Each DFT/aug-cc-pVDZ combination produces optimized equilibrium geometries for 21 small vdW complexes of organic molecules (up to four non-hydrogen atoms total) that agree with high-level CCSD(T)/CBS reference geometries to within ±0.1 Å for the averages of the center-of-mass displacement and the mean least root-mean-squared displacement. The DFT/aug-cc-pVDZ combinations are also able to reproduce the optimal center-of-mass displacements interpolated from CCSD(T)/CBS radial potential energy surfaces in both NBC7x and HBC6 test sets to within ±0.1 Å. We therefore conclude that each of these denisty functional methods, together with the aug-cc-pVDZ basis set, is suitable for producing equilibrium geometries of generic nonbonded complexes.
Collapse
Affiliation(s)
- Dominic A Sirianni
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry , Georgia Institute of Technology , Atlanta , Georgia 30332-0400 , United States
| | - Asem Alenaizan
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry , Georgia Institute of Technology , Atlanta , Georgia 30332-0400 , United States
| | - Daniel L Cheney
- Molecular Structure and Design , Bristol-Myers Squibb Company , P.O. Box 5400, Princeton , New Jersey 08543 , United States
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering , Georgia Institute of Technology , Atlanta , Georgia 30332-0400 , United States
| |
Collapse
|
39
|
Smith DGA, Burns LA, Sirianni DA, Nascimento DR, Kumar A, James AM, Schriber JB, Zhang T, Zhang B, Abbott AS, Berquist EJ, Lechner MH, Cunha LA, Heide AG, Waldrop JM, Takeshita TY, Alenaizan A, Neuhauser D, King RA, Simmonett AC, Turney JM, Schaefer HF, Evangelista FA, DePrince AE, Crawford TD, Patkowski K, Sherrill CD. Psi4NumPy: An Interactive Quantum Chemistry Programming Environment for Reference Implementations and Rapid Development. J Chem Theory Comput 2018; 14:3504-3511. [DOI: 10.1021/acs.jctc.8b00286] [Citation(s) in RCA: 73] [Impact Index Per Article: 12.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Daniel G. A. Smith
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| | - Lori A. Burns
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| | - Dominic A. Sirianni
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| | - Daniel R. Nascimento
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, United States
| | - Ashutosh Kumar
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Andrew M. James
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Jeffrey B. Schriber
- Department of Chemistry, Emory University, Atlanta, Georgia 30322, United States
| | - Tianyuan Zhang
- Department of Chemistry, Emory University, Atlanta, Georgia 30322, United States
| | - Boyi Zhang
- Center for Computational Quantum Chemistry, University of Georgia, Athens, Georgia 30602, United States
| | - Adam S. Abbott
- Center for Computational Quantum Chemistry, University of Georgia, Athens, Georgia 30602, United States
| | - Eric J. Berquist
- University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - Marvin H. Lechner
- Department of Chemistry, Technical University of Munich, 80333 Munich, Germany
| | - Leonardo A. Cunha
- The Technical Institute of Aeronautics, São José dos Campos, 12228-900, Brazil
| | - Alexander G. Heide
- Department of Chemistry, Bethel University, St. Paul, Minnesota 55112, United States
| | - Jonathan M. Waldrop
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849, United States
| | - Tyler Y. Takeshita
- Department of Chemistry, University of California Berkeley, Berkeley, California 94720, United States
| | - Asem Alenaizan
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| | - Daniel Neuhauser
- Department of Chemistry and Biochemistry, University of California, Los Angeles, California 90095, United States
| | - Rollin A. King
- Department of Chemistry, Bethel University, St. Paul, Minnesota 55112, United States
| | - Andrew C. Simmonett
- National Institutes of Health - National Heart, Lung and Blood Institute, Laboratory of Computational Biology, 5635 Fishers Lane, T-900 Suite, Rockville, Maryland 20852, United States
| | - Justin M. Turney
- Center for Computational Quantum Chemistry, University of Georgia, Athens, Georgia 30602, United States
| | - Henry F. Schaefer
- Center for Computational Quantum Chemistry, University of Georgia, Athens, Georgia 30602, United States
| | | | - A. Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, United States
| | - T. Daniel Crawford
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Konrad Patkowski
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849, United States
| | - C. David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| |
Collapse
|
40
|
Affiliation(s)
- Ryan M. Richard
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| | - Brandon W. Bakr
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| | - C. David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| |
Collapse
|
41
|
Brahmachari U, Gonthier JF, Sherrill CD, Barry BA. Chloride Maintains a Protonated Internal Water Network in the Photosynthetic Oxygen Evolving Complex. J Phys Chem B 2017; 121:10327-10337. [DOI: 10.1021/acs.jpcb.7b08358] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
|
42
|
Moore KB, Sadeghian K, Sherrill CD, Ochsenfeld C, Schaefer HF. C-H···O Hydrogen Bonding. The Prototypical Methane-Formaldehyde System: A Critical Assessment. J Chem Theory Comput 2017; 13:5379-5395. [PMID: 29039941 DOI: 10.1021/acs.jctc.7b00753] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
Distinguishing the functionality of C-H···O hydrogen bonds (HBs) remains challenging, because their properties are difficult to quantify reliably. Herein, we present a study of the model methane-formaldehyde complex (MFC). Six stationary points on the MFC potential energy surface (PES) were obtained at the CCSD(T)/ANO2 level. The CCSDT(Q)/CBS interaction energies of the conformers range from only -1.12 kcal mol-1 to -0.33 kcal mol-1, denoting a very flat PES. Notably, only the lowest energy stationary point (MFC1) corresponds to a genuine minimum, whereas all other stationary points-including the previously studied ideal case of ae(C-H···O) = 180°-exhibit some degree of freedom that leads to MFC1. Despite the flat PES, we clearly see that the HB properties of MFC1 align with those of the prototypical water dimer O-H···O HB. Each HB property generally becomes less prominent in the higher-energy conformers. Only the MFC1 conformer prominently exhibits (1) elongated C-H donor bonds, (2) attractive C-H···O═C interactions, (3) n(O) → σ*(C-H) hyperconjugation, (4) critical points in the electron density from Bader's method and from the noncovalent interactions method, (5) positively charged donor hydrogen, and (6) downfield NMR chemical shifts and nonzero 2J(CM-HM···OF) coupling constants. Based on this research, some issues merit further study. The flat PES hinders reliable determinations of the HB-induced shifts of the C-H stretches; a similarly difficult challenge is observed for the experiment. The role of charge transfer in HBs remains an intriguing open question, although our BLW and NBO computations suggest that it is relevant to the C-H···O HB geometries. These issues notwithstanding, the prominence of the HB properties in MFC1 serves as clear evidence that the MFC is predominantly bound by a C-H···O HB.
Collapse
Affiliation(s)
- Kevin B Moore
- Center for Computational Quantum Chemistry, University of Georgia , Athens, Georgia 30602, United States
| | - Keyarash Sadeghian
- Department of Chemistry, Ludwig-Maximilians University (LMU) , Munich D-81377, Germany
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology , Atlanta, Georgia 30332, United States
| | - Christian Ochsenfeld
- Department of Chemistry, Ludwig-Maximilians University (LMU) , Munich D-81377, Germany
| | - Henry F Schaefer
- Center for Computational Quantum Chemistry, University of Georgia , Athens, Georgia 30602, United States
| |
Collapse
|
43
|
Parrish RM, Burns LA, Smith DGA, Simmonett AC, DePrince AE, Hohenstein EG, Bozkaya U, Sokolov AY, Di Remigio R, Richard RM, Gonthier JF, James AM, McAlexander HR, Kumar A, Saitow M, Wang X, Pritchard BP, Verma P, Schaefer HF, Patkowski K, King RA, Valeev EF, Evangelista FA, Turney JM, Crawford TD, Sherrill CD. Psi4 1.1: An Open-Source Electronic Structure Program Emphasizing Automation, Advanced Libraries, and Interoperability. J Chem Theory Comput 2017; 13:3185-3197. [PMID: 28489372 PMCID: PMC7495355 DOI: 10.1021/acs.jctc.7b00174] [Citation(s) in RCA: 742] [Impact Index Per Article: 106.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Psi4 is an ab initio electronic structure program providing methods such as Hartree-Fock, density functional theory, configuration interaction, and coupled-cluster theory. The 1.1 release represents a major update meant to automate complex tasks, such as geometry optimization using complete-basis-set extrapolation or focal-point methods. Conversion of the top-level code to a Python module means that Psi4 can now be used in complex workflows alongside other Python tools. Several new features have been added with the aid of libraries providing easy access to techniques such as density fitting, Cholesky decomposition, and Laplace denominators. The build system has been completely rewritten to simplify interoperability with independent, reusable software components for quantum chemistry. Finally, a wide range of new theoretical methods and analyses have been added to the code base, including functional-group and open-shell symmetry adapted perturbation theory, density-fitted coupled cluster with frozen natural orbitals, orbital-optimized perturbation and coupled-cluster methods (e.g., OO-MP2 and OO-LCCD), density-fitted multiconfigurational self-consistent field, density cumulant functional theory, algebraic-diagrammatic construction excited states, improvements to the geometry optimizer, and the "X2C" approach to relativistic corrections, among many other improvements.
Collapse
Affiliation(s)
- Robert M Parrish
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology , Atlanta, Georgia 30332-0400, United States
| | - Lori A Burns
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology , Atlanta, Georgia 30332-0400, United States
| | - Daniel G A Smith
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology , Atlanta, Georgia 30332-0400, United States
| | - Andrew C Simmonett
- National Institutes of Health , National Heart, Lung and Blood Institute, Laboratory of Computational Biology, 5635 Fishers Lane, T-900 Suite, Rockville, Maryland 20852, United States
| | - A Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University , Tallahassee, Florida 32306-4390, United States
| | - Edward G Hohenstein
- Department of Chemistry and Biochemistry, The City College of New York , New York, New York 10031, United States
| | - Uğur Bozkaya
- Department of Chemistry, Hacettepe University , Ankara 06800, Turkey
| | - Alexander Yu Sokolov
- Division of Chemistry and Chemical Engineering, California Institute of Technology , Pasadena, California 91125, United States
| | - Roberto Di Remigio
- Department of Chemistry, Centre for Theoretical and Computational Chemistry, UiT, The Arctic University of Norway , N-9037 Tromsø, Norway
| | - Ryan M Richard
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology , Atlanta, Georgia 30332-0400, United States
| | - Jérôme F Gonthier
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology , Atlanta, Georgia 30332-0400, United States
| | - Andrew M James
- Department of Chemistry, Virginia Tech , Blacksburg, Virginia 24061, United States
| | - Harley R McAlexander
- Department of Chemistry, Virginia Tech , Blacksburg, Virginia 24061, United States
| | - Ashutosh Kumar
- Department of Chemistry, Virginia Tech , Blacksburg, Virginia 24061, United States
| | - Masaaki Saitow
- Department of Chemistry and Research Center for Smart Molecules, Rikkyo University , 3-34-1 Nishi-ikebukuro, Toshima-ku, Tokyo 171-8501, Japan
| | - Xiao Wang
- Department of Chemistry, Virginia Tech , Blacksburg, Virginia 24061, United States
| | - Benjamin P Pritchard
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology , Atlanta, Georgia 30332-0400, United States
| | - Prakash Verma
- Department of Chemistry, Emory University , Atlanta, Georgia 30322, United States
| | - Henry F Schaefer
- Center for Computational Quantum Chemistry, University of Georgia , Athens, Georgia 30602, United States
| | - Konrad Patkowski
- Department of Chemistry and Biochemistry, Auburn University , Auburn, Alabama 36849, United States
| | - Rollin A King
- Department of Chemistry, Bethel University , St. Paul, Minnesota 55112, United States
| | - Edward F Valeev
- Department of Chemistry, Virginia Tech , Blacksburg, Virginia 24061, United States
| | | | - Justin M Turney
- Center for Computational Quantum Chemistry, University of Georgia , Athens, Georgia 30602, United States
| | - T Daniel Crawford
- Department of Chemistry, Virginia Tech , Blacksburg, Virginia 24061, United States
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology , Atlanta, Georgia 30332-0400, United States
| |
Collapse
|
44
|
Parrish RM, Sitkoff DF, Cheney DL, Sherrill CD. The Surprising Importance of Peptide Bond Contacts in Drug–Protein Interactions. Chemistry 2017; 23:7887-7890. [DOI: 10.1002/chem.201701031] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2017] [Indexed: 01/08/2023]
Affiliation(s)
- Robert M. Parrish
- Center for Computational Molecular Science and Technology School of Chemistry and Biochemistry School of Computational Science and Engineering Georgia Institute of Technology Atlanta GA 30332-0400 USA
| | - Doree F. Sitkoff
- Molecular Structure and Design Bristol-Myers Squibb Company 311 Pennington-Rocky Hill Road Pennington NJ 08534 USA
| | - Daniel L. Cheney
- Molecular Structure and Design Bristol-Myers Squibb Company 311 Pennington-Rocky Hill Road Pennington NJ 08534 USA
| | - C. David Sherrill
- Center for Computational Molecular Science and Technology School of Chemistry and Biochemistry School of Computational Science and Engineering Georgia Institute of Technology Atlanta GA 30332-0400 USA
| |
Collapse
|
45
|
Gonthier JF, Sherrill CD. Density-fitted open-shell symmetry-adapted perturbation theory and application to π-stacking in benzene dimer cation and ionized DNA base pair steps. J Chem Phys 2017; 145:134106. [PMID: 27782424 DOI: 10.1063/1.4963385] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023] Open
Abstract
Symmetry-Adapted Perturbation Theory (SAPT) is one of the most popular approaches to energy component analysis of non-covalent interactions between closed-shell systems, yielding both accurate interaction energies and meaningful interaction energy components. In recent years, the full open-shell equations for SAPT up to second-order in the intermolecular interaction and zeroth-order in the intramolecular correlation (SAPT0) were published [P. S. Zuchowski et al., J. Chem. Phys. 129, 084101 (2008); M. Hapka et al., ibid. 137, 164104 (2012)]. Here, we utilize density-fitted electron repulsion integrals to produce an efficient computational implementation. This approach is used to examine the effect of ionization on π-π interactions. For the benzene dimer radical cation, comparison against reference values indicates a good performance for open-shell SAPT0, except in cases with substantial charge transfer. For π stacking between hydrogen-bonded pairs of nucleobases, dispersion interactions still dominate binding, in spite of the creation of a positive charge.
Collapse
Affiliation(s)
- Jérôme F Gonthier
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - C David Sherrill
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| |
Collapse
|
46
|
Bozkaya U, Sherrill CD. Analytic energy gradients for the coupled-cluster singles and doubles method with the density-fitting approximation. J Chem Phys 2017; 144:174103. [PMID: 27155621 DOI: 10.1063/1.4948318] [Citation(s) in RCA: 30] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
An efficient implementation is presented for analytic gradients of the coupled-cluster singles and doubles (CCSD) method with the density-fitting approximation, denoted DF-CCSD. Frozen core terms are also included. When applied to a set of alkanes, the DF-CCSD analytic gradients are significantly accelerated compared to conventional CCSD for larger molecules. The efficiency of our DF-CCSD algorithm arises from the acceleration of several different terms, which are designated as the "gradient terms": computation of particle density matrices (PDMs), generalized Fock-matrix (GFM), solution of the Z-vector equation, formation of the relaxed PDMs and GFM, back-transformation of PDMs and GFM to the atomic orbital (AO) basis, and evaluation of gradients in the AO basis. For the largest member of the alkane set (C10H22), the computational times for the gradient terms (with the cc-pVTZ basis set) are 2582.6 (CCSD) and 310.7 (DF-CCSD) min, respectively, a speed up of more than 8-folds. For gradient related terms, the DF approach avoids the usage of four-index electron repulsion integrals. Based on our previous study [U. Bozkaya, J. Chem. Phys. 141, 124108 (2014)], our formalism completely avoids construction or storage of the 4-index two-particle density matrix (TPDM), using instead 2- and 3-index TPDMs. The DF approach introduces negligible errors for equilibrium bond lengths and harmonic vibrational frequencies.
Collapse
Affiliation(s)
- Uğur Bozkaya
- Department of Chemistry, Hacettepe University, Ankara 06800, Turkey
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| |
Collapse
|
47
|
Dolgounitcheva O, Díaz-Tinoco M, Zakrzewski VG, Richard RM, Marom N, Sherrill CD, Ortiz JV. Correction to Accurate Ionization Potentials and Electron Affinities of Acceptor Molecules IV: Electron-Propagator Methods. J Chem Theory Comput 2016; 13:389-391. [DOI: 10.1021/acs.jctc.6b01180] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
48
|
Sirianni DA, Burns LA, Sherrill CD. Comparison of Explicitly Correlated Methods for Computing High-Accuracy Benchmark Energies for Noncovalent Interactions. J Chem Theory Comput 2016; 13:86-99. [DOI: 10.1021/acs.jctc.6b00797] [Citation(s) in RCA: 43] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Dominic A. Sirianni
- Center for Computational
Molecular Science and Technology, School of Chemistry and Biochemistry,
School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| | - Lori A. Burns
- Center for Computational
Molecular Science and Technology, School of Chemistry and Biochemistry,
School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| | - C. David Sherrill
- Center for Computational
Molecular Science and Technology, School of Chemistry and Biochemistry,
School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, United States
| |
Collapse
|
49
|
Abstract
Since the original fitting of Grimme's DFT-D3 damping parameters, the number and quality of benchmark interaction energies has increased significantly. Here, conventional benchmark sets, which focus on minimum-orientation radial curves at the expense of angular diversity, are augmented by new databases such as side chain-side chain interactions (SSI), which are composed of interactions gleaned from crystal data and contain no such minima-focused bias. Moreover, some existing databases such as S22×5 are extended to shorter intermolecular separations. This improved DFT-D3 training set provides a balanced description of distances, covers the entire range of interaction types, and at 1526 data points is far larger than the original training set of 130. The results are validated against a new collection of 6773 data points and demonstrate that the effect of refitting the damping parameters ranges from no change in accuracy (LC-ωPBE-D3) to an almost 2-fold decrease in average error (PBE-D3).
Collapse
Affiliation(s)
- Daniel G A Smith
- Department of Chemistry and Biochemistry, Auburn University , Auburn, Alabama 36849, United States
| | - Lori A Burns
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology , Atlanta, Georgia 30332-0400, United States
| | - Konrad Patkowski
- Department of Chemistry and Biochemistry, Auburn University , Auburn, Alabama 36849, United States
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry and School of Computational Science and Engineering, Georgia Institute of Technology , Atlanta, Georgia 30332-0400, United States
| |
Collapse
|
50
|
Bakr BW, Sherrill CD. Analysis of transition state stabilization by non-covalent interactions in the Houk-List model of organocatalyzed intermolecular Aldol additions using functional-group symmetry-adapted perturbation theory. Phys Chem Chem Phys 2016; 18:10297-308. [PMID: 27020417 DOI: 10.1039/c5cp07281f] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Rational design of catalysts would be aided by a better understanding of how non-covalent interactions stabilize transition states. Here, we apply the newly-developed Functional-Group Symmetry-Adapted Perturbation Theory (F-SAPT) to quantify non-covalent interactions in transition states of the proline-catalyzed intermolecular aldol reaction between benzaldehyde and cyclohexanone, according to the Houk-List mechanism [Bahmanyar et al., J. Am. Chem. Soc., 2003, 125, 2475]. A recent re-examination of this organocatalytic reaction by Rzepa and co-workers [Armstrong et al., Chem. Sci., 2014, 5, 2057] used electron density analysis to identify three key non-covalent interactions thought to influence stereoselectivity: (1) a favorable electrostatic interaction (originally identified by Houk and List) between the NCH(δ+) group of the enamine intermediate and the (δ-)O[double bond, length as m-dash]C of benzaldehyde; (2) a C-H/π interaction between the cyclohexene group of the enamine intermediate and the benzaldehyde phenyl ring; (3) a stabilizing contact between an ortho-hydrogen of the phenyl and an oxygen of the carboxylic acid group of the enamine. These three interactions have been directly computed using F-SAPT, which confirms the stabilizing interaction between an ortho-hydrogen and the carboxylic acid in the (S,S) and (R,S) transition state stereoisomers. F-SAPT analysis also finds stabilizing dispersion and electrostatic interactions due to a C-H/π interaction between the cyclohexene and phenyl groups in the (S,S) and (R,R) transition states. However, unfavorable exchange-repulsion cancels the attractive terms that favor these stereoisomers. Surprisingly, the interaction thought to be most important for stereoselectivity, the NCH(δ+)(δ-)O[double bond, length as m-dash]C interaction, is actually found to be repulsive due to the negative charge on the nitrogen. Hence, our results indicate that geometric analysis and/or density-based analysis does not necessarily produce a reliable picture of non-covalent stabilization. As confirmed by high-level coupled-cluster computations, intermolecular interaction energies are strongest for the (R,R) transition states, which are not the experimentally favored products. This suggests that at least for this reaction, stereoselectivity is also strongly dependent on the energy required to distort the reacting molecules into the transition state geometry.
Collapse
Affiliation(s)
- Brandon W Bakr
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, GA 30332-0400, USA.
| | | |
Collapse
|