NMR studies of organic polymorphs and solvates

Spiers Memorial Lecture: NMR crystallography

Chemical function is directly related to the spatial arrangement of atoms. Consequently, the determination of atomic-level three-dimensional structures has transformed molecular and materials science over the past 60 years. In this context, solid-state NMR has emerged to become the method of choice for atomic-level characterization of complex materials in powder form. In the following we present an overview of current methods for chemical shift driven NMR crystallography, illustrated with applications to complex materials.

Exploring the crystallisation of aspirin in a confined porous material using solid-state nuclear magnetic resonance

In this study, nuclear magnetic resonance (NMR) is used to investigate the crystallisation behaviour of aspirin within a mesoporous SBA-15 silica material. The potential of dynamic nuclear polarisation (DNP) experiments is also investigated using specifically designed porous materials that incorporate polarising agents within their walls. The formation of the metastable crystalline form II is observed when crystallisation occurs within the pores of the mesoporous structure. Conversely, bulk crystallisation yields the most stable form, namely form I, of aspirin. Remarkably, the metastable form II remains trapped within the pores of mesoporous SBA-15 silica material even 30 days after impregnation, underscoring its persistent stability within this confined environment.

Tale of Two Polymorphs: Investigating the Structural Differences and Dynamic Relationship between Nirmatrelvir Solid Forms (Paxlovid)

Two anhydrous polymorphs of the novel antiviral medicine nirmatrelvir were discovered during the development of Paxlovid, Pfizer's oral Covid-19 treatment. A comprehensive experimental and computational approach was necessary to distinguish the two closely related polymorphs, herein identified as Forms 1 and 4. This approach paired experimental methods, including powder X-ray diffraction and single-crystal X-ray diffraction, solid-state experimental methods, thermal analysis, solid-state nuclear magnetic resonance and Raman spectroscopy with computational investigations comprising crystal structure prediction, Gibbs free energy calculations, and molecular dynamics simulations of the polymorphic transition. Forms 1 and 4 were ultimately determined to be enantiotropically related polymorphs with Form 1 being the stable form above the transition temperature of ∼17 °C and designated as the nominated form for drug development. The work described in this paper shows the importance of using highly specialized orthogonal approaches to elucidate the subtle differences in structure and properties of similar solid-state forms. This synergistic approach allowed for unprecedented speed in bringing Paxlovid to patients in record time amidst the pandemic.

An unusual ionic cocrystal of ponatinib hydrochloride: characterization by single-crystal X-ray diffraction and ultra-high field NMR spectroscopy

This study describes the discovery of a unique ionic cocrystal of the active pharmaceutical ingredient (API) ponatinib hydrochloride (pon·HCl), and characterization using single-crystal X-ray diffraction (SCXRD) and solid-state NMR (SSNMR) spectroscopy. Pon·HCl is a multicomponent crystal that features an unusual stoichiometry, with an asymmetric unit containing both monocations and dications of the ponatinib molecule, three water molecules, and three chloride ions. Structural features include (i) a charged imidazopyridazine moiety that forms a hydrogen bond between the ponatinib monocations and dications and (ii) a chloride ion that does not feature hydrogen bonds involving any organic moiety, instead being situated in a "square" arrangement with three water molecules. Multinuclear SSNMR, featuring high and ultra-high fields up to 35.2 T, provides the groundwork for structural interpretation of complex multicomponent crystals in the absence of diffraction data. A 13C CP/MAS spectrum confirms the presence of two crystallographically distinct ponatinib molecules, whereas 1D 1H and 2D 1H-1H DQ-SQ spectra identify and assign the unusually deshielded imidazopyridazine proton. 1D 35Cl spectra obtained at multiple fields confirm the presence of three distinct chloride ions, with density functional theory calculations providing key relationships between the SSNMR spectra and H⋯Cl- hydrogen bonding arrangements. A 2D 35Cl → 1H D-RINEPT spectrum confirms the spatial proximities between the chloride ions, water molecules, and amine moieties. This all suggests future application of multinuclear SSNMR at high and ultra-high fields to the study of complex API solid forms for which SCXRD data are unavailable, with potential application to heterogeneous mixtures or amorphous solid dispersions.

Exploiting solid-state dynamic nuclear polarization NMR spectroscopy to establish the spatial distribution of polymorphic phases in a solid material

