NMR crystallography of campho[2,3-c]pyrazole (Z' = 6): combining high-resolution 1H-13C solid-state MAS NMR spectroscopy and GIPAW chemical-shift calculations

An NMR crystallographic characterisation of solid (+)-usnic acid

We use a combination of one- and two-dimensional solid-state nuclear magnetic resonance (NMR) spectroscopy and density functional theory (DFT) calculations to obtain a full assignment of the 1H and 13C signals for solid (+)-usnic acid, which contains two molecules in the asymmetric unit. By combining through-space 1H-1H correlation data with computation it is possible to assign signals not just to the same molecules (relative assignment) but to assign the signals to specific crystallographic molecules (absolute assignment). Variable-temperature measurements reveal that there is some variation in many of the 13C chemical shifts with temperature, likely arising from varying populations of different tautomeric forms of the molecule. The NMR spectrum of crystalline (+)-usnic acid is then compared with that of ground Usnea dasopoga lichen (the source material of the usnic acid). The abundance of usnic acid is so great in the lichen that this natural product can be observed directly in the NMR spectrum without further purification. This natural sample of usnic acid appears to have the same crystalline form as that in the pure commercial sample.

Polymorph Identification for Flexible Molecules: Linear Regression Analysis of Experimental and Calculated Solution- and Solid-State NMR Data

The Δδ regression approach of Blade et al. [ J. Phys. Chem. A 2020, 124(43), 8959-8977] for accurately discriminating between solid forms using a combination of experimental solution- and solid-state NMR data with density functional theory (DFT) calculation is here extended to molecules with multiple conformational degrees of freedom, using furosemide polymorphs as an exemplar. As before, the differences in measured 1H and 13C chemical shifts between solution-state NMR and solid-state magic-angle spinning (MAS) NMR (Δδexperimental) are compared to those determined by gauge-including projector augmented wave (GIPAW) calculations (Δδcalculated) by regression analysis and a t-test, allowing the correct furosemide polymorph to be precisely identified. Monte Carlo random sampling is used to calculate solution-state NMR chemical shifts, reducing computation times by avoiding the need to systematically sample the multidimensional conformational landscape that furosemide occupies in solution. The solvent conditions should be chosen to match the molecule's charge state between the solution and solid states. The Δδ regression approach indicates whether or not correlations between Δδexperimental and Δδcalculated are statistically significant; the approach is differently sensitive to the popular root mean squared error (RMSE) method, being shown to exhibit a much greater dynamic range. An alternative method for estimating solution-state NMR chemical shifts by approximating the measured solution-state dynamic 3D behavior with an ensemble of 54 furosemide crystal structures (polymorphs and cocrystals) from the Cambridge Structural Database (CSD) was also successful in this case, suggesting new avenues for this method that may overcome its current dependency on the prior determination of solution dynamic 3D structures.

Discovering the Solid-State Secrets of Lorlatinib by NMR Crystallography: To Hydrogen Bond or not to Hydrogen Bond

Lorlatinib is an active pharmaceutical ingredient (API) used in the treatment of lung cancer. Here, an NMR crystallography analysis is presented whereby the single-crystal X-ray diffraction structure (CSD: 2205098) determination is complemented by multinuclear (1H, 13C, 14/15N, 19F) magic-angle spinning (MAS) solid-state NMR and gauge-including projector augmented wave (GIPAW) calculation of NMR chemical shifts. Lorlatinib crystallises in the P21 space group, with two distinct molecules in the asymmetric unit cell, Z' = 2. Three of the four NH2 hydrogen atoms form intermolecular hydrogen bonds, N30-H…N15 between the two distinct molecules and N30-H…O2 between two equivalent molecules. This is reflected in one of the NH21H chemical shifts being significantly lower, 4.0 ppm compared to 7.0 ppm. Two-dimensional 1H-13C, 14N-1H and 1H (double-quantum, DQ)-1H (single-quantum, SQ) MAS NMR spectra are presented. The 1H resonances are assigned and specific HH proximities corresponding to the observed DQ peaks are identified. The resolution enhancement at a 1H Larmor frequency of 1 GHz as compared to 500 or 600 MHz is demonstrated.

Exploring the Potential of Multinuclear Solid-State 1 H, 13 C, and 35 Cl Magnetic Resonance To Characterize Static and Dynamic Disorder in Pharmaceutical Hydrochlorides

Crystallographic disorder, whether static or dynamic, can be detrimental to the physical and chemical stability, ease of crystallization and dissolution rate of an active pharmaceutical ingredient. Disorder can result in a loss of manufacturing control leading to batch-to-batch variability and can lengthen the process of structural characterization. The range of NMR active nuclei makes solid-state NMR a unique technique for gaining nucleus-specific information about crystallographic disorder. Here, we explore the use of high-field ³⁵ Cl solid-state NMR at 23.5 T to characterize both static and dynamic crystallographic disorder: specifically, dynamic disorder occurring in duloxetine hydrochloride (1), static disorder in promethazine hydrochloride (2), and trifluoperazine dihydrochloride (3). In all structures, the presence of crystallographic disorder was confirmed by ¹³ C cross-polarization magic-angle spinning (CPMAS) NMR and supported by GIPAW-DFT calculations, and in the case of 3, ¹ H solid-state NMR provided additional confirmation. Applying ³⁵ Cl solid-state NMR to these compounds, we show that higher magnetic fields are beneficial for resolving the crystallographic disorder in 1 and 3, while broad spectral features were observed in 2 even at higher fields. Combining the data obtained from ¹ H, ¹³ C, and ³⁵ Cl NMR, we show that 3 exhibits a unique case of disorder involving the ⁺ N-H hydrogen positions of the piperazinium ring, driving the chloride anions to occupy three distinct sites.

Resolving the Chemical Formula of Nesquehonite via NMR Crystallography, DFT Computation, and Complementary Neutron Diffraction

Nesquehonite is a magnesium carbonate mineral relevant to carbon sequestration envisioned for carbon capture and storage of CO₂ . Its chemical formula remains controversial today, assigned as either a hydrated magnesium carbonate [MgCO₃  ⋅ 3H₂ O], or a hydroxy- hydrated- magnesium bicarbonate [Mg(HCO₃ )OH ⋅ 2H₂ O]. The resolution of this controversy is central to understanding this material's thermodynamic, phase, and chemical behavior. In an NMR crystallography study, using rotational-echo double-resonance ¹³ C{¹ H} (REDOR), ¹³ C-¹ H distances are determined with precision, and the combination of ¹³ C static NMR lineshapes and density functional theory (DFT) calculations are used to model different H atomic coordinates. [MgCO₃  ⋅ 3H₂ O] is found to be accurate, and evidence from neutron powder diffraction bolsters these assignments. Refined H positions can help understand how H-bonding stabilizes this structure against dehydration to MgCO₃ . More broadly, these results illustrate the power of NMR crystallography as a technique for resolving questions where X-ray diffraction is inconclusive.

Combining heteronuclear correlation NMR with spin-diffusion to detect relayed Cl-H-H and N-H-H proximities in molecular solids

Analysis of short-to-intermediate range intermolecular interactions offers a great way of characterizing the solid-state organization of small molecules and materials. This can be achieved by two-dimensional (2D) homo- and heteronuclear correlation NMR spectroscopy, for example, by carrying out experiments at high magnetic fields in conjunction with fast magic-angle spinning (MAS) techniques. But, detecting 2D peaks for heteronuclear dipolar coupled spin pairs separated by greater than 3 Å is not always straightforward, particularly when low-gamma quadrupolar nuclei are involved. Here, we present a 2D correlation NMR experiment that combines the advantages of heteronuclear-multiple quantum coherence (HMQC) and proton-based spin-diffusion (SD) pulse sequences using radio-frequency-driven-recouping (RFDR) to probe inter and intramolecular ¹H-X (X = ¹⁴N, ³⁵Cl) interactions. This experiment can be used to acquire 2D ¹H{X}-HMQC filtered ¹H-¹H correlation as well as 2D ¹H-X HMQC spectra. Powder forms of dopamine·HCl and l-histidine·HCl·H₂O are characterized at high fields (21.1 T and 18.8 T) with fast MAS (60 kHz) using the 2D HMQC-SD-RFDR approach. Solid-state NMR results are complemented with NMR crystallography analyses using the gauge-including projector augmented wave (GIPAW) approach. For histidine·HCl·H₂O, 2D peaks associated with ¹⁴N-¹H-¹H and ³⁵Cl-¹H-¹H distances of up to 4.4 and 3.9 Å have been detected. This is further corroborated by the observation of 2D peaks corresponding to ¹⁴N-¹H-¹H and ³⁵Cl-¹H-¹H distances of up to 4.2 and 3.7 Å in dopamine·HCl, indicating the suitability of the HMQC-SD-RFDR experiments for detecting medium-range proximities in molecular solids.

Complete resonance assignment of a pharmaceutical drug at natural isotopic abundance from DNP-Enhanced solid-state NMR

Solid-state dynamic nuclear polarization enhanced magic angle spinning (DNP-MAS) NMR measurements coupled with density functional theory (DFT) calculations enable the full resonance assignment of a complex pharmaceutical drug molecule without the need for isotopic enrichment. DNP dramatically enhances the NMR signals, thereby making possible previously intractable two-dimensional correlation NMR spectra at natural abundance. Using inputs from DFT calculations, herein we describe a significant improvement to the structure elucidation process for complex organic molecules. Further, we demonstrate that a series of two-dimensional correlation experiments, including ¹⁵N-¹³C TEDOR, ¹³C-¹³C INADEQUATE/SARCOSY, ¹⁹F-¹³C HETCOR, and ¹H-¹³C HETCOR, can be obtained at natural isotopic abundance within reasonable experiment times, thus enabling a complete resonance assignment of sitagliptin, a pharmaceutical used for the treatment of type 2 diabetes.

Solid-State Study of the Structure, Dynamics, and Thermal Processes of Safinamide Mesylate─A New Generation Drug for the Treatment of Neurodegenerative Diseases

Safinamide mesylate (SM), the pure active pharmaceutical ingredient (API) recently used in Parkinson disease treatment, recrystallized employing water–ethanol mixture of solvents (vol/vol 1:9) gives a different crystallographic form compared to SM in Xadago tablets. Pure SM crystallizes as a hemihydrate in the monoclinic system with the P2₁ space group. Its crystal and molecular structure were determined by means of cryo X-ray crystallography at 100 K. SM in the Xadago tablet exists in anhydrous form in the orthorhombic crystallographic system with the P2₁2₁2₁ space group. The water migration and thermal processes in the crystal lattice were monitored by solid-state NMR spectroscopy, differential scanning calorimetry, and thermogravimetric analysis. SM in Xadago in the high-humidity environment undergoes phase transformation to the P2₁ form which can be easily reversed just by heating up to 80 °C. For the commercial form of the API, there is also a reversible thermal transformation observed between Z′ = 1 ↔ Z′ = 3 crystallographic forms in the 0–20 °C temperature range. Analysis of molecular motion in the crystal lattice proves that the observed conformational polymorphism is forced by intramolecular dynamics. All above-mentioned processes were analyzed and described employing the NMR crystallography approach with the support of advanced theoretical calculations.

A toolbox for improving the workflow of NMR crystallography

NMR crystallography is a powerful tool with applications in structural characterization and crystal structure verification, to name two. However, applying this tool presents several challenges, especially for industrial users, in terms of consistency, workflow, time consumption, and the requirement for a high level of understanding of experimental solid-state NMR and GIPAW-DFT calculations. Here, we have developed a series of fully parameterized scripts for use in Materials Studio and TopSpin, based on the .magres file format, with a focus on organic molecules (e.g. pharmaceuticals), improving efficiency, robustness, and workflow. We separate these tools into three major categories: performing the DFT calculations, extracting & visualizing the results, and crystallographic modelling. These scripts will rapidly submit fully parameterized CASTEP jobs, extract data from the calculations, assist in visualizing the results, and expedite the process of structural modelling. Accompanied with these tools is a description on their functionality, documentation on how to get started and use the scripts, and links to video tutorials for guiding new users. Through the use of these tools, we hope to facilitate NMR crystallography and to harmonize the process across users.

Single-Crystal X-ray and Solid-State NMR Characterisation of AND-1184 and Its Hydrochloride Form

In this study, we report on a structural investigation of AND-1184, with the chemical name N-[3-[4-(6-fluoro-1,2-benzoxazol-3-yl)piperidin-1-yl]propyl]-3-methylbenzenesulfonamide (MBS), and its hydrochloride form (MBS^HCl); AND-1184 is a potential API for the treatment of dementia. The single-crystal X-ray investigation of both forms results in monoclinic crystal systems with P2₁/c and C2/c symmetry for MBS and MBS^HCl, respectively. This solid-state NMR study, combined with quantum-chemical calculations, allowed us to assign all ¹³C and most ¹H signals. The MBS structure was defined as a completely rigid system without significant dynamic behaviours, whereas MBS^HCl exhibited limited dynamic motion of the aromatic part of the molecule.

Accurate fragment-based 51-V chemical shift predictions in molecular crystals

Nuclear magnetic resonance (NMR) spectroscopy plays a crucial role in determining molecular structure for complex biological and pharmaceutical compounds. NMR investigations are increasingly reliant on computation for mapping spectral features to chemical structures. Here we benchmark the accuracy of fragment-based ⁵¹V chemical shielding tensor calculations using a training set comprised of 10 biologically and pharmaceutically relevant oxovanadium complexes. Using our self-consistent reproduction of the Madelung potential (SCRMP) electrostatic embedding model, we demonstrate comparable performance between fragment methods and computationally demanding cluster-based techniques. Specifically, fragment methods employing hybrid density functionals are capable of reproducing the experimental ⁵¹V isotropic chemical shifts with a training set rms error of ~9 ​ppm, representing a 20% improvement over traditional plane wave techniques. We provide training set-derived linear regression models for mapping the absolute shieldings obtained from computation to the experimentally determined chemical shifts using four common density functionals; PBE0, B3LYP, PBE, and BLYP. Finally, we establish the utility of fragment methods and the reported regression parameters examining four oxovanadium structures excluded from the training set including the tetracoordinate oxovanadium silicate [Formula: see text] , VO¹⁵NGlySalbz which contains redox-active ligands, and the solid-state form of the common ⁵¹V NMR reference compound VOCl₃.

Pseudocyclic Form of 4-Hydroxypyrrolidine-2-carboxanilide Podands with Trioxyethylene Chain: Modeling, Conformational Search, and NMR Analysis

The 4-hydroxypyrrolidine-2-carboxanilide podand salt demonstrates catalytic activity in asymmetric Biginelli reaction. The systematic search for prevalent conformational state of the cation was carried out by computer simulations in combination with one- and two-dimensional NMR experiments. For that purpose, we proposed a novel algorithm for the generation and selection of conformers based on molecular dynamics and clustering in the space of principal components. The search had found an important trend of the podand to form a pseudocyclic structure with a horseshoe-shaped conformation of the oligooxyethylene fragment. This conformation is stabilized by different types of intramolecular hydrogen bonds between the acidic and basic centers of the two 4-hydroxypyrrolidine-2-carboxanilide residuals (branches). The proposed approach had made it possible to identify the major structural factors, providing a correlation between the calculated and experimental chemical shifts of hydrogen atoms in the ¹H NMR spectra of the protonated podand.

Synergy of Solid-State NMR, Single-Crystal X-ray Diffraction, and Crystal Structure Prediction Methods: A Case Study of Teriflunomide (TFM)

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 ¹H-¹³C cross-polarization (CP) magic-angle spinning (MAS) solid-state NMR, while the crystallographic system is preserved (triclinic). Two-dimensional ¹³C-¹H and ¹H-¹H 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.

Modeling Small Structural and Environmental Differences in Solids with 15 N NMR Chemical Shift Tensors

The ability to theoretically predict accurate NMR chemical shifts in solids is increasingly important due to the role such shifts play in selecting among proposed model structures. Herein, two theoretical methods are evaluated for their ability to assign ¹⁵ N shifts from guanosine dihydrate to one of the two independent molecules present in the lattice. The NMR data consist of ¹⁵ N shift tensors from 10 resonances. Analysis using periodic boundary or fragment methods consider a benchmark dataset to estimate errors and predict uncertainties of 5.6 and 6.2 ppm, respectively. Despite this high accuracy, only one of the five sites were confidently assigned to a specific molecule of the asymmetric unit. This limitation is not due to negligible differences in experimental data, as most sites exhibit differences of >6.0 ppm between pairs of resonances representing a given position. Instead, the theoretical methods are insufficiently accurate to make assignments at most positions.

A structural analysis of 2,5-diaryl-4H-2,4-dihydro-3H-1,2,4-triazol-3-ones: NMR in the solid state, X-ray crystallography, and GIPAW calculations

The ¹ H, ¹³ C, ¹⁵ N, and ¹⁹ F nuclear magnetic resonance (NMR) spectra of 11 2,5-diaryl-2,4-dihydro-3H-1,2,4-triazol-3-ones have been acquired in DMSO-d₆ solution and the ¹³ C, ¹⁵ N, and ¹⁹ F NMR spectra have also been acquired in the solid state (solid-state nuclear magnetic resonance [SSNMR] and magic angle spinning [MAS]). The X-ray structures of Compounds 3, 5, and 6 have been determined by X-ray diffraction. Theoretical calculations at the gauge-independent atomic orbital (GIAO)/B3LYP/6-311++G(d,p) level have provided a set of 321 chemical shifts that were compared with 310 experimental values in DMSO-d₆ . To obtain good agreements, some effects need to be included. The SSNMR chemical shifts have been compared with gauge-including projector-augmented wave (GIPAW) calculations and with the heavy atom-light atom (HALA) effects.

Hydrogen-Mediated Noncovalent Interactions in Solids: What Can NMR Crystallography Tell About?

Hydrogen atoms play a crucial role in the aggregation of organic (bio)molecules through diverse number of noncovalent interactions that they mediate, such as electrostatic in proton transfer systems, hydrogen bonding, and CH-π interactions, to mention only the most prominent. To identify and adequately describe such low-energy interactions, increasingly sensitive methods have been developed over time, among which quantum chemical computations have witnessed impressive advances in recent years. For reaching the present state-of-the-art, computations had to rely on a pool of relevant experimental data, needed at least for validation, if not also for other purposes. In the case of molecular crystals, the best illustration for the synergy between computations and experiment is given by the so-called NMR crystallography approach. Originally designed to increase the confidence level in crystal structure determination of organic compounds from powders, NMR crystallography is able now to offer also a wealth of information regarding the noncovalent interactions that drive molecules to pack in a given crystalline pattern or another. This is particularly true for the noncovalent interactions which depend on the exact location of labile hydrogen atoms in the system: in such cases, NMR crystallography represents a valuable characterization tool, in some cases complementing even the standard single-crystal X-ray diffraction technique. A concise introduction in the field is made in this mini-review, which is aimed at providing a comprehensive picture with respect to the current accuracy level reached by NMR crystallography in the characterization of hydrogen-mediated noncovalent interactions in organic solids. Different types of practical applications are illustrated with the example of molecular crystals studied by our research group, but references to other representative developments reported in the literature are also made. By summarizing the major concepts and methodological progresses, the present work is also intended to be a guide to the practical potential of this relatively recent analytical tool for the scientists working in areas where crystal engineering represents the main approach for rational design of novel materials.

Investigating discrepancies between experimental solid-state NMR and GIPAW calculation: NC-N 13C and OH⋯O 1H chemical shifts in pyridinium fumarates and their cocrystals

An NMR crystallography analysis is presented for four solid-state structures of pyridine fumarates and their cocrystals, using crystal structures deposited in the Cambridge Crystallographic Data Centre, CCDC. Experimental one-dimensional one-pulse ¹H and ¹³C cross-polarisation (CP) magic-angle spinning (MAS) nuclear magnetic resonance (NMR) and two-dimensional ¹⁴N-¹H heteronuclear multiple-quantum coherence MAS NMR spectra are compared with gauge-including projector augmented wave (GIPAW) calculations of the ¹H and ¹³C chemical shifts and the ¹⁴N shifts that additionally depend on the quadrupolar interaction. Considering the high ppm (>10 ​ppm) ¹H resonances, while there is good agreement (within 0.4 ​ppm) between experiment and GIPAW calculation for the hydrogen-bonded NH moieties, the hydrogen-bonded fumaric acid OH resonances are 1.2-1.9 ​ppm higher in GIPAW calculation as compared to experiment. For the cocrystals of a salt and a salt formed by 2-amino-5-methylpyridinium and 2-amino-6-methylpyridinium ions, a large discrepancy of 4.2 and 5.9 ​ppm between experiment and GIPAW calculation is observed for the quaternary ring carbon ¹³C resonance that is directly bonded to two nitrogens (in the ring and in the amino group). By comparison, there is excellent agreement (within 0.2 ​ppm) for the quaternary ring carbon ¹³C resonance directly bonded to the ring nitrogen for the salt and cocrystal of a salt formed by 2,6-lutidinium and 2,5-lutidinium, respectively.

A GIPAW versus GIAO-ZORA-SO study of 13C and 15N CPMAS NMR chemical shifts of aromatic and heterocyclic bromo derivatives

Theoretical simulation of NMR parameters in compounds bearing heavy atoms generally requires the application of relativistic corrections. We report herein the theoretical characterization of ¹³C and ¹⁵N CPMAS NMR of known bromo-derivative crystals by using both the GIPAW and the combined GIAO-ZORA-SO approximation methods. Several statistical analyses were performed to compare both approaches, with non-relativistic GIPAW method being more useful to predict the ¹³C and ¹⁵N chemical shifts. The problem of applying GIPAW to crystal structures showing static or dynamic crystalline disorder of the special class resulting in half-protons will be discussed in detail.

NMR crystallography of molecular organics

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.

A Bayesian approach to NMR crystal structure determination

Nuclear Magnetic Resonance (NMR) spectroscopy is particularly well suited to determine the structure of molecules and materials in powdered form. Structure determination usually proceeds by finding the best match between experimentally observed NMR chemical shifts and those of candidate structures. Chemical shifts for the candidate configurations have traditionally been computed by electronic-structure methods, and more recently predicted by machine learning. However, the reliability of the determination depends on the errors in the predicted shifts. Here we propose a Bayesian framework for determining the confidence in the identification of the experimental crystal structure, based on knowledge of the typical errors in the electronic structure methods. We demonstrate the approach on the determination of the structures of six organic molecular crystals. We critically assess the reliability of the structure determinations, facilitated by the introduction of a visualization of the similarity between candidate configurations in terms of their chemical shifts and their structures. We also show that the commonly used values for the errors in calculated ¹³C shifts are underestimated, and that more accurate, self-consistently determined uncertainties make it possible to use ¹³C shifts to improve the accuracy of structure determinations. Finally, we extend the recently-developed ShiftML model to render it more efficient, accurate, and, most importantly, to evaluate the uncertainties in its predictions. By quantifying the confidence in structure determinations based on ShiftML predictions we further substantiate that it provides a valid replacement for first-principles calculations in NMR crystallography.

Understanding the formation of apremilast cocrystals

Apremilast (APR), an anti-psoriatic agent, easily forms isostructural cocrystals and solvates with aromatic entities, often disobeying at the same time Kitaigorodsky's rule as to the saturation of possible hydrogen-bonding sites. In this paper the reasons for this peculiar behavior are investigated, employing a joint experimental and theoretical approach. This includes the design of cocrystals with coformers having a high propensity towards the formation of both aromatic-aromatic and hydrogen-bonding interactions, determination of their structure, using solid-state NMR spectroscopy and X-ray crystallography, as well as calculations of stabilization energies of formation of the obtained cocrystals, followed by crystal structure prediction calculations and solubility measurements. The findings indicate that the stabilization energies of cocrystal formation are positive in all cases, which results from strain in the APR conformation in these crystal forms. On the other hand, solubility measurements show that the Gibbs free energy of formation of the apremilast:picolinamide cocrystal is negative, suggesting that the formation of the studied cocrystals is entropy driven. This entropic stabilization is associated with the disorder observed in almost all known cocrystals and solvates of APR.

Rapid Structure Determination of Molecular Solids Using Chemical Shifts Directed by Unambiguous Prior Constraints

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 ⟨r_av⟩ = 0.176 Å, which corresponds to an average equivalent displacement parameter U_eq = 0.0103 Ų.

The use of variable temperature 13 C solid-state MAS NMR and GIPAW DFT calculations to explore the dynamics of diethylcarbamazine citrate

Experimental ¹³ C solid-state magic-angle spinning (MAS) Nuclear Magnetic Resonance (NMR) as well as Density-Functional Theory (DFT) gauge-including projector augmented wave (GIPAW) calculations were used to probe disorder and local mobility in diethylcarbamazine citrate, (DEC)⁺ (citrate)^- . This compound has been used as the first option drug for the treatment of filariasis, a disease endemic in tropical countries and caused by adult worms of Wuchereria bancrofti, which is transmitted by mosquitoes. We firstly present 2D ¹³ C─¹ H dipolar-coupling-mediated heteronuclear correlation spectra recorded at moderate spinning frequency, to explore the intermolecular interaction between DEC and citrate molecules. Secondly, we investigate the dynamic behavior of (DEC)⁺ (citrate)^- by varying the temperature and correlating the experimental MAS NMR results with DFT GIPAW calculations that consider two (DEC)⁺ conformers (in a 70:30 ratio) for crystal structures determined at 293 and 235 K. Solid-state NMR provides insights on slow exchange dynamics revealing conformational changes involving particularly the DEC ethyl groups.

An NMR crystallography investigation of furosemide

This paper presents an NMR crystallography study of three polymorphs of furosemide. Experimental magic-angle spinning (MAS) solid-state NMR spectra are reported for form I of furosemide, and these are assigned using density-functional theory (DFT)-based gauge-including projector augmented wave (GIPAW) calculations. Focusing on the three known polymorphs, we examine the changes to the NMR parameters due to crystal packing effects. We use a recently developed formalism to visualise which regions are responsible for the chemical shielding of particular sites and hence understand the variation in NMR parameters between the three polymorphs.

NMR Crystallography: Evaluation of Hydrogen Positions in Hydromagnesite by 13 C{1 H} REDOR Solid-State NMR and Density Functional Theory Calculation of Chemical Shielding Tensors

Solid-state NMR measurements coupled with density functional theory (DFT) calculations demonstrate how hydrogen positions can be refined in a crystalline system. The precision afforded by rotational-echo double-resonance (REDOR) NMR to interrogate ¹³ C-¹ H distances is exploited along with DFT determinations of the ¹³ C tensor of carbonates (CO₃ ^2- ). Nearby ¹ H nuclei perturb the axial symmetry of the carbonate sites in the hydrated carbonate mineral, hydromagnesite [4 MgCO₃ ⋅Mg(OH)₂ ⋅4 H₂ O]. A match between the calculated structure and solid-state NMR was found by testing multiple semi-local and dispersion-corrected DFT functionals and applying them to optimize atom positions, starting from X-ray diffraction (XRD)-determined atomic coordinates. This was validated by comparing calculated to experimental ¹³ C{¹ H} REDOR and ¹³ C chemical shift anisotropy (CSA) tensor values. The results show that the combination of solid-state NMR, XRD, and DFT can improve structure refinement for hydrated materials.

An XRD and NMR crystallographic investigation of the structure of 2,6-lutidinium hydrogen fumarate

A crystallographic study highlighting the benefits of a combined XRD and NMR approach in investigating both stability and variation within an organic multicomponent crystal.

Chemical shifts in molecular solids by machine learning

Due to their strong dependence on local atonic environments, NMR chemical shifts are among the most powerful tools for strucutre elucidation of powdered solids or amorphous materials. Unfortunately, using them for structure determination depends on the ability to calculate them, which  comes at the cost of high accuracy first-principles calculations. Machine learning has recently emerged as a way to overcome the need for quantum chemical calculations, but for chemical shifts in solids it is hindered by the chemical and combinatorial space spanned by molecular solids, the strong dependency of chemical shifts on their environment, and the lack of an experimental database of shifts. We propose a machine learning method based on local environments to accurately predict chemical shifts of molecular solids and their polymorphs to within DFT accuracy. We also demonstrate that the trained model is able to determine, based on the match between experimentally measured and ML-predicted shifts, the structures of cocaine and the drug 4-[4-(2-adamantylcarbamoyl)-5-tert-butylpyrazol-1-yl]benzoic acid.

Solid-state nuclear magnetic resonance combined with quantum chemical shift predictions is limited by high computational cost. Here, the authors use machine learning based on local atomic environments to predict experimental chemical shifts in molecular solids with accuracy similar to density functional theory.

19F-NMR Diastereotopic Signals in Two N-CHF₂ Derivatives of (4S,7R)-7,8,8-Trimethyl-4,5,6,7-tetrahydro-4,7-methano-2H-indazole

In this paper, we report the anisochrony of the fluorine atoms of a CHF₂ group when linked to a pyrazole ring. The pyrazole is part of (4S,7R)-7,8,8-trimethyl-4,5,6,7-tetrahydro-4,7-methano-2H-indazole also known as (4S,7R)-campho[2,3-c]pyrazole, which has two stereogenic centers. Gauge-Independent Atomic Orbital (GIAO)/Becke, 3-parameter, Lee-Yang-Parr (B3LYP)/6-311++G(d,f) calculated 19F chemical shifts of the minimum energy conformations satisfactorily agree with the experimental data. The energy differences between minima need to consider solvent effects (continuum model) to be satisfactorily reproduced.

Ab initio random structure searching of organic molecular solids: assessment and validation against experimental data

This paper explores the capability of using the DFT-D ab initio random structure searching (AIRSS) method to generate crystal structures of organic molecular materials, focusing on a system (m-aminobenzoic acid; m-ABA) that is known from experimental studies to exhibit abundant polymorphism. Within the structural constraints selected for the AIRSS calculations (specifically, centrosymmetric structures with Z = 4 for zwitterionic m-ABA molecules), the method is shown to successfully generate the two known polymorphs of m-ABA (form III and form IV) that have these structural features. We highlight various issues that are encountered in comparing crystal structures generated by AIRSS to experimental powder X-ray diffraction (XRD) data and solid-state magic-angle spinning (MAS) NMR data, demonstrating successful fitting for some of the lowest energy structures from the AIRSS calculations against experimental low-temperature powder XRD data for known polymorphs of m-ABA, and showing that comparison of computed and experimental solid-state NMR parameters allows different hydrogen-bonding motifs to be discriminated.

Probing intermolecular interactions in a diethylcarbamazine citrate salt by fast MAS 1H solid-state NMR spectroscopy and GIPAW calculations

Fast magic-angle spinning (MAS) NMR is used to probe intermolecular interactions in a diethylcarbamazine salt, that is widely used as a treatment against adult worms of Wuchereria bancrofti which cause a common disease in tropical countries named filariasis. Specifically, a dihydrogen citrate salt that has improved thermal stability and solubility as compared to the free form is studied. One-dimensional ¹H, ¹³C and ¹⁵N and two-dimensional ¹H-¹³C and ¹⁴N-¹H heteronuclear correlation NMR experiments under moderate and fast MAS together with GIPAW (CASTEP) calculations enable the assignment of the ¹H, ¹³C and ¹⁴N/¹⁵N resonances. A two-dimensional ¹H-¹H double-quantum (DQ) -single-quantum (SQ) MAS spectrum recorded with BaBa recoupling at 60kHz MAS identifies specific proton-proton proximities associated with citrate-citrate and citrate-diethylcarbamazine intermolecular interactions.

Determination of a complex crystal structure in the absence of single crystals: analysis of powder X-ray diffraction data, guided by solid-state NMR and periodic DFT calculations, reveals a new 2'-deoxyguanosine structural motif

Derivatives of guanine exhibit diverse supramolecular chemistry, with a variety of distinct hydrogen-bonding motifs reported in the solid state, including ribbons and quartets, which resemble the G-quadruplex found in nucleic acids with sequences rich in guanine. Reflecting this diversity, the solid-state structural properties of 3',5'-bis-O-decanoyl-2'-deoxyguanosine, reported in this paper, reveal a hydrogen-bonded guanine ribbon motif that has not been observed previously for 2'-deoxyguanosine derivatives. In this case, structure determination was carried out directly from powder XRD data, representing one of the most challenging organic molecular structures (a 90-atom molecule) that has been solved to date by this technique. While specific challenges were encountered in the structure determination process, a successful outcome was achieved by augmenting the powder XRD analysis with information derived from solid-state NMR data and with dispersion-corrected periodic DFT calculations for structure optimization. The synergy of experimental and computational methodologies demonstrated in the present work is likely to be an essential feature of strategies to further expand the application of powder XRD as a technique for structure determination of organic molecular materials of even greater complexity in the future.

Single-crystal X-ray diffraction and NMR crystallography of a 1:1 cocrystal of dithianon and pyrimethanil

A single-crystal X-ray diffraction structure of a 1:1 cocrystal of two fungicides, namely dithianon (DI) and pyrimethanil (PM), is reported [systematic name: 5,10-dioxo-5H,10H-naphtho[2,3-b][1,4]dithiine-2,3-dicarbonitrile-4,6-dimethyl-N-phenylpyrimidin-2-amine (1/1), C14H4N2O2S2·C12H13N2]. Following an NMR crystallography approach, experimental solid-state magic angle spinning (MAS) NMR spectra are presented together with GIPAW (gauge-including projector augmented wave) calculations of NMR chemical shieldings. Specifically, experimental 1H and 13C chemical shifts are determined from two-dimensional 1H-13C MAS NMR correlation spectra recorded with short and longer contact times so as to probe one-bond C-H connectivities and longer-range C...H proximities, whereas H...H proximities are identified in a 1H double-quantum (DQ) MAS NMR spectrum. The performing of separate GIPAW calculations for the full periodic crystal structure and for isolated molecules allows the determination of the change in chemical shift upon going from an isolated molecule to the full crystal structure. For the 1H NMR chemical shifts, changes of 3.6 and 2.0 ppm correspond to intermolecular N-H...O and C-H...O hydrogen bonding, while changes of -2.7 and -1.5 ppm are due to ring current effects associated with C-H...π interactions. Even though there is a close intermolecular S...O distance of 3.10 Å, it is of note that the molecule-to-crystal chemical shifts for the involved sulfur or oxygen nuclei are small.

Ab initio crystal structure prediction of magnesium (poly)sulfides and calculation of their NMR parameters

Ab initio prediction of sensible crystal structures can be regarded as a crucial task in the quickly-developing methodology of NMR crystallography. In this contribution, an evolutionary algorithm was used for the prediction of magnesium (poly)sulfide crystal structures with various compositions. The employed approach successfully identified all three experimentally detected forms of MgS, i.e. the stable rocksalt form and the metastable wurtzite and zincblende forms. Among magnesium polysulfides with a higher content of sulfur, the most probable structure with the lowest formation energy was found to be MgS2, exhibiting a modified rocksalt structure, in which S2− anions were replaced by S2 2− dianions. Magnesium polysulfides with even larger fractions of sulfur were not predicted to be stable. For the lowest-energy structures, 25Mg quadrupolar coupling constants and chemical shift parameters were calculated using the density functional theory approach. The calculated NMR parameters could be well rationalized by the symmetries of the local magnesium environments, by the coordination of magnesium cations and by the nature of the surrounding anions. In the future, these parameters could serve as a reference for the experimentally determined 25Mg NMR parameters of magnesium sulfide species.

Enhanced NMR Discrimination of Pharmaceutically Relevant Molecular Crystal Forms through Fragment-Based Ab Initio Chemical Shift Predictions

Chemical shift prediction plays an important role in the determination or validation of crystal structures with solid-state nuclear magnetic resonance (NMR) spectroscopy. One of the fundamental theoretical challenges lies in discriminating variations in chemical shifts resulting from different crystallographic environments. Fragment-based electronic structure methods provide an alternative to the widely used plane wave gauge-including projector augmented wave (GIPAW) density functional technique for chemical shift prediction. Fragment methods allow hybrid density functionals to be employed routinely in chemical shift prediction, and we have recently demonstrated appreciable improvements in the accuracy of the predicted shifts when using the hybrid PBE0 functional instead of generalized gradient approximation (GGA) functionals like PBE. Here, we investigate the solid-state ¹³C and ¹⁵N NMR spectra for multiple crystal forms of acetaminophen, phenobarbital, and testosterone. We demonstrate that the use of the hybrid density functional instead of a GGA provides both higher accuracy in the chemical shifts and increased discrimination among the different crystallographic environments. Finally, these results also provide compelling evidence for the transferability of the linear regression parameters mapping predicted chemical shieldings to chemical shifts that were derived in an earlier study.

Crystallographic and Dynamic Aspects of Solid-State NMR Calibration Compounds: Towards ab Initio NMR Crystallography

The excellent results of dispersion-corrected density functional theory (DFT-D) calculations for static systems have been well established over the past decade. The introduction of dynamics into DFT-D calculations is a target, especially for the field of molecular NMR crystallography. Four (13) C ss-NMR calibration compounds are investigated by single-crystal X-ray diffraction, molecular dynamics and DFT-D calculations. The crystal structure of 3-methylglutaric acid is reported. The rotator phases of adamantane and hexamethylbenzene at room temperature are successfully reproduced in the molecular dynamics simulations. The calculated (13) C chemical shifts of these compounds are in excellent agreement with experiment, with a root-mean-square deviation of 2.0 ppm. It is confirmed that a combination of classical molecular dynamics and DFT-D chemical shift calculation improves the accuracy of calculated chemical shifts.

The role of solid-state nuclear magnetic resonance in crystal engineering

This Highlight article discusses the role of solid-state NMR spectroscopy in crystal engineering with the aid of several examples from the literature.

Polymorphs of Theophylline Characterized by DNP Enhanced Solid-State NMR

We show how dynamic nuclear polarization (DNP) enhanced solid-state NMR spectroscopy can be used to characterize polymorphs and solvates of organic solids. We applied DNP to three polymorphs and one hydrated form of the asthma drug molecule theophylline. For some forms of theophylline, sample grinding and impregnation with the radical-containing solution, which are necessary to prepare the samples for DNP, were found to induce polymorphic transitions or desolvation between some forms. We present protocols for sample preparation for solid-state magic-angle spinning (MAS) DNP experiments that avoid the polymorphic phase transitions in theophylline. These protocols include cryogrinding, grinding under inert atmosphere, and the appropriate choice of the impregnating liquid. By applying these procedures, we subsequently demonstrate that two-dimensional correlation experiments, such as (1)H-(13)C and (1)H-(15)N HETCOR or (13)C-(13)C INADEQUATE, can be obtained at natural isotopic abundance in reasonable times, thus enabling more advanced structural characterization of polymorphs.

Converging nuclear magnetic shielding calculations with respect to basis and system size in protein systems

Ab initio chemical shielding calculations greatly facilitate the interpretation of nuclear magnetic resonance (NMR) chemical shifts in biological systems, but the large sizes of these systems requires approximations in the chemical models used to represent them. Achieving good convergence in the predicted chemical shieldings is necessary before one can unravel how other complex structural and dynamical factors affect the NMR measurements. Here, we investigate how to balance trade-offs between using a better basis set or a larger cluster model for predicting the chemical shieldings of the substrates in two representative examples of protein-substrate systems involving different domains in tryptophan synthase: the N-(4'-trifluoromethoxybenzoyl)-2-aminoethyl phosphate (F9) ligand which binds in the α active site, and the 2-aminophenol quinonoid intermediate formed in the β active site. We first demonstrate that a chemically intuitive three-layer, locally dense basis model that uses a large basis on the substrate, a medium triple-zeta basis to describe its hydrogen-bonding partners and/or surrounding van der Waals cavity, and a crude basis set for more distant atoms provides chemical shieldings in good agreement with much more expensive large basis calculations. Second, long-range quantum mechanical interactions are important, and one can accurately estimate them as a small-basis correction to larger-basis calculations on a smaller cluster. The combination of these approaches enables one to perform density functional theory NMR chemical shift calculations in protein systems that are well-converged with respect to both basis set and cluster size.

An NMR crystallography study of the hemihydrate of 2', 3'-O-isopropylidineguanosine

An NMR crystallography study of the hemihydrate of 2', 3'-O-isopropylidineguanosine (Gace) is presented, together with powder X-ray diffraction and thermogravimetric analysis. (1)H double-quantum and (14)N-(1)H HMQC spectra recorded at 850MHz and 75kHz MAS (using a JEOL 1mm probe) are presented together with a (1)H-(13)C refocused INEPT spectrum recorded at 500MHz and 12.5kHz MAS using eDUMBO-122(1)H homonuclear decoupling. NMR chemical shieldings are calculated using the GIPAW (gauge-including projector augmented wave) method; good two-dimensional agreement between calculation and experiment is observed for (13)C and (1)H chemical shifts for directly bonded CH and CH3 peaks. There are two Gace molecules in the asymmetric unit cell: differences in specific (1)H chemical shifts are rationalised in terms of the strength of CH-π and intermolecular hydrogen bonding interactions.

NMR crystallography of ezetimibe co-crystals

An efficient, simplified protocol for solvent-drop assisted co-crystal preparation of ezetimibe (a drug for the treatment of primary hypercholesterolemia) with both imidazole and l-proline has been derived. The structures of the white powders were successfully solved via "NMR crystallography" combining solid-state NMR, powder X-ray diffraction and DFT chemical shift computations. Detailed insights into the likely crystallization mechanism were obtained from competition experiments, where efficient co-crystallization was feasible using ezetimibe monohydrate as precursor indicating that the crystal water acts as "molecular catalyst". It was also found that co-crystallization of imidazole is favored over l-proline, thus suggesting a clear preference of neutral hydrogen bonds compared to charge-assisted motifs.