NMR Crystallography of the Polymorphs of Metergoline

1H isotropic chemical shift metrics for NMR crystallography of powdered molecular organics

Hydrogen magnetic shielding values from gauge including projector augmented wave (GIPAW) density functional theory (DFT) calculations, when combined with experimental solid-state 1H nuclear magnetic resonance (NMR) chemical shift data collected on powdered microcrystalline organics, have been used to perform various crystal structure characterization tasks (e.g., refinements, verifications, determinations). These tasks fall under the umbrella of 'NMR crystallography'. In several instances, an isotropic 1H chemical shift (δiso) root-mean-squared deviation (RMSD) metric has been applied during these studies (including the first de novo crystal structure determination: M. Baias, J.-N. Dumez, P. H. Svensson, S. Schantz, G. M. Day and L. Emsley, De Novo Determination of the Crystal Structure of a Large Drug Molecule by Crystal Structure Prediction-Based Powder NMR Crystallography, J. Am. Chem. Soc., 2013, 135, 17501-17507). While it is assumed that the 1H δiso RMSD metrics are converged, our study probes the robustness of these metrics. Specifically, we consider how the structure of the δiso(1H) RMSD metric varies depending on: (i) selected GIPAW DFT input parameters; (ii) the number of fitting parameters used during linear mapping; and (iii) the GIPAW DFT computational software. These δiso(1H) RMSD metrics were produced from a set of 24 benchmark crystal structures (428 crystallographically unique hydrogen atom environments). Interestingly, we find that the δiso(1H) RMSD metric structures are very robust to substantial degradation in the quality of the GIPAW DFT computations, which is unexplored in the NMR crystallography literature as prior studies focus on convergence rather than divergence. We then briefly consider the impact of our findings using the structure determination of thymol as an illustrative example and our results strongly suggest that if δiso(1H) RMSD metrics are being used, then the GIPAW DFT computations can be performed much more efficiently than at present. Overall, this should allow for more efficient NMR crystallography characterization tasks of important materials that contain 1H nuclei, such as organic pharmaceuticals.

Anisotropy, segmental dynamics and polymorphism of crystalline biogenic carboxylic acids

Carboxylic acids of the Krebs cycle possess invaluable biochemical significance. Still, there are severe gaps in the availability of thermodynamic and crystallographic data, as well as ambiguities prevailing in the literature on the thermodynamic characterization and polymorph ranking. Providing an unambiguous description of the structure, thermodynamics and polymorphism of their neat crystalline phases requires a complex multidisciplinary approach. This work presents results of an extensive investigation of the structural anisotropy of the thermal expansion and local dynamics within these crystals, obtained from a beneficial cooperation of NMR crystallography and ab initio calculations of non-covalent interactions. The observed structural anisotropy and spin-lattice relaxation times are traced to large spatial variations in the strength of molecular interactions in the crystal lattice, especially in the orientation of the hydrogen bonds. A completely resolved crystal structure for oxaloacetic acid is reported for the first time. Thanks to multi-instrumental calorimetric effort, this work clarifies phase behavior, determines third-law entropies of the crystals, and states definitive polymorph ranking for succinic and fumaric acids. These thermodynamic observations are then interpreted in terms of first-principles quasi-harmonic calculations of cohesive properties. A sophisticated model capturing electronic, thermal, and configurational-entropic effects on the crystal structure approaches captures the subtle Gibbs energy differences governing polymorph ranking for succinic and fumaric acids, representing another success story of computational chemistry.

Combining X-ray and NMR Crystallography to Explore the Crystallographic Disorder in Salbutamol Oxalate

Salbutamol is an active pharmaceutical ingredient commonly used to treat respiratory distress and is listed by the World Health Organization as an essential medicine. Here, we establish the crystal structure of its oxalate form, salbutamol oxalate, and explore the nature of its crystallographic disorder by combined X-ray crystallography and ¹³C cross-polarization (CP) magic-angle spinning (MAS) solid-state NMR. The *C-OH chiral center of salbutamol (note that the crystal structures are a racemic mixture of the two enantiomers of salbutamol) is disordered over two positions, and the tert-butyl group is rotating rapidly, as revealed by ¹³C solid-state NMR. The impact of crystallization conditions on the disorder was investigated, finding variations in the occupancy ratio of the *C-OH chiral center between single crystals and a consistency across samples in the bulk powder. Overall, this work highlights the contrast between investigating crystallographic disorder by X-ray diffraction and solid-state NMR experiment, and gauge-including projector-augmented-wave (GIPAW) density functional theory (DFT) calculations, with their combined use, yielding an improved understanding of the nature of the crystallographic disorder between the local (i.e., as viewed by NMR) and longer-range periodic (i.e., as viewed by diffraction) scale.

Enantiotropy of Simvastatin as a Result of Weakened Interactions in the Crystal Lattice: Entropy-Driven Double Transitions and the Transient Modulated Phase as Seen by Solid-State NMR Spectroscopy

