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.

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 35 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 13 C cross-polarization magic-angle spinning (CPMAS) NMR and supported by GIPAW-DFT calculations, and in the case of 3, 1 H solid-state NMR provided additional confirmation. Applying 35 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 1 H, 13 C, and 35 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.

Characterization of Amorphous Solid Dispersions and Identification of Low Levels of Crystallinity by Transmission Electron Microscopy.

Amorphous solid dispersions (ASDs) are used to increase the solubility of oral medicines by kinetically stabilizing the more soluble amorphous phase of an active pharmaceutical ingredient with a suitable amorphous polymer. Low levels of a crystalline material in an ASD can negatively impact the desired dissolution properties of the drug. Characterization techniques such as powder X-ray diffraction (pXRD), differential scanning calorimetry (DSC), and Fourier transform infrared spectroscopy (FTIR) are often used to detect and measure any crystallinity within ASDs. These techniques are unable to detect or quantify very low levels because they have limits of detection typically in the order of 1-5%. Herein, an ASD of felodipine (FEL) and polyvinylpyrrolidone/vinyl acetate copolymer (PVP/VA) prepared via a hot melt extrusion (HME) in a mass ratio of 30:70 was characterized using a range of techniques. No signs of residual crystallinity were found by pXRD, DSC, or FTIR. However, transmission electron microscopy (TEM) did identify two areas containing crystals at the edges of milled particles from a total of 55 examined. Both crystalline areas contained Cl Kα X-ray peaks when measured by energy-dispersive X-ray spectroscopy, confirming the presence of FEL (due to the presence of Cl atoms in FEL and not in PVP/VA). Further analysis was carried out by TEM using conical dark field (DF) imaging of a HME ASD of 50:50 FEL-PVP/VA to provide insights into the recrystallization process that occurs at the edges of particles during accelerated ageing conditions in an atmosphere of 75% relative humidity. Multiple metastable polymorphs of recrystallized FEL could be identified by selected area electron diffraction (SAED), predominately form II and the more stable form I. Conical DF imaging was also successful in spatially resolving and sizing crystals. This work highlights the potential for TEM-based techniques to improve the limit of detection of crystallinity in ASDs, while also providing insights into transformation pathways by identifying the location, size, and form of any crystallization that might occur on storage. This opens up the possibility of providing an enhanced understanding of a drug product's stability and performance.

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.

NMR and Thermal Studies for the Characterization of Mass Transport and Phase Separation in Paracetamol/Copovidone Hot-Melt Extrusion Formulations.

