Conformations in Solution and in Solid-State Polymorphs: Correlating Experimental and Calculated Nuclear Magnetic Resonance Chemical Shifts for Tolfenamic Acid

A new approach for quantitively assessing putative crystal structures with applications in crystal structure prediction (CSP) is introduced that is based upon experimental solution- and magic-angle spinning (MAS) solid-state NMR data and density functional theory (DFT) calculation. For the specific case of tolfenamic acid (TFA), we consider experimental solution-state NMR for a range of solvents, experimental MAS NMR of polymorphs I and II, and DFT calculations for four polymorphs. The change in NMR chemical shift observed in passing from the solution state to the solid state (Delta delta(Experimental)) is calculated as the difference between H-1 and C-13 experimental solid-state chemical shifts for each polymorphic form (delta(Solid expt)) and the corresponding solution-state NMR chemical shifts (delta(Solution expt)). Separately, we use the gauge-included projector augmented wave (GIPAW) method to calculate the NMR chemical shifts for each form (delta(Solid calc)) and for TFA in solution (delta(Solution calc)) using the dynamic 3D solution conformational ensemble determined from NMR spectroscopy. The calculated change in passing from the solution state to the solid state (Delta delta(Calculated)) is then calculated as the difference of delta(Solid calc) and delta(Solution calc). Regression analysis for Delta delta(Calculated )against Delta delta(Experimental) followed by a t-test for statistical significance provides a robust quantitative assessment. We show that this assessment clearly identifies the correct polymorph, i.e., when comparing Delta delta(Experimental )based on the experimental MAS NMR chemical shifts of form I or II with Delta delta(Calculated )based on calculated chemical shifts for polymorphs I, II, III, and IV. Complementarity to the established approach of comparing delta(Solid expt) to delta(Solid calc) is explored. We further show that our approach is applicable if there are no solid-state crystal structure data. Specifically, delta(Solid calc) in Delta delta(Calculated) is replaced by the chemical shift for an isolated molecule with a specific conformation. Sampling conformations at specific 1S degrees angle values and comparing them against experimental C-13 chemical shift data for forms I and H identifies matching narrow ranges of conformations, successfully predicting the conformation of tolfenamic acid in each form. This methodology can therefore be used in crystal structure prediction to both reduce the initial conformational search space and also quantitatively assess subsequent putative structures to reliably and unambiguously identify the correct structure.