Solid-state DNP NMR can enhance the ability to detect minor amounts of solid phases within heterogenous materials. Here we demonstrate that NMR contrast based on the transport of DNP-enhanced polarization can be exploited in the challenging case of early detection of a small amount of a minor polymorphic phase within a major polymorph, and we show that this approach can yield quantitative information on the spatial distribution of the two polymorphs. We focus on the detection of a minor amount (<4%) of polymorph III of m-aminobenzoic acid within a powder sample of polymorph I at natural isotopic abundance. Based on proposed models of the spatial distribution of the two polymorphs, simulations of 1H spin diffusion allow NMR data to be calculated for each model as a function of particle size and the relative amounts of the polymorphs. A comparison between simulated and experimental NMR data allows the model(s) best representing the spatial distribution of the polymorphs in the system to be established.

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.

Advanced Solid-State NMR Analysis of Two Crystal Forms of Ranitidine Hydrochloride: Detection of 1H-14N Intra-/Intermolecular Correlations

Understanding the characteristics of crystal polymorphism of active pharmaceutical ingredients and analyzing them with high sensitivity is important for quality of drug products, appropriate characterization strategies, and appropriate screening and selection processes. However, there are few methods to measure intra- and intermolecular correlations in crystals other than X-ray crystallography, for which it is sometimes difficult to obtain suitable single crystals. Recently, solid-state NMR has been recognized as a straightforward method for measuring molecular correlations. In this study, we selected ranitidine hydrochloride, which is known to exist in two forms, 1 and 2, as the model drug and investigated each form using solid-state NMR. In conducting the analysis, rotating the sample tube, which had a 1-mm inner diameter, increased the solid-state NMR resolution at 70 kHz. The ¹H-¹⁴N dipolar-based heteronuclear multiple quantum coherence (D-HMQC) analysis revealed the intermolecular correlation of Form 1 between the N atom of the nitro group and a proton of the furan moiety, which were closer than those of the intramolecular correlation reported using single X-ray crystal analysis. Thus, ¹H-¹⁴N D-HMQC analysis could be useful for characterizing intermolecular interaction in ranitidine hydrochloride crystals. In addition, we reassigned the ¹³C solid-state NMR signals of ranitidine hydrochloride according to the liquid-state and multiple solid-state NMR experiments.

Non-Covalent Reactions Supporting Antiviral Development

Viruses are the current big enemy of the world's healthcare systems. As the small infector causes various deadly diseases, from influenza and HIV to COVID-19, the virus continues to evolve from one type to its mutants. Therefore, the development of antivirals demands tremendous attention and resources for drug researchers around the world. Active pharmaceutical ingredients (API) development includes discovering new drug compounds and developing existing ones. However, to innovate a new antiviral takes a very long time to test its safety and effectiveness, from structure modeling to synthesis, and then requires various stages of clinical trials. Meanwhile, developing the existing API can be more efficient because it reduces many development stages. One approach in this effort is to modify the solid structures to improve their physicochemical properties and enhance their activity. This review discusses antiviral multicomponent systems under the research phase and has been marketed. The discussion includes the types of antivirals, their counterpart compound, screening, manufacturing methods, multicomponent systems yielded, characterization methods, physicochemical properties, and their effects on their pharmacological activities. It is hoped that the opportunities and challenges of solid antiviral drug modifications can be drawn in this review as important information for further antiviral development.

E-selective Semi-hydrogenation of Alkynes under Mild Conditions by a Diruthenium Hydride Complex

The synthesis, characterization and catalytic activity of a new class of diruthenium hydrido carbonyl complexes bound to the tBu PNNP expanded pincer ligand is described. Reacting tBu PNNP with two equiv of RuHCl(PPh3 )3 (CO) at 140 °C produces an insoluble air-stable complex, which was structurally characterized as [Ru2 (tBu PNNP)H(μ-H)Cl(μ-Cl)(CO)2 ] (1) using solid-state NMR, IR and X-ray absorption spectroscopies and follow-up reactivity. A reaction with KOtBu results in deprotonation of a methylene linker to produce [Ru2 (tBu PNNP* )H(μ-H)(μ-OtBu)(CO)2 ] (3) featuring a partially dearomatized naphthyridine core. This enables metal-ligand cooperative activation of H2 analogous to the mononuclear analogue, [Ru(tBu PNP*)H(CO)]. In contrast to the mononuclear system, the bimetallic analogue 3 catalyzes the E-selective semi-hydrogenation of alkynes at ambient temperature and atmospheric H2 pressure with good functional group tolerance. Monitoring the semi-hydrogenation of diphenylacetylene by 1 H NMR spectroscopy shows the intermediacy of Z-stilbene, which is subsequently isomerized to the E-isomer. Initial findings into the mode of action of this system are provided, including the spectroscopic characterization of a polyhydride intermediate and the isolation of a deactivated species with a partially hydrogenated naphthyridine backbone.

