1
|
Chen WK, Fang WH, Cui G. Extending multi-layer energy-based fragment method for excited-state calculations of large covalently bonded fragment systems. J Chem Phys 2023; 158:044110. [PMID: 36725521 DOI: 10.1063/5.0129458] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023] Open
Abstract
Recently, we developed a low-scaling Multi-Layer Energy-Based Fragment (MLEBF) method for accurate excited-state calculations and nonadiabatic dynamics simulations of nonbonded fragment systems. In this work, we extend the MLEBF method to treat covalently bonded fragment ones. The main idea is cutting a target system into many fragments according to chemical properties. Fragments with dangling bonds are first saturated by chemical groups; then, saturated fragments, together with the original fragments without dangling bonds, are grouped into different layers. The accurate total energy expression is formulated with the many-body energy expansion theory, in combination with the inclusion-exclusion principle that is used to delete the contribution of chemical groups introduced to saturate dangling bonds. Specifically, in a two-layer MLEBF model, the photochemically active and inert layers are calculated with high-level and efficient electronic structure methods, respectively. Intralayer and interlayer energies can be truncated at the two- or three-body interaction level. Subsequently, through several systems, including neutral and charged covalently bonded fragment systems, we demonstrate that MLEBF can provide accurate ground- and excited-state energies and gradients. Finally, we realize the structure, conical intersection, and path optimizations by combining our MLEBF program with commercial and free packages, e.g., ASE and SciPy. These developments make MLEBF a practical and reliable tool for studying complex photochemical and photophysical processes of large nonbonded and bonded fragment systems.
Collapse
Affiliation(s)
- Wen-Kai Chen
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Wei-Hai Fang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Ganglong Cui
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| |
Collapse
|
2
|
Hellmers J, Hedegård ED, König C. Fragmentation-Based Decomposition of a Metalloenzyme-Substrate Interaction: A Case Study for a Lytic Polysaccharide Monooxygenase. J Phys Chem B 2022; 126:5400-5412. [PMID: 35833656 DOI: 10.1021/acs.jpcb.2c02883] [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
We present a novel decomposition scheme for electronic interaction energies based on the flexible formulation of fragmentation schemes through fragment combination ranges (FCRs; J. Chem. Phys., 2021, 155, 164105). We devise a clear additive decomposition with contribution of nondisjoint fragments and correction terms for overlapping fragments and apply this scheme to the metalloenzyme-substrate complex of a lytic polysaccharide monooxygenase (LPMO) with an oligosaccharide. By this, we further illustrate the straightforward adaptability of the FCR-based schemes to novel systems. Our calculations suggest that the description of the electronic structure is a larger error source than the fragmentation scheme. In particular, we find a large impact of the basis set size on the interaction energies. Still, the introduction of three-body interaction terms in the fragmentation setup improves the agreement to the supermolecular reference. Yet, the qualitative results for the decomposition scheme with two-body terms only largely agree within the investigated electronic-structure approaches and basis sets, which are B97-3c, DFT (TPSS and B3LYP), and MP2 methods. The overlap contributions are found to be small, allowing analysis of the interaction energy into individual amino acid residues: We find a particularly strong interaction between the substrate and the LPMO copper active site.
Collapse
Affiliation(s)
- Janine Hellmers
- Institute of Physical Chemistry and Electrochemistry, Leibniz University Hannover, 30167 Hannover, Germany
| | - Erik Donovan Hedegård
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, 5230 Odense, Denmark
| | - Carolin König
- Institute of Physical Chemistry and Electrochemistry, Leibniz University Hannover, 30167 Hannover, Germany
| |
Collapse
|
3
|
Jindal S, Hsu PJ, Phan HT, Tsou PK, Kuo JL. Capturing the potential energy landscape of large size molecular clusters from atomic interactions up to a 4-body system using deep learning. Phys Chem Chem Phys 2022; 24:27263-27276. [DOI: 10.1039/d2cp04441b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
We propose a new method that utilizes the database of stable conformers and borrow the fragmentation concept of many-body-expansion (MBE) methods in ab initio methods to train a deep-learning machine learning (ML) model using SchNet.
Collapse
Affiliation(s)
- Shweta Jindal
- Institute of Atomic and Molecular Sciences, Academia Sinica, No. 1 Roosevelt Road, Section 4, Daan District, Taipei City 10617, Taiwan
| | - Po-Jen Hsu
- Institute of Atomic and Molecular Sciences, Academia Sinica, No. 1 Roosevelt Road, Section 4, Daan District, Taipei City 10617, Taiwan
| | - Huu Trong Phan
- Institute of Atomic and Molecular Sciences, Academia Sinica, No. 1 Roosevelt Road, Section 4, Daan District, Taipei City 10617, Taiwan
- Molecular Science and Technology Program, Taiwan International Graduate Program, Academia Sinica, Taipei, 11529, Taiwan
- Department of Chemistry, National Tsing Hua University, Hsinchu 30013, Taiwan
| | - Pei-Kang Tsou
- Institute of Atomic and Molecular Sciences, Academia Sinica, No. 1 Roosevelt Road, Section 4, Daan District, Taipei City 10617, Taiwan
| | - Jer-Lai Kuo
- Institute of Atomic and Molecular Sciences, Academia Sinica, No. 1 Roosevelt Road, Section 4, Daan District, Taipei City 10617, Taiwan
- Department of Chemistry, National Tsing Hua University, Hsinchu 30013, Taiwan
- Molecular Science and Technology, National Taiwan University, Section 4, Daan District, Taipei City 10617, Taiwan
| |
Collapse
|
4
|
Liu KY, Carter-Fenk K, Herbert JM. Self-consistent charge embedding at very low cost, with application to symmetry-adapted perturbation theory. J Chem Phys 2019; 151:031102. [DOI: 10.1063/1.5111869] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Affiliation(s)
- Kuan-Yu Liu
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Kevin Carter-Fenk
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - John M. Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| |
Collapse
|
5
|
Xu M, He X, Zhu T, Zhang JZH. A Fragment Quantum Mechanical Method for Metalloproteins. J Chem Theory Comput 2019; 15:1430-1439. [PMID: 30620584 DOI: 10.1021/acs.jctc.8b00966] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
An accurate energy calculation of metalloprotein is of crucial importance and also a theoretical challenge. In this work, a metal molecular fractionation with conjugate caps (metal-MFCC) approach is developed for efficient linear-scaling quantum calculation of potential energy and atomic forces of metalloprotein. In this approach, the potential energy of a given protein is calculated by a linear combination of potential energies of the neighboring residues, two-body interaction energy between non-neighboring residues that are spatially in close contact and the potential energy of the metal binding group. The calculation of each fragment is embedded in a field of point charges representing the remaining protein environment. Numerical studies were carried out to check the performance of this method, and the calculated potential energies and atomic forces all show excellent agreement with the full system calculations at the M06-2X/6-31G(d) level. By combining the energy calculation with molecular dynamic simulation, we performed an ab initio structural optimization for a zinc finger protein with high efficiency. The present metal-MFCC approach is linear-scaling with a low prefactor and trivially parallelizable. The individual fragment typically contains about 50 atoms, and it is thus possible to be calculated at higher levels of the quantum chemistry method. This fragment method can be routinely applied to perform structural optimization and ab initio molecular dynamic simulation for metalloproteins of any size.
Collapse
Affiliation(s)
- Mingyuan Xu
- Shanghai Engineering Research Center of Molecular Therapeutics & New Drug Development, School of Chemistry and Molecular Engineering , East China Normal University , Shanghai , 200062 , China
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics & New Drug Development, School of Chemistry and Molecular Engineering , East China Normal University , Shanghai , 200062 , China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai , Shanghai , 200062 , China
| | - Tong Zhu
- Shanghai Engineering Research Center of Molecular Therapeutics & New Drug Development, School of Chemistry and Molecular Engineering , East China Normal University , Shanghai , 200062 , China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai , Shanghai , 200062 , China
| | - John Z H Zhang
- Shanghai Engineering Research Center of Molecular Therapeutics & New Drug Development, School of Chemistry and Molecular Engineering , East China Normal University , Shanghai , 200062 , China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai , Shanghai , 200062 , China.,Department of Chemistry , New York University , New York 10003 , United States
| |
Collapse
|
6
|
Yuan D, Li Y, Li W, Li S. Structures and properties of large supramolecular coordination complexes predicted with the generalized energy-based fragmentation method. Phys Chem Chem Phys 2018; 20:28894-28902. [DOI: 10.1039/c8cp05548c] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The generalized energy-based fragmentation (GEBF) method has been extended to facilitate ab initio calculations of large supramolecular coordination complexes.
Collapse
Affiliation(s)
- Dandan Yuan
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing 210023
| | - Yunzhi Li
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing 210023
| | - Wei Li
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing 210023
| | - Shuhua Li
- School of Chemistry and Chemical Engineering
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education
- Institute of Theoretical and Computational Chemistry
- Nanjing University
- Nanjing 210023
| |
Collapse
|
7
|
Liu J, Qi LW, Zhang JZH, He X. Fragment Quantum Mechanical Method for Large-Sized Ion–Water Clusters. J Chem Theory Comput 2017; 13:2021-2034. [DOI: 10.1021/acs.jctc.7b00149] [Citation(s) in RCA: 43] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Affiliation(s)
- Jinfeng Liu
- State
Key Laboratory of Natural Medicines, Department of Basic Medicine
and Clinical Pharmacy, China Pharmaceutical University, Nanjing, 210009, China
| | - Lian-Wen Qi
- State
Key Laboratory of Natural Medicines, Department of Basic Medicine
and Clinical Pharmacy, China Pharmaceutical University, Nanjing, 210009, China
| | - John Z. H. Zhang
- School
of Chemistry and Molecular Engineering, East China Normal University, Shanghai, 200062, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, 200062, China
- Department
of Chemistry, New York University, New York, New York 10003, United States
| | - Xiao He
- School
of Chemistry and Molecular Engineering, East China Normal University, Shanghai, 200062, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, 200062, China
| |
Collapse
|
8
|
Raghavachari K, Saha A. Accurate Composite and Fragment-Based Quantum Chemical Models for Large Molecules. Chem Rev 2015; 115:5643-77. [PMID: 25849163 DOI: 10.1021/cr500606e] [Citation(s) in RCA: 183] [Impact Index Per Article: 20.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Krishnan Raghavachari
- Department of Chemistry, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Arjun Saha
- Department of Chemistry, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| |
Collapse
|
9
|
Affiliation(s)
- Michael A Collins
- †Research School of Chemistry, Australian National University, Canberra, ACT 0200, Australia
| | - Ryan P A Bettens
- ‡Department of Chemistry, National University of Singapore, 3 Science Drive 3, Singapore 117543, Singapore
| |
Collapse
|
10
|
Wang B, Yang KR, Xu X, Isegawa M, Leverentz HR, Truhlar DG. Quantum mechanical fragment methods based on partitioning atoms or partitioning coordinates. Acc Chem Res 2014; 47:2731-8. [PMID: 24841937 DOI: 10.1021/ar500068a] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Conspectus The development of more efficient and more accurate ways to represent reactive potential energy surfaces is a requirement for extending the simulation of large systems to more complex systems, longer-time dynamical processes, and more complete statistical mechanical sampling. One way to treat large systems is by direct dynamics fragment methods. Another way is by fitting system-specific analytic potential energy functions with methods adapted to large systems. Here we consider both approaches. First we consider three fragment methods that allow a given monomer to appear in more than one fragment. The first two approaches are the electrostatically embedded many-body (EE-MB) expansion and the electrostatically embedded many-body expansion of the correlation energy (EE-MB-CE), which we have shown to yield quite accurate results even when one restricts the calculations to include only electrostatically embedded dimers. The third fragment method is the electrostatically embedded molecular tailoring approach (EE-MTA), which is more flexible than EE-MB and EE-MB-CE. We show that electrostatic embedding greatly improves the accuracy of these approaches compared with the original unembedded approaches. Quantum mechanical fragment methods share with combined quantum mechanical/molecular mechanical (QM/MM) methods the need to treat a quantum mechanical fragment in the presence of the rest of the system, which is especially challenging for those parts of the rest of the system that are close to the boundary of the quantum mechanical fragment. This is a delicate matter even for fragments that are not covalently bonded to the rest of the system, but it becomes even more difficult when the boundary of the quantum mechanical fragment cuts a bond. We have developed a suite of methods for more realistically treating interactions across such boundaries. These methods include redistributing and balancing the external partial atomic charges and the use of tuned fluorine atoms for capping dangling bonds, and we have shown that they can greatly improve the accuracy. Finally we present a new approach that goes beyond QM/MM by combining the convenience of molecular mechanics with the accuracy of fitting a potential function to electronic structure calculations on a specific system. To make the latter practical for systems with a large number of degrees of freedom, we developed a method to interpolate between local internal-coordinate fits to the potential energy. A key issue for the application to large systems is that rather than assigning the atoms or monomers to fragments, we assign the internal coordinates to reaction, secondary, and tertiary sets. Thus, we make a partition in coordinate space rather than atom space. Fits to the local dependence of the potential energy on tertiary coordinates are arrayed along a preselected reaction coordinate at a sequence of geometries called anchor points; the potential energy function is called an anchor points reactive potential. Electrostatically embedded fragment methods and the anchor points reactive potential, because they are based on treating an entire system by quantum mechanical electronic structure methods but are affordable for large and complex systems, have the potential to open new areas for accurate simulations where combined QM/MM methods are inadequate.
Collapse
Affiliation(s)
- Bo Wang
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, 207 Pleasant St. SE, Minneapolis, Minnesota 55455-0431, United States
| | - Ke R. Yang
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, 207 Pleasant St. SE, Minneapolis, Minnesota 55455-0431, United States
| | - Xuefei Xu
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, 207 Pleasant St. SE, Minneapolis, Minnesota 55455-0431, United States
| | - Miho Isegawa
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, 207 Pleasant St. SE, Minneapolis, Minnesota 55455-0431, United States
| | - Hannah R. Leverentz
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, 207 Pleasant St. SE, Minneapolis, Minnesota 55455-0431, United States
| | - Donald G. Truhlar
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, 207 Pleasant St. SE, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
11
|
Ouyang JF, Cvitkovic MW, Bettens RPA. Trouble with the Many-Body Expansion. J Chem Theory Comput 2014; 10:3699-707. [DOI: 10.1021/ct500396b] [Citation(s) in RCA: 41] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- John F. Ouyang
- Department of Chemistry, National University of Singapore, 3 Science
Drive 3, Singapore 117543
| | - Milan W. Cvitkovic
- Department of Chemistry, National University of Singapore, 3 Science
Drive 3, Singapore 117543
| | - Ryan P. A. Bettens
- Department of Chemistry, National University of Singapore, 3 Science
Drive 3, Singapore 117543
| |
Collapse
|
12
|
Saha A, Raghavachari K. Dimers of Dimers (DOD): A New Fragment-Based Method Applied to Large Water Clusters. J Chem Theory Comput 2013; 10:58-67. [DOI: 10.1021/ct400472v] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Arjun Saha
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Krishnan Raghavachari
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| |
Collapse
|
13
|
Ramabhadran RO, Raghavachari K. Extrapolation to the Gold-Standard in Quantum Chemistry: Computationally Efficient and Accurate CCSD(T) Energies for Large Molecules Using an Automated Thermochemical Hierarchy. J Chem Theory Comput 2013; 9:3986-94. [PMID: 26592394 DOI: 10.1021/ct400465q] [Citation(s) in RCA: 44] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Abstract
The CCSD(T) method is known as the gold-standard in quantum chemistry and has been the method of choice in quantum chemistry for over 20 years to obtain accurate bond energies and molecular properties. Its computational cost formally scales as the seventh power of the size of the system and can be prohibitive for large molecules. As part of our efforts to reduce the computational cost of the CCSD(T) method yet retain its accuracy, we present a simple, efficient, and user-friendly protocol to extrapolate to CCSD(T) energies in conjunction with MP2 energies. The method is based on the automated error-canceling thermochemical hierarchy previously developed by us called the Connectivity-Based Hierarchy (CBH). For a test set containing 30 diverse nonaromatic organic molecules and biomonomers, we obtain highly accurate extrapolated CCSD(T) energies (with a mean absolute error of only 0.2-0.3 kcal/mol with different basis-set). Additionally, the work also features the successful extrapolation to CCSD energies using a similar protocol.
Collapse
Affiliation(s)
| | - Krishnan Raghavachari
- Department of Chemistry, Indiana University , Bloomington, Indiana 47405, United States
| |
Collapse
|