1
|
DFT and TD-DFT study of hydrogen bonded complexes of aspartic acid and n water (n = 1 and 2). J Mol Model 2023; 29:94. [PMID: 36905452 DOI: 10.1007/s00894-023-05500-z] [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: 01/18/2023] [Accepted: 03/02/2023] [Indexed: 03/12/2023]
Abstract
CONTEXT Hydrogen bonds (HB) influence the conformational preferences of biomolecules and their optical and electronic properties. The directional interaction of molecules of water can be a prototype to understand the effects of HBs on biomolecules. Among the neurotransmitters (NT), L-aspartic acid (ASP) stands out due to its importance in health and as a precursor of several biomolecules. As it presents different functional groups and readily forms inter- and intramolecular HBs, ASP can be considered a prototype for understanding the behavior of NTs when interacting by HB with other substances. Although several theoretical studies have been performed in the past on isolated ASP and its formed complexes with water, both in gas and liquid phases, using DFT and TD-DFT formalisms, these works did not perform large basis set calculations or study electronic transitions of ASP-water complexes. We investigated the HB interactions in complexes of ASP and water molecules. The results show that the interactions between the carboxylic groups of ASP with water molecules, forming cyclic structures with two HBs, lead to more stable and less polar complexes than other conformers formed between water and the NH2 group. It was observed that there is a relationship between the deviation in the UV-Vis absorption band of the ASP and the interactions of water with the HOMO and LUMO orbitals with the stabilization/destabilization of the S1 state to the S0 of the complexes. However, in some cases, such as 1:1 complex ASP-W2, this analysis may be inaccurate due to small changes in ΔE. METHODS We studied the landscapes of the ground state surface of different conformers of isolated L-ASP and the L-ASP-(H2O)n complexes (n = 1 and 2) using the DFT formalism, with the B3LYP functional, and six different basis sets: 6-31 + + G(d,p), 6-311 + + G(d,p), D95 + + (d,p), D95V + + (d,p), cc-pVDZ, and, cc-pVTZ basis sets. The cc-pVTZ basis set provides the minimum energy of all conformers, and therefore, we performed the analysis with this basis set. We evaluated the stabilization of the ASP and complexes using the minimum ground state energy, corrected by the zero point energy and the interaction energy between the ASP and the water molecules. We also calculated the vertical electronic transitions S1 ← S0, and their properties using the TD-DFT formalism at B3LYP/cc-pVTZ level with the optimized geometries for S0 state with the same basis set. For the analysis of the vertical transitions of isolated ASP and the ASP-(H2O)n complexes, we calculated the electrostatic energy in the S0 and S1 states. We performed the calculations with the Gaussian 09 software package. We used the VMD software package to visualize the geometries and shapes of the molecule and complexes.
Collapse
|
2
|
Nakai H, Kobayashi M, Yoshikawa T, Seino J, Ikabata Y, Nishimura Y. Divide-and-Conquer Linear-Scaling Quantum Chemical Computations. J Phys Chem A 2023; 127:589-618. [PMID: 36630608 DOI: 10.1021/acs.jpca.2c06965] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
Fragmentation and embedding schemes are of great importance when applying quantum-chemical calculations to more complex and attractive targets. The divide-and-conquer (DC)-based quantum-chemical model is a fragmentation scheme that can be connected to embedding schemes. This feature article explains several DC-based schemes developed by the authors over the last two decades, which was inspired by the pioneering study of DC self-consistent field (SCF) method by Yang and Lee (J. Chem. Phys. 1995, 103, 5674-5678). First, the theoretical aspects of the DC-based SCF, electron correlation, excited-state, and nuclear orbital methods are described, followed by the two-component relativistic theory, quantum-mechanical molecular dynamics simulation, and the introduction of three programs, including DC-based schemes. Illustrative applications confirmed the accuracy and feasibility of the DC-based schemes.
Collapse
Affiliation(s)
- Hiromi Nakai
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo169-8555, Japan.,Waseda Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo169-8555, Japan
| | - Masato Kobayashi
- Department of Chemistry, Faculty of Science, Hokkaido University, Kita 10 Nishi 8, Kita-ku, Sapporo, Hokkaido060-0810, Japan.,Institute for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Kita 21 Nishi 10, Kita-ku, Sapporo, Hokkaido001-0021, Japan
| | - Takeshi Yoshikawa
- Faculty of Pharmaceutical Sciences, Toho University, 2-2-1 Miyama, Funabashi, Chiba274-8510, Japan
| | - Junji Seino
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo169-8555, Japan.,Waseda Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo169-8555, Japan
| | - Yasuhiro Ikabata
- Information and Media Center, Toyohashi University of Technology, 1-1 Hibarigaoka, Tempaku-cho, Toyohashi, Aichi441-8580, Japan.,Department of Computer Science and Engineering, Toyohashi University of Technology, 1-1 Hibarigaoka, Tempaku-cho, Toyohashi, Aichi441-8580, Japan
| | - Yoshifumi Nishimura
- Waseda Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo169-8555, Japan
| |
Collapse
|
3
|
Yan S, Wang B, Lin H. Tracking the Delocalized Proton in Concerted Proton Transfer in Bulk Water. J Chem Theory Comput 2023; 19:448-459. [PMID: 36630655 DOI: 10.1021/acs.jctc.2c01097] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
A solvated proton in water is often characterized as a charge or structural defect, and it is important to track its evolution on-the-fly in certain dynamics simulations. Previously, we introduced the proton indicator, a pseudo-atom, whose position approximates the location of the excess proton modeled as a structural defect. The proton indicator generally yields a smooth trajectory of a hydrated proton diffusing in aqueous solutions, including in the events of stepwise proton transfer via the Grotthuss mechanism; however, the proton indicator did not perform well in the events of concerted proton transfer, for which it occasionally yielded large position displacements between two successive time steps. To overcome this hurdle, we develop a new algorithm of a proton indicator with greatly enhanced performance for concerted proton transfer in bulk water. A protocol is proposed to exhaustively explore the hydrogen-bonding network of the water wires over which the excess proton is delocalized and to properly account for the contributions of the water molecules in this network as the geometry evolves. The new proton indicator (called Indicator 2.0) is assessed in dynamics simulations of an excess proton in bulk water and in specially constructed model systems of more complex architectures. The results demonstrate that the new indicator yields a smooth trajectory in both stepwise and concerted proton transfers.
Collapse
Affiliation(s)
- Shengheng Yan
- State Key Laboratory of Physical Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen360015P. R. China
| | - Binju Wang
- State Key Laboratory of Physical Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen360015P. R. China
| | - Hai Lin
- Department of Chemistry, CB 194, University of Colorado Denver, P.O. Box 173364, Denver, Colorado80217, United States
| |
Collapse
|
4
|
Muñoz-Santiburcio D. Accurate diffusion coefficients of the excess proton and hydroxide in water via extensive ab initio simulations with different schemes. J Chem Phys 2022; 157:024504. [PMID: 35840376 DOI: 10.1063/5.0093958] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Despite its simple molecular formula, obtaining an accurate in silico description of water is far from straightforward. Many of its very peculiar properties are quite elusive, and in particular, obtaining good estimations of the diffusion coefficients of the solvated proton and hydroxide at a reasonable computational cost has been an unsolved challenge until now. Here, I present extensive results of several unusually long ab initio molecular dynamics (MD) simulations employing different combinations of the Born-Oppenheimer and second-generation Car-Parrinello MD propagation methods with different ensembles (NVE and NVT) and thermostats, which show that these methods together with the RPBE-D3 functional provide a very accurate estimation of the diffusion coefficients of the solvated H3O+ and OH- ions, together with an extremely accurate description of several properties of neutral water (such as the structure of the liquid and its diffusion and shear viscosity coefficients). In addition, I show that the estimations of DH3O+ and DOH- depend dramatically on the simulation length, being necessary to reach timescales in the order of hundreds of picoseconds to obtain reliable results.
Collapse
Affiliation(s)
- Daniel Muñoz-Santiburcio
- CIC nanoGUNE BRTA, Tolosa Hiribidea 76, 20018 San Sebastián, Spain and Instituto de Fusión Nuclear "Guillermo Velarde," Universidad Politécnica de Madrid, C/ José Gutiérrez Abascal 2, 28006 Madrid, Spain
| |
Collapse
|
5
|
Liu R, Zhang C, Liang X, Liu J, Wu X, Chen M. Structural and Dynamic Properties of Solvated Hydroxide and Hydronium Ions in Water from Ab Initio Modeling. J Chem Phys 2022; 157:024503. [DOI: 10.1063/5.0094944] [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/14/2022] Open
Abstract
Predicting the asymmetric structure and dynamics of solvated hydroxide and hydronium in water has been a challenging task from ab initio molecular dynamics (AIMD). The difficulty mainly comes from a lack of accurate and efficient exchange-correlation functional in elucidating the amphiphilic nature and the ubiquitous proton transfer behaviors of the two ions. By adopting the strongly-constrained and appropriately normed (SCAN) meta-GGA functional in AIMD simulations, we systematically examine the amphiphilic properties, the solvation structures, the electronic structures, and the dynamic properties of the two water ions. In particular, we compare these results to those predicted by the PBE0-TS functional, which is an accurate yet computationally more expensive exchange-correlation functional. We demonstrate that the general-purpose SCAN functional provides a reliable choice in describing the two water ions. Specifically, in the SCAN picture of water ions, the appearance of the fourth and fifth hydrogen bonds near hydroxide stabilizes the pot-like shape solvation structure and suppresses the structural diffusion, while the hydronium stably donates three hydrogen bonds to its neighbors. We apply a detailed analysis of the proton transfer mechanism of the two ions and find the two ions exhibit substantially different proton transfer patterns. Our AIMD simulations indicate hydroxide diffuses slower than hydronium in water, which is consistent with the experiments.
Collapse
Affiliation(s)
| | | | | | | | - Xifan Wu
- Physics, Temple University, United States of America
| | - Mohan Chen
- College of Engineering, Peking University, China
| |
Collapse
|
6
|
Charzewski Ł, Krzyśko KA, Lesyng B. Exploring Covalent Docking Mechanisms of Boron-Based Inhibitors to Class A, C and D β-Lactamases Using Time-dependent Hybrid QM/MM Simulations. Front Mol Biosci 2021; 8:633181. [PMID: 34434961 PMCID: PMC8380965 DOI: 10.3389/fmolb.2021.633181] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2020] [Accepted: 07/19/2021] [Indexed: 11/13/2022] Open
Abstract
Recently, molecular covalent docking has been extensively developed to design new classes of inhibitors that form chemical bonds with their biological targets. This strategy for the design of such inhibitors, in particular boron-based inhibitors, holds great promise for the vast family of β-lactamases produced, inter alia, by Gram-negative antibiotic-resistant bacteria. However, the description of covalent docking processes requires a quantum-mechanical approach, and so far, only a few studies of this type have been presented. This study accurately describes the covalent docking process between two model inhibitors - representing two large families of inhibitors based on boronic-acid and bicyclic boronate scaffolds, and three β-lactamases which belong to the A, C, and D classes. Molecular fragments containing boron can be converted from a neutral, trigonal, planar state with sp2 hybridization to the anionic, tetrahedral sp3 state in a process sometimes referred to as morphing. This study applies multi-scale modeling methods, in particular, the hybrid QM/MM approach which has predictive power reaching well beyond conventional molecular modeling. Time-dependent QM/MM simulations indicated several structural changes and geometric preferences, ultimately leading to covalent docking processes. With current computing technologies, this approach is not computationally expensive, can be used in standard molecular modeling and molecular design works, and can effectively support experimental research which should allow for a detailed understanding of complex processes important to molecular medicine. In particular, it can support the rational design of covalent boron-based inhibitors for β-lactamases as well as for many other enzyme systems of clinical relevance, including SARS-CoV-2 proteins.
Collapse
Affiliation(s)
| | | | - Bogdan Lesyng
- Department of Biophysics, Faculty of Physics, University of Warsaw, Warsaw, Poland
| |
Collapse
|
7
|
Ono J, Imai M, Nishimura Y, Nakai H. Hydroxide Ion Carrier for Proton Pumps in Bacteriorhodopsin: Primary Proton Transfer. J Phys Chem B 2020; 124:8524-8539. [DOI: 10.1021/acs.jpcb.0c05507] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Junichi Ono
- Waseda Research Institute for Science and Engineering (WISE), Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo 169-8555, Japan
- Elements Strategy Initiative for Catalysts & Batteries (ESICB), Kyoto University, 1-30 Goryo-Ohara, Nishikyo-ku, Kyoto 615-8245, Japan
| | - Minori Imai
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo 169-8555, Japan
| | - Yoshifumi Nishimura
- Waseda Research Institute for Science and Engineering (WISE), Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo 169-8555, Japan
| | - Hiromi Nakai
- Waseda Research Institute for Science and Engineering (WISE), Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo 169-8555, Japan
- Elements Strategy Initiative for Catalysts & Batteries (ESICB), Kyoto University, 1-30 Goryo-Ohara, Nishikyo-ku, Kyoto 615-8245, Japan
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo 169-8555, Japan
| |
Collapse
|
8
|
Ono J, Nakai H. Weighted histogram analysis method for multiple short-time metadynamics simulations. Chem Phys Lett 2020. [DOI: 10.1016/j.cplett.2020.137384] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
|