1
|
Zwitterionic or Not? Fast and Reliable Structure Determination by Combining Crystal Structure Prediction and Solid-State NMR. Molecules 2023; 28:molecules28041876. [PMID: 36838863 PMCID: PMC9966216 DOI: 10.3390/molecules28041876] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/16/2023] [Revised: 02/06/2023] [Accepted: 02/10/2023] [Indexed: 02/18/2023] Open
Abstract
When it comes to crystal structure determination, computational approaches such as Crystal Structure Prediction (CSP) have gained more and more attention since they offer some insight on how atoms and molecules are packed in the solid state, starting from only very basic information without diffraction data. Furthermore, it is well known that the coupling of CSP with solid-state NMR (SSNMR) greatly enhances the performance and the accuracy of the predictive method, leading to the so-called CSP-NMR crystallography (CSP-NMRX). In this paper, we present the successful application of CSP-NMRX to determine the crystal structure of three structural isomers of pyridine dicarboxylic acid, namely quinolinic, dipicolinic and dinicotinic acids, which can be in a zwitterionic form, or not, in the solid state. In a first step, mono- and bidimensional SSNMR spectra, i.e., 1H Magic-Angle Spinning (MAS), 13C and 15N Cross Polarisation Magic-Angle Spinning (CPMAS), 1H Double Quantum (DQ) MAS, 1H-13C HETeronuclear CORrelation (HETCOR), were used to determine the correct molecular structure (i.e., zwitterionic or not) and the local molecular arrangement; at the end, the RMSEs between experimental and computed 1H and 13C chemical shifts allowed the selection of the correct predicted structure for each system. Interestingly, while quinolinic and dipicolinic acids are zwitterionic and non-zwitterionic, respectively, in the solid state, dinicotinic acid exhibits in its crystal structure a "zwitterionic-non-zwitterionic continuum state" in which the proton is shared between the carboxylic moiety and the pyridinic nitrogen. Very refined SSNMR experiments were carried out, i.e., 14N-1H Phase-Modulated (PM) pulse and Rotational-Echo Saturation-Pulse Double-Resonance (RESPDOR), to provide an accurate N-H distance value confirming the hybrid nature of the molecule. The CSP-NMRX method showed a remarkable match between the selected structures and the experimental ones. The correct molecular input provided by SSNMR reduced the number of CSP calculations to be performed, leading to different predicted structures, while RMSEs provided an independent parameter with respect to the computed energy for the selection of the best candidate.
Collapse
|
2
|
Schlesinger C, Fitterer A, Buchsbaum C, Habermehl S, Chierotti MR, Nervi C, Schmidt MU. Ambiguous structure determination from powder data: four different structural models of 4,11-di-fluoro-quinacridone with similar X-ray powder patterns, fit to the PDF, SSNMR and DFT-D. IUCRJ 2022; 9:406-424. [PMID: 35844476 PMCID: PMC9252154 DOI: 10.1107/s2052252522004237] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/10/2022] [Accepted: 04/20/2022] [Indexed: 05/31/2023]
Abstract
Four different structural models, which all fit the same X-ray powder pattern, were obtained in the structure determination of 4,11-di-fluoro-quinacridone (C20H10N2O2F2) from unindexed X-ray powder data by a global fit. The models differ in their lattice parameters, space groups, Z, Z', molecular packing and hydrogen bond patterns. The molecules form a criss-cross pattern in models A and B, a layer structure built from chains in model C and a criss-cross arrangement of dimers in model D. Nevertheless, all models give a good Rietveld fit to the experimental powder pattern with acceptable R-values. All molecular geometries are reliable, except for model D, which is slightly distorted. All structures are crystallochemically plausible, concerning density, hydrogen bonds, intermolecular distances etc. All models passed the checkCIF test without major problems; only in model A a missed symmetry was detected. All structures could have probably been published, although 3 of the 4 structures were wrong. The investigation, which of the four structures is actually the correct one, was challenging. Six methods were used: (1) Rietveld refinements, (2) fit of the crystal structures to the pair distribution function (PDF) including the refinement of lattice parameters and atomic coordinates, (3) evaluation of the colour, (4) lattice-energy minimizations with force fields, (5) lattice-energy minimizations by two dispersion-corrected density functional theory methods, and (6) multinuclear CPMAS solid-state NMR spectroscopy (1H, 13C, 19F) including the comparison of calculated and experimental chemical shifts. All in all, model B (perhaps with some disorder) can probably be considered to be the correct one. This work shows that a structure determination from limited-quality powder data may result in totally different structural models, which all may be correct or wrong, even if they are chemically sensible and give a good Rietveld refinement. Additionally, the work is an excellent example that the refinement of an organic crystal structure can be successfully performed by a fit to the PDF, and the combination of computed and experimental solid-state NMR chemical shifts can provide further information for the selection of the most reliable structure among several possibilities.
Collapse
Affiliation(s)
- Carina Schlesinger
- Institute of Inorganic and Analytical Chemistry, Johann Wolfgang Goethe University, Max-von-Laue-Straße 7, 60438 Frankfurt am Main, Germany
| | - Arnd Fitterer
- Institute of Inorganic and Analytical Chemistry, Johann Wolfgang Goethe University, Max-von-Laue-Straße 7, 60438 Frankfurt am Main, Germany
| | - Christian Buchsbaum
- Institute of Inorganic and Analytical Chemistry, Johann Wolfgang Goethe University, Max-von-Laue-Straße 7, 60438 Frankfurt am Main, Germany
| | - Stefan Habermehl
- Institute of Inorganic and Analytical Chemistry, Johann Wolfgang Goethe University, Max-von-Laue-Straße 7, 60438 Frankfurt am Main, Germany
| | - Michele R. Chierotti
- Department of Chemistry and NIS centre, University of Torino, V. Giuria 7, Torino 10125, Italy
| | - Carlo Nervi
- Department of Chemistry and NIS centre, University of Torino, V. Giuria 7, Torino 10125, Italy
| | - Martin U. Schmidt
- Institute of Inorganic and Analytical Chemistry, Johann Wolfgang Goethe University, Max-von-Laue-Straße 7, 60438 Frankfurt am Main, Germany
| |
Collapse
|
3
|
Bravetti F, Bordignon S, Alig E, Eisenbeil D, Fink L, Nervi C, Gobetto R, Schmidt MU, Chierotti MR. Solid-State NMR-Driven Crystal Structure Prediction of Molecular Crystals: The Case of Mebendazole. Chemistry 2021; 28:e202103589. [PMID: 34962330 DOI: 10.1002/chem.202103589] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2021] [Indexed: 11/06/2022]
Abstract
Among all possible NMR crystallography approaches for crystal-structure determination, crystal structure prediction - NMR crystallography (CSP-NMRX) has recently turned out to be a powerful method. In the latter, the original procedure exploited solid-state NMR (SSNMR) information during the final steps of the prediction. In particular, it used the comparison of computed and experimental chemical shifts for the selection of the correct crystal packing. Still, the prediction procedure, generally carried out with DFT methods, may require important computational resources and be quite time-consuming, especially if there are no available constraints to use at the initial stage. Herein, the successful application of this combined prediction method, which exploits NMR information also in the input step to reduce the search space of the predictive algorithm, is presented. Herein, this method was applied on mebendazole, which is characterized by desmotropism. The use of SSNMR data as constraints for the selection of the right tautomer and the determination of the number of independent molecules in the unit cell led to a considerably faster process, reducing the number of calculations to be performed. In this way, the crystal packing was successfully predicted for the three known phases of mebendazole. To evaluate the quality of the predicted structures, these were compared to the experimental ones. The crystal structure of phase B of mebendazole, in particular, was determined de novo by powder diffraction and is presented for the first time in this paper.
Collapse
Affiliation(s)
- Federica Bravetti
- Department of Chemistry, Università degli Studi di Torino, via Pietro Giuria 7, 10125, Torino, Italy
| | - Simone Bordignon
- Department of Chemistry, Università degli Studi di Torino, via Pietro Giuria 7, 10125, Torino, Italy
| | - Edith Alig
- Institute of Inorganic and Analytical Chemistry, Goethe University, Max-von-Laue-Strasse 7, 60438, Frankfurt am Main, Germany
| | - Daniel Eisenbeil
- Institute of Inorganic and Analytical Chemistry, Goethe University, Max-von-Laue-Strasse 7, 60438, Frankfurt am Main, Germany
| | - Lothar Fink
- Institute of Inorganic and Analytical Chemistry, Goethe University, Max-von-Laue-Strasse 7, 60438, Frankfurt am Main, Germany
| | - Carlo Nervi
- Department of Chemistry, Università degli Studi di Torino, via Pietro Giuria 7, 10125, Torino, Italy
| | - Roberto Gobetto
- Department of Chemistry, Università degli Studi di Torino, via Pietro Giuria 7, 10125, Torino, Italy
| | - Martin U Schmidt
- Institute of Inorganic and Analytical Chemistry, Goethe University, Max-von-Laue-Strasse 7, 60438, Frankfurt am Main, Germany
| | - Michele R Chierotti
- Department of Chemistry, Università degli Studi di Torino, via Pietro Giuria 7, 10125, Torino, Italy
| |
Collapse
|
4
|
Mathew R, Stevensson B, Edén M. Refined Structures of O-Phospho-l-serine and Its Calcium Salt by New Multinuclear Solid-State NMR Crystallography Methods. J Phys Chem B 2021; 125:10985-11004. [PMID: 34553936 PMCID: PMC8503883 DOI: 10.1021/acs.jpcb.1c05587] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2021] [Indexed: 12/17/2022]
Abstract
O-phospho-l-serine (Pser) and its Ca salt, Ca[O-phospho-l-serine]·H2O (CaPser), play important roles for bone mineralization and were recently also proposed to account for the markedly improved bone-adhesive properties of Pser-doped calcium phosphate-based cements for biomedical implants. However, the hitherto few proposed structural models of Pser and CaPser were obtained by X-ray diffraction, thereby leaving the proton positions poorly defined. Herein, we refine the Pser and CaPser structures by density functional theory (DFT) calculations and contrast them with direct interatomic-distance constraints from two-dimensional (2D) nuclear magnetic resonance (NMR) correlation experimentation at fast magic-angle spinning (MAS), encompassing double-quantum-single-quantum (2Q-1Q) 1H NMR along with heteronuclear 13C{1H} and 31P{1H} correlation NMR experiments. The Pser and CaPser structures before and after refinements by DFT were validated against sets of NMR-derived effective 1H-1H, 1H-31P, and 1H-13C distances, which confirmed the improved accuracy of the refined structures. Each distance set was derived from one sole 2D NMR experiment applied to a powder without isotopic enrichment. The distances were extracted without invoking numerical spin-dynamics simulations or approximate phenomenological models. We highlight the advantages and limitations of the new distance-extraction procedure. Isotropic 1H, 13C, and 31P chemical shifts obtained by DFT calculations using the gauge including projector augmented wave (GIPAW) method agreed very well with the experimental results. We discuss the isotropic and anisotropic 13C and 31P chemical-shift parameters in relation to the previous literature, where most data on CaPser are reported herein for the first time.
Collapse
Affiliation(s)
- Renny Mathew
- Department of Materials and Environmental
Chemistry, Stockholm University, SE-106 91 Stockholm, Sweden
| | - Baltzar Stevensson
- Department of Materials and Environmental
Chemistry, Stockholm University, SE-106 91 Stockholm, Sweden
| | - Mattias Edén
- Department of Materials and Environmental
Chemistry, Stockholm University, SE-106 91 Stockholm, Sweden
| |
Collapse
|
5
|
Pawlak T, Sudgen I, Bujacz G, Iuga D, Brown SP, Potrzebowski MJ. Synergy of Solid-State NMR, Single-Crystal X-ray Diffraction, and Crystal Structure Prediction Methods: A Case Study of Teriflunomide (TFM). CRYSTAL GROWTH & DESIGN 2021; 21:3328-3343. [PMID: 34267599 PMCID: PMC8273857 DOI: 10.1021/acs.cgd.1c00123] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 02/01/2021] [Revised: 04/27/2021] [Indexed: 06/13/2023]
Abstract
In this work, for the first time, we present the X-ray diffraction crystal structure and spectral properties of a new, room-temperature polymorph of teriflunomide (TFM), CSD code 1969989. As revealed by DSC, the low-temperature TFM polymorph recently reported by Gunnam et al. undergoes a reversible thermal transition at -40 °C. This reversible process is related to a change in Z' value, from 2 to 1, as observed by variable-temperature 1H-13C cross-polarization (CP) magic-angle spinning (MAS) solid-state NMR, while the crystallographic system is preserved (triclinic). Two-dimensional 13C-1H and 1H-1H double-quantum MAS NMR spectra are consistent with the new room-temperature structure, including comparison with GIPAW (gauge-including projector augmented waves) calculated NMR chemical shifts. A crystal structure prediction procedure found both experimental teriflunomide polymorphs in the energetic global minimum region. Differences between the polymorphs are seen for the torsional angle describing the orientation of the phenyl ring relative to the planarity of the TFM molecule. In the low-temperature structure, there are two torsion angles of 4.5 and 31.9° for the two Z' = 2 molecules, while in the room-temperature structure, there is disorder that is modeled with ∼50% occupancy between torsion angles of -7.8 and 28.6°. These observations are consistent with a broad energy minimum as revealed by DFT calculations. PISEMA solid-state NMR experiments show a reduction in the C-H dipolar coupling in comparison to the static limit for the aromatic CH moieties of 75% and 51% at 20 and 40 °C, respectively, that is indicative of ring flips at the higher temperature. Our study shows the power of combining experiments, namely DSC, X-ray diffraction, and MAS NMR, with DFT calculations and CSP to probe and understand the solid-state landscape, and in particular the role of dynamics, for pharmaceutical molecules.
Collapse
Affiliation(s)
- Tomasz Pawlak
- Centre
of Molecular and Macromolecular Studies, Polish Academy of Sciences, Sienkiewicza 112, 90-363 Lodz, Poland
| | - Isaac Sudgen
- Molecular
Systems Engineering Group, Centre for Process Systems Engineering,
Department of Chemical Engineering, Imperial
College London, London SW7 2AZ, U.K.
| | - Grzegorz Bujacz
- Institute
of Molecular and Industrial Biotechnology, Lodz University of Technology, Stefanowskiego 4/10, 90-924, Lodz, Poland
| | - Dinu Iuga
- Department
of Physics, University of Warwick, Coventry CV4 7AL, U.K.
| | - Steven P. Brown
- Department
of Physics, University of Warwick, Coventry CV4 7AL, U.K.
| | - Marek J. Potrzebowski
- Centre
of Molecular and Macromolecular Studies, Polish Academy of Sciences, Sienkiewicza 112, 90-363 Lodz, Poland
| |
Collapse
|
6
|
Hodgkinson P. NMR crystallography of molecular organics. PROGRESS IN NUCLEAR MAGNETIC RESONANCE SPECTROSCOPY 2020; 118-119:10-53. [PMID: 32883448 DOI: 10.1016/j.pnmrs.2020.03.001] [Citation(s) in RCA: 65] [Impact Index Per Article: 16.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/02/2019] [Revised: 02/25/2020] [Accepted: 03/13/2020] [Indexed: 06/11/2023]
Abstract
Developments of NMR methodology to characterise the structures of molecular organic structures are reviewed, concentrating on the previous decade of research in which density functional theory-based calculations of NMR parameters in periodic solids have become widespread. With a focus on demonstrating the new structural insights provided, it is shown how "NMR crystallography" has been used in a spectrum of applications from resolving ambiguities in diffraction-derived structures (such as hydrogen atom positioning) to deriving complete structures in the absence of diffraction data. As well as comprehensively reviewing applications, the different aspects of the experimental and computational techniques used in NMR crystallography are surveyed. NMR crystallography is seen to be a rapidly maturing subject area that is increasingly appreciated by the wider crystallographic community.
Collapse
Affiliation(s)
- Paul Hodgkinson
- Department of Chemistry, Durham University, Stockton Road, Durham DH1 3LE, UK.
| |
Collapse
|
7
|
Dudek MK, Paluch P, Śniechowska J, Nartowski KP, Day GM, Potrzebowski MJ. Crystal structure determination of an elusive methanol solvate – hydrate of catechin using crystal structure prediction and NMR crystallography. CrystEngComm 2020. [DOI: 10.1039/d0ce00452a] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Abstract
A useful short-cut was developed to limit the number of molecular conformations that need to be regarded in crystal structure prediction calculations, which led to the crystal structure determination of new methanol solvate – hydrate of catechin.
Collapse
Affiliation(s)
- Marta K. Dudek
- Centre of Molecular and Macromolecular Studies of Polish Academy of Sciences
- 90-363 Lodz
- Poland
| | - Piotr Paluch
- Centre of Molecular and Macromolecular Studies of Polish Academy of Sciences
- 90-363 Lodz
- Poland
| | - Justyna Śniechowska
- Centre of Molecular and Macromolecular Studies of Polish Academy of Sciences
- 90-363 Lodz
- Poland
| | - Karol P. Nartowski
- Department of Drug Form Technology
- Wroclaw Medical University
- 50-556 Wroclaw
- Poland
| | - Graeme M. Day
- Computational Systems Chemistry
- School of Chemistry
- University of Southampton
- UK
| | - Marek J. Potrzebowski
- Centre of Molecular and Macromolecular Studies of Polish Academy of Sciences
- 90-363 Lodz
- Poland
| |
Collapse
|
8
|
Hofstetter A, Balodis M, Paruzzo FM, Widdifield CM, Stevanato G, Pinon AC, Bygrave PJ, Day GM, Emsley L. Rapid Structure Determination of Molecular Solids Using Chemical Shifts Directed by Unambiguous Prior Constraints. J Am Chem Soc 2019; 141:16624-16634. [PMID: 31117663 PMCID: PMC7540916 DOI: 10.1021/jacs.9b03908] [Citation(s) in RCA: 33] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Abstract
NMR-based crystallography approaches involving the combination of crystal structure prediction methods, ab initio calculated chemical shifts and solid-state NMR experiments are powerful methods for crystal structure determination of microcrystalline powders. However, currently structural information obtained from solid-state NMR is usually included only after a set of candidate crystal structures has already been independently generated, starting from a set of single-molecule conformations. Here, we show with the case of ampicillin that this can lead to failure of structure determination. We propose a crystal structure determination method that includes experimental constraints during conformer selection. In order to overcome the problem that experimental measurements on the crystalline samples are not obviously translatable to restrict the single-molecule conformational space, we propose constraints based on the analysis of absent cross-peaks in solid-state NMR correlation experiments. We show that these absences provide unambiguous structural constraints on both the crystal structure and the gas-phase conformations, and therefore can be used for unambiguous selection. The approach is parametrized on the crystal structure determination of flutamide, flufenamic acid, and cocaine, where we reduce the computational cost by around 50%. Most importantly, the method is then shown to correctly determine the crystal structure of ampicillin, which would have failed using current methods because it adopts a high-energy conformer in its crystal structure. The average positional RMSE on the NMR powder structure is ⟨rav⟩ = 0.176 Å, which corresponds to an average equivalent displacement parameter Ueq = 0.0103 Å2.
Collapse
Affiliation(s)
- Albert Hofstetter
- Institut des Sciences et Ingénierie Chimiques , École Polytechnique Fédérale de Lausanne (EPFL) , 1015 Lausanne , Switzerland
| | - Martins Balodis
- Institut des Sciences et Ingénierie Chimiques , École Polytechnique Fédérale de Lausanne (EPFL) , 1015 Lausanne , Switzerland
| | - Federico M Paruzzo
- Institut des Sciences et Ingénierie Chimiques , École Polytechnique Fédérale de Lausanne (EPFL) , 1015 Lausanne , Switzerland
| | - Cory M Widdifield
- Department of Chemistry, Mathematics and Science Center , Oakland University , 146 Library Drive , Rochester , Michigan 48309-4479 , United States
| | - Gabriele Stevanato
- Institut des Sciences et Ingénierie Chimiques , École Polytechnique Fédérale de Lausanne (EPFL) , 1015 Lausanne , Switzerland
| | - Arthur C Pinon
- Institut des Sciences et Ingénierie Chimiques , École Polytechnique Fédérale de Lausanne (EPFL) , 1015 Lausanne , Switzerland
| | - Peter J Bygrave
- School of Chemistry , University of Southampton , Highfield , Southampton SO17 1BJ , United Kingdom
| | - Graeme M Day
- School of Chemistry , University of Southampton , Highfield , Southampton SO17 1BJ , United Kingdom
| | - Lyndon Emsley
- Institut des Sciences et Ingénierie Chimiques , École Polytechnique Fédérale de Lausanne (EPFL) , 1015 Lausanne , Switzerland
| |
Collapse
|