1
|
Ballesteros F, Lao KU. Analysis of two overlapping fragmentation approaches in density matrix construction: GMBE-DM vs. ADMA. Phys Chem Chem Phys 2024; 26:4386-4394. [PMID: 38236152 DOI: 10.1039/d3cp05759c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2024]
Abstract
In this study, we conduct a comparative analysis of two density matrix construction methods: the generalized many-body expansion for building density matrices (GMBE-DM) based on the set-theoretical principle of inclusion/exclusion and the adjustable density matrix assembler (ADMA) based on the Mulliken-Mezey ansatz. We apply these methods to various noncovalent clusters, including water clusters, ion-water clusters, and ion-pair clusters, using both small 6-31G(d) and large def2-TZVPPD basis sets. Our findings reveal that the GMBE-DM method, particularly when combined with the purification scheme and truncation at the one-body level [GMBE(1)-DM-P], exhibits superior performance across all test systems and basis sets. In contrast, all ADMA set of methods show reasonable results only with small and compact basis sets. For example, GMBE(1)-DM-P outperforms the best ADMA method by at least 4 and 16 times with small and large basis sets, respectively, in the case of (H2O)N=6-55. This highlights the significance of the basis set choice for ADMA, which is even more critical than the fragmentation scheme, such as the size of subsystems, while GMBE-DM consistently produces accurate results irrespective of the chosen basis set. Consequently, the efficient and robust GMBE(1)-DM-P approach is recommended as a fragmentation method for generating accurate absolute and relative energies across different binding patterns and basis sets for noncovalent clusters.
Collapse
Affiliation(s)
| | - Ka Un Lao
- Department of Chemistry, Virginia Commonwealth University, Richmond, VA, USA.
| |
Collapse
|
2
|
Vornweg JR, Wolter M, Jacob CR. A simple and consistent quantum-chemical fragmentation scheme for proteins that includes two-body contributions. J Comput Chem 2023; 44:1634-1644. [PMID: 37171574 DOI: 10.1002/jcc.27114] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2023] [Revised: 03/28/2023] [Accepted: 03/30/2023] [Indexed: 05/13/2023]
Abstract
The Molecular Fractionation with Conjugate Caps (MFCC) method is a popular fragmentation method for the quantum-chemical treatment of proteins. However, it does not account for interactions between the amino acid fragments, such as intramolecular hydrogen bonding. Here, we present a combination of the MFCC fragmentation scheme with a second-order many-body expansion (MBE) that consistently accounts for all fragment-fragment, fragment-cap, and cap-cap interactions, while retaining the overall simplicity of the MFCC scheme with its chemically meaningful fragments. We show that with the resulting MFCC-MBE(2) scheme, the errors in the total energies of selected polypeptides and proteins can be reduced by up to one order of magnitude and relative energies of different protein conformers can be predicted accurately.
Collapse
Affiliation(s)
- Johannes R Vornweg
- Institute of Physical and Theoretical Chemistry, Technische Universität Braunschweig, Braunschweig, Germany
| | - Mario Wolter
- Institute of Physical and Theoretical Chemistry, Technische Universität Braunschweig, Braunschweig, Germany
| | - Christoph R Jacob
- Institute of Physical and Theoretical Chemistry, Technische Universität Braunschweig, Braunschweig, Germany
| |
Collapse
|
3
|
Nakai H, Kobayashi M, Yoshikawa T, Seino J, Ikabata Y, Nishimura Y. Divide-and-Conquer Linear-Scaling Quantum Chemical Computations. J Phys Chem A 2023; 127:589-618. [PMID: 36630608 DOI: 10.1021/acs.jpca.2c06965] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
Fragmentation and embedding schemes are of great importance when applying quantum-chemical calculations to more complex and attractive targets. The divide-and-conquer (DC)-based quantum-chemical model is a fragmentation scheme that can be connected to embedding schemes. This feature article explains several DC-based schemes developed by the authors over the last two decades, which was inspired by the pioneering study of DC self-consistent field (SCF) method by Yang and Lee (J. Chem. Phys. 1995, 103, 5674-5678). First, the theoretical aspects of the DC-based SCF, electron correlation, excited-state, and nuclear orbital methods are described, followed by the two-component relativistic theory, quantum-mechanical molecular dynamics simulation, and the introduction of three programs, including DC-based schemes. Illustrative applications confirmed the accuracy and feasibility of the DC-based schemes.
Collapse
Affiliation(s)
- Hiromi Nakai
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo169-8555, Japan.,Waseda Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo169-8555, Japan
| | - Masato Kobayashi
- Department of Chemistry, Faculty of Science, Hokkaido University, Kita 10 Nishi 8, Kita-ku, Sapporo, Hokkaido060-0810, Japan.,Institute for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Kita 21 Nishi 10, Kita-ku, Sapporo, Hokkaido001-0021, Japan
| | - Takeshi Yoshikawa
- Faculty of Pharmaceutical Sciences, Toho University, 2-2-1 Miyama, Funabashi, Chiba274-8510, Japan
| | - Junji Seino
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo169-8555, Japan.,Waseda Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo169-8555, Japan
| | - Yasuhiro Ikabata
- Information and Media Center, Toyohashi University of Technology, 1-1 Hibarigaoka, Tempaku-cho, Toyohashi, Aichi441-8580, Japan.,Department of Computer Science and Engineering, Toyohashi University of Technology, 1-1 Hibarigaoka, Tempaku-cho, Toyohashi, Aichi441-8580, Japan
| | - Yoshifumi Nishimura
- Waseda Research Institute for Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo169-8555, Japan
| |
Collapse
|
4
|
Ballesteros F, Lao KU. Accelerating the Convergence of Self-Consistent Field Calculations Using the Many-Body Expansion. J Chem Theory Comput 2021; 18:179-191. [PMID: 34881906 DOI: 10.1021/acs.jctc.1c00765] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
The balance between cost-effective and sufficiently accurate methods represents the proverbial "promised land" for quantum chemistry calculations. The burden thus falls upon theoretical and computational chemists to provide such alternatives to mitigate the issues that arise from the employ of finite computing resources. In this paper, we attempt to demonstrate the importance of the quality of the initial guess for the self-consistent field (SCF) calculation when considering cost reduction techniques. We broach this challenge by using the many body expansion (MBE) to yield high quality density matrices (DMs) which, in turn, are applied as an SCF initial guess. The MBE-DM approaches combined with purification schemes and distance-based cutoff schemes can serve as initial guesses to reduce the SCF cycles necessary for convergence or derive energy directly through one Fock build. To this end, four unique types of clusters including water clusters, fluoride anion water clusters, sodium cation water clusters, and ammonium-bisulfate salt clusters have been used to test the performance of MBE-DM where its truncation at three-body expansion, MBE(3)-DM, shows vast improvement for those four clusters with reductions in the number of SCF cycles up to 40% as compared with the traditional superposition of atomic densities (SAD) guess. Other types of typical initial guesses, superposition of atomic potentials (SAP) and basis set projection (BSP), perform much worse than MBE-DM and SAD. In addition, the MBE-DM shows consistency across an array of fragment types irrespective of charges, size, level of theory, and basis set selection. Through MBE(3)-DM with the distance cutoff and the average purification scheme, the energy can be obtained directly with a mere 3.2 mH of the mean absolute deviation (MAD) for (H2O)N=6-55 which is at least 73 times better than the energy prediction using the typical initial guesses (SAD, SAP, and BSP). The corresponding MAD per monomer is only 0.14 mH which reaches the threshold of the "dynamical accuracy". The promising results of the methods outlined in this paper not only indicate two direct routes for computational cost reduction but also lay the possible foundation for composite techniques (i.e., ab initio sampling) that make best use of near converged values as their starting point.
Collapse
Affiliation(s)
- Francisco Ballesteros
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284, United States
| | - Ka Un Lao
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284, United States
| |
Collapse
|
5
|
Reinholdt P, Jørgensen FK, Kongsted J, Olsen JMH. Polarizable Density Embedding for Large Biomolecular Systems. J Chem Theory Comput 2020; 16:5999-6006. [PMID: 32991163 DOI: 10.1021/acs.jctc.0c00763] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
We present an efficient and robust fragment-based quantum-classical embedding model capable of accurately capturing effects from complex environments such as proteins and nucleic acids. This is realized by combining the molecular fractionation with conjugate caps (MFCC) procedure with the polarizable density embedding (PDE) model at the level of Fock matrix construction. The PDE contributions to the Fock matrix of the core region are constructed using the local molecular basis of the individual fragments rather than the supermolecular basis of the entire system. Thereby, we avoid complications associated with the application of the MFCC procedure on environment quantities such as electronic densities and molecular-orbital energies. Moreover, the computational cost associated with solving self-consistent field (SCF) equations of the core region remains unchanged from that of purely classical polarized embedding models. We analyze the performance of the resulting model in terms of the reproduction of the electrostatic potential of an insulin monomer protein and further in the context of solving problems related to electron spill-out. Finally, we showcase the model for the calculation of one- and two-photon properties of the Nile red molecule in a protein environment. Based on our analyses, we find that the combination of the MFCC approach with the PDE model is an efficient, yet accurate approach for calculating molecular properties of molecules embedded in structured biomolecular environments.
Collapse
Affiliation(s)
- Peter Reinholdt
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| | - Frederik Kamper Jørgensen
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| | - Jacob Kongsted
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| | - Jógvan Magnus Haugaard Olsen
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, UiT The Arctic University of Norway, N-9037 Tromsø, Norway.,Department of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark
| |
Collapse
|
6
|
Abstract
Since the introduction of the fragment molecular orbital method 20 years ago, fragment-based approaches have occupied a small but growing niche in quantum chemistry. These methods decompose a large molecular system into subsystems small enough to be amenable to electronic structure calculations, following which the subsystem information is reassembled in order to approximate an otherwise intractable supersystem calculation. Fragmentation sidesteps the steep rise (with respect to system size) in the cost of ab initio calculations, replacing it with a distributed cost across numerous computer processors. Such methods are attractive, in part, because they are easily parallelizable and therefore readily amenable to exascale computing. As such, there has been hope that distributed computing might offer the proverbial "free lunch" in quantum chemistry, with the entrée being high-level calculations on very large systems. While fragment-based quantum chemistry can count many success stories, there also exists a seedy underbelly of rarely acknowledged problems. As these methods begin to mature, it is time to have a serious conversation about what they can and cannot be expected to accomplish in the near future. Both successes and challenges are highlighted in this Perspective.
Collapse
Affiliation(s)
- John M Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| |
Collapse
|
7
|
Sun CL, Liu LP, Tian F, Ding F, Wang LW. Charge-patching method for the calculation of electronic structure of polypeptides. Phys Chem Chem Phys 2018; 20:23301-23310. [DOI: 10.1039/c8cp01803k] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Based on the CPM method, the charge densities of polypeptides can be generated and their electronic structure can be further calculated.
Collapse
Affiliation(s)
- Chang-Liang Sun
- Center of Physical Chemistry Test
- Shenyang University of Chemical Technology
- Shenyang 110142
- People's Republic of China
- Materials Science Division
| | - Li-Ping Liu
- Materials Science Division
- Lawrence Berkeley National Laboratory
- Berkeley
- USA
- School of Physics
| | - Fubo Tian
- Materials Science Division
- Lawrence Berkeley National Laboratory
- Berkeley
- USA
- College of Physics
| | - Fu Ding
- Center of Physical Chemistry Test
- Shenyang University of Chemical Technology
- Shenyang 110142
- People's Republic of China
| | - Lin-Wang Wang
- Materials Science Division
- Lawrence Berkeley National Laboratory
- Berkeley
- USA
| |
Collapse
|
8
|
|
9
|
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
|
10
|
Jin X, Zhang JZH, He X. Full QM Calculation of RNA Energy Using Electrostatically Embedded Generalized Molecular Fractionation with Conjugate Caps Method. J Phys Chem A 2017; 121:2503-2514. [PMID: 28264557 DOI: 10.1021/acs.jpca.7b00859] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
In this study, the electrostatically embedded generalized molecular fractionation with conjugate caps (concaps) method (EE-GMFCC) was employed for efficient linear-scaling quantum mechanical (QM) calculation of total energies of RNAs. In the EE-GMFCC approach, the total energy of RNA is calculated by taking a proper combination of the QM energy of each nucleotide-centric fragment with large caps or small caps (termed EE-GMFCC-LC and EE-GMFCC-SC, respectively) deducted by the energies of concaps. The two-body QM interaction energy between non-neighboring ribonucleotides which are spatially in close contact are also taken into account for the energy calculation. Numerical studies were carried out to calculate the total energies of a number of RNAs using the EE-GMFCC-LC and EE-GMFCC-SC methods at levels of the Hartree-Fock (HF) method, density functional theory (DFT), and second-order many-body perturbation theory (MP2), respectively. The results show that the efficiency of the EE-GMFCC-SC method is about 3 times faster than the EE-GMFCC-LC method with minimal accuracy sacrifice. The EE-GMFCC-SC method is also applied for relative energy calculations of 20 different conformers of two RNA systems using HF and DFT, respectively. Both single-point and relative energy calculations demonstrate that the EE-GMFCC method has deviations from the full system results of only a few kcal/mol.
Collapse
Affiliation(s)
- Xinsheng Jin
- School of Chemistry and Molecular Engineering, East China Normal University , Shanghai 200062, 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
|
11
|
Sahu N, Gadre SR. Accurate vibrational spectra via molecular tailoring approach: a case study of water clusters at MP2 level. J Chem Phys 2016; 142:014107. [PMID: 25573553 DOI: 10.1063/1.4905004] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In spite of the recent advents in parallel algorithms and computer hardware, high-level calculation of vibrational spectra of large molecules is still an uphill task. To overcome this, significant effort has been devoted to the development of new algorithms based on fragmentation methods. The present work provides the details of an efficient and accurate procedure for computing the vibrational spectra of large clusters employing molecular tailoring approach (MTA). The errors in the Hessian matrix elements and dipole derivatives arising due to the approximation nature of MTA are reduced by grafting the corrections from a smaller basis set. The algorithm has been tested out for obtaining vibrational spectra of neutral and charged water clusters at Møller-Plesset second order level of theory, and benchmarking them against the respective full calculation (FC) and/or experimental results. For (H2O)16 clusters, the estimated vibrational frequencies are found to differ by a maximum of 2 cm(-1) with reference to the corresponding FC values. Unlike the FC, the MTA-based calculations including grafting procedure can be performed on a limited hardware, yet take a fraction of the FC time. The present methodology, thus, opens a possibility of the accurate estimation of the vibrational spectra of large molecular systems, which is otherwise impossible or formidable.
Collapse
Affiliation(s)
- Nityananda Sahu
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208 016, India
| | - Shridhar R Gadre
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208 016, India
| |
Collapse
|
12
|
Liu J, Herbert JM. Pair-Pair Approximation to the Generalized Many-Body Expansion: An Alternative to the Four-Body Expansion for ab Initio Prediction of Protein Energetics via Molecular Fragmentation. J Chem Theory Comput 2016; 12:572-84. [PMID: 26730608 DOI: 10.1021/acs.jctc.5b00955] [Citation(s) in RCA: 43] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We introduce a "pair-pair" approximation to the generalized many-body expansion (pp-GMBE) as an approximation to a traditional four-body expansion, the latter of which is accurate but quickly becomes numerically unstable and ultimately intractable as the number of "bodies" (fragments) increases. The pp-GMBE method achieves a good balance between accuracy and efficiency by defining significant fragment pairs and then fragment quartets. An efficient fragmentation scheme is introduced for proteins such that the largest subsystems contain about 60 atoms. Application of the pp-GMBE method to proteins with as many as 70 amino acids (1142 atoms) reveals that pp-GMBE energies are quite faithful to those obtained when the same level of density functional theory is applied to the entire macromolecule. When combined with embedding charges obtained from natural population analysis, the pp-GMBE approach affords absolute energies that differ by 1-3 kcal/mol from full supersystem results, but it yields conformational energy profiles that are practically indistinguishable from the supersystem calculation at the same level of theory.
Collapse
Affiliation(s)
- Jie Liu
- Department of Chemistry and Biochemistry, The Ohio State University , Columbus, Ohio 43210, United States
| | - John M Herbert
- Department of Chemistry and Biochemistry, The Ohio State University , Columbus, Ohio 43210, United States
| |
Collapse
|
13
|
Liu J, Zhang JZH, He X. Fragment quantum chemical approach to geometry optimization and vibrational spectrum calculation of proteins. Phys Chem Chem Phys 2016; 18:1864-75. [DOI: 10.1039/c5cp05693d] [Citation(s) in RCA: 52] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Geometry optimization and vibrational spectra (infrared and Raman spectra) calculations of proteins are carried out by a quantum chemical approach using the EE-GMFCC (electrostatically embedded generalized molecular fractionation with conjugate caps) method (J. Phys. Chem. A, 2013, 117, 7149).
Collapse
Affiliation(s)
- Jinfeng Liu
- State Key Laboratory of Precision Spectroscopy
- Institute of Theoretical and Computational Science
- College of Chemistry and Molecular Engineering
- East China Normal University
- Shanghai
| | - John Z. H. Zhang
- State Key Laboratory of Precision Spectroscopy
- Institute of Theoretical and Computational Science
- College of Chemistry and Molecular Engineering
- East China Normal University
- Shanghai
| | - Xiao He
- State Key Laboratory of Precision Spectroscopy
- Institute of Theoretical and Computational Science
- College of Chemistry and Molecular Engineering
- East China Normal University
- Shanghai
| |
Collapse
|
14
|
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
|
15
|
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
|
16
|
Liu J, Wang X, Zhang JZH, He X. Calculation of protein–ligand binding affinities based on a fragment quantum mechanical method. RSC Adv 2015. [DOI: 10.1039/c5ra20185c] [Citation(s) in RCA: 38] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022] Open
Abstract
An efficient fragment-based quantum mechanical method has been successfully applied for reliable prediction of protein–ligand binding affinities.
Collapse
Affiliation(s)
- Jinfeng Liu
- State Key Laboratory of Precision Spectroscopy
- Institute of Theoretical and Computational Science
- College of Chemistry and Molecular Engineering
- East China Normal University
- Shanghai
| | - Xianwei Wang
- Center for Optics & Optoelectronics Research
- College of Science
- Zhejiang University of Technology
- Hangzhou
- China
| | - John Z. H. Zhang
- State Key Laboratory of Precision Spectroscopy
- Institute of Theoretical and Computational Science
- College of Chemistry and Molecular Engineering
- East China Normal University
- Shanghai
| | - Xiao He
- State Key Laboratory of Precision Spectroscopy
- Institute of Theoretical and Computational Science
- College of Chemistry and Molecular Engineering
- East China Normal University
- Shanghai
| |
Collapse
|
17
|
He X, Zhu T, Wang X, Liu J, Zhang JZH. Fragment quantum mechanical calculation of proteins and its applications. Acc Chem Res 2014; 47:2748-57. [PMID: 24851673 DOI: 10.1021/ar500077t] [Citation(s) in RCA: 145] [Impact Index Per Article: 14.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Conspectus The desire to study molecular systems that are much larger than what the current state-of-the-art ab initio or density functional theory methods could handle has naturally led to the development of novel approximate methods, including semiempirical approaches, reduced-scaling methods, and fragmentation methods. The major computational limitation of ab initio methods is the scaling problem, because the cost of ab initio calculation scales nth power or worse with system size. In the past decade, the fragmentation approach based on chemical locality has opened a new door for developing linear-scaling quantum mechanical (QM) methods for large systems and for applications to large molecular systems such as biomolecules. The fragmentation approach is highly attractive from a computational standpoint. First, the ab initio calculation of individual fragments can be conducted almost independently, which makes it suitable for massively parallel computations. Second, the electron properties, such as density and energy, are typically combined in a linear fashion to reproduce those for the entire molecular system, which makes the overall computation scale linearly with the size of the system. In this Account, two fragmentation methods and their applications to macromolecules are described. They are the electrostatically embedded generalized molecular fractionation with conjugate caps (EE-GMFCC) method and the automated fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) approach. The EE-GMFCC method is developed from the MFCC approach, which was initially used to obtain accurate protein-ligand QM interaction energies. The main idea of the MFCC approach is that a pair of conjugate caps (concaps) is inserted at the location where the subsystem is divided by cutting the chemical bond. In addition, the pair of concaps is fused to form molecular species such that the overcounted effect from added concaps can be properly removed. By introducing the electrostatic embedding field in each fragment calculation and two-body interaction energy correction on top of the MFCC approach, the EE-GMFCC method is capable of accurately reproducing the QM molecular properties (such as the dipole moment, electron density, and electrostatic potential), the total energy, and the electrostatic solvation energy from full system calculations for proteins. On the other hand, the AF-QM/MM method was used for the efficient QM calculation of protein nuclear magnetic resonance (NMR) parameters, including the chemical shift, chemical shift anisotropy tensor, and spin-spin coupling constant. In the AF-QM/MM approach, each amino acid and all the residues in its vicinity are automatically assigned as the QM region through a distance cutoff for each residue-centric QM/MM calculation. Local chemical properties of the central residue can be obtained from individual QM/MM calculations. The AF-QM/MM approach precisely reproduces the NMR chemical shifts of proteins in the gas phase from full system QM calculations. Furthermore, via the incorporation of implicit and explicit solvent models, the protein NMR chemical shifts calculated by the AF-QM/MM method are in excellent agreement with experimental values. The applications of the AF-QM/MM method may also be extended to more general biological systems such as DNA/RNA and protein-ligand complexes.
Collapse
Affiliation(s)
- Xiao He
- State
Key Laboratory of Precision Spectroscopy, Institute of Theoretical
and Computational Science, East China Normal University, Shanghai 200062, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| | - Tong Zhu
- State
Key Laboratory of Precision Spectroscopy, Institute of Theoretical
and Computational Science, East China Normal University, Shanghai 200062, China
| | - Xianwei Wang
- State
Key Laboratory of Precision Spectroscopy, Institute of Theoretical
and Computational Science, East China Normal University, Shanghai 200062, China
| | - Jinfeng Liu
- State
Key Laboratory of Precision Spectroscopy, Institute of Theoretical
and Computational Science, East China Normal University, Shanghai 200062, China
| | - John Z. H. Zhang
- State
Key Laboratory of Precision Spectroscopy, Institute of Theoretical
and Computational Science, East China Normal University, Shanghai 200062, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| |
Collapse
|
18
|
Jia X, Wang X, Liu J, Zhang JZH, Mei Y, He X. An improved fragment-based quantum mechanical method for calculation of electrostatic solvation energy of proteins. J Chem Phys 2013; 139:214104. [DOI: 10.1063/1.4833678] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
19
|
Wang X, Liu J, Zhang JZH, He X. Electrostatically embedded generalized molecular fractionation with conjugate caps method for full quantum mechanical calculation of protein energy. J Phys Chem A 2013; 117:7149-61. [PMID: 23452268 DOI: 10.1021/jp400779t] [Citation(s) in RCA: 84] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
An electrostatically embedded generalized molecular fractionation with conjugate caps (EE-GMFCC) method is developed for efficient linear-scaling quantum mechanical (QM) calculation of protein energy. This approach is based on our previously proposed GMFCC/MM method (He; et al. J. Chem. Phys. 2006, 124, 184703), In this EE-GMFCC scheme, the total energy of protein is calculated by taking a linear combination of the QM energy of the neighboring residues and the two-body QM interaction energy between non-neighboring residues that are spatially in close contact. All the fragment calculations are embedded in a field of point charges representing the remaining protein environment, which is the major improvement over our previous GMFCC/MM approach. Numerical studies are carried out to calculate the total energies of 18 real three-dimensional proteins of up to 1142 atoms using the EE-GMFCC approach at the HF/6-31G* level. The overall mean unsigned error of EE-GMFCC for the 18 proteins is 2.39 kcal/mol with reference to the full system HF/6-31G* energies. The EE-GMFCC approach is also applied for proteins at the levels of the density functional theory (DFT) and second-order many-body perturbation theory (MP2), also showing only a few kcal/mol deviation from the corresponding full system result. The EE-GMFCC method is linear-scaling with a low prefactor, trivially parallel, and can be readily applied to routinely perform structural optimization of proteins and molecular dynamics simulation with high level ab initio electronic structure theories.
Collapse
Affiliation(s)
- Xianwei Wang
- State Key Laboratory of Precision Spectroscopy and Department of Physics, Institute of Theoretical and Computational Science, East China Normal University, Shanghai 200062, China
| | | | | | | |
Collapse
|
20
|
Hua S, Li W, Li S. The Generalized Energy-Based Fragmentation Approach with an Improved Fragmentation Scheme: Benchmark Results and Illustrative Applications. Chemphyschem 2012; 14:108-15. [DOI: 10.1002/cphc.201200867] [Citation(s) in RCA: 33] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2012] [Indexed: 11/09/2022]
|
21
|
Faver JC, Zheng Z, Merz KM. Statistics-based model for basis set superposition error correction in large biomolecules. Phys Chem Chem Phys 2012; 14:7795-9. [DOI: 10.1039/c2cp23715f] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
22
|
Faver JC, Zheng Z, Merz KM. Model for the fast estimation of basis set superposition error in biomolecular systems. J Chem Phys 2011; 135:144110. [PMID: 22010701 PMCID: PMC3212865 DOI: 10.1063/1.3641894] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2011] [Accepted: 09/01/2011] [Indexed: 11/14/2022] Open
Abstract
Basis set superposition error (BSSE) is a significant contributor to errors in quantum-based energy functions, especially for large chemical systems with many molecular contacts such as folded proteins and protein-ligand complexes. While the counterpoise method has become a standard procedure for correcting intermolecular BSSE, most current approaches to correcting intramolecular BSSE are simply fragment-based analogues of the counterpoise method which require many (two times the number of fragments) additional quantum calculations in their application. We propose that magnitudes of both forms of BSSE can be quickly estimated by dividing a system into interacting fragments, estimating each fragment's contribution to the overall BSSE with a simple statistical model, and then propagating these errors throughout the entire system. Such a method requires no additional quantum calculations, but rather only an analysis of the system's interacting fragments. The method is described herein and is applied to a protein-ligand system, a small helical protein, and a set of native and decoy protein folds.
Collapse
Affiliation(s)
- John C Faver
- Quantum Theory Project, The University of Florida, 2328 New Physics Building, P.O. Box 118435, Gainesville, Florida 32611-8435, USA
| | | | | |
Collapse
|
23
|
Gordon MS, Fedorov DG, Pruitt SR, Slipchenko LV. Fragmentation Methods: A Route to Accurate Calculations on Large Systems. Chem Rev 2011; 112:632-72. [DOI: 10.1021/cr200093j] [Citation(s) in RCA: 836] [Impact Index Per Article: 64.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Mark S. Gordon
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames Iowa 50011, United States
| | - Dmitri G. Fedorov
- Nanosystem Research Institute, National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan
| | - Spencer R. Pruitt
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames Iowa 50011, United States
| | - Lyudmila V. Slipchenko
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, United States
| |
Collapse
|
24
|
Le HA, Bettens RPA. Distributed Multipoles and Energies of Flexible Molecules. J Chem Theory Comput 2011; 7:921-30. [DOI: 10.1021/ct100683u] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Hai-Anh Le
- 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
|
25
|
Jiang N, Tan RX, Ma J. Simulations of Solid-State Vibrational Circular Dichroism Spectroscopy of (S)-Alternarlactam by Using Fragmentation Quantum Chemical Calculations. J Phys Chem B 2011; 115:2801-13. [DOI: 10.1021/jp110152q] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Nan Jiang
- School of Chemistry and Chemical Engineering, Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry of MOE, Nanjing University, Nanjing, Jiangsu, 210093, People’s Republic of China
| | - Ren Xiang Tan
- Institute of Functional Biomolecules, State Key Laboratory of Pharmaceutical Biotechnology, Nanjing University, Nanjing, Jiangsu 210093, People’s Republic of China
| | - Jing Ma
- School of Chemistry and Chemical Engineering, Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry of MOE, Nanjing University, Nanjing, Jiangsu, 210093, People’s Republic of China
| |
Collapse
|
26
|
Yang Z, Hua S, Hua W, Li S. Low-Lying Structures and Stabilities of Large Water Clusters: Investigation Based on the Combination of the AMOEBA Potential and Generalized Energy-Based Fragmentation Approach. J Phys Chem A 2010; 114:9253-61. [DOI: 10.1021/jp1038267] [Citation(s) in RCA: 39] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Zhen Yang
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of the Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210093, People’s Republic of China
| | - Shugui Hua
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of the Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210093, People’s Republic of China
| | - Weijie Hua
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of the Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210093, People’s Republic of China
| | - Shuhua Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of the Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210093, People’s Republic of China
| |
Collapse
|
27
|
Hua S, Hua W, Li S. An Efficient Implementation of the Generalized Energy-Based Fragmentation Approach for General Large Molecules. J Phys Chem A 2010; 114:8126-34. [DOI: 10.1021/jp103074f] [Citation(s) in RCA: 82] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Shugui Hua
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210093, People’s Republic of China, and Department of Theoretical Chemistry, School of Biotechnology, Royal Institute of Technology, S-106 91 Stockholm, Sweden
| | - Weijie Hua
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210093, People’s Republic of China, and Department of Theoretical Chemistry, School of Biotechnology, Royal Institute of Technology, S-106 91 Stockholm, Sweden
| | - Shuhua Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210093, People’s Republic of China, and Department of Theoretical Chemistry, School of Biotechnology, Royal Institute of Technology, S-106 91 Stockholm, Sweden
| |
Collapse
|
28
|
Mullin JM, Roskop LB, Pruitt SR, Collins MA, Gordon MS. Systematic fragmentation method and the effective fragment potential: an efficient method for capturing molecular energies. J Phys Chem A 2010; 113:10040-9. [PMID: 19739681 DOI: 10.1021/jp9036183] [Citation(s) in RCA: 61] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The systematic fragmentation method fragments a large molecular system into smaller pieces, in such a way as to greatly reduce the computational cost while retaining nearly the accuracy of the parent ab initio electronic structure method. In order to attain the desired (sub-kcal/mol) accuracy, one must properly account for the nonbonded interactions between the separated fragments. Since, for a large molecular species, there can be a great many fragments and therefore a great many nonbonded interactions, computations of the nonbonded interactions can be very time-consuming. The present work explores the efficacy of employing the effective fragment potential (EFP) method to obtain the nonbonded interactions since the EFP method has been shown previously to capture nonbonded interactions with an accuracy that is often comparable to that of second-order perturbation theory. It is demonstrated that for nonbonded interactions that are not high on the repulsive wall (generally >2.7 A), the EFP method appears to be a viable approach for evaluating the nonbonded interactions. The efficacy of the EFP method for this purpose is illustrated by comparing the method to ab initio methods for small water clusters, the ZOVGAS molecule, retinal, and the alpha-helix. Using SFM with EFP for nonbonded interactions yields an error of 0.2 kcal/mol for the retinal cis-trans isomerization and a mean error of 1.0 kcal/mol for the isomerization energies of five small (120-170 atoms) alpha-helices.
Collapse
|
29
|
Abstract
The ability to perform ab initio electronic structure calculations that scales linearly with the system size is one of the central aims in theoretical chemistry. In this study, the implementation of the divide-and-conquer (DC) algorithm, an algorithm with the potential to aid the achievement of true linear scaling within Hartree-Fock (HF) theory is revisited. Standard HF calculations solve the Roothaan-Hall equations for the whole system; in the DC-HF approach, the diagonalization of the Fock matrix is carried out on smaller subsystems. The DC algorithm for HF calculations was validated on polyglycines, polyalanines and eleven real three-dimensional proteins of up to 608 atoms in this work. We also found that a fragment-based initial guess using molecular fractionation with conjugated caps (MFCC) method significantly reduces the number of SCF cycles and even is capable of achieving convergence for some globular proteins where the simple superposition of atomic densities (SAD) initial guess fails.
Collapse
Affiliation(s)
- Xiao He
- Department of Chemistry and the Quantum Theory Project, 2328 New Physics Building, P.O. Box 118435, University of Florida, Gainesville, Florida 32611-8435
| | | |
Collapse
|
30
|
Roy RK, Saha S. Studies of regioselectivity of large molecular systems using DFT based reactivity descriptors. ACTA ACUST UNITED AC 2010. [DOI: 10.1039/b811052m] [Citation(s) in RCA: 77] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
31
|
Leverentz HR, Truhlar DG. Electrostatically Embedded Many-Body Approximation for Systems of Water, Ammonia, and Sulfuric Acid and the Dependence of Its Performance on Embedding Charges. J Chem Theory Comput 2009; 5:1573-84. [DOI: 10.1021/ct900095d] [Citation(s) in RCA: 56] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Hannah R. Leverentz
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431
| | - Donald G. Truhlar
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431
| |
Collapse
|
32
|
Gordon MS, Mullin JM, Pruitt SR, Roskop LB, Slipchenko LV, Boatz JA. Accurate Methods for Large Molecular Systems. J Phys Chem B 2009; 113:9646-63. [DOI: 10.1021/jp811519x] [Citation(s) in RCA: 170] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Mark S. Gordon
- Department of Chemistry, Iowa State University, Ames, Iowa 50011, Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, and Space and Missile Propulsion Division, Air Force Research Laboratory, AFRL/RZS, 10 East Saturn Boulevard, Edward AFB, California 93524
| | - Jonathan M. Mullin
- Department of Chemistry, Iowa State University, Ames, Iowa 50011, Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, and Space and Missile Propulsion Division, Air Force Research Laboratory, AFRL/RZS, 10 East Saturn Boulevard, Edward AFB, California 93524
| | - Spencer R. Pruitt
- Department of Chemistry, Iowa State University, Ames, Iowa 50011, Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, and Space and Missile Propulsion Division, Air Force Research Laboratory, AFRL/RZS, 10 East Saturn Boulevard, Edward AFB, California 93524
| | - Luke B. Roskop
- Department of Chemistry, Iowa State University, Ames, Iowa 50011, Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, and Space and Missile Propulsion Division, Air Force Research Laboratory, AFRL/RZS, 10 East Saturn Boulevard, Edward AFB, California 93524
| | - Lyudmila V. Slipchenko
- Department of Chemistry, Iowa State University, Ames, Iowa 50011, Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, and Space and Missile Propulsion Division, Air Force Research Laboratory, AFRL/RZS, 10 East Saturn Boulevard, Edward AFB, California 93524
| | - Jerry A. Boatz
- Department of Chemistry, Iowa State University, Ames, Iowa 50011, Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, and Space and Missile Propulsion Division, Air Force Research Laboratory, AFRL/RZS, 10 East Saturn Boulevard, Edward AFB, California 93524
| |
Collapse
|
33
|
Hua W, Fang T, Li W, Yu JG, Li S. Geometry Optimizations and Vibrational Spectra of Large Molecules from a Generalized Energy-Based Fragmentation Approach. J Phys Chem A 2008; 112:10864-72. [DOI: 10.1021/jp8026385] [Citation(s) in RCA: 112] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Weijie Hua
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210093, P. R. China, and College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
| | - Tao Fang
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210093, P. R. China, and College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
| | - Wei Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210093, P. R. China, and College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
| | - Jian-Guo Yu
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210093, P. R. China, and College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
| | - Shuhua Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210093, P. R. China, and College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China
| |
Collapse
|
34
|
Jiang N, Ma J. Conformational Simulations of Aqueous Solvated α-Conotoxin GI and Its Single Disulfide Analogues Using a Polarizable Force Field Model. J Phys Chem A 2008; 112:9854-67. [DOI: 10.1021/jp8029693] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Affiliation(s)
- Nan Jiang
- School of Chemistry and Chemical Engineering, Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry of MOE, Nanjing University, Nanjing, 210093, People’s Republic of China
| | - Jing Ma
- School of Chemistry and Chemical Engineering, Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry of MOE, Nanjing University, Nanjing, 210093, People’s Republic of China
| |
Collapse
|
35
|
Li H, Li W, Li S, Ma J. Fragmentation-Based QM/MM Simulations: Length Dependence of Chain Dynamics and Hydrogen Bonding of Polyethylene Oxide and Polyethylene in Aqueous Solutions. J Phys Chem B 2008; 112:7061-70. [DOI: 10.1021/jp800777e] [Citation(s) in RCA: 27] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Affiliation(s)
- Hui Li
- School of Chemistry and Chemical Engineering, Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry of MOE, Nanjing University, Nanjing 210093, P. R. China
| | - Wei Li
- School of Chemistry and Chemical Engineering, Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry of MOE, Nanjing University, Nanjing 210093, P. R. China
| | - Shuhua Li
- School of Chemistry and Chemical Engineering, Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry of MOE, Nanjing University, Nanjing 210093, P. R. China
| | - Jing Ma
- School of Chemistry and Chemical Engineering, Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry of MOE, Nanjing University, Nanjing 210093, P. R. China
| |
Collapse
|
36
|
Hatfield MPD, Palermo NY, Csontos J, Murphy RF, Lovas S. Evaluation of methods to cap molecular fragments in calculating energies of interaction in avian pancreatic polypeptide. INTERNATIONAL JOURNAL OF QUANTUM CHEMISTRY 2008; 108:1017-1021. [PMID: 18985167 PMCID: PMC2577377 DOI: 10.1002/qua.21553] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/27/2023]
Abstract
The accuracy of the determination of the energy of interaction between Phe20 and the Pro5-Thr6-Tyr7-Pro8 complex inside the hydrophobic core of avian pancreatic polypeptide was investigated using three capping strategies for molecular fractionation with conjugated caps and DFT quantum chemical calculations at the BHandHLYP/cc-pVTZ level of theory. The most accurate determination resulted from acetylation of the alpha-amino group combined with methyl amidation of the alpha-carbonyl group with relative deviations less than 10%. Combinations of hydrogenation of the alpha-amino group with the replacement of the alpha-carbonyl group with a hydrogen and the hydrogenation of the alpha-amino group with methylation of the alpha-carbonyl group were less accurate, leading to relative deviations up to 35%. Choice of capping methods depends on the structural features of the polypeptide system, the desired accuracy and the available computational resources.
Collapse
Affiliation(s)
- Marcus P. D. Hatfield
- Nebraska Wesleyan University, Lincoln, NE 68504, U.S.A
- Department of Biomedical Sciences, Creighton University, Omaha, NE 68178, U.S.A
| | - Nicholas Y. Palermo
- Department of Biomedical Sciences, Creighton University, Omaha, NE 68178, U.S.A
| | - József Csontos
- Department of Biomedical Sciences, Creighton University, Omaha, NE 68178, U.S.A
| | - Richard F. Murphy
- Department of Biomedical Sciences, Creighton University, Omaha, NE 68178, U.S.A
| | - Sándor Lovas
- Department of Biomedical Sciences, Creighton University, Omaha, NE 68178, U.S.A
| |
Collapse
|
37
|
Li S, Li W. Fragment energy approach to Hartree–Fock calculations of macromolecules. ACTA ACUST UNITED AC 2008. [DOI: 10.1039/b703896h] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
38
|
Abstract
The nuclear magnetic shielding tensor is a molecular property that can be computed from first principles. In this work we show that by utilizing the fragmentation approach, one is able to accurately compute this property for a large class of molecules. This is of great significance because the computational expense required in the evaluation of the shielding tensor for all nuclei in a large molecule is now subject to near linear scaling. On the basis of previous studies and this work, it is also very likely that all molecular properties that can be expressed as derivatives of the total energy of the system are also amenable to accurate evaluation via fragmentation. If only the chemical shifts for nuclei in a small part of a large molecule are of interest, then only those molecular fragments containing those nuclei need to have their shielding tensors evaluated. Further, the fragmentation approach allows one to construct a database of molecular fragments that could, in principle, be used in the NMR characterization of molecules and at the same time provide possible three-dimensional representations of these molecules.
Collapse
Affiliation(s)
- Adrian M Lee
- Department of Chemistry, National University of Singapore, 3 Science Drive 3, Singapore 117543
| | | |
Collapse
|
39
|
Li W, Li S, Jiang Y. Generalized Energy-Based Fragmentation Approach for Computing the Ground-State Energies and Properties of Large Molecules. J Phys Chem A 2007; 111:2193-9. [PMID: 17388268 DOI: 10.1021/jp067721q] [Citation(s) in RCA: 221] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We present a generalized energy-based fragmentation (GEBF) approach for approximately predicting the ground-state energies and molecular properties of large molecules, especially those charged and polar molecules. In this approach, the total energy (or properties) of a large molecule can be approximately obtained from energy (or properties) calculations on various small subsystems, each of which is constructed to contain a certain fragment and its local surroundings within a given distance. In the quantum chemistry calculation of a given subsystem, those distant atoms (outside this subsystem) are modeled as background point charges at the corresponding nuclear centers. This treatment allows long-range electrostatic interaction and polarization effects between distant fragments to be taken into account approximately, which are very important for polar and charged molecules. We also propose a new fragmentation scheme for constructing subsystems. Our test calculations at the Hartree-Fock and second-order Møller-Plesser perturbation theory levels demonstrate that the approach could yield satisfactory ground-state energies, the dipole moments, and static polarizabilities for polar and charged molecules such as water clusters and proteins.
Collapse
Affiliation(s)
- Wei Li
- School of Chemistry and Chemical Engineering, Institute of Theoretical and Computational Chemistry, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Nanjing University, Nanjing, 210093, P. R. China.
| | | | | |
Collapse
|