1
|
Suda K, Yokogawa D. Vibrational Self-Consistent Field (VSCF) and Post-VSCF Method Calculations Combined with the Reference Interaction Site Model Self-Consistent Field Method Coupled with the Constrained Spatial Electron Density Distribution: Applications to NaHCOO in Aqueous Phase. J Chem Theory Comput 2024; 20:4885-4892. [PMID: 38815984 DOI: 10.1021/acs.jctc.4c00360] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/01/2024]
Abstract
Investigating vibrational behavior in solution is crucial for understanding molecular dynamics within a solvent environment. Notably, the analysis of Raman spectra for molecules in solution is important owing to its ability to unveil intricate solute-solvent interactions. Previous studies have effectively employed frequency calculations utilizing the reference interaction site model self-consistent field method in conjunction with constrained spatial electron density distribution (RISM-SCF-cSED) to understand molecular vibrations in solution, primarily focusing on fundamental vibrational modes. However, the oversight of overtones and combination tones in these studies prompted us to combine the vibrational self-consistent field (VSCF) and vibrational second-order Mo̷ller-Plesset perturbation (VMP2) methods with RISM-SCF-cSED to address these aspects theoretically. Illustrating the efficacy of this integrated approach, we computed the Raman spectra of sodium formate (NaHCOO) in water, revealing the necessity of accounting for molecular anharmonicity in solution vibrational analysis. Our findings underscore the potency of VSCF and VMP2 in conjunction with RISM-SCF-cSED as a robust theoretical framework for such calculations.
Collapse
Affiliation(s)
- Kayo Suda
- Graduate School of Arts and Sciences, The University of Tokyo, 3-8-1 Komaba, Meguro-ku, Tokyo 153-8902, Japan
| | - Daisuke Yokogawa
- Graduate School of Arts and Sciences, The University of Tokyo, 3-8-1 Komaba, Meguro-ku, Tokyo 153-8902, Japan
| |
Collapse
|
2
|
Imamura K, Yokogawa D, Sato H. Recent developments and applications of reference interaction site model self-consistent field with constrained spatial electron density (RISM-SCF-cSED): A hybrid model of quantum chemistry and integral equation theory of molecular liquids. J Chem Phys 2024; 160:050901. [PMID: 38341702 DOI: 10.1063/5.0190116] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2023] [Accepted: 01/04/2024] [Indexed: 02/13/2024] Open
Abstract
The significance of solvent effects in electronic structure calculations has long been noted, and various methods have been developed to consider this effect. The reference interaction site model self-consistent field with constrained spatial electron density (RISM-SCF-cSED) is a hybrid model that combines the integral equation theory of molecular liquids with quantum chemistry. This method can consider the statistically convergent solvent distribution at a significantly lower cost than molecular dynamics simulations. Because the RISM theory explicitly considers the solvent structure, it performs well for systems where hydrogen bonds are formed between the solute and solvent molecules, which is a challenge for continuum solvent models. Taking advantage of being founded on the variational principle, theoretical developments have been made in calculating various properties and incorporating electron correlation effects. In this review, we organize the theoretical aspects of RISM-SCF-cSED and its distinctions from other hybrid methods involving integral equation theories. Furthermore, we carefully present its progress in terms of theoretical developments and recent applications.
Collapse
Affiliation(s)
- Kosuke Imamura
- Department of Molecular Engineering, Graduate School of Engineering, Kyoto University, Kyoto 615-8510, Japan
| | - Daisuke Yokogawa
- Graduate School of Arts and Science, The University of Tokyo, Komaba, Meguro-ku, Tokyo 153-8902, Japan
| | - Hirofumi Sato
- Department of Molecular Engineering, Graduate School of Engineering, Kyoto University, Kyoto 615-8510, Japan
- Fukui Institute for Fundamental Chemistry, Kyoto University, Kyoto 606-8103, Japan
| |
Collapse
|
3
|
Yokogawa D, Suda K. Interpretable Attribution Assignment for Octanol-Water Partition Coefficient. J Phys Chem B 2023; 127:7004-7010. [PMID: 37498912 DOI: 10.1021/acs.jpcb.3c02740] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/29/2023]
Abstract
With the increasing development of machine learning models, their credibility has become an important issue. In chemistry, attribution assignment is gaining relevance when it comes to designing molecules and debugging models. However, attention has only been paid to which atoms are important in the prediction and not to whether the attribution is reasonable. In this study, we developed a graph neural network model, a highly interpretable attribution model in chemistry, and modified the integrated gradients method. The credibility of our approach was confirmed by predicting the octanol-water partition coefficient (logP) and evaluating the three metrics (accuracy, consistency, and stability) in the attribution assignment.
Collapse
Affiliation(s)
- Daisuke Yokogawa
- Graduate School of Arts and Sciences, The University of Tokyo, 3-8-1 Komaba Meguro-ku, Tokyo 153-8902, Japan
| | - Kayo Suda
- Graduate School of Arts and Sciences, The University of Tokyo, 3-8-1 Komaba Meguro-ku, Tokyo 153-8902, Japan
| |
Collapse
|
4
|
Kramarova EP, Borisevich SS, Khamitov EM, Korlyukov AA, Dorovatovskii PV, Shagina AD, Mineev KS, Tarasenko DV, Novikov RA, Lagunin AA, Boldyrev I, Ezdoglian AA, Karpechenko NY, Shmigol TA, Baukov YI, Negrebetsky VV. Pyridine Carboxamides Based on Sulfobetaines: Design, Reactivity, and Biological Activity. Molecules 2022; 27:7542. [PMID: 36364369 PMCID: PMC9658115 DOI: 10.3390/molecules27217542] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2022] [Revised: 10/20/2022] [Accepted: 10/27/2022] [Indexed: 03/09/2024] Open
Abstract
The synthesis of the products of the 1,3-propanesultone ring opening during its interaction with amides of pyridinecarboxylic acids has been carried out. The dependence of the yield of the reaction products on the position (ortho-, meta-, para-) of the substituent in the heteroaromatic fragment and temperature condition was revealed. In contrast to the meta- and para-substituted substrates, the reaction involving ortho-derivatives at the boiling point of methanol unexpectedly led to the formation of a salt. On the basis of spectroscopic, X-Ray, and quantum-chemical calculation data, a model of the transition-state, as well as a mechanism for this alkylation reaction of pyridine carboxamides with sultone were proposed in order to explain the higher yields obtained with the nicotinamide and its N-methyl analog compared to ortho or meta parents. Based on the analysis of ESP maps, the positions of the binding sites of reagents with a potential complexing agent in space were determined. The in silico evaluation of possible biological activity showed that the synthetized compounds revealed some promising pharmacological effects and low acute toxicity.
Collapse
Affiliation(s)
- Eugene P. Kramarova
- Department of Medical Chemistry and Toxicology, Pirogov National Research Medical University, Ministry of Health of Russia, 117997 Moscow, Russia
| | - Sophia S. Borisevich
- Laboratory of Physical Chemistry, Ufa Institute of Chemistry, Ufa Federal Research Center, Russian Academy of Sciences, 450071 Ufa, Russia
| | - Edward M. Khamitov
- Laboratory of Physical Chemistry, Ufa Institute of Chemistry, Ufa Federal Research Center, Russian Academy of Sciences, 450071 Ufa, Russia
| | - Alexander A. Korlyukov
- Department of Medical Chemistry and Toxicology, Pirogov National Research Medical University, Ministry of Health of Russia, 117997 Moscow, Russia
| | | | - Anastasia D. Shagina
- Department of Medical Chemistry and Toxicology, Pirogov National Research Medical University, Ministry of Health of Russia, 117997 Moscow, Russia
| | - Konstantin S. Mineev
- Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry, Russian Academy of Sciences, 117997 Moscow, Russia
| | - Dmitri V. Tarasenko
- Department of Medical Chemistry and Toxicology, Pirogov National Research Medical University, Ministry of Health of Russia, 117997 Moscow, Russia
| | - Roman A. Novikov
- N.D. Zelinsky Institute of Organic Chemistry, Russian Academy of Sciences, 119991 Moscow, Russia
| | - Alexey A. Lagunin
- Department of Medical Chemistry and Toxicology, Pirogov National Research Medical University, Ministry of Health of Russia, 117997 Moscow, Russia
| | - Ivan Boldyrev
- Shemyakin-Ovchinnikov Institute of Bioorganic Chemistry, Russian Academy of Sciences, 117997 Moscow, Russia
| | - Aiarpi A. Ezdoglian
- Department of Medical Chemistry and Toxicology, Pirogov National Research Medical University, Ministry of Health of Russia, 117997 Moscow, Russia
| | - Natalia Yu. Karpechenko
- Department of Medical Chemistry and Toxicology, Pirogov National Research Medical University, Ministry of Health of Russia, 117997 Moscow, Russia
| | - Tatiana A. Shmigol
- Department of Medical Chemistry and Toxicology, Pirogov National Research Medical University, Ministry of Health of Russia, 117997 Moscow, Russia
| | - Yuri I. Baukov
- Department of Medical Chemistry and Toxicology, Pirogov National Research Medical University, Ministry of Health of Russia, 117997 Moscow, Russia
| | - Vadim V. Negrebetsky
- Department of Medical Chemistry and Toxicology, Pirogov National Research Medical University, Ministry of Health of Russia, 117997 Moscow, Russia
| |
Collapse
|
5
|
Yokogawa D, Suda K. Analytical second derivatives of the free energy in solution by the reference interaction site model self-consistent field explicitly including constrained spatial electron density distribution. J Chem Phys 2021; 155:204102. [PMID: 34852465 DOI: 10.1063/5.0067248] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The application of analytical derivative methods to solution systems is important because several chemical reactions occur in solution. The reference interaction site model (RISM) is one of the solvation theories used to study solution systems and has shown good performance, especially in the polar solvent systems. Although the analytical first derivative based on the RISM coupled with quantum methods (RISM-SCF) has already been derived, the analytical second derivative has not been proposed yet. Therefore, in this study, the analytical second derivative was derived using RISM-SCF explicitly including constrained spatial electron density distribution (RISM-SCF-cSED). The performance of this method was validated with the Hessian calculations of formaldehyde and para-nitroaniline in solution, and the results demonstrated that the method accurately calculated frequency values at a small computational cost.
Collapse
Affiliation(s)
- Daisuke Yokogawa
- Graduate School of Arts and Sciences, The University of Tokyo, 3-8-1 Komaba, Meguro-ku, Tokyo 153-8902, Japan
| | - Kayo Suda
- Graduate School of Arts and Sciences, The University of Tokyo, 3-8-1 Komaba, Meguro-ku, Tokyo 153-8902, Japan
| |
Collapse
|