New salts of teriflunomide (TFM) - Single crystal X-ray and solid state NMR investigation

New salts of teriflunomide TFM (drug approved for Multiple Sclerosis treatment) with inorganic counterions: lithium (TFM_Li), sodium (TFM_Na), potassium (TFM_K), rubidium (TFM_Rb), caesium (TFM_Cs) and ammonium (TFM_NH₄) were prepared and investigated employing solid state NMR Spectroscopy, Powder X-ray Diffraction PXRD and Single Crystal X-ray Diffraction (SC XRD). Crystal and molecular structures of three salts: TFM_Na (CCDC: 2173257), TFM_Cs (CCDC: 2165288) and TFM_NH₄ (CCDC: 2165281) were determined and deposited. Compared to the native TFM, for all crystalline salt structures, a conformational change of the teriflunomide molecule involving about 180-degree rotation of the end group, forming an intramolecular hydrogen bond N-H⋯O is observed. By applying a complementary multi-technique approach, employing 1D and 2D solid state MAS NMR techniques, single and powder X-ray diffraction measurements, as well as the DFT-based GIPAW calculations of NMR chemical shifts for TFM_Na and TFM_Cs allowed to propose structural features of TFM_Li for which it was not possible to obtain adequate material for single crystal X-Ray measurement.

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.

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.

Multicomponent crystals of an artemisinin derivative and cinchona alkaloids for use as antimalarial drugs

Multicomponent crystals of an artemisinin derivative and cinchona alkaloids were produced, combining two major types of antimalaria drugs with unique hydrogen bond interactions. These salts demonstrate a new category of antimalarial pharmaceuticals.

Two-Photon Absorbing AIEgens: Influence of Stereoconfiguration on Their Crystallinity and Spectroscopic Properties and Applications in Bioimaging

This paper aims at designing chromophores with efficient aggregation-induced emission (AIE) properties for two-photon fluorescence microscopy (2PFM), which is one of the best-suited types of microscopy for the imaging of living organisms or thick biological tissues. Tetraphenylethylene (TPE) derivatives are common building blocks in the design of chromophores with efficient AIE properties. Therefore, in this study, extended TPE AIEgens specifically optimized for two-photon absorption (2PA) are synthesized and the resulting (E/Z) isomers are separated using chromatography on chiral supports. Comparative characterization of the AIE properties is performed on the pure (Z) and (E) isomers and the mixture, allowing us, in combination with powder X-ray diffraction and solid-state NMR, to document a profound impact of crystallinity on solid-state fluorescence properties. In particular, we show that stereopure AIEgens form aggregates of superior crystallinity, which in turn exhibit a higher fluorescence quantum yield compared to diastereoisomers mixtures. Preparation of stereopure organic nanoparticles affords very bright fluorescent contrast agents, which are then used for cellular and intravital two-photon microscopy on human breast cancer cells and on zebrafish embryos.

Selective detection of active pharmaceutical ingredients in tablet formulations using solid-state NMR spectroscopy

Atomic-level characterization of active pharmaceutical ingredients (API) is crucial in pharmaceutical industry because APIs play an important role in physicochemical properties of drug formulations. However, the analysis of targeted APIs in intact tablet formulations is less straightforward due to the coexistence of excipients as major components and different APIs at dilute concentrations (often below 10 ​wt% loading). Although solid-state (ss) NMR spectroscopy is widely used to investigate short-range order, polymorphism, and pseudo-polymorphism in neat pharmaceutical compounds, the analysis of complex drug formulations is often limited by overlapped signals that originate from structurally different APIs and excipients. In particular, such examples are frequently encountered in the analysis of ¹H ssNMR spectra of pharmaceutical formulations. While the high-resolution in ¹H ssNMR spectra can be attained by, for example, high magnetic fields accompanied by fast magic-angle spinning (MAS) approaches, the spectral complexity associated with the mixtures of compounds hinders the accurate determination of chemical shifts and through-space proximities. Here we propose a fast MAS (70 ​kHz) NMR experiment for the selective detection of ¹H signals associated with an API from a severely overlapped NMR spectrum of a tablet formulation. Spectral simplification is achieved by combining (i) symmetry-based dipolar recoupling (SR4₁²) rotational-echo saturation-pulse double-resonance (RESPDOR) with phase-modulate (PM) saturation pulses, (ii) radio frequency-driven recoupling (RFDR), and (iii) double-quantum excitation using Back-to-Back (BaBa) pulse sequence elements. First, ¹H sites in close proximities to ¹⁴N nuclei of an API are excited using a PM-S-RESPDOR sequence, and simultaneously, the other unwanted ¹H signals of excipients are suppressed. Then, ¹H magnetization transfer to adjacent ¹H sites in the API is achieved by spin diffusion process using a RFDR sequence, which polarizes to ¹H sites within the crystalline API regions of the drug formulation. Next, a PM-S-RESPDOR-RFDR sequence is combined with a Back-to-Back (BaBa) sequence to elucidate local-structures and ¹H-¹H proximities of the API in a dosage form. The PM-S-RESPDOR-RFDR-BaBa experiment is employed in one- (1D) and two-dimensional (2D) versions to selectively detect the ¹H ssNMR spectrum of l-cysteine (10.6 ​wt% or 0.11 ​mg) in a commercial formulation, and compared with the spectra of neat l-cysteine recorded using a standard BaBa experiment. The 2D ¹H double-quantum-single-quantum (DQ-SQ) spectrum of the API (l-cysteine)-detected pharmaceutical tablet is in good agreement with the 2D ¹H DQ-SQ spectrum obtained from the pure API molecule. Furthermore, the sensitivity and robustness of the experiment is examined by selectively detecting ¹H{¹⁴N} signals in an amino acid salt, l-histidine^.H₂O^.HCl.

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.

