1
|
Gupta AK, Maier S, Thapa B, Raghavachari K. Toward Post-Hartree-Fock Accuracy for Protein-Ligand Affinities Using the Molecules-in-Molecules Fragmentation-Based Method. J Chem Theory Comput 2024; 20:2774-2785. [PMID: 38530869 DOI: 10.1021/acs.jctc.3c01293] [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: 03/28/2024]
Abstract
The complexity and size of large molecular systems, such as protein-ligand complexes, pose computational challenges for accurate post-Hartree-Fock calculations. This study delivers a thorough benchmarking of the Molecules-in-Molecules (MIM) method, presenting a clear and accessible strategy for layer/theory selections in post-Hartree-Fock computations on substantial molecular systems, notably protein-ligand complexes. An approach is articulated, enabling augmented computational efficiency by strategically canceling out common subsystem energy terms between complexes and proteins within the supermolecular equation. Employing DLPNO-based post-Hartree-Fock methods in conjunction with the three-layer MIM method (MIM3), this study demonstrates the achievement of protein-ligand binding energies with remarkable accuracy (errors <1 kcal mol-1), while significantly reducing computational costs. Furthermore, noteworthy correlations between theoretically computed interaction energies and their experimental equivalents were observed, with R2 values of approximately 0.90 and 0.78 for CDK2 and BZT-ITK sets, respectively, thus validating the efficacy of the MIM method in calculating binding energies. By highlighting the crucial role of diffuse or small Pople-style basis sets in the middle layer for reducing energy errors, this work provides valuable insights and practical methodologies for interaction energy computations in large molecular complexes and opens avenues for their application across a diverse range of molecular systems.
Collapse
Affiliation(s)
- Ankur K Gupta
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Sarah Maier
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Bishnu Thapa
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Krishnan Raghavachari
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| |
Collapse
|
2
|
Pereira RW, Ramabhadran RO. Accurate Computation of Aqueous p Kas of Biologically Relevant Organic Acids: Overcoming the Challenges Posed by Multiple Conformers, Tautomeric Equilibria, and Disparate Functional Groups with the Fully Black-Box p K-Yay Method. J Phys Chem A 2023; 127:9121-9138. [PMID: 37862610 DOI: 10.1021/acs.jpca.3c02977] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2023]
Abstract
The use of static electronic structure calculations to compute solution-phase pKas offers a great advantage in that a macroscopic bulk property could be computed via microscopic computations involving very few molecules. There are various sources of errors in the quantum chemical calculations though. Overcoming these errors to accurately compute pKas of a plethora of acids is an active area of research in physical chemistry pursued by both computational as well as experimental chemists. We recently developed the pK-Yay method in our attempt to accurately compute aqueous pKas of strong and weak acids. The method is fully black-box, computationally inexpensive, and is very easy for even a nonexpert to use. However, the method was thus far tested on very few molecules (only 16 in all). Herein, in order to assess the future applicability of pK-Yay, we study the effect of multiple conformers, the presence of tautomers under equilibrium, and the impact of a wide variety of functional groups (derivatives of acetic acid with substituents at various positions, dicarboxylic acids, aromatic carboxylic acids, amines and amides, phenols and thiols, and fluorine bearing organic acids). Starting with more than 1000 conformers and tautomers, this study establishes that overall errors of ∼ 1.0 pKa units are routinely obtained for a majority of the molecules. Larger errors are noted in cases where multiple charges, intramolecular hydrogen bonding, and several ionizable functional groups are simultaneously present. An important conclusion to emerge from this work is that, the computed pKas are insensitive (difference <0.5) to whether we consider multiple conformers/tautomers or only choose the most stable conformer/tautomer. Further, pK-Yay captures the stereoelectronic effects arising due to differing axial vs equatorial pattern, and is useful to predict the dominant acid-base equilibrium in a system featuring several equilibria. Overall, pK-Yay may be employed in several chemical applications featuring organic molecules and biomonomers.
Collapse
Affiliation(s)
- Roshni W Pereira
- Department of Chemistry, Indian Institute of Science Education and Research (IISER) Tirupati, Andhra Pradesh 517507, India
- Centre for Atomic Molecular Optical Sciences and Technology (CAMOST), Tirupati, Andhra Pradesh 517507, India
| | - Raghunath O Ramabhadran
- Department of Chemistry, Indian Institute of Science Education and Research (IISER) Tirupati, Andhra Pradesh 517507, India
- Centre for Atomic Molecular Optical Sciences and Technology (CAMOST), Tirupati, Andhra Pradesh 517507, India
| |
Collapse
|
3
|
Olive LN, Dornshuld EV, Schaefer HF, Tschumper GS. Competition between Solvent···Solvent and Solvent···Solute Interactions in the Microhydration of the Tetrafluoroborate Anion, BF 4-(H 2O) n=1,2,3,4. J Phys Chem A 2023; 127:8806-8820. [PMID: 37774368 DOI: 10.1021/acs.jpca.3c04014] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/01/2023]
Abstract
This study systematically examines the interactions of the tetrafluoroborate anion (BF4-) with up to four water molecules (BF4-(H2O)n=1,2,3,4). Full geometry optimizations and subsequent harmonic vibrational frequency computations are performed using a variety of density functional theory (DFT) methods (B3LYP, B3LYP-D3BJ, and M06-2X) and the MP2 ab initio method with a triple-ζ correlation consistent basis set augmented with diffuse functions on all non-hydrogen atoms (cc-pVTZ for H and aug-cc-pVTZ for B, O, and F; denoted as haTZ). Optimized structures and harmonic vibrational frequencies were also obtained with the CCSD(T) ab initio method and the haTZ basis set for the mono- and dihydrate (n = 1, 2) structures. The 2-body:Many-body (2b:Mb) technique, in which CCSD(T) computations capture the 1- and 2-body contributions to the interactions and MP2 computations recover all higher-order contributions, was used to extend these demanding computations to the tri- and tetrahydrate (n = 3, 4) systems. Four, five, and eight new stationary points have been identified for the di-, tri-, and tetrahydrate systems, respectively. The global minimum of the monohydrate adopts a symmetric double ionic hydrogen bond motif with C2v symmetry and an electronic dissociation energy of 13.17 kcal mol-1 at the CCSD(T)/haTZ level of theory. This strong solvent···solute interaction, however, competes with solute···solute interactions in the lowest-energy BF4-(H2O)n=2,3,4 minima that are not seen in the other di-, tri-, or tetrahydrate minima. The latter interactions help increase the 2b:Mb dissociation energies to more than 26, 41, and 51 kcal mol-1 for n = 2, 3, and 4, respectively. Structures that form hydrogen bonds between the solvating water molecules also exhibit the largest shifts in the harmonic OH stretching frequencies for the waters of hydration. These shifts can exceed -280 cm-1 relative to an isolated H2O molecule at the 2b:Mb/haTZ level of theory.
Collapse
Affiliation(s)
- Laura N Olive
- Center for Computational Quantum Chemistry, University of Georgia, Athens, Georgia 30602, United States
| | - Eric V Dornshuld
- Department of Chemistry, Mississippi State University, Mississippi State, Mississippi 39762, United States
| | - Henry F Schaefer
- Center for Computational Quantum Chemistry, University of Georgia, Athens, Georgia 30602, United States
| | - Gregory S Tschumper
- Department of Chemistry and Biochemistry, University of Mississippi, University, Mississippi 38677, United States
| |
Collapse
|
4
|
Tripathy V, Raghavachari K. Fragment-based models for dissociation of strong acids in water: Electrostatic embedding minimizes the dependence on the fragmentation schemes. J Chem Phys 2023; 159:124106. [PMID: 38127382 DOI: 10.1063/5.0164089] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2023] [Accepted: 08/28/2023] [Indexed: 12/23/2023] Open
Abstract
Fragmentation methods such as MIM (Molecules-in-Molecules) provide a route to accurately model large systems and have been successful in predicting their structures, energies, and spectroscopic properties. However, their use is often limited to systems at equilibrium due to the inherent complications in the choice of fragments in systems away from equilibrium. Furthermore, the presence of charges resulting from any heterolytic bond breaking may increase the fragmentation error. We have previously suggested EE-MIM (Electrostatically Embedded Molecules-In-Molecules) as a method to mitigate the errors resulting from the missing long-range interactions in molecular clusters in equilibrium. Here, we show that the same method can be applied to improve the performance of MIM to solve the longstanding problem of dependency of the fragmentation energy error on the choice of the fragmentation scheme. We chose four widely used acid dissociation reactions (HCl, HClO4, HNO3, and H2SO4) as test cases due to their importance in chemical processes and complex reaction potential energy surfaces. Electrostatic embedding improves the performance at both one and two-layer MIM as shown by lower EE-MIM1 and EE-MIM2 errors. The EE-MIM errors are also demonstrated to be less dependent on the choice of the fragmentation scheme by analyzing the variation in fragmentation energy at the points with more than one possible fragmentation scheme (points where the fragmentation scheme changes). EE-MIM2 with M06-2X as the low-level resulted in a variation of less than 1 kcal/mol for all the cases and 1 kJ/mol for all but three cases, rendering our method fragmentation scheme-independent for acid dissociation processes.
Collapse
Affiliation(s)
- Vikrant Tripathy
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, USA
| | | |
Collapse
|
5
|
Sahu N, Khire SS, Gadre SR. Combining fragmentation method and high-performance computing: Geometry optimization and vibrational spectra of proteins. J Chem Phys 2023; 159:044309. [PMID: 37522406 DOI: 10.1063/5.0149572] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2023] [Accepted: 07/12/2023] [Indexed: 08/01/2023] Open
Abstract
Exploring the structures and spectral features of proteins with advanced quantum chemical methods is an uphill task. In this work, a fragment-based molecular tailoring approach (MTA) is appraised for the CAM-B3LYP/aug-cc-pVDZ-level geometry optimization and vibrational infrared (IR) spectra calculation of ten real proteins containing up to 407 atoms and 6617 basis functions. The use of MTA and the inherently parallel nature of the fragment calculations enables a rapid and accurate calculation of the IR spectrum. The applicability of MTA to optimize the protein geometry and evaluate its IR spectrum employing a polarizable continuum model with water as a solvent is also showcased. The typical errors in the total energy and IR frequencies computed by MTA vis-à-vis their full calculation (FC) counterparts for the studied protein are 5-10 millihartrees and 5 cm-1, respectively. Moreover, due to the independent execution of the fragments, large-scale parallelization can also be achieved. With increasing size and level of theory, MTA shows an appreciable advantage in computer time as well as memory and disk space requirement over the corresponding FCs. The present study suggests that the geometry optimization and IR computations on the biomolecules containing ∼1000 atoms and/or ∼15 000 basis functions using MTA and HPC facility can be clearly envisioned in the near future.
Collapse
Affiliation(s)
- Nityananda Sahu
- Theoretische Chemie, Philipps-Universität Marburg, 35032 Marburg, Germany
| | - Subodh S Khire
- RIKEN Center for Computational Science, Kobe 650-0047, Japan
| | - Shridhar R Gadre
- Departments of Scientific Computing, Modelling & Simulation and Chemistry, Savitribai Phule Pune University, Pune 411007, India
| |
Collapse
|
6
|
Bowling PE, Broderick DR, Herbert JM. Fragment-Based Calculations of Enzymatic Thermochemistry Require Dielectric Boundary Conditions. J Phys Chem Lett 2023; 14:3826-3834. [PMID: 37061921 DOI: 10.1021/acs.jpclett.3c00533] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/19/2023]
Abstract
Electronic structure calculations on enzymes require hundreds of atoms to obtain converged results, but fragment-based approximations offer a cost-effective solution. We present calculations on enzyme models containing 500-600 atoms using the many-body expansion, comparing to benchmarks in which the entire enzyme-substrate complex is described at the same level of density functional theory. When the amino acid fragments contain ionic side chains, the many-body expansion oscillates under vacuum boundary conditions but rapid convergence is restored using low-dielectric boundary conditions. This implies that full-system calculations in the gas phase are inappropriate benchmarks for assessing errors in fragment-based approximations. A three-body protocol retains sub-kilocalorie per mole fidelity with respect to a supersystem calculation, as does a two-body calculation combined with a full-system correction at a low-cost level of theory. These protocols pave the way for application of high-level quantum chemistry to large systems via rigorous, ab initio treatment of many-body polarization.
Collapse
Affiliation(s)
- Paige E Bowling
- Biophysics Graduate Program, The Ohio State University, Columbus, Ohio 43210, United States
- Department of Chemistry & Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| | - Dustin R Broderick
- Department of Chemistry & Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| | - John M Herbert
- Biophysics Graduate Program, The Ohio State University, Columbus, Ohio 43210, United States
- Department of Chemistry & Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| |
Collapse
|
7
|
Gupta AK, Gamoke BC, Raghavachari K. Interaction-Deletion: A Composite Energy Method for the Optimization of Molecular Systems Selectively Removing Specific Nonbonded Interactions. J Phys Chem A 2021; 125:4668-4682. [PMID: 34014684 DOI: 10.1021/acs.jpca.1c02918] [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
The complex interactions between different portions of a large molecule can be challenging to analyze through traditional electronic structure calculations. Moreover, standard methods cannot easily quantify the physical consequences of individual pairwise interactions inside a molecule. By creating a set of molecular fragments, we propose a composite energy method to explore changes in a molecule caused by removing selected nonbonded interactions between different molecular portions. Energies and forces are easily obtained with this composite approach, allowing geometry optimizations that lead to chemically meaningful structures that describe how the omitted interactions contribute to changes in the local geometrical minima. We illustrate the application of our new hybrid scheme by computing the influence of intramolecular hydrogen-bonding interactions in two small molecules: 1,6-(tG+G+TG+G+g-)-hexanediol and a cyclic analogue, cis-1,4-cyclohexanediol. The resulting structural and energetic changes are interpreted to yield key physical insights and quantify concepts such as "preparation energy" or "reorganization energy". We demonstrate that the composite method can be extended to larger molecular systems by showing its application on a Si(100) surface model containing interactions between dissociated ammonia molecules on adjacent surface dimers. The scheme's efficacy is also tested by applying it to systems having multiple intramolecular interactions, viz., 310-polyglycine and H+GPGG. Furthermore, the cooperative nature of intramolecular hydrogen bonds is explored by using interaction-deletion in 2-nitrobenzene-1,3-diol.
Collapse
Affiliation(s)
- Ankur K Gupta
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Benjamin C Gamoke
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Krishnan Raghavachari
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| |
Collapse
|
8
|
Noffke BW, Beckett D, Li LS, Raghavachari K. Aromatic Fragmentation Based on a Ring Overlap Scheme: An Algorithm for Large Polycyclic Aromatic Hydrocarbons Using the Molecules-in-Molecules Fragmentation-Based Method. J Chem Theory Comput 2020; 16:2160-2171. [DOI: 10.1021/acs.jctc.9b00566] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Affiliation(s)
- Benjamin W. Noffke
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Daniel Beckett
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Liang-shi Li
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Krishnan Raghavachari
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| |
Collapse
|
9
|
Liu KY, Herbert JM. Energy-Screened Many-Body Expansion: A Practical Yet Accurate Fragmentation Method for Quantum Chemistry. J Chem Theory Comput 2019; 16:475-487. [PMID: 31765559 DOI: 10.1021/acs.jctc.9b01095] [Citation(s) in RCA: 41] [Impact Index Per Article: 8.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
We introduce an implementation of the truncated many-body expansion, MBE(n), in which the n-body corrections are screened using the effective fragment potential force field, and only those that exceed a specified energy threshold are computed at a quantum-mechanical level of theory. This energy-screened MBE(n) approach is tested at the n = 3 level for a sequence of water clusters, (H2O)N=6-34. A threshold of 0.25 kJ/mol eliminates more than 80% of the subsystem electronic structure calculations and is even more efficacious in that respect than is distance-based screening. Even so, the energy-screened MBE(3) method is faithful to a full-system quantum chemistry calculation to within 1-2 kJ/mol/monomer, even in good quality basis sets such as aug-cc-pVTZ. These errors can be reduced by means of a two-layer approach that involves a Hartree-Fock calculation for the entire cluster. Such a correction proves to be necessary in order to obtain accurate relative energies for conformational isomers of (H2O)20, but the cost of a full-system Hartree-Fock calculation remains smaller than the cost of three-body subsystem calculations at correlated levels of theory. At the level of second-order Møller-Plesset perturbation theory (MP2), a screened MBE(3) calculation plus a full-system Hartree-Fock calculation is less expensive than a full-system MP2 calculation starting at N = 12 water molecules. This is true even if all MBE(3) subsystem calculations are performed on a single 40-core compute node, i.e., without significant parallelization. Energy-screened MBE(n) thus provides a fragment-based method that is accurate, stable in large basis sets, and low in cost, even when the latter is measured in aggregate computer time.
Collapse
Affiliation(s)
- Kuan-Yu 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
|
10
|
Sadhukhan T, Beckett D, Thapa B, Raghavachari K. Coupling Constants, High Spin, and Broken Symmetry States of Organic Radicals: an Assessment of the Molecules-in-Molecules Fragmentation-Based Method. J Chem Theory Comput 2019; 15:5998-6009. [DOI: 10.1021/acs.jctc.9b00563] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Affiliation(s)
- Tumpa Sadhukhan
- Department of Chemistry, Indiana University, Bloomington, 47405 Indiana, United States
| | - Daniel Beckett
- Department of Chemistry, Indiana University, Bloomington, 47405 Indiana, United States
| | - Bishnu Thapa
- Department of Chemistry, Indiana University, Bloomington, 47405 Indiana, United States
| | - Krishnan Raghavachari
- Department of Chemistry, Indiana University, Bloomington, 47405 Indiana, United States
| |
Collapse
|
11
|
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
|
12
|
Kaliakin DS, Fedorov DG, Alexeev Y, Varganov SA. Locating Minimum Energy Crossings of Different Spin States Using the Fragment Molecular Orbital Method. J Chem Theory Comput 2019; 15:6074-6084. [DOI: 10.1021/acs.jctc.9b00641] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Affiliation(s)
- Danil S. Kaliakin
- Department of Chemistry, University of Nevada, Reno, 1664 N. Virginia Street, Reno, Nevada 89557-0216, United States
| | - Dmitri G. Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba 305-8568, Japan
| | - Yuri Alexeev
- Computational Science Division and Argonne Leadership Computing Facility, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - Sergey A. Varganov
- Department of Chemistry, University of Nevada, Reno, 1664 N. Virginia Street, Reno, Nevada 89557-0216, United States
| |
Collapse
|
13
|
Fedorov DG. Solvent Screening in Zwitterions Analyzed with the Fragment Molecular Orbital Method. J Chem Theory Comput 2019; 15:5404-5416. [PMID: 31461277 DOI: 10.1021/acs.jctc.9b00715] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Based on induced solvent charges, a new model of solvent screening is developed in the framework of the fragment molecular orbital combined with the polarizable continuum model. The developed model is applied to analyze interactions in a prototypical zwitterionic system, sodium chloride in water, and it is shown that the large underestimation of the interaction in the original solvent screening based on local charges is successfully corrected. The model is also applied to a complex of the Trp-cage (PDB: 1L2Y ) miniprotein with an anionic ligand, and the physical factors determined protein-ligand binding in solution are unraveled.
Collapse
Affiliation(s)
- Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat) , National Institute of Advanced Industrial Science and Technology (AIST) , Central 2, Umezono 1-1-1 , Tsukuba 305-8568 , Japan
| |
Collapse
|
14
|
Fedorov DG, Brekhov A, Mironov V, Alexeev Y. Molecular Electrostatic Potential and Electron Density of Large Systems in Solution Computed with the Fragment Molecular Orbital Method. J Phys Chem A 2019; 123:6281-6290. [DOI: 10.1021/acs.jpca.9b04936] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Dmitri G. Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba, 305-8568, Japan
| | - Anton Brekhov
- Department of Chemistry, Lomonosov Moscow State University, Moscow, 119991, Russian Federation
| | - Vladimir Mironov
- Department of Chemistry, Lomonosov Moscow State University, Moscow, 119991, Russian Federation
| | - Yuri Alexeev
- Argonne Leadership Computing Facility and Computational Science Division, Argonne National Laboratory, Argonne, Illinois, 60439, United States
| |
Collapse
|
15
|
Gupta AK, Thapa B, Raghavachari K. Exploring Reaction Energy Profiles Using the Molecules-in-Molecules Fragmentation-Based Approach. J Chem Theory Comput 2019; 15:3991-4002. [PMID: 31181886 DOI: 10.1021/acs.jctc.9b00152] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The Molecules-in-Molecules (MIM) fragmentation-based approach has been successfully used in previous studies to obtain the energies, optimized geometries, and spectroscopic properties of large molecular systems. The present work delineates a protocol to study the potential energy profiles for multistep chemical reactions using the MIM methodology. In a complex multistep chemical reaction, the fragmentation scheme needs to be changed as the reacting species transition into a new reaction step, resulting in a discontinuity in the potential energy curve of the reaction. In our approach, the fragmentation scheme for a particular step in a reaction is chosen on the basis of the nature of the bonding changes associated with that step. Thus, the reactant, transition state, and product are treated consistently throughout the reaction step, leading to an accurate energy barrier for that step. The discontinuity now occurs in describing the energies of reaction intermediates at the transition point between two reaction steps that are treated by two different fragmentation schemes. To address this issue, we propose a systematic procedure for obtaining continuous potential energy curves that are least shifted from their initial positions. The corrected MIM potential energy curves are continuous with activation energies preserved. Following this approach, energy profiles of complex reactions involving large molecular species can be obtained at high levels of theory with a reasonable computational cost.
Collapse
Affiliation(s)
- Ankur Kumar Gupta
- Department of Chemistry , Indiana University , Bloomington , Indiana 47405 , United States
| | - Bishnu Thapa
- Department of Chemistry , Indiana University , Bloomington , Indiana 47405 , United States
| | - Krishnan Raghavachari
- Department of Chemistry , Indiana University , Bloomington , Indiana 47405 , United States
| |
Collapse
|
16
|
Sexton TM, Tschumper GS. 2-body:Many-body QM:QM study of structures, energetics, and vibrational frequencies for microhydrated halide ions. Mol Phys 2018. [DOI: 10.1080/00268976.2018.1554827] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Affiliation(s)
- Thomas More Sexton
- Department of Chemistry and Biochemistry, University of Mississippi, University, MS, USA
| | - Gregory S. Tschumper
- Department of Chemistry and Biochemistry, University of Mississippi, University, MS, USA
| |
Collapse
|
17
|
Debnath S, Sengupta A, Jose KVJ, Raghavachari K. Fragment-Based Approaches for Supramolecular Interaction Energies: Applications to Foldamers and Their Complexes with Anions. J Chem Theory Comput 2018; 14:6226-6239. [DOI: 10.1021/acs.jctc.8b00525] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Sibali Debnath
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Arkajyoti Sengupta
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - K. V. Jovan Jose
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Krishnan Raghavachari
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| |
Collapse
|
18
|
Nakata H, Fedorov DG. Analytic second derivatives for the efficient electrostatic embedding in the fragment molecular orbital method. J Comput Chem 2018; 39:2039-2050. [DOI: 10.1002/jcc.25360] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2018] [Revised: 04/27/2018] [Accepted: 04/29/2018] [Indexed: 01/09/2023]
Affiliation(s)
- Hiroya Nakata
- Department of Fundamental Technology Research; Research and Development Center Kagoshima, Kyocera, 1-4 Kokubu Yamashita-cho; Kirishima-shi Kagoshima, 899-4312 Japan
| | - Dmitri G. Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono; Tsukuba Ibaraki, 305-8568 Japan
| |
Collapse
|
19
|
Empirical corrections and pair interaction energies in the fragment molecular orbital method. Chem Phys Lett 2018. [DOI: 10.1016/j.cplett.2018.06.025] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
|
20
|
Fedorov DG. Analysis of solute-solvent interactions using the solvation model density combined with the fragment molecular orbital method. Chem Phys Lett 2018. [DOI: 10.1016/j.cplett.2018.05.002] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
|
21
|
Thapa B, Beckett D, Jovan Jose KV, Raghavachari K. Assessment of Fragmentation Strategies for Large Proteins Using the Multilayer Molecules-in-Molecules Approach. J Chem Theory Comput 2018; 14:1383-1394. [DOI: 10.1021/acs.jctc.7b01198] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Affiliation(s)
- Bishnu Thapa
- Department of Chemistry, Indiana University, Bloomington 47405, Indiana, United States
| | - Daniel Beckett
- Department of Chemistry, Indiana University, Bloomington 47405, Indiana, United States
| | - K. V. Jovan Jose
- Department of Chemistry, Indiana University, Bloomington 47405, Indiana, United States
| | - Krishnan Raghavachari
- Department of Chemistry, Indiana University, Bloomington 47405, Indiana, United States
| |
Collapse
|
22
|
Nishimoto Y, Fedorov DG. Adaptive frozen orbital treatment for the fragment molecular orbital method combined with density-functional tight-binding. J Chem Phys 2018; 148:064115. [DOI: 10.1063/1.5012935] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Affiliation(s)
- Yoshio Nishimoto
- Fukui Institute for Fundamental Chemistry, Kyoto University, 34-4 Takano Nishihiraki-cho, Sakyo-ku, Kyoto 606-8103, Japan
| | - Dmitri G. Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan
| |
Collapse
|
23
|
Fedorov DG, Kitaura K. Pair Interaction Energy Decomposition Analysis for Density Functional Theory and Density-Functional Tight-Binding with an Evaluation of Energy Fluctuations in Molecular Dynamics. J Phys Chem A 2018; 122:1781-1795. [DOI: 10.1021/acs.jpca.7b12000] [Citation(s) in RCA: 37] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Dmitri G. Fedorov
- Research
Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba 305-8568, Japan
| | - Kazuo Kitaura
- Advanced
Institute for Computational Science (AICS), RIKEN, 7-1-26 Minatojima-Minami-Machi,
Chuo-ku, Kobe, Hyogo 650-0047, Japan
- Fukui
Institute for Fundamental Chemistry, Kyoto University, Takano-Nishihiraki-cho
34-4, Sakyou-ku, Kyoto 606-8103, Japan
| |
Collapse
|
24
|
Jacobs AD, Jovan Jose KV, Horness R, Raghavachari K, Thielges MC, Clemmer DE. Cooperative Formation of Icosahedral Proline Clusters from Dimers. JOURNAL OF THE AMERICAN SOCIETY FOR MASS SPECTROMETRY 2018; 29:95-102. [PMID: 29127569 PMCID: PMC6884317 DOI: 10.1007/s13361-017-1833-6] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/29/2017] [Revised: 09/28/2017] [Accepted: 10/01/2017] [Indexed: 06/07/2023]
Abstract
Ion mobility spectrometry-mass spectrometry and Fourier transform infrared spectroscopy (FTIR) techniques were combined with quantum chemical calculations to examine the origin of icosahedral clusters of the amino acid proline. When enantiopure proline solutions are electrosprayed (using nanospray) from 100 mM ammonium acetate, only three peaks are observed in the mass spectrum across a concentration range of five orders of magnitude: a monomer [Pro+H]+ species, favored from 0.001 to 0.01 mM proline concentrations; a dimer [2Pro+H]+ species, the most abundant species for proline concentrations above 0.01 mM; and, the dimer and dodecamer [12Pro+2H]2+ for 1.0 mM and more concentrated proline solutions. Electrospraying racemic D/L-proline solutions from 100 mM ammonium acetate leads to a monomer at low proline concentrations (0.001 to 0.1 mM), and a dimer at higher concentrations (>0.09 mM), as well as a very small population of 8 to 15 Pro clusters that comprise <0.1% of the total ion signals even at the highest proline concentration. Solution FTIR studies show unique features that increase in intensity in the enantiopure proline solutions, consistent with clustering, presumably from the icosahedral geometry in bulk solution. When normalized for the total proline, these results are indicative of a cooperative formation of the enantiopure 12Pro species from 2Pro. Graphical Abstract.
Collapse
Affiliation(s)
- Alexander D Jacobs
- Department of Chemistry, Indiana University, Bloomington, IN, 47405, USA
| | - K V Jovan Jose
- Department of Chemistry, Indiana University, Bloomington, IN, 47405, USA
| | - Rachel Horness
- Department of Chemistry, Indiana University, Bloomington, IN, 47405, USA
| | | | - Megan C Thielges
- Department of Chemistry, Indiana University, Bloomington, IN, 47405, USA
| | - David E Clemmer
- Department of Chemistry, Indiana University, Bloomington, IN, 47405, USA.
| |
Collapse
|
25
|
|
26
|
Jose KVJ, Raghavachari K. Fragment-Based Approach for the Evaluation of NMR Chemical Shifts for Large Biomolecules Incorporating the Effects of the Solvent Environment. J Chem Theory Comput 2017; 13:1147-1158. [DOI: 10.1021/acs.jctc.6b00922] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Affiliation(s)
- K. V. Jovan Jose
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Krishnan Raghavachari
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| |
Collapse
|
27
|
Sahu N, Gadre SR. Vibrational infrared and Raman spectra of polypeptides: Fragments-in-fragments within molecular tailoring approach. J Chem Phys 2017; 144:114113. [PMID: 27004868 DOI: 10.1063/1.4943966] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Abstract
The present work reports the calculation of vibrational infrared (IR) and Raman spectra of large molecular systems employing molecular tailoring approach (MTA). Further, it extends the grafting procedure for the accurate evaluation of IR and Raman spectra of large molecular systems, employing a new methodology termed as Fragments-in-Fragments (FIF), within MTA. Unlike the previous MTA-based studies, the accurate estimation of the requisite molecular properties is achieved without performing any full calculations (FC). The basic idea of the grafting procedure is implemented by invoking the nearly basis-set-independent nature of the MTA-based error vis-à-vis the respective FCs. FIF has been tested out for the estimation of the above molecular properties for three isomers, viz., β-strand, 310- and α-helix of acetyl(alanine)nNH2 (n = 10, 15) polypeptides, three conformers of doubly protonated gramicidin S decapeptide and trpzip2 protein (PDB id: 1LE1), respectively, employing BP86/TZVP, M06/6-311G**, and M05-2X/6-31G** levels of theory. For most of the cases, a maximum difference of 3 cm(-1) is achieved between the grafted-MTA frequencies and the corresponding FC values. Further, a comparison of the BP86/TZVP level IR and Raman spectra of α-helical (alanine)20 and its N-deuterated derivative shows an excellent agreement with the existing experimental spectra. In view of the requirement of only MTA-based calculations and the ability of FIF to work at any level of theory, the current methodology provides a cost-effective solution for obtaining accurate spectra of large molecular systems.
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
|
28
|
Jose KVJ, Raghavachari K. Molecules-in-molecules fragment-based method for the calculation of chiroptical spectra of large molecules: Vibrational circular dichroism and Raman optical activity spectra of alanine polypeptides. Chirality 2016; 28:755-768. [DOI: 10.1002/chir.22651] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2016] [Revised: 08/29/2016] [Accepted: 08/30/2016] [Indexed: 12/17/2022]
Affiliation(s)
- K. V. Jovan Jose
- Department of Chemistry; Indiana University; Bloomington Indiana USA
| | | |
Collapse
|
29
|
Meitei OR, Heßelmann A. On the Stability of Cyclophane Derivates Using a Molecular Fragmentation Method. Chemphyschem 2016; 17:3863-3874. [PMID: 27653807 DOI: 10.1002/cphc.201600942] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2016] [Indexed: 11/11/2022]
Affiliation(s)
- Oinam Romesh Meitei
- Lehrstuhl für Theoretische Chemie; Universität Erlangen-Nürnberg; Egerlandstr. 3 D-91058 Erlangen Germany
| | - Andreas Heßelmann
- Lehrstuhl für Theoretische Chemie; Universität Erlangen-Nürnberg; Egerlandstr. 3 D-91058 Erlangen Germany
| |
Collapse
|
30
|
Fedorov DG, Kitaura K. Subsystem Analysis for the Fragment Molecular Orbital Method and Its Application to Protein-Ligand Binding in Solution. J Phys Chem A 2016; 120:2218-31. [PMID: 26949816 DOI: 10.1021/acs.jpca.6b00163] [Citation(s) in RCA: 47] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A subsystem analysis is derived incorporating interfragment interactions into the fragment properties, such as energies or charges. The relative stabilities of three alanine isomers, the α-helix, the β-turn, and the extended form are studied and the differences in fragment properties are elucidated. The analysis is further elaborated for studies of binding energies. The binding of the Trp-cage protein (PDB: 1L2Y ) to two ligands is studied in detail. Binding energies defined for each fragment can be used as a convenient descriptor for analyzing contributions to binding in solution.
Collapse
Affiliation(s)
- Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST) , Central 2, Umezono 1-1-1, Tsukuba, 305-8568, Japan
| | - Kazuo Kitaura
- Graduate School of System Informatics, Kobe University , 1-1 Rokkodai-cho, Nada-ku, Kobe 657-8501, Japan
| |
Collapse
|
31
|
Jovan Jose KV, Raghavachari K. Raman Optical Activity Spectra for Large Molecules through Molecules-in-Molecules Fragment-Based Approach. J Chem Theory Comput 2016; 12:585-94. [DOI: 10.1021/acs.jctc.5b01127] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- K. V. Jovan Jose
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Krishnan Raghavachari
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| |
Collapse
|
32
|
Jovan Jose K, Raghavachari K. Molecules-in-molecules fragment-based method for the evaluation of Raman spectra of large molecules. Mol Phys 2015. [DOI: 10.1080/00268976.2015.1074744] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
33
|
Jose KVJ, Beckett D, Raghavachari K. Vibrational Circular Dichroism Spectra for Large Molecules through Molecules-in-Molecules Fragment-Based Approach. J Chem Theory Comput 2015; 11:4238-47. [DOI: 10.1021/acs.jctc.5b00647] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- K. V. Jovan Jose
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Daniel Beckett
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| | - Krishnan Raghavachari
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| |
Collapse
|
34
|
Sahu N, Khire SS, Gadre SR. Structures, energetics and vibrational spectra of (H2O)32clusters: a journey from model potentials to correlated theory. Mol Phys 2015. [DOI: 10.1080/00268976.2015.1062150] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|