1
|
Zwolenik A, Tchoń D, Makal A. Evolution of structure and spectroscopic properties of a new 1,3-diacetylpyrene polymorph with temperature and pressure. IUCRJ 2024; 11:519-527. [PMID: 38727170 PMCID: PMC11220879 DOI: 10.1107/s2052252524003634] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/04/2024] [Accepted: 04/22/2024] [Indexed: 07/04/2024]
Abstract
A new polymorph of 1,3-diacetylpyrene has been obtained from its melt and thoroughly characterized using single-crystal X-ray diffraction, steady-state UV-Vis spectroscopy and periodic density functional theory calculations. Experimental studies covered the temperature range from 90 to 390 K and the pressure range from atmospheric to 4.08 GPa. Optimal sample placement in a diamond anvil cell according to our previously presented methodology ensured over 80% data coverage up to 0.8 Å for a monoclinic sample. Unrestrained Hirshfeld atom refinement of the high-pressure crystal structures was successful and anharmonic behavior of carbonyl oxygen atoms was observed. Unlike the previously characterized polymorph, the structure of 2°AP-β is based on infinite π-stacks of antiparallel 2°AP molecules. 2°AP-β displays piezochromism and piezofluorochromism which are directly related to the variation in interplanar distances within the π-stacking. The importance of weak intermolecular interactions is reflected in the substantial negative thermal expansion coefficient of -55.8 (57) MK-1 in the direction of C-H...O interactions.
Collapse
Affiliation(s)
- A. Zwolenik
- Biological and Chemical Research Centre, Faculty of ChemistryUniversity of WarsawŻwirki i Wigury 10102-089WarszawaPoland
| | - D. Tchoń
- Biological and Chemical Research Centre, Faculty of ChemistryUniversity of WarsawŻwirki i Wigury 10102-089WarszawaPoland
- Molecular Biophysics and Integrated Bioimaging DivisionLawrence Berkeley National Laboratory,1 Cyclotron RoadBerkeleyCA94720USA
| | - A. Makal
- Biological and Chemical Research Centre, Faculty of ChemistryUniversity of WarsawŻwirki i Wigury 10102-089WarszawaPoland
| |
Collapse
|
2
|
Kumar A, Jha KK, Olech B, Goral T, Malinska M, Woźniak K, Dominiak PM. TAAM refinement on high-resolution experimental and simulated 3D ED/MicroED data for organic molecules. Acta Crystallogr C Struct Chem 2024; 80:264-277. [PMID: 38934273 PMCID: PMC11225613 DOI: 10.1107/s2053229624005357] [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/17/2023] [Accepted: 06/05/2024] [Indexed: 06/28/2024] Open
Abstract
3D electron diffraction (3D ED), or microcrystal electron diffraction (MicroED), has become an alternative technique for determining the high-resolution crystal structures of compounds from sub-micron-sized crystals. Here, we considered L-alanine, α-glycine and urea, which are known to form good-quality crystals, and collected high-resolution 3D ED data on our in-house TEM instrument. In this study, we present a comparison of independent atom model (IAM) and transferable aspherical atom model (TAAM) kinematical refinement against experimental and simulated data. TAAM refinement on both experimental and simulated data clearly improves the model fitting statistics (R factors and residual electrostatic potential) compared to IAM refinement. This shows that TAAM better represents the experimental electrostatic potential of organic crystals than IAM. Furthermore, we compared the geometrical parameters and atomic displacement parameters (ADPs) resulting from the experimental refinements with the simulated refinements, with the periodic density functional theory (DFT) calculations and with published X-ray and neutron crystal structures. The TAAM refinements on the 3D ED data did not improve the accuracy of the bond lengths between the non-H atoms. The experimental 3D ED data provided more accurate H-atom positions than the IAM refinements on the X-ray diffraction data. The IAM refinements against 3D ED data had a tendency to lead to slightly longer X-H bond lengths than TAAM, but the difference was statistically insignificant. Atomic displacement parameters were too large by tens of percent for L-alanine and α-glycine. Most probably, other unmodelled effects were causing this behaviour, such as radiation damage or dynamical scattering.
Collapse
Affiliation(s)
- Anil Kumar
- Biological and Chemical Research Centre Faculty of Chemistry University of Warsaw, ul Żwirki i Wigury 101 02-089 Warszawa Poland
| | - Kunal Kumar Jha
- Biological and Chemical Research Centre Faculty of Chemistry University of Warsaw, ul Żwirki i Wigury 101 02-089 Warszawa Poland
- Centre of New Technologies University of Warsaw, ul S Banacha 2c 02-097 Warszawa Poland
| | - Barbara Olech
- Biological and Chemical Research Centre Faculty of Chemistry University of Warsaw, ul Żwirki i Wigury 101 02-089 Warszawa Poland
- Centre of New Technologies University of Warsaw, ul S Banacha 2c 02-097 Warszawa Poland
| | - Tomasz Goral
- Biological and Chemical Research Centre Faculty of Chemistry University of Warsaw, ul Żwirki i Wigury 101 02-089 Warszawa Poland
- Centre of New Technologies University of Warsaw, ul S Banacha 2c 02-097 Warszawa Poland
| | - Maura Malinska
- Faculty of Chemistry University of Warsaw, Pasteura 1 02-093 Warszawa Poland
| | - Krzysztof Woźniak
- Centre of New Technologies University of Warsaw, ul S Banacha 2c 02-097 Warszawa Poland
- Faculty of Chemistry University of Warsaw, Pasteura 1 02-093 Warszawa Poland
| | - Paulina Maria Dominiak
- Biological and Chemical Research Centre Faculty of Chemistry University of Warsaw, ul Żwirki i Wigury 101 02-089 Warszawa Poland
| |
Collapse
|
3
|
Hawash M, Qneibi M, Natsheh H, Mohammed NH, Hamda LA, Kumar A, Olech B, Dominiak PM, Bdir S, Bdair M. Evaluating the Neuroprotective Potential of Novel Benzodioxole Derivatives in Parkinson's Disease via AMPA Receptor Modulation. ACS Chem Neurosci 2024; 15:2334-2349. [PMID: 38747411 DOI: 10.1021/acschemneuro.4c00163] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/06/2024] Open
Abstract
Parkinson's disease (PD) is a significant health issue because it gradually damages the nervous system. α-Amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) receptors play a significant role in the development of PD. The current investigation employed hybrid benzodioxole-propanamide (BDZ-P) compounds to get information on AMPA receptors, analyze their biochemical and biophysical properties, and assess their neuroprotective effects. Examining the biophysical characteristics of all the subunits of the AMPA receptor offers insights into the impact of BDZ-P on the desensitization and deactivation rate. It demonstrates a partial improvement in the locomotor capacities in a mouse model of Parkinson's disease. In addition, the in vivo experiment assessed the locomotor activity by utilizing the open-field test. Our findings demonstrated that BDZ-P7 stands out with its remarkable potency, inhibiting the GluA2 subunit nearly 8-fold with an IC50 of 3.03 μM, GluA1/2 by 7.5-fold with an IC50 of 3.14 μM, GluA2/3 by nearly 7-fold with an IC50 of 3.19 μM, and GluA1 by 6.5-fold with an IC50 of 3.2 μM, significantly impacting the desensitization and deactivation rate of the AMPA receptor. BDZ-P7 showed an in vivo impact of partially reinstating locomotor abilities in a mouse model of PD. The results above suggest that the BDZ-P7 compounds show great promise as top contenders for the development of novel neuroprotective therapies.
Collapse
Affiliation(s)
- Mohammed Hawash
- Department of Pharmacy, Faculty of Medicine and Health Sciences, An-Najah National University, P403, Nablus 00970, Palestine
| | - Mohammad Qneibi
- Department of Biomedical Sciences, Faculty of Medicine and Health Sciences, An-Najah National University, P403, Nablus 00970, Palestine
| | - Hiba Natsheh
- Department of Pharmacy, Faculty of Medicine and Health Sciences, An-Najah National University, P403, Nablus 00970, Palestine
| | - Noor Haj Mohammed
- Department of Biomedical Sciences, Faculty of Medicine and Health Sciences, An-Najah National University, P403, Nablus 00970, Palestine
| | - Lubaba Abu Hamda
- Department of Biomedical Sciences, Faculty of Medicine and Health Sciences, An-Najah National University, P403, Nablus 00970, Palestine
| | - Anil Kumar
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, 02-089 Warsaw, Poland
| | - Barbara Olech
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, 02-089 Warsaw, Poland
- Centre of New Technologies, University of Warsaw, ul. S. Banacha 2c, 02-097 Warsaw, Poland
| | - Paulina Maria Dominiak
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, 02-089 Warsaw, Poland
| | - Sosana Bdir
- Department of Biomedical Sciences, Faculty of Medicine and Health Sciences, An-Najah National University, P403, Nablus 00970, Palestine
| | - Mohammad Bdair
- Department of Biomedical Sciences, Faculty of Medicine and Health Sciences, An-Najah National University, P403, Nablus 00970, Palestine
| |
Collapse
|
4
|
Chodkiewicz M, Patrikeev L, Pawlędzio S, Woźniak K. Transferable Hirshfeld atom model for rapid evaluation of aspherical atomic form factors. IUCRJ 2024; 11:249-259. [PMID: 38446457 PMCID: PMC10916294 DOI: 10.1107/s2052252524001507] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/28/2023] [Accepted: 02/14/2024] [Indexed: 03/07/2024]
Abstract
Form factors based on aspherical models of atomic electron density have brought great improvement in the accuracies of hydrogen atom parameters derived from X-ray crystal structure refinement. Today, two main groups of such models are available, the banks of transferable atomic densities parametrized using the Hansen-Coppens multipole model which allows for rapid evaluation of atomic form factors and Hirshfeld atom refinement (HAR)-related methods which are usually more accurate but also slower. In this work, a model that combines the ideas utilized in the two approaches is tested. It uses atomic electron densities based on Hirshfeld partitions of electron densities, which are precalculated and stored in a databank. This model was also applied during the refinement of the structures of five small molecules. A comparison of the resulting hydrogen atom parameters with those derived from neutron diffraction data indicates that they are more accurate than those obtained with the Hansen-Coppens based databank, and only slightly less accurate than those obtained with a version of HAR that neglects the crystal environment. The advantage of using HAR becomes more noticeable when the effects of the environment are included. To speed up calculations, atomic densities were represented by multipole expansion with spherical harmonics up to l = 7, which used numerical radial functions (a different approach to that applied in the Hansen-Coppens model). Calculations of atomic form factors for the small protein crambin (at 0.73 Å resolution) took only 68 s using 12 CPU cores.
Collapse
Affiliation(s)
- Michał Chodkiewicz
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury 101, Warszawa 02-089, Poland
| | - Leonid Patrikeev
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury 101, Warszawa 02-089, Poland
| | - Sylwia Pawlędzio
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury 101, Warszawa 02-089, Poland
- Neutron Scattering Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831, USA
| | - Krzysztof Woźniak
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury 101, Warszawa 02-089, Poland
| |
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
|
Woińska M, Hoser AA, Chodkiewicz ML, Woźniak K. Enhancing hydrogen positions in X-ray structures of transition metal hydride complexes with dynamic quantum crystallography. IUCRJ 2024; 11:45-56. [PMID: 37990870 PMCID: PMC10833390 DOI: 10.1107/s205225252300951x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/25/2023] [Accepted: 10/31/2023] [Indexed: 11/23/2023]
Abstract
Hirshfeld atom refinement (HAR) is a method which enables the user to obtain more accurate positions of hydrogen atoms bonded to light chemical elements using X-ray data. When data quality permits, this method can be extended to hydrogen-bonded transition metals (TMs), as in hydride complexes. However, addressing hydrogen thermal motions with HAR, particularly in TM hydrides, presents a challenge. At the same time, proper description of thermal vibrations can be vital for determining hydrogen positions correctly. In this study, we employ tools such as SHADE3 and Normal Mode Refinement (NoMoRe) to estimate anisotropic displacement parameters (ADPs) for hydrogen atoms during HAR and IAM refinements performed for seven structures of TM (Fe, Ni, Cr, Nb, Rh and Os) and metalloid (Sb) hydride complexes for which both the neutron and the X-ray structures have been determined. A direct comparison between neutron and HAR/SHADE3/NoMoRe ADPs reveals that the similarity between neutron hydrogen ADPs and those estimated with NoMoRe or SHADE3 is significantly higher than when hydrogen ADPs are refined with HAR. Regarding TM-H bond lengths, traditional HAR exhibits a slight advantage over the other methods. However, combining NoMoRe/SHADE3 with HAR results in a minor decrease in agreement with neutron TM-H bond lengths. For the Cr complex, for which high-resolution X-ray data were collected, an investigation of resolution-related effects was possible.
Collapse
Affiliation(s)
- Magdalena Woińska
- Biological and Chemical Research Centre, Chemistry Department, University of Warsaw, Żwirki i Wigury 101, Warsaw 02-089, Poland
| | - Anna A. Hoser
- Biological and Chemical Research Centre, Chemistry Department, University of Warsaw, Żwirki i Wigury 101, Warsaw 02-089, Poland
| | - Michał L. Chodkiewicz
- Biological and Chemical Research Centre, Chemistry Department, University of Warsaw, Żwirki i Wigury 101, Warsaw 02-089, Poland
| | - Krzysztof Woźniak
- Biological and Chemical Research Centre, Chemistry Department, University of Warsaw, Żwirki i Wigury 101, Warsaw 02-089, Poland
| |
Collapse
|
7
|
Hawash M, Al-Smadi D, Kumar A, Olech B, Dominiak PM, Jaradat N, Antari S, Mohammed S, Nasasrh A, Abualhasan M, Musa A, Suboh S, Çapan İ, Qneibi M, Natsheh H. Characterization and Investigation of Novel Benzodioxol Derivatives as Antidiabetic Agents: An In Vitro and In Vivo Study in an Animal Model. Biomolecules 2023; 13:1486. [PMID: 37892167 PMCID: PMC10604990 DOI: 10.3390/biom13101486] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2023] [Revised: 09/30/2023] [Accepted: 10/03/2023] [Indexed: 10/29/2023] Open
Abstract
In this study, we synthesized benzodioxol carboxamide derivatives and investigated their antidiabetic potential. The synthesized compounds (Ia-Ic and IIa-IId) underwent characterization via HRMS, 1H-, 13CAPT-NMR, and MicroED. Their efficacy against α-amylase was assessed in vitro, while MTS assays were employed to gauge cytotoxicity across cancer and normal cell lines. Additionally, the antidiabetic impact of compound IIc was evaluated in vivo using a streptozotocin-induced diabetic mice model. Notably, IIa and IIc displayed potent α-amylase inhibition (IC50 values of 0.85 and 0.68 µM, respectively) while exhibiting a negligible effect on the Hek293t normal cell line (IC50 > 150 µM), suggesting their safety. Compound IId demonstrated significant activity against four cancer cell lines (26-65 µM). In vivo experiments revealed that five doses of IIc substantially reduced mice blood glucose levels from 252.2 mg/dL to 173.8 mg/dL in contrast to the control group. The compelling in vitro anticancer efficacy of IIc and its safety for normal cells underscores the need for further in vivo assessment of this promising compound. This research highlights the potential of benzodioxol derivatives as candidates for the future development of synthetic antidiabetic drugs.
Collapse
Affiliation(s)
- Mohammed Hawash
- Department of Pharmacy, Faculty of Medicine and Health Sciences, An-Najah National University, Nablus 00970, Palestine; (N.J.); (S.A.); (S.M.); (A.N.); (M.A.)
| | - Derar Al-Smadi
- Department of Chemistry, Faculty of Sciences, An-Najah National University, Nablus 00970, Palestine;
| | - Anil Kumar
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, 02-089 Warsaw, Poland; (A.K.); (B.O.); (P.M.D.)
| | - Barbara Olech
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, 02-089 Warsaw, Poland; (A.K.); (B.O.); (P.M.D.)
- Centre of New Technologies, University of Warsaw, ul. S. Banacha 2c, 02-097 Warsaw, Poland
| | - Paulina Maria Dominiak
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, 02-089 Warsaw, Poland; (A.K.); (B.O.); (P.M.D.)
| | - Nidal Jaradat
- Department of Pharmacy, Faculty of Medicine and Health Sciences, An-Najah National University, Nablus 00970, Palestine; (N.J.); (S.A.); (S.M.); (A.N.); (M.A.)
| | - Sarah Antari
- Department of Pharmacy, Faculty of Medicine and Health Sciences, An-Najah National University, Nablus 00970, Palestine; (N.J.); (S.A.); (S.M.); (A.N.); (M.A.)
| | - Sarah Mohammed
- Department of Pharmacy, Faculty of Medicine and Health Sciences, An-Najah National University, Nablus 00970, Palestine; (N.J.); (S.A.); (S.M.); (A.N.); (M.A.)
| | - Ala’a Nasasrh
- Department of Pharmacy, Faculty of Medicine and Health Sciences, An-Najah National University, Nablus 00970, Palestine; (N.J.); (S.A.); (S.M.); (A.N.); (M.A.)
| | - Murad Abualhasan
- Department of Pharmacy, Faculty of Medicine and Health Sciences, An-Najah National University, Nablus 00970, Palestine; (N.J.); (S.A.); (S.M.); (A.N.); (M.A.)
| | - Ahmed Musa
- Department of Biomedical Sciences, Faculty of Medicine and Health Sciences, An-Najah National University, Nablus 00970, Palestine; (A.M.); (S.S.); (M.Q.)
| | - Shorooq Suboh
- Department of Biomedical Sciences, Faculty of Medicine and Health Sciences, An-Najah National University, Nablus 00970, Palestine; (A.M.); (S.S.); (M.Q.)
| | - İrfan Çapan
- Department of Material and Material Processing Technologies, Technical Sciences Vocational College, Gazi University, 06560 Ankara, Turkey;
- Basic and Engineering Sciences Central Laboratory Application and Research Center (GUTMAM), Gazi University, 06500 Ankara, Turkey
| | - Mohammad Qneibi
- Department of Biomedical Sciences, Faculty of Medicine and Health Sciences, An-Najah National University, Nablus 00970, Palestine; (A.M.); (S.S.); (M.Q.)
| | - Hiba Natsheh
- Department of Pharmacy, Faculty of Medicine and Health Sciences, An-Najah National University, Nablus 00970, Palestine; (N.J.); (S.A.); (S.M.); (A.N.); (M.A.)
| |
Collapse
|
8
|
Woińska M, Pawlędzio S, Chodkiewicz ML, Woźniak K. Hirshfeld Atom Refinement of Metal-Organic Complexes: Treatment of Hydrogen Atoms Bonded to Transition Metals. J Phys Chem A 2023; 127:3020-3035. [PMID: 36947670 PMCID: PMC10084459 DOI: 10.1021/acs.jpca.2c06998] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/24/2023]
Abstract
Hydrogen positions in hydrides play a key role in hydrogen storage materials and high-temperature superconductors. Our recently published study of five crystal structures of transition-metal-bound hydride complexes showed that using aspherical atomic scattering factors for Hirshfeld atom refinement (HAR) resulted in a systematic elongation of metal-hydrogen bonds compared to using spherical scattering factors with the Independent Atom Model (IAM). Even though only standard-resolution X-ray data was used, for the highest-quality data, we obtained excellent agreement between the X-ray and the neutron-derived bond lengths. We present an extended version of this study including 10 crystal structures of metal-organic complexes containing hydrogen atoms bonded to transition-metal atoms for which both X-ray and neutron data are available. The neutron structures were used as a benchmark, and the X-ray structures were refined by applying Hirshfeld atom refinement using various basis sets and DFT functionals in order to investigate the influence of the technical aspects on the length of metal-hydrogen bonds. The result of including relativistic effects in the Hamiltonian and using a cluster of multipoles simulating interactions with a crystal environment during wave function calculations was examined. The effect of the data quality on the final result was also evaluated. The study confirms that a high quality of experimental data is the key factor allowing us to obtain significant improvement in transition metal (TM)-hydrogen bond lengths from HAR in comparison with the IAM. Individual adjustments and better choices of the basis set can improve hydrogen positions. Average differences between TM-H bond lengths obtained with various DFT functionals upon including relativistic effects or between double-ζ and triple-ζ basis sets were not statistically significant. However, if all bonds formed by H atoms were considered, significant differences caused by different refinement strategies were observed. Finally, we examined the refinement of atomic thermal motions. Anisotropic refinement of hydrogen thermal motions with HAR was feasible only in some cases, and isotropically refined hydrogen thermal motions were in similar agreement with neutron values whether obtained with HAR or with the IAM.
Collapse
Affiliation(s)
- Magdalena Woińska
- Biological and Chemical Research Centre, Chemistry Department, University of Warsaw, Żwirki i Wigury 101, 02-089 Warszawa, Poland
| | - Sylwia Pawlędzio
- Biological and Chemical Research Centre, Chemistry Department, University of Warsaw, Żwirki i Wigury 101, 02-089 Warszawa, Poland
| | - Michał L Chodkiewicz
- Biological and Chemical Research Centre, Chemistry Department, University of Warsaw, Żwirki i Wigury 101, 02-089 Warszawa, Poland
| | - Krzysztof Woźniak
- Biological and Chemical Research Centre, Chemistry Department, University of Warsaw, Żwirki i Wigury 101, 02-089 Warszawa, Poland
| |
Collapse
|
9
|
Jha KK, Kleemiss F, Chodkiewicz ML, Dominiak PM. Aspherical atom refinements on X-ray data of diverse structures including disordered and covalent organic framework systems: a time-accuracy trade-off. J Appl Crystallogr 2023; 56:116-127. [PMID: 36777135 PMCID: PMC9901929 DOI: 10.1107/s1600576722010883] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2022] [Accepted: 11/13/2022] [Indexed: 12/24/2022] Open
Abstract
Aspherical atom refinement is the key to achieving accurate structure models, displacement parameters, hydrogen-bond lengths and analysis of weak interactions, amongst other examples. There are various quantum crystallographic methods to perform aspherical atom refinement, including Hirshfeld atom refinement (HAR) and transferable aspherical atom model (TAAM) refinement. Both HAR and TAAM have their limitations and advantages, the former being more accurate and the latter being faster. With the advent of non-spherical atoms in Olex2 (NoSpherA2), it is now possible to overcome some limitations, like treating disorder, twinning and network structures, in aspherical refinements using HAR, TAAM or both together. TAAM refinement in NoSpherA2 showed significant improvement in refinement statistics compared with independent atom model (IAM) refinements on a diverse set of X-ray diffraction data. The sensitivity of TAAM towards poor data quality and disorder was observed in terms of higher refinement statistics for such structures. A comparison of IAM with TAAM and HAR in NoSpherA2 indicated that the time taken by TAAM refinements was of the same order of magnitude as that taken by IAM, while in HAR the time taken using a minimal basis set was 50 times higher than for IAM and rapidly increased with increasing size of the basis sets used. The displacement parameters for hydrogen and non-hydrogen atoms were very similar in both HAR and TAAM refinements. The hydrogen-bond lengths were slightly closer to neutron reference values in the case of HAR with higher basis sets than in TAAM. To benefit from the advantages of each method, a new hybrid refinement approach has been introduced, allowing a combination of IAM, HAR and TAAM in one structure refinement. Refinement of coordination complexes involving metal-organic compounds and network structures such as covalent organic frameworks and metal-organic frameworks is now possible in a hybrid mode such as IAM-TAAM or HAR-TAAM, where the metal atoms are treated via either the IAM or HAR method and the organic part via TAAM, thus reducing the computational costs without compromising the accuracy. Formal charges on the metal and ligand can also be introduced in hybrid-mode refinement.
Collapse
Affiliation(s)
- Kunal Kumar Jha
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, Warsaw, 02-089, Poland
| | - Florian Kleemiss
- Fakultät für Chemie und Pharmazie, Universität Regensburg, Universitätstrasse 31, Regensburg, Bayern 93053, Germany
| | - Michał Leszek Chodkiewicz
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, Warsaw, 02-089, Poland
| | - Paulina Maria Dominiak
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, Warsaw, 02-089, Poland
| |
Collapse
|
10
|
Landeros-Rivera B, Ramírez-Palma D, Cortés-Guzmán F, Dominiak PM, Contreras-García J. How do density functionals affect the Hirshfeld atom refinement? Phys Chem Chem Phys 2023; 25:12702-12711. [PMID: 36644944 DOI: 10.1039/d2cp04098k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
In this work, the effect of mixing different amounts of Hartree-Fock (HF) exchange with hybrid density functionals applied to the Hirshfeld atom refinement (HAR) of urea and oxalic acid dihydrate is explored. Together, the influence of using different basis sets, methods (including MP2 and HF) and cluster sizes (to model bulk effects) is studied. The results show that changing the amount of HF exchange, no matter the level of theory, has an impact almost exclusively on the H atom refinement parameters. Contrary to pure quantum mechanical calculations where good geometries are obtained with intermediate HF exchange mixtures, in the HAR the best match with neutron diffraction reference values is not necessarily found for these admixtures. While the non-hydrogen covalent bond lengths are insensitive to the combination of method or basis set employed, the X-H bond lengths always increase proportionally to the HF exchange for the analysed systems. This outcome is opposite to what is normally observed from geometry optimisations, i.e., shorter bonds are obtained with greater HF exchange. Additionally, the thermal ellipsoids tend to shrink with larger HF exchange, especially for the H atoms involved in strong hydrogen bonding. Thus, it may be the case that the development of density functionals or basis sets suitable for quantum crystallography should take a different path than those fitted for quantum chemistry calculations.
Collapse
Affiliation(s)
| | - David Ramírez-Palma
- Departamento de Química, Centro de Investigación y de Estudios Avanzados del IPN (CINVESTAV), Avenida IPN 2508, Col. San Pedro Zacatenco, Ciudad de México, 07360, Mexico
| | - Fernando Cortés-Guzmán
- Universidad Nacional Autónoma de México, Instituto de Química, Ciudad Universitaria, Ciudad de México, 04510, Mexico
| | - Paulina M Dominiak
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Poland
| | | |
Collapse
|
11
|
Chodkiewicz ML, Gajda R, Lavina B, Tkachev S, Prakapenka VB, Dera P, Wozniak K. Accurate crystal structure of ice VI from X-ray diffraction with Hirshfeld atom refinement. IUCRJ 2022; 9:573-579. [PMID: 36071798 PMCID: PMC9438488 DOI: 10.1107/s2052252522006662] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 05/03/2022] [Accepted: 06/28/2022] [Indexed: 06/15/2023]
Abstract
Water is an essential chemical compound for living organisms, and twenty of its different crystal solid forms (ices) are known. Still, there are many fundamental problems with these structures such as establishing the correct positions and thermal motions of hydrogen atoms. The list of ice structures is not yet complete as DFT calculations have suggested the existence of additional and - to date - unknown phases. In many ice structures, neither neutron diffraction nor DFT calculations nor X-ray diffraction methods can easily solve the problem of hydrogen atom disorder or accurately determine their anisotropic displacement parameters (ADPs). Here, accurate crystal structures of H2O, D2O and mixed (50%H2O/50%D2O) ice VI obtained by Hirshfeld atom refinement (HAR) of high-pressure single-crystal synchrotron and laboratory X-ray diffraction data are presented. It was possible to obtain O-H/D bond lengths and ADPs for disordered hydrogen atoms which are in good agreement with the corresponding single-crystal neutron diffraction data. These results show that HAR combined with X-ray diffraction can compete with neutron diffraction in detailed studies of polymorphic forms of ice and crystals of other hydrogen-rich compounds. As neutron diffraction is relatively expensive, requires larger crystals which can be difficult to obtain and access to neutron facilities is restricted, cheaper and more accessible X-ray measurements combined with HAR can facilitate the verification of the existing ice polymorphs and the quest for new ones.
Collapse
Affiliation(s)
- Michal L. Chodkiewicz
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury, Warszawa 02-089, Poland
| | - Roman Gajda
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury, Warszawa 02-089, Poland
| | - Barbara Lavina
- Advanced Photon Source, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, IL 60439, USA
| | - Sergey Tkachev
- Advanced Photon Source, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, IL 60439, USA
| | - Vitali B. Prakapenka
- Hawai’i Institute of Geophysics and Planetology, Université d’hawaï à mānoa, 1680 East-West Road, Honolulu, HI 96822, USA
| | - Przemyslaw Dera
- Hawai’i Institute of Geophysics and Planetology, Université d’hawaï à mānoa, 1680 East-West Road, Honolulu, HI 96822, USA
| | - Krzysztof Wozniak
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury, Warszawa 02-089, Poland
| |
Collapse
|
12
|
Jha KK, Gruza B, Sypko A, Kumar P, Chodkiewicz ML, Dominiak PM. Multipolar Atom Types from Theory and Statistical Clustering (MATTS) Data Bank: Restructurization and Extension of UBDB. J Chem Inf Model 2022; 62:3752-3765. [PMID: 35943747 PMCID: PMC9400107 DOI: 10.1021/acs.jcim.2c00144] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
A fast and accurate operational model of electron density
is crucial
in many scientific disciplines including crystallography, molecular
biology, pharmaceutical, and structural chemistry. In quantum crystallography,
the aspherical refinement of crystal structures is becoming increasingly
popular because of its accurate description in terms of physically
meaningful properties. The transferable aspherical atom model (TAAM)
is quick and precise, though it requires a robust algorithm for atom
typing and coverage of the most popular atom types present in small
organic molecules. Thus, the University at Buffalo Databank (UBDB)
has been renamed to the Multipolar Atom Types from Theory and Statistical
clustering (MATTS) data bank, broadened, restructured, and implemented
into the software DiSCaMB with 651 atom types obtained from 2316 small-molecule
crystal structures containing C, H, N, O, P, S, F, Cl, and Br atoms.
MATTS2021 data bank now covers most of the small molecules, peptides,
RNA, DNA, and some frequently occurring cations and anions in biological,
pharmaceutical, and organic materials, including the majority of known
crystal structures composed of the above elements. The multipole model
parameters (Pval, κ, κ′, Plm) obtained for different
atom types were greatly influenced by neighboring atom types, hybridization,
geometrical strain in the ring system, and charges on the molecule.
Contrary to previous findings, the atoms showing variable oxidation
states and ions deviate from the linear dependence of monopole-derived
charges on the expansion–contraction κ parameter.
Collapse
Affiliation(s)
- Kunal Kumar Jha
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, 02-089 Warszawa, Poland
| | - Barbara Gruza
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, 02-089 Warszawa, Poland
| | - Aleksandra Sypko
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, 02-089 Warszawa, Poland
| | - Prashant Kumar
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, 02-089 Warszawa, Poland
| | - Michał Leszek Chodkiewicz
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, 02-089 Warszawa, Poland
| | - Paulina Maria Dominiak
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, 02-089 Warszawa, Poland
| |
Collapse
|
13
|
Kulik M, Chodkiewicz ML, Dominiak PM. Theoretical 3D electron diffraction electrostatic potential maps of proteins modeled with a multipolar pseudoatom data bank. Acta Crystallogr D Struct Biol 2022; 78:1010-1020. [PMID: 35916225 PMCID: PMC9344478 DOI: 10.1107/s2059798322005836] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2021] [Accepted: 05/31/2022] [Indexed: 11/10/2022] Open
Abstract
The availability of atomic resolution experimental maps of electrostatic potential from 3D electron diffraction (3D ED) extends the possibility of investigating the electrostatic potential beyond the determination of non-H-atom positions. However, accurate tools to calculate this potential for macromolecules, without the use of expensive quantum calculations, are lacking. The University at Buffalo Data Bank (UBDB) gathers atom types that can be used to calculate accurate electrostatic potential maps via structure-factor calculations. Here, the transferable aspherical atom model (TAAM) is applied with UBDB to investigate theoretically obtained electrostatic potential maps of lysozyme and proteinase K, and compare them with experimental maps from 3D ED. UBDB better reproduces the molecular electrostatic potential of molecules within their entire volume compared with the neutral spherical models used in the popular independent atom model (IAM). Additionally, the theoretical electron-density maps of the studied proteins are shown and compared with the electrostatic potential maps. The atomic displacement parameters (B factors) may affect the electrostatic potential maps in a different way than in the case of electron-density maps. The computational method presented in this study could potentially facilitate the interpretation of the less resolved regions of cryo-electron microscopy density maps and pave the way for distinguishing between different ions/water molecules in the active sites of macromolecules in high-resolution structures, which is of interest for drug-design purposes.
Collapse
Affiliation(s)
- Marta Kulik
- Biological and Chemical Research Centre, Faculty of Chemistry, University of Warsaw, Zwirki i Wigury 101, 02-089 Warsaw, Poland
| | - Michał Leszek Chodkiewicz
- Biological and Chemical Research Centre, Faculty of Chemistry, University of Warsaw, Zwirki i Wigury 101, 02-089 Warsaw, Poland
| | - Paulina Maria Dominiak
- Biological and Chemical Research Centre, Faculty of Chemistry, University of Warsaw, Zwirki i Wigury 101, 02-089 Warsaw, Poland
| |
Collapse
|
14
|
Pawlędzio S, Malinska M, Kleemiss F, Grabowsky S, Woźniak K. Influence of modelling disorder on Hirshfeld atom refinement results of an organo-gold(I) compound. IUCRJ 2022; 9:497-507. [PMID: 35844484 PMCID: PMC9252150 DOI: 10.1107/s2052252522005309] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 03/21/2022] [Accepted: 05/19/2022] [Indexed: 06/15/2023]
Abstract
Details of the validation of disorder modelling with Hirshfeld atom refinement (HAR) for a previously investigated organo-gold(I) compound are presented here. The impact of refining disorder on HAR results is discussed using an analysis of the differences of dynamic structure factors. These dynamic structure factor differences are calculated from thermally smeared quantum mechanical electron densities based on wavefunctions that include or exclude electron correlation and relativistic effects. When disorder is modelled, the electron densities stem from a weighted superposition of two (or more) different conformers. Here this is shown to impact the relative importance of electron correlation and relativistic effect estimates expressed by the structure factor magnitudes. The role of disorder modelling is also compared with the effect of the treatment of hydrogen anisotropic displacement parameter (ADP) values and atomic anharmonicity of the gold atom. The analysis of ADP values of gold and disordered carbon atoms showed that the effect of disorder significantly altered carbon ADP values and did not influence those of the gold atom.
Collapse
Affiliation(s)
- Sylwia Pawlędzio
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury 101, Warsaw 02-089, Poland
| | - Maura Malinska
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury 101, Warsaw 02-089, Poland
| | - Florian Kleemiss
- Faculty for Chemistry und Pharmacy, University of Regensburg, Universitätsstraße 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, Warsaw 02-089, Poland
| |
Collapse
|
15
|
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
|
16
|
Chodkiewicz M, Pawlędzio S, Woińska M, Woźniak K. Fragmentation and transferability in Hirshfeld atom refinement. IUCRJ 2022; 9:298-315. [PMID: 35371499 PMCID: PMC8895009 DOI: 10.1107/s2052252522000690] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/18/2021] [Accepted: 01/19/2022] [Indexed: 05/06/2023]
Abstract
Hirshfeld atom refinement (HAR) is one of the most effective methods for obtaining accurate structural parameters for hydrogen atoms from X-ray diffraction data. Unfortunately, it is also relatively computationally expensive, especially for larger molecules due to wavefunction calculations. Here, a fragmentation approach has been tested as a remedy for this problem. It gives an order of magnitude improvement in computation time for larger organic systems and is a few times faster for metal-organic systems at the cost of only minor differences in the calculated structural parameters when compared with the original HAR calculations. Fragmentation was also applied to polymeric and disordered systems where it provides a natural solution to problems that arise when HAR is applied. The concept of fragmentation is closely related to the transferable aspherical atom model (TAAM) and allows insight into possible ways to improve TAAM. Hybrid approaches combining fragmentation with the transfer of atomic densities between chemically similar atoms have been tested. An efficient handling of intermolecular interactions was also introduced for calculations involving fragmentation. When applied in fragHAR (a fragmentation approach for polypeptides) as a replacement for the original approach, it allowed for more efficient calculations. All of the calculations were performed with a locally modified version of Olex2 combined with a development version of discamb2tsc and ORCA. Care was taken to efficiently use the power of multicore processors by simple implementation of load-balancing, which was found to be very important for lowering computational time.
Collapse
Affiliation(s)
- Michał Chodkiewicz
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury 101, Warszawa 02-089, Poland
- Correspondence e-mail: ,
| | - Sylwia Pawlędzio
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury 101, Warszawa 02-089, Poland
| | - Magdalena Woińska
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury 101, Warszawa 02-089, Poland
| | - Krzysztof Woźniak
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury 101, Warszawa 02-089, Poland
- Correspondence e-mail: ,
| |
Collapse
|
17
|
Tchoń D, Makal A. Maximizing completeness in single-crystal high-pressure diffraction experiments: phase transitions in 2°AP. IUCRJ 2021; 8:1006-1017. [PMID: 34804552 PMCID: PMC8562673 DOI: 10.1107/s2052252521009532] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 07/22/2021] [Accepted: 09/14/2021] [Indexed: 06/13/2023]
Abstract
Sufficiently high completeness of diffraction data is necessary to correctly determine the space group, observe solid-state structural transformations or investigate charge density distribution under pressure. Regrettably, experiments performed at high pressure in a diamond anvil cell (DAC) yield inherently incomplete datasets. The present work systematizes the combined influence of radiation wavelength, DAC opening angle and sample orientation in a DAC on the completeness of diffraction data collected in a single-crystal high-pressure (HP) experiment with the help of dedicated software. In particular, the impact of the sample orientation on the achievable data completeness is quantified and proved to be substantial. Graphical guides for estimating the most beneficial sample orientation depending on the sample Laue class and assuming a few commonly used experimental setups are proposed. The usefulness of these guides has been tested in the case of luminescent 1,3-diacetylpyrene, suspected to undergo transitions from the α phase (Pnma) to the γ phase (Pn21 a) and δ phase (P1121/a) under pressure. Effective sample orientation has ensured over 90% coverage even for the monoclinic system and enabled unrestrained structure refinements and access to complete systematic extinction patterns.
Collapse
Affiliation(s)
- D. Tchoń
- Biological and Chemical Research Centre, Faculty of Chemistry, University of Warsaw, Żwirki i Wigury 101, 02-089 Warszawa, Poland
| | - A. Makal
- Biological and Chemical Research Centre, Faculty of Chemistry, University of Warsaw, Żwirki i Wigury 101, 02-089 Warszawa, Poland
| |
Collapse
|
18
|
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
|
19
|
Woińska M, Chodkiewicz ML, Woźniak K. Towards accurate and precise positions of hydrogen atoms bonded to heavy metal atoms. Chem Commun (Camb) 2021; 57:3652-3655. [PMID: 33870351 DOI: 10.1039/d0cc07661a] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A comparison of five X-ray structures of transition-metal-bound hydride complexes, successfully refined using Hirshfeld Atom Refinement (HAR) against low resolution X-ray diffraction data (including the positions and ADPs of all hydrogen atoms), with neutron structures shows that using aspherical atomic scattering factors instead of spherical ones results in systematic elongation of metal-hydrogen bonds, which in the case of the highest-quality data leads to excellent agreement of the X-ray and the neutron-derived bond lengths.
Collapse
Affiliation(s)
- Magdalena Woińska
- Biological and Chemical Research Centre, Chemistry Department, University of Warsaw, Żwirki i Wigury 101, 02-089 Warszawa, Poland.
| | | | | |
Collapse
|
20
|
Chodkiewicz ML, Woińska M, Woźniak K. Hirshfeld atom like refinement with alternative electron density partitions. IUCRJ 2020; 7:1199-1215. [PMID: 33209330 PMCID: PMC7642787 DOI: 10.1107/s2052252520013603] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/26/2020] [Accepted: 10/12/2020] [Indexed: 05/03/2023]
Abstract
Hirshfeld atom refinement is one of the most successful methods for the accurate determination of structural parameters for hydrogen atoms from X-ray diffraction data. This work introduces a generalization of the method [generalized atom refinement (GAR)], consisting of the application of various methods of partitioning electron density into atomic contributions. These were tested on three organic structures using the following partitions: Hirshfeld, iterative Hirshfeld, iterative stockholder, minimal basis iterative stockholder and Becke. The effects of partition choice were also compared with those caused by other factors such as quantum chemical methodology, basis set, representation of the crystal field and a combination of these factors. The differences between the partitions were small in terms of R factor (e.g. much smaller than for refinements with different quantum chemistry methods, i.e. Hartree-Fock and coupled cluster) and therefore no single partition was clearly the best in terms of experimental data reconstruction. In the case of structural parameters the differences between the partitions are comparable to those related to the choice of other factors. We have observed the systematic effects of the partition choice on bond lengths and ADP values of polar hydrogen atoms. The bond lengths were also systematically influenced by the choice of electron density calculation methodology. This suggests that GAR-derived structural parameters could be systematically improved by selecting an optimal combination of the partition and quantum chemistry method. The results of the refinements were compared with those of neutron diffraction experiments. This allowed a selection of the most promising partition methods for further optimization of GAR settings, namely the Hirshfeld, iterative stockholder and minimal basis iterative stockholder.
Collapse
Affiliation(s)
- Michał Leszek Chodkiewicz
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury 101, Warszawa, 02-089 Warszawa, Poland
| | - Magdalena Woińska
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury 101, Warszawa, 02-089 Warszawa, Poland
| | - Krzysztof Woźniak
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury 101, Warszawa, 02-089 Warszawa, Poland
| |
Collapse
|
21
|
Sanjuan-Szklarz WF, Woińska M, Domagała S, Dominiak PM, Grabowsky S, Jayatilaka D, Gutmann M, Woźniak K. On the accuracy and precision of X-ray and neutron diffraction results as a function of resolution and the electron density model. IUCRJ 2020; 7:920-933. [PMID: 32939284 PMCID: PMC7467170 DOI: 10.1107/s2052252520010441] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/29/2020] [Accepted: 07/28/2020] [Indexed: 05/20/2023]
Abstract
X-ray diffraction is the main source of three-dimensional structural information. In total, more than 1.5 million crystal structures have been refined and deposited in structural databanks (PDB, CSD and ICSD) to date. Almost 99.7% of them were obtained by approximating atoms as spheres within the independent atom model (IAM) introduced over a century ago. In this study, X-ray datasets for single crystals of hydrated α-oxalic acid were refined using several alternative electron density models that abandon the crude spherical approximation: the multipole model (MM), the transferable aspherical atom model (TAAM) and the Hirshfeld atom refinement (HAR) model as a function of the resolution of X-ray data. The aspherical models (MM, TAAM, HAR) give far more accurate and precise single-crystal X-ray results than IAM, sometimes identical to results obtained from neutron diffraction and at low resolution. Hence, aspherical approaches open new routes for improving existing structural information collected over the last century.
Collapse
Affiliation(s)
- W. Fabiola Sanjuan-Szklarz
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury, Warszawa, Poland
| | - Magdalena Woińska
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury, Warszawa, Poland
| | - Sławomir Domagała
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury, Warszawa, Poland
| | - Paulina M. Dominiak
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury, Warszawa, Poland
| | - Simon Grabowsky
- Department of Chemistry and Biochemistry, University of Bern, Freiestrasse 3, 3012 Bern, Switzerland
| | - Dylan Jayatilaka
- School of Molecular Sciences, University of Western Australia, 35 Stirling Highway, Perth WA 6009, Australia
| | - Matthias Gutmann
- Rutherford Appleton Laboratory, ISIS Facility, Chilton, Didcot, Oxfordshire OX11 OQX, United Kingdom
| | - Krzysztof Woźniak
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, Żwirki i Wigury, Warszawa, Poland
| |
Collapse
|
22
|
Jha KK, Gruza B, Kumar P, Chodkiewicz ML, Dominiak PM. TAAM: a reliable and user friendly tool for hydrogen-atom location using routine X-ray diffraction data. ACTA CRYSTALLOGRAPHICA SECTION B, STRUCTURAL SCIENCE, CRYSTAL ENGINEERING AND MATERIALS 2020; 76:296-306. [PMID: 32831250 DOI: 10.1107/s2052520620002917] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/05/2019] [Accepted: 03/02/2020] [Indexed: 06/11/2023]
Abstract
Hydrogen is present in almost all of the molecules in living things. It is very reactive and forms bonds with most of the elements, terminating their valences and enhancing their chemistry. X-ray diffraction is the most common method for structure determination. It depends on scattering of X-rays from electron density, which means the single electron of hydrogen is difficult to detect. Generally, neutron diffraction data are used to determine the accurate position of hydrogen atoms. However, the requirement for good quality single crystals, costly maintenance and the limited number of neutron diffraction facilities means that these kind of results are rarely available. Here it is shown that the use of Transferable Aspherical Atom Model (TAAM) instead of Independent Atom Model (IAM) in routine structure refinement with X-ray data is another possible solution which largely improves the precision and accuracy of X-H bond lengths and makes them comparable to averaged neutron bond lengths. TAAM, built from a pseudoatom databank, was used to determine the X-H bond lengths on 75 data sets for organic molecule crystals. TAAM parametrizations available in the modified University of Buffalo Databank (UBDB) of pseudoatoms applied through the DiSCaMB software library were used. The averaged bond lengths determined by TAAM refinements with X-ray diffraction data of atomic resolution (dmin ≤ 0.83 Å) showed very good agreement with neutron data, mostly within one single sample standard deviation, much like Hirshfeld atom refinement (HAR). Atomic displacements for both hydrogen and non-hydrogen atoms obtained from the refinements systematically differed from IAM results. Overall TAAM gave better fits to experimental data of standard resolution compared to IAM. The research was accompanied with development of software aimed at providing user-friendly tools to use aspherical atom models in refinement of organic molecules at speeds comparable to routine refinements based on spherical atom model.
Collapse
Affiliation(s)
- Kunal Kumar Jha
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, Warszawa, 02-089, Poland
| | - Barbara Gruza
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, Warszawa, 02-089, Poland
| | - Prashant Kumar
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, Warszawa, 02-089, Poland
| | - Michal Leszek Chodkiewicz
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, Warszawa, 02-089, Poland
| | - Paulina Maria Dominiak
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, Warszawa, 02-089, Poland
| |
Collapse
|
23
|
Gruza B, Chodkiewicz ML, Krzeszczakowska J, Dominiak PM. Refinement of organic crystal structures with multipolar electron scattering factors. Acta Crystallogr A Found Adv 2020; 76:92-109. [PMID: 31908353 PMCID: PMC8127334 DOI: 10.1107/s2053273319015304] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2019] [Accepted: 11/13/2019] [Indexed: 12/14/2022] Open
Abstract
A revolution in resolution is occurring now in electron microscopy arising from the development of methods for imaging single particles at cryogenic temperatures and obtaining electron diffraction data from nanocrystals of small organic molecules or macromolecules. Near-atomic or even atomic resolution of molecular structures can be achieved. The basis of these methods is the scattering of an electron beam due to the electrostatic potential of the sample. To analyse these high-quality experimental data, it is necessary to use appropriate atomic scattering factors. The independent atom model (IAM) is commonly used although various more advanced models, already known from X-ray diffraction, can also be applied to enhance the analysis. In this study a comparison is presented of IAM and TAAM (transferable aspherical atom model), the latter with the parameters of the Hansen-Coppens multipole model transferred from the University at Buffalo Databank (UBDB). By this method, TAAM takes into account the fact that atoms in molecules are partially charged and are not spherical. Structure refinements were performed on a carbamazepine crystal using electron structure-factor amplitudes determined experimentally [Jones et al. (2018). ACS Cent. Sci. 4, 1587-1592] or modelled with theoretical quantum-mechanical methods. The results show the possibilities and limitations of the TAAM method when applied to electron diffraction. Among others, the method clearly improves model fitting statistics, when compared with IAM, and allows for reliable refinement of atomic thermal parameters. The improvements are more pronounced with poorer-resolution diffraction data.
Collapse
Affiliation(s)
- Barbara Gruza
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, Warsaw, 02-089, Poland
| | - Michał Leszek Chodkiewicz
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, Warsaw, 02-089, Poland
| | - Joanna Krzeszczakowska
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, Warsaw, 02-089, Poland
| | - Paulina Maria Dominiak
- Biological and Chemical Research Centre, Department of Chemistry, University of Warsaw, ul. Żwirki i Wigury 101, Warsaw, 02-089, Poland
| |
Collapse
|
24
|
Cachau RE, Zhu J, Nicklaus MC. The upcoming subatomic resolution revolution. Curr Opin Struct Biol 2019; 58:53-58. [PMID: 31233975 DOI: 10.1016/j.sbi.2019.05.013] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2019] [Revised: 05/12/2019] [Accepted: 05/13/2019] [Indexed: 10/26/2022]
Abstract
Subatomic resolution macromolecular crystallography has been revealing the most fascinating details of macromolecular structures for many years. This most extreme form of macromolecular crystallography is going through rapid changes. A new generation of superbrilliant X-ray sources and detectors is facilitating the rapid acquisition of high-quality datasets. Equally important, a new breed of methods and highly integrated advanced computational tools for structure refinement and analysis is poised to change the way we use subatomic resolution data and reposition high-resolution macromolecular crystallography in medicinal chemistry studies. Subatomic resolution macromolecular crystallography may soon be a routine source of detailed molecular information besides precise geometries, including binding energies and other chemical descriptors, opening new possibilities of application.
Collapse
Affiliation(s)
- Raul E Cachau
- Advanced Biomedical Computational Science, Frederick National Laboratory for Cancer Research, Leidos Biomedical Inc., Frederick, MD 21702, USA.
| | - Jianghai Zhu
- Advanced Biomedical Computational Science, Frederick National Laboratory for Cancer Research, Leidos Biomedical Inc., Frederick, MD 21702, USA
| | - Marc C Nicklaus
- Chemical Biology Laboratory, National Cancer Institute, Frederick, MD 21702, USA
| |
Collapse
|