Genarris: Random generation of molecular crystal structures and fast screening with a Harris approximation

We present Genarris, a Python package that performs configuration space screening for molecular crystals of rigid molecules by random sampling with physical constraints. For fast energy evaluations, Genarris employs a Harris approximation, whereby the total density of a molecular crystal is constructed via superposition of single molecule densities. Dispersion-inclusive density functional theory is then used for the Harris density without performing a self-consistency cycle. Genarris uses machine learning for clustering, based on a relative coordinate descriptor developed specifically for molecular crystals, which is shown to be robust in identifying packing motif similarity. In addition to random structure generation, Genarris offers three workflows based on different sequences of successive clustering and selection steps: the "Rigorous" workflow is an exhaustive exploration of the potential energy landscape, the "Energy" workflow produces a set of low energy structures, and the "Diverse" workflow produces a maximally diverse set of structures. The latter is recommended for generating initial populations for genetic algorithms. Here, the implementation of Genarris is reported and its application is demonstrated for three test cases.

DNP-enhanced solid-state NMR spectroscopy of active pharmaceutical ingredients

Solid-state NMR spectroscopy has become a valuable tool for the characterization of both pure and formulated active pharmaceutical ingredients (APIs). However, NMR generally suffers from poor sensitivity that often restricts NMR experiments to nuclei with favorable properties, concentrated samples, and acquisition of one-dimensional (1D) NMR spectra. Here, we review how dynamic nuclear polarization (DNP) can be applied to routinely enhance the sensitivity of solid-state NMR experiments by one to two orders of magnitude for both pure and formulated APIs. Sample preparation protocols for relayed DNP experiments and experiments on directly doped APIs are detailed. Numerical spin diffusion models illustrate the dependence of relayed DNP enhancements on the relaxation properties and particle size of the solids and can be used for particle size determination when the other factors are known. We then describe the advanced solid-state NMR experiments that have been enabled by DNP and how they provide unique insight into the molecular and macroscopic structure of APIs. For example, with large sensitivity gains provided by DNP, natural isotopic abundance, ¹³ C-¹³ C double-quantum single-quantum homonuclear correlation NMR spectra of pure APIs can be routinely acquired. DNP also enables solid-state NMR experiments with unreceptive quadrupolar nuclei such as ² H, ¹⁴ N, and ³⁵ Cl that are commonly found in APIs. Applications of DNP-enhanced solid-state NMR spectroscopy for the molecular level characterization of low API load formulations such as commercial tablets and amorphous solid dispersions are described. Future perspectives for DNP-enhanced solid-state NMR experiments on APIs are briefly discussed.

Triboluminescence from Pharmaceutical Formulations

Triboluminescence (TL) is shown to enable selective detection of trace crystallinity within nominally amorphous solid dispersions (ASDs). ASDs are increasingly used for the preparation of pharmaceutical formulations, the physical stability of which can be negatively impacted by trace crystallinity introduced during manufacturing or storage. In the present study, TL measurements of a model ASD consisting of griseofulvin in polyethylene glycol produced limits of detection of 140 ppm. Separate studies of the particle size dependence of sucrose crystals and the dependence on polymorphism in clopidogrel bisulfate particles are both consistent with a mechanism for TL closely linked to the piezoelectric response of the crystalline fraction. Whereas disordered polymeric materials cannot support piezoelectric activity, molecular crystals produced from homochiral molecules adopt crystal structures that are overwhelmingly symmetry-allowed for piezoelectricity. Consequently, TL may provide a broadly applicable and simple experimental route for sensitive detection of trace crystallinity within nominally amorphous materials.