In crystalline molecular solids, in the absence of strong intermolecular interactions, entropy-driven processes play a key role in the formation of dynamically modulated transient phases. Specifically, in crystalline simvastatin, the observed fully reversible enantiotropic behavior is associated with multiple order–disorder transitions: upon cooling, the dynamically disordered high-temperature polymorphic Form I is transformed to the completely ordered low-temperature polymorphic Form III via the intermediate (transient) modulated phase II. This behavior is associated with a significant reduction in the kinetic energy of the rotating and flipping ester substituents, as well as a decrease in structural ordering into two distinct positions. In transient phase II, the conventional three-dimensional structure is modulated by periodic distortions caused by cooperative conformation exchange of the ester substituent between the two states, which is enabled by weakened hydrogen bonding. Based on solid-state NMR data analysis, the mechanism of the enantiotropic phase transition and the presence of the transient modulated phase are documented.

Polymorphic Forms of Valinomycin Investigated by NMR Crystallography

A dodecadepsipeptide valinomycin (VLM) has been most recently reported to be a potential anti-coronavirus drug that could be efficiently produced on a large scale. It is thus of importance to study solid-phase forms of VLM in order to be able to ensure its polymorphic purity in drug formulations. The previously available solid-state NMR (SSNMR) data are combined with the plane-wave DFT computations in the NMR crystallography framework. Structural/spectroscopical predictions (the PBE functional/GIPAW method) are obtained to characterize four polymorphs of VLM. Interactions which confer a conformational stability to VLM molecules in these crystalline forms are described in detail. The way how various structural factors affect the values of SSNMR parameters is thoroughly analyzed, and several SSNMR markers of the respective VLM polymorphs are identified. The markers are connected to hydrogen bonding effects upon the corresponding (13C/15N/1H) isotropic chemical shifts of (CO, Namid, Hamid, Hα) VLM backbone nuclei. These results are expected to be crucial for polymorph control of VLM and in probing its interactions in dosage forms.

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.

Monitoring the Site-Specific Solid-State NMR Data in Oligopeptides

Reliable values of the solid-state NMR (SSNMR) parameters together with precise structural data specific for a given amino acid site in an oligopeptide are needed for the proper interpretation of measurements aiming at an understanding of oligopeptides' function. The periodic density functional theory (DFT)-based computations of geometries and SSNMR chemical shielding tensors (CSTs) of solids are shown to be accurate enough to support the SSNMR investigations of suitably chosen models of oriented samples of oligopeptides. This finding is based on a thorough comparison between the DFT and experimental data for a set of tripeptides with both ¹³C_α and ¹⁵N_amid CSTs available from the single-crystal SSNMR measurements and covering the three most common secondary structural elements of polypeptides. Thus, the ground is laid for a quantitative description of local spectral parameters of crystalline oligopeptides, as demonstrated for the backbone ¹⁵N_amid nuclei of samarosporin I, which is a pentadecapeptide (composed of five classical and ten nonproteinogenic amino acids) featuring a strong antimicrobial activity.

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.

A theoretical NMR study of polymorphism in crystal structures of azoles and benzazoles

The NMR chemical shifts of two azoles and one benzazole whose crystal structures present polymorphism have been computed using the GIPAW approach. ¹⁵ N and ¹³ C nuclei have been studied. Statistical analysis of the computed ¹³ C and ¹⁵ N chemical shifts indicates that the GIPAW chemical shifts reproduce with a high degree of accuracy those experimentally reported. This methodology can be used to identify other polymorphic crystal structures.

Exploring Accuracy Limits of Predictions of the 1H NMR Chemical Shielding Anisotropy in the Solid State

The 1H chemical shielding anisotropy (CSA) is an NMR parameter that is exquisitely sensitive to the local environment of protons in crystalline systems, but it is difficult to obtain it experimentally due to the need to concomitantly suppress other anisotropic interactions in the solid-state NMR (SSNMR) pulse sequences. The SSNMR measurements of the 1H CSA are particularly challenging if the fast magic-angle-spinning (MAS) is applied. It is thus important to confront the results of both the single-crystal (SC) and fast-MAS experiments with their theoretical counterparts. Here the plane-waves (PW) DFT calculations have been carried out using two functionals in order to precisely characterize the structures and the 1H NMR chemical shielding tensors (CSTs) of the solid forms of maleic, malonic, and citric acids, and of L-histidine hydrochloride monohydrate. The level of agreement between the PW DFT and either SC or fast-MAS SSNMR 1H CSA data has been critically compared. It has been found that for the eigenvalues of the 1H CSTs provided by the fast-MAS measurements, an accuracy limit of current PW DFT predictions is about two ppm in terms of the standard deviation of the linear regression model, and sources of this error have been thoroughly discussed.