51
|
Chirality-induced spin polarization places symmetry constraints on biomolecular interactions. Proc Natl Acad Sci U S A 2017; 114:2474-2478. [PMID: 28228525 DOI: 10.1073/pnas.1611467114] [Citation(s) in RCA: 121] [Impact Index Per Article: 17.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
Noncovalent interactions between molecules are key for many biological processes. Necessarily, when molecules interact, the electronic charge in each of them is redistributed. Here, we show experimentally that, in chiral molecules, charge redistribution is accompanied by spin polarization. We describe how this spin polarization adds an enantioselective term to the forces, so that homochiral interaction energies differ from heterochiral ones. The spin polarization was measured by using a modified Hall effect device. An electric field that is applied along the molecules causes charge redistribution, and for chiral molecules, a Hall voltage is measured that indicates the spin polarization. Based on this observation, we conjecture that the spin polarization enforces symmetry constraints on the biorecognition process between two chiral molecules, and we describe how these constraints can lead to selectivity in the interaction between enantiomers based on their handedness. Model quantum chemistry calculations that rigorously enforce these constraints show that the interaction energy for methyl groups on homochiral molecules differs significantly from that found for heterochiral molecules at van der Waals contact and shorter (i.e., ∼0.5 kcal/mol at 0.26 nm).
Collapse
|
52
|
Hakulinen R, Puranen S. Probabilistic model to treat flexibility in molecular contacts. Mol Phys 2016. [DOI: 10.1080/00268976.2016.1225129] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
Affiliation(s)
- Riku Hakulinen
- Structural Bioinformatics Laboratory/Biochemistry, Faculty of Science and Engineering, Åbo Akademi University, Turku, Finland
| | - Santeri Puranen
- Department of Computer Science, Aalto University, Espoo, Finland
- Department of Computer Science, University of Helsinki, Helsinki, Finland
- Helsinki Institute for Information Technology HIIT, Helsinki, Finland
| |
Collapse
|
53
|
Mereghetti P, Maccari G, Spampinato GLB, Tozzini V. Optimization of Analytical Potentials for Coarse-Grained Biopolymer Models. J Phys Chem B 2016; 120:8571-9. [PMID: 27150459 DOI: 10.1021/acs.jpcb.6b02555] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
The increasing trend in the recent literature on coarse grained (CG) models testifies their impact in the study of complex systems. However, the CG model landscape is variegated: even considering a given resolution level, the force fields are very heterogeneous and optimized with very different parametrization procedures. Along the road for standardization of CG models for biopolymers, here we describe a strategy to aid building and optimization of statistics based analytical force fields and its implementation in the software package AsParaGS (Assisted Parameterization platform for coarse Grained modelS). Our method is based on the use and optimization of analytical potentials, optimized by targeting internal variables statistical distributions by means of the combination of different algorithms (i.e., relative entropy driven stochastic exploration of the parameter space and iterative Boltzmann inversion). This allows designing a custom model that endows the force field terms with a physically sound meaning. Furthermore, the level of transferability and accuracy can be tuned through the choice of statistical data set composition. The method-illustrated by means of applications to helical polypeptides-also involves the analysis of two and three variable distributions, and allows handling issues related to the FF term correlations. AsParaGS is interfaced with general-purpose molecular dynamics codes and currently implements the "minimalist" subclass of CG models (i.e., one bead per amino acid, Cα based). Extensions to nucleic acids and different levels of coarse graining are in the course.
Collapse
Affiliation(s)
- Paolo Mereghetti
- Center for Nanotechnology and Innovation @NEST, Istituto Italiano di Tecnologia , Piazza San Silvestro 12, 56127 Pisa, Italy
| | - Giuseppe Maccari
- Center for Nanotechnology and Innovation @NEST, Istituto Italiano di Tecnologia , Piazza San Silvestro 12, 56127 Pisa, Italy
| | | | - Valentina Tozzini
- NEST, Istituto Nanoscienze - CNR and Scuola Normale Superiore , Piazza San Silvestro 12, 56127 Pisa, Italy
| |
Collapse
|
54
|
Abstract
Path-integral free energy perturbation (PI-FEP) theory is presented to directly determine the ratio of quantum mechanical partition functions of different isotopologs in a single simulation. Furthermore, a double averaging strategy is used to carry out the practical simulation, separating the quantum mechanical path integral exactly into two separate calculations, one corresponding to a classical molecular dynamics simulation of the centroid coordinates, and another involving free-particle path-integral sampling over the classical, centroid positions. An integrated centroid path-integral free energy perturbation and umbrella sampling (PI-FEP/UM, or simply, PI-FEP) method along with bisection sampling was summarized, which provides an accurate and fast convergent method for computing kinetic isotope effects for chemical reactions in solution and in enzymes. The PI-FEP method is illustrated by a number of applications, to highlight the computational precision and accuracy, the rule of geometrical mean in kinetic isotope effects, enhanced nuclear quantum effects in enzyme catalysis, and protein dynamics on temperature dependence of kinetic isotope effects.
Collapse
Affiliation(s)
- J Gao
- Theoretical Chemistry Institute, Jilin University, Changchun, Jilin Province, PR China; University of Minnesota, Minneapolis, MN, United States.
| |
Collapse
|
55
|
Cisneros G, Wikfeldt KT, Ojamäe L, Lu J, Xu Y, Torabifard H, Bartók AP, Csányi G, Molinero V, Paesani F. Modeling Molecular Interactions in Water: From Pairwise to Many-Body Potential Energy Functions. Chem Rev 2016; 116:7501-28. [PMID: 27186804 PMCID: PMC5450669 DOI: 10.1021/acs.chemrev.5b00644] [Citation(s) in RCA: 272] [Impact Index Per Article: 34.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2015] [Indexed: 12/17/2022]
Abstract
Almost 50 years have passed from the first computer simulations of water, and a large number of molecular models have been proposed since then to elucidate the unique behavior of water across different phases. In this article, we review the recent progress in the development of analytical potential energy functions that aim at correctly representing many-body effects. Starting from the many-body expansion of the interaction energy, specific focus is on different classes of potential energy functions built upon a hierarchy of approximations and on their ability to accurately reproduce reference data obtained from state-of-the-art electronic structure calculations and experimental measurements. We show that most recent potential energy functions, which include explicit short-range representations of two-body and three-body effects along with a physically correct description of many-body effects at all distances, predict the properties of water from the gas to the condensed phase with unprecedented accuracy, thus opening the door to the long-sought "universal model" capable of describing the behavior of water under different conditions and in different environments.
Collapse
Affiliation(s)
| | - Kjartan Thor Wikfeldt
- Science
Institute, University of Iceland, VR-III, 107, Reykjavik, Iceland
- Department
of Physics, Albanova, Stockholm University, S-106 91 Stockholm, Sweden
| | - Lars Ojamäe
- Department
of Chemistry, Linköping University, SE-581 83 Linköping, Sweden
| | - Jibao Lu
- Department
of Chemistry, The University of Utah, Salt Lake City, Utah 84112-0850, United States
| | - Yao Xu
- Lehrstuhl
Physikalische Chemie II, Ruhr-Universität
Bochum, 44801 Bochum, Germany
| | - Hedieh Torabifard
- Department
of Chemistry, Wayne State University, Detroit, Michigan 48202, United States
| | - Albert P. Bartók
- Engineering
Laboratory, University of Cambridge, Trumpington Street, Cambridge CB21PZ, United Kingdom
| | - Gábor Csányi
- Engineering
Laboratory, University of Cambridge, Trumpington Street, Cambridge CB21PZ, United Kingdom
| | - Valeria Molinero
- Department
of Chemistry, The University of Utah, Salt Lake City, Utah 84112-0850, United States
| | - Francesco Paesani
- Department
of Chemistry and Biochemistry, University
of California San Diego, La Jolla, California 92093, United States
| |
Collapse
|
56
|
Loco D, Polack É, Caprasecca S, Lagardère L, Lipparini F, Piquemal JP, Mennucci B. A QM/MM Approach Using the AMOEBA Polarizable Embedding: From Ground State Energies to Electronic Excitations. J Chem Theory Comput 2016; 12:3654-61. [DOI: 10.1021/acs.jctc.6b00385] [Citation(s) in RCA: 109] [Impact Index Per Article: 13.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Daniele Loco
- Dipartimento
di Chimica e Chimica Industriale, Università di Pisa, via G. Moruzzi
13, I-56124 Pisa, Italy
| | - Étienne Polack
- Sorbonne Universités,
UPMC Univ. Paris 06, UMR 7616, Laboratoire de Chimie Théorique, F-75005, Paris, France
- Sorbonne Universités,
UPMC Univ. Paris 06, UMR 7598, Laboratoire Jacques-Louis Lions, F-75005, Paris, France
| | - Stefano Caprasecca
- Dipartimento
di Chimica e Chimica Industriale, Università di Pisa, via G. Moruzzi
13, I-56124 Pisa, Italy
| | - Louis Lagardère
- Sorbonne Universités,
UPMC Univ. Paris 06, UMR 7616, Laboratoire de Chimie Théorique, F-75005, Paris, France
- Sorbonne Universités,
UPMC Univ. Paris 06, Institut du Calcul et de la Simulation, F-75005, Paris, France
| | - Filippo Lipparini
- Institut
für Physikalische Chemie, Universität Mainz, Duesbergweg 10-14, D-55128 Mainz, Germany
| | - Jean-Philip Piquemal
- Sorbonne Universités,
UPMC Univ. Paris 06, UMR 7616, Laboratoire de Chimie Théorique, F-75005, Paris, France
- CNRS, UMR 7598 and 7616, F-75005, Paris, France
- Department
of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA
| | - Benedetta Mennucci
- Dipartimento
di Chimica e Chimica Industriale, Università di Pisa, via G. Moruzzi
13, I-56124 Pisa, Italy
| |
Collapse
|
57
|
Giese TJ, Panteva MT, Chen H, York DM. Multipolar Ewald methods, 2: applications using a quantum mechanical force field. J Chem Theory Comput 2016; 11:451-61. [PMID: 25691830 PMCID: PMC4325604 DOI: 10.1021/ct500799g] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2014] [Indexed: 01/08/2023]
Abstract
A fully quantum mechanical force field (QMFF) based on a modified “divide-and-conquer” (mDC) framework is applied to a series of molecular simulation applications, using a generalized Particle Mesh Ewald method extended to multipolar charge densities. Simulation results are presented for three example applications: liquid water, p-nitrophenylphosphate reactivity in solution, and crystalline N,N-dimethylglycine. Simulations of liquid water using a parametrized mDC model are compared to TIP3P and TIP4P/Ew water models and experiment. The mDC model is shown to be superior for cluster binding energies and generally comparable for bulk properties. Examination of the dissociative pathway for dephosphorylation of p-nitrophenylphosphate shows that the mDC method evaluated with the DFTB3/3OB and DFTB3/OPhyd semiempirical models bracket the experimental barrier, whereas DFTB2 and AM1/d-PhoT QM/MM simulations exhibit deficiencies in the barriers, the latter for which is related, in part, to the anomalous underestimation of the p-nitrophenylate leaving group pKa. Simulations of crystalline N,N-dimethylglycine are performed and the overall structure and atomic fluctuations are compared with the experiment and the general AMBER force field (GAFF). The QMFF, which was not parametrized for this application, was shown to be in better agreement with crystallographic data than GAFF. Our simulations highlight some of the application areas that may benefit from using new QMFFs, and they demonstrate progress toward the development of accurate QMFFs using the recently developed mDC framework.
Collapse
Affiliation(s)
- Timothy J Giese
- Center for Integrative Proteomics Research, BioMaPS Institute for Quantitative Biology and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854-8087, United States
| | | | | | | |
Collapse
|
58
|
Fast calculation of molecular total energy with ABEEMσπ/MM method – For some series of organic molecules and peptides. Chem Phys 2016. [DOI: 10.1016/j.chemphys.2016.03.004] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
59
|
Christensen A, Kubař T, Cui Q, Elstner M. Semiempirical Quantum Mechanical Methods for Noncovalent Interactions for Chemical and Biochemical Applications. Chem Rev 2016; 116:5301-37. [PMID: 27074247 PMCID: PMC4867870 DOI: 10.1021/acs.chemrev.5b00584] [Citation(s) in RCA: 246] [Impact Index Per Article: 30.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2015] [Indexed: 12/28/2022]
Abstract
Semiempirical (SE) methods can be derived from either Hartree-Fock or density functional theory by applying systematic approximations, leading to efficient computational schemes that are several orders of magnitude faster than ab initio calculations. Such numerical efficiency, in combination with modern computational facilities and linear scaling algorithms, allows application of SE methods to very large molecular systems with extensive conformational sampling. To reliably model the structure, dynamics, and reactivity of biological and other soft matter systems, however, good accuracy for the description of noncovalent interactions is required. In this review, we analyze popular SE approaches in terms of their ability to model noncovalent interactions, especially in the context of describing biomolecules, water solution, and organic materials. We discuss the most significant errors and proposed correction schemes, and we review their performance using standard test sets of molecular systems for quantum chemical methods and several recent applications. The general goal is to highlight both the value and limitations of SE methods and stimulate further developments that allow them to effectively complement ab initio methods in the analysis of complex molecular systems.
Collapse
Affiliation(s)
- Anders
S. Christensen
- Department
of Chemistry and Theoretical Chemistry Institute, University of Wisconsin—Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Tomáš Kubař
- Institute of Physical
Chemistry & Center for Functional Nanostructures and Institute of Physical
Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
| | - Qiang Cui
- Department
of Chemistry and Theoretical Chemistry Institute, University of Wisconsin—Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Marcus Elstner
- Institute of Physical
Chemistry & Center for Functional Nanostructures and Institute of Physical
Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
| |
Collapse
|
60
|
Cao Z, Voth GA. The multiscale coarse-graining method. XI. Accurate interactions based on the centers of charge of coarse-grained sites. J Chem Phys 2016; 143:243116. [PMID: 26723601 DOI: 10.1063/1.4933249] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023] Open
Abstract
It is essential to be able to systematically construct coarse-grained (CG) models that can efficiently and accurately reproduce key properties of higher-resolution models such as all-atom. To fulfill this goal, a mapping operator is needed to transform the higher-resolution configuration to a CG configuration. Certain mapping operators, however, may lose information related to the underlying electrostatic properties. In this paper, a new mapping operator based on the centers of charge of CG sites is proposed to address this issue. Four example systems are chosen to demonstrate this concept. Within the multiscale coarse-graining framework, CG models that use this mapping operator are found to better reproduce the structural correlations of atomistic models. The present work also demonstrates the flexibility of the mapping operator and the robustness of the force matching method. For instance, important functional groups can be isolated and emphasized in the CG model.
Collapse
Affiliation(s)
- Zhen Cao
- Department of Chemistry, James Franck Institute, Institute for Biophysical Dynamics, and Computation Institute, The University of Chicago, 5735 S Ellis Ave., Chicago, Illinois 60637, USA
| | - Gregory A Voth
- Department of Chemistry, James Franck Institute, Institute for Biophysical Dynamics, and Computation Institute, The University of Chicago, 5735 S Ellis Ave., Chicago, Illinois 60637, USA
| |
Collapse
|
61
|
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
|
62
|
Abstract
The design of optimal light-harvesting (supra)molecular systems and materials is one of the most challenging frontiers of science. Theoretical methods and computational models play a fundamental role in this difficult task, as they allow the establishment of structural blueprints inspired by natural photosynthetic organisms that can be applied to the design of novel artificial light-harvesting devices. Among theoretical strategies, the application of quantum chemical tools represents an important reality that has already reached an evident degree of maturity, although it still has to show its real potentials. This Review presents an overview of the state of the art of this strategy, showing the actual fields of applicability but also indicating its current limitations, which need to be solved in future developments.
Collapse
Affiliation(s)
- Carles Curutchet
- Departament de Fisicoquímica, Facultat de Farmàcia, Universitat de Barcelona , Av. Joan XXIII s/n, 08028 Barcelona, Spain
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, University of Pisa , via G. Moruzzi 13, 56124 Pisa, Italy
| |
Collapse
|
63
|
|
64
|
Nishimoto Y, Fedorov DG. The fragment molecular orbital method combined with density-functional tight-binding and the polarizable continuum model. Phys Chem Chem Phys 2016; 18:22047-61. [DOI: 10.1039/c6cp02186g] [Citation(s) in RCA: 49] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The electronic gap in proteins is analyzed in detail, and it is shown that FMO-DFTB/PCM is efficient and accurate in describing the molecular structure of proteins in solution.
Collapse
Affiliation(s)
- Yoshio Nishimoto
- Fukui Institute for Fundamental Chemistry
- Kyoto University
- 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)
- Tsukuba
- Japan
| |
Collapse
|
65
|
Hagler AT. Quantum Derivative Fitting and Biomolecular Force Fields: Functional Form, Coupling Terms, Charge Flux, Nonbond Anharmonicity, and Individual Dihedral Potentials. J Chem Theory Comput 2015; 11:5555-72. [DOI: 10.1021/acs.jctc.5b00666] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Affiliation(s)
- A. T. Hagler
- Department of Chemistry University of Massachusetts, Amherst, Massachusetts 01003, United States
- Shifa Biopharm, Shifa Biomedical Corporation, Malvern, Pennsylvania 19355, United States
| |
Collapse
|
66
|
Liu J, Zhu T, Wang X, He X, Zhang JZH. Quantum Fragment Based ab Initio Molecular Dynamics for Proteins. J Chem Theory Comput 2015; 11:5897-905. [PMID: 26642993 DOI: 10.1021/acs.jctc.5b00558] [Citation(s) in RCA: 51] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Developing ab initio molecular dynamics (AIMD) methods for practical application in protein dynamics is of significant interest. Due to the large size of biomolecules, applying standard quantum chemical methods to compute energies for dynamic simulation is computationally prohibitive. In this work, a fragment based ab initio molecular dynamics approach is presented for practical application in protein dynamics study. In this approach, the energy and forces of the protein are calculated by a recently developed electrostatically embedded generalized molecular fractionation with conjugate caps (EE-GMFCC) method. For simulation in explicit solvent, mechanical embedding is introduced to treat protein interaction with explicit water molecules. This AIMD approach has been applied to MD simulations of a small benchmark protein Trpcage (with 20 residues and 304 atoms) in both the gas phase and in solution. Comparison to the simulation result using the AMBER force field shows that the AIMD gives a more stable protein structure in the simulation, indicating that quantum chemical energy is more reliable. Importantly, the present fragment-based AIMD simulation captures quantum effects including electrostatic polarization and charge transfer that are missing in standard classical MD simulations. The current approach is linear-scaling, trivially parallel, and applicable to performing the AIMD simulation of proteins with a large size.
Collapse
Affiliation(s)
- Jinfeng Liu
- State Key Laboratory of Precision Spectroscopy, Institute of Theoretical and Computational Science, East China Normal University , Shanghai 200062, China
| | - Tong Zhu
- 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
| | - Xianwei Wang
- Center for Optics & Optoelectronics Research, College of Science, Zhejiang University of Technology , Hangzhou, Zhejiang 310023, China
| | - 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
| | - 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.,Department of Chemistry, New York University , New York, New York 10003, United States
| |
Collapse
|
67
|
Ghosh S, Parui S, Jana B, Bhattacharyya K. Ionic liquid induced dehydration and domain closure in lysozyme: FCS and MD simulation. J Chem Phys 2015; 143:125103. [DOI: 10.1063/1.4931974] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Shirsendu Ghosh
- Department of Physical Chemistry, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700 032, India
| | - Sridip Parui
- Department of Physical Chemistry, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700 032, India
| | - Biman Jana
- Department of Physical Chemistry, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700 032, India
| | - Kankan Bhattacharyya
- Department of Physical Chemistry, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700 032, India
| |
Collapse
|
68
|
Nakata H, Fedorov DG, Nagata T, Kitaura K, Nakamura S. Simulations of Chemical Reactions with the Frozen Domain Formulation of the Fragment Molecular Orbital Method. J Chem Theory Comput 2015; 11:3053-64. [DOI: 10.1021/acs.jctc.5b00277] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Hiroya Nakata
- Department
of Biomolecular Engineering, Tokyo Institute of Technology, 4259 Nagatsuta-cho, Midori-ku, Yokohama, Kanagawa 226-8501, Japan
- Research
Cluster for Innovation, RIKEN, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
- Japan Society for the Promotion of Science, Kojimachi Business Center Building, 5-3-1 Kojimachi, Chiyoda-ku, Tokyo 102-0083, Japan
| | - Dmitri G. Fedorov
- Nanosystem
Research Institute, National Institute of Advanced Industrial Science and Technology, 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan
| | - Takeshi Nagata
- Nanosystem
Research Institute, National Institute of Advanced Industrial Science and Technology, 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan
- Graduate
School of System Informatics, Kobe University, 1-1 Rokkodai-cho, Nada-ku, Kobe 657-8501, Japan
| | - Kazuo Kitaura
- Graduate
School of System Informatics, Kobe University, 1-1 Rokkodai-cho, Nada-ku, Kobe 657-8501, Japan
| | - Shinichiro Nakamura
- Research
Cluster for Innovation, RIKEN, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| |
Collapse
|
69
|
Giese T, Panteva MT, Chen H, York DM. Multipolar Ewald methods, 1: theory, accuracy, and performance. J Chem Theory Comput 2015; 11:436-50. [PMID: 25691829 PMCID: PMC4325605 DOI: 10.1021/ct5007983] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2014] [Indexed: 11/29/2022]
Abstract
The Ewald, Particle Mesh Ewald (PME), and Fast Fourier–Poisson (FFP) methods are developed for systems composed of spherical multipole moment expansions. A unified set of equations is derived that takes advantage of a spherical tensor gradient operator formalism in both real space and reciprocal space to allow extension to arbitrary multipole order. The implementation of these methods into a novel linear-scaling modified “divide-and-conquer” (mDC) quantum mechanical force field is discussed. The evaluation times and relative force errors are compared between the three methods, as a function of multipole expansion order. Timings and errors are also compared within the context of the quantum mechanical force field, which encounters primary errors related to the quality of reproducing electrostatic forces for a given density matrix and secondary errors resulting from the propagation of the approximate electrostatics into the self-consistent field procedure, which yields a converged, variational, but nonetheless approximate density matrix. Condensed-phase simulations of an mDC water model are performed with the multipolar PME method and compared to an electrostatic cutoff method, which is shown to artificially increase the density of water and heat of vaporization relative to full electrostatic treatment.
Collapse
Affiliation(s)
- Timothy
J. Giese
- Center for Integrative Proteomics
Research, BioMaPS Institute for Quantitative Biology and Department
of Chemistry and Chemical Biology, Rutgers
University, Piscataway, New Jersey 08854-8087, United States
| | - Maria T. Panteva
- Center for Integrative Proteomics
Research, BioMaPS Institute for Quantitative Biology and Department
of Chemistry and Chemical Biology, Rutgers
University, Piscataway, New Jersey 08854-8087, United States
| | - Haoyuan Chen
- Center for Integrative Proteomics
Research, BioMaPS Institute for Quantitative Biology and Department
of Chemistry and Chemical Biology, Rutgers
University, Piscataway, New Jersey 08854-8087, United States
| | - Darrin M. York
- Center for Integrative Proteomics
Research, BioMaPS Institute for Quantitative Biology and Department
of Chemistry and Chemical Biology, Rutgers
University, Piscataway, New Jersey 08854-8087, United States
| |
Collapse
|
70
|
Omer A, Suryanarayanan V, Selvaraj C, Singh SK, Singh P. Explicit Drug Re-positioning: Predicting Novel Drug-Target Interactions of the Shelved Molecules with QM/MM Based Approaches. ADVANCES IN PROTEIN CHEMISTRY AND STRUCTURAL BIOLOGY 2015; 100:89-112. [PMID: 26415842 DOI: 10.1016/bs.apcsb.2015.07.001] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
With the demand to enhance the speed of the drug discovery process there has been an increased usage of computational approaches in drug discovery studies. However because of their probabilistic outcomes, the challenge is to exactly mimic the natural environment which can provide the exact charge polarization effect while estimating the binding energy between protein and ligand. There has been a large number of scoring functions from simple one to the complex one available for estimating binding energy. The quantum mechanics/molecular mechanics (QM/MM) hybrid approach has been the preferred choice of interest since last decade for modeling reactions in biomolecular systems. The application of QM/MM approach has been expanded right from rescoring the already known complexes and depicting the correct position of some novel molecule to ranking a large number of molecules. It is expected that the application of QM/MM-based scoring will grow in all areas of drug discovery. However, the most promising area will be its application in repositioning, that is, assigning novel functions or targets to the already existing drugs, as this would stop the rising attrition rates as well as reduce the overall time and cost of drug discovery procedure.
Collapse
Affiliation(s)
- Ankur Omer
- Division of Toxicology, CSIR-Central Drug Research Institute, Lucknow, Uttar Pradesh, India; Academy of Scientific & Innovative Research (AcSIR), New Delhi, India
| | - Venkatesan Suryanarayanan
- Department of Bioinformatics, Computer Aided Drug Design and Molecular Modeling Lab, Science Block, Alagappa University, Karaikudi, Tamil Nadu, India
| | - Chandrabose Selvaraj
- Department of Bioinformatics, Computer Aided Drug Design and Molecular Modeling Lab, Science Block, Alagappa University, Karaikudi, Tamil Nadu, India
| | - Sanjeev Kumar Singh
- Department of Bioinformatics, Computer Aided Drug Design and Molecular Modeling Lab, Science Block, Alagappa University, Karaikudi, Tamil Nadu, India.
| | - Poonam Singh
- Division of Toxicology, CSIR-Central Drug Research Institute, Lucknow, Uttar Pradesh, India; Academy of Scientific & Innovative Research (AcSIR), New Delhi, India.
| |
Collapse
|
71
|
Matthews RP, Welton T, Hunt PA. Hydrogen bonding and π–π interactions in imidazolium-chloride ionic liquid clusters. Phys Chem Chem Phys 2015; 17:14437-53. [DOI: 10.1039/c5cp00459d] [Citation(s) in RCA: 95] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The importance of 1° and 2° hydrogen-bonding and anion–π+ interactions for ionic liquid structuring.
Collapse
Affiliation(s)
| | - Tom Welton
- Department of Chemistry
- Imperial College London
- London
- UK
| | | |
Collapse
|
72
|
Ramabhadran RO, Raghavachari K. The successful merger of theoretical thermochemistry with fragment-based methods in quantum chemistry. Acc Chem Res 2014; 47:3596-604. [PMID: 25393551 DOI: 10.1021/ar500294s] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
CONSPECTUS: Quantum chemistry and electronic structure theory have proven to be essential tools to the experimental chemist, in terms of both a priori predictions that pave the way for designing new experiments and rationalizing experimental observations a posteriori. Translating the well-established success of electronic structure theory in obtaining the structures and energies of small chemical systems to increasingly larger molecules is an exciting and ongoing central theme of research in quantum chemistry. However, the prohibitive computational scaling of highly accurate ab initio electronic structure methods poses a fundamental challenge to this research endeavor. This scenario necessitates an indirect fragment-based approach wherein a large molecule is divided into small fragments and is subsequently reassembled to compute its energy accurately. In our quest to further reduce the computational expense associated with the fragment-based methods and overall enhance the applicability of electronic structure methods to large molecules, we realized that the broad ideas involved in a different area, theoretical thermochemistry, are transferable to the area of fragment-based methods. This Account focuses on the effective merger of these two disparate frontiers in quantum chemistry and how new concepts inspired by theoretical thermochemistry significantly reduce the total number of electronic structure calculations needed to be performed as part of a fragment-based method without any appreciable loss of accuracy. Throughout, the generalized connectivity based hierarchy (CBH), which we developed to solve a long-standing problem in theoretical thermochemistry, serves as the linchpin in this merger. The accuracy of our method is based on two strong foundations: (a) the apt utilization of systematic and sophisticated error-canceling schemes via CBH that result in an optimal cutting scheme at any given level of fragmentation and (b) the use of a less expensive second layer of electronic structure method to recover all the missing long-range interactions in the parent large molecule. Overall, the work featured here dramatically decreases the computational expense and empowers the execution of very accurate ab initio calculations (gold-standard CCSD(T)) on large molecules and thereby facilitates sophisticated electronic structure applications to a wide range of important chemical problems.
Collapse
Affiliation(s)
| | - Krishnan Raghavachari
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| |
Collapse
|
73
|
Lao KU, Herbert JM. Accurate and Efficient Quantum Chemistry Calculations for Noncovalent Interactions in Many-Body Systems: The XSAPT Family of Methods. J Phys Chem A 2014; 119:235-52. [DOI: 10.1021/jp5098603] [Citation(s) in RCA: 72] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Ka Un Lao
- 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
|