Does Z' equal 1 or 2? Enhanced powder NMR crystallography verification of a disordered room temperature crystal structure of a p38 inhibitor for chronic obstructive pulmonary disease

The crystal structure of the Form A polymorph of N-cyclopropyl-3-fluoro-4-methyl-5-[3-[[1-[2-[2-(methylamino)ethoxy]phenyl]cyclopropyl]amino]-2-oxo-pyrazin-1-yl]benzamide (i.e., AZD7624), determined using single-crystal X-ray diffraction (scXRD) at 100 K, contains two molecules in the asymmetric unit (Z' = 2) and has regions of local static disorder. This substance has been in phase IIa drug development trials for the treatment of chronic obstructive pulmonary disease, a disease which affects over 300 million people and contributes to nearly 3 million deaths annually. While attempting to verify the crystal structure using nuclear magnetic resonance crystallography (NMRX), we measured ¹³C solid-state NMR (SSNMR) spectra at 295 K that appeared consistent with Z' = 1 rather than Z' = 2. To understand this surprising observation, we used multinuclear SSNMR (¹H, ¹³C, ¹⁵N), gauge-including projector augmented-wave density functional theory (GIPAW DFT) calculations, crystal structure prediction (CSP), and powder XRD (pXRD) to determine the room temperature crystal structure. Due to the large size of AZD7624 (ca. 500 amu, 54 distinct ¹³C environments for Z' = 2), static disorder at 100 K, and (as we show) dynamic disorder at ambient temperatures, NMR spectral assignment was a challenge. We introduce a method to enhance confidence in NMR assignments by comparing experimental ¹³C isotropic chemical shifts against site-specific DFT-calculated shift distributions established using CSP-generated crystal structures. The assignment and room temperature NMRX structure determination process also included measurements of ¹³C shift tensors and the observation of residual dipolar coupling between ¹³C and ¹⁴N. CSP generated ca. 90 reasonable candidate structures (Z' = 1 and Z' = 2), which when coupled with GIPAW DFT results, room temperature pXRD, and the assigned SSNMR data, establish Z' = 2 at room temperature. We find that the polymorphic Form A of AZD7624 is maintained at room temperature, although dynamic disorder is present on the NMR timescale. Of the CSP-generated structures, 2 are found to be fully consistent with the SSNMR and pXRD data; within this pair, they are found to be structurally very similar (RMSD₁₆ = 0.30 Å). We establish that the CSP structure in best agreement with the NMR data possesses the highest degree of structural similarity with the scXRD-determined structure (RMSD₁₆ = 0.17 Å), and has the lowest DFT-calculated energy amongst all CSP-generated structures with Z' = 2.

Enhancing the resolution of 1H and 13C solid-state NMR spectra by reduction of anisotropic bulk magnetic susceptibility broadening

We demonstrate that natural isotopic abundance 2D heteronuclear correlation (HETCOR) solid-state NMR spectra can be used to significantly reduce or eliminate the broadening of ¹H and ¹³C solid-state NMR spectra of organic solids due to anisotropic bulk magnetic susceptibility (ABMS). ABMS often manifests in solids with aromatic groups, such as active pharmaceutical ingredients (APIs), and inhomogeneously broadens the NMR peaks of all nuclei in the sample. Inhomogeneous peaks with full widths at half maximum (FWHM) of ∼1 ppm typically result from ABMS broadening and the low spectral resolution impedes the analysis of solid-state NMR spectra. ABMS broadening of solid-state NMR spectra has previously been eliminated using 2D multiple-quantum correlation experiments, or by performing NMR experiments on diluted materials or single crystals. However, these experiments are often infeasible due to their poor sensitivity and/or provide limited gains in resolution. 2D ¹H-¹³C HETCOR experiments have previously been applied to reduce susceptibility broadening in paramagnetic solids and we show that this strategy can significantly reduce ABMS broadening in diamagnetic organic solids. Comparisons of 1D solid-state NMR spectra and ¹H and ¹³C solid-state NMR spectra obtained from 2D ¹H-¹³C HETCOR NMR spectra show that the HETCOR spectrum directly increases resolution by a factor of 1.5 to 8. The direct gain in resolution is determined by the ratio of the inhomogeneous ¹³C/¹H linewidth to the homogeneous ¹H linewidth, with the former depending on the magnitude of the ABMS broadening and the strength of the applied field and the latter on the efficiency of homonuclear decoupling. The direct gains in resolution obtained using the 2D HETCOR experiments are better than that obtained by dilution. For solids with long proton longitudinal relaxation times, dynamic nuclear polarization (DNP) was applied to enhance sensitivity and enable the acquisition of 2D ¹H-¹³C HETCOR NMR spectra. 2D ¹H-¹³C HETCOR experiments were applied to resolve and partially assign the NMR signals of the form I and form II polymorphs of aspirin in a sample containing both forms. These findings have important implications for ultra-high field NMR experiments, optimization of decoupling schemes and assessment of the fundamental limits on the resolution of solid-state NMR spectra.

