1
|
Malaspina LA, Genoni A, Jayatilaka D, Turner MJ, Sugimoto K, Nishibori E, Grabowsky S. The advanced treatment of hydrogen bonding in quantum crystallography. J Appl Crystallogr 2021; 54:718-729. [PMID: 34188611 PMCID: PMC8202034 DOI: 10.1107/s1600576721001126] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2020] [Accepted: 01/31/2021] [Indexed: 11/16/2022] Open
Abstract
Although hydrogen bonding is one of the most important motifs in chemistry and biology, H-atom parameters are especially problematic to refine against X-ray diffraction data. New developments in quantum crystallography offer a remedy. This article reports how hydrogen bonds are treated in three different quantum-crystallographic methods: Hirshfeld atom refinement (HAR), HAR coupled to extremely localized molecular orbitals and X-ray wavefunction refinement. Three different compound classes that form strong intra- or intermolecular hydrogen bonds are used as test cases: hydrogen maleates, the tripeptide l-alanyl-glycyl-l-alanine co-crystallized with water, and xylitol. The differences in the quantum-mechanical electron densities underlying all the used methods are analysed, as well as how these differences impact on the refinement results.
Collapse
Affiliation(s)
- Lorraine A. Malaspina
- Universität Bern, Departement für Chemie, Biochemie und Pharmazie, Freiestrasse 3, 3012 Bern, Switzerland
- Universität Bremen, Fachbereich 2 – Biologie/Chemie, Institut für Anorganische Chemie und Kristallographie, Leobener Strasse 3, 28359 Bremen, Germany
| | - Alessandro Genoni
- Université de Lorraine and CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019, 1 Boulevard Arago, 57078 Metz, France
| | - Dylan Jayatilaka
- The University of Western Australia, School of Molecular Sciences, 35 Stirling Highway, Perth, WA 6009, Australia
| | - Michael J. Turner
- The University of Western Australia, School of Molecular Sciences, 35 Stirling Highway, Perth, WA 6009, Australia
| | - Kunihisa Sugimoto
- Japan Synchrotron Radiation Research Institute/Diffraction and Scattering Division, 1-1-1 Kouto, Sayo-cho, Sayo-gun, Hyogo 679-5198, Japan
- Institute for Integrated Cell-Material Sciences (iCeMS), Kyoto University, Yoshida-Ushinomiya-cho, Sakyo-ku, Kyoto 606-8501, Japan
| | - Eiji Nishibori
- Department of Physics, Faculty of Pure and Applied Sciences, Tsukuba Research Center for Energy Materials Science (TREMS), University of Tsukuba, Tsukuba, Japan
| | - Simon Grabowsky
- Universität Bern, Departement für Chemie, Biochemie und Pharmazie, Freiestrasse 3, 3012 Bern, Switzerland
- Universität Bremen, Fachbereich 2 – Biologie/Chemie, Institut für Anorganische Chemie und Kristallographie, Leobener Strasse 3, 28359 Bremen, Germany
| |
Collapse
|
2
|
Kleemiss F, Wieduwilt EK, Hupf E, Shi MW, Stewart SG, Jayatilaka D, Turner MJ, Sugimoto K, Nishibori E, Schirmeister T, Schmidt TC, Engels B, Grabowsky S. Similarities and Differences between Crystal and Enzyme Environmental Effects on the Electron Density of Drug Molecules. Chemistry 2021; 27:3407-3419. [PMID: 33090581 PMCID: PMC7898524 DOI: 10.1002/chem.202003978] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2020] [Indexed: 01/28/2023]
Abstract
The crystal interaction density is generally assumed to be a suitable measure of the polarization of a low-molecular weight ligand inside an enzyme, but this approximation has seldomly been tested and has never been quantified before. In this study, we compare the crystal interaction density and the interaction electrostatic potential for a model compound of loxistatin acid (E64c) with those inside cathepsin B, in solution, and in vacuum. We apply QM/MM calculations and experimental quantum crystallography to show that the crystal interaction density is indeed very similar to the enzyme interaction density. Less than 0.1 e are shifted between these two environments in total. However, this difference has non-negligible consequences for derived properties.
Collapse
Affiliation(s)
- Florian Kleemiss
- Department 2 – Biology/Chemistry, Institute of Inorganic Chemistry and CrystallographyUniversity of BremenLeobener Str. 3 and 7, 28359 BremenGermany
- Department of Chemistry and BiochemistryUniversity of BernFreiestrasse 3, 3012 BernSwitzerland
| | - Erna K. Wieduwilt
- Department 2 – Biology/Chemistry, Institute of Inorganic Chemistry and CrystallographyUniversity of BremenLeobener Str. 3 and 7, 28359 BremenGermany
- Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019Université de Lorraine & CNRSBoulevard Arago, 57078 MetzFrance
| | - Emanuel Hupf
- Department 2 – Biology/Chemistry, Institute of Inorganic Chemistry and CrystallographyUniversity of BremenLeobener Str. 3 and 7, 28359 BremenGermany
| | - Ming W. Shi
- School of Molecular SciencesUniversity of Western Australia35 Stirling Highway, Perth WA 6009Australia
| | - Scott G. Stewart
- School of Molecular SciencesUniversity of Western Australia35 Stirling Highway, Perth WA 6009Australia
| | - Dylan Jayatilaka
- School of Molecular SciencesUniversity of Western Australia35 Stirling Highway, Perth WA 6009Australia
| | - Michael J. Turner
- School of Molecular SciencesUniversity of Western Australia35 Stirling Highway, Perth WA 6009Australia
| | - Kunihisa Sugimoto
- Japan Synchrotron Radiation Research InstituteSPring-81-1-1 Kouto, Sayo-cho, Sayo-gun, Hyogo 679-5198Japan
- Institute for Integrated Cell-Materials SciencesKyoto UniversityYoshida-Ushinomiya-cho, Sakyo-ku, Kyoto 606-8501Japan
| | - Eiji Nishibori
- Division of Physics, Faculty of Pure and Applied Sciences, Tsukuba Research Center for Energy Materials ScienceUniversity of TsukubaTsukubaJapan
| | - Tanja Schirmeister
- Institute of Pharmaceutical and Biomedical SciencesJohannes-Gutenberg University MainzStaudingerweg 5, 55128 MainzGermany
| | - Thomas C. Schmidt
- Institute for Physical and Theoretical ChemistryJulius-Maximilians-University WürzburgEmil-Fischer-Str. 42, 97074 WürzburgGermany
| | - Bernd Engels
- Institute for Physical and Theoretical ChemistryJulius-Maximilians-University WürzburgEmil-Fischer-Str. 42, 97074 WürzburgGermany
| | - Simon Grabowsky
- Department 2 – Biology/Chemistry, Institute of Inorganic Chemistry and CrystallographyUniversity of BremenLeobener Str. 3 and 7, 28359 BremenGermany
- Department of Chemistry and BiochemistryUniversity of BernFreiestrasse 3, 3012 BernSwitzerland
| |
Collapse
|
3
|
Woinska M, Wanat M, Taciak P, Pawinski T, Minor W, Wozniak K. Energetics of interactions in the solid state of 2-hydroxy-8- X-quinoline derivatives ( X = Cl, Br, I, S-Ph): comparison of Hirshfeld atom, X-ray wavefunction and multipole refinements. IUCRJ 2019; 6:868-883. [PMID: 31576220 PMCID: PMC6760436 DOI: 10.1107/s2052252519007358] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/07/2019] [Accepted: 05/20/2019] [Indexed: 05/20/2023]
Abstract
In this work, two methods of high-resolution X-ray data refinement: multipole refinement (MM) and Hirshfeld atom refinement (HAR) - together with X-ray wavefunction refinement (XWR) - are applied to investigate the refinement of positions and anisotropic thermal motion of hydrogen atoms, experiment-based reconstruction of electron density, refinement of anharmonic thermal vibrations, as well as the effects of excluding the weakest reflections in the refinement. The study is based on X-ray data sets of varying quality collected for the crystals of four quinoline derivatives with Cl, Br, I atoms and the -S-Ph group as substituents. Energetic investigations are performed, comprising the calculation of the energy of intermolecular interactions, cohesive and geometrical relaxation energy. The results obtained for experimentally derived structures are verified against the values calculated for structures optimized using dispersion-corrected periodic density functional theory. For the high-quality data sets (the Cl and -S-Ph compounds), both MM and XWR could be successfully used to refine the atomic displacement parameters and the positions of hydrogen atoms; however, the bond lengths obtained with XWR were more precise and closer to the theoretical values. In the application to the more challenging data sets (the Br and I compounds), only XWR enabled free refinement of hydrogen atom geometrical parameters, nevertheless, the results clearly showed poor data quality. For both refinement methods, the energy values (intermolecular interactions, cohesive and relaxation) calculated for the experimental structures were in similar agreement with the values associated with the optimized structures - the most significant divergences were observed when experimental geometries were biased by poor data quality. XWR was found to be more robust in avoiding incorrect distortions of the reconstructed electron density as a result of data quality issues. Based on the problem of anharmonic thermal motion refinement, this study reveals that for the most correct interpretation of the obtained results, it is necessary to use the complete data set, including the weak reflections in order to draw conclusions.
Collapse
Affiliation(s)
- Magdalena Woinska
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Zwirki i Wigury 101, Warszawa 02-089, Poland
- Department of Molecular Physiology and Biological Physics, University of Virginia, Charlottesville, VA, USA
| | - Monika Wanat
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Zwirki i Wigury 101, Warszawa 02-089, Poland
- College of Inter-Faculty Individual Studies in Mathematics and Natural Sciences (MISMaP), University of Warsaw, Stefana Banacha 2C, Warszawa 02-097, Poland
| | - Przemyslaw Taciak
- Department of Drug Chemistry, Faculty of Pharmacy, Medicinal University of Warsaw, Stefana Banacha 1, Warszawa 02-091, Poland
- Department of Pharmacodynamics, Centre of Preclinical Research and Technology, Faculty of Pharmacy, Medical University of Warsaw, Stefana Banacha 1, Warszawa 02-091, Poland
| | - Tomasz Pawinski
- Department of Drug Chemistry, Faculty of Pharmacy, Medicinal University of Warsaw, Stefana Banacha 1, Warszawa 02-091, Poland
| | - Wladek Minor
- Department of Molecular Physiology and Biological Physics, University of Virginia, Charlottesville, VA, USA
- Department of Chemistry, University of Warsaw, Pasteura 1, Warszawa 02-093, Poland
| | - Krzysztof Wozniak
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Zwirki i Wigury 101, Warszawa 02-089, Poland
| |
Collapse
|
4
|
Grabowsky S, Genoni A, Bürgi HB. Quantum crystallography. Chem Sci 2017; 8:4159-4176. [PMID: 28878872 PMCID: PMC5576428 DOI: 10.1039/c6sc05504d] [Citation(s) in RCA: 97] [Impact Index Per Article: 13.9] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2016] [Accepted: 03/03/2017] [Indexed: 12/12/2022] Open
Abstract
Approximate wavefunctions can be improved by constraining them to reproduce observations derived from diffraction and scattering experiments. Conversely, charge density models, incorporating electron-density distributions, atomic positions and atomic motion, can be improved by supplementing diffraction experiments with quantum chemically calculated, tailor-made electron densities (form factors). In both cases quantum chemistry and diffraction/scattering experiments are combined into a single, integrated tool. The development of quantum crystallographic research is reviewed. Some results obtained by quantum crystallography illustrate the potential and limitations of this field.
Collapse
Affiliation(s)
- Simon Grabowsky
- Universität Bremen , Fachbereich 2 - Biologie/Chemie , Institut für Anorganische Chemie und Kristallographie , Leobener Str. NW2 , 28359 Bremen , Germany .
| | - Alessandro Genoni
- CNRS , Laboratoire SRSMC , UMR 7565 , Vandoeuvre-lès-Nancy , F-54506 , France
- Université de Lorraine , Laboratoire SRSMC , UMR 7565 , Vandoeuvre-lès-Nancy , F-54506 , France .
| | - Hans-Beat Bürgi
- Universität Bern , Departement für Chemie und Biochemie , Freiestr. 3 , CH-3012 Bern , Switzerland .
- Universität Zürich , Institut für Chemie , Winterthurerstrasse 190 , CH-8057 Zürich , Switzerland
| |
Collapse
|
5
|
Shi MW, Stewart SG, Sobolev AN, Dittrich B, Schirmeister T, Luger P, Hesse M, Chen Y, Spackman PR, Spackman MA, Grabowsky S. Approaching an experimental electron density model of the biologically active
trans
‐epoxysuccinyl amide group—Substituent effects vs. crystal packing. J PHYS ORG CHEM 2017. [DOI: 10.1002/poc.3683] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Ming W. Shi
- School of Chemistry and Biochemistry The University of Western Australia Perth WA Australia
| | - Scott G. Stewart
- School of Chemistry and Biochemistry The University of Western Australia Perth WA Australia
| | - Alexandre N. Sobolev
- School of Chemistry and Biochemistry The University of Western Australia Perth WA Australia
| | - Birger Dittrich
- Anorganische Chemie und Strukturchemie Heinrich‐Heine‐Universität Düsseldorf Düsseldorf Germany
| | - Tanja Schirmeister
- Institut für Pharmazie und Biochemie Johannes‐Gutenberg‐Universität Mainz Mainz Germany
| | - Peter Luger
- Institut für Chemie und Biochemie, Anorganische Chemie Freie Universität Berlin Berlin Germany
| | - Malte Hesse
- Fachbereich 2—Biologie/Chemie, Institut für Anorganische Chemie und Kristallographie Universität Bremen Bremen Germany
| | - Yu‐Sheng Chen
- ChemMatCARS The University of Chicago Argonne IL USA
| | - Peter R. Spackman
- School of Chemistry and Biochemistry The University of Western Australia Perth WA Australia
| | - Mark A. Spackman
- School of Chemistry and Biochemistry The University of Western Australia Perth WA Australia
| | - Simon Grabowsky
- School of Chemistry and Biochemistry The University of Western Australia Perth WA Australia
- Fachbereich 2—Biologie/Chemie, Institut für Anorganische Chemie und Kristallographie Universität Bremen Bremen Germany
| |
Collapse
|
6
|
|
7
|
Fucke K, McIntyre GJ, Lemée-Cailleau MH, Wilkinson C, Edwards AJ, Howard JAK, Steed JW. Insights into the crystallisation process from anhydrous, hydrated and solvated crystal forms of diatrizoic acid. Chemistry 2015; 21:1036-47. [PMID: 25370384 DOI: 10.1002/chem.201404693] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2014] [Indexed: 11/11/2022]
Abstract
Diatrizoic acid (DTA), a clinically used X-ray contrast agent, crystallises in two hydrated, three anhydrous and nine solvated solid forms, all of which have been characterised by X-ray crystallography. Single-crystal neutron structures of DTA dihydrate and monosodium DTA tetrahydrate have been determined. All of the solid-state structures have been analysed using partial atomic charges and hardness algorithm (PACHA) calculations. Even though in general all DTA crystal forms reveal similar intermolecular interactions, the overall crystal packing differs considerably from form to form. The water of the dihydrate is encapsulated between a pair of host molecules, which calculations reveal to be an extraordinarily stable motif. DTA presents functionalities that enable hydrogen and halogen bonding, and whilst an extended hydrogen-bonding network is realised in all crystal forms, halogen bonding is not present in the hydrated crystal forms. This is due to the formation of a hydrogen-bonding network based on individual enclosed water squares, which is not amenable to the concomitant formation of halogen bonds. The main interaction in the solvates involves the carboxylic acid, which corroborates the hypothesis that this strong interaction is the last one to be broken during the crystal desolvation and nucleation process.
Collapse
Affiliation(s)
- Katharina Fucke
- School of Medicine, Pharmacy and Health, Durham University, University Boulevard, Stockton-on-Tees, TS17 6BH (United Kingdom).
| | | | | | | | | | | | | |
Collapse
|
8
|
Shi MW, Sobolev AN, Schirmeister T, Engels B, Schmidt TC, Luger P, Mebs S, Dittrich B, Chen YS, Bąk JM, Jayatilaka D, Bond CS, Turner MJ, Stewart SG, Spackman MA, Grabowsky S. Electrostatic complementarity in pseudoreceptor modeling based on drug molecule crystal structures: the case of loxistatin acid (E64c). NEW J CHEM 2015. [DOI: 10.1039/c4nj01503g] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
A combination of pseudoreceptor modeling and electrostatic complementarity maps properties of a native pocket for an enzyme ligand.
Collapse
|
9
|
Capelli SC, Bürgi HB, Dittrich B, Grabowsky S, Jayatilaka D. Hirshfeld atom refinement. IUCRJ 2014; 1:361-79. [PMID: 25295177 PMCID: PMC4174878 DOI: 10.1107/s2052252514014845] [Citation(s) in RCA: 171] [Impact Index Per Article: 17.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/02/2014] [Accepted: 06/24/2014] [Indexed: 05/20/2023]
Abstract
Hirshfeld atom refinement (HAR) is a method which determines structural parameters from single-crystal X-ray diffraction data by using an aspherical atom partitioning of tailor-made ab initio quantum mechanical molecular electron densities without any further approximation. Here the original HAR method is extended by implementing an iterative procedure of successive cycles of electron density calculations, Hirshfeld atom scattering factor calculations and structural least-squares refinements, repeated until convergence. The importance of this iterative procedure is illustrated via the example of crystalline ammonia. The new HAR method is then applied to X-ray diffraction data of the dipeptide Gly-l-Ala measured at 12, 50, 100, 150, 220 and 295 K, using Hartree-Fock and BLYP density functional theory electron densities and three different basis sets. All positions and anisotropic displacement parameters (ADPs) are freely refined without constraints or restraints - even those for hydrogen atoms. The results are systematically compared with those from neutron diffraction experiments at the temperatures 12, 50, 150 and 295 K. Although non-hydrogen-atom ADPs differ by up to three combined standard uncertainties (csu's), all other structural parameters agree within less than 2 csu's. Using our best calculations (BLYP/cc-pVTZ, recommended for organic molecules), the accuracy of determining bond lengths involving hydrogen atoms from HAR is better than 0.009 Å for temperatures of 150 K or below; for hydrogen-atom ADPs it is better than 0.006 Å(2) as judged from the mean absolute X-ray minus neutron differences. These results are among the best ever obtained. Remarkably, the precision of determining bond lengths and ADPs for the hydrogen atoms from the HAR procedure is comparable with that from the neutron measurements - an outcome which is obtained with a routinely achievable resolution of the X-ray data of 0.65 Å.
Collapse
Affiliation(s)
- Silvia C. Capelli
- Department of Chemistry, University of Milan, Via Golgi 19, 20133 Milan, Italy
| | - Hans-Beat Bürgi
- Department of Chemistry and Biochemistry, Universität Bern, Freiestrasse 3, 3012 Bern, Switzerland
- Department of Chemistry, Universität Zürich, Winterthurstrasse 190, 8057 Zürich, Switzerland
| | - Birger Dittrich
- Institut für Anorganische und Angewandte Chemie, Universität Hamburg, Martin-Luther-King-Platz 6, 20146 Hamburg, Germany
- Institut für Anorganische Chemie, Georg-August-Universität Göttingen, Tammannstr. 4, 37077 Göttingen, Germany
| | - Simon Grabowsky
- School of Chemistry and Biochemistry, The University of Western Australia, 35 Stirling Highway, Crawley WA 6009, Australia
| | - Dylan Jayatilaka
- School of Chemistry and Biochemistry, The University of Western Australia, 35 Stirling Highway, Crawley WA 6009, Australia
| |
Collapse
|