1
|
Genoni A, Martín Pendás Á. Critical assessment of the x-ray restrained wave function approach: Advantages, drawbacks, and perspectives for density functional theory and periodic ab initio calculations. J Chem Phys 2024; 160:234108. [PMID: 38899684 DOI: 10.1063/5.0213247] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2024] [Accepted: 06/03/2024] [Indexed: 06/21/2024] Open
Abstract
The x-ray restrained wave function (XRW) method is a quantum crystallographic technique to extract wave functions compatible with experimental x-ray diffraction data. The approach looks for wave functions that minimize the energies of the investigated systems and also reproduce sets of x-ray structure factors. Given the strict relationship between x-ray structure factors and electron distributions, the strategy practically allows determining wave functions that correspond to given (usually experimental) electron densities. In this work, the capabilities of the XRW approach were further tested. The aim was to evaluate whether the XRW technique could serve as a tool for suggesting new exchange-correlation functionals for density functional theory or refining existing ones. Additionally, the ability of the method to address the influences of the crystalline environment was also assessed. The outcomes of XRW computations were thus compared to those of traditional gas-phase, embedding quantum mechanics/molecular mechanics, and fully periodic calculations. The results revealed that, irrespective of the initial conditions, the XRW computations practically yield a consensus electron density, in contrast to the currently employed density functional approximations (DFAs), which tend to give a too large range of electron distributions. This is encouraging in view of exploiting the XRW technique to develop improved functionals. Conversely, the calculations also emphasized that the XRW method is limited in its ability to effectively address the influences of the crystalline environment. This underscores the need for a periodic XRW technique, which would allow further untangling the shortcomings of DFAs from those inherent to the XRW approach.
Collapse
Affiliation(s)
- Alessandro Genoni
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019, 1 Boulevard Arago, 57078 Metz, France
| | - Ángel Martín Pendás
- Departamento de Química Física y Analítica, Facultad de Química, Universidad de Oviedo, Avenida Julian Clavería 8, 33006 Oviedo, Spain
| |
Collapse
|
2
|
Hasecke L, Mata RA. Optimization of Quantum Nuclei Positions with the Adaptive Nuclear-Electronic Orbital Approach. J Phys Chem A 2024; 128:3205-3211. [PMID: 38619054 PMCID: PMC11056972 DOI: 10.1021/acs.jpca.4c00096] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2024] [Revised: 03/08/2024] [Accepted: 03/08/2024] [Indexed: 04/16/2024]
Abstract
The use of multicomponent methods has become increasingly popular over the last years. Under this framework, nuclei (commonly protons) are treated quantum mechanically on the same footing as the electronic structure problem. Under the use of atomic-centered orbitals, this can lead to some complications as the ideal location of the nuclear basis centers must be optimized. In this contribution, we propose a straightforward approach to determine the position of such centers within the self-consistent cycle of a multicomponent calculation, making use of individual proton charge centroids. We test the method on model systems including the water dimer, a protonated water tetramer, and a porphine system. Comparing to numerical gradient calculations, the adaptive nuclear-electronic orbital (NEO) procedure is able to converge the basis centers to within a few cents of an Ångström and with less than 0.1 kcal/mol differences in absolute energies. This is achieved in one single calculation and with a small added computational effort of up to 80% compared to a regular NEO- self-consistent field run. An example application for the human transketolase proton wire is also provided.
Collapse
Affiliation(s)
- Lukas Hasecke
- Institute of Physical Chemistry, University
of Göttingen, Tammannstrasse 6, 37077 Göttingen, Germany
| | - Ricardo A. Mata
- Institute of Physical Chemistry, University
of Göttingen, Tammannstrasse 6, 37077 Göttingen, Germany
| |
Collapse
|
3
|
Meurer F, Kleemiss F, Riesinger C, Balázs G, Vuković V, Shenderovich IG, Jelsch C, Bodensteiner M. Probing the Isolobal Relation between Cp'''NiP 3 and White Phosphorus by Experimental Charge Density Analysis. Chemistry 2024; 30:e202303762. [PMID: 38277228 DOI: 10.1002/chem.202303762] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2023] [Revised: 12/20/2023] [Accepted: 01/26/2024] [Indexed: 01/28/2024]
Abstract
An in-depth analysis of the description of bonding within Cp'''Ni-cyclo-P3 (Cp'''=1,2,4-tri-tert-butylcyclopentadienyl, [Ni]P3) employing X-ray diffraction based multipolar modeling, density functional theory (DFT) as well as an "experimental wavefunction" obtained from X-ray restrained wavefunction (XRW) fitting is presented. The results are compared to DFT calculations on white phosphorus - an isolobal analogue to [Ni]P3. A complementary bonding analysis shows insights into the reactivity of [Ni]P3. The isolobal principle is reflected in every aspect of our analysis and the employed methods seamlessly predict the differences in reactivity of [Ni]P3 and P4. Crystallographic modeling, solid-state NMR, and DFT calculations describe the dynamic behavior of the cyclo-P3 unit in the title molecule.
Collapse
Affiliation(s)
- Florian Meurer
- Faculty for Chemistry and Pharmacy, University of Regensburg, Universitätsstraße 31, 93053, Regensburg, Germany
| | - Florian Kleemiss
- Faculty for Chemistry and Pharmacy, University of Regensburg, Universitätsstraße 31, 93053, Regensburg, Germany
- Institute for Inorganic Chemistry, RWTH Aachen, Landoltweg 1a, 52074, Aachen, Germany
| | - Christoph Riesinger
- Faculty for Chemistry and Pharmacy, University of Regensburg, Universitätsstraße 31, 93053, Regensburg, Germany
| | - Gábor Balázs
- Faculty for Chemistry and Pharmacy, University of Regensburg, Universitätsstraße 31, 93053, Regensburg, Germany
| | - Vedran Vuković
- Faculty for Chemistry and Pharmacy, University of Regensburg, Universitätsstraße 31, 93053, Regensburg, Germany
- Cristallographie, Résonance Magnétique et Modélisations CNRS, UMR 7036, Institut Jean Barriol, CNRS, Université de Lorraine BP 70239, F54506, Vandoeuvre-lès-Nancy CEDEX, France
| | - Ilya G Shenderovich
- Faculty for Chemistry and Pharmacy, University of Regensburg, Universitätsstraße 31, 93053, Regensburg, Germany
| | - Christian Jelsch
- Cristallographie, Résonance Magnétique et Modélisations CNRS, UMR 7036, Institut Jean Barriol, CNRS, Université de Lorraine BP 70239, F54506, Vandoeuvre-lès-Nancy CEDEX, France
| | - Michael Bodensteiner
- Faculty for Chemistry and Pharmacy, University of Regensburg, Universitätsstraße 31, 93053, Regensburg, Germany
| |
Collapse
|
4
|
Singh A, Avinash K, Malaspina LA, Banoo M, Alhameedi K, Jayatilaka D, Grabowsky S, Thomas SP. Dynamic Covalent Bonds in the Ebselen Class of Antioxidants Probed by X-ray Quantum Crystallography. Chemistry 2024; 30:e202303384. [PMID: 38126954 DOI: 10.1002/chem.202303384] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2023] [Revised: 12/19/2023] [Accepted: 12/20/2023] [Indexed: 12/23/2023]
Abstract
Dynamic bonds are essential structural ingredients of dynamic covalent chemistry that involve reversible cleavage and formation of bonds. Herein, we explore the electronic characteristics of Se-N bonds in the organo-selenium antioxidant ebselen and its derivatives for their propensity to function as dynamic covalent bonds by employing high-resolution X-ray quantum crystallography and complementary computational studies. An analysis of the experimentally reconstructed X-ray wavefunctions reveals the salient electronic features of the Se-N bonds with very low electron density localized at the bonding region and a positive Laplacian value at the bond critical point. Bond orders and percentage covalency and ionicity estimated from the X-ray wavefunctions, along with localized orbital locator (LOL) and electron localization function (ELF) analyses show that the Se-N bond is unique in its closed shell-like features, despite being a covalent bond. Time-dependent DFT calculations simulate the cleavage of Se-N bonds in ebselen in the excited state, further substantiating their nature as dynamic bonds.
Collapse
Affiliation(s)
- Ashi Singh
- Department of Chemistry, Indian Institute of Technology Delhi, New Delhi, 110016, India
| | - Kiran Avinash
- Department of Chemistry, Indian Institute of Technology Delhi, New Delhi, 110016, India
| | - Lorraine A Malaspina
- University of Bern, Department of Chemistry, Biochemistry and Pharmaceutical Sciences, Freiestrasse 3, 3012, Bern, Switzerland
| | - Masoumeh Banoo
- Department of Chemistry, Indian Institute of Technology Delhi, New Delhi, 110016, India
| | - Khidhir Alhameedi
- School of Molecular Sciences, University of Western Australia, Perth, WA, 6009, Australia
- Faculty of Education for Pure Sciences -, University of Kerbala, Karbala, Iraq
| | - Dylan Jayatilaka
- School of Molecular Sciences, University of Western Australia, Perth, WA, 6009, Australia
| | - Simon Grabowsky
- University of Bern, Department of Chemistry, Biochemistry and Pharmaceutical Sciences, Freiestrasse 3, 3012, Bern, Switzerland
| | - Sajesh P Thomas
- Department of Chemistry, Indian Institute of Technology Delhi, New Delhi, 110016, India
| |
Collapse
|
5
|
Pawlędzio S, Ziemniak M, Trzybiński D, Arhangelskis M, Makal A, Woźniak K. Influence of N-protonation on electronic properties of acridine derivatives by quantum crystallography. RSC Adv 2024; 14:5340-5350. [PMID: 38348299 PMCID: PMC10859733 DOI: 10.1039/d3ra08081a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/26/2023] [Accepted: 01/30/2024] [Indexed: 02/15/2024] Open
Abstract
Applications of 9-aminoacridine (9aa) and its derivatives span fields such as chemistry, biology, and medicine, including anticancer and antimicrobial activities. Protonation of such molecules can alter their bioavailability as weakly basic drugs like aminoacridines exhibit reduced solubility at high pH levels potentially limiting their effectiveness in patients with elevated gastric pH. In this study, we analyse the influence of protonation on the electronic characteristics of the molecular organic crystals of 9-aminoacridine. The application of quantum crystallography, including aspherical atom refinement, has enriched the depiction of electron density in the studied systems and non-covalent interactions, providing more details than previous studies. Our experimental results, combined with a topological analysis of the electron density and its Laplacian, provided detailed descriptions of how protonation changes the electron density distribution around the amine group and water molecule, concurrently decreasing the electron density at bond critical points of N/O-H bonds. Protonation also alters the molecular architecture of the systems under investigation. This is reflected in different proportions of the N⋯H and O⋯H intermolecular contacts for the neutral and protonated forms. Periodic DFT calculations of the cohesive energies of the crystal lattice, as well as computed interaction energies between molecules in the crystal, confirm that protonation stabilises the crystal structure due to a positive synergy between strong halogen and hydrogen bonds. Our findings highlight the potential of quantum crystallography in predicting crystal structure properties and point to its possible applications in developing new formulations for poorly soluble drugs.
Collapse
Affiliation(s)
- Sylwia Pawlędzio
- Neutron Scattering Division, Oak Ridge National Laboratory Oak Ridge TN 37831 USA
- Department of Chemistry, Biological and Chemical Research Centre, University of Warsaw Żwirki i Wigury 101 02-093 Warszawa Poland
| | - Marcin Ziemniak
- Department of Chemistry, Biological and Chemical Research Centre, University of Warsaw Żwirki i Wigury 101 02-093 Warszawa Poland
| | - Damian Trzybiński
- Department of Chemistry, Biological and Chemical Research Centre, University of Warsaw Żwirki i Wigury 101 02-093 Warszawa Poland
| | - Mihails Arhangelskis
- Department of Chemistry, Biological and Chemical Research Centre, University of Warsaw Żwirki i Wigury 101 02-093 Warszawa Poland
| | - Anna Makal
- Department of Chemistry, Biological and Chemical Research Centre, University of Warsaw Żwirki i Wigury 101 02-093 Warszawa Poland
| | - Krzysztof Woźniak
- Department of Chemistry, Biological and Chemical Research Centre, University of Warsaw Żwirki i Wigury 101 02-093 Warszawa Poland
| |
Collapse
|
6
|
Priyadharsan RR, Timothy RA, Thomas JM, Jeyakumar TC, Rajaram R, Louis H. Investigating the structure, bonding, and energy decomposition analysis of group 10 transition metal carbonyls with substituted terminal germanium chalcogenides [M(CO) 3GeX] (M = Ni, Pd, and Pt; X = O, S, Se, and Te) complexes: insight from first-principles calculations. J Mol Model 2023; 29:344. [PMID: 37847395 DOI: 10.1007/s00894-023-05745-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2023] [Accepted: 10/02/2023] [Indexed: 10/18/2023]
Abstract
CONTEXT This research focused on the theoretical investigation of transition metal carbonyls [M(CO)4] coordinated with terminal germanium chalcogenides complexes [M(CO)3GeX], where M represents Ni, Pd, and Pt and X represents O, S, Se, and Te labeled 1-15. While the notable complexes M(CO)4 (where M = Ni, Pd, Pt) numbered 1, 6, and 11 are of significance, substituting one of the CO ligands in 1, 6, and 11 with a GeX ligand (where X = O, S, Se, or Te) result in substituted complexes (2-5, 7-10, and 11-15). Substituting of the CO ligand slightly alters these bond angles. Specifically, the ∠CMC bond angles for [Ni] complexes range from 111.9° to 112.2°, for [Pd] complexes from 111.4° to 111.7°, and for [Pt] complexes from 112.4° to 112.8°. These findings indicate a minor deviation from the tetrahedral geometry due to the influence of the new GeX ligand. Similarly, there is a slight change in the geometry of the metal complexes, where the ∠GeMC angles for [Ni] complexes are between 106.7° and 106.9°, for [Pd] complexes between 107.2° and 107.5°, and for [Pt] complexes between 105.9° and 106.4°. Comparing among the substituted GeX complexes, those containing GeTe exhibit a higher natural bond orbital (NBO) contribution from the Ge atom compared to the M atom. Consequently, based on the above observations, it can be inferred that GeX acts as an effective sigma donor in contrast to carbonyl compounds. Results of energy decomposition analysis (EDA) for the M-CO bond in 1, 6, and 11 and for the M-GeX bond in the other [M(CO)3(GeX)] complexes where M = Ni, Pd and Pt. The percentage contribution of ΔEelstat and ΔEorb shows a relatively identical behavior for all ligands in case of each metal complexes. METHODS Density functional theory (DFT) calculations were conducted using the B3LYP/gen/6-31G*/LanL2DZ level of theory to examine transition metal carbonyls [M(CO)4] coordinated with terminal germanium chalcogenides complexes [M(CO)3GeX], where M represents Ni, Pd, and Pt, and X represents O, S, Se, and Te labeled 1-15 utilized through the use of Gaussian 09W and GaussView 6.0.16 software packages. Post-processing computational code such as multi-wave function was employed for results analysis and visualization.
Collapse
Affiliation(s)
- R Rameshbabu Priyadharsan
- PG & Research Department of Chemistry, The American College (Autonomous), Madurai, Tamil Nadu, India
| | - Rawlings A Timothy
- Computational and Bio-Simulation Research Group, University of Calabar, Calabar, Nigeria
| | - Jisha Mary Thomas
- Department of Chemistry, Pondicherry University, Puducherry, 605014, India
| | | | - Rajendran Rajaram
- Department of Chemistry, Madanapalle Institute of Technology and Science, Angallu (V), Madanapalle, Andhra Pradesh, 517325, India
| | - Hitler Louis
- Computational and Bio-Simulation Research Group, University of Calabar, Calabar, Nigeria.
- Centre for Herbal Pharmacology and Environmental Sustainability, Chettinad Hospital and Research Institute, Chettinad Academy of Research and Education, Kelambakkam, 603103, Tamil Nadu, India.
| |
Collapse
|
7
|
Genoni A. Elucidating the nature of chemical bonds in a coordination compound through quantum crystallographic techniques. ACTA CRYSTALLOGRAPHICA SECTION B, STRUCTURAL SCIENCE, CRYSTAL ENGINEERING AND MATERIALS 2023; 79:253-254. [PMID: 37561073 DOI: 10.1107/s2052520623006364] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/21/2023] [Accepted: 07/21/2023] [Indexed: 08/11/2023]
Abstract
Investigations simultaneously involving multiple techniques of quantum crystallography could be very useful to prove the consistency of obtained results or to highlight different facets of the same scientific phenomenon or problem. Pinto et al. [Acta Cryst. (2023), B79, 282-296] exploit three different quantum crystallographic techniques (Hansen & Coppens multipole model refinement, QTAIM analysis of the electron density, and Hirshfeld atom refinement) to characterize the nature of chemical bonds and of intra/intermolecular interactions in an organometallic compound.
Collapse
Affiliation(s)
- Alessandro Genoni
- Université de Lorraine and CNRS, Laboratoire de Physique et Chimie Théoriques, 1 Boulevard Arago, Metz, F-57070, France
| |
Collapse
|
8
|
Fokt I, Cybulski M, Skora S, Pająk B, Ziemniak M, Woźniak K, Zielinski R, Priebe W. d-Glucose- and d-mannose-based antimetabolites. Part 4: Facile synthesis of mono- and di-acetates of 2-deoxy-d-glucose prodrugs as potentially useful antimetabolites. Carbohydr Res 2023; 531:108861. [PMID: 37356236 DOI: 10.1016/j.carres.2023.108861] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2023] [Revised: 05/29/2023] [Accepted: 06/05/2023] [Indexed: 06/27/2023]
Abstract
2-Deoxy-d-glucose (2-DG), a compound known to interfere with d-glucose and d-mannose metabolism, has been tested as a potential anticancer and antiviral agent. Preclinical and clinical studies focused on 2-DG have highlighted several limitations related to 2-DG drug-like properties, such as poor pharmacokinetic properties. To overcome this problem, we proposed design and synthesis of novel 2-DG prodrugs that subsequently could be tested using a variety of biochemical and molecular methods. We narrowed here our focus to esters of 2-DG as potential prodrugs based on the hypothesis that ubiquitous esterases will regenerate 2-DG, leading to increased circulation time of drug and adequate organ and tumor penetration. Testing this hypothesis in vitro and, especially, in vivo requires significant amounts of respective pure mono- and previously unknown di-acetylated water-soluble derivatives of 2-DG. Development of their efficient and practical method of synthesis was imperative. We describe novel facile and scalable syntheses of seven selectively acetylated water-soluble derivatives of 2-DG and present a detailed 1H and 13C NMR analysis of all final products. X-ray diffraction analysis has been performed for compound WP1122 that was selected for detailed preclinical and subsequent clinical evaluation as potential anticancer or antiviral agent.
Collapse
Affiliation(s)
- Izabela Fokt
- Department of Experimental Therapeutics, The University of Texas MD Anderson Cancer Center, Houston, TX, 77030, USA
| | - Marcin Cybulski
- Department of Experimental Therapeutics, The University of Texas MD Anderson Cancer Center, Houston, TX, 77030, USA; Lukasiewicz-Industrial Chemistry Institute, Rydygiera 8, 01-793, Warsaw, Poland
| | - Stanisław Skora
- Department of Experimental Therapeutics, The University of Texas MD Anderson Cancer Center, Houston, TX, 77030, USA
| | - Beata Pająk
- Independent Laboratory of Genetics and Molecular Biology, Kaczkowski Military Institute of Hygiene and Epidemiology, Kozielska 4, 01-163, Warsaw, Poland
| | - Marcin Ziemniak
- Department of Chemistry, Biological and Chemical Research Centre, University of Warsaw, Żwirki I Wigury 101, 02-089, Warsaw, Poland
| | - Krzysztof Woźniak
- Department of Chemistry, Biological and Chemical Research Centre, University of Warsaw, Żwirki I Wigury 101, 02-089, Warsaw, Poland
| | - Rafal Zielinski
- Department of Experimental Therapeutics, The University of Texas MD Anderson Cancer Center, Houston, TX, 77030, USA
| | - Waldemar Priebe
- Department of Experimental Therapeutics, The University of Texas MD Anderson Cancer Center, Houston, TX, 77030, USA.
| |
Collapse
|
9
|
Wuttig M, Schön CF, Lötfering J, Golub P, Gatti C, Raty JY. Revisiting the Nature of Chemical Bonding in Chalcogenides to Explain and Design their Properties. ADVANCED MATERIALS (DEERFIELD BEACH, FLA.) 2023; 35:e2208485. [PMID: 36456187 DOI: 10.1002/adma.202208485] [Citation(s) in RCA: 11] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/15/2022] [Revised: 10/31/2022] [Indexed: 05/19/2023]
Abstract
Quantum chemical bonding descriptors have recently been utilized to design materials with tailored properties. Their usage to facilitate a quantitative description of bonding in chalcogenides as well as the transition between different bonding mechanisms is reviewed. More importantly, these descriptors can also be employed as property predictors for several important material characteristics, including optical and transport properties. Hence, these quantum chemical bonding descriptors can be utilized to tailor material properties of chalcogenides relevant for thermoelectrics, photovoltaics, and phase-change memories. Relating material properties to bonding mechanisms also shows that there is a class of materials, which are characterized by unconventional properties such as a pronounced anharmonicity, a large chemical bond polarizability, and strong optical absorption. This unusual property portfolio is attributed to a novel bonding mechanism, fundamentally different from ionic, metallic, and covalent bonding, which is called "metavalent." In the concluding section, a number of promising research directions are sketched, which explore the nature of the property changes upon changing bonding mechanism and extend the concept of quantum chemical property predictors to more complex compounds.
Collapse
Affiliation(s)
- Matthias Wuttig
- I. Institute of Physics, Physics of Novel Materials, RWTH Aachen University, 52056, Aachen, Germany
- Jülich-Aachen Research Alliance (JARA FIT and JARA HPC), RWTH Aachen University, 52056, Aachen, Germany
- PGI 10 (Green IT), Forschungszentrum Jülich GmbH, 52428, Jülich, Germany
| | - Carl-Friedrich Schön
- I. Institute of Physics, Physics of Novel Materials, RWTH Aachen University, 52056, Aachen, Germany
| | - Jakob Lötfering
- I. Institute of Physics, Physics of Novel Materials, RWTH Aachen University, 52056, Aachen, Germany
| | - Pavlo Golub
- Department of Theoretical Chemistry, J. Heyrovský Institute of Physical Chemistry, Dolejškova 2155/3, Prague 8, 182 23, Czech Republic
| | - Carlo Gatti
- CNR-SCITEC, Istituto di Scienze e Tecnologie Chimiche "Giulio Natta", sezione di via Golgi, via Golgi 19, Milano, 20133, Italy
| | - Jean-Yves Raty
- CESAM B5, Université de Liège, Sart-Tilman, B4000, Belgium
| |
Collapse
|
10
|
Wieduwilt EK, Boto RA, Macetti G, Laplaza R, Contreras-García J, Genoni A. Extracting Quantitative Information at Quantum Mechanical Level from Noncovalent Interaction Index Analyses. J Chem Theory Comput 2023; 19:1063-1079. [PMID: 36656682 DOI: 10.1021/acs.jctc.2c01092] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
The noncovalent interaction (NCI) index is nowadays a well-known strategy to detect NCIs in molecular systems. Even though it initially provided only qualitative descriptions, the technique has been recently extended to also extract quantitative information. To accomplish this task, integrals of powers of the electron distribution were considered, with the requirement that the overall electron density can be clearly decomposed as sum of distinct fragment contributions to enable the definition of the (noncovalent) integration region. So far, this was done by only exploiting approximate promolecular electron densities, which are given by the sum of spherically averaged atomic electron distributions and thus represent too crude approximations. Therefore, to obtain more quantum mechanically (QM) rigorous results from NCI index analyses, in this work, we propose to use electron densities obtained through the transfer of extremely localized molecular orbitals (ELMOs) or through the recently developed QM/ELMO embedding technique. Although still approximate, the electron distributions resulting from the abovementioned methods are fully QM and, above all, are again partitionable into subunit contributions, which makes them completely suitable for the NCI integral approach. Therefore, we benchmarked the integrals resulting from NCI index analyses (both those based on the promolecular densities and those based on ELMO electron distributions) against interaction energies computed at a high quantum chemical level (in particular, at the coupled cluster level). The performed test calculations have indicated that the NCI integrals based on ELMO electron densities outperform the promolecular ones. Furthermore, it was observed that the novel quantitative NCI-(QM/)ELMO approach can be also profitably exploited both to characterize and evaluate the strength of specific interactions between ligand subunits and protein residues in protein-ligand complexes and to follow the evolution of NCIs along trajectories of molecular dynamics simulations. Although further methodological improvements are still possible, the new quantitative ELMO-based technique could be already exploited in situations in which fast and reliable assessments of NCIs are crucial, such as in computational high-throughput screenings for drug discovery.
Collapse
Affiliation(s)
- Erna K Wieduwilt
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019, 1 Boulevard Arago, Metz F-57078, France
| | - Roberto A Boto
- Laboratoire de Chimie Théorique (LCT), UMR 7616, Sorbonne Université & CNRS, 4 Place Jussieu, Paris F-75005, France
| | - Giovanni Macetti
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019, 1 Boulevard Arago, Metz F-57078, France
| | - Rubén Laplaza
- Laboratoire de Chimie Théorique (LCT), UMR 7616, Sorbonne Université & CNRS, 4 Place Jussieu, Paris F-75005, France
| | - Julia Contreras-García
- Laboratoire de Chimie Théorique (LCT), UMR 7616, Sorbonne Université & CNRS, 4 Place Jussieu, Paris F-75005, France
| | - Alessandro Genoni
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019, 1 Boulevard Arago, Metz F-57078, France
| |
Collapse
|
11
|
Macetti G, Genoni A. Introduction of a weighting scheme for the X-ray restrained wavefunction approach: advantages and drawbacks. Acta Crystallogr A Found Adv 2023; 79:25-40. [PMID: 36601761 DOI: 10.1107/s2053273322010221] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/04/2022] [Accepted: 10/23/2022] [Indexed: 11/22/2022] Open
Abstract
In a quite recent study [Genoni et al. (2017). IUCrJ, 4, 136-146], it was observed that the X-ray restrained wavefunction (XRW) approach allows a more efficient and larger capture of electron correlation effects on the electron density if high-angle reflections are not considered in the calculations. This is due to the occurrence of two concomitant effects when one uses theoretical X-ray diffraction data corresponding to a single-molecule electron density in a large unit cell: (i) the high-angle reflections are generally much more numerous than the low- and medium-angle ones, and (ii) they are already very well described at unrestrained level. Nevertheless, since high-angle data also contain important information that should not be disregarded, it is not advisable to neglect them completely. For this reason, based on the results of the previous investigation, this work introduces a weighting scheme for XRW calculations to up-weight the contribution of low- and medium-angle reflections, and, at the same time, to reasonably down-weight the importance of the high-angle data. The proposed strategy was tested through XRW computations with both theoretical and experimental structure-factor amplitudes. The tests have shown that the new weighting scheme works optimally if it is applied with theoretically generated X-ray diffraction data, while it is not advantageous when traditional experimental X-ray diffraction data (even of very high resolution) are employed. This also led to the conclusion that the use of a specific external parameter λJ for each resolution range might not be a suitable strategy to adopt in XRW calculations exploiting experimental X-ray data as restraints.
Collapse
Affiliation(s)
- Giovanni Macetti
- Université de Lorraine and CNRS, Laboratoire de Physique et Chimie Théoriques, 1 Boulevard Arago, Metz, F-57078, France
| | - Alessandro Genoni
- Université de Lorraine and CNRS, Laboratoire de Physique et Chimie Théoriques, 1 Boulevard Arago, Metz, F-57078, France
| |
Collapse
|
12
|
Hill NDD, Lilienthal E, Bender CO, Boeré RT. Accurate Crystal Structures of C 12H 9CN, C 12H 8(CN) 2, and C 16H 11CN Valence Isomers Using Nonspherical Atomic Scattering Factors. J Org Chem 2022; 87:16213-16229. [DOI: 10.1021/acs.joc.2c01706] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
Affiliation(s)
- Nathan D. D. Hill
- Department of Chemistry and Biochemistry, University of Lethbridge, 4401 University Dr. W, Lethbridge, AB, Canada T1K 3M4
- The Canadian Centre for Research in Advanced Fluorine Technologies (C-CRAFT), University of Lethbridge, 4401 University Dr. W, Lethbridge, AB, Canada T1K 3M4
| | - Elaura Lilienthal
- Department of Chemistry and Biochemistry, University of Lethbridge, 4401 University Dr. W, Lethbridge, AB, Canada T1K 3M4
- The Canadian Centre for Research in Advanced Fluorine Technologies (C-CRAFT), University of Lethbridge, 4401 University Dr. W, Lethbridge, AB, Canada T1K 3M4
| | - Christopher O. Bender
- Department of Chemistry and Biochemistry, University of Lethbridge, 4401 University Dr. W, Lethbridge, AB, Canada T1K 3M4
| | - René T. Boeré
- Department of Chemistry and Biochemistry, University of Lethbridge, 4401 University Dr. W, Lethbridge, AB, Canada T1K 3M4
- The Canadian Centre for Research in Advanced Fluorine Technologies (C-CRAFT), University of Lethbridge, 4401 University Dr. W, Lethbridge, AB, Canada T1K 3M4
| |
Collapse
|
13
|
Leherte L, Bodart L, Wouters J, Vercauteren DP. Description of non-covalent interactions in benzyl chalcocyanate crystals from smoothed Cromer-Mann electron density distribution functions. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2022; 34:494003. [PMID: 36223781 DOI: 10.1088/1361-648x/ac99c9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/25/2022] [Accepted: 10/12/2022] [Indexed: 06/16/2023]
Abstract
A well-known method to characterize non-covalent interactions consists in the topological analysis of electron density distribution (EDD) functions, complemented by the search for minima in the reduced density gradient (RDG) distributions. Here, we characterize intermolecular interactions occurring in crystals of benzyl chalcocyanate compounds through bond critical points (BCP) of the promolecular electron density (ED) built from the crystallographic Cromer-Mann parameters, at several smoothing levelst. The trajectories formed by thet-dependent BCP locations are interpreted in terms of the intermolecular interactions occurring within the crystal arrangements. Chalcogen…nitro BCPs are clearly present in the unsmoothed EDDs but are annihilated astincreases, while chalcogen…chalcogen BCPs appear and are among the only BCPs left at the highest smoothing level. The chalcogen bonds are differentiated from the other chalcogen interactions through the linear chalcogen…BCP…nitro geometry at low smoothing level and their more negative Laplacian values. The annihilation of CPs can be followed by the apparition of a RDG minimum, associated with a very weak interaction. Along the BCP trajectories, the Laplacian shows a progressive concentration of the ED in the intermolecular space within the crystals and adopts the most negative values at the shortest atom…atom separations. At the termination point of a BCP trajectory, the drastic increase of the ellipticity value illustrates the flattening of the EDD.
Collapse
Affiliation(s)
- Laurence Leherte
- Laboratory of Structural Biological Chemistry, Unit of Theoretical and Structural Physical Chemistry, Department of Chemistry, NAmur Research Institute for LIfe Sciences (NARILIS), Namur Institute of Structured Matter (NISM), NAmur MEdicine and Drug Innovation Center (NAMEDIC), University of Namur, Rue de Bruxelles 61, B-5000 Namur, Belgium
| | - Laurie Bodart
- Laboratory of Structural Biological Chemistry, Unit of Theoretical and Structural Physical Chemistry, Department of Chemistry, NAmur Research Institute for LIfe Sciences (NARILIS), Namur Institute of Structured Matter (NISM), NAmur MEdicine and Drug Innovation Center (NAMEDIC), University of Namur, Rue de Bruxelles 61, B-5000 Namur, Belgium
| | - Johan Wouters
- Laboratory of Structural Biological Chemistry, Unit of Theoretical and Structural Physical Chemistry, Department of Chemistry, NAmur Research Institute for LIfe Sciences (NARILIS), Namur Institute of Structured Matter (NISM), NAmur MEdicine and Drug Innovation Center (NAMEDIC), University of Namur, Rue de Bruxelles 61, B-5000 Namur, Belgium
| | - Daniel P Vercauteren
- Laboratory of Structural Biological Chemistry, Unit of Theoretical and Structural Physical Chemistry, Department of Chemistry, NAmur Research Institute for LIfe Sciences (NARILIS), Namur Institute of Structured Matter (NISM), NAmur MEdicine and Drug Innovation Center (NAMEDIC), University of Namur, Rue de Bruxelles 61, B-5000 Namur, Belgium
| |
Collapse
|
14
|
Electron density is not spherical: the many applications of the transferable aspherical atom model. Comput Struct Biotechnol J 2022; 20:6237-6243. [DOI: 10.1016/j.csbj.2022.10.018] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2022] [Revised: 10/12/2022] [Accepted: 10/13/2022] [Indexed: 11/20/2022] Open
|
15
|
Shteingolts SA, Stash AI, Tsirelson VG, Fayzullin RR. Real-Space Interpretation of Interatomic Charge Transfer and Electron Exchange Effects by Combining Static and Kinetic Potentials and Associated Vector Fields. Chemistry 2022; 28:e202200985. [PMID: 35638164 DOI: 10.1002/chem.202200985] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2022] [Indexed: 11/09/2022]
Abstract
Intricate behaviour of one-electron potentials from the Euler equation for electron density and corresponding gradient force fields in crystals was studied. Channels of locally enhanced kinetic potential and corresponding saddle Lagrange points were found between chemically bonded atoms. Superposition of electrostatic ϕ e s r and kinetic ϕ k r potentials and electron density ρ r allowed partitioning any molecules and crystals into atomic ρ - and potential-based ϕ -basins; ϕ k -basins explicitly account for the electron exchange effect, which is missed for ϕ e s -ones. Phenomena of interatomic charge transfer and related electron exchange were explained in terms of space gaps between zero-flux surfaces of ρ - and ϕ -basins. The gap between ϕ e s - and ρ -basins represents the charge transfer, while the gap between ϕ k - and ρ -basins is a real-space manifestation of sharing the transferred electrons caused by the static exchange and kinetic effects as a response against the electron transfer. The regularity describing relative positions of ρ -, ϕ e s -, and ϕ k - basin boundaries between interacting atoms was proposed. The position of ϕ k -boundary between ϕ e s - and ρ -ones within an electron occupier atom determines the extent of transferred electron sharing. The stronger an H⋅⋅⋅O hydrogen bond is, the deeper hydrogen atom's ϕ k -basin penetrates oxygen atom's ρ -basin, while for covalent bonds a ϕ k -boundary closely approaches a ϕ e s -one indicating almost complete sharing of the transferred electrons. In the case of ionic bonds, the same region corresponds to electron pairing within the ρ -basin of an electron occupier atom.
Collapse
Affiliation(s)
- Sergey A Shteingolts
- Arbuzov Institute of Organic and Physical Chemistry, FRC Kazan Scientific Center, Russian Academy of Sciences, 8 Arbuzov Street, Kazan, 420088, Russian Federation
| | - Adam I Stash
- A.N. Nesmeyanov Institute of Organoelement Compounds, Russian Academy of Sciences, 28 Vavilov Street, Moscow, 119991, Russian Federation
| | - Vladimir G Tsirelson
- D.I. Mendeleev University of Chemical Technology, 9 Miusskaya Square, Moscow, 125047, Russian Federation.,South Ural State University, 76 Lenin Avenue, Chelyabinsk, 454080, Russian Federation
| | - Robert R Fayzullin
- Arbuzov Institute of Organic and Physical Chemistry, FRC Kazan Scientific Center, Russian Academy of Sciences, 8 Arbuzov Street, Kazan, 420088, Russian Federation
| |
Collapse
|
16
|
Genoni A. On the termination of the X-ray constrained wavefunction procedure: reformulation of the method for an unequivocal determination of λ. Acta Crystallogr A Found Adv 2022; 78:302-308. [DOI: 10.1107/s2053273322003746] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2022] [Accepted: 04/05/2022] [Indexed: 11/10/2022] Open
Abstract
The X-ray constrained/restrained wavefunction (XCW/XRW) approach of quantum crystallography is revisited by introducing the stationary condition of the Jayatilaka functional with respect to the Lagrange multiplier λ. The theoretical derivation has unequivocally shown that the right value of λ is a maximum stationary point of the functional to optimize, thus enabling the solution of the longstanding problem of establishing the point at which to halt the XCW/XRW procedure. Based on the new finding, a reformulation of the X-ray constrained wavefunction algorithm is proposed and its implementation is envisaged. In addition to relying on more solid mathematical grounds, the new variant of the method will be intrinsically more physically meaningful, allowing a straightforward evaluation of the highest level of confidence with which the experimental X-ray diffraction data can be possibly reproduced.
Collapse
|
17
|
Hoser AA, Rekis T, Madsen AØ. Dynamics and disorder: on the stability of pyrazinamide polymorphs. ACTA CRYSTALLOGRAPHICA SECTION B STRUCTURAL SCIENCE, CRYSTAL ENGINEERING AND MATERIALS 2022; 78:416-424. [PMID: 35695115 PMCID: PMC9254588 DOI: 10.1107/s2052520622004577] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/11/2022] [Accepted: 04/29/2022] [Indexed: 11/10/2022]
Abstract
The enantiotropic relationship between the four polymorphs of pyrazinamide is analyzed by means of accurate X-ray diffraction measurements, normal-mode refinement and periodic DFT calculations. This article focuses on the structure and relative stability of four pyrazinamide polymorphs. New single crystal X-ray diffraction data collected for all forms at 10 K and 122 K are presented. By combining periodic abinitio DFT calculations with normal-mode refinement against X-ray diffraction data, both enthalpic and entropic contributions to the free energy of all polymorphs are calculated. On the basis of the estimated free energies, the stability order of the polymorphs as a function of temperature and the corresponding solid state phase transition temperatures are anticipated. It can be concluded that the α and γ forms have higher vibrational entropy than that of the β and δ forms and therefore they are significantly more stabilized at higher temperatures. Due to the entropy which arises from the disorder in γ form, it overcomes form α and is the most stable form at temperatures above ∼500 K. Our findings are in qualitative agreement with the experimental calorimetry results.
Collapse
|
18
|
Davidson ML, Grabowsky S, Jayatilaka D. X-ray constrained wavefunctions based on Hirshfeld atoms. I. Method and review. ACTA CRYSTALLOGRAPHICA SECTION B, STRUCTURAL SCIENCE, CRYSTAL ENGINEERING AND MATERIALS 2022; 78:312-332. [PMID: 35695105 DOI: 10.1107/s2052520622004097] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/07/2022] [Accepted: 04/18/2022] [Indexed: 06/15/2023]
Abstract
The X-ray constrained wavefunction (XCW) procedure for obtaining an experimentally reconstructed wavefunction from X-ray diffraction data is reviewed. The two-center probability distribution model used to perform nuclear-position averaging in the original paper [Grimwood & Jayatilaka (2001). Acta Cryst. A57, 87-100] is carefully distinguished from the newer one-center probability distribution model. In the one-center model, Hirshfeld atoms are used, and the Hirshfeld atom based X-ray constrained wavefunction (HA-XCW) procedure is described for the first time, as well as its efficient implementation. In this context, the definition of the related X-ray wavefunction refinement (XWR) method is refined. The key halting problem for the XCW method - the procedure by which one determines when overfitting has occurred - is named and work on it reviewed.
Collapse
Affiliation(s)
- Max L Davidson
- School of Molecular Sciences, The University of Western Australia, 35 Stirling Highway, Crawley 6009, Western Australia, Australia
| | - Simon Grabowsky
- Departement für Chemie, Biochemie und Pharmazie, Universität Bern, Freiestrasse 3, 3012 Bern, Switzerland
| | - Dylan Jayatilaka
- School of Molecular Sciences, The University of Western Australia, 35 Stirling Highway, Crawley 6009, Western Australia, Australia
| |
Collapse
|
19
|
Pal R, Jelsch C, Momma K, Grabowsky S. π-Hole bonding in a new co-crystal hydrate of gallic acid and pyrazine: static and dynamic charge density analysis. ACTA CRYSTALLOGRAPHICA SECTION B, STRUCTURAL SCIENCE, CRYSTAL ENGINEERING AND MATERIALS 2022; 78:231-246. [PMID: 35411861 PMCID: PMC9004022 DOI: 10.1107/s2052520622001457] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 10/02/2021] [Accepted: 02/07/2022] [Indexed: 06/14/2023]
Abstract
A new cocrystal hydrate of gallic acid with pyrazine (4GA, Py, 4H2O; GA4PyW4) was obtained and characterized by single crystal X-ray diffraction. In addition to structure determination, experimental charge density analysis was carried out in terms of Multipole Modelling (MP), X-ray wavefunction refinement (XWR) and maximum entropy method (MEM). As a part of XWR, the structural refinement via Hirshfeld atom refinement was carried out and resulted in O-H bond lengths close to values from neutron diffraction. A systematic comparison of molecular conformations and aromatic interactions in this new cocrystal hydrate was performed with other existing polymorphs of gallic acid. In GA4PyW4, the two symmetry-independent gallic acid molecules have a syn COOH orientation and form the common (COOH)2 dimeric synthon. The carboxyl C atom displays the characteristics of π-holes with electropositive regions above and below the molecular plane and engages in acceptor-donor interactions with oxygen atoms of acidic O-H groups and phenol groups of neighbouring gallic acid molecules. The signature of the π-hole was identified from experimental charge density analysis, both in static density maps in MP and XWR as well as dynamic density in MEM, but it cannot be pinned down to a specific atom-atom interaction. This study presents the first comparison between an XWR and a MEM experimental electron-density determination.
Collapse
Affiliation(s)
- Rumpa Pal
- Faculty of Pure and Applied Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki 305-8571, Japan
- Institute of Inorganic Chemistry and Crystallography, Department 2 – Biology/Chemistry, University of Bremen, Leobener Str. 3, 28359 Bremen, Germany
| | | | - Koichi Momma
- National Museum of Nature and Science, 4-1-1 Amakubo, Tsukuba, Ibaraki, Japan
| | - Simon Grabowsky
- Institute of Inorganic Chemistry and Crystallography, Department 2 – Biology/Chemistry, University of Bremen, Leobener Str. 3, 28359 Bremen, Germany
- Department of Chemistry, Biochemistry and Pharmaceutical Sciences, University of Bern, Freiestrasse 3, 3012 Bern, Switzerland
| |
Collapse
|
20
|
Pal R, Jayatilaka D, Nishibori E. Structure Factors and Charge Density Description of Aluminum: A Quantum Crystallographic Study. J Phys Chem A 2022; 126:2042-2049. [PMID: 35316043 DOI: 10.1021/acs.jpca.1c10730] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Static structure factors and charge density for metallic aluminum were investigated by periodic calculations using atom-centered Gaussian-type basis sets with the Perdew-Burke-Ernzerhof (PBE) functional implemented in the CRYSTAL14 package and X-ray constrained wave function (XCW) fitting. The effects of additional diffuse d and f basis functions on structure factors were compared with synchrotron powder X-ray diffraction and quantitative convergent electron beam diffraction data. Changes in structure factors from an independent atom model at 022, 113, and 222 reflections introduced d and f basis functions similar to those of the experimental data. The XCW fitting was applied to different sizes of aluminum clusters. The charge density features for a 50-atom cluster clearly demonstrated electron accumulation at tetrahedral sites and electron depletion at octahedral sites. The resolution dependence of the XCW study suggests that structure factors of the five lowest resolution reflections with 0.1% accuracy were indispensable for determining the detailed bonding description in the case of metallic aluminum.
Collapse
Affiliation(s)
- Rumpa Pal
- Faculty of Pure and Applied Sciences, University of Tsukuba, Tsukuba, Ibaraki 305-8571, Japan
| | - Dylan Jayatilaka
- School of Molecular Sciences, University of Western Australia, 35 Stirling Highway, Perth, Washington 6009, Australia
| | - Eiji Nishibori
- Faculty of Pure and Applied Sciences, University of Tsukuba, Tsukuba, Ibaraki 305-8571, Japan.,Tsukuba Research Center for Energy Materials Science, University of Tsukuba, Tsukuba 305-8571, Japan
| |
Collapse
|
21
|
Pawlȩdzio S, Malinska M, Kleemiss F, Grabowsky S, Woźniak K. Aurophilic Interactions Studied by Quantum Crystallography. Inorg Chem 2022; 61:4235-4239. [PMID: 35230099 PMCID: PMC8924918 DOI: 10.1021/acs.inorgchem.1c03333] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
This is the first
use of a wave-function-based crystallographic
method to characterize aurophilic interactions from X-ray diffraction
data. Theoretical calculations previously suggested the importance
of electron correlation and dispersion forces, but no influence of
relativistic corrections to the Au...Au interaction energy was found.
In this study, we confirm the importance of relativistic corrections
in the characterization of aurophilic interactions in addition to
electron correlation and dispersion. Hirshfeld
atom refinement was used to characterize aurophilic
interactions from X-ray diffraction data. An intermediate closed-shell
type of aurophilic interaction with some features of covalency was
identified when both electron correlation and relativistic corrections
were applied. Relativistic correction changes the electron density
distribution more than electron correlation. Relativistic effects
strongly dominate the metal core region also in the direction of the
noncovalent interactions and all of the valence and bonding regions
with regard to the Au···Au interaction.
Collapse
Affiliation(s)
- Sylwia Pawlȩdzio
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury 101, 02-089 Warszawa, Poland
| | - Maura Malinska
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury 101, 02-089 Warszawa, Poland
| | - Florian Kleemiss
- Department of Chemistry, Biochemistry and Pharmaceutical Sciences, University of Bern, Freiestrasse 3, 3012 Bern, Switzerland.,Faculty for Chemistry und Pharmacy, University of Regensburg, Universitätsstrasse 31, 93053 Regensburg, Germany
| | - Simon Grabowsky
- Department of Chemistry, Biochemistry and Pharmaceutical Sciences, University of Bern, Freiestrasse 3, 3012 Bern, Switzerland
| | - Krzysztof Woźniak
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury 101, 02-089 Warszawa, Poland
| |
Collapse
|
22
|
Matta CF, Huang L, Massa L. Quantum Crystallography:
N
‐Representability Big and Small**. Isr J Chem 2021. [DOI: 10.1002/ijch.202100108] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
- Chérif F. Matta
- Department of Chemistry and Physics Mount Saint Vincent University Halifax Nova Scotia Canada B3M 2J6
- Department of Chemistry Dalhousie University Halifax Nova Scotia Canada B3H 4J3
- Department of Chemistry Saint Mary's University Halifax Nova Scotia Canada B3H 3C3
- Département de chimie Université Laval, Québec Québec Canada G1V 0A6
| | - Lulu Huang
- Department of Chemistry Hunter College City University of New York NY 10065 USA
| | - Lou Massa
- Department of Chemistry Hunter College City University of New York NY 10065 USA
- Departments of Chemistry and Physics Graduate Center City University of New York NY 10016 USA
| |
Collapse
|
23
|
Wenger JS, Wang X, Johnstone TC. H-Atom Assignment and Sb-O Bonding of [Mes 3SbOH][O 3SPh] Confirmed by Neutron Diffraction, Multipole Modeling, and Hirshfeld Atom Refinement. Inorg Chem 2021; 60:16048-16052. [PMID: 34661394 DOI: 10.1021/acs.inorgchem.1c02229] [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
Neutron wavelength-resolved Laue diffraction experiments permit accurate refinement of the H-atom positions and anisotropic displacement parameters of [Mes3SbOH][O3SPh]. A multipole-based charge density refinement and a topological analysis of the refined electron density were also performed. Hirshfeld atom refinement (HAR) recovers the neutron-determined H-atom parameters, and the quantum-mechanical electron density used in HAR recovers the electron density topology from the refined multipole model. These results confirm that [Mes3SbOH][O3SPh] does indeed feature a hydroxystibonium cation with a nominal Sb-O single bond and not a stibine oxide with an Sb=O/Sb+-O- bond.
Collapse
Affiliation(s)
- John S Wenger
- Department of Chemistry and Biochemistry, University of California Santa Cruz, Santa Cruz, California 95064, United States
| | - Xiaoping Wang
- Neutron Scattering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| | - Timothy C Johnstone
- Department of Chemistry and Biochemistry, University of California Santa Cruz, Santa Cruz, California 95064, United States
| |
Collapse
|
24
|
Midgley L, Bourhis LJ, Dolomanov OV, Grabowsky S, Kleemiss F, Puschmann H, Peyerimhoff N. Vanishing of the atomic form factor derivatives in non-spherical structural refinement - a key approximation scrutinized in the case of Hirshfeld atom refinement. ACTA CRYSTALLOGRAPHICA A-FOUNDATION AND ADVANCES 2021; 77:519-533. [PMID: 34726630 DOI: 10.1107/s2053273321009086] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/10/2021] [Accepted: 09/02/2021] [Indexed: 11/10/2022]
Abstract
When calculating derivatives of structure factors, there is one particular term (the derivatives of the atomic form factors) that will always be zero in the case of tabulated spherical atomic form factors. What happens if the form factors are non-spherical? The assumption that this particular term is very close to zero is generally made in non-spherical refinements (for example, implementations of Hirshfeld atom refinement or transferable aspherical atom models), unless the form factors are refinable parameters (for example multipole modelling). To evaluate this general approximation for one specific method, a numerical differentiation was implemented within the NoSpherA2 framework to calculate the derivatives of the structure factors in a Hirshfeld atom refinement directly as accurately as possible, thus bypassing the approximation altogether. Comparing wR2 factors and atomic parameters, along with their uncertainties from the approximate and numerically differentiating refinements, it turns out that the impact of this approximation on the final crystallographic model is indeed negligible.
Collapse
Affiliation(s)
- Laura Midgley
- Department of Mathematical Sciences, Durham University, Upper Mountjoy Campus, Stockton Road, Durham, DH1 3LE, United Kingdom
| | - Luc J Bourhis
- Bruker, 4 Allée Lorentz, Champs-sur-Marne, 77447 Marne-la-Vallée cedex 2, France
| | - Oleg V Dolomanov
- OlexSys Ltd, Department of Chemistry, Durham University, South Road, Durham, DH1 3LE, United Kingdom
| | - Simon Grabowsky
- Departement für Chemie, Biochemie und Pharmazie, Universität Bern, Freiestrasse 3, 3012 Bern, Switzerland
| | - Florian Kleemiss
- Departement für Chemie, Biochemie und Pharmazie, Universität Bern, Freiestrasse 3, 3012 Bern, Switzerland
| | - Horst Puschmann
- OlexSys Ltd, Department of Chemistry, Durham University, South Road, Durham, DH1 3LE, United Kingdom
| | - Norbert Peyerimhoff
- Department of Mathematical Sciences, Durham University, Upper Mountjoy Campus, Stockton Road, Durham, DH1 3LE, United Kingdom
| |
Collapse
|
25
|
Saunders LK, Pallipurath AR, Gutmann MJ, Nowell H, Zhang N, Allan DR. A quantum crystallographic approach to short hydrogen bonds. CrystEngComm 2021; 23:6180-6190. [PMID: 34588923 PMCID: PMC8436739 DOI: 10.1039/d1ce00355k] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2021] [Accepted: 07/27/2021] [Indexed: 11/21/2022]
Abstract
In this work we use high-resolution synchrotron X-ray diffraction for electron density mapping, in conjunction with ab initio modelling, to study short O-H⋯O and O+-H⋯O- hydrogen bonds whose behaviour is known to be tuneable by temperature. The short hydrogen bonds have donor-acceptor distances in the region of 2.45 Å and are formed in substituted urea and organic acid molecular complexes of N,N'-dimethylurea oxalic acid 2 : 1 (1), N,N-dimethylurea 2,4-dinitrobenzoate 1 : 1 (2) and N,N-dimethylurea 3,5-dinitrobenzoic acid 2 : 2 (3). From the combined analyses, these complexes are found to fall within the salt-cocrystal continuum and exhibit short hydrogen bonds that can be characterised as both strong and electrostatic (1, 3) or very strong with a significant covalent contribution (2). An additional charge assisted component is found to be important in distinguishing the relatively uncommon O-H⋯O pseudo-covalent interaction from a typical strong hydrogen bond. The electron density is found to be sensitive to the extent of static proton transfer, presenting it as a useful parameter in the study of the salt-cocrystal continuum. From complementary calculated hydrogen atom potentials, we attribute changes in proton position to the molecular environment. Calculated potentials also show zero barrier to proton migration, forming an 'energy slide' between the donor and acceptor atoms. The better fundamental understanding of the short hydrogen bond in the 'zone of fluctuation' presented in a salt-cocrystal continuum, enabled by studies like this, provide greater insight into their related properties and can have implications in the regulation of pharmaceutical materials.
Collapse
Affiliation(s)
- Lucy K Saunders
- Diamond Light Source, Harwell Science and Innovation Campus Didcot OX11 0DE UK
| | - Anuradha R Pallipurath
- School of Chemical and Process Engineering, University of Leeds Leeds LS2 9JT UK
- Research Complex at Harwell Didcot Oxfordshire OX11 0DE UK
- EPSRC Centre for Innovative Manufacturing in Continuous Manufacturing and Advanced Crystallization, University of Strathclyde G1 1RD UK
| | - Matthias J Gutmann
- ISIS Pulsed Muon and Neutron Source, Rutherford Appleton Laboratory, Harwell Oxford Didcot OX11 0QX UK
| | - Harriott Nowell
- Diamond Light Source, Harwell Science and Innovation Campus Didcot OX11 0DE UK
| | - Ningjin Zhang
- Chemistry, Faculty of Natural and Environmental Sciences, Highfield Campus, University of Southampton Southampton SO17 1HE UK
| | - David R Allan
- Diamond Light Source, Harwell Science and Innovation Campus Didcot OX11 0DE UK
| |
Collapse
|
26
|
Macetti G, Genoni A. Initial Maximum Overlap Method for Large Systems by the Quantum Mechanics/Extremely Localized Molecular Orbital Embedding Technique. J Chem Theory Comput 2021; 17:4169-4182. [PMID: 34196174 DOI: 10.1021/acs.jctc.1c00388] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Quantum chemistry offers a large variety of methods to treat excited states. Many of them are based on a multireference wave function ansatz and are therefore characterized by an intrinsic complexity and high computational costs. To overcome these drawbacks and also some limitations of simpler single-reference approaches (e.g., configuration interaction singles and time-dependent density functional theory), the single-determinant Δself-consistent field-initial maximum overlap method (ΔSCF-IMOM) has been proposed. This strategy substitutes the aufbau principle with a criterion that occupies molecular orbitals at successive SCF iterations on the basis of their maximum overlap with a proper set of guess orbitals for the target excited state. In this way, it prevents the SCF process to collapse to the ground state wave function and provides excited state single Slater determinant solutions to the SCF equations. Here, we propose to extend the applicability of the IMOM to the treatment of localized excited states of large systems. To accomplish this task, we coupled it with the QM/ELMO (quantum mechanics/extremely localized molecular orbitals) strategy, a quantum mechanical embedding method in which the most chemically relevant part of the system is treated with traditional quantum chemical approaches, while the rest is described by extremely localized molecular orbitals transferred from recently constructed libraries or proper model molecules. After presenting the theoretical foundations of the new IMOM/ELMO technique, in this paper, we will show and discuss the results of preliminary test calculations carried out on both model systems (i.e., decanoic acid, decene, decapentaene, and solvated acrolein) and a system of biological interest (flavin mononucleotide in the flavodoxin protein). We observed that, for localized excited states, the new IMOM/ELMO method provides reliable results, and it reproduces the outcomes of fully IMOM calculations within the chemical accuracy threshold (i.e., 0.043 eV) by including only a limited number of atoms in the QM region. Furthermore, the first application of our embedding technique to a larger biological system gave completely plausible results in line with those obtained through more traditional quantum mechanical methods, thus opening the possibility of using the new approach in future investigations of photobiology problems.
Collapse
Affiliation(s)
- Giovanni Macetti
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019, 1 Boulevard Arago, F-57078 Metz, France
| | - Alessandro Genoni
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019, 1 Boulevard Arago, F-57078 Metz, France
| |
Collapse
|
27
|
Gruene T, Holstein JJ, Clever GH, Keppler B. Establishing electron diffraction in chemical crystallography. Nat Rev Chem 2021; 5:660-668. [PMID: 37118416 DOI: 10.1038/s41570-021-00302-4] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 06/07/2021] [Indexed: 02/06/2023]
Abstract
The emerging field of 3D electron diffraction (3D ED) opens new opportunities for structure determination from sub-micrometre-sized crystals. Although the foundations of this technology emerged earlier, the past decade has seen developments in cryo-electron microscopy and (X-ray) crystallography that particularly enable the widespread use of 3D ED. This Perspective describes to chemists and chemical crystallographers just how similar electron and X-ray diffraction are and discusses their complementary aspects. We wish to establish 3D ED in the broader chemistry community, such that electron crystallography becomes a common part of the analytical chemistry toolkit. With a suitable instrument at their disposal, every skilled crystallographer can quickly learn to perform structure determinations using 3D ED.
Collapse
|
28
|
Milašinović V, Molčanov K, Krawczuk A, Bogdanov NE, Zakharov BA, Boldyreva EV, Jelsch C, Kojić-Prodić B. Charge density studies of multicentre two-electron bonding of an anion radical at non-ambient temperature and pressure. IUCRJ 2021; 8:644-654. [PMID: 34258012 PMCID: PMC8256703 DOI: 10.1107/s2052252521005273] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/30/2020] [Accepted: 05/18/2021] [Indexed: 06/13/2023]
Abstract
The variation of charge density of two-electron multicentre bonding (pancake bonding) between semi-quinone radicals with pressure and temperature was studied on a salt of 5,6-di-chloro-2,3-di-cyano-semi-quinone radical anion (DDQ) with 4-cyano-N-methyl-pyridinium cation (4-CN) using the Transferable Aspheric Atom Model (TAAM) refinement. The pancake-bonded radical dimers are stacked by non-bonding π-interactions. With rising pressure, the covalent character of interactions between radicals increases, and above 2.55 GPa, the electron density indicates multicentric covalent interactions throughout the stack. The experimental charge densities were verified and corroborated by periodic DFT computations. The TAAM approach has been tested and validated for atomic resolution data measured at ambient pressure; this work shows this approach can also be applied to diffraction data obtained at pressures up to several gigapascals.
Collapse
Affiliation(s)
- Valentina Milašinović
- Department of Physical Chemistry, Rudjer Bošković Institute, Bijenička 54, Zagreb 10000, Croatia
| | - Krešimir Molčanov
- Department of Physical Chemistry, Rudjer Bošković Institute, Bijenička 54, Zagreb 10000, Croatia
| | - Anna Krawczuk
- Institute of Inorganic Chemistry, University of Göttingen, Tammannstrasse 4, Göttingen 37077, Germany
- Faculty of Chemistry, Jagiellonian University in Krakow, Gronostajowa 2, Krakow 30-387, Poland
| | - Nikita E. Bogdanov
- Boreskov Institute of Catalysis, SB RAS, Lavrentiev Avenue 5, Novosibirsk 630090, Russian Federation
- Novosibirsk State University, Pirogova Street 2, Novosibirsk 630090 Russian Federation
| | - Boris A. Zakharov
- Boreskov Institute of Catalysis, SB RAS, Lavrentiev Avenue 5, Novosibirsk 630090, Russian Federation
- Novosibirsk State University, Pirogova Street 2, Novosibirsk 630090 Russian Federation
| | - Elena V. Boldyreva
- Boreskov Institute of Catalysis, SB RAS, Lavrentiev Avenue 5, Novosibirsk 630090, Russian Federation
- Novosibirsk State University, Pirogova Street 2, Novosibirsk 630090 Russian Federation
| | - Christian Jelsch
- CRM2, CNRS, UMR 7036, Université de Lorraine, BP 70239 Nancy, France
| | - Biserka Kojić-Prodić
- Department of Physical Chemistry, Rudjer Bošković Institute, Bijenička 54, Zagreb 10000, Croatia
| |
Collapse
|
29
|
Macetti G, Genoni A. Three-Layer Multiscale Approach Based on Extremely Localized Molecular Orbitals to Investigate Enzyme Reactions. J Phys Chem A 2021; 125:6013-6027. [PMID: 34190569 DOI: 10.1021/acs.jpca.1c05040] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Quantum mechanics/molecular mechanics (QM/MM) calculations are widely used embedding techniques to computationally investigate enzyme reactions. In most QM/MM computations, the quantum mechanical region is treated through density functional theory (DFT), which offers the best compromise between chemical accuracy and computational cost. Nevertheless, to obtain more accurate results, one should resort to wave function-based methods, which however lead to a much larger computational cost already for relatively small QM subsystems. To overcome this drawback, we propose the coupling of our QM/ELMO (quantum mechanics/extremely localized molecular orbital) approach with molecular mechanics, thus introducing the three-layer QM/ELMO/MM technique. The QM/ELMO strategy is an embedding method in which the chemically relevant part of the system is treated at the quantum mechanical level, while the rest is described through frozen ELMOs. Since the QM/ELMO method reproduces results of fully QM computations within chemical accuracy and with a much lower computational effort, it can be considered a suitable strategy to extend the range of applicability and accuracy of the QM/MM scheme. In this paper, other than briefly presenting the theoretical bases of the QM/ELMO/MM technique, we will also discuss its validation on the well-tested deprotonation of acetyl coenzyme A by aspartate in citrate synthase.
Collapse
Affiliation(s)
- Giovanni Macetti
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019, 1 Boulevard Arago, F-57078 Metz, France
| | - Alessandro Genoni
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019, 1 Boulevard Arago, F-57078 Metz, France
| |
Collapse
|
30
|
Wanat M, Malinska M, Hoser AA, Woźniak K. Further Validation of Quantum Crystallography Approaches. Molecules 2021; 26:molecules26123730. [PMID: 34207308 PMCID: PMC8233966 DOI: 10.3390/molecules26123730] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/12/2021] [Revised: 06/14/2021] [Accepted: 06/14/2021] [Indexed: 11/16/2022] Open
Abstract
Quantum crystallography is a fast-developing multidisciplinary area of crystallography. In this work, we analyse the influence of different charge density models (i.e., the multipole model (MM), Hirshfeld atom refinement (HAR), and the transferable aspherical atom model (TAAM)), modelling of the thermal motion of hydrogen atoms (anisotropic, isotropic, and with the aid of SHADE or NoMoRe), and the type of radiation used (Mo Kα and Cu Kα) on the final results. To achieve this aim, we performed a series of refinements against X-ray diffraction data for three model compounds and compared their final structures, geometries, shapes of ADPs, and charge density distributions. Our results were also supported by theoretical calculations that enabled comparisons of the lattice energies of these structures. It appears that geometrical parameters are better described (closer to the neutron values) when HAR is used; however, bonds to H atoms more closely match neutron values after MM or TAAM refinement. Our analysis shows the superiority of the NoMoRe method in the description of H-atom ADPs. Moreover, the shapes of the ADPs of H atoms, as well as their electron density distributions, were better described with low-resolution Cu Kα data in comparison to low-resolution Mo Kα data.
Collapse
Affiliation(s)
- Monika Wanat
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, 101 Żwirki i Wigury, 02-089 Warszawa, Poland; (M.W.); (M.M.); (A.A.H.)
- College of Inter-Faculty Individual Studies in Mathematics and Natural Sciences (MISMaP), University of Warsaw, 2C Stefana Banacha, 02-097 Warszawa, Poland
| | - Maura Malinska
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, 101 Żwirki i Wigury, 02-089 Warszawa, Poland; (M.W.); (M.M.); (A.A.H.)
| | - Anna A. Hoser
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, 101 Żwirki i Wigury, 02-089 Warszawa, Poland; (M.W.); (M.M.); (A.A.H.)
| | - Krzysztof Woźniak
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, 101 Żwirki i Wigury, 02-089 Warszawa, Poland; (M.W.); (M.M.); (A.A.H.)
- Correspondence:
| |
Collapse
|
31
|
Malaspina LA, Genoni A, Grabowsky S. lamaGOET: an interface for quantum crystallography. J Appl Crystallogr 2021; 54:987-995. [PMID: 34188618 PMCID: PMC8202027 DOI: 10.1107/s1600576721002545] [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: 09/11/2020] [Accepted: 03/08/2021] [Indexed: 11/24/2022] Open
Abstract
In quantum crystallography, theoretical calculations and crystallographic refinements are closely intertwined. This means that the employed software must be able to perform both quantum-mechanical calculations and crystallographic least-squares refinements. So far, the program Tonto is the only one able to do that. The lamaGOET interface described herein deals with this issue since it interfaces dedicated quantum-chemical software (the widely used Gaussian package and the specialized ELMOdb program) with the refinement capabilities of Tonto. Three different flavours of quantum-crystallographic refinements of the dipetide glycyl-l-threonine dihydrate are presented to showcase the capabilities of lamaGOET: Hirshfeld atom refinement (HAR), HAR-ELMO, namely HAR coupled with extremely localized molecular orbitals, and X-ray constrained wavefunction fitting.
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
| | - 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
|
32
|
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
|
33
|
Peng D, Nakashima PNH. Measuring Density Functional Parameters from Electron Diffraction Patterns. PHYSICAL REVIEW LETTERS 2021; 126:176402. [PMID: 33988394 DOI: 10.1103/physrevlett.126.176402] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/17/2020] [Revised: 01/13/2021] [Accepted: 03/15/2021] [Indexed: 06/12/2023]
Abstract
We integrate density functional theory (DFT) into quantitative convergent-beam electron diffraction (QCBED) to create a synergy between experiment and theory called QCBED-DFT. This synergy resides entirely in the electron density which, in real materials, gives rise to the experimental CBED patterns used by QCBED-DFT to refine DFT model parameters. We use it to measure the Hubbard energy U for two strongly correlated electron systems, NiO and CeB_{6} (U=7.4±0.6 eV for d orbitals in NiO and U=3.0±0.6 eV for f orbitals in CeB_{6}), and the boron position parameter x for CeB_{6} (x=0.1992±0.0003). In verifying our measurements, we demonstrate an accuracy test for any modeled electron density.
Collapse
Affiliation(s)
- Ding Peng
- Department of Materials Science and Engineering, Monash University, Victoria 3800, Australia
| | - Philip N H Nakashima
- Department of Materials Science and Engineering, Monash University, Victoria 3800, Australia
| |
Collapse
|
34
|
Shteingolts SA, Stash AI, Tsirelson VG, Fayzullin RR. Orbital-Free Quantum Crystallographic View on Noncovalent Bonding: Insights into Hydrogen Bonds, π⋅⋅⋅π and Reverse Electron Lone Pair⋅⋅⋅π Interactions. Chemistry 2021; 27:7789-7809. [PMID: 33769620 DOI: 10.1002/chem.202005497] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/28/2020] [Indexed: 01/16/2023]
Abstract
A detailed analysis of a complete set of the local potentials that appear in the Euler equation for electron density is carried out for noncovalent interactions in the crystal of a uracil derivative using experimental X-ray charge density. The interplay between the quantum theory of atoms in molecules and crystals and the local potentials and corresponding inner-crystal electronic forces of electrostatic and kinetic origin is explored. Partitioning of crystal space into atomic basins and atomic-like potential basins led us to the definite description of interatomic interaction and charge transfer. Novel physically grounded bonding descriptors derived within the orbital-free quantum crystallography provided the detailed examination of π-stacking and intricate C=O⋅⋅⋅π interactions and nonclassical hydrogen bonds present in the crystal. The donor-acceptor character of these interactions is revealed by analysis of Pauli and von Weizsäcker potentials together with well-known functions, e. g., deformation electron density and electron localization function. In this way, our analysis throws light on aspects of these closed-shell interactions hitherto hidden from the description.
Collapse
Affiliation(s)
- Sergey A Shteingolts
- Arbuzov Institute of Organic and Physical Chemistry FRC Kazan Scientific Center, Russian Academy of Sciences, 8 Arbuzov Street, Kazan, 420088, Russian Federation
| | - Adam I Stash
- A.N. Nesmeyanov Institute of Organoelement Compounds Russian Academy of Sciences, 28 Vavilov Street, Moscow, 119991, Russian Federation
| | - Vladimir G Tsirelson
- D.I. Mendeleev University of Chemical Technology, 9 Miusskaya Square, Moscow, 125047, Russian Federation
| | - Robert R Fayzullin
- Arbuzov Institute of Organic and Physical Chemistry FRC Kazan Scientific Center, Russian Academy of Sciences, 8 Arbuzov Street, Kazan, 420088, Russian Federation
| |
Collapse
|
35
|
Bergmann J, Oksanen E, Ryde U. Quantum-refinement studies of the bidentate ligand of V‑nitrogenase and the protonation state of CO-inhibited Mo‑nitrogenase. J Inorg Biochem 2021; 219:111426. [PMID: 33756394 DOI: 10.1016/j.jinorgbio.2021.111426] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2020] [Revised: 02/19/2021] [Accepted: 03/10/2021] [Indexed: 10/21/2022]
Abstract
Nitrogenase is the only enzyme that can cleave the triple bond in N2, making nitrogen available to plants (although the enzyme itself is strictly microbial). It has been studied extensively with both experimental and computational methods, but many details of the reaction mechanism are still unclear. X-ray crystallography is the main source of structural information for biomacromolecules, but it has problems to discern hydrogen atoms or to distinguish between elements with the same number of electrons. These problems can sometimes be alleviated by introducing quantum chemical calculations in the refinement, providing information about the ideal structure (in the same way as the empirical restraints used in standard crystallographic refinement) and comparing different interpretations of the structure with normal crystallographic and quantum mechanical quality measures. We have performed such quantum-refinement calculations to address two important issues for nitrogenase. First, we show that the bidentate ligand of the active-site FeV cluster in V‑nitrogenase is carbonate, rather than bicarbonate or nitrate. Second, we study the CO-inhibited structure of Mo‑nitrogenase. CO binds to a reduced and protonated state of the enzyme by replacing one of the sulfide ions (S2B) in the active-site FeMo cluster. We examined if it is possible to deduce from the crystal structure the location of the protons. Our results indicates that the crystal structure is best modelled as fully deprotonated.
Collapse
Affiliation(s)
- Justin Bergmann
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Esko Oksanen
- European Spallation Source ESS ERIC, Lund, Sweden
| | - Ulf Ryde
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden.
| |
Collapse
|
36
|
Macetti G, Wieduwilt EK, Genoni A. QM/ELMO: A Multi-Purpose Fully Quantum Mechanical Embedding Scheme Based on Extremely Localized Molecular Orbitals. J Phys Chem A 2021; 125:2709-2726. [DOI: 10.1021/acs.jpca.0c11450] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Affiliation(s)
- Giovanni Macetti
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019, 1 Boulevard Arago, F-57078 Metz, France
| | - Erna K. Wieduwilt
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019, 1 Boulevard Arago, F-57078 Metz, France
| | - Alessandro Genoni
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019, 1 Boulevard Arago, F-57078 Metz, France
| |
Collapse
|
37
|
Shteingolts SA, Saifina AF, Saifina LF, Semenov VE, Fukin GK, Fayzullin RR. X-ray charge density study of the 6-methyluracil derivative in the crystal: Revealing, consequences, and multipole refinement of minor static disorder. J Mol Struct 2021. [DOI: 10.1016/j.molstruc.2020.129724] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
|
38
|
Dittrich B. On modelling disordered crystal structures through restraints from molecule-in-cluster computations, and distinguishing static and dynamic disorder. IUCRJ 2021; 8:305-318. [PMID: 33708406 PMCID: PMC7924241 DOI: 10.1107/s2052252521000531] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 10/02/2020] [Accepted: 01/14/2021] [Indexed: 06/12/2023]
Abstract
Distinguishing disorder into static and dynamic based on multi-temperature X-ray or neutron diffraction experiments is the current state of the art, but is only descriptive, not predictive. Here, several disordered structures are revisited from the Cambridge Crystallographic Data Center 'drug subset', the Cambridge Structural Database and own earlier work, where experimental intensities of Bragg diffraction data were available. Using the molecule-in-cluster approach, structures with distinguishable conformations were optimized separately, as extracted from available or generated disorder models of the respective disordered crystal structures. Re-combining these 'archetype structures' by restraining positional and constraining displacement parameters for conventional least-squares refinement, based on the optimized geometries, then often achieves a superior fit to the experimental diffraction data compared with relying on experimental information alone. It also simplifies and standardizes disorder refinement. Ten example structures were analysed. It is observed that energy differences between separate disorder conformations are usually within a small energy window of RT (T = crystallization temperature). Further computations classify disorder into static or dynamic, using single experiments performed at one single temperature, and this was achieved for propionamide.
Collapse
Affiliation(s)
- Birger Dittrich
- Novartis Campus, Novartis Pharma AG, Postfach, Basel, CH-4002, Switzerland
- Mathematisch-Naturwissenschaftliche Fakultät, Universität Zürich, Winterthurerstrasse 190, Zürich, CH-8057, Switzerland
| |
Collapse
|
39
|
Crystallographic and Computational Electron Density of d x2-y2 Orbitals of Azo-Schiff Base Metal Complexes Using Conventional Programs. Molecules 2021; 26:molecules26030551. [PMID: 33494463 PMCID: PMC7865803 DOI: 10.3390/molecules26030551] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/24/2020] [Revised: 01/12/2021] [Accepted: 01/20/2021] [Indexed: 11/23/2022] Open
Abstract
The crystal structures of two azobenzene derivative Schiff base metal complexes (new C44H40CuN6O2 of P-1 and known C44H38MnN6O7 of P21/c abbreviated as Cu and Mn, respectively) were (re-)determined experimentally using conventional X-ray analysis to obtain electron density using a PLATON program. Cu affords a four-coordinated square planar geometry, while Mn affords a hexa-coordinated distorted octahedral geometry whose apical sites are occupied by an acetate ion and water ligands, which are associated with hydrogen bonds. The π-π or CH-π and hydrogen bonding intermolecular interactions were found in both crystals, which were also analyzed using a Hirshfeld surface analysis program. To compare these results with experimental results, a density functional theory (DFT) calculation was also carried out based on the crystal structures to obtain calculated electron density using a conventional Gaussian program. These results revealed that the axial Mn-O coordination bonds of Mn were relatively weaker than the in-plane M-N or M-O coordination bonds.
Collapse
|
40
|
Wieduwilt EK, Macetti G, Genoni A. Climbing Jacob's Ladder of Structural Refinement: Introduction of a Localized Molecular Orbital-Based Embedding for Accurate X-ray Determinations of Hydrogen Atom Positions. J Phys Chem Lett 2021; 12:463-471. [PMID: 33369421 DOI: 10.1021/acs.jpclett.0c03421] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
The positions of hydrogen atoms in molecules are fundamental in many aspects of chemistry. Nevertheless, most molecular structures are obtained from refinements of X-ray data exploiting the independent atom model (IAM), which uses spherical atomic densities and provides bond lengths involving hydrogen atoms that are too short compared to the neutron reference values. To overcome the IAM shortcomings, the wave function-based Hirshfeld atom refinement (HAR) method has been recently proposed, emerging as a promising strategy able to give element-hydrogen bond distances in excellent agreement with the neutron ones in terms of accuracy and precision. In this Letter, we propose a significant improvement of HAR based on the idea of describing the crystal environment explicitly in the underlying wave function calculation through a quantum mechanical embedding strategy that exploits extremely localized molecular orbitals. Test-bed refinements on a crystal structure characterized by strong intermolecular interactions are also discussed.
Collapse
Affiliation(s)
- Erna K Wieduwilt
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019, 1 Boulevard Arago, F-57078 Metz, France
| | - Giovanni Macetti
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019, 1 Boulevard Arago, F-57078 Metz, France
| | - Alessandro Genoni
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019, 1 Boulevard Arago, F-57078 Metz, France
| |
Collapse
|
41
|
Wieduwilt EK, Boisson JC, Terraneo G, Hénon E, Genoni A. A Step toward the Quantification of Noncovalent Interactions in Large Biological Systems: The Independent Gradient Model-Extremely Localized Molecular Orbital Approach. J Chem Inf Model 2021; 61:795-809. [PMID: 33444021 DOI: 10.1021/acs.jcim.0c01188] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
The independent gradient model (IGM) is a recent electron density-based computational method that enables to detect and quantify covalent and noncovalent interactions. When applied to large systems, the original version of the technique still relies on promolecular electron densities given by the sum of spherically averaged atomic electron distributions, which leads to approximate evaluations of the inter- and intramolecular interactions occurring in systems of biological interest. To overcome this drawback and perform IGM analyses based on quantum mechanically rigorous electron densities also for macromolecular systems, we coupled the IGM approach with the recently constructed libraries of extremely localized molecular orbitals (ELMOs) that allow fast and reliable reconstructions of polypeptide and protein electron densities. The validation tests performed on small polypeptides and peptide dimers have shown that the novel IGM-ELMO strategy provides results that are systematically closer to the fully quantum mechanical ones and outperforms the IGM method based on the crude promolecular approximation, but still keeping a quite low computational cost. The results of the test calculations carried out on proteins have also confirmed the trends observed for the IGM analyses conducted on small systems. This makes us envisage the future application of the novel IGM-ELMO approach to unravel complicated noncovalent interaction networks (e.g., in protein-protein contacts) or to rationally design new drugs through molecular docking calculations and virtual high-throughput screenings.
Collapse
Affiliation(s)
- Erna K Wieduwilt
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques, UMR CNRS 7019, 1 Boulevard Arago, Metz F-57078, France
| | - Jean-Charles Boisson
- CReSTIC EA 3804, Université de Reims Champagne-Ardenne, Moulin de la Housse, Reims Cedex 02 BP39, F-51687, France
| | - Giancarlo Terraneo
- Laboratory of Supramolecular and Bio-Nanomaterials (SupraBioNanoLab), Department of Chemistry, Materials and Chemical Engineering "Giulio Natta", Politecnico di Milano, Via L. Mancinelli 7, Milan I-20131, Italy
| | - Eric Hénon
- Institut de Chimie Moléculaire de Reims UMR CNRS 7312, Université de Reims Champagne-Ardenne, Moulin de la Housse, Reims Cedex 02 BP39, F-51687, France
| | - Alessandro Genoni
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques, UMR CNRS 7019, 1 Boulevard Arago, Metz F-57078, France
| |
Collapse
|
42
|
Podhorský M, Bučinský L, Jayatilaka D, Grabowsky S. HgH 2 meets relativistic quantum crystallography. How to teach relativity to a non-relativistic wavefunction. ACTA CRYSTALLOGRAPHICA A-FOUNDATION AND ADVANCES 2021; 77:54-66. [PMID: 33399131 DOI: 10.1107/s2053273320014837] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/30/2020] [Accepted: 11/09/2020] [Indexed: 11/10/2022]
Abstract
The capability of X-ray constrained wavefunction (XCW) fitting to introduce relativistic effects into a non-relativistic wavefunction is tested. It is quantified how much of the reference relativistic effects can be absorbed in the non-relativistic XCW calculation when fitted against relativistic structure factors of a model HgH2 molecule. Scaling of the structure-factor sets to improve the agreement statistics is found to introduce a significant systematic error into the XCW fitting of relativistic effects.
Collapse
Affiliation(s)
- Michal Podhorský
- Institute of Physical Chemistry and Chemical Physics FCHPT, Slovak University of Technology, Radlinskeho 9, Bratislava SK-812 37, Slovakia
| | - Lukáš Bučinský
- Institute of Physical Chemistry and Chemical Physics FCHPT, Slovak University of Technology, Radlinskeho 9, Bratislava SK-812 37, Slovakia
| | - Dylan Jayatilaka
- School of Molecular Sciences, The University of Western Australia, 35 Stirling Hwy, Perth WA 6009, Australia
| | - Simon Grabowsky
- Department of Chemistry and Biochemistry, University of Bern, Freiestrasse 3, 3012 Bern, Switzerland
| |
Collapse
|
43
|
Quantum mechanics/extremely localized molecular orbital embedding technique: Theoretical foundations and further validation. ADVANCES IN QUANTUM CHEMISTRY 2021. [DOI: 10.1016/bs.aiq.2021.05.004] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
44
|
Affiliation(s)
- Piero Macchi
- Department, Chemistry, Materials and Chemical Engineering, Politecnico di Milano, Milano, Italy
- Center for Nano Science and Technology CNST@polimi, Italian Institute of Technology, Milano, Italy
| |
Collapse
|
45
|
Macetti G, Genoni A. Quantum Mechanics/Extremely Localized Molecular Orbital Embedding Strategy for Excited States: Coupling to Time-Dependent Density Functional Theory and Equation-of-Motion Coupled Cluster. J Chem Theory Comput 2020; 16:7490-7506. [PMID: 33241930 DOI: 10.1021/acs.jctc.0c00956] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The QM/ELMO (quantum mechanics/extremely localized molecular orbital) method is a recently developed embedding technique in which the most important region of the system under examination is treated at fully quantum mechanical level, while the rest is described by means of transferred and frozen extremely localized molecular orbitals. In this paper, we propose the first application of the QM/ELMO approach to the investigation of excited states, and, in particular, we present the coupling of the QM/ELMO philosophy with Time-Dependent Density Functional Theory (TDDFT) and Equation-of-Motion Coupled Cluster with single and double substitutions (EOM-CCSD). The proposed TDDFT/ELMO and EOM-CCSD/ELMO strategies underwent a series of preliminary tests that were already considered for the validation of other embedding methods for excited states. The obtained results showed that the novel techniques allow the accurate description of localized excitations in large systems by only including a relatively small number of atoms in the region treated at fully quantum chemical level. Furthermore, for TDDFT/ELMO, it was also observed that (i) the method enables to avoid the presence of artificial low-lying charge-transfer states that may affect traditional TDDFT calculations, even using functionals that do not take into account long-range corrections, and (ii) the novel approach can be also successfully exploited to investigate local electronic transitions in quite large systems (e.g., reduced model of the Green Fluorescent Protein), and the accuracy of the results can be improved by including a sufficient number of chemically crucial fragments/residues in the quantum mechanical region. Finally, concerning EOM-CCSD/ELMO, it was also seen that, despite the quite crude approximation of an embedding potential given by frozen extremely localized molecular orbitals, the new strategy is able to satisfactorily account for the effects of the environment. This work paves the way to further extensions of the QM/ELMO philosophy for the study of local excitations in extended systems, suggesting the coupling of the QM/ELMO approach with other quantum chemical strategies for excited states, from the simplest ΔSCF techniques to the most advanced and computationally expensive multireferences methods.
Collapse
Affiliation(s)
- Giovanni Macetti
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019, 1 Boulevard Arago, F-57078 Metz, France
| | - Alessandro Genoni
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019, 1 Boulevard Arago, F-57078 Metz, France
| |
Collapse
|
46
|
Kleemiss F, Dolomanov OV, Bodensteiner M, Peyerimhoff N, Midgley L, Bourhis LJ, Genoni A, Malaspina LA, Jayatilaka D, Spencer JL, White F, Grundkötter-Stock B, Steinhauer S, Lentz D, Puschmann H, Grabowsky S. Accurate crystal structures and chemical properties from NoSpherA2. Chem Sci 2020; 12:1675-1692. [PMID: 34163928 PMCID: PMC8179328 DOI: 10.1039/d0sc05526c] [Citation(s) in RCA: 114] [Impact Index Per Article: 28.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2020] [Accepted: 11/06/2020] [Indexed: 12/21/2022] Open
Abstract
The relationship between the structure and the properties of a drug or material is a key concept of chemistry. Knowledge of the three-dimensional structure is considered to be of such importance that almost every report of a new chemical compound is accompanied by an X-ray crystal structure - at least since the 1970s when diffraction equipment became widely available. Crystallographic software of that time was restricted to very limited computing power, and therefore drastic simplifications had to be made. It is these simplifications that make the determination of the correct structure, especially when it comes to hydrogen atoms, virtually impossible. We have devised a robust and fast system where modern chemical structure models replace the old assumptions, leading to correct structures from the model refinement against standard in-house diffraction data using no more than widely available software and desktop computing power. We call this system NoSpherA2 (Non-Spherical Atoms in Olex2). We explain the theoretical background of this technique and demonstrate the far-reaching effects that the improved structure quality that is now routinely available can have on the interpretation of chemical problems exemplified by five selected examples.
Collapse
Affiliation(s)
- Florian Kleemiss
- Universität Bern, Departement für Chemie und Biochemie Freiestrasse 3 3012 Bern Switzerland
| | | | - Michael Bodensteiner
- Universität Regensburg, Fakultät für Chemie und Pharmazie, Universitätsstr. 31 93053 Regensburg Germany
| | - Norbert Peyerimhoff
- Durham University, Department of Mathematical Sciences South Road Durham DH1 3LE UK
| | - Laura Midgley
- Durham University, Department of Mathematical Sciences South Road Durham DH1 3LE UK
| | - Luc J Bourhis
- Bruker France 4 Allée Lorentz, Champs-sur-Marne 77447 Marne-la-Vallée cedex 2 France
| | - Alessandro Genoni
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019 1 Boulevard Arago 57078 Metz France
| | - Lorraine A Malaspina
- Universität Bern, Departement für Chemie und Biochemie Freiestrasse 3 3012 Bern Switzerland
| | - Dylan Jayatilaka
- University of Western Australia, School of Molecular Sciences 35 Stirling Highway WA 6009 Perth Australia
| | - John L Spencer
- Victoria University of Wellington, School of Chemical and Physical Sciences Wellington 6012 New Zealand
| | - Fraser White
- Rigaku Europe SE Hugenottenallee 167 63263 Neu-Isenburg Germany
| | - Bernhard Grundkötter-Stock
- Freie Universität Berlin, Institut für Chemie und Biochemie Anorganische Chemie, Fabeckstr. 34/36 14195 Berlin Germany
| | - Simon Steinhauer
- Freie Universität Berlin, Institut für Chemie und Biochemie Anorganische Chemie, Fabeckstr. 34/36 14195 Berlin Germany
| | - Dieter Lentz
- Freie Universität Berlin, Institut für Chemie und Biochemie Anorganische Chemie, Fabeckstr. 34/36 14195 Berlin Germany
| | | | - Simon Grabowsky
- Universität Bern, Departement für Chemie und Biochemie Freiestrasse 3 3012 Bern Switzerland
| |
Collapse
|
47
|
Cao L, Ryde U. Quantum refinement with multiple conformations: application to the P-cluster in nitrogenase. Acta Crystallogr D Struct Biol 2020; 76:1145-1156. [PMID: 33135685 PMCID: PMC7604908 DOI: 10.1107/s2059798320012917] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2020] [Accepted: 09/21/2020] [Indexed: 11/17/2022] Open
Abstract
X-ray crystallography is the main source of atomistic information on the structure of proteins. Normal crystal structures are obtained as a compromise between the X-ray scattering data and a set of empirical restraints that ensure chemically reasonable bond lengths and angles. However, such restraints are not always available or accurate for nonstandard parts of the structure, for example substrates, inhibitors and metal sites. The method of quantum refinement, in which these empirical restraints are replaced by quantum-mechanical (QM) calculations, has previously been suggested for small but interesting parts of the protein. Here, this approach is extended to allow for multiple conformations in the QM region by performing separate QM calculations for each conformation. This approach is shown to work properly and leads to improved structures in terms of electron-density maps and real-space difference density Z-scores. It is also shown that the quality of the structures can be gauged using QM strain energies. The approach, called ComQumX-2QM, is applied to the P-cluster in two different crystal structures of the enzyme nitrogenase, i.e. an Fe8S7Cys6 cluster, used for electron transfer. One structure is at a very high resolution (1.0 Å) and shows a mixture of two different oxidation states, the fully reduced PN state (Fe82+, 20%) and the doubly oxidized P2+ state (80%). In the original crystal structure the coordinates differed for only two iron ions, but here it is shown that the two states also show differences in other atoms of up to 0.7 Å. The second structure is at a more modest resolution, 2.1 Å, and was originally suggested to show only the one-electron oxidized state, P1+. Here, it is shown that it is rather a 50/50% mixture of the P1+ and P2+ states and that many of the Fe-Fe and Fe-S distances in the original structure were quite inaccurate (by up to 0.8 Å). This shows that the new ComQumX-2QM approach can be used to sort out what is actually seen in crystal structures with dual conformations and to give locally improved coordinates.
Collapse
Affiliation(s)
- Lili Cao
- Department of Theoretical Chemistry, Lund University, PO Box 124, 221 00 Lund, Sweden
| | - Ulf Ryde
- Department of Theoretical Chemistry, Lund University, PO Box 124, 221 00 Lund, Sweden
| |
Collapse
|
48
|
Sadjadi S, Matta CF, Hamilton IP. Bonding and metastability for Group 12 dications. J Comput Chem 2020; 42:40-49. [PMID: 33063900 DOI: 10.1002/jcc.26431] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/06/2020] [Revised: 09/25/2020] [Accepted: 09/28/2020] [Indexed: 11/06/2022]
Abstract
Electronic structure and bonding properties of the Group 12 dications M2 2+ (M = Zn, Cd, Hg) are investigated and electron density-derived quantities are used to characterize the metastability of these species. Of particular interest are the complementary descriptions afforded by the Laplacian of the electron density ∇2 ρ(r) and the one-electron Bohm quantum potential (Q = ∇ 2 ρ r / 2 ρ r ) along the bond path. Further, properties derived from the pair density including the localization-delocalization matrices (LDMs) and the interacting quantum atoms (IQA) energies are analyzed within the framework of the quantum theory of atoms in molecules (QTAIM). From the crossing points of the singlet (ground) and triplet (excited) potential energy curves, the barriers for dissociation (BFD) are estimated to be 25.2 kcal/mol (1.09 eV) for Zn2 2+ , 22.8 kcal/mol (0.99 eV) for Cd2 2+ , and 26.4 kcal/mol (1.14 eV) for Hg2 2+ . For comparison and benchmarking purposes, the case of N2 2+ is considered as a texbook example of metastability. At the equilibrium geometries, LDMs, which are used here as an electronic fingerprinting tool, discriminate and group together Group 12 M2 2+ from its isoelectronic Group 11 M2 . While "classical" bonding indices are inconclusive in establishing regions of metastability in the bonding, it is shown that the one-electron Bohm quantum potential is promising in this regard.
Collapse
Affiliation(s)
- SeyedAbdolreza Sadjadi
- Department of Physics, Faculty of Science, Laboratory for Space Research, The University of Hong Kong, Hong Kong SAR, China
| | - Chérif F Matta
- Department of Chemistry and Physics, Mount Saint Vincent University, Halifax, Nova Scotia, Canada.,Department of Chemistry, Dalhousie University, Halifax, Nova Scotia, Canada
| | - Ian P Hamilton
- Department of Chemistry and Biochemistry, Wilfrid Laurier University, Waterloo, Ontario, Canada
| |
Collapse
|
49
|
Gatti C. Looking at local classical and quantum forces in stable crystals using multipole-model refined electron densities and orbital-free DFT approximations. ACTA CRYSTALLOGRAPHICA SECTION B, STRUCTURAL SCIENCE, CRYSTAL ENGINEERING AND MATERIALS 2020; 76:724-726. [PMID: 33017305 DOI: 10.1107/s2052520620012895] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Two distinct approaches, that of energy and that of force, are adopted in quantum mechanics to get insights on chemical processes. In the second one, the net forces acting on the electrons and nuclei in a system (Ehrnefest and Hellmann-Feynman forces, respectively) are determined and a local version of the approach, in terms of force density fields rather than forces, has also been proposed for electrons. This is the path followed by Tsirelson & Stash (2020) in this issue of Acta Crystallographica Section B, to study for the first time the spatial distribution of the electronic forces of different nature acting in stable crystals. Interestingly, by relying on approximations taken from orbital-free DFT, all components of the inner-crystal force can be easily retrieved from multipole-model refined experimental electron densities and their derivatives. No less important is that these calculations are becoming easily doable for any X-ray density crystallographer thanks to a new version of the computer program WinXPRO, purposely developed in the study which is discussed in this commentary.
Collapse
Affiliation(s)
- Carlo Gatti
- CNR-SCITEC, via Golgi section, via Golgi 19, Milano, Italy 20134, Italy
| |
Collapse
|
50
|
Ricardi N, Ernst M, Macchi P, Wesolowski TA. Embedding-theory-based simulations using experimental electron densities for the environment. Acta Crystallogr A Found Adv 2020; 76:571-579. [PMID: 32869754 PMCID: PMC7459768 DOI: 10.1107/s2053273320008062] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2020] [Accepted: 06/16/2020] [Indexed: 11/26/2022] Open
Abstract
The basic idea of frozen-density embedding theory (FDET) is the constrained minimization of the Hohenberg-Kohn density functional EHK[ρ] performed using the auxiliary functional E_{v_{AB}}^{\rm FDET}[\Psi _A, \rho _B], where ΨA is the embedded NA-electron wavefunction and ρB(r) is a non-negative function in real space integrating to a given number of electrons NB. This choice of independent variables in the total energy functional E_{v_{AB}}^{\rm FDET}[\Psi _A, \rho _B] makes it possible to treat the corresponding two components of the total density using different methods in multi-level simulations. The application of FDET using ρB(r) reconstructed from X-ray diffraction data for a molecular crystal is demonstrated for the first time. For eight hydrogen-bonded clusters involving a chromophore (represented as ΨA) and the glycylglycine molecule [represented as ρB(r)], FDET is used to derive excitation energies. It is shown that experimental densities are suitable for use as ρB(r) in FDET-based simulations.
Collapse
Affiliation(s)
- Niccolò Ricardi
- Department of Physical Chemistry, University of Geneva, 30, Quai Ernest-Ansermet, CH-1211 Genève 4, Switzerland
| | - Michelle Ernst
- University of Bern, Freiestraße 3, 3012 Bern, Switzerland
| | - Piero Macchi
- Department of Chemistry, Materials and Chemical Engineering, Polytechnic of Milan, via Mancinelli 7, Milano 20131, Italy
| | - Tomasz Adam Wesolowski
- Department of Physical Chemistry, University of Geneva, 30, Quai Ernest-Ansermet, CH-1211 Genève 4, Switzerland
| |
Collapse
|