Rotator side chains trigger cooperative transition for shape and function memory effect in organic semiconductors

Martensitic transition is a solid-state phase transition involving cooperative movement of atoms, mostly studied in metallurgy. The main characteristics are low transition barrier, ultrafast kinetics, and structural reversibility. They are rarely observed in molecular crystals, and hence the origin and mechanism are largely unexplored. Here we report the discovery of martensitic transition in single crystals of two different organic semiconductors. In situ microscopy, single-crystal X-ray diffraction, Raman and nuclear magnetic resonance spectroscopy, and molecular simulations combined indicate that the rotating bulky side chains trigger cooperative transition. Cooperativity enables shape memory effect in single crystals and function memory effect in thin film transistors. We establish a molecular design rule to trigger martensitic transition in organic semiconductors, showing promise for designing next-generation smart multifunctional materials.

Martensitic transition is commonly seen in steel and shape memory alloys but rarely in organic materials. Chung et al. discover martensitic transitions in organic electronics and utilize it in designing field-effect transistors, leading to shape memory effect that in return modifies charge transport properties.

Screening polymorphism in a Ni(ii) metal–organic framework: experimental observations, Hirshfeld surface analyses and DFT studies

A new polymorphic form of diaquaiminodiacetatonickel(ii) is structurally characterized by a detailed analysis of Hirshfeld surface and further analysed by DFT and QTAIM calculations.

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.

Reactions of Antitumor Active Dirhodium(II) Tetraacetate Rh2(CH3COO)4 with Cysteine and Its Derivatives

