1
|
Li S, Xie BB, Yin BW, Liu L, Shen L, Fang WH. Construction of Highly Accurate Machine Learning Potential Energy Surfaces for Excited-State Dynamics Simulations Based on Low-Level Data Sets. J Phys Chem A 2024; 128:5516-5524. [PMID: 38954640 DOI: 10.1021/acs.jpca.4c02028] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/04/2024]
Abstract
Machine learning is capable of effectively predicting the potential energies of molecules in the presence of high-quality data sets. Its application in the construction of ground- and excited-state potential energy surfaces is attractive to accelerate nonadiabatic molecular dynamics simulations of photochemical reactions. Because of the huge computational cost of excited-state electronic structure calculations, the construction of a high-quality data set becomes a bottleneck. In the present work, we first built two data sets. One was obtained from surface hopping dynamics simulations at the semiempirical OM2/MRCI level. Another was extracted from the dynamics trajectories at the CASSCF level, which was reported previously. The ground- and excited-state potential energy surfaces of ethylene-bridged azobenzene at the CASSCF computational level were constructed based on the former low-level data set. Although non-neural network machine learning methods can achieve good or modest performance during the training process, only neural network models provide reliable predictions on the latter external test data set. The BPNN and SchNet combined with the Δ-ML scheme and the force term in the loss functions are recommended for dynamics simulations. Then, we performed excited-state dynamics simulations of the photoisomerization of ethylene-bridged azobenzene on machine learning potential energy surfaces. Compared with the lifetimes of the first excited state (S1) estimated at different computational levels, our results on the E isomer are in good agreement with the high-level estimation. However, the overestimation of the Z isomer is unimproved. It suggests that smaller errors during the training process do not necessarily translate to more accurate predictions on high-level potential energies or better performance on nonadiabatic dynamics simulations, at least in the present case.
Collapse
Affiliation(s)
- Shuai Li
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
| | - Bin-Bin Xie
- Hangzhou Institute of Advanced Studies, Zhejiang Normal University, Hangzhou 311231, Zhejiang, P. R. China
| | - Bo-Wen Yin
- Hangzhou Institute of Advanced Studies, Zhejiang Normal University, Hangzhou 311231, Zhejiang, P. R. China
| | - Lihong Liu
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
| | - Lin Shen
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
- Yantai-Jingshi Institute of Material Genome Engineering, Yantai 265505, Shandong, P. R. China
| | - Wei-Hai Fang
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
- Shandong Laboratory of Yantai Advanced Materials and Green Manufacturing, Yantai 264006, Shandong, P. R. China
| |
Collapse
|
2
|
Martyka M, Jankowska J. New insights into the photocyclization reaction of a popular diarylethene switch: a nonadiabatic molecular dynamics study. Phys Chem Chem Phys 2024; 26:13383-13394. [PMID: 38646878 DOI: 10.1039/d3cp06256b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/23/2024]
Abstract
Diarylethene (DAE) molecular switches have continued to attract the attention of researchers for over 20 years. Their remarkable photophysical properties endow them with countless applications in photonics and molecular technologies. However, despite extensive experimental and theoretical research, the mechanism of DAE photoswitching is not yet fully rationalized. In this work, we investigate the ring closure dynamics of a popular DAE switch, 1,2-bis(3-methyl-5-phenyl-2 thienyl)perfluorocyclopentene (PT), using nonadiabatic molecular dynamics (NAMD) simulations. Employing the fewest switches surface hopping protocol, along with the semi-empirical multireference ODM2/MRCI-SD method, we investigate possible reaction pathways for this photoprocess, as well as their timescales and resulting photoproducts. Furthermore, using a dynamic configuration-space sampling procedure, we elucidate the role of triplet states in the photocyclization of PT, supporting available experimental data for the closely related DMPT molecule, which indicate an ultrafast intersystem crossing (ISC) transition competing with the singlet-driven photoswitching reaction. Our findings not only corroborate experimental studies on DAE switches, but also provide new mechanistic insights into the potential use in the rational design of DAE switches tailored for specific technological applications.
Collapse
Affiliation(s)
- Mikołaj Martyka
- Faculty of Chemistry, University of Warsaw, Pasteura 1, Warsaw, 02-093, Poland.
- Interdisciplinary Doctoral School, University of Warsaw, Dobra 56/66, Warsaw, 00-312, Poland
| | - Joanna Jankowska
- Faculty of Chemistry, University of Warsaw, Pasteura 1, Warsaw, 02-093, Poland.
| |
Collapse
|
3
|
Jorner K, Pollice R, Lavigne C, Aspuru-Guzik A. Ultrafast Computational Screening of Molecules with Inverted Singlet-Triplet Energy Gaps Using the Pariser-Parr-Pople Semiempirical Quantum Chemistry Method. J Phys Chem A 2024; 128:2445-2456. [PMID: 38485448 PMCID: PMC10983003 DOI: 10.1021/acs.jpca.3c06357] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2023] [Revised: 02/23/2024] [Accepted: 02/23/2024] [Indexed: 04/04/2024]
Abstract
Molecules with an inverted energy gap between their first singlet and triplet excited states have promising applications in the next generation of organic light-emitting diode (OLED) materials. Unfortunately, such molecules are rare, and only a handful of examples are currently known. High-throughput virtual screening could assist in finding novel classes of these molecules, but current efforts are hampered by the high computational cost of the required quantum chemical methods. We present a method based on the semiempirical Pariser-Parr-Pople theory augmented by perturbation theory and show that it reproduces inverted gaps at a fraction of the cost of currently employed excited-state calculations. Our study paves the way for ultrahigh-throughput virtual screening and inverse design to accelerate the discovery and development of this new generation of OLED materials.
Collapse
Affiliation(s)
- Kjell Jorner
- Institute
of Chemical and Bioengineering, Department of Chemistry and Applied
Biosciences, ETH Zürich, Vladimir-Prelog-Weg 1, Zürich CH-8093, Switzerland
- Department
of Chemistry and Chemical Engineering, Chalmers
University of Technology, Kemigården 4, Gothenburg SE-41258, Sweden
- Chemical
Physics Theory Group, Department of Chemistry, University of Toronto, 80 St. George Street, Toronto M5S 3H6, Canada
- Department
of Computer Science, University of Toronto, 40 St. George Street, Toronto M5S 2E4, Canada
| | - Robert Pollice
- Chemical
Physics Theory Group, Department of Chemistry, University of Toronto, 80 St. George Street, Toronto M5S 3H6, Canada
- Department
of Computer Science, University of Toronto, 40 St. George Street, Toronto M5S 2E4, Canada
- Stratingh
Institute for Chemistry, University of Groningen, Nijenborgh 4, Groningen 9747, AG, The Netherlands
| | - Cyrille Lavigne
- Chemical
Physics Theory Group, Department of Chemistry, University of Toronto, 80 St. George Street, Toronto M5S 3H6, Canada
- Department
of Computer Science, University of Toronto, 40 St. George Street, Toronto M5S 2E4, Canada
| | - Alán Aspuru-Guzik
- Chemical
Physics Theory Group, Department of Chemistry, University of Toronto, 80 St. George Street, Toronto M5S 3H6, Canada
- Department
of Computer Science, University of Toronto, 40 St. George Street, Toronto M5S 2E4, Canada
- Department
of Chemical Engineering & Applied Chemistry, University of Toronto, 200 College Street, Toronto M5S 3E5, Canada
- Department
of Materials Science & Engineering, University of Toronto, 184 College Street, Toronto M5S 3E4, Canada
- Vector
Institute for Artificial Intelligence, 661 University Ave. Suite 710, Toronto M5G 1M1, Canada
- Lebovic
Fellow, Canadian Institute for Advanced
Research (CIFAR), 661
University Avenue, Toronto M5G 1M1, Canada
- Acceleration
Consortium, University of Toronto, 700 University Avenue, Toronto M5G 1Z5, Canada
| |
Collapse
|
4
|
Salvadori G, Mazzeo P, Accomasso D, Cupellini L, Mennucci B. Deciphering Photoreceptors Through Atomistic Modeling from Light Absorption to Conformational Response. J Mol Biol 2024; 436:168358. [PMID: 37944793 DOI: 10.1016/j.jmb.2023.168358] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2023] [Revised: 10/28/2023] [Accepted: 11/02/2023] [Indexed: 11/12/2023]
Abstract
In this review, we discuss the successes and challenges of the atomistic modeling of photoreceptors. Throughout our presentation, we integrate explanations of the primary methodological approaches, ranging from quantum mechanical descriptions to classical enhanced sampling methods, all while providing illustrative examples of their practical application to specific systems. To enhance the effectiveness of our analysis, our primary focus has been directed towards the examination of applications across three distinct photoreceptors. These include an example of Blue Light-Using Flavin (BLUF) domains, a bacteriophytochrome, and the orange carotenoid protein (OCP) employed by cyanobacteria for photoprotection. Particular emphasis will be placed on the pivotal role played by the protein matrix in fine-tuning the initial photochemical event within the embedded chromophore. Furthermore, we will investigate how this localized perturbation initiates a cascade of events propagating from the binding pocket throughout the entire protein structure, thanks to the intricate network of interactions between the chromophore and the protein.
Collapse
Affiliation(s)
- Giacomo Salvadori
- Department of Chemistry and Industrial Chemistry, University of Pisa, 56124 Pisa, Italy
| | - Patrizia Mazzeo
- Department of Chemistry and Industrial Chemistry, University of Pisa, 56124 Pisa, Italy
| | - Davide Accomasso
- Department of Chemistry and Industrial Chemistry, University of Pisa, 56124 Pisa, Italy
| | - Lorenzo Cupellini
- Department of Chemistry and Industrial Chemistry, University of Pisa, 56124 Pisa, Italy
| | - Benedetta Mennucci
- Department of Chemistry and Industrial Chemistry, University of Pisa, 56124 Pisa, Italy
| |
Collapse
|
5
|
Dral PO, Ge F, Hou YF, Zheng P, Chen Y, Barbatti M, Isayev O, Wang C, Xue BX, Pinheiro Jr M, Su Y, Dai Y, Chen Y, Zhang L, Zhang S, Ullah A, Zhang Q, Ou Y. MLatom 3: A Platform for Machine Learning-Enhanced Computational Chemistry Simulations and Workflows. J Chem Theory Comput 2024; 20:1193-1213. [PMID: 38270978 PMCID: PMC10867807 DOI: 10.1021/acs.jctc.3c01203] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2023] [Revised: 12/29/2023] [Accepted: 01/03/2024] [Indexed: 01/26/2024]
Abstract
Machine learning (ML) is increasingly becoming a common tool in computational chemistry. At the same time, the rapid development of ML methods requires a flexible software framework for designing custom workflows. MLatom 3 is a program package designed to leverage the power of ML to enhance typical computational chemistry simulations and to create complex workflows. This open-source package provides plenty of choice to the users who can run simulations with the command-line options, input files, or with scripts using MLatom as a Python package, both on their computers and on the online XACS cloud computing service at XACScloud.com. Computational chemists can calculate energies and thermochemical properties, optimize geometries, run molecular and quantum dynamics, and simulate (ro)vibrational, one-photon UV/vis absorption, and two-photon absorption spectra with ML, quantum mechanical, and combined models. The users can choose from an extensive library of methods containing pretrained ML models and quantum mechanical approximations such as AIQM1 approaching coupled-cluster accuracy. The developers can build their own models using various ML algorithms. The great flexibility of MLatom is largely due to the extensive use of the interfaces to many state-of-the-art software packages and libraries.
Collapse
Affiliation(s)
- Pavlo O. Dral
- State
Key Laboratory of Physical Chemistry of Solid Surfaces, College of
Chemistry and Chemical Engineering, and Innovation Laboratory for
Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
- Fujian
Provincial Key Laboratory of Theoretical and Computational Chemistry, Xiamen, Fujian 361005, China
| | - Fuchun Ge
- State
Key Laboratory of Physical Chemistry of Solid Surfaces, College of
Chemistry and Chemical Engineering, and Innovation Laboratory for
Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
- Fujian
Provincial Key Laboratory of Theoretical and Computational Chemistry, Xiamen, Fujian 361005, China
| | - Yi-Fan Hou
- State
Key Laboratory of Physical Chemistry of Solid Surfaces, College of
Chemistry and Chemical Engineering, and Innovation Laboratory for
Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
- Fujian
Provincial Key Laboratory of Theoretical and Computational Chemistry, Xiamen, Fujian 361005, China
| | - Peikun Zheng
- State
Key Laboratory of Physical Chemistry of Solid Surfaces, College of
Chemistry and Chemical Engineering, and Innovation Laboratory for
Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
- Fujian
Provincial Key Laboratory of Theoretical and Computational Chemistry, Xiamen, Fujian 361005, China
| | - Yuxinxin Chen
- State
Key Laboratory of Physical Chemistry of Solid Surfaces, College of
Chemistry and Chemical Engineering, and Innovation Laboratory for
Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
- Fujian
Provincial Key Laboratory of Theoretical and Computational Chemistry, Xiamen, Fujian 361005, China
| | - Mario Barbatti
- Aix
Marseille University, CNRS, ICR, Marseille 13013, France
- Institut
Universitaire de France, Paris 75231, France
| | - Olexandr Isayev
- Department
of Chemistry, Carnegie Mellon University, Pittsburgh, Pennsylvania15213, United States
| | - Cheng Wang
- State
Key Laboratory of Physical Chemistry of Solid Surfaces, College of
Chemistry and Chemical Engineering, and Innovation Laboratory for
Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
- iChem, Xiamen University, Xiamen, Fujian 361005, China
| | - Bao-Xin Xue
- State
Key Laboratory of Physical Chemistry of Solid Surfaces, College of
Chemistry and Chemical Engineering, and Innovation Laboratory for
Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
- Fujian
Provincial Key Laboratory of Theoretical and Computational Chemistry, Xiamen, Fujian 361005, China
| | - Max Pinheiro Jr
- Aix
Marseille University, CNRS, ICR, Marseille 13013, France
| | - Yuming Su
- State
Key Laboratory of Physical Chemistry of Solid Surfaces, College of
Chemistry and Chemical Engineering, and Innovation Laboratory for
Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
- iChem, Xiamen University, Xiamen, Fujian 361005, China
| | - Yiheng Dai
- State
Key Laboratory of Physical Chemistry of Solid Surfaces, College of
Chemistry and Chemical Engineering, and Innovation Laboratory for
Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
- iChem, Xiamen University, Xiamen, Fujian 361005, China
| | - Yangtao Chen
- State
Key Laboratory of Physical Chemistry of Solid Surfaces, College of
Chemistry and Chemical Engineering, and Innovation Laboratory for
Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
- iChem, Xiamen University, Xiamen, Fujian 361005, China
| | - Lina Zhang
- State
Key Laboratory of Physical Chemistry of Solid Surfaces, College of
Chemistry and Chemical Engineering, and Innovation Laboratory for
Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
- Fujian
Provincial Key Laboratory of Theoretical and Computational Chemistry, Xiamen, Fujian 361005, China
| | - Shuang Zhang
- State
Key Laboratory of Physical Chemistry of Solid Surfaces, College of
Chemistry and Chemical Engineering, and Innovation Laboratory for
Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
- Fujian
Provincial Key Laboratory of Theoretical and Computational Chemistry, Xiamen, Fujian 361005, China
| | - Arif Ullah
- School
of Physics and Optoelectronic Engineering, Anhui University, Hefei230601, China
| | - Quanhao Zhang
- State
Key Laboratory of Physical Chemistry of Solid Surfaces, College of
Chemistry and Chemical Engineering, and Innovation Laboratory for
Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
- Fujian
Provincial Key Laboratory of Theoretical and Computational Chemistry, Xiamen, Fujian 361005, China
| | - Yanchi Ou
- State
Key Laboratory of Physical Chemistry of Solid Surfaces, College of
Chemistry and Chemical Engineering, and Innovation Laboratory for
Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen, Fujian 361005, China
- Fujian
Provincial Key Laboratory of Theoretical and Computational Chemistry, Xiamen, Fujian 361005, China
| |
Collapse
|
6
|
Pieri E, Weingart O, Huix-Rotllant M, Ledentu V, Garavelli M, Ferré N. Modeling pH-Dependent Biomolecular Photochemistry. J Chem Theory Comput 2024; 20:842-855. [PMID: 38198619 DOI: 10.1021/acs.jctc.3c00980] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2024]
Abstract
The tuning mechanism of pH can be extremely challenging to model computationally in complex biological systems, especially with respect to the photochemical properties. This article reports a protocol aimed at modeling pH-dependent photodynamics using a combination of constant-pH molecular dynamics and semiclassical nonadiabatic molecular dynamics simulations. With retinal photoisomerization in Anabaena sensory rhodopsin (ASR) as a testbed, we show that our protocol produces pH-dependent photochemical properties, such as the isomerization quantum yield or decay rates. We decompose our results into single-titrated residue contributions, identifying some key tuning amino acids. Additionally, we assess the validity of the single protonation state picture to represent the system at a given pH and propose the most populated protein charge state as a compromise between cost and accuracy.
Collapse
Affiliation(s)
- Elisa Pieri
- Aix-Marseille Univ, CNRS, Institut de Chimie Radicalaire, 13013 Marseille, France
| | - Oliver Weingart
- Faculty of Mathematics and Natural Sciences, Institute for Theoretical and Computational Chemistry, Heinrich Heine University, Universitätsstr. 1, 40225 Düsseldorf, Germany
| | - Miquel Huix-Rotllant
- Aix-Marseille Univ, CNRS, Institut de Chimie Radicalaire, 13013 Marseille, France
| | - Vincent Ledentu
- Aix-Marseille Univ, CNRS, Institut de Chimie Radicalaire, 13013 Marseille, France
| | - Marco Garavelli
- Dipartimento di Chimica Industriale "Toso Montanari", Università degli Studi di Bologna, Viale del Risorgimento, 4, 40136 Bologna, Italy
| | - Nicolas Ferré
- Aix-Marseille Univ, CNRS, Institut de Chimie Radicalaire, 13013 Marseille, France
| |
Collapse
|
7
|
Janoš J, Slavíček P. What Controls the Quality of Photodynamical Simulations? Electronic Structure Versus Nonadiabatic Algorithm. J Chem Theory Comput 2023; 19:8273-8284. [PMID: 37939301 PMCID: PMC10688183 DOI: 10.1021/acs.jctc.3c00908] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2023] [Revised: 10/18/2023] [Accepted: 10/19/2023] [Indexed: 11/10/2023]
Abstract
The field of nonadiabatic dynamics has matured over the last decade with a range of algorithms and electronic structure methods available at the moment. While the community currently focuses more on developing and benchmarking new nonadiabatic dynamics algorithms, the underlying electronic structure controls the outcome of nonadiabatic simulations. Yet, the electronic-structure sensitivity analysis is typically neglected. In this work, we present a sensitivity analysis of the nonadiabatic dynamics of cyclopropanone to electronic structure methods and nonadiabatic dynamics algorithms. In particular, we compare wave function-based CASSCF, FOMO-CASCI, MS- and XMS-CASPT2, density-functional REKS, and semiempirical MRCI-OM3 electronic structure methods with the Landau-Zener surface hopping, fewest switches surface hopping, and ab initio multiple spawning with informed stochastic selection algorithms. The results clearly demonstrate that the electronic structure choice significantly influences the accuracy of nonadiabatic dynamics for cyclopropanone even when the potential energy surfaces exhibit qualitative and quantitative similarities. Thus, selecting the electronic structure solely on the basis of the mapping of potential energy surfaces can be misleading. Conversely, we observe no discernible differences in the performance of the nonadiabatic dynamics algorithms across the various methods. Based on the above results, we discuss the present-day practice in computational photodynamics.
Collapse
Affiliation(s)
- Jiří Janoš
- Department of Physical Chemistry, University of Chemistry and Technology, Technická 5, 16628 Prague 6, Czech Republic
| | - Petr Slavíček
- Department of Physical Chemistry, University of Chemistry and Technology, Technická 5, 16628 Prague 6, Czech Republic
| |
Collapse
|
8
|
Budiutama G, Li R, Manzhos S, Ihara M. Hybrid Density Functional Tight Binding (DFTB)─Molecular Mechanics Approach for a Low-Cost Expansion of DFTB Applicability. J Chem Theory Comput 2023. [PMID: 37450317 DOI: 10.1021/acs.jctc.3c00310] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/18/2023]
Abstract
The density functional-based tight binding (DFTB) method has seen a rise in adoption for materials modeling, as it offers significant improvement in scalability with accuracy comparable to the density functional theory (DFT) when good parameterizations exist. The cost reduction in DFTB compared to DFT is achieved by the pre-parameterization of the elements of the Hamiltonian matrix as well as the repulsion potential between all pairs of atoms. Parameterization for new systems with accuracies competitive with DFT in specific applications requires specialized manpower and computational resources. This prevents the application of the DFTB method to systems for which it was not parameterized. In this work, we explore an approach to address the problem of missing parameters of DFTB by modeling the interactions with missing Slater-Koster parameters with an interatomic interaction potential. When the distance between two atoms modeled at the force-field level is sufficiently large, the approach results in accurate structural and electronic properties. The resulting calculation is therefore a hybrid between DFTB and molecular mechanics, a pure DFTB for atoms with a complete set of interatomic parameterizations, and a mix between DFTB and molecular mechanics for atoms with a missing interatomic parameterization. The approach is expected to be particularly useful for hybrid materials and interfaces. The method is tested on the examples of 2D materials, mixed oxides, and a large-scale calculation of an oxide-oxide interface.
Collapse
Affiliation(s)
- Gekko Budiutama
- School of Materials and Chemical Technology, Tokyo Institute of Technology, Ookayama 2-12-1, Meguro-ku, Tokyo 152-8552 Japan
| | - Ruicheng Li
- School of Materials and Chemical Technology, Tokyo Institute of Technology, Ookayama 2-12-1, Meguro-ku, Tokyo 152-8552 Japan
| | - Sergei Manzhos
- School of Materials and Chemical Technology, Tokyo Institute of Technology, Ookayama 2-12-1, Meguro-ku, Tokyo 152-8552 Japan
| | - Manabu Ihara
- School of Materials and Chemical Technology, Tokyo Institute of Technology, Ookayama 2-12-1, Meguro-ku, Tokyo 152-8552 Japan
| |
Collapse
|
9
|
Wang Z, Neese F. Development of NOTCH, an all-electron, beyond-NDDO semiempirical method: Application to diatomic molecules. J Chem Phys 2023; 158:2889026. [PMID: 37154284 DOI: 10.1063/5.0141686] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/07/2023] [Accepted: 04/18/2023] [Indexed: 05/10/2023] Open
Abstract
In this work, we develop a new semiempirical method, dubbed NOTCH (Natural Orbital Tied Constructed Hamiltonian). Compared to existing semiempirical methods, NOTCH is less empirical in its functional form as well as parameterization. Specifically, in NOTCH, (1) the core electrons are treated explicitly; (2) the nuclear-nuclear repulsion term is calculated analytically, without any empirical parameterization; (3) the contraction coefficients of the atomic orbital (AO) basis depend on the coordinates of the neighboring atoms, which allows the size of AOs to depend on the molecular environment, despite the fact that a minimal basis set is used; (4) the one-center integrals of free atoms are derived from scalar relativistic multireference equation-of-motion coupled cluster calculations instead of empirical fitting, drastically reducing the number of necessary empirical parameters; (5) the (AA|AB) and (AB|AB)-type two-center integrals are explicitly included, going beyond the neglect of differential diatomic overlap approximation; and (6) the integrals depend on the atomic charges, effectively mimicking the "breathing" of AOs when the atomic charge varies. For this preliminary report, the model has been parameterized for the elements H-Ne, giving only 8 empirical global parameters. Preliminary results on the ionization potentials, electron affinities, and excitation energies of atoms and diatomic molecules, as well as the equilibrium geometries, vibrational frequencies dipole moments, and bond dissociation energies of diatomic molecules, show that the accuracy of NOTCH rivals or exceeds those of popular semiempirical methods (including PM3, PM7, OM2, OM3, GFN-xTB, and GFN2-xTB) as well as the cost-effective ab initio method Hartree-Fock-3c.
Collapse
Affiliation(s)
- Zikuan Wang
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Frank Neese
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
10
|
Ma J, Zhao D, Yu L, Jiang C, Lan Z, Li F. Simultaneously improving the efficiencies of photo- and thermal isomerization of an oxindole-based light-driven molecular rotary motor by a structural redesign. Phys Chem Chem Phys 2023; 25:12800-12809. [PMID: 37129050 DOI: 10.1039/d3cp00559c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/03/2023]
Abstract
We designed a novel highly efficient light-driven molecular rotary motor theoretically by using electronic structure calculations and nonadiabatic dynamics simulations, and it showed excellent performance for both photo- and thermal isomerization processes simultaneously. By the small structural modification based on 3-(2,7-dimethyl-2,3-dihydro-1H-inden-1-ylidene)-1-methylindolin-2-one (DDIYM) synthesized by Feringa et al. recently, an oxindole-based light-driven molecular rotary motor, 3-(1,5-dimethyl-4,5-dihydrocyclopenta[b]pyrrol-6(1H)-ylidene)-1-methylindolin-2-one (DDPYM), is proposed, which displays a significant electronic push-pull character and weak steric hindrance for double-bond isomerization. The newly designed motor DDPYM shows a remarkable improvement of the quantum yield for both EP → ZM and ZP → EM photoisomerization processes, compared to the original motor DDIYM. Furthermore, the rotary motion in photoisomerization processes of DDPYM behaves more like a pure axial rotational motion approximately, while that of DDIYM is an obvious precessional motion. The weakness of the steric hindrance reduces the energy barriers of the thermal helix EM → EP and ZM → ZP inversion steps, and would accelerate two ground-state isomerization steps significantly. Our results confirm the feasibility of simultaneously improving the efficiencies of photo- and thermal isomerization of oxindole-based light-driven molecular rotary motors and this design idea sheds light on the future development of more efficient molecular motors.
Collapse
Affiliation(s)
- Jianzheng Ma
- Ministry of Education Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, Shaanxi Province Key Laboratory of Quantum Information and Quantum Optoelectronic Devices, School of Physics, Xi'an Jiaotong University, Xian 710049, China.
| | - Di Zhao
- Ministry of Education Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, Shaanxi Province Key Laboratory of Quantum Information and Quantum Optoelectronic Devices, School of Physics, Xi'an Jiaotong University, Xian 710049, China.
| | - Le Yu
- Key Laboratory of Synthetic and Natural Functional Molecule Chemistry of Ministry of Education, College of Chemistry & Materials Science, Shaanxi Key Laboratory of Physico-Inorganic Chemistry, Northwest University, Xi'an 710127, China
| | - Chenwei Jiang
- Ministry of Education Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, Shaanxi Province Key Laboratory of Quantum Information and Quantum Optoelectronic Devices, School of Physics, Xi'an Jiaotong University, Xian 710049, China.
| | - Zhenggang Lan
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety & MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, China; School of Environment, South China Normal University, Guangzhou 510006, China.
| | - Fuli Li
- Ministry of Education Key Laboratory for Nonequilibrium Synthesis and Modulation of Condensed Matter, Shaanxi Province Key Laboratory of Quantum Information and Quantum Optoelectronic Devices, School of Physics, Xi'an Jiaotong University, Xian 710049, China.
| |
Collapse
|
11
|
Martyka M, Jankowska J. Nonadiabatic molecular dynamics study of a complete photoswitching cycle for a full-size diarylethene system. J Photochem Photobiol A Chem 2023. [DOI: 10.1016/j.jphotochem.2022.114513] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
|
12
|
Ong AWW, Cao SY, Kwek LC. An improved parameterization procedure for NDDO-descendant semi-empirical methods. J Mol Model 2023; 29:118. [PMID: 36977949 PMCID: PMC10050048 DOI: 10.1007/s00894-023-05499-3] [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/30/2022] [Accepted: 03/02/2023] [Indexed: 03/30/2023]
Abstract
CONCEPT MNDO-based semi-empirical methods in quantum chemistry have found widespread application in the modelling of large and complex systems. A method for the analytic evaluation of first and second derivatives of molecular properties against semi-empirical parameters in MNDO-based NDDO-descendant models is presented, and the resultant parameter Hessian is compared against the approximant currently used in parameterization for the PMx models. METHODS As a proof of concept, the exact parameter Hessian is employed in a limited reparameterization of MNDO for the elements C, H, N, O and F using 1206 molecules for reference data (heats of formation, ionization energies, dipole moments and reference geometries). The correctness of our MNDO implementation was verified by comparing the calculated molecular properties with the MOPAC program.
Collapse
Affiliation(s)
- Adrian Wee Wen Ong
- NUS High School of Mathematics and Science, 20 Clementi Avenue 1, 129957, Singapore, Singapore.
- Centre for Quantum Technologies, National University of Singapore, 117543, Singapore, Singapore.
| | - Steve Yueran Cao
- NUS High School of Mathematics and Science, 20 Clementi Avenue 1, 129957, Singapore, Singapore
- Centre for Quantum Technologies, National University of Singapore, 117543, Singapore, Singapore
| | - Leong Chuan Kwek
- Centre for Quantum Technologies, National University of Singapore, 117543, Singapore, Singapore
- MajuLab, CNRS-UNS-NUS-NTU International Joint Research Unit, UMI 3654, 117543, Singapore, Singapore
- National Institute of Education, Nanyang Technological University, 1 Nanyang Walk, 637616, Singapore, Singapore
- School of Electrical and Electronic Engineering Block S2.1, 50 Nanyang Avenue, 639798, Singapore, Singapore
| |
Collapse
|
13
|
Tichý O, Pederzoli M, Pittner J, Burda JV. Vertical Excitation Energies and Lifetimes of the Two Lowest Singlet Excited States of Cytosine, 5-Aza-cytosine, and the Triazine Family: Quantum Mechanics-Molecular Mechanics Studies. J Chem Theory Comput 2023; 19:1976-1985. [PMID: 36961980 PMCID: PMC10100535 DOI: 10.1021/acs.jctc.2c01262] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/26/2023]
Abstract
A swarm of semi-classical quantum mechanics/molecular mechanics molecular-dynamics simulations where OM2/MNDO is combined with the Gromacs program for consideration of explicit water is performed, solving the time-dependent Schrödinger equation in each step of the trajectories together with the Tully's fewest switches algorithm. Within this stochastic treatment, time dependent probabilities of the three lowest electronic states are determined. The fact that nucleobases are quickly deactivated is confirmed in the cytosine case where our best lifetime estimation is τ1=0.82 ps for the model with 100 water molecules with the SPCE force field and a time step of 0.1 fs. Lifetimes of the remaining molecules are visibly longer: 5-azacytosine, 2,4-diamino-1,3,5-triazine (DT), and 2,4,6-triamino-1,3,5-triazine (TT) molecules have an S1 → S0 de-excitation time of slightly above 10 ps. The lifetimes of the triazine family increases with the increasing number of exocyclic amino groups, that is, s-triazine < 2-amino-1,3,5-triazine < DT < TT. This can be explained by a higher mobility of the carbon-bonded hydrogen atoms in comparison with heavier amino groups since their movement is slowed down due to a substantially higher mass than hydrogen atoms, which can easier reach the out-of-plane positions required in the conical intersection structures. Moreover, bulkier NH2 ligands suffer due to greater friction caused by the surrounding water environment. These mechanical aspects caused a change in the explored lifetime dependences in comparison with our previous gas-phase study.
Collapse
Affiliation(s)
- Ondřej Tichý
- Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, 121 16 Prague 2, Czech Republic
| | - Marek Pederzoli
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences, Dolejškova 3, 182 23 Prague 8, Czech Republic
| | - Jiří Pittner
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences, Dolejškova 3, 182 23 Prague 8, Czech Republic
| | - Jaroslav V Burda
- Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, 121 16 Prague 2, Czech Republic
| |
Collapse
|
14
|
Fatková K, Cajzl R, Burda JV. The vertical excitation energies and a lifetime of the two lowest singlet excited states of the conjugated polyenes from C2 to C22: Ab initio, DFT, and semiclassical MNDO-MD simulations. J Comput Chem 2023; 44:777-787. [PMID: 36444915 DOI: 10.1002/jcc.27040] [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: 07/12/2022] [Revised: 10/04/2022] [Accepted: 10/17/2022] [Indexed: 12/03/2022]
Abstract
Electronic excited states in the series of polyene molecules were explored. Optimal ground-state geometry was used for the evaluation of vertical excitation energies. Results of a chosen set of functionals were compared to post-HF methods (EOM-CCSD, NEVPT2, CASPT2, and MRCI). In addition, the semiempirical OM2/MNDO method using MRCISD computational level was confronted with the above-mentioned techniques. Despite the fact that the first excited state has a significant double-excitation character some functionals were able to qualitatively determine the correct state order (where the lowest excited state has a A g - character). The most successful functionals in transition energies predictions were PBE, TPSS and BLYP in Tamm-Dancoff approach (TDA), which had the smallest root-mean-square deviation (RMSD) scoring towards the experimental values. Regarding RMSD scoring, the OM2/MNDO method performed fairly well, too. Besides absorption spectra, lifetimes of the first two excited states were estimated based on a stochastic approach exploring a swarm of OM2/MNDO hopping dynamics using the Tully fewest switch algorithm for each molecule. The longest lifetime of the first excited state (S1 ) was found for decapentaene (about 5 ps). Further elongation of the conjugated chain caused a mild decrease of this value to ca 1.5 ps for docosaundecaene.
Collapse
Affiliation(s)
- Kateřina Fatková
- Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Prague, Czech Republic
| | - Radim Cajzl
- Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Prague, Czech Republic
| | - Jaroslav V Burda
- Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Prague, Czech Republic
| |
Collapse
|
15
|
Oyeniyi E, Akin-Ojo O, Popoola OO. oeINDO: Efficient determination of excitation energies and UV–Vis absorption spectra of nano-sized Zn, Cd, S and their complexes. COMPUT THEOR CHEM 2023. [DOI: 10.1016/j.comptc.2023.114096] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/13/2023]
|
16
|
Bosia F, Zheng P, Vaucher A, Weymuth T, Dral PO, Reiher M. Ultra-fast semi-empirical quantum chemistry for high-throughput computational campaigns with Sparrow. J Chem Phys 2023; 158:054118. [PMID: 36754821 DOI: 10.1063/5.0136404] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023] Open
Abstract
Semi-empirical quantum chemical approaches are known to compromise accuracy for the feasibility of calculations on huge molecules. However, the need for ultrafast calculations in interactive quantum mechanical studies, high-throughput virtual screening, and data-driven machine learning has shifted the emphasis toward calculation runtimes recently. This comes with new constraints for the software implementation as many fast calculations would suffer from a large overhead of the manual setup and other procedures that are comparatively fast when studying a single molecular structure, but which become prohibitively slow for high-throughput demands. In this work, we discuss the effect of various well-established semi-empirical approximations on calculation speed and relate this to data transfer rates from the raw-data source computer to the results of the visualization front end. For the former, we consider desktop computers, local high performance computing, and remote cloud services in order to elucidate the effect on interactive calculations, for web and cloud interfaces in local applications, and in world-wide interactive virtual sessions. The models discussed in this work have been implemented into our open-source software SCINE Sparrow.
Collapse
Affiliation(s)
- Francesco Bosia
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Peikun Zheng
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, Department of Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| | - Alain Vaucher
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Thomas Weymuth
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Pavlo O Dral
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, Department of Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| | - Markus Reiher
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| |
Collapse
|
17
|
Fang Y, Huang H, Lin K, Xu C, Gu FL, Lan Z. The impact of different geometrical restrictions on the nonadiabatic photoisomerization of biliverdin chromophores. Phys Chem Chem Phys 2022; 24:26190-26199. [PMID: 36278817 DOI: 10.1039/d2cp02941c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
The photoisomerization mechanism of the chromophore of bacterial biliverdin (BV) phytochromes is explored via nonadiabatic dynamics simulation by using the on-the-fly trajectory surface-hopping method at the semi-empirical OM2/MRCI level. Particularly, the current study focuses on the influence of geometrical constrains on the nonadiabatic photoisomerization dynamics of the BV chromophore. Here a rather simplified approach is employed in the nonadiabatic dynamics to capture the features of geometrical constrains, which adds mechanical restrictions to the specific moieties of the BV chromophore. This simplified method provides a rather quick approach to examine the influence of geometrical restrictions on photoisomerization. As expected, different constrains bring distinctive influences on the photoisomerization mechanism of the BV chromophore, giving either strong or minor modification of both involved reaction channels and excited-state lifetimes after the constrains are added in different ring moieties. These observations not only contribute to the primary understanding of the role of the spatial restriction caused by biological environments in photoinduced dynamics of the BV chromophore, but also provide useful ideas for the artificial regulation of the photoisomerization reaction channels of phytochrome proteins.
Collapse
Affiliation(s)
- Yuan Fang
- School of Chemistry, South China Normal University, Guangzhou 510006, China
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, P. R. China.
| | - Haiyi Huang
- School of Chemistry, South China Normal University, Guangzhou 510006, China
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, P. R. China.
| | - Kunni Lin
- School of Chemistry, South China Normal University, Guangzhou 510006, China
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, P. R. China.
| | - Chao Xu
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, P. R. China.
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| | - Feng Long Gu
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, P. R. China.
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| | - Zhenggang Lan
- MOE Key Laboratory of Environmental Theoretical Chemistry, South China Normal University, Guangzhou 510006, P. R. China.
- SCNU Environmental Research Institute, Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| |
Collapse
|
18
|
Pracht P, Bannwarth C. Fast Screening of Minimum Energy Crossing Points with Semiempirical Tight-Binding Methods. J Chem Theory Comput 2022; 18:6370-6385. [PMID: 36121838 DOI: 10.1021/acs.jctc.2c00578] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The investigation of photochemical processes is a highly active field in computational chemistry. One research direction is the automated exploration and identification of minimum energy conical intersection (MECI) geometries. However, due to the immense technical effort required to calculate nonadiabatic potential energy landscapes, the routine application of such computational protocols is severely limited. In this study, we will discuss the prospect of combining adiabatic potential energy surfaces from semiempirical quantum mechanical calculations with specialized confinement potential and metadynamics simulations to identify S0/T1 minimum energy crossing point (MECP) geometries. It is shown that MECPs calculated at the GFN2-xTB level can provide suitable approximations to high-level S0/S1ab initio conical intersection geometries at a fraction of the computational cost. Reference MECIs of benzene are studied to illustrate the basic concept. An example application of the presented protocol is demonstrated for a set of photoswitch molecules.
Collapse
Affiliation(s)
- Philipp Pracht
- Institute of Physical Chemistry, RWTH Aachen University, Melatener Str. 20, 52056Aachen, Germany
| | - Christoph Bannwarth
- Institute of Physical Chemistry, RWTH Aachen University, Melatener Str. 20, 52056Aachen, Germany
| |
Collapse
|
19
|
Wang S, Kind T, Bremer PL, Tantillo DJ, Fiehn O. Beyond the Ground State: Predicting Electron Ionization Mass Spectra Using Excited-State Molecular Dynamics. J Chem Inf Model 2022; 62:4403-4410. [DOI: 10.1021/acs.jcim.2c00597] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Shunyang Wang
- West Coast Metabolomics Center, UC Davis Genome Center, University of California, 451 Health Sciences Drive, Davis, California 95616, United States
- Department of Chemistry, University of California, 1 Shields Avenue, Davis, California 95616, United States
| | - Tobias Kind
- West Coast Metabolomics Center, UC Davis Genome Center, University of California, 451 Health Sciences Drive, Davis, California 95616, United States
| | - Parker Ladd Bremer
- West Coast Metabolomics Center, UC Davis Genome Center, University of California, 451 Health Sciences Drive, Davis, California 95616, United States
- Department of Chemistry, University of California, 1 Shields Avenue, Davis, California 95616, United States
| | - Dean J. Tantillo
- Department of Chemistry, University of California, 1 Shields Avenue, Davis, California 95616, United States
| | - Oliver Fiehn
- West Coast Metabolomics Center, UC Davis Genome Center, University of California, 451 Health Sciences Drive, Davis, California 95616, United States
| |
Collapse
|
20
|
Sun Q, Gieseking RLM. Parametrization of the PM7 Semiempirical Quantum Mechanical Method for Silver Nanoclusters. J Phys Chem A 2022; 126:6558-6569. [PMID: 36082665 DOI: 10.1021/acs.jpca.2c05782] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Semiempirical quantum mechanical methods (SEQMs) are widely used in computational chemistry because of their low computational cost, but their accuracy depends on the quality of the parameters. The neglect of diatomic differential overlap method PM7 is among the few SEQMs that contain parameters for Ag, but the experimental reference data was insufficient to obtain reliable parameters in the original parametrization. In this work, we reparametrize the PM7 parameters for Ag to accurately reproduce the ground-state potential energy surfaces of Ag clusters. Since little experimental data is available, we use reference data obtained from the ab initio method CCSD(T). The resulting parameters significantly reduce the errors in binding energies, energies required to displace clusters along their normal modes, and relative energies of isomers compared to the default PM7 Ag parameters.
Collapse
Affiliation(s)
- Qiwei Sun
- Department of Chemistry, Brandeis University, 415 South Street, Waltham, Massachusetts 02453, United States
| | - Rebecca L M Gieseking
- Department of Chemistry, Brandeis University, 415 South Street, Waltham, Massachusetts 02453, United States
| |
Collapse
|
21
|
Mocci F, de Villiers Engelbrecht L, Olla C, Cappai A, Casula MF, Melis C, Stagi L, Laaksonen A, Carbonaro CM. Carbon Nanodots from an In Silico Perspective. Chem Rev 2022; 122:13709-13799. [PMID: 35948072 PMCID: PMC9413235 DOI: 10.1021/acs.chemrev.1c00864] [Citation(s) in RCA: 17] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Abstract
Carbon nanodots (CNDs) are the latest and most shining rising stars among photoluminescent (PL) nanomaterials. These carbon-based surface-passivated nanostructures compete with other related PL materials, including traditional semiconductor quantum dots and organic dyes, with a long list of benefits and emerging applications. Advantages of CNDs include tunable inherent optical properties and high photostability, rich possibilities for surface functionalization and doping, dispersibility, low toxicity, and viable synthesis (top-down and bottom-up) from organic materials. CNDs can be applied to biomedicine including imaging and sensing, drug-delivery, photodynamic therapy, photocatalysis but also to energy harvesting in solar cells and as LEDs. More applications are reported continuously, making this already a research field of its own. Understanding of the properties of CNDs requires one to go to the levels of electrons, atoms, molecules, and nanostructures at different scales using modern molecular modeling and to correlate it tightly with experiments. This review highlights different in silico techniques and studies, from quantum chemistry to the mesoscale, with particular reference to carbon nanodots, carbonaceous nanoparticles whose structural and photophysical properties are not fully elucidated. The role of experimental investigation is also presented. Hereby, we hope to encourage the reader to investigate CNDs and to apply virtual chemistry to obtain further insights needed to customize these amazing systems for novel prospective applications.
Collapse
Affiliation(s)
- Francesca Mocci
- Department
of Chemical and Geological Sciences, University
of Cagliari, I-09042 Monserrato, Italy,
| | | | - Chiara Olla
- Department
of Physics, University of Cagliari, I-09042 Monserrato, Italy
| | - Antonio Cappai
- Department
of Physics, University of Cagliari, I-09042 Monserrato, Italy
| | - Maria Francesca Casula
- Department
of Mechanical, Chemical and Materials Engineering, University of Cagliari, Via Marengo 2, IT 09123 Cagliari, Italy
| | - Claudio Melis
- Department
of Physics, University of Cagliari, I-09042 Monserrato, Italy
| | - Luigi Stagi
- Department
of Chemistry and Pharmacy, Laboratory of Materials Science and Nanotechnology, University of Sassari, Via Vienna 2, 07100 Sassari, Italy
| | - Aatto Laaksonen
- Department
of Chemical and Geological Sciences, University
of Cagliari, I-09042 Monserrato, Italy,Department
of Materials and Environmental Chemistry, Arrhenius Laboratory, Stockholm University, SE-106 91 Stockholm, Sweden,State Key
Laboratory of Materials-Oriented and Chemical Engineering, Nanjing Tech University, Nanjing 210009, P. R. China,Centre
of Advanced Research in Bionanoconjugates and Biopolymers, PetruPoni Institute of Macromolecular Chemistry, Aleea Grigore Ghica-Voda 41A, 700487 Iasi, Romania,Division
of Energy Science, Energy Engineering, Luleå
University of Technology, Luleå 97187, Sweden,
| | | |
Collapse
|
22
|
Different timescales during ultrafast stilbene isomerization in the gas and liquid phases revealed using time-resolved photoelectron spectroscopy. Nat Chem 2022; 14:1126-1132. [PMID: 35953643 PMCID: PMC7613649 DOI: 10.1038/s41557-022-01012-0] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2020] [Accepted: 06/28/2022] [Indexed: 11/08/2022]
Abstract
Directly contrasting ultrafast excited-state dynamics in the gas and liquid phases is crucial to understanding the influence of complex environments. Previous studies have often relied on different spectroscopic observables, rendering direct comparisons challenging. Here, we apply extreme-ultraviolet time-resolved photoelectron spectroscopy to both gaseous and liquid cis-stilbene, revealing the coupled electronic and nuclear dynamics that underlie its isomerization. Our measurements track the excited-state wave packets from excitation along the complete reaction path to the final products. We observe coherent excited-state vibrational dynamics in both phases of matter that persist to the final products, enabling the characterization of the branching space of the S1-S0 conical intersection. We observe a systematic lengthening of the relaxation timescales in the liquid phase and a red shift of the measured excited-state frequencies that is most pronounced for the complex reaction coordinate. These results characterize in detail the influence of the liquid environment on both electronic and structural dynamics during a complete photochemical transformation.
Collapse
|
23
|
Menezes F, Popowicz GM. ULYSSES: An Efficient and Easy to Use Semiempirical Library for C+. J Chem Inf Model 2022; 62:3685-3694. [PMID: 35930308 DOI: 10.1021/acs.jcim.2c00757] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We introduce ULYSSES, a user-friendly and robust C++ library for semiempirical quantum chemical calculations. In its current version, ULYSSES is equipped with a large set of different semiempirical models, most of which are based on the Neglect of Diatomic Differential Overlap (NDDO) approximation. Empirical corrections for dispersion and hydrogen bonding are available for most methods, so that higher quality is achieved in the calculation of energies of nonbonded complexes. The library is furthermore equipped with geometry optimization, as well as modules for calculating molecular properties of general interest. Ideal gas thermodynamics is available and allows single structure as well as conformer (multistructure) averaged properties to be calculated. We offer the possibility to use several vibrational partition functions according to the nature of interactions being studied: for covalent systems, the traditional harmonic oscillator approximation is available; for nonbonded complexes, we systematically extended the partition function proposed by Grimme for all thermodynamic functions. The library is also capable of running Born-Oppenheimer molecular dynamics.
Collapse
Affiliation(s)
- Filipe Menezes
- Institute of Structural Biology, Helmholtz Zentrum Muenchen, Ingolstaedter Landstrasse 1, 85764 Neuherberg, Germany
| | - Grzegorz M Popowicz
- Institute of Structural Biology, Helmholtz Zentrum Muenchen, Ingolstaedter Landstrasse 1, 85764 Neuherberg, Germany
| |
Collapse
|
24
|
Deep learning of dynamically responsive chemical Hamiltonians with semiempirical quantum mechanics. Proc Natl Acad Sci U S A 2022; 119:e2120333119. [PMID: 35776544 PMCID: PMC9271210 DOI: 10.1073/pnas.2120333119] [Citation(s) in RCA: 19] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022] Open
Abstract
Machine learning is revolutionizing computational chemistry by greatly reducing the computational difficulty of many simulations performed by computational chemists while maintaining accuracies of 1 kcal/mol or better. A major challenge in this field is addressing the poor extensibility and transferability of conventional machine-learning (ML) models, which result in degraded accuracy when applying these models to large or new chemical systems. To build a more general and interpretable model, we incorporate a quantum chemistry framework into the deep neural network, resulting in an interpretable Hamiltonian-based model with markedly high training efficiency. We validate this method on multiple large biochemical molecules by predicting various properties with consistently high accuracies, indicating the model is both extensible and transferable. Conventional machine-learning (ML) models in computational chemistry learn to directly predict molecular properties using quantum chemistry only for reference data. While these heuristic ML methods show quantum-level accuracy with speeds several orders of magnitude faster than traditional quantum chemistry methods, they suffer from poor extensibility and transferability; i.e., their accuracy degrades on large or new chemical systems. Incorporating quantum chemistry frameworks into the ML models directly solves this problem. Here we take the structure of semiempirical quantum mechanics (SEQM) methods to construct dynamically responsive Hamiltonians. SEQM methods use empirical parameters fitted to experimental properties to construct reduced-order Hamiltonians, facilitating much faster calculations than ab initio methods but with compromised accuracy. By replacing these static parameters with machine-learned dynamic values inferred from the local environment, we greatly improve the accuracy of the SEQM methods. Trained on molecular energies and atomic forces, these dynamically generated Hamiltonian parameters show a strong correlation with atomic hybridization and bonding. Trained with only about 60,000 small organic molecular conformers, the resulting model retains interpretability, extensibility, and transferability when testing on much larger chemical systems and predicting various molecular properties. Overall, this work demonstrates the virtues of incorporating physics-based descriptions with ML to develop models that are simultaneously accurate, transferable, and interpretable.
Collapse
|
25
|
Wang Y, Seritan S, Lahana D, Ford JE, Valentini A, Hohenstein EG, Martínez TJ. InteraChem: Exploring Excited States in Virtual Reality with Ab Initio Interactive Molecular Dynamics. J Chem Theory Comput 2022; 18:3308-3317. [PMID: 35649124 DOI: 10.1021/acs.jctc.2c00005] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
InteraChem is an ab initio interactive molecular dynamics (AI-IMD) visualizer that leverages recent advances in virtual reality hardware and software, as well as the graphical processing unit (GPU)-accelerated TeraChem electronic structure package, in order to render quantum chemistry in real time. We introduce the exploration of electronically excited states via AI-IMD using the floating occupation molecular orbital-complete active space configuration interaction method. The optimization tools in InteraChem enable identification of excited state minima as well as minimum energy conical intersections for further characterization of excited state chemistry in small- to medium-sized systems. We demonstrate that finite-temperature Hartree-Fock theory is an efficient method to perform ground state AI-IMD. InteraChem allows users to track electronic properties such as molecular orbitals and bond order in real time, resulting in an interactive visualization tool that aids in the interpretation of excited state chemistry data and makes quantum chemistry more accessible for both research and educational purposes.
Collapse
Affiliation(s)
- Yuanheng Wang
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Stefan Seritan
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Dean Lahana
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Jason E Ford
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Alessio Valentini
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Edward G Hohenstein
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Todd J Martínez
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| |
Collapse
|
26
|
Salazar E, Reinink S, Faraji S. Providing theoretical insight into the role of symmetry in the photoisomerization mechanism of a non-symmetric dithienylethene photoswitch. Phys Chem Chem Phys 2022; 24:11592-11602. [PMID: 35531648 PMCID: PMC9116444 DOI: 10.1039/d2cp00550f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Dithienylethene (DTE) molecular photoswitches have shown to be excellent candidates in the design of efficient optoelectronic devices, due to their high photoisomerization quantum yield (QY), for which symmetry is suggested to play a crucial role. Here, we present a theoretical study on the photochemistry of a non-symmetric dithienylethene photoswitch, with a special emphasis on the effect of asymmetric substitution on the photocyclization and photoreversion mechanisms. We used the Spin-Flip Time Dependent Density Functional Theory (SF-TDDFT) method to locate and characterize the main structures (conical intersections and minima) of the ground state and the first two excited states, S1 and S2, along the ring-opening/closure reaction coordinate of the photocyclization and photoreversion processes, and to identify the important coordinates governing the radiationless decay pathways. Our results suggest that while the main features that characterize the photoisomerization of symmetric DTEs are also present for the photoisomerization of the non-symmetric DTE, the lower energy barrier on S1 along the cycloreversion reaction speaks in favor of a more efficient and therefore a higher cycloreversion QY for the non-symmetric DTEs, making them a better candidate for molecular optoelectronic devices than their symmetric counterparts.
Collapse
Affiliation(s)
- Edison Salazar
- Theoretical Chemistry, Zernike Institute for Advanced Materials, University of GroningenNijenborgh 49747 AG GroningenThe Netherlands
| | - Suzanne Reinink
- Theoretical Chemistry, Zernike Institute for Advanced Materials, University of GroningenNijenborgh 49747 AG GroningenThe Netherlands
| | - Shirin Faraji
- Theoretical Chemistry, Zernike Institute for Advanced Materials, University of GroningenNijenborgh 49747 AG GroningenThe Netherlands
| |
Collapse
|
27
|
Daub CD, Zakai I, Valiev R, Salo VT, Gerber RB, Kurtén T. Energy transfer, pre-reactive complex formation and recombination reactions during the collision of peroxy radicals. Phys Chem Chem Phys 2022; 24:10033-10043. [PMID: 35415732 DOI: 10.1039/d1cp04720e] [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/15/2022]
Abstract
In this paper we study collisions between polyatomic radicals - an important process in fields ranging from biology to combustion. Energy transfer, formation of intermediate complexes and recombination reactions are treated, with applications to peroxy radicals in atmospheric chemistry. Multi-reference perturbation theory, supplemented by coupled-cluster calculations, describes the potential energy surfaces with high accuracy, including the interaction of singlet and triplet spin states during radical recombination. Our multi-reference molecular dynamics (MD) trajectories on methyl peroxy radicals confirm the reaction mechanism postulated in earlier studies. Specifically, they show that if suitable pre-reactive complexes are formed, they will rapidly lead to the formation and subsequent decomposition of tetroxide intermediates. However, generating multi-reference MD trajectories is exceedingly computationally demanding, and we cannot adequately sample the whole conformational space. To answer this challenge, we promote the use of a novel simplified semi-empirical MD methodology. It assumes the collision is governed by two states, a singlet (S0) and a triplet (T1) state. The method predicts differences between collisions on S0 and T1 surfaces, and qualitatively includes not only pre-reactive complex formation, but also recombination processes such as tetroxide formation. Finally, classical MD simulations using force-fields for non-reactive collisions are employed to generate thousands of collision trajectories, to verify that the semi-empirical method is sampling collisions adequately, and to carry out preliminary investigations of larger systems. For systems with low activation energies, the experimental rate coefficient is surprisingly well reproduced by simply multiplying the gas-kinetic collision rate by the simulated probability for long-lived complex formation.
Collapse
Affiliation(s)
| | - Itai Zakai
- Department of Chemistry, Hebrew University of Jerusalem, Jerusalem, Israel.
| | - Rashid Valiev
- Department of Chemistry, University of Helsinki, P.O. Box 55, Helsinki 00014, Finland.
| | - Vili-Taneli Salo
- Department of Chemistry, University of Helsinki, P.O. Box 55, Helsinki 00014, Finland.
| | - R Benny Gerber
- Department of Chemistry, Hebrew University of Jerusalem, Jerusalem, Israel. .,Department of Chemistry, University of California Irvine, Irvine, CA 92697, USA
| | - Theo Kurtén
- Department of Chemistry, University of Helsinki, P.O. Box 55, Helsinki 00014, Finland.
| |
Collapse
|
28
|
Zheng P, Yang W, Wu W, Isayev O, Dral PO. Toward Chemical Accuracy in Predicting Enthalpies of Formation with General-Purpose Data-Driven Methods. J Phys Chem Lett 2022; 13:3479-3491. [PMID: 35416675 DOI: 10.1021/acs.jpclett.2c00734] [Citation(s) in RCA: 14] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Enthalpies of formation and reaction are important thermodynamic properties that have a crucial impact on the outcome of chemical transformations. Here we implement the calculation of enthalpies of formation with a general-purpose ANI-1ccx neural network atomistic potential. We demonstrate on a wide range of benchmark sets that both ANI-1ccx and our other general-purpose data-driven method AIQM1 approach the coveted chemical accuracy of 1 kcal/mol with the speed of semiempirical quantum mechanical methods (AIQM1) or faster (ANI-1ccx). It is remarkably achieved without specifically training the machine learning parts of ANI-1ccx or AIQM1 on formation enthalpies. Importantly, we show that these data-driven methods provide statistical means for uncertainty quantification of their predictions, which we use to detect and eliminate outliers and revise reference experimental data. Uncertainty quantification may also help in the systematic improvement of such data-driven methods.
Collapse
Affiliation(s)
- Peikun Zheng
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, Department of Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| | - Wudi Yang
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, Department of Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| | - Wei Wu
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, Department of Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| | - Olexandr Isayev
- Department of Chemistry, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States
| | - Pavlo O Dral
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, Department of Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| |
Collapse
|
29
|
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
|
30
|
Ríos-García M, Fernández B, Rodríguez-Otero J, Cabaleiro-Lago EM, Vázquez SA. The PM6-FGC Method: Improved Corrections for Amines and Amides. Molecules 2022; 27:molecules27051678. [PMID: 35268779 PMCID: PMC8924896 DOI: 10.3390/molecules27051678] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2022] [Revised: 02/23/2022] [Accepted: 03/01/2022] [Indexed: 12/10/2022] Open
Abstract
Recently, we reported a new approach to develop pairwise analytical corrections to improve the description of noncovalent interactions, by approximate methods of electronic structures, such as semiempirical quantum mechanical (SQM) methods. In particular, and as a proof of concept, we used the PM6 Hamiltonian and we named the method PM6-FGC, where the FGC acronym, corresponding to Functional Group Corrections, emphasizes the idea that the corrections work for specific functional groups rather than for individual atom pairs. The analytical corrections were derived from fits to B3LYP-D3/def2-TZVP (reference). PM6 interaction energy differences, evaluated for a reduced set of small bimolecular complexes, were chosen as representatives of saturated hydrocarbons, carboxylic, amine and, tentatively, amide functional groups. For the validation, the method was applied to several complexes of well-known databases, as well as to complexes of diglycine and dialanine, assuming the transferability of amine group corrections to amide groups. The PM6-FGC method showed great potential but revealed significant inaccuracies for the description of some interactions involving the –NH2 group in amines and amides, caused by the inadequate selection of the model compound used to represent these functional groups (an NH3 molecule). In this work, methylamine and acetamide are used as representatives of amine and amide groups, respectively. This new selection leads to significant improvements in the calculation of noncovalent interactions in the validation set.
Collapse
|
31
|
Tichý O, Burda JV. Estimation of electron absorption spectra and lifetime of the two lowest singlet excited states of pyrimidine nucleobases and their derivatives. J Mol Struct 2022. [DOI: 10.1016/j.molstruc.2021.131863] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
32
|
Zheng P, Zubatyuk R, Wu W, Isayev O, Dral PO. Artificial intelligence-enhanced quantum chemical method with broad applicability. Nat Commun 2021; 12:7022. [PMID: 34857738 PMCID: PMC8640006 DOI: 10.1038/s41467-021-27340-2] [Citation(s) in RCA: 34] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2021] [Accepted: 11/10/2021] [Indexed: 11/08/2022] Open
Abstract
High-level quantum mechanical (QM) calculations are indispensable for accurate explanation of natural phenomena on the atomistic level. Their staggering computational cost, however, poses great limitations, which luckily can be lifted to a great extent by exploiting advances in artificial intelligence (AI). Here we introduce the general-purpose, highly transferable artificial intelligence-quantum mechanical method 1 (AIQM1). It approaches the accuracy of the gold-standard coupled cluster QM method with high computational speed of the approximate low-level semiempirical QM methods for the neutral, closed-shell species in the ground state. AIQM1 can provide accurate ground-state energies for diverse organic compounds as well as geometries for even challenging systems such as large conjugated compounds (fullerene C60) close to experiment. This opens an opportunity to investigate chemical compounds with previously unattainable speed and accuracy as we demonstrate by determining geometries of polyyne molecules-the task difficult for both experiment and theory. Noteworthy, our method's accuracy is also good for ions and excited-state properties, although the neural network part of AIQM1 was never fitted to these properties.
Collapse
Affiliation(s)
- Peikun Zheng
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, Department of Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, 361005, China
| | - Roman Zubatyuk
- Department of Chemistry, Carnegie Mellon University, Pittsburgh, PA, 15213, USA
| | - Wei Wu
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, Department of Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, 361005, China
| | - Olexandr Isayev
- Department of Chemistry, Carnegie Mellon University, Pittsburgh, PA, 15213, USA.
| | - Pavlo O Dral
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, Department of Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, 361005, China.
| |
Collapse
|
33
|
Rational prioritization strategy allows the design of macrolide derivatives that overcome antibiotic resistance. Proc Natl Acad Sci U S A 2021; 118:2113632118. [PMID: 34750269 PMCID: PMC8609559 DOI: 10.1073/pnas.2113632118] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 09/24/2021] [Indexed: 11/21/2022] Open
Abstract
Due to the development of resistance against commonly used antibiotics, new derivatives that avoid resistance mechanisms need to be developed. To address this problem, a rational prioritization strategy is outlined for macrolide antibiotics. Candidates are screened based on their solubility, membrane permeability, and binding affinity using a tiered optimization approach of free energy simulations and quantum mechanics/molecular mechanics calculations. After prioritization by computational methods, the best candidates are evaluated experimentally. The strategy creates a targeted substance library that is highly enriched in compounds with antibacterial activity. This allows for faster iterations in the development of new antibiotic derivatives. Antibiotic resistance is a major threat to global health; this problem can be addressed by the development of new antibacterial agents to keep pace with the evolutionary adaptation of pathogens. Computational approaches are essential tools to this end since their application enables fast and early strategical decisions in the drug development process. We present a rational design approach, in which acylide antibiotics were screened based on computational predictions of solubility, membrane permeability, and binding affinity toward the ribosome. To assess our design strategy, we tested all candidates for in vitro inhibitory activity and then evaluated them in vivo with several antibiotic-resistant strains to determine minimal inhibitory concentrations. The predicted best candidate is synthetically more accessible, exhibits higher solubility and binding affinity to the ribosome, and is up to 56 times more active against resistant pathogens than telithromycin. Notably, the best compounds designed by us show activity, especially when combined with the membrane-weakening drug colistin, against Acinetobacter baumanii, Pseudomonas aeruginosa, and Escherichia coli, which are the three most critical targets from the priority list of pathogens of the World Health Organization.
Collapse
|
34
|
Hohenstein EG, Yu JK, Bannwarth C, List NH, Paul AC, Folkestad SD, Koch H, Martínez TJ. Predictions of Pre-edge Features in Time-Resolved Near-Edge X-ray Absorption Fine Structure Spectroscopy from Hole-Hole Tamm-Dancoff-Approximated Density Functional Theory. J Chem Theory Comput 2021; 17:7120-7133. [PMID: 34623139 DOI: 10.1021/acs.jctc.1c00478] [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/16/2022]
Abstract
Time-resolved near-edge X-ray absorption fine structure (TR-NEXAFS) spectroscopy is a powerful technique for studying photochemical reaction dynamics with femtosecond time resolution. In order to avoid ambiguity in TR-NEXAFS spectra from nonadiabatic dynamics simulations, core- and valence-excited states must be evaluated on equal footing and those valence states must also define the potential energy surfaces used in the nonadiabatic dynamics simulation. In this work, we demonstrate that hole-hole Tamm-Dancoff-approximated density functional theory (hh-TDA) is capable of directly simulating TR-NEXAFS spectroscopies. We apply hh-TDA to the excited-state dynamics of acrolein. We identify two pre-edge features in the oxygen K-edge TR-NEXAFS spectrum associated with the S2 (ππ*) and S1 (nπ*) excited states. We show that these features can be used to follow the internal conversion dynamics between the lowest three electronic states of acrolein. Due to the low, O(N2) apparent computational complexity of hh-TDA and our GPU-accelerated implementation, this method is promising for the simulation of pre-edge features in TR-NEXAFS spectra of large molecules and molecules in the condensed phase.
Collapse
Affiliation(s)
- Edward G Hohenstein
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Jimmy K Yu
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States.,Biophysics Program, Stanford University, Stanford, California 94305, United States
| | - Christoph Bannwarth
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Nanna Holmgaard List
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Alexander C Paul
- Department of Chemistry, Norwegian University of Science and Technology, NTNU, 7491 Trondheim, Norway
| | - Sarai D Folkestad
- Department of Chemistry, Norwegian University of Science and Technology, NTNU, 7491 Trondheim, Norway
| | - Henrik Koch
- Department of Chemistry, Norwegian University of Science and Technology, NTNU, 7491 Trondheim, Norway.,Scuola Normale Superiore, Piazza dei Cavaleri 7, 56126 Pisa, Italy
| | - Todd J Martínez
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| |
Collapse
|
35
|
Liu J, Lan Z, Yang J. An efficient implementation of spin-orbit coupling within the framework of semiempirical orthogonalization-corrected methods for ultrafast intersystem crossing dynamics. Phys Chem Chem Phys 2021; 23:22313-22323. [PMID: 34591049 DOI: 10.1039/d1cp03477d] [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/21/2022]
Abstract
We implement spin-orbit coupling (SOC) within the framework of semiempirical orthogonalization-corrected methods (OMx). The excited-state wavefunction is generated from configuration interaction with single excitations (CIS). The SOC Hamiltonian in terms of the one-electron Breit-Pauli operator with effective nuclear charges is adopted in this work. Benchmark calculations show that SOCs evaluated using the OMx/CIS method agree very well with those obtained from time-dependent density functional theory. As a particularly attractive application, we incorporate SOCs between singlet and triplet states into Tully's fewest switches surface hopping algorithm to enable excited-state nonadiabatic dynamics simulations, treating internal conversion and intersystem crossing on an equal footing. This semiempirical dynamics simulation approach is applied to investigate ultrafast intersystem crossing processes in core-substituted naphthalenediimides.
Collapse
Affiliation(s)
- Jie Liu
- Hefei National Laboratory for Physical Sciences at Microscale, University of Science and Technology of China, Hefei, Anhui 230026, China.
| | - Zhenggang Lan
- Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety and MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| | - Jinlong Yang
- Hefei National Laboratory for Physical Sciences at Microscale, University of Science and Technology of China, Hefei, Anhui 230026, China.
| |
Collapse
|
36
|
Pérez-Tabero S, Fernández B, Cabaleiro-Lago EM, Martínez-Núñez E, Vázquez SA. New Approach for Correcting Noncovalent Interactions in Semiempirical Quantum Mechanical Methods: The Importance of Multiple-Orientation Sampling. J Chem Theory Comput 2021; 17:5556-5567. [PMID: 34424696 PMCID: PMC8486165 DOI: 10.1021/acs.jctc.1c00365] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
![]()
A new
approach is presented to improve the performance of semiempirical
quantum mechanical (SQM) methods in the description of noncovalent
interactions. To show the strategy, the PM6 Hamiltonian was selected,
although, in general, the procedure can be applied to other semiempirical
Hamiltonians and to different methodologies. A set of small molecules
were selected as representative of various functional groups, and
intermolecular potential energy curves (IPECs) were evaluated for
the most relevant orientations of interacting molecular pairs. Then,
analytical corrections to PM6 were derived from fits to B3LYP-D3/def2-TZVP
reference–PM6 interaction energy differences. IPECs provided
by the B3LYP-D3/def2-TZVP combination of the electronic structure
method and basis set were chosen as the reference because they are
in excellent agreement with CCSD(T)/aug-cc-pVTZ curves for the studied
systems. The resulting method, called PM6-FGC (from functional group
corrections), significantly improves the performance of PM6 and shows
the importance of including a sufficient number of orientations of
the interacting molecules in the reference data set in order to obtain
well-balanced descriptions.
Collapse
Affiliation(s)
- Sergio Pérez-Tabero
- Departamento de Química Física, Facultade de Química, Universidade de Santiago de Compostela, Santiago de Compostela 15782, Spain
| | - Berta Fernández
- Departamento de Química Física, Facultade de Química, Universidade de Santiago de Compostela, Santiago de Compostela 15782, Spain
| | - Enrique M Cabaleiro-Lago
- Departamento de Química Física, Facultade de Química, Universidade de Santiago de Compostela, Santiago de Compostela 15782, Spain
| | - Emilio Martínez-Núñez
- Departamento de Química Física, Facultade de Química, Universidade de Santiago de Compostela, Santiago de Compostela 15782, Spain
| | - Saulo A Vázquez
- Departamento de Química Física, Facultade de Química, Universidade de Santiago de Compostela, Santiago de Compostela 15782, Spain
| |
Collapse
|
37
|
Bondanza M, Jacquemin D, Mennucci B. Excited States of Xanthophylls Revisited: Toward the Simulation of Biologically Relevant Systems. J Phys Chem Lett 2021; 12:6604-6612. [PMID: 34251826 PMCID: PMC8311646 DOI: 10.1021/acs.jpclett.1c01929] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2021] [Accepted: 07/08/2021] [Indexed: 06/13/2023]
Abstract
Xanthophylls are a class of oxygen-containing carotenoids, which play a fundamental role in light-harvesting pigment-protein complexes and in many photoresponsive proteins. The complexity of the manifold of the electronic states and the large sensitivity to the environment still prevent a clear and coherent interpretation of their photophysics and photochemistry. In this Letter, we compare cutting-edge ab initio methods (CC3 and DMRG/NEVPT2) with time-dependent DFT and semiempirical CI (SECI) on model keto-carotenoids and show that SECI represents the right compromise between accuracy and computational cost to be applied to real xanthophylls in their biological environment. As an example, we investigate canthaxanthin in the orange carotenoid protein and show that the conical intersections between excited states and excited-ground states are mostly determined by the effective bond length alternation coordinate, which is significantly tuned by the protein through geometrical constraints and electrostatic effects.
Collapse
Affiliation(s)
- Mattia Bondanza
- Dipartimento
di Chimica e Chimica Industriale, University
of Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| | - Denis Jacquemin
- Université
de Nantes, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| | - Benedetta Mennucci
- Dipartimento
di Chimica e Chimica Industriale, University
of Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| |
Collapse
|
38
|
Zubatiuk T, Nebgen B, Lubbers N, Smith JS, Zubatyuk R, Zhou G, Koh C, Barros K, Isayev O, Tretiak S. Machine learned Hückel theory: Interfacing physics and deep neural networks. J Chem Phys 2021; 154:244108. [PMID: 34241371 DOI: 10.1063/5.0052857] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The Hückel Hamiltonian is an incredibly simple tight-binding model known for its ability to capture qualitative physics phenomena arising from electron interactions in molecules and materials. Part of its simplicity arises from using only two types of empirically fit physics-motivated parameters: the first describes the orbital energies on each atom and the second describes electronic interactions and bonding between atoms. By replacing these empirical parameters with machine-learned dynamic values, we vastly increase the accuracy of the extended Hückel model. The dynamic values are generated with a deep neural network, which is trained to reproduce orbital energies and densities derived from density functional theory. The resulting model retains interpretability, while the deep neural network parameterization is smooth and accurate and reproduces insightful features of the original empirical parameterization. Overall, this work shows the promise of utilizing machine learning to formulate simple, accurate, and dynamically parameterized physics models.
Collapse
Affiliation(s)
- Tetiana Zubatiuk
- Department of Chemistry, Mellon College of Science, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA
| | - Benjamin Nebgen
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87544, USA
| | - Nicholas Lubbers
- Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - Justin S Smith
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87544, USA
| | - Roman Zubatyuk
- Department of Chemistry, Mellon College of Science, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA
| | - Guoqing Zhou
- Department of Physics and Astronomy, University of Southern California, Los Angeles, California 90089, USA
| | - Christopher Koh
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87544, USA
| | - Kipton Barros
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87544, USA
| | - Olexandr Isayev
- Department of Chemistry, Mellon College of Science, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA
| | - Sergei Tretiak
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87544, USA
| |
Collapse
|
39
|
Jankowska J, Martyka M, Michalski M. Photo-cycloreversion mechanism in diarylethenes revisited: A multireference quantum-chemical study at the ODM2/MRCI level. J Chem Phys 2021; 154:204305. [PMID: 34241185 DOI: 10.1063/5.0045830] [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/14/2022] Open
Abstract
Photoswitchable diarylethenes (DAEs), over years of intense fundamental and applied research, have been established among the most commonly chosen molecular photoswitches, often employed as controlling units in molecular devices and smart materials. At the same time, providing reliable explanation for their photophysical behavior, especially the mechanism of the photo-cycloreversion transformation, turned out to be a highly challenging task. Herein, we investigate this mechanism in detail by means of multireference semi-empirical quantum chemistry calculations, allowing, for the first time, for a balanced treatment of the static and dynamic correlation effects, both playing a crucial role in DAE photochemistry. In the course of our study, we find the second singlet excited state of double electronic-excitation character to be the key to understanding the nature of the photo-cycloreversion transformation in DAE molecular photoswitches.
Collapse
Affiliation(s)
- J Jankowska
- Faculty of Chemistry, University of Warsaw, Warsaw 02-093, Poland
| | - M Martyka
- Faculty of Chemistry, University of Warsaw, Warsaw 02-093, Poland
| | - M Michalski
- Faculty of Chemistry, University of Warsaw, Warsaw 02-093, Poland
| |
Collapse
|
40
|
Feng M, Gilson MK. Mechanistic analysis of light-driven overcrowded alkene-based molecular motors by multiscale molecular simulations. Phys Chem Chem Phys 2021; 23:8525-8540. [PMID: 33876015 PMCID: PMC8102045 DOI: 10.1039/d0cp06685k] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
We analyze light-driven overcrowded alkene-based molecular motors, an intriguing class of small molecules that have the potential to generate MHz-scale rotation rates. The full rotation process is simulated at multiple scales by combining quantum surface-hopping molecular dynamics (MD) simulations for the photoisomerization step with classical MD simulations for the thermal helix inversion step. A Markov state analysis resolves conformational substates, their interconversion kinetics, and their roles in the motor's rotation process. Furthermore, motor performance metrics, including rotation rate and maximal power output, are computed to validate computations against experimental measurements and to inform future designs. Lastly, we find that to correctly model these motors, the force field must be optimized by fitting selected parameters to reference quantum mechanical energy surfaces. Overall, our simulations yield encouraging agreement with experimental observables such as rotation rates, and provide mechanistic insights that may help future designs.
Collapse
Affiliation(s)
- Mudong Feng
- Department of Chemistry and Biochemistry, University of California, San Diego, 9500 Gilman Drive, La Jolla, 92093, USA.
| | | |
Collapse
|
41
|
Rocha MND, Alves DR, Marinho MM, Morais SMD, Marinho ES. Virtual Screening of Citrus Flavonoid Tangeretin: A Promising Pharmacological Tool for the Treatment and Prevention of Zika fever and COVID-19. JOURNAL OF COMPUTATIONAL BIOPHYSICS AND CHEMISTRY 2021. [DOI: 10.1142/s2737416521500137] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
It is of great importance for the pharmaceutical industry to find therapeutic substances extracted from natural sources, which are abundant, obtained with low costs and presenting the antiviral potential for the treatment of Zika virus (ZIKV) and COVID-19. Tangeretin (TAN) is a citrus polymethoxyflavone from Citrus reticulata peel oil with known antiviral activities, whose physico-chemical properties are not reported. The present study aimed to investigate by a theoretical screening of electronic, structural properties and pharmacodynamic and pharmacokinetic parameters that characterize TAN as a therapeutic drug in the treatment and prevention of zika fever and COVID-19. The molecule reached its minimum energy-forming state of [Formula: see text]795.85747[Formula: see text]kJ/mol and the HOMO and LUMO boundary orbitals reactivity descriptors suggest that the compound is stable and does not tend to be reactive in intermolecular interactions. The ligand connects to the NS1 ZIKV receptor with strong H-bond interactions, also connects with the NS5 ZIKV receptor in a competitive effect with the SAM inhibitor and acts in a supplementary effect with the N3 inhibitor and the BRT drug in the Mpro SARS-CoV-2 receptor. The properties of ADMET shows that the compound suffers few amounts of drug alterations because it inhibits the metabolic enzymes CYP2C9 and CYP3A4 and penetrates the central nervous system, without accumulation of drug residues in the blood or in the lumen in the gastrointestinal tract, without risk of toxicity to the patient. With the results obtained, it is possible to identify TAN as a promising pharmacological tool for the treatment and prevention of Zika fever and COVID-19.
Collapse
Affiliation(s)
- Matheus Nunes da Rocha
- Group of Theoretical Chemistry and Electrochemical, FAFIDAM, Ceará State University, Limoeiro do Norte, Ceará, Brazil
| | - Daniela Ribeiro Alves
- Animal Health Research Center, Ceará State University, Campus Itaperi, Fortaleza, Ceará, Brazil
| | - Marcia Machado Marinho
- Iguatu Faculty of Education, Science and Letters/FECLI, State University of Ceará Iguatu, Ceará, Brazil
| | - Selene Maia de Morais
- Animal Health Research Center, Ceará State University, Campus Itaperi, Fortaleza, Ceará, Brazil
| | - Emmanuel Silva Marinho
- Group of Theoretical Chemistry and Electrochemical, FAFIDAM, Ceará State University, Limoeiro do Norte, Ceará, Brazil
| |
Collapse
|
42
|
Krettler CA, Thallinger GG. A map of mass spectrometry-based in silico fragmentation prediction and compound identification in metabolomics. Brief Bioinform 2021; 22:6184408. [PMID: 33758925 DOI: 10.1093/bib/bbab073] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/26/2020] [Revised: 01/29/2021] [Accepted: 02/12/2021] [Indexed: 12/27/2022] Open
Abstract
Metabolomics, the comprehensive study of the metabolome, and lipidomics-the large-scale study of pathways and networks of cellular lipids-are major driving forces in enabling personalized medicine. Complicated and error-prone data analysis still remains a bottleneck, however, especially for identifying novel metabolites. Comparing experimental mass spectra to curated databases containing reference spectra has been the gold standard for identification of compounds, but constructing such databases is a costly and time-demanding task. Many software applications try to circumvent this process by utilizing cutting-edge advances in computational methods-including quantum chemistry and machine learning-and simulate mass spectra by performing theoretical, so called in silico fragmentations of compounds. Other solutions concentrate directly on experimental spectra and try to identify structural properties by investigating reoccurring patterns and the relationships between them. The considerable progress made in the field allows recent approaches to provide valuable clues to expedite annotation of experimental mass spectra. This review sheds light on individual strengths and weaknesses of these tools, and attempts to evaluate them-especially in view of lipidomics, when considering complex mixtures found in biological samples as well as mass spectrometer inter-instrument variability.
Collapse
Affiliation(s)
- Christoph A Krettler
- Institute of Biomedical Informatics, Graz University of Technology, Stremayrgasse 16/I, 8010, Graz, Austria.,Omics Center Graz, BioTechMed-Graz, Stiftingtalstrasse 24, 8010, Graz, Austria
| | - Gerhard G Thallinger
- Institute of Biomedical Informatics, Graz University of Technology, Stremayrgasse 16/I, 8010, Graz, Austria.,Omics Center Graz, BioTechMed-Graz, Stiftingtalstrasse 24, 8010, Graz, Austria
| |
Collapse
|
43
|
Gieseking RLM. A new release of MOPAC incorporating the INDO/S semiempirical model with CI excited states. J Comput Chem 2021; 42:365-378. [PMID: 33227163 DOI: 10.1002/jcc.26455] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2020] [Revised: 09/28/2020] [Accepted: 11/05/2020] [Indexed: 11/10/2022]
Abstract
The semiempirical INDO/S Hamiltonian is incorporated into a new release of MOPAC2016. The MOPAC2016 software package has long been at the forefront of semiempirical quantum chemical methods (SEQMs) for small molecules, proteins, and solids and until this release has included only NDDO-type SEQMs. The new code enables the calculation of excited states using the INDO/S Hamiltonian combined with a configuration interaction (CI) approach using single excitations (CIS), single and double excitations (CISD), or multiple reference determinants (MRCI) where reference determinants are generated using a complete active space (CAS) approach. The capacity to perform excited-state calculations beyond the CIS level makes INDO/CI one of the few low-cost computational methods capable of accurately modeling states with substantial double-excitation character. Solvent corrections to the ground-state and excited-state energies can be computed using the COSMO implicit solvent model, incorporating state-specific corrections to the excited states based on the solvent refractive index. This code produces physically reasonable electronic structures, absorption spectra, and solvatochromic shifts at low computational costs for systems up to hundreds of atoms, and for both organic molecules and metal clusters.
Collapse
|
44
|
Jankowska J, Góra RW. Ultrafast nonradiative deactivation of photoexcited 8-oxo-hypoxanthine: a nonadiabatic molecular dynamics study. Phys Chem Chem Phys 2021; 23:1234-1241. [PMID: 33355573 DOI: 10.1039/d0cp05271j] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
In the scientific endeavor to understand the chemical origins of life, the photochemistry of the smallest life building blocks, nucleobases, has been a constant object of focus and intense research. Here, we report the results of the first theoretical study on the photo-properties of an 8-oxo-hypoxanthine molecule, the chromophore of 8-oxo-inosine, which is relevant to the recently proposed, prebiotically plausible synthetic routes to the formation of purine- and pyrimidine-nucleotides. With ab initio and semi-empirical OM2/MRCI quantum-chemistry calculations, we predict a strong photostability of the 8-oxo-hypoxanthine system and see the origin of this effect in ultrafast nonradiative relaxation through puckering of the 6-membered heterocyclic ring.
Collapse
Affiliation(s)
- Joanna Jankowska
- Faculty of Chemistry, University of Warsaw, 02-093 Warsaw, Poland.
| | - Robert W Góra
- Department of Physical and Quantum Chemistry, Wroclaw University of Science and Technology, Wybrzeże Wyspiańskiego 27, 50-370, Wrocław, Poland.
| |
Collapse
|
45
|
Bursch M, Hansen A, Pracht P, Kohn JT, Grimme S. Theoretical study on conformational energies of transition metal complexes. Phys Chem Chem Phys 2021; 23:287-299. [PMID: 33336657 DOI: 10.1039/d0cp04696e] [Citation(s) in RCA: 43] [Impact Index Per Article: 14.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Abstract
Conformational energies are an important chemical property for which a performance assessment of theoretical methods is mandatory. Existing benchmark sets are often limited to biochemical or main group element containing molecules, while organometallic systems are generally less studied. A key problem herein is to routinely generate conformers for these molecules due to their complexity and manifold of possible coordination patterns. In this study we used our recently published CREST protocol [Pracht et al., Phys. Chem. Chem. Phys., 2020, 22, 7169-7192] to generate conformer ensembles for a variety of 40 challenging transition metal containing molecules, which were then used to form a comprehensive conformational energy benchmark set termed TMCONF40. Several low-cost semiempirical, density functional theory (DFT) and force-field methods were compared to high level DLPNO-CCSD(T1) and double-hybrid DFT reference values. Close attention was paid to the energetic ordering of the conformers in the statistical evaluation. With respect to the double-hybrid references, both tested low-cost composite DFT methods produce high Pearson correlation coefficients of rp,mean,B97-3c//B97-3c = 0.922 and rp,mean,PBEh-3c//B97-3c = 0.890, with mean absolute deviations close to or below 1 kcal mol-1. This good performance also holds for a comparison to DLPNO-CCSD(T1) reference energies for a smaller subset termed TMCONF5. Based on DFT geometries, the GFNn-xTB methods yield reasonable Pearson correlation coefficients of rp,mean,GFN1-xTB//B97-3c = 0.617 (MADmean = 2.15 kcal mol-1) and rp,mean,GFN2-xTB//B97-3c = 0.567 (MADmean = 2.68 kcal mol-1), outperforming the widely used PMx methods on the TMCONF40 test set. Employing the low-cost composite DFT method B97-3c on GFN2-xTB geometries yields an slightly improved correlation of rp,mean,B97-3c//GFN2-xTB = 0.632. Furthermore, for 68% of the investigated complexes at least one low-energy conformer was found that is more stable than the respective crystal structure conformation, which signals the importance of conformational studies. General recommendations for the application of the CREST protocol and DFT methods for transition metal conformational energies are given.
Collapse
Affiliation(s)
- Markus Bursch
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstr. 4, 53115 Bonn, Germany.
| | | | | | | | | |
Collapse
|
46
|
Tang D, Shen L, Fang WH. Evaluation of mixed quantum-classical molecular dynamics on cis-azobenzene photoisomerization. Phys Chem Chem Phys 2021; 23:13951-13964. [PMID: 34142685 DOI: 10.1039/d1cp01374b] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The quantitative prediction of nonadiabatic transitions between different electronic states is important to understand ultrafast processes in photochemistry. A variety of mixed quantum-classical molecular dynamics methods such as surface hopping and Ehrenfest mean-field have been developed. However, how to choose an appropriate one from a wide diversity of dynamics algorithms to study a realistic photochemical process is still unclear. In this work, we implemented 30 combinations of different mixed quantum-classical dynamics methods, including 24 surface hopping models with 8 decoherence corrections and 3 momentum rescaling strategies as well as 6 mean-field models. Then we performed numerical investigations by simulating the photoisomerization of cis-azobenzene combined with on-the-fly electronic structure calculations. Predictions of the S1 lifetime and the quantum yield of the photoproduct using different models are distinct. Surface hopping is more robust than mean-field in our test system. Moreover, the choice of momentum rescaling methods in surface hopping brings more significant changes than decoherence corrections, while a large discrepancy between simulation results with different mean-field algorithms has been observed.
Collapse
Affiliation(s)
- Diandong Tang
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China.
| | - Lin Shen
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China.
| | - Wei-Hai Fang
- Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China.
| |
Collapse
|
47
|
Janoš J, Madea D, Mahvidi S, Mujawar T, Švenda J, Suchan J, Slavíček P, Klán P. Conformational Control of the Photodynamics of a Bilirubin Dipyrrinone Subunit: Femtosecond Spectroscopy Combined with Nonadiabatic Simulations. J Phys Chem A 2020; 124:10457-10471. [PMID: 33283519 DOI: 10.1021/acs.jpca.0c08945] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
The photochemistry of bilirubin has been extensively studied due to its importance in the phototherapy of hyperbilirubinemia. In the present work, we investigated the ultrafast photodynamics of a bilirubin dipyrrinone subunit, vinylneoxanthobilirubic acid methyl ester. The photoisomerization and photocyclization reactions of its (E) and (Z) isomers were studied using femtosecond transient absorption spectroscopy and by multireference electronic structure theory, where the nonadiabatic dynamics was modeled with a Landau-Zener surface hopping technique. The following picture has emerged from the combined theoretical and experimental approach. Upon excitation, dipyrrinone undergoes a very fast vibrational relaxation, followed by an internal conversion on a picosecond time scale. The internal conversion leads either to photoisomerization or regeneration of the starting material. Further relaxation dynamics on the order of tens of picoseconds was observed in the ground state. The nonadiabatic simulations revealed a strong conformational control of the photodynamics. The ultrafast formation of a cyclic photochemical product from a less-populated conformer of the studied subunit was predicted by our calculations. We discuss the relevance of the present finding for the photochemistry of native bilirubin. The work has also pointed to the limits of semiclassical nonadiabatic simulations for simulating longer photochemical processes, probably due to the zero-point leakage issue.
Collapse
Affiliation(s)
- Jiří Janoš
- Department of Physical Chemistry, University of Chemistry and Technology, Technická 5, Prague 6, Czech Republic
| | - Dominik Madea
- Department of Chemistry, Faculty of Science, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic.,RECETOX, Faculty of Science, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic
| | - Sadegh Mahvidi
- Department of Chemistry, Faculty of Science, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic.,RECETOX, Faculty of Science, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic
| | - Taufiqueahmed Mujawar
- Department of Chemistry, Faculty of Science, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic
| | - Jakub Švenda
- Department of Chemistry, Faculty of Science, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic
| | - Jiří Suchan
- Department of Physical Chemistry, University of Chemistry and Technology, Technická 5, Prague 6, Czech Republic
| | - Petr Slavíček
- Department of Physical Chemistry, University of Chemistry and Technology, Technická 5, Prague 6, Czech Republic
| | - Petr Klán
- Department of Chemistry, Faculty of Science, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic.,RECETOX, Faculty of Science, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic
| |
Collapse
|
48
|
König G, Riniker S. On the faithfulness of molecular mechanics representations of proteins towards quantum-mechanical energy surfaces. Interface Focus 2020; 10:20190121. [PMID: 33184586 DOI: 10.1098/rsfs.2019.0121] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 08/06/2020] [Indexed: 12/11/2022] Open
Abstract
Force fields based on molecular mechanics (MM) are the main computational tool to study the relationship between protein structure and function at the molecular level. To validate the quality of such force fields, high-level quantum-mechanical (QM) data are employed to test their capability to reproduce the features of all major conformational substates of a series of blocked amino acids. The phase-space overlap between MM and QM is quantified in terms of the average structural reorganization energies over all energy minima. Here, the structural reorganization energy is the MM potential-energy difference between the structure of the respective QM energy minimum and the structure of the closest MM energy minimum. Thus, it serves as a measure for the relative probability of visiting the QM minimum during an MM simulation. We evaluate variants of the AMBER, CHARMM, GROMOS and OPLS biomolecular force fields. In addition, the two blocked amino acids alanine and serine are used to demonstrate the dependence of the measured agreement on the QM method, the phase, and the conformational preferences. Blocked serine serves as an example to discuss possible improvements of the force fields, such as including polarization with Drude particles, or using tailored force fields. The results show that none of the evaluated force fields satisfactorily reproduces all energy minima. By decomposing the average structural reorganization energies in terms of individual energy terms, we can further assess the individual weaknesses of the parametrization strategies of each force field. The dominant problem for most force fields appears to be the van der Waals parameters, followed to a lesser degree by dihedral and bonded terms. Our results show that performing a simple QM energy optimization from an MM-optimized structure can be a first test of the validity of a force field for a particular target molecule.
Collapse
Affiliation(s)
- Gerhard König
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany.,Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Sereina Riniker
- Laboratory of Physical Chemistry, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
49
|
Neugebauer H, Bohle F, Bursch M, Hansen A, Grimme S. Benchmark Study of Electrochemical Redox Potentials Calculated with Semiempirical and DFT Methods. J Phys Chem A 2020; 124:7166-7176. [DOI: 10.1021/acs.jpca.0c05052] [Citation(s) in RCA: 28] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Hagen Neugebauer
- Mulliken Center for Theoretical Chemistry, Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| | - Fabian Bohle
- Mulliken Center for Theoretical Chemistry, Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| | - Markus Bursch
- Mulliken Center for Theoretical Chemistry, Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| |
Collapse
|
50
|
Pecina A, Eyrilmez SM, Köprülüoğlu C, Miriyala VM, Lepšík M, Fanfrlík J, Řezáč J, Hobza P. SQM/COSMO Scoring Function: Reliable Quantum-Mechanical Tool for Sampling and Ranking in Structure-Based Drug Design. Chempluschem 2020; 85:2362-2371. [PMID: 32609421 DOI: 10.1002/cplu.202000120] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2020] [Revised: 05/27/2020] [Indexed: 12/17/2022]
Abstract
Quantum mechanical (QM) methods have been gaining importance in structure-based drug design where a reliable description of protein-ligand interactions is of utmost significance. However, strategies i. e. QM/MM, fragmentation or semiempirical (SQM) methods had to be pursued to overcome the unfavorable scaling of QM methods. Various SQM-based approaches have significantly contributed to the accuracy of docking and improvement of lead compounds. Parametrizations of SQM and implicit solvent methods in our laboratory have been instrumental to obtain a reliable SQM-based scoring function. The experience gained in its application for activity ranking of ligands binding to tens of protein targets resulted in setting up a faster SQM/COSMO scoring approach, which outperforms standard scoring methods in native pose identification for two dozen protein targets with ten thousand poses. Recently, SQM/COSMO was effectively applied in a proof-of-concept study of enrichment in virtual screening. Due to its superior performance, feasibility and chemical generality, we propose the SQM/COSMO approach as an efficient tool in structure-based drug design.
Collapse
Affiliation(s)
- Adam Pecina
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic
| | - Saltuk M Eyrilmez
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Palacky University, 771 46, Olomouc, Czech Republic
| | - Cemal Köprülüoğlu
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Palacky University, 771 46, Olomouc, Czech Republic
| | - Vijay Madhav Miriyala
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic
| | - Martin Lepšík
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic
| | - Jindřich Fanfrlík
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic
| | - Jan Řezáč
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic
| | - Pavel Hobza
- Institute of Organic Chemistry, and Biochemistry of Czech Academy of Sciences, Flemingovo namesti 2, 166 10, Prague, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Palacky University, 771 46, Olomouc, Czech Republic
| |
Collapse
|