The formulation of drug/polymer amorphous solid dispersions (ASDs) is one of the most successful strategies for improving the oral bioavailability of poorly soluble active pharmaceutical ingredients (APIs). Hot-melt extrusion (HME) is one method for preparing ASDs that is growing in importance in the pharmaceutical industry, but there are still substantial gaps in our understanding regarding the dynamics of drug dissolution and dispersion in viscous polymers and the physical stability of the final formulations. Furthermore, computational models have been built to predict optimal processing conditions, but they are limited by the lack of experimental data for key mass transport parameters, such as the diffusion coefficient. The work presented here reports direct measurements of API diffusion in pharmaceutical polymer melts, using high-temperature pulsed-field gradient NMR. The diffusion coefficient of a model drug/polymer system (paracetamol/copovidone) was determined for different drug loadings and at temperatures relevant to the HME process. The mechanisms of the diffusion process are also explored with the Stokes-Einstein and Arrhenius models. The results show that diffusivity is linked exponentially to temperature. Furthermore, this study includes rheological characterization, differential scanning calorimetry (DSC), and 1H ssNMR T1 and T1ρ measurements to give additional insights into the physical state, phase separation, and API/polymer interactions in paracetamol/copovidone ASD formulations.

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 (ΔδExperimental) is calculated as the difference between 1H and 13C experimental solid-state chemical shifts for each polymorphic form (δSolid expt) and the corresponding solution-state NMR chemical shifts (δSolution expt). Separately, we use the gauge-included projector augmented wave (GIPAW) method to calculate the NMR chemical shifts for each form (δSolid calc) and for TFA in solution (δ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 (ΔδCalculated) is then calculated as the difference of δSolid calc and δSolution calc. Regression analysis for ΔδCalculated against Δδ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 ΔδExperimental based on the experimental MAS NMR chemical shifts of form I or II with ΔδCalculated based on calculated chemical shifts for polymorphs I, II, III, and IV. Complementarity to the established approach of comparing δSolid expt to δSolid calc is explored. We further show that our approach is applicable if there are no solid-state crystal structure data. Specifically, δSolid calc in ΔδCalculated is replaced by the chemical shift for an isolated molecule with a specific conformation. Sampling conformations at specific 15° angle values and comparing them against experimental 13C chemical shift data for forms I and II 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.

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 1H and 13C cross-polarisation (CP) magic-angle spinning (MAS) nuclear magnetic resonance (NMR) and two-dimensional 14N-1H heteronuclear multiple-quantum coherence MAS NMR spectra are compared with gauge-including projector augmented wave (GIPAW) calculations of the 1H and 13C chemical shifts and the 14N shifts that additionally depend on the quadrupolar interaction. Considering the high ppm (>10 â€‹ppm) 1H 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 13C 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 13C 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.

5-amino-2-methylpyridinium hydrogen fumarate: An XRD and NMR crystallography analysis.

Single-crystal X-ray diffraction structures of the 5-amino-2-methylpyridinium hydrogen fumarate salt have been solved at 150 and 300 K (CCDC 1952142 and 1952143). A base-acid-base-acid ring is formed through pyridinium-carboxylate and amine-carboxylate hydrogen bonds that hold together chains formed from hydrogen-bonded hydrogen fumarate ions. 1 H and 13 C chemical shifts as well as 14 N shifts that additionally depend on the quadrupolar interaction are determined by experimental magic angle spinning (MAS) solid-state nuclear magnetic resonance (NMR) and gauge-including projector-augmented wave (GIPAW) calculation. Two-dimensional homonuclear 1 H-1 H double-quantum (DQ) MAS and heteronuclear 1 H-13 C and 14 N-1 H spectra are presented. Only small differences of up to 0.1 and 0.6 ppm for 1 H and 13 C are observed between GIPAW calculations starting with the two structures solved at 150 and 300 K (after geometry optimisation of atomic positions, but not unit cell parameters). A comparison of GIPAW-calculated 1 H chemical shifts for isolated molecules and the full crystal structures is indicative of hydrogen bonding strength.

Toward Developing a Predictive Approach To Assess Electron Beam Instability during Transmission Electron Microscopy of Drug Molecules.

During drug development control of polymorphism, particle properties and impurities are critical for ensuring a good quality, reproducible, and safe medicine. A wide variety of analytical techniques are employed in demonstrating the regulators control over the drug substance and product manufacturing, storage, and supply. Transmission electron microscopy (TEM) offers the opportunity to analyze in detail pharmaceutical systems at a length scale and limit of detection not readily achieved by many traditional techniques. However, the use of TEM as a characterization tool for drug development is uncommon due to possible damage caused by the electron beam. This work outlines the development of a model, using molecular descriptors, to predict the electron beam stability of active pharmaceutical ingredients (API). For a given set of conditions and a particular imaging or analytical mode, the total number of electrons per unit area, which causes observable damage to a sample in the TEM, can be defined as the critical fluence ( CF). Here the CF of 20 poorly water-soluble APIs were measured using selected area electron diffraction. Principal component analysis was used to select the most influential molecular descriptors on CF, which were shown to be descriptors involving the degree of conjugation, the number of hydrogen bond donors and acceptors, and the number of rotatable bonds. These were used to generate several multiple linear regression models. The model that provided the best fit to the measured CF data included the ratio of the number of conjugated carbons to nonconjugated carbons, the ratio of the number of hydrogen bond donors to acceptors, and the ratio of the number of hydrogen bond acceptors to donors. Using this model, the CF of the majority of the compounds was predicted within ±2 e-/Å2. Molecules with no hydrogen bond acceptors did not fit the model accurately possibly due to the limited sample size or the influence of other parameters not included in this model, such as intermolecular bond energies. The model presented can be used to support pharmaceutical development by quickly assessing the stability of other poorly soluble drugs in TEM. Provided that the model suggests that the API is relatively stable to electron irradiation, TEM offers the prospect of determining the presence of crystalline material at low levels at length scales and limits of detection unobtainable by other techniques. This is particularly so for amorphous solid dispersions.

Indomethacin-Kollidon VA64 Extrudates: A Mechanistic Study of pH-Dependent Controlled Release.

Because of its weakly acidic nature (pKa of 4.5), indomethacin presents an aqueous solubility that significantly increases when changing from acidic to neutral/alkaline pH (1.5 µg/mL at pH 1.2 and 105.2 µg/mL at pH 7.4). We have therefore investigated the impact of the dissolution medium pH on the dissolution performance of indomethacin:Kollidon VA64 extrudates. The impact of the drug loading on the dissolution properties of these systems was also examined (5%, 15%, 30%, 50%, 70%, and 90% drug loading). Time-resolved Raman spectroscopy along with in-line UV-vis spectrophotometry was employed to directly relate changes in dissolution behavior to physicochemical changes that occur to the extrudate during the test. The dissolution tests were performed in pH 2 HCl (to mimic the stomach conditions), and this was then switched during the experiment to pH 6.8 phosphate buffer (to simulate the poststomach conditions). The rotating disc dissolution rate test was also used to simultaneously measure the dissolution rate of both the drug and the polymer. We found that in pH 2 HCl buffer, for the 15% or higher drug-loaded extrudates, Kollidon VA64 preferentially dissolves from the exterior of the compact leaving an amorphous drug-rich hydrophobic shell, which, similarly to an enteric coating, inhibits the drug release. The in situ formation of an enteric coating has been previously hypothesized, and this has been the first time that is directly observed in a pH-variable dissolution test. The dissolution medium switch to pH 6.8 phosphate buffer, due to the large increase of the aqueous solubility of indomethacin at this pH, leads to rapid dissolution of the material forming the coating and therefore total drug release. In contrast, the 5% extrudate is fully hydrated and quickly dissolves at low pH pointing to a dissolution performance dependent on highly water-soluble Kollidon VA64.

Monitoring the Dissolution Mechanisms of Amorphous Bicalutamide Solid Dispersions via Real-Time Raman Mapping.

Real-time in situ Raman mapping has been employed to monitor, during dissolution, the crystallization transitions of amorphous bicalutamide formulated as a molecular dispersion in a copovidone VA64 matrix. The dissolution performance was also investigated using the rotating disc dissolution rate methodology, which allows simultaneous determination of the dissolution rate of both active ingredient and polymer. The dissolution behavior of two bicalutamide:copovidone VA64 dispersion formulations, containing 5% (w/w) and 50% (w/w) bicalutamide, respectively, was investigated, with the aim of exploring the effect of increasing the bicalutamide loading on the dissolution performance. Spatially time-resolved Raman maps generated using multivariate curve resolution indicated the simultaneous transformation of amorphous bicalutamide present in the 50% drug-loaded extrudate into metastable polymorphic form II and low-energy polymorphic form I. Fitting a kinetic model and spatially correlating the data extracted from the Raman maps also allowed us to understand the re-crystallization mechanisms by which the low-energy form I appears. Form I was shown to crystallize mainly directly from the amorphous solid dispersion, with crystallization from the metastable form II being a minor contribution.

Investigating the Dissolution Performance of Amorphous Solid Dispersions Using Magnetic Resonance Imaging and Proton NMR.

We have investigated the dissolution performance of amorphous solid dispersions of poorly water-soluble bicalutamide in a Kollidon VA64 polymeric matrix as a function of the drug loading (5% vs. 30% bicalutamide). A combined suite of state-of-the-art analytical techniques were employed to obtain a clear picture of the drug release, including an integrated magnetic resonance imaging UV-Vis flow cell system and 1H-NMR. Off-line 1H-NMR was used for the first time to simultaneously measure the dissolution profiles and rates of both the drug and the polymer from a solid dispersion. MRI and 1H-NMR data showed that the 5% drug loading compact erodes linearly, and that bicalutamide and Kollidon VA64 are released at approximately the same rate from the molecular dispersion. For the 30% extrudate, data indicated a slower water ingress into the compact which corresponds to a slower dissolution rate of both bicalutamide and Kollidon VA64.

Proton NMR: a new tool for understanding dissolution.

We present the use of (1)H NMR as a new measurement approach for improving understanding of the dissolution of pharmaceutical tablets. NMR has benefits over the conventional UV measurement approach in respect to much greater analyte selectivity and the ability to detect non-UV-absorbing species such as sugars. We used an in-line flow cell and water suppression experiments to determine the release profiles of three drug substances and lactose from the same tablet. Dissolution was performed in a pharmacopieal dissolution system with a standard protic buffer. NMR was shown to give high selectivity with each analyte having a well-resolved signal and sufficient sensitivity to determine the full release profile of even a compound present at only 5 mg in the tablet. The in-line flow cell gives excellent quality NMR spectra having little impact on peak shape. Dissolution of all the drug substances and lactose were determined to proceed at the same relative rates.

The Global Characterisation of a Drug-Dendrimer Conjugate - PEGylated poly-lysine Dendrimer.

The recent emergence of drug-dendrimer conjugates within pharmaceutical industry research and development introduces a range of challenges for analytical and measurement science. These molecules are very high molecular weight (100-200kDa) with a significant degree of structural complexity. The characteristics and quality attributes that require understanding and definition, and impact efficacy and safety, are diverse. They relate to the intact conjugate, the various building blocks of these complex systems and the level of the free and bound active pharmaceutical ingredient (API). From an analytical and measurement science perspective, this necessitates the measurement of the molecular weight, impurity characterisation, the quantitation of the number of conjugated versus free API molecules, the determination of the impurity profiles of the building blocks, primary structure and both particle size and morphology. Here we report the first example of a global characterisation of a drug-dendrimer conjugate - PEGylated poly-lysine dendrimer currently under development (AZD0466). The impact of the wide variety of analytical and measurement techniques on the overall understanding of this complex molecular entity is discussed, with the relative capabilities of the various approaches compared. The results of this study are an essential platform for the research and development of the future generations of related dendrimer-based medicines.

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 13C 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 13C 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.

Measuring the Rate of In-vitro Drug Release From Polymeric Nanoparticles by 19F Solution State NMR Spectroscopy.

We report what we believe is the first use of 19F NMR spectroscopy to directly measure in-vitro release (IVR) from polymeric nanoparticles (PNPs). Using 19F NMR we selectively measured IVR of AZD2811 from PNPs. Due to rapid nuclear relaxation in solid-like environments only AZD2811 in solution is detected, and physical separation from the PNPs isn't required. The NMR approach and ultra-centrifugation/UHPLC were shown to be equivalent. The selectivity of 19F NMR means it is readily applied to complex IVR media such as recombinant human serum albumin (rHSA).

A curious case of dynamic disorder in pyrrolidine rings elucidated by NMR crystallography.

A pharmaceutical exhibits differing dynamics in crystallographically distinct pyrrolidine rings despite being nearly related by symmetry, with one performing ring inversions while the other is constrained to torsional librations. Using 13C solid-state magic-angle spinning (MAS) NMR and DFT calculations, we show that this contrast originates from C-HH-C close contacts and less efficient C-Hπ intermolecular interactions observed in the transition state of the constrained pyrrolidine ring, highlighting the influence of the crystallographic environment on the molecular motion.

Structural characterisation of amorphous solid dispersions via metropolis matrix factorisation of pair distribution function data.

We measure the X-ray pair distribution functions (PDFs) of a series of felodipine:copovidone amorphous solid dispersions. Using a newly-developed Metropolis Matrix Factorisation (MMF) algorithm we extract from these data the PDF of the amorphous felodipine component in isolation. Our MMF analysis allows quantification of the degree of drug crystallinity in each sample, and structural characterisation of the amorphous drug via its PDF. Comparison with atomistic simulations reveals that the (in)accessibility of conformational rotamers distinguishes amorphous and crystalline felodipine, in turn suggesting design routes for stabilising the amorphous form. We discuss the conceptual importance of our results in the context of characterising not only amorphous pharmaceuticals, but complex mixtures in general.

Probing intermolecular hydrogen bonding in sibenadet hydrochloride polymorphs by high-resolution (1) H double-quantum solid-state NMR spectroscopy.

Molecular packing in two polymorphs of sibenadet hydrochloride (AR-C68397AA, Viozan™) is investigated using a combined experimental (1) H double-quantum (DQ) solid-state magic-angle spinning nuclear magnetic resonance and computational (gauge including projected augmented wave calculation of chemical shifts) approach. For Form I, NH-NH and NH-OH (1) H DQ peaks are observed corresponding to nearest distances of 2.62 and 2.87 Å, respectively, for the intermolecular hydrogen-bonding arrangement in the single-crystal X-ray diffraction structure. The same (1) H DQ peaks at the same (1) H chemical shifts are observed for Form II, for which there is no single-crystal diffraction structure, indicating the same intermolecular hydrogen-bonding arrangement of the benzothiazolone moieties as in Form I. (1) H DQ build-up (as a function of the DQ recoupling time) curves are presented for the resolved NH-NH and NH-OH DQ peaks for the two polymorphs. For Form I, the ratio of the maximum intensity for the NH-OH and NH-NH DQ peaks is in excellent agreement with the ratio of the summed squares of the H-H dipolar couplings, as determined using H-H distances from the crystal structure up to 4 Å. Small differences in the (1) H DQ build-up behaviour for the two polymorphs are attributed to differences in the longer-range NH-OH distances associated with different inter-layer arrangements.