We have combined results from several spectroscopic techniques to investigate the aerobic reactions of Rh2(AcO)4 (AcO- = CH3COO-) with l-cysteine (H2Cys) and its derivatives d-penicillamine (3,3'-dimethylcysteine, H2Pen), with steric hindrance at the thiol group, and N-acetyl-l-cysteine (H2NAC), with its amino group blocked. Previous investigations have shown that antitumor active dirhodium(II) carboxylates may irreversibly inhibit enzymes containing a thiol group at or near their active sites. Also, cysteine, the only thiol-containing proteinogenic amino acid, interacts in vivo with this class of antitumor compounds, but structural information on the products of such reactions is lacking. In the present study, the reactions of Rh2(AcO)4 and H2L were carried out in aqueous solutions at the pH of mixing (acidic) and at physiological pH, using the different mole ratios 1:2, 1:4, and 1:6, which resulted in the same products in increasing yields. Electrospray ionization mass spectrometry (ESI-MS) indicates formation of dimeric [RhIII 2Pen4]2- or oligomeric {RhIII 2L4} n (L = Cys, NAC) complexes with bridging thiolate groups. Analyses of Rh K edge extended X-ray absorption fine structure (EXAFS) data reveal 3-4 Rh-S and 2-3 Rh-(N/O) bonds around six-coordinated Rh(III) ions at mean distances of 2.33 ± 0.02 and 2.09 ± 0.02 Å, respectively. In the N-acetyl-l-cysteine compound, the RhIII···RhIII distance 3.10 ± 0.02 Å obtained from the EXAFS spectrum supports trithiolate bridges between the Rh(III) ions, as was also found when using glutathione as ligand. In the cysteine and penicillamine complexes, double thiolate bridges join the Rh(III) ions, with the nonbridging Cys2- and Pen2- ligands in tridentate chelating (S,N,O) mode, which is consistent with the ΔδC = 7.3-8.4 ppm shift of the COO- signal in their carbon-13 cross polarization magic angle spinning (CPMAS) NMR spectra. For the penicillamine complex, the 2475.6 eV peak in its S K edge X-ray absorption near edge structure (XANES) spectrum shows partial oxidation, probably caused by peroxide generated from reduction of dissolved O2, of thiolato to sulfenato (S=O) groups, which were also identified by ESI-MS for all three {RhIII 2L4} n compounds.

Pharmaceutical cocrystals, salts and polymorphs: Advanced characterization techniques

The main goal of a novel drug development is to obtain it with optimal physiochemical, pharmaceutical and biological properties. Pharmaceutical companies and scientists modify active pharmaceutical ingredients (APIs), which often are cocrystals, salts or carefully selected polymorphs, to improve the properties of a parent drug. To find the best form of a drug, various advanced characterization methods should be used. In this review, we have described such analytical methods, dedicated to solid drug forms. Thus, diffraction, spectroscopic, thermal and also pharmaceutical characterization methods are discussed. They all are necessary to study a solid API in its intrinsic complexity from bulk down to the molecular level, gain information on its structure, properties, purity and possible transformations, and make the characterization efficient, comprehensive and complete. Furthermore, these methods can be used to monitor and investigate physical processes, involved in the drug development, in situ and in real time. The main aim of this paper is to gather information on the current advancements in the analytical methods and highlight their pharmaceutical relevance.

Exploring Systematic Discrepancies in DFT Calculations of Chlorine Nuclear Quadrupole Couplings

Previous studies have revealed significant discrepancies between density functional theory (DFT)-calculated and experimental nuclear quadrupolar coupling constants (C_Q) for chlorine atoms, particularly in ionic solids. Various aspects of the computations are systematically investigated here, including the choice of the DFT functional, basis set convergence, and geometry optimization protocol. The effects of fast (fs) time-scale dynamics are probed using molecular dynamics (MD) and nuclear quantum effects (NQEs) are considered using path-integral MD calculations. It is shown that the functional choice is the most important factor related to improving the accuracy of the quadrupolar coupling calculations, and that functionals beyond the generalized gradient approximation (GGA) level, such as hybrid and meta-GGA functionals, are required for good correlations with experiment. The influence of molecular dynamics and NQEs is less important than the functional choice in the studied systems. A method which involves scaling the calculated quadrupolar coupling constant is proposed here; its application leads to good agreement with experimental data.

Terahertz spectra of l-phenylalanine and its monohydrate

The low-frequency vibrational property of l-phenylalanine (l-Phe) and l-phenylalanine monohydrate (l-Phe·H₂O) has been investigated by terahertz time-domain spectroscopy (THz-TDS) at room and low temperature ranging from 0.5 to 4.5THz. Distinctive THz absorption spectra of the two compounds were observed. Density functional theory (DFT) calculations based on the crystal structures have been performed to simulate the vibrational modes of l-Phe and l-Phe·H₂O and the results agree well with the experimental observations. The study indicates that the characterized features of l-Phe mainly originate from the collective vibration of molecules. And the characterized features of l-Phe·H₂O mainly come from hydrogen bond interactions between l-Phe and water molecules. l-Phe and l-Phe·H₂O were also verified by differential scanning calorimetry and thermogravimetry (DSC-TG) and powder X-ray diffraction (PXRD) examinations.

Advanced solid-state NMR spectroscopy of natural organic matter

Solid-state NMR is essential for the characterization of natural organic matter (NOM) and is gaining importance in geosciences and environmental sciences. This review is intended to highlight advanced solid-state NMR techniques, especially a systematic approach to NOM characterization, and their applications to the study of NOM. We discuss some basics of how to acquire high-quality and quantitative solid-state ¹³C NMR spectra, and address some common technical mistakes that lead to unreliable spectra of NOM. The identification of specific functional groups in NOM, primarily based on ¹³C spectral-editing techniques, is described and the theoretical background of some recently-developed spectral-editing techniques is provided. Applications of solid-state NMR to investigating nitrogen (N) in NOM are described, focusing on limitations of the widely used ¹⁵N CP/MAS experiment and the potential of improved advanced NMR techniques for characterizing N forms in NOM. Then techniques used for identifying proximities, heterogeneities and domains are reviewed, and some examples provided. In addition, NMR techniques for studying segmental dynamics in NOM are reviewed. We also briefly discuss applications of solid-state NMR to NOM from various sources, including soil organic matter, aquatic organic matter, organic matter in atmospheric particulate matter, carbonaceous meteoritic organic matter, and fossil fuels. Finally, examples of NMR-based structural models and an outlook are provided.

Positional Variance in NMR Crystallography

We propose a method to quantify positional uncertainties in crystal structures determined by chemical-shift-based NMR crystallography. The method combines molecular dynamics simulations and density functional theory calculations with experimental and computational chemical shift uncertainties. In this manner we find the average positional accuracy as well as the isotropic and anisotropic positional accuracy associated with each atom in a crystal structure determined by NMR crystallography. The approach is demonstrated on the crystal structures of cocaine, flutamide, flufenamic acid, the K salt of penicillin G, and form 4 of the drug 4-[4-(2-adamantylcarbamoyl)-5-tert-butylpyrazol-1-yl]benzoic acid (AZD8329). We find that, for the crystal structure of cocaine, the uncertainty corresponds to a positional RMSD of 0.17 Å. This is a factor of 2.5 less than for single-crystal X-ray-diffraction-based structure determination.

`NMR Crystallization': in-situ NMR techniques for time-resolved monitoring of crystallization processes

Solid-state NMR spectroscopy is a well-established and versatile technique for studying the structural and dynamic properties of solids, and there is considerable potential to exploit the power and versatility of solid-state NMR for in-situ studies of chemical processes. However, a number of technical challenges are associated with adapting this technique for in-situ studies, depending on the process of interest. Recently, an in-situ solid-state NMR strategy for monitoring the evolution of crystallization processes has been developed and has proven to be a promising approach for identifying the sequence of distinct solid forms present as a function of time during crystallization from solution, and for the discovery of new polymorphs. The latest development of this technique, called `CLASSIC' NMR, allows the simultaneous measurement of both liquid-state and solid-state NMR spectra as a function of time, thus yielding complementary information on the evolution of both the liquid phase and the solid phase during crystallization from solution. This article gives an overview of the range of NMR strategies that are currently available for in-situ studies of crystallization processes, with examples of applications that highlight the potential of these strategies to deepen our understanding of crystallization phenomena.

Solid state NMR studies of gels derived from low molecular mass gelators

Since its invention more than six decades ago, nuclear magnetic resonance (NMR) spectroscopy has evolved as an inevitable part of chemical as well as structural analysis of small molecules, polymers, biomaterials and hybrid materials. In the solution state, due to the increased viscosity of complex viscoelastic fluids such as gels, liquid crystals and other soft materials, the rate of molecular tumbling is reduced, which in turn affects the chemical shift anisotropy, dipolar and quadrupolar interactions. As a consequence the solution state NMR spectra show broad lines, and therefore, extracting detailed structural information is a challenging task. In this context, solid state (SS) NMR has the ability to distinguish between a minute amount of polymorphic forms, conformational changes, and the number of non-equivalent molecules in an asymmetric unit of a crystal lattice, and to provide both qualitative as well as quantitative analytical data with a short-range order. Therefore, SS NMR has continued to evolve as an indispensable tool for structural analysis and gave birth to a new field called NMR crystallography. Solid state cross polarization (CP) and high resolution (HR) magic angle spinning (MAS) NMR spectroscopy has been used to study weak interactions in polymer gels. However, the application of SS NMR spectroscopy to study gels derived from low molecular weight gelators has been limited until recently. In this review, we will focus on the importance of solid state NMR spectroscopy in understanding and elucidating the structure of supramolecular gels derived from low molecular weight gelators with selected examples.

Structure elucidation of a complex CO2-based organic framework material by NMR crystallography

A three-dimensional structural model of a complex CO2-based organic framework made from high molecular weight, self-assembled, flexible and multi-functional oligomeric constituents has been determined de novo by solid-state NMR including DNP-enhanced experiments. The complete assignment of the 15N, 13C and 1H resonances was obtained from a series of two-dimensional through space and through bond correlation experiments. MM-QM calculations were used to generate different model structures for the material which were then evaluated by comparing multiple experimental and calculated NMR parameters. Both NMR and powder X-ray diffraction were evaluated as tools to determine the packing by crystal modelling, and at the level of structural modelling used here PXRD was found not to be a useful complement. The structure determined reveals a highly optimised H-bonding network that explains the unusual selectivity of the self-assembly process which generates the material. The NMR crystallography approach used here should be applicable for the structure determination of other complex solid materials.

Characterising the role of water in sildenafil citrate by NMR crystallography

A combination of solid-state NMR techniques, including ¹³C/¹H correlation, ²H magic-angle spinning NMR and first principles calculation are employed to characterise the role of water in different hydration states of sildenafil citrate.

Natural abundance 14N and 15N solid-state NMR of pharmaceuticals and their polymorphs

¹⁴N and ¹⁵N solid-state NMR at natural abundance are used in tandem for the investigation of pharmaceuticals and their polymorphs.

35Cl dynamic nuclear polarization solid-state NMR of active pharmaceutical ingredients

In this work, we show how to obtain efficient dynamic nuclear polarization (DNP) enhanced ³⁵Cl solid-state NMR (SSNMR) spectra at 9.4 T and demonstrate how they can be used to characterize the molecular-level structure of hydrochloride salts of active pharmaceutical ingredients (APIs) in both bulk and low wt% API dosage forms.