1
|
Deng Z, Liu C, Li Z, Zhang Y. An efficient method by combining different basis sets and SAPT levels. J Comput Chem 2024; 45:1936-1944. [PMID: 38703182 DOI: 10.1002/jcc.27386] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2023] [Revised: 04/15/2024] [Accepted: 04/18/2024] [Indexed: 05/06/2024]
Abstract
In symmetry-adapted perturbation theory (SAPT), accurate calculations on non-covalent interaction (NCI) for large complexes with more than 50 atoms are time-consuming using large basis sets. More efficient ones with smaller basis sets usually result in poor prediction in terms of dispersion and overall energies. In this study, we propose two composite methods with baseline calculated at SAPT2/aug-cc-pVDZ and SAPT2/aug-cc-pVTZ with dispersion term corrected at SAPT2+ level using bond functions and smaller basis set with δ MP2 corrections respectively. Benchmark results on representative NCI data sets, such as S22, S66, and so forth, show significant improvements on the accuracy compared to the original SAPT Silver standard and comparable to SAPT Gold standard in some cases with much less computational cost.
Collapse
Affiliation(s)
- Zhihao Deng
- Beijing StoneWise Technology Co Ltd., Beijing, China
| | - Chang Liu
- Beijing StoneWise Technology Co Ltd., Beijing, China
| | - Zhongwei Li
- Yantai Gogetter Technology Co Ltd., Yantai, China
| | | |
Collapse
|
2
|
Szczygiel M, Derewenda U, Scheiner S, Minor W, Derewenda ZS. A structural role for tryptophan in proteins, and the ubiquitous Trp C δ1-H...O=C (backbone) hydrogen bond. Acta Crystallogr D Struct Biol 2024; 80:551-562. [PMID: 38941144 PMCID: PMC11220837 DOI: 10.1107/s2059798324005515] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2024] [Accepted: 06/09/2024] [Indexed: 06/29/2024] Open
Abstract
Tryptophan is the most prominent amino acid found in proteins, with multiple functional roles. Its side chain is made up of the hydrophobic indole moiety, with two groups that act as donors in hydrogen bonds: the Nϵ-H group, which is a potent donor in canonical hydrogen bonds, and a polarized Cδ1-H group, which is capable of forming weaker, noncanonical hydrogen bonds. Due to adjacent electron-withdrawing moieties, C-H...O hydrogen bonds are ubiquitous in macromolecules, albeit contingent on the polarization of the donor C-H group. Consequently, Cα-H groups (adjacent to the carbonyl and amino groups of flanking peptide bonds), as well as the Cϵ1-H and Cδ2-H groups of histidines (adjacent to imidazole N atoms), are known to serve as donors in hydrogen bonds, for example stabilizing parallel and antiparallel β-sheets. However, the nature and the functional role of interactions involving the Cδ1-H group of the indole ring of tryptophan are not well characterized. Here, data mining of high-resolution (r ≤ 1.5 Å) crystal structures from the Protein Data Bank was performed and ubiquitous close contacts between the Cδ1-H groups of tryptophan and a range of electronegative acceptors were identified, specifically main-chain carbonyl O atoms immediately upstream and downstream in the polypeptide chain. The stereochemical analysis shows that most of the interactions bear all of the hallmarks of proper hydrogen bonds. At the same time, their cohesive nature is confirmed by quantum-chemical calculations, which reveal interaction energies of 1.5-3.0 kcal mol-1, depending on the specific stereochemistry.
Collapse
Affiliation(s)
- Michal Szczygiel
- Department of Molecular Physiology and Biological PhysicsUniversity of Virginia1340 Jefferson Park AvenueCharlottesvilleVA22908-0736USA
| | - Urszula Derewenda
- Department of Molecular Physiology and Biological PhysicsUniversity of Virginia1340 Jefferson Park AvenueCharlottesvilleVA22908-0736USA
| | - Steve Scheiner
- Department of Chemistry and BiochemistryUtah State UniversityLoganUtahUSA
| | - Wladek Minor
- Department of Molecular Physiology and Biological PhysicsUniversity of Virginia1340 Jefferson Park AvenueCharlottesvilleVA22908-0736USA
| | - Zygmunt S. Derewenda
- Department of Molecular Physiology and Biological PhysicsUniversity of Virginia1340 Jefferson Park AvenueCharlottesvilleVA22908-0736USA
| |
Collapse
|
3
|
Scheiner S, Amonov A. Types of noncovalent bonds within complexes of thiazole with CF 4 and SiF 4. Phys Chem Chem Phys 2024; 26:6127-6137. [PMID: 38299682 DOI: 10.1039/d4cp00057a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2024]
Abstract
The five-membered heteroaromatic thiazole molecule contains a number of electron-rich regions that could attract an electrophile, namely the N and S lone pairs that lie in the molecular plane, and π-system areas above the plane. The possibility of each of these sites engaging in a tetrel bond (TB) with CF4 and SiF4, as well as geometries that encompass a CH⋯F H-bond, was explored via DFT calculations. There are a number of minima that occur in the pairing of thiazole with CF4 that are very close in energy, but these complexes are weakly bound by less than 2 kcal mol-1 and the presence of a true TB is questionable. The inclusion of zero-point vibrational energies alters the energetic ordering, which is further modified when entropic effects are added. The preferred geometry would thus be sensitive to the temperature of an experiment. Replacement of CF4 by SiF4 leaves intact most of the configurations, and their tight energetic clustering, the ordering of which is again altered as the temperature rises. But there is one exception in that by far the most tightly bound complex involves a strong Si⋯N TB between SiF4 and the lone pair of the thiazole N, with an interaction energy of 30 kcal mol-1. Even accounting for its high deformation energy and entropic considerations, this structure remains as clearly the most stable at any temperature.
Collapse
Affiliation(s)
- Steve Scheiner
- Department of Chemistry and Biochemistry Utah State University Logan, Utah 84322-0300, USA.
| | - Akhtam Amonov
- Department of Optics and Spectroscopy, Institute of Engineering Physics Samarkand State University 140104, University blv. 15, Samarkand, Uzbekistan
| |
Collapse
|
4
|
Ali Z, Kukhta T, Trant JF, Sharma P. An Atlas of the base inter-RNA stacks involved in bacterial translation. Biophys Chem 2024; 305:107144. [PMID: 38061282 DOI: 10.1016/j.bpc.2023.107144] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2023] [Accepted: 11/28/2023] [Indexed: 01/03/2024]
Abstract
Nucleobase-specific noncovalent interactions play a crucial role in translation. Herein, we provide a comprehensive analysis of the stacks between different RNA components in the crystal structures of the bacterial ribosome caught at different translation stages. Analysis of tRNA||rRNA stacks reveals distinct behaviour; both the A-and E-site tRNAs exhibit unique stacking patterns with 23S rRNA bases, while P-site tRNAs stack with 16S rRNA bases. Furthermore, E-site stacks exhibit diverse face orientations and ring topologies-rare for inter-chain RNA interactions-with higher average interaction energies than A or P-site stacks. This suggests that stacking may be essential for stabilizing tRNA progression through the E-site. Additionally, mRNA||rRNA stacks reveal other geometries, which depend on the tRNA binding site, whereas 16S rRNA||23S rRNA stacks highlight the importance of specific bases in maintaining the integrity of the translational complex by linking the two rRNAs. Furthermore, tRNA||mRNA stacks exhibit distinct geometries and energetics at the E-site, indicating their significance during tRNA translocation and elimination. Overall, both A and E-sites display a more diverse distribution of inter-RNA stacks compared to the P-site. Stacking interactions in the active ribosome are not simply accidental byproducts of biochemistry but are likely invoked to compensate and support the integrity and dynamics of translation.
Collapse
MESH Headings
- RNA, Ribosomal, 16S/chemistry
- RNA, Ribosomal, 16S/genetics
- RNA, Ribosomal, 16S/metabolism
- RNA, Ribosomal, 23S/chemistry
- RNA, Ribosomal, 23S/genetics
- RNA, Ribosomal, 23S/metabolism
- Ribosomes/chemistry
- RNA, Transfer/chemistry
- RNA, Transfer/genetics
- RNA, Transfer/metabolism
- RNA, Messenger/genetics
- RNA, Messenger/analysis
- RNA, Messenger/metabolism
- Nucleic Acid Conformation
Collapse
Affiliation(s)
- Zakir Ali
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh 160014, India
| | - Teagan Kukhta
- Department of Chemistry and Biochemistry, University of Windsor, 401 Sunset Ave., Windsor, ON N9B 3P4, Canada
| | - John F Trant
- Department of Chemistry and Biochemistry, University of Windsor, 401 Sunset Ave., Windsor, ON N9B 3P4, Canada; We-Spark Health Institute, 401 Sunset Ave., Windsor, ON N9B 3P4, Canada; Binary Star Research Services, LaSalle, ON N9J 3X8, Canada.
| | - Purshotam Sharma
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh 160014, India; Department of Chemistry and Biochemistry, University of Windsor, 401 Sunset Ave., Windsor, ON N9B 3P4, Canada.
| |
Collapse
|
5
|
Thürlemann M, Riniker S. Hybrid classical/machine-learning force fields for the accurate description of molecular condensed-phase systems. Chem Sci 2023; 14:12661-12675. [PMID: 38020395 PMCID: PMC10646964 DOI: 10.1039/d3sc04317g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2023] [Accepted: 10/24/2023] [Indexed: 12/01/2023] Open
Abstract
Electronic structure methods offer in principle accurate predictions of molecular properties, however, their applicability is limited by computational costs. Empirical methods are cheaper, but come with inherent approximations and are dependent on the quality and quantity of training data. The rise of machine learning (ML) force fields (FFs) exacerbates limitations related to training data even further, especially for condensed-phase systems for which the generation of large and high-quality training datasets is difficult. Here, we propose a hybrid ML/classical FF model that is parametrized exclusively on high-quality ab initio data of dimers and monomers in vacuum but is transferable to condensed-phase systems. The proposed hybrid model combines our previous ML-parametrized classical model with ML corrections for situations where classical approximations break down, thus combining the robustness and efficiency of classical FFs with the flexibility of ML. Extensive validation on benchmarking datasets and experimental condensed-phase data, including organic liquids and small-molecule crystal structures, showcases how the proposed approach may promote FF development and unlock the full potential of classical FFs.
Collapse
Affiliation(s)
- Moritz Thürlemann
- Department of Chemistry and Applied Biosciences, ETH Zürich Vladimir-Prelog-Weg 2 Zürich 8093 Switzerland
| | - Sereina Riniker
- Department of Chemistry and Applied Biosciences, ETH Zürich Vladimir-Prelog-Weg 2 Zürich 8093 Switzerland
| |
Collapse
|
6
|
Amonov A, Scheiner S. Relation between Halogen Bond Strength and IR and NMR Spectroscopic Markers. Molecules 2023; 28:7520. [PMID: 38005241 PMCID: PMC10673387 DOI: 10.3390/molecules28227520] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2023] [Revised: 11/04/2023] [Accepted: 11/08/2023] [Indexed: 11/26/2023] Open
Abstract
The relationship between the strength of a halogen bond (XB) and various IR and NMR spectroscopic quantities is assessed through DFT calculations. Three different Lewis acids place a Br or I atom on a phenyl ring; each is paired with a collection of N and O bases of varying electron donor power. The weakest of the XBs display a C-X bond contraction coupled with a blue shift in the associated frequency, whereas the reverse trends occur for the stronger bonds. The best correlations with the XB interaction energy are observed with the NMR shielding of the C atom directly bonded to X and the coupling constants involving the C-X bond and the C-H/F bond that lies ortho to the X substituent, but these correlations are not accurate enough for the quantitative assessment of energy. These correlations tend to improve as the Lewis acid becomes more potent, which makes for a wider range of XB strengths.
Collapse
Affiliation(s)
- Akhtam Amonov
- Department of Optics and Spectroscopy, Engineering Physics Institute, Samarkand State University, University blv. 15, Samarkand 140104, Uzbekistan;
| | - Steve Scheiner
- Department of Chemistry and Biochemistry, Utah State University, Logan, UT 84322-0300, USA
| |
Collapse
|
7
|
Su X, Li S, Yan L. N-H Bond Activation Catalyzed by an Anderson-Type Polyoxometalate-Based Compound: Key Role of Transition-Metal Heteroatom. Inorg Chem 2023; 62:15673-15679. [PMID: 37708077 DOI: 10.1021/acs.inorgchem.3c02462] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/16/2023]
Abstract
Polyoxometalates (POMs) have a broad array of applied platforms with well-characterized catalysis to achieve N-H bond activation. Herein, the mechanism of the Anderson-type POM-based catalyst [FeIIIMoVI6O18{(OCH2)3CNH2}2]3- ([TrisFeIIIMoVI6O18]3-, Tris = {(OCH2)3CNH2}2) for the N-H bond activation of hydrazine (PhHNNHPh) was investigated by density functional theory calculations. The results reveal that [TrisFeIIIMoVI6O18]3- as the active species is responsible for the continuous abstraction of two electrons and two protons of PhHNNHPh via a proton-coupled electron transfer pathway, resulting in the activation of two N-H bonds in PhHNNHPh and thus the product PhNNPh. H2O2 acts as an oxidant to regulate catalyst regeneration. Based on the proposed catalytic mechanism, the key role of the heteroatom FeIII in [TrisFeIIIMoVI6O18]3- was disclosed. The d-orbital of FeIII in [TrisFeIIIMoVI6O18]3- acts as an electron receptor to promote the electron transfer (ET) in the rate-determining step (RDS) of the catalytic cycle. The substitution of the heteroatom FeIII of [TrisFeIIIMoVI6O18]3- with CoIII, RuIII, or MnIII is expected to improve the catalytic activity for several reasons: (i) the unoccupied molecular orbitals of POM-based compounds containing CoIII or RuIII are low, which is beneficial for the ET of RDS; (ii) For N-H bond activation catalyzed by the MnIII-containing POM-based compound, the transition state of RDS is stable because the d-orbital of its active site is half-filled, which results in a low free-energy barrier.
Collapse
Affiliation(s)
- Xiaofang Su
- Collaborative Innovation Centre of Henan Province for Green Manufacturing of Fine Chemicals, Key Laboratory of Green Chemical Media and Reactions, Ministry of Education, School of Chemistry and Chemical Engineering, Henan Normal University, Xinxiang 453007, Henan, P. R. China
| | - Shujun Li
- Collaborative Innovation Centre of Henan Province for Green Manufacturing of Fine Chemicals, Key Laboratory of Green Chemical Media and Reactions, Ministry of Education, School of Chemistry and Chemical Engineering, Henan Normal University, Xinxiang 453007, Henan, P. R. China
| | - Likai Yan
- Key Laboratory of Polyoxometalate Science of Ministry of Education, Institute of Functional Material Chemistry, Faculty of Chemistry, Northeast Normal University, Changchun 130018, P. R. China
| |
Collapse
|
8
|
Wang Y, Xu C, Li Z, Barati Farimani A. Denoise Pretraining on Nonequilibrium Molecules for Accurate and Transferable Neural Potentials. J Chem Theory Comput 2023; 19:5077-5087. [PMID: 37390120 PMCID: PMC10413865 DOI: 10.1021/acs.jctc.3c00289] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2023] [Indexed: 07/02/2023]
Abstract
Recent advances in equivariant graph neural networks (GNNs) have made deep learning amenable to developing fast surrogate models to expensive ab initio quantum mechanics (QM) approaches for molecular potential predictions. However, building accurate and transferable potential models using GNNs remains challenging, as the data are greatly limited by the expensive computational costs and level of theory of QM methods, especially for large and complex molecular systems. In this work, we propose denoise pretraining on nonequilibrium molecular conformations to achieve more accurate and transferable GNN potential predictions. Specifically, atomic coordinates of sampled nonequilibrium conformations are perturbed by random noises, and GNNs are pretrained to denoise the perturbed molecular conformations which recovers the original coordinates. Rigorous experiments on multiple benchmarks reveal that pretraining significantly improves the accuracy of neural potentials. Furthermore, we show that the proposed pretraining approach is model-agnostic, as it improves the performance of different invariant and equivariant GNNs. Notably, our models pretrained on small molecules demonstrate remarkable transferability, improving performance when fine-tuned on diverse molecular systems, including different elements, charged molecules, biomolecules, and larger systems. These results highlight the potential for leveraging denoise pretraining approaches to build more generalizable neural potentials for complex molecular systems.
Collapse
Affiliation(s)
- Yuyang Wang
- Department
of Mechanical Engineering, Carnegie Mellon
University, Pittsburgh, Pennsylvania 15213, United States
- Machine
Learning Department, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States
| | - Changwen Xu
- Department
of Mechanical Engineering, Carnegie Mellon
University, Pittsburgh, Pennsylvania 15213, United States
| | - Zijie Li
- Department
of Mechanical Engineering, Carnegie Mellon
University, Pittsburgh, Pennsylvania 15213, United States
| | - Amir Barati Farimani
- Department
of Mechanical Engineering, Carnegie Mellon
University, Pittsburgh, Pennsylvania 15213, United States
- Machine
Learning Department, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States
- Department
of Materials Science and Engineering, Carnegie
Mellon University, Pittsburgh, Pennsylvania 15213, United States
- Department
of Chemical Engineering, Carnegie Mellon
University, Pittsburgh, Pennsylvania 15213, United States
| |
Collapse
|
9
|
Ali Z, Kaur S, Kukhta T, Abu-Saleh AAAA, Jhunjhunwala A, Mitra A, Trant JF, Sharma P. Structural Mapping of the Base Stacks Containing Post-transcriptionally Modified Bases in RNA. J Phys Chem B 2023. [PMID: 37369074 DOI: 10.1021/acs.jpcb.3c02268] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/29/2023]
Abstract
Post-transcriptionally modified bases play vital roles in many biochemical processes involving RNA. Analysis of the non-covalent interactions associated with these bases in RNA is crucial for providing a more complete understanding of the RNA structure and function; however, the characterization of these interactions remains understudied. To address this limitation, we present a comprehensive analysis of base stacks involving all crystallographic occurrences of the most biologically relevant modified bases in a large dataset of high-resolution RNA crystal structures. This is accompanied by a geometrical classification of the stacking contacts using our established tools. Coupled with quantum chemical calculations and an analysis of the specific structural context of these stacks, this provides a map of the stacking conformations available to modified bases in RNA. Overall, our analysis is expected to facilitate structural research on altered RNA bases.
Collapse
Affiliation(s)
- Zakir Ali
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh 160014, India
| | - Sarabjeet Kaur
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh 160014, India
- Surface Chemistry and Catalysis: Characterisation and Application Team (COK-KAT), Leuven (Arenberg) Celestijnenlaan 200f─Box 2461, 3001 Leuven, Belgium
| | - Teagan Kukhta
- Department of Chemistry and Biochemistry, University of Windsor, 401 Sunset Avenue, Windsor, Ontario N9B 3P4, Canada
| | - Abd Al-Aziz A Abu-Saleh
- Department of Chemistry and Biochemistry, University of Windsor, 401 Sunset Avenue, Windsor, Ontario N9B 3P4, Canada
- Binary Star Research Services, LaSalle, Ontario N9J 3X8, Canada
| | - Ayush Jhunjhunwala
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology (IIIT-H), Gachibowli, Hyderabad, Telangana 500032, India
| | - Abhijit Mitra
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology (IIIT-H), Gachibowli, Hyderabad, Telangana 500032, India
| | - John F Trant
- Department of Chemistry and Biochemistry, University of Windsor, 401 Sunset Avenue, Windsor, Ontario N9B 3P4, Canada
- Binary Star Research Services, LaSalle, Ontario N9J 3X8, Canada
| | - Purshotam Sharma
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh 160014, India
- Department of Chemistry and Biochemistry, University of Windsor, 401 Sunset Avenue, Windsor, Ontario N9B 3P4, Canada
| |
Collapse
|
10
|
Abstract
Quantum calculations study the potential of an intramolecular H-bond between the halogen atom (X) of a halobenzene and a substituent placed ortho to it, to amplify the ability of X to engage in a halogen bond (XB) with a Lewis base. H-bonding substituents NH2, CH2CH2OH, CH2OH, OH, and COOH were added to halobenzenes (X = Cl, Br, I). The amino group had little effect, but those containing OH increased the CX···N XB energy to a NH3 nucleophile by about 0.5 kcal/mol; the increment associated with COOH is larger, nearly 2 kcal/mol. These energy increments were approximately doubled if two such H-bonding substituents are present. Combining a pair of ortho COOH groups with an electron-withdrawing NO2 group in the para position has a particularly large effect, raising the XB energy by about 4 kcal/mol, which can amount to as much as a 4-fold magnification.
Collapse
Affiliation(s)
- Steve Scheiner
- Department of Chemistry and Biochemistry, Utah State University, Logan, Utah 84322-0300, United States
| |
Collapse
|
11
|
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: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [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
|
12
|
Jaegers NR, Iglesia E. Theoretical Assessment of the Mechanism and Active Sites in Alkene Dimerization on Ni Monomers Grafted onto Aluminosilicates: (Ni-OH) + Centers and C-C Coupling Mediated by Lewis Acid-Base Pairs. J Am Chem Soc 2023; 145:6349-6361. [PMID: 36914428 DOI: 10.1021/jacs.2c13487] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/16/2023]
Abstract
Ni-based solids are effective catalysts for alkene dimerization, but the nature of active centers and identity and kinetic relevance of bound species and elementary reactions remain speculative and based on organometallic chemistry. Ni centers grafted onto ordered MCM-41 mesopores lead to well-defined monomers that are rendered stable by the presence of an intrapore nonpolar liquid, thus enabling accurate experimental inquiries and indirect evidence for grafted (Ni-OH)+ monomers. Density functional theory (DFT) treatments presented here confirm the plausible involvement of pathways and active centers not previously considered as mediators of high turnover rates for C2-C4 alkenes at cryogenic temperatures. (Ni-OH)+ species act as Lewis acid-base pairs that stabilize C-C coupling transition states by polarizing two alkenes in opposite directions via concerted interactions with the O and H atoms in these pairs. DFT-derived activation barriers for ethene dimerization (59 kJ mol-1) are similar to measured values (46 ± 5 kJ mol-1) and the weak binding of ethene on (Ni-OH)+ is consistent with kinetic trends that require sites to remain essentially bare at subambient temperatures and high alkene pressures (1-15 bar). DFT treatments of classical metallacycle and Cossee-Arlman dimerization routes (Ni+ and Ni2+-H grafted onto Al-MCM-41, respectively) show that such sites bind ethene strongly and lead to saturation coverages, in contradiction with observed kinetic trends. These C-C coupling routes at acid-base pairs in (Ni-OH)+ differ from molecular catalysts in (i) the type of elementary steps; (ii) the nature of active centers; and (iii) their catalytic competence at subambient temperatures without requiring co-catalysts or activators.
Collapse
Affiliation(s)
- Nicholas R Jaegers
- Department of Chemical and Biomolecular Engineering, University of California at Berkeley, Berkeley, California 94720, United States
| | - Enrique Iglesia
- Department of Chemical and Biomolecular Engineering, University of California at Berkeley, Berkeley, California 94720, United States
| |
Collapse
|
13
|
Zhao H, Tang Z, He M, Yang X, Lai S, An K, Han S, Qu Z, Zhou W, Wang Z. Effect of oxygen functional groups on competitive adsorption of benzene and water on carbon materials: Density functional theory study. THE SCIENCE OF THE TOTAL ENVIRONMENT 2023; 863:160772. [PMID: 36513224 DOI: 10.1016/j.scitotenv.2022.160772] [Citation(s) in RCA: 10] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/20/2022] [Revised: 11/16/2022] [Accepted: 12/04/2022] [Indexed: 06/17/2023]
Abstract
It is important to study the effect of oxygen-containing functional groups on the competitive adsorption mechanism of benzene and water on the surface of carbon materials, and to directional modification of activated carbon to improve its selective adsorption of benzene in air. In this study, the adsorption characteristics of benzene and water on original and linked ester, carboxyl, hydroxyl, carbon materials linked by ether groups were calculated by quantum chemical simulation based on density functional theory. The types and proportions of weak interactions in the adsorption process were calculated by energy decomposition analysis, and the adsorption mechanism of carbon materials for water and benzene was described. The influence and contribution of oxygen-containing functional groups on the adsorption of benzene and water were further analyzed by van der Waals potential and electrostatic potential, respectively, so as to determine the difference in the adsorption effect of different types of oxygen-containing functional groups on the two molecules. It was found that the carboxyl group has a great influence on the hydrophilicity of carbon materials, and the electrostatic potential distribution before and after linking the carboxyl group changed significantly. Therefore, they can attract each other with water through hydrogen bonds and occupy the surface adsorption sites of carbon materials, thereby inhibiting the adsorption of benzene on carbon materials. On the contrary, due to its hydrophobic properties, the ether group will free up adsorption space for the adsorption of benzene on the surface of the carbon material, which is beneficial to the adsorption of benzene. The adsorption experiments were carried out, and the results were consistent with the simulation. This study provides an idea for preparing efficient carbonaceous adsorbent of benzene and reducing benzene pollution in industry.
Collapse
Affiliation(s)
- Haiqian Zhao
- School of Mechanical Science and Engineering, Northeast Petroleum University, Daqing 163318, China
| | - Ziyu Tang
- School of Civil and Architectural Engineering, Northeast Petroleum University, Daqing 163318, China
| | - Mingqi He
- School of Mechanical Science and Engineering, Northeast Petroleum University, Daqing 163318, China
| | - Xue Yang
- School of Mechanical Science and Engineering, Northeast Petroleum University, Daqing 163318, China
| | - Shiwei Lai
- School of Civil and Architectural Engineering, Northeast Petroleum University, Daqing 163318, China
| | - Kaibo An
- School of Civil and Architectural Engineering, Northeast Petroleum University, Daqing 163318, China
| | - Shuaishuai Han
- School of Civil and Architectural Engineering, Northeast Petroleum University, Daqing 163318, China
| | - Zhibin Qu
- School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, China
| | - Wei Zhou
- School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, China
| | - Zhonghua Wang
- School of Civil and Architectural Engineering, Northeast Petroleum University, Daqing 163318, China.
| |
Collapse
|
14
|
Xie J, Xi Z, Yang Z, Zhang X, Yuan H, Yang Y, Ni L, He M. Computational Approach to Understanding the Structures, Properties, and Supramolecular Chemistry of Pagoda[n]arenes. J Mol Struct 2023. [DOI: 10.1016/j.molstruc.2023.135073] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
|
15
|
Kadhim MM, Ihsan Mahmood Helmi Al-Bayati A, Taha A, Abdullaha SA, Jassim Khalil M, Mahdi Rheima A, Hachim SK. Pristine, Co, Rh, and, Ir decorated aluminum phosphide nano-sheet semiconductors as viable chemical sensors for cathinone. INORG CHEM COMMUN 2023. [DOI: 10.1016/j.inoche.2022.110305] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
|
16
|
Kříž K, Schmidt L, Andersson AT, Walz MM, van der Spoel D. An Imbalance in the Force: The Need for Standardized Benchmarks for Molecular Simulation. J Chem Inf Model 2023; 63:412-431. [PMID: 36630710 PMCID: PMC9875315 DOI: 10.1021/acs.jcim.2c01127] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2022] [Indexed: 01/12/2023]
Abstract
Force fields (FFs) for molecular simulation have been under development for more than half a century. As with any predictive model, rigorous testing and comparisons of models critically depends on the availability of standardized data sets and benchmarks. While such benchmarks are rather common in the fields of quantum chemistry, this is not the case for empirical FFs. That is, few benchmarks are reused to evaluate FFs, and development teams rather use their own training and test sets. Here we present an overview of currently available tests and benchmarks for computational chemistry, focusing on organic compounds, including halogens and common ions, as FFs for these are the most common ones. We argue that many of the benchmark data sets from quantum chemistry can in fact be reused for evaluating FFs, but new gas phase data is still needed for compounds containing phosphorus and sulfur in different valence states. In addition, more nonequilibrium interaction energies and forces, as well as molecular properties such as electrostatic potentials around compounds, would be beneficial. For the condensed phases there is a large body of experimental data available, and tools to utilize these data in an automated fashion are under development. If FF developers, as well as researchers in artificial intelligence, would adopt a number of these data sets, it would become easier to compare the relative strengths and weaknesses of different models and to, eventually, restore the balance in the force.
Collapse
Affiliation(s)
- Kristian Kříž
- Department
of Cell and Molecular Biology, Uppsala University, Box 596, SE-75124Uppsala, Sweden
| | - Lisa Schmidt
- Faculty
of Biosciences, University of Heidelberg, Heidelberg69117, Germany
| | - Alfred T. Andersson
- Department
of Cell and Molecular Biology, Uppsala University, Box 596, SE-75124Uppsala, Sweden
| | - Marie-Madeleine Walz
- Department
of Cell and Molecular Biology, Uppsala University, Box 596, SE-75124Uppsala, Sweden
| | - David van der Spoel
- Department
of Cell and Molecular Biology, Uppsala University, Box 596, SE-75124Uppsala, Sweden
| |
Collapse
|
17
|
Negi I, Jangra R, Gharu A, Trant JF, Sharma P. Guanidinium–amino acid hydrogen-bonding interactions in protein crystal structures: implications for guanidinium-induced protein denaturation. Phys Chem Chem Phys 2023; 25:857-869. [DOI: 10.1039/d2cp04943k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Structural analysis of guanidinium–amino acid interaction pairs in protein crystal structures is coupled with an effective scheme for classifying the optimized pairs, to gain understanding of the guanidinium:protein hydrogen bonding modes.
Collapse
Affiliation(s)
- Indu Negi
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh, 160014, India
| | - Raman Jangra
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh, 160014, India
| | - Amit Gharu
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh, 160014, India
| | - John F. Trant
- Department of Chemistry and Biochemistry, University of Windsor, 401 Sunset Ave. Windsor ON, N9B 3P4, Canada
- Binary Star Research Services, LaSalle, ON, N9J 3 X 8, Canada
| | - Purshotam Sharma
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh, 160014, India
- Department of Chemistry and Biochemistry, University of Windsor, 401 Sunset Ave. Windsor ON, N9B 3P4, Canada
| |
Collapse
|
18
|
Theoretical insight into the acidity and cooperativity effect of the LLM-105∙(HNO3)2 system. J Mol Model 2022; 28:401. [DOI: 10.1007/s00894-022-05376-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2022] [Accepted: 11/04/2022] [Indexed: 11/28/2022]
|
19
|
Domagała M, Jabłoński M, Dubis AT, Zabel M, Pfitzner A, Palusiak M. Testing of Exchange-Correlation Functionals of DFT for a Reliable Description of the Electron Density Distribution in Organic Molecules. Int J Mol Sci 2022; 23:ijms232314719. [PMID: 36499046 PMCID: PMC9740346 DOI: 10.3390/ijms232314719] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2022] [Revised: 11/08/2022] [Accepted: 11/18/2022] [Indexed: 11/27/2022] Open
Abstract
Researchers carrying out calculations using the DFT method face the problem of the correct choice of the exchange-correlation functional to describe the quantities they are interested in. This article deals with benchmark calculations aimed at testing various exchange-correlation functionals in terms of a reliable description of the electron density distribution in molecules. For this purpose, 30 functionals representing all rungs of Jacob's Ladder are selected and then the values of some QTAIM-based parameters are compared with their reference equivalents obtained at the CCSD/aug-cc-pVTZ level of theory. The presented results show that the DFT method undoubtedly has the greatest problems with a reliable description of the electron density distribution in multiple strongly polar bonds, such as C=O, and bonds associated with large electron charge delocalization. The performance of the tested functionals turned out to be unsystematic. Nevertheless, in terms of a reliable general description of QTAIM-based parameters, the M11, SVWN, BHHLYP, M06-HF, and, to a slightly lesser extent, also BLYP, B3LYP, and X3LYP functionals turned out to be the worst. It is alarming to find the most popular B3LYP functional in this group. On the other hand, in the case of the electron density at the bond critical point, being the most important QTAIM-based parameter, the M06-HF functional is especially discouraged due to the very poor description of the C=O bond. On the contrary, the VSXC, M06-L, SOGGA11-X, M06-2X, MN12-SX, and, to a slightly lesser extent, also TPSS, TPSSh, and B1B95 perform well in this respect. Particularly noteworthy is the overwhelming performance of double hybrids in terms of reliable values of bond delocalization indices. The results show that there is no clear improvement in the reliability of describing the electron density distribution with climbing Jacob's Ladder, as top-ranked double hybrids are also, in some cases, able to produce poor values compared to CCSD.
Collapse
Affiliation(s)
- Małgorzata Domagała
- Faculty of Chemistry, University of Lodz, Pomorska 163/165, 90-236 Lodz, Poland
| | - Mirosław Jabłoński
- Faculty of Chemistry, Nicolaus Copernicus University in Toruń, Gagarina 7, 87-100 Toruń, Poland
- Correspondence: ; Tel.: +48-056-611-4695
| | - Alina T. Dubis
- Faculty of Chemistry, University of Białystok, Ciołkowskiego 1K, 15-245 Białystok, Poland
| | - Manfred Zabel
- Institute of Inorganic Chemistry, University of Regensburg, Universitätsstr. 31, 93040 Regensburg, Germany
| | - Arno Pfitzner
- Institute of Inorganic Chemistry, University of Regensburg, Universitätsstr. 31, 93040 Regensburg, Germany
| | - Marcin Palusiak
- Faculty of Chemistry, University of Lodz, Pomorska 163/165, 90-236 Lodz, Poland
| |
Collapse
|
20
|
Su J, Wang L, Lu Q, Li M, Guan Y, Wang C, Fu S. Alkali resistance mechanism of cyano-containing heterocyclic disperse dyes. J Mol Struct 2022. [DOI: 10.1016/j.molstruc.2022.133438] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
21
|
A DFT study on the tautomerization of vitamin B3 (niacin). COMPUT THEOR CHEM 2022. [DOI: 10.1016/j.comptc.2022.113800] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
22
|
Biswas A, Mallik BS. Ionic Dynamics and Vibrational Spectral Diffusion of a Protic Alkylammonium Ionic Salt through Intrinsic Cationic N-H Vibrational Probe from FPMD Simulations. J Phys Chem A 2022; 126:5134-5147. [PMID: 35900106 DOI: 10.1021/acs.jpca.2c03387] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We employed density functional theory (DFT)-based molecular dynamics simulations to explore the structure, dynamics, and spectral properties of the protic ionic entity trimethylammonium chloride (TMACl). Structural investigations include calculating the site-site radial distribution functions (RDFs), the distribution of constituent cations and anions in three-dimensional space, and combined distribution functions of the hydrogen-bonded pair RDF versus angle, revealing the structural characteristics of the ionic solvation and the intermolecular interactions within ions. Further, we determined the instantaneous vibrational stretching frequencies of the intrinsic N-H stretch probe modes by applying the time-series wavelet method. The associated ionic dynamics within the protic ionic compound were investigated by calculating the time-evolution of the fluctuating frequencies and the frequency-time correlation functions (FFCFs). The time scale related to the local structural relaxation process and the average hydrogen bond lifetime, ion cage dynamics, and mean squared displacement were investigated. The faster decay component of the FFCFs, depicting the intermolecular motion of intact hydrogen bonds in TMACl, is 0.07 ps for the Perdew-Burke-Ernzerhof (PBE)-based simulation and 0.06 ps for the PBE-D2 representation. The slower time scale of the longer picosecond decay time component of PBE and PBE-D2 representations are 3.13 and 2.87 ps, respectively. These picosecond time scales represent more significant fluctuations of the hydrogen-bonding partners in the ionic entity and hydrogen-bond jump events accompanied by large angular jumps. The longest picosecond time scales represent structural relaxation, including large angular jumps and ion-pair dynamics. Also, ion cage lifetimes correlate with the slowest time scale of the associated dynamics of vibrational spectral diffusion despite the type of DFT functional. This study benchmarks DFT treatments of the exchange-correlation functional with and without the van der Waals (vdW) dispersion correction scheme. The inclusion of vdW interactions to the PBE functional represents a less structured state of the ionic entity and faster dynamics of the molecular motions relative to the one predicted by the PBE system. All the results illustrate the necessity of accurately describing the Coulomb interactions, vdW dispersive interactive forces, and localized hydrogen bonds required to sustain the energetic balance in this ionic salt.
Collapse
Affiliation(s)
- Aritri Biswas
- Department of Chemistry, Indian Institute of Technology Hyderabad, Sangareddy 502284, Telangana, India
| | - Bhabani S Mallik
- Department of Chemistry, Indian Institute of Technology Hyderabad, Sangareddy 502284, Telangana, India
| |
Collapse
|
23
|
Yamabe S, Tsuchida N, Yamazaki S. A DFT study on the degradation mechanism of vitamin B2. FOOD CHEMISTRY: MOLECULAR SCIENCES 2022; 4:100080. [PMID: 35415686 PMCID: PMC8991990 DOI: 10.1016/j.fochms.2022.100080] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/27/2021] [Revised: 01/18/2022] [Accepted: 01/22/2022] [Indexed: 11/26/2022]
Abstract
Elementary processes of degradation from riboflavin were found by DFT calculations. Photochemical reaction courses in the lowest triplet spin state were elucidated. Base-catalyzed degradation paths from formylmethylflavin were determined. All the transition states of cleavage of C-C and C-N covalent bonds were determined.
Degradation reaction paths starting from riboflavin (RF) were investigated using DFT (density functional theory) as the first attempt to reveal their elementary processes. Photochemical reactions were followed in the lowest triplet spin state, “(T)”. Two intermediates [Int1(T) and Int2(T)] were found in the course, RF(T) → FMF (7,8-dimethyl-10-formylflavin, T). From FMF(T), there are two degradation channels. Release of ketene(T) and carbon monoxide leads to LC (lumichrome, S0) and LF (lumiflavin, T), respectively. The base-catalyzed (ground state) degradation of FMF was investigated with HO–(H2O)3. The Grotthuss-type proton transfer along hydrogen bonds controlled the degradation reaction. All the transition states of cleavage of C—C and C—N covalent bonds were determined, and the degradation mechanism was clarified.
Collapse
|
24
|
Yamabe S, Tsuchida N, Yamazaki S. A DFT study of the active role of the phosphate group of an internal aldimine in a transamination reaction. Org Biomol Chem 2022; 20:5334-5341. [PMID: 35748359 DOI: 10.1039/d2ob00913g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A transamination reaction from an internal aldimine ([PLP]) and (S)-alanine to pyridoxamine phosphate (PMP) and pyruvic acid was investigated by DFT calculations. As [PLP], a model where the lysine (-Lys) part was approximated by -CH[-NH-C(O)-CH3]-C(O)-NH-CH3 was adopted. (H2O)4 was also included to trace reaction paths involving proton transfers. 13 elementary processes were obtained. For (the external aldimine → quinoid), (quinoid → ketimine) and (ketimine → carbinol amine) processes, the water dimer was found to connect a phosphate-group oxygen with the moving proton. The connection promoted the Grotthuss-type proton transfer in transition states. It was revealed that the phosphate group is not a mere substituent but has the central role in the transfer.
Collapse
Affiliation(s)
- Shinichi Yamabe
- Department of Chemistry, Nara University of Education, Takabatake-cho, Nara 630-8528, Japan.
| | - Noriko Tsuchida
- Department of Liberal Arts, Faculty of Medicine, Saitama Medical University, 38 Morohongo, Moroyama-machi, Iruma-gun, Saitama 350-0495, Japan
| | - Shoko Yamazaki
- Department of Chemistry, Nara University of Education, Takabatake-cho, Nara 630-8528, Japan.
| |
Collapse
|
25
|
Prasad VK, Otero-de-la-Roza A, DiLabio GA. Small-Basis Set Density-Functional Theory Methods Corrected with Atom-Centered Potentials. J Chem Theory Comput 2022; 18:2913-2930. [PMID: 35412817 DOI: 10.1021/acs.jctc.2c00036] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Density functional theory (DFT) is currently the most popular method for modeling noncovalent interactions and thermochemistry. The accurate calculation of noncovalent interaction energies, reaction energies, and barrier heights requires choosing an appropriate functional and, typically, a relatively large basis set. Deficiencies of the density-functional approximation and the use of a limited basis set are the leading sources of error in the calculation of noncovalent and thermochemical properties in molecular systems. In this article, we present three new DFT methods based on the BLYP, M06-2X, and CAM-B3LYP functionals in combination with the 6-31G* basis set and corrected with atom-centered potentials (ACPs). ACPs are one-electron potentials that have the same form as effective-core potentials, except they do not replace any electrons. The ACPs developed in this work are used to generate energy corrections to the underlying DFT/basis-set method such that the errors in predicted chemical properties are minimized while maintaining the low computational cost of the parent methods. ACPs were developed for the elements H, B, C, N, O, F, Si, P, S, and Cl. The ACP parameters were determined using an extensive training set of 118655 data points, mostly of complete basis set coupled-cluster level quality. The target molecular properties for the ACP-corrected methods include noncovalent interaction energies, molecular conformational energies, reaction energies, barrier heights, and bond separation energies. The ACPs were tested first on the training set and then on a validation set of 42567 additional data points. We show that the ACP-corrected methods can predict the target molecular properties with accuracy close to complete basis set wavefunction theory methods, but at a computational cost of double-ζ DFT methods. This makes the new BLYP/6-31G*-ACP, M06-2X/6-31G*-ACP, and CAM-B3LYP/6-31G*-ACP methods uniquely suited to the calculation of noncovalent, thermochemical, and kinetic properties in large molecular systems.
Collapse
Affiliation(s)
- Viki Kumar Prasad
- Department of Chemistry, University of British Columbia, Okanagan, 3247 University Way, Kelowna, British Columbia V1V 1V7, Canada
| | - Alberto Otero-de-la-Roza
- Departamento de Química Física y Analítica, Facultad de Química, Universidad de Oviedo, MALTA Consolider Team, Oviedo E-33006, Spain
| | - Gino A DiLabio
- Department of Chemistry, University of British Columbia, Okanagan, 3247 University Way, Kelowna, British Columbia V1V 1V7, Canada
| |
Collapse
|
26
|
Interaction of simple amino acids (glycine, α-alanine, β-alanine and L-valine) with germatranol hydrate. J Mol Struct 2022. [DOI: 10.1016/j.molstruc.2021.132245] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
27
|
Prasad VK, Otero-de-la-Roza A, DiLabio GA. Fast and Accurate Quantum Mechanical Modeling of Large Molecular Systems Using Small Basis Set Hartree-Fock Methods Corrected with Atom-Centered Potentials. J Chem Theory Comput 2022; 18:2208-2232. [PMID: 35313106 DOI: 10.1021/acs.jctc.1c01128] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
There has been significant interest in developing fast and accurate quantum mechanical methods for modeling large molecular systems. In this work, by utilizing a machine learning regression technique, we have developed new low-cost quantum mechanical approaches to model large molecular systems. The developed approaches rely on using one-electron Gaussian-type functions called atom-centered potentials (ACPs) to correct for the basis set incompleteness and the lack of correlation effects in the underlying minimal or small basis set Hartree-Fock (HF) methods. In particular, ACPs are proposed for ten elements common in organic and bioorganic chemistry (H, B, C, N, O, F, Si, P, S, and Cl) and four different base methods: two minimal basis sets (MINIs and MINIX) plus a double-ζ basis set (6-31G*) in combination with dispersion-corrected HF (HF-D3/MINIs, HF-D3/MINIX, HF-D3/6-31G*) and the HF-3c method. The new ACPs are trained on a very large set (73 832 data points) of noncovalent properties (interaction and conformational energies) and validated additionally on a set of 32 048 data points. All reference data are of complete basis set coupled-cluster quality, mostly CCSD(T)/CBS. The proposed ACP-corrected methods are shown to give errors in the tenths of a kcal/mol range for noncovalent interaction energies and up to 2 kcal/mol for molecular conformational energies. More importantly, the average errors are similar in the training and validation sets, confirming the robustness and applicability of these methods outside the boundaries of the training set. In addition, the performance of the new ACP-corrected methods is similar to complete basis set density functional theory (DFT) but at a cost that is orders of magnitude lower, and the proposed ACPs can be used in any computational chemistry program that supports effective-core potentials without modification. It is also shown that ACPs improve the description of covalent and noncovalent bond geometries of the underlying methods and that the improvement brought about by the application of the ACPs is directly related to the number of atoms to which they are applied, allowing the treatment of systems containing some atoms for which ACPs are not available. Overall, the ACP-corrected methods proposed in this work constitute an alternative accurate, economical, and reliable quantum mechanical approach to describe the geometries, interaction energies, and conformational energies of systems with hundreds to thousands of atoms.
Collapse
Affiliation(s)
- Viki Kumar Prasad
- Department of Chemistry, University of British Columbia, Okanagan, 3247 University Way, Kelowna, British Columbia, Canada V1V 1V7
| | - Alberto Otero-de-la-Roza
- MALTA Consolider Team, Departamento de Química Física y Analítica, Facultad de Química, Universidad de Oviedo, E-33006 Oviedo, Spain
| | - Gino A DiLabio
- Department of Chemistry, University of British Columbia, Okanagan, 3247 University Way, Kelowna, British Columbia, Canada V1V 1V7
| |
Collapse
|
28
|
Priyadarsini A, Mallik BS. Structure and rotational dynamics of water around hydrogen peroxide. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2021.118054] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
29
|
Mohamed SY, Monge-Palacios M, Giri BR, Khaled F, Liu D, Farooq A, Sarathy SM. The Effect of Hydrogen Bonding on the Reactivity of OH Radicals with Prenol and Isoprenol: A Shock Tube and Multi-Structural Torsional Variational Transition State Theory Study. Phys Chem Chem Phys 2022; 24:12601-12620. [DOI: 10.1039/d2cp00737a] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The presence of two functional groups (OH and double bond) in C5 methyl-substituted enols (i.e., isopentenols), such as 3-methyl-2-buten-1-ol (prenol) and 3-methyl-3-buten-1-ol (isoprenol), makes them excellent biofuel candidates as fuel...
Collapse
|
30
|
Fan D, Chen L, Wang C, Yin S, Mo Y. Inter-anion chalcogen bonds: Are they anti-electrostatic in nature? J Chem Phys 2021; 155:234302. [PMID: 34937369 DOI: 10.1063/5.0076872] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
Inter-anion hydrogen and halogen bonds have emerged as counterintuitive linkers and inspired us to expand the range of this unconventional bonding pattern. Here, the inter-anion chalcogen bond (IAChB) was proposed and theoretically analyzed in a series of complexes formed by negatively charged bidentate chalcogen bond donors with chloride anions. The kinetic stability of IAChB was evidenced by the minima on binding energy profiles and further supported by ab initio molecular dynamic simulations. The block-localized wave function (BLW) method and its subsequent energy decomposition (BLW-ED) approach were employed to elucidate the physical origin of IAChB. While all other energy components vary monotonically as anions get together, the electrostatic interaction behaves exceptionally as it experiences a Coulombic repulsion barrier. Before reaching the barrier, the electrostatic repulsion increases with the shortening Ch⋯Cl- distance as expected from classical electrostatics. However, after passing the barrier, the electrostatic repulsion decreases with the Ch⋯Cl- distance shortening and subsequently turns into the most favorable trend among all energy terms at short ranges, representing a dominating force for the kinetic stability of inter-anions. For comparison, all energy components exhibit the same trends and vary monotonically in the conventional counterparts where donors are neutral. By comparing inter-anions and their conventional counterparts, we found that only the electrostatic energy term is affected by the extra negative charge. Remarkably, the distinctive (nonmonotonic) electrostatic energy profiles were reproduced using quantum mechanical-based atomic multipoles, suggesting that the crucial electrostatic interaction in IAChB can be rationalized within the classical electrostatic theory just like conventional non-covalent interactions.
Collapse
Affiliation(s)
- Dan Fan
- Key Laboratory for Macromolecular Science of Shaanxi Province, School of Chemistry & Chemical Engineering, Shaanxi Normal University, Xi'an 710119, China
| | - Li Chen
- Key Laboratory for Macromolecular Science of Shaanxi Province, School of Chemistry & Chemical Engineering, Shaanxi Normal University, Xi'an 710119, China
| | - Changwei Wang
- Key Laboratory for Macromolecular Science of Shaanxi Province, School of Chemistry & Chemical Engineering, Shaanxi Normal University, Xi'an 710119, China
| | - Shiwei Yin
- Key Laboratory for Macromolecular Science of Shaanxi Province, School of Chemistry & Chemical Engineering, Shaanxi Normal University, Xi'an 710119, China
| | - Yirong Mo
- Department of Nanoscience, Joint School of Nanoscience and Nanoengineering, University of North Carolina at Greensboro, Greensboro, North Carolina 27401, USA
| |
Collapse
|
31
|
Sparrow ZM, Ernst BG, Joo PT, Lao KU, DiStasio RA. NENCI-2021. I. A large benchmark database of non-equilibrium non-covalent interactions emphasizing close intermolecular contacts. J Chem Phys 2021; 155:184303. [PMID: 34773949 DOI: 10.1063/5.0068862] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this work, we present NENCI-2021, a benchmark database of ∼8000 Non-Equilibirum Non-Covalent Interaction energies for a large and diverse selection of intermolecular complexes of biological and chemical relevance. To meet the growing demand for large and high-quality quantum mechanical data in the chemical sciences, NENCI-2021 starts with the 101 molecular dimers in the widely used S66 and S101 databases and extends the scope of these works by (i) including 40 cation-π and anion-π complexes, a fundamentally important class of non-covalent interactions that are found throughout nature and pose a substantial challenge to theory, and (ii) systematically sampling all 141 intermolecular potential energy surfaces (PESs) by simultaneously varying the intermolecular distance and intermolecular angle in each dimer. Designed with an emphasis on close contacts, the complexes in NENCI-2021 were generated by sampling seven intermolecular distances along each PES (ranging from 0.7× to 1.1× the equilibrium separation) and nine intermolecular angles per distance (five for each ion-π complex), yielding an extensive database of 7763 benchmark intermolecular interaction energies (Eint) obtained at the coupled-cluster with singles, doubles, and perturbative triples/complete basis set [CCSD(T)/CBS] level of theory. The Eint values in NENCI-2021 span a total of 225.3 kcal/mol, ranging from -38.5 to +186.8 kcal/mol, with a mean (median) Eint value of -1.06 kcal/mol (-2.39 kcal/mol). In addition, a wide range of intermolecular atom-pair distances are also present in NENCI-2021, where close intermolecular contacts involving atoms that are located within the so-called van der Waals envelope are prevalent-these interactions, in particular, pose an enormous challenge for molecular modeling and are observed in many important chemical and biological systems. A detailed symmetry-adapted perturbation theory (SAPT)-based energy decomposition analysis also confirms the diverse and comprehensive nature of the intermolecular binding motifs present in NENCI-2021, which now includes a significant number of primarily induction-bound dimers (e.g., cation-π complexes). NENCI-2021 thus spans all regions of the SAPT ternary diagram, thereby warranting a new four-category classification scheme that includes complexes primarily bound by electrostatics (3499), induction (700), dispersion (1372), or mixtures thereof (2192). A critical error analysis performed on a representative set of intermolecular complexes in NENCI-2021 demonstrates that the Eint values provided herein have an average error of ±0.1 kcal/mol, even for complexes with strongly repulsive Eint values, and maximum errors of ±0.2-0.3 kcal/mol (i.e., ∼±1.0 kJ/mol) for the most challenging cases. For these reasons, we expect that NENCI-2021 will play an important role in the testing, training, and development of next-generation classical and polarizable force fields, density functional theory approximations, wavefunction theory methods, and machine learning based intra- and inter-molecular potentials.
Collapse
Affiliation(s)
- Zachary M Sparrow
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Brian G Ernst
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Paul T Joo
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Ka Un Lao
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Robert A DiStasio
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| |
Collapse
|
32
|
Jabłoński M, Krygowski TM. Dependence of the substituent energy on the level of theory. J Comput Chem 2021; 42:2079-2088. [PMID: 34432304 DOI: 10.1002/jcc.26738] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2021] [Revised: 06/24/2021] [Accepted: 07/22/2021] [Indexed: 01/08/2023]
Abstract
Most often, the substituent effects are described using rather troublesome Hammett constants. Quite recently, it has been proposed to use the so-called substituent energy, which is based on total energies of the X-substituted polycyclic aromatic hydrocarbon and phenyl. This article concerns the influence of the applied level of theory (i.e., both the basis set and the method) on the determined values of the substituent energies. For this purpose, the energies of the NH2 and NO2 groups in 16 unique positions of naphthalene, anthracene, tetracene, phenanthrene, and pyrene were calculated using 10 different basis sets and 23 various exchange-correlation functionals representing the entire Jacob's Ladder, from local, through gradient- and meta-gradient-corrected, to hybrid and double-hybrid ones. Additionally, using even larger 6-311++G(2df,2p) basis set, the energies of NH2 , NO2 , CN, and Cl were also computed. Both the basis set and the method used have little effect on the substituent energy if the substituent is in the benzene-like position. On the contrary, the effect of the level of theory is pronounced especially in the case of the most spatially crowded 4-substituted phenanthrene. Substituent energies have been shown to be very useful theoretical parameters describing the proximity effect in the substituted derivatives of polycyclic aromatic hydrocarbons.
Collapse
Affiliation(s)
- Mirosław Jabłoński
- Faculty of Chemistry, Nicolaus Copernicus University in Toruń, Toruń, Poland
| | | |
Collapse
|
33
|
Yamabe S, Tsuchida N, Yamazaki S. How Is the Oxidation Related to the Tautomerization in Vitamin B9? J Phys Chem A 2021; 125:9346-9354. [PMID: 34663066 DOI: 10.1021/acs.jpca.1c06678] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The relationship between the lactim-lactam tautomerization and the free-radical scavenging reaction in vitamin B9 [folic acid (FA)] was investigated by density functional theory calculations. 6-Methylpterin was also adopted for the detailed analyses of various reaction paths. For pterin, the transition state of the tautomerization with two water molecules (n = 2) was calculated to be of the lowest activation energy. The proton-transfer circuit of n = 2 is retained (not broken) even with the addition of outer water molecules, n = 2 + 2, 2 + 4, 2 + 8, and 2 + 14. At the oxidation of the system composed of 6-methylpterin + (H2O)2 + HO•, the radical character of HO• is directly transmitted to the pterin ring along with the C-O → H → O → H → O → H → OH proton transfer. The patterns of the electron transfer (pterin ring → OX•) and the concomitant proton transfer via the water dimer were commonly obtained for the oxidant (OX•) = HO•, Cl3C-O2•, N3•, or SO4-•. The hydrogen atom transfer mechanism was ruled out. Two conformations of the puckered form with the -C(═O)-OH···N intramolecular hydrogen bonds of FA were found to have the stability similar to that of the linear conformer. Both the tautomerization and the oxidation were calculated to occur competitively in the three conformers.
Collapse
Affiliation(s)
- Shinichi Yamabe
- Department of Chemistry, Nara University of Education, Takabatake-cho, Nara 630-8528, Japan
| | - Noriko Tsuchida
- Department of Liberal Arts, Faculty of Medicine, Saitama Medical University, 38 Morohongo, Moroyama-machi, Iruma-gun, Saitama 350-0495, Japan
| | - Shoko Yamazaki
- Department of Chemistry, Nara University of Education, Takabatake-cho, Nara 630-8528, Japan
| |
Collapse
|
34
|
Zare S, Qomi MJA. Reactive force fields for aqueous and interfacial magnesium carbonate formation. Phys Chem Chem Phys 2021; 23:23106-23123. [PMID: 34617078 DOI: 10.1039/d1cp02627e] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
We develop Mg/C/O/H ReaxFF parameter sets for two environments: an aqueous force field for magnesium ions in solution and an interfacial force field for minerals and mineral-water interfaces. Since magnesium is highly ionic, we choose to fix the magnesium charge and model its interaction with C/O/H through Coulomb, Lennard-Jones, and Buckingham potentials. We parameterize the forcefields against several crystal structures, including brucite, magnesite, magnesia, magnesium hydride, and magnesium carbide, as well as Mg2+ water binding energies for the aqueous forcefield. Then, we test the forcefield for other magnesium-containing crystals, solvent separated and contact ion-pairs and single-molecule/multilayer water adsorption energies on mineral surfaces. We also apply the forcefield to the forsterite-water and brucite-water interface that contains a bicarbonate ion. We observe that a long-range proton transfer mechanism deprotonates the bicarbonate ion to carbonate at the interface. Free energy calculations show that carbonate can attach to the magnesium surface with an energy barrier of about 0.22 eV, consistent with the free energy required for aqueous Mg-CO3 ion pairing. Also, the diffusion constant of the hydroxide ions in the water layers formed on the forsterite surface are shown to be anisotropic and heterogeneous. These findings can help explain the experimentally observed fast nucleation and growth of magnesite at low temperature at the mineral-water-CO2 interface in water-poor conditions.
Collapse
Affiliation(s)
- Siavash Zare
- Department of Civil and Environmental Engineering, University of California, Irvine, CA, USA.
| | | |
Collapse
|
35
|
Le Bras L, Dory YL, Champagne B. Computational prediction of the supramolecular self-assembling properties of organic molecules: the role of conformational flexibility of amide moieties. Phys Chem Chem Phys 2021; 23:20453-20465. [PMID: 34498627 DOI: 10.1039/d1cp02675e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Two families of organic molecules with different backbones have been considered. The first family is based on a macrolactam-like unit that is constrained in a particular conformation. The second family is composed by a substituted central phenyl that allows a larger mobility for its substituents. They have however a common feature, three amide moieties (within the cycle for the macrolactam-like molecule and as substituents for the phenyl) that permit hydrogen bonding when molecules are stacked. In this study we propose a computational protocol to unravel the ability of the different families to self-assemble into organic nanotubes. Starting from the monomer and going towards larger assemblies like dimers, trimers, and pentamers we applied the different protocols to rationalize the behavior of the different assemblies. Both structures and thermodynamics were investigated to give a complete picture of the process. Thanks to the combination of a quantum mechanics approach and molecular dynamics simulations along with the use of tailored tools (non covalent interaction visualization) and techniques (umbrella sampling), we have been able to differentiate the two families and highlight the best candidate for self-assembling purposes.
Collapse
Affiliation(s)
- Laura Le Bras
- Unité de Chimie Physique Théorique et Structurale, Chemistry Department, Namur Institute of Structured Matter, University of Namur, Belgium.
| | - Yves L Dory
- Laboratoire de Synthèse Supramoléculaire, Département de Chimie, Institut de Pharmacologie, Université de Sherbrooke, Sherbrooke, Québec, Canada
| | - Benoît Champagne
- Unité de Chimie Physique Théorique et Structurale, Chemistry Department, Namur Institute of Structured Matter, University of Namur, Belgium.
| |
Collapse
|
36
|
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] [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
|
37
|
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] [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
|
38
|
Malloum A, Conradie J. Accurate binding energies of ammonia clusters and benchmarking of hybrid DFT functionals. COMPUT THEOR CHEM 2021. [DOI: 10.1016/j.comptc.2021.113236] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
39
|
Yamabe S, Tsuchida N, Yamazaki S. A density functional theory study of the hydride shift in the Eschweiler–Clarke reaction. J PHYS ORG CHEM 2021. [DOI: 10.1002/poc.4253] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Affiliation(s)
- Shinichi Yamabe
- Department of Chemistry Nara University of Education Nara Japan
| | - Noriko Tsuchida
- Department of Liberal Arts, Faculty of Medicine Saitama Medical University Saitama Japan
| | - Shoko Yamazaki
- Department of Chemistry Nara University of Education Nara Japan
| |
Collapse
|
40
|
Yamabe S, Tsuchida N, Yamazaki S. How is vitamin B1 oxidized to thiochrome? Elementary processes revealed by a DFT study. Org Biomol Chem 2021; 19:4529-4536. [PMID: 33929469 DOI: 10.1039/d1ob00677k] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
The oxidation reaction of thiamine (vitamin B1) to thiochrome was investigated by DFT calculations. Three reaction systems, [A] thiamine + methyl peroxy radical + (H2O)8, [B] thiamine + cyanogen bromide + HO-(H2O)8 and [C] thiamine + mercury(ii) chloride + HO-(H2O)8, were investigated. wB97X-D/6-311+G** for [A] and [B] and wB97X-D/SDD&6-311(+)G** for [C] geometry optimizations were carried out with the solvent effect (water). The effect is of the self-consistent reaction field (SCRF) with the polarizable continuum model (PCM). In [A], the H3C-O2˙ adduct of thiamine undergoes simultaneous cleavage of the C-H and O-O bonds, leading to a very stable 2(3H)-thiazolone intermediate. The same intermediate was obtained after the cleavage of the C-H and O-H bonds of the HO adduct of thiamine in [B] and [C]. After the formation of the key intermediate, the N-protonated thiochrome was afforded via three steps. In reflection of the water-soluble character of vitamin B1, proton transfers along hydrogen bonds of the water cluster enhance those steps.
Collapse
Affiliation(s)
- Shinichi Yamabe
- Department of Chemistry, Nara University of Education, Takabatake-cho, Nara 630-8528, Japan.
| | - Noriko Tsuchida
- Department of Liberal Arts, Faculty of Medicine, Saitama Medical University, 38 Morohongo, Moroyama-machi, Iruma-gun, Saitama 350-0495, Japan
| | - Shoko Yamazaki
- Department of Chemistry, Nara University of Education, Takabatake-cho, Nara 630-8528, Japan.
| |
Collapse
|
41
|
Cooperation/Competition between Halogen Bonds and Hydrogen Bonds in Complexes of 2,6-Diaminopyridines and X-CY3 (X = Cl, Br; Y = H, F). Symmetry (Basel) 2021. [DOI: 10.3390/sym13050766] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Abstract
The DFT calculations have been performed on a series of two-element complexes formed by substituted 2,6-diaminopyridine (R−PDA) and pyridine (R−Pyr) with X−CY3 molecules (where X = Cl, Br and Y = H, F). The primary aim of this study was to examine the intermolecular hydrogen and halogen bonds in the condition of their mutual coexistence. Symmetry/antisymmetry of the interrelation between three individual interactions is addressed. It appears that halogen bonds play the main role in the stabilization of the structures of the selected systems. However, the occurrence of one or two hydrogen bonds was associated with the favourable geometry of the complexes. Moreover, the impact of different substituent groups attached in the para position to the aromatic ring of the 2,6-diaminopyridine and pyridine on the character of the intermolecular hydrogen and halogen bonds was examined. The results indicate that the presence of electron-donating substituents strengthens the bonds. In turn, the presence of electron-withdrawing substituents reduces the strength of halogen bonds. Additionally, when hydrogen and halogen bonds lose their leading role in the complex formation, the nonspecific electrostatic interactions between dipole moments take their place. Analysis was based on geometric, energetic, and topological parameters of the studied systems.
Collapse
|
42
|
Assessing the Intrinsic Strengths of Ion–Solvent and Solvent–Solvent Interactions for Hydrated Mg2+ Clusters. INORGANICS 2021. [DOI: 10.3390/inorganics9050031] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022] Open
Abstract
Information resulting from a comprehensive investigation into the intrinsic strengths of hydrated divalent magnesium clusters is useful for elucidating the role of aqueous solvents on the Mg2+ ion, which can be related to those in bulk aqueous solution. However, the intrinsic Mg–O and intermolecular hydrogen bond interactions of hydrated magnesium ion clusters have yet to be quantitatively measured. In this work, we investigated a set of 17 hydrated divalent magnesium clusters by means of local vibrational mode force constants calculated at the ωB97X-D/6-311++G(d,p) level of theory, where the nature of the ion–solvent and solvent–solvent interactions were interpreted from topological electron density analysis and natural population analysis. We found the intrinsic strength of inner shell Mg–O interactions for [Mg(H2O)n]2+ (n = 1–6) clusters to relate to the electron density at the bond critical point in Mg–O bonds. From the application of a secondary hydration shell to [Mg(H2O)n]2+ (n = 5–6) clusters, stronger Mg–O interactions were observed to correspond to larger instances of charge transfer between the lp(O) orbitals of the inner hydration shell and the unfilled valence shell of Mg. As the charge transfer between water molecules of the first and second solvent shell increased, so did the strength of their intermolecular hydrogen bonds (HBs). Cumulative local vibrational mode force constants of explicitly solvated Mg2+, having an outer hydration shell, reveal a CN of 5, rather than a CN of 6, to yield slightly more stable configurations in some instances. However, the cumulative local mode stretching force constants of implicitly solvated Mg2+ show the six-coordinated cluster to be the most stable. These results show that such intrinsic bond strength measures for Mg–O and HBs offer an effective way for determining the coordination number of hydrated magnesium ion clusters.
Collapse
|
43
|
Kodrycka M, Patkowski K. Efficient Density-Fitted Explicitly Correlated Dispersion and Exchange Dispersion Energies. J Chem Theory Comput 2021; 17:1435-1456. [PMID: 33606539 DOI: 10.1021/acs.jctc.0c01158] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The leading-order dispersion and exchange-dispersion terms in symmetry-adapted perturbation theory (SAPT), Edisp(20) and Eexch-disp(20), suffer from slow convergence to the complete basis set limit. To alleviate this problem, explicitly correlated variants of these corrections, Edisp(20)-F12 and Eexch-disp(20)-F12, have been proposed recently. However, the original formalism (M., Kodrycka , J. Chem. Theory Comput. 2019, 15, 5965-5986), while highly successful in terms of improving convergence, was not competitive to conventional orbital-based SAPT in terms of computational efficiency due to the need to manipulate several kinds of two-electron integrals. In this work, we eliminate this need by decomposing all types of two-electron integrals using robust density fitting. We demonstrate that the error of the density fitting approximation is negligible when standard auxiliary bases such as aug-cc-pVXZ/MP2FIT are employed. The new implementation allowed us to study all complexes in the A24 database in basis sets up to aug-cc-pV5Z, and the Edisp(20)-F12 and Eexch-disp(20)-F12 values exhibit vastly improved basis set convergence over their conventional counterparts. The well-converged Edisp(20)-F12 and Eexch-disp(20)-F12 numbers can be substituted for conventional Edisp(20) and Eexch-disp(20) ones in a calculation of the total SAPT interaction energy at any level (SAPT0, SAPT2+3, ...). We show that the addition of F12 terms does not improve the accuracy of low-level SAPT treatments. However, when the theory errors are minimized in high-level SAPT approaches such as SAPT2+3(CCD)δMP2, the reduction of basis set incompleteness errors thanks to the F12 treatment substantially improves the accuracy of small-basis calculations.
Collapse
Affiliation(s)
- Monika Kodrycka
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849, United States
| | - Konrad Patkowski
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849, United States
| |
Collapse
|
44
|
Firme CL. Local potential energy density model (LPE): Applications and limitations to quantify intra/intermolecular interactions. COMPUT THEOR CHEM 2021. [DOI: 10.1016/j.comptc.2021.113143] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
45
|
Yamabe S, Tsuchida N, Yamazaki S. A density functional theory study of the reaction mechanism of formation of phenolphthalein and fluorescein. J PHYS ORG CHEM 2021. [DOI: 10.1002/poc.4136] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Affiliation(s)
- Shinichi Yamabe
- Department of Chemistry Nara University of Education Nara Japan
| | - Noriko Tsuchida
- Department of Liberal Arts, Faculty of Medicine Saitama Medical University Saitama Japan
| | - Shoko Yamazaki
- Department of Chemistry Nara University of Education Nara Japan
| |
Collapse
|
46
|
Waldrop JM, Patkowski K. Nonapproximated third-order exchange induction energy in symmetry-adapted perturbation theory. J Chem Phys 2021; 154:024103. [PMID: 33445897 DOI: 10.1063/5.0035050] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The exchange terms in symmetry-adapted perturbation theory (SAPT) are normally calculated within the so-called S2 or single exchange approximation, which approximates the all-electron antisymmetrizer by interchanges of at most one electron pair between the interacting molecules. This approximation is typically very accurate at the van der Waals minimum separation and at larger intermolecular distances but begins to deteriorate at short range. Nonapproximated expressions for the second-order SAPT exchange corrections have been derived some time ago by Schäffer and Jansen [Mol. Phys. 111, 2570 (2013)]. In this work, we extend Schäffer and Jansen's formalism to derive and implement a nonapproximated expression for the third-order exchange-induction correction. Numerical tests on several representative noncovalent databases show that the S2 approximation underestimates the exchange-induction contributions in both second and third orders. This underestimation is very similar in relative terms, but the larger absolute values of the third-order exchange-induction effects, and their near complete cancellation with the corresponding induction energies, make the third-order errors more severe. In the worst-case scenario of interactions involving ions, the breakdown of the S2 approximation can result in a qualitatively wrong, attractive character of SAPT total energies at short range {as first observed by Lao and Herbert [J. Phys. Chem. A 116, 3042 (2012)]}. As expected, the inclusion of the full third-order exchange-induction energy in place of its S2-approximated counterpart restores the correct, repulsive short-range behavior of the SAPT potential energy curves computed through the third order.
Collapse
Affiliation(s)
- Jonathan M Waldrop
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849, USA
| | - Konrad Patkowski
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849, USA
| |
Collapse
|
47
|
Spectral Probe for Electron Transfer and Addition Reactions of Azide Radicals with Substituted Quinoxalin-2-Ones in Aqueous Solutions. Int J Mol Sci 2021; 22:ijms22020633. [PMID: 33435233 PMCID: PMC7828026 DOI: 10.3390/ijms22020633] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2020] [Revised: 01/03/2021] [Accepted: 01/07/2021] [Indexed: 12/30/2022] Open
Abstract
The azide radical (N3●) is one of the most important one-electron oxidants used extensively in radiation chemistry studies involving molecules of biological significance. Generally, it was assumed that N3● reacts in aqueous solutions only by electron transfer. However, there were several reports indicating the possibility of N3● addition in aqueous solutions to organic compounds containing double bonds. The main purpose of this study was to find an experimental approach that allows a clear assignment of the nature of obtained products either to its one-electron oxidation or its addition products. Radiolysis of water provides a convenient source of one-electron oxidizing radicals characterized by a very broad range of reduction potentials. Two inorganic radicals (SO4●-, CO3●-) and Tl2+ ions with the reduction potentials higher, and one radical (SCN)2●- with the reduction potential slightly lower than the reduction potential of N3● were selected as dominant electron-acceptors. Transient absorption spectra formed in their reactions with a series of quinoxalin-2-one derivatives were confronted with absorption spectra formed from reactions of N3● with the same series of compounds. Cases, in which the absorption spectra formed in reactions involving N3● differ from the absorption spectra formed in the reactions involving other one-electron oxidants, strongly indicate that N3● is involved in the other reaction channel such as addition to double bonds. Moreover, it was shown that high-rate constants of reactions of N3● with quinoxalin-2-ones do not ultimately prove that they are electron transfer reactions. The optimized structures of the radical cations (7-R-3-MeQ)●+, radicals (7-R-3-MeQ)● and N3● adducts at the C2 carbon atom in pyrazine moiety and their absorption spectra are reasonably well reproduced by density functional theory quantum mechanics calculations employing the ωB97XD functional combined with the Dunning's aug-cc-pVTZ correlation-consistent polarized basis sets augmented with diffuse functions.
Collapse
|
48
|
Biswas A, Mallik BS. Conformational dynamics of aqueous hydrogen peroxide from first principles molecular dynamics simulations. Phys Chem Chem Phys 2020; 22:28286-28296. [PMID: 33295373 DOI: 10.1039/d0cp05451h] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
We performed first principles molecular dynamics simulations of a relatively dilute aqueous hydrogen peroxide (H2O2) solution to examine its structural alterations and relevant dynamics upon solvation. The internal rotation of the OH groups about the O-O bond facilitates the flexible structure of H2O2. Structural calculations reveal dihedral angle fluctuations in the aqueous solution. Water molecules make stronger hydrogen bonds through the hydrogen atom of the solute than the oxygen atom leading to distinct hydrogen bonding configurations inside the first solvation shell. Time-dependent dihedral angle alterations result in conformational changes and the normalized dihedral angle distribution plot displays characteristic peaks at ∼100-120° and ∼230°, illustrating various conformational states. Within the simulation time, flexibility-induced interconversion of hydrogen peroxide gives rise to several cisoid and transoid conformers. In this study, we examine the relative population of the associated conformational states and the lifetime of the cisoid and transoid conformers from the torsion angle variations. We also determine the free energy landscape of the rotational isomerization process in H2O2 and explore two distinct energy barriers during such interconversion.
Collapse
Affiliation(s)
- Aritri Biswas
- Department of Chemistry, Indian Institute of Technology Hyderabad, Sangareddy - 502285, Telangana, India.
| | | |
Collapse
|
49
|
Samala NR, Grinberg I. Tuning of ORR activity through the stabilization of the adsorbates by hydrogen bonding with substituent groups. Phys Chem Chem Phys 2020; 22:27811-27817. [PMID: 33245314 DOI: 10.1039/d0cp04478d] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
Metallocorroles and metalloporphyrins (M-N-C) are some of the best alternative molecular catalysts for the replacement of the expensive platinum-group metals (PGM) in oxygen reduction reaction (ORR) catalysis in polymer electrolyte membrane (PEM) fuel cells. To date, Co-based corroles have shown the best performance, but still suffer from the poor stability and the toxicity of the Co metal. Mn-N-C are more stable than Co-N-C, and are also less reactive towards peroxide formation. In this work, using first-principles density functional theory calculations, we study the improvement of the Mn-based corrole ORR activity by exploiting hydrogen bonding with substituent groups to modify the adsorption energies of the ORR intermediates and obtain higher onset potential (Vonset) values. We found that by using phenyl acetic acid as a substituent, Vonset increased from 0.54 V for the unsubstituted corrole to ∼0.9 V which is competitive with the Vonset of the Co-based corroles. Our results suggest that hydrogen bonding with substituent groups should be considered in the analysis and design of the reactivity of active sites in non-PGM ORR catalysts.
Collapse
|
50
|
Delgado AAA, Sethio D, Munar I, Aviyente V, Kraka E. Local vibrational mode analysis of ion-solvent and solvent-solvent interactions for hydrated Ca 2+ clusters. J Chem Phys 2020; 153:224303. [PMID: 33317306 DOI: 10.1063/5.0034765] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
Hydrated calcium ion clusters have received considerable attention due to their essential role in biological processes such as bone development, hormone regulation, blood coagulation, and neuronal signaling. To better understand the biological role of the cation, the interactions between the Ca2+ ions and water molecules have been frequently investigated. However, a quantitative measure for the intrinsic Ca-O (ion-solvent) and intermolecular hydrogen bond (solvent-solvent) interactions has been missing so far. Here, we report a topological electron density analysis and a natural population analysis to analyze the nature of these interactions for a set of 14 hydrated calcium clusters via local mode stretching force constants obtained at the ωB97X-D/6-311++G(d,p) level of theory. The results revealed that the strength of inner Ca-O interactions for Ca(H2O)n 2+ (n = 1-8) clusters correlates with the electron density. The application of a second hydration shell to Ca(H2O)n 2+ (n = 6-8) clusters resulted in stronger Ca-O interactions where a larger electron charge transfer between lp(O) of the first hydration shell and the lower valence of Ca prevailed. The strength of the intermolecular hydrogen bonds, formed between the first and second hydration shells, became stronger when the charge transfers between hydrogen bond (HB) donors and HB acceptors were enhanced. From the local mode stretching force constants of implicitly and explicitly solvated Ca2+, we found the six-coordinated cluster to possess the strongest stabilizations, and these results prove that the intrinsic bond strength measures for Ca-O and hydrogen bond interactions form new effective tools to predict the coordination number for the hydrated calcium ion clusters.
Collapse
Affiliation(s)
- Alexis A A Delgado
- Department of Chemistry, Computational and Theoretical Chemistry Group (CATCO), Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, USA
| | - Daniel Sethio
- Department of Chemistry - BMC, Uppsala University, Husargatan 3, 75237 Uppsala, Sweden
| | - Ipek Munar
- Department of Chemistry, Boǧaziçi University, Bebek 34342, Istanbul, Turkey
| | - Viktorya Aviyente
- Department of Chemistry, Boǧaziçi University, Bebek 34342, Istanbul, Turkey
| | - Elfi Kraka
- Department of Chemistry, Computational and Theoretical Chemistry Group (CATCO), Southern Methodist University, 3215 Daniel Avenue, Dallas, Texas 75275-0314, USA
| |
Collapse
|