Narrowing down the conformational space with solid-state NMR in crystal structure prediction of linezolid cocrystals

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.

Atomic-level structure determination of amorphous molecular solids by NMR

Structure determination of amorphous materials remains challenging, owing to the disorder inherent to these materials. Nuclear magnetic resonance (NMR) powder crystallography is a powerful method to determine the structure of molecular solids, but disorder leads to a high degree of overlap between measured signals, and prevents the unambiguous identification of a single modeled periodic structure as representative of the whole material. Here, we determine the atomic-level ensemble structure of the amorphous form of the drug AZD4625 by combining solid-state NMR experiments with molecular dynamics (MD) simulations and machine-learned chemical shifts. By considering the combined shifts of all 1H and 13C atomic sites in the molecule, we determine the structure of the amorphous form by identifying an ensemble of local molecular environments that are in agreement with experiment. We then extract and analyze preferred conformations and intermolecular interactions in the amorphous sample in terms of the stabilization of the amorphous form of the drug.

On the Use of Deuterated Organic Solvents without TMS to Report 1H/13C NMR Spectral Data of Organic Compounds: Current State of the Method, Its Pitfalls and Benefits, and Related Issues

The quite popular, simple but imperfect method of referencing NMR spectra to residual ¹H and ¹³C signals of TMS-free deuterated organic solvents (referred to as Method A) is critically discussed for six commonly used NMR solvents with respect to their δ_H and δ_C values that exist in the literature. Taking into account the most reliable data, it was possible to recommend 'best' δ_X values for such secondary internal standards. The position of these reference points on the δ scale strongly depends on the concentration and type of analyte under study and the solvent medium used. For some solvents, chemically induced shifts (CISs) of residual ¹H lines were considered, also taking into account the formation of 1:1 molecular complexes (for CDCl₃). Typical potential errors that can occur as a result of improper application of Method A are considered in detail. An overview of all found δ_X values adopted by users of this method revealed a discrepancy of up to 1.9 ppm in δ_C reported for CDCl₃, most likely caused by the CIS mentioned above. The drawbacks of Method A are discussed in relation to the classical use of an internal standard (Method B), two 'instrumental' schemes in which Method A is often implicitly applied, that is, the default Method C using ²H lock frequencies and Method D based on Ξ values, recommended by the IUPAC but only occasionally used for ¹H/¹³C spectra, and external referencing (Method E). Analysis of current needs and opportunities for NMR spectrometers led to the conclusion that, for the most accurate application of Method A, it is necessary to (a) use dilute solutions in a single NMR solvent and (b) to report δ_X data applied for the reference ¹H/¹³C signals to the nearest 0.001/0.01 ppm to ensure the precise characterization of new synthesized or isolated organic systems, especially those with complex or unexpected structures. However, the use of TMS in Method B is strongly recommended in all such cases.

Zwitterionic or Not? Fast and Reliable Structure Determination by Combining Crystal Structure Prediction and Solid-State NMR

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., ¹H Magic-Angle Spinning (MAS), ¹³C and ¹⁵N Cross Polarisation Magic-Angle Spinning (CPMAS), ¹H Double Quantum (DQ) MAS, ¹H-¹³C 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 ¹H and ¹³C 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., ¹⁴N-¹H 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.