Enforcing Ostwald's rule of stages: Isolation of paracetamol forms III and II

Polymorphic transitions in flufenamic acid-trehalose composites

The combination of poorly-soluble drugs with small molecule co-formers to generate amorphous solid dispersions (ASDs) has great potential to improve dissolution rate and kinetic solubility, and thus increase the bioavailability of these active ingredients. However, such ASDs are known to be unstable and to crystallise upon storage or heating. In this work, we explore the crystallisation of flufenamic acid (FFA) from ASDs prepared with trehalose. FFA-trehalose mixtures were prepared at a range of w/w composition ratios, heated to melting and crash cooled to form ASDs. They were then subject to a further heat/cool cycle, which was monitored by simultaneous differential scanning calorimetry - X-ray diffraction to observe the phase changes occurring. These varied with the composition of the blend. Upon short-term storage, formulations with low trehalose contents (FFA:trehalose 5:1 w/w) recrystallised into form I FFA, while higher trehalose contents crystallised to FFA form IV. When heated, all FFA trehalose combinations ultimately recrystallised into form I before melting. Upon a second cooling cycle, systems with low trehalose content (FFA:trehalose 5:1 w/w) recrystallised into form IV, while higher trehalose contents led to FFA form I. It is thus clear that even with a single excipient it is possible to control the crystallisation pathway through judicious choice of the formulation parameters.

Exploring the role of crystal habit in the Ostwald rule of stages

The crystallization of calcium carbonate is shown to be dictated by the Ostwald rule of stages (ORS), for high relative initial supersaturations (S CaCO 3 = [ C a 2 + ] [ CO 3 2 - ] / K SP, Calcite > 2500 ), under sweet (carbon dioxide saturated) and anoxic (oxygen depleted) solution conditions. Rhombohedral calcite crystals emerge after the sequential crystallization and dissolution of the metastable polymorphs: vaterite (snowflake-shaped) and aragonite (needle-shaped). However, the presence of certain cations, which can form trigonal carbonates (e.g. Fe2+ and Ni2+), in concentrations as low as 1.5 mM, triggers the emergence of calcite crystals, with a star-shaped crystal habit, first. These star-shaped crystals dissolve to yield needle-shaped aragonite crystals, which in turn dissolve to give the rhombohedral calcite crystals. The star-shaped crystals, formed at high SCaCO3 , possess higher surface free energy (therefore higher apparent solubility) than their rhombohedral counterparts. This sequence of dissolution and recrystallization demonstrates that the ORS does not only drive the crystal towards its thermodynamically most stable polymorph but also towards its most stable crystal habit.

Drug-Polymer Interactions in Acetaminophen/Hydroxypropylmethylcellulose Acetyl Succinate Amorphous Solid Dispersions Revealed by Multidimensional Multinuclear Solid-State NMR Spectroscopy

The bioavailability of insoluble crystalline active pharmaceutical ingredients (APIs) can be enhanced by formulation as amorphous solid dispersions (ASDs). One of the key factors of ASD stabilization is the formation of drug-polymer interactions at the molecular level. Here, we used a range of multidimensional and multinuclear nuclear magnetic resonance (NMR) experiments to identify these interactions in amorphous acetaminophen (paracetamol)/hydroxypropylmethylcellulose acetyl succinate (HPMC-AS) ASDs at various drug loadings. At low drug loading (<20 wt %), we showed that 1H-13C through-space heteronuclear correlation experiments identify proximity between aromatic protons in acetaminophen with cellulose backbone protons in HPMC-AS. We also show that 14N-1H heteronuclear multiple quantum coherence (HMQC) experiments are a powerful approach in probing spatial interactions in amorphous materials and establish the presence of hydrogen bonds (H-bond) between the amide nitrogen of acetaminophen with the cellulose ring methyl protons in these ASDs. In contrast, at higher drug loading (40 wt %), no acetaminophen/HPMC-AS spatial proximity was identified and domains of recrystallization of amorphous acetaminophen into its crystalline form I, the most thermodynamically stable polymorph, and form II are identified. These results provide atomic scale understanding of the interactions in the acetaminophen/HPMC-AS ASD occurring via H-bond interactions.

Τhe role of surface energy in the apparent solubility of two different calcite crystal habits

The interplay between polymorphism and facet-specific surface energy on the dissolution of crystals is examined in this work. It is shown that, using cationic additives, it is possible to produce star-shaped calcite crystals at very high supersaturations. In crystallization processes following the Ostwald rule of stages these star-shaped crystals appear to have higher solubility than both their rhombohedral counterparts and needle-shaped aragonite crystals. The vapour pressures of vaterite, aragonite, star-shaped calcite and rhombohedral calcite crystals are measured using thermogravimetric analysis and the corresponding enthalpies of melting are obtained. Using inverse gas chromatography, the surface energy of the aforementioned crystals is measured as well and the surface energy of the main crystal facets is calculated. Combining the effect of facet-specific surface energies and the enthalpies of melting on a modified version of the classical solubility equation for regular solutions, it is proved that the star-shaped calcite crystals can indeed have higher apparent solubility than aragonitecrystals.

Mechanochemistry: A Green Approach in the Preparation of Pharmaceutical Cocrystals

Mechanochemistry is considered an alternative attractive greener approach to prepare diverse molecular compounds and has become an important synthetic tool in different fields (e.g., physics, chemistry, and material science) since is considered an ecofriendly procedure that can be carried out under solvent free conditions or in the presence of minimal quantities of solvent (catalytic amounts). Being able to substitute, in many cases, classical solution reactions often requiring significant amounts of solvents. These sustainable methods have had an enormous impact on a great variety of chemistry fields, including catalysis, organic synthesis, metal complexes formation, preparation of multicomponent pharmaceutical solid forms, etc. In this sense, we are interested in highlighting the advantages of mechanochemical methods on the obtaining of pharmaceutical cocrystals. Hence, in this review, we describe and discuss the relevance of mechanochemical procedures in the formation of multicomponent solid forms focusing on pharmaceutical cocrystals. Additionally, at the end of this paper, we collect a chronological survey of the most representative scientific papers reporting the mechanochemical synthesis of cocrystals.

Selective formation of form II paracetamol through the assistance of paracetamol co-crystals as templates in a solution

The arrangement of paracetamol templated by the structure of 1 : 1 co-crystal of paracetamol–maleic acid in solution phase facilitating the nucleation of metastable form II paracetamol.

Challenging the Ostwald rule of stages in mechanochemical cocrystallisation

Mechanochemistry provides an efficient, but still poorly understood route to synthesize and screen for polymorphs of organic solids. We present a hitherto unexplored effect of the milling assembly on the polymorphic outcome of mechanochemical cocrystallisation, tentatively related to the efficiency of mechanical energy transfer to the milled sample. Previous work on mechanochemical cocrystallisation has established that introducing liquid or polymer additives to milling systems can be used to direct polymorphic behavior, leading to extensive studies how the amount and nature of grinding additive affect reaction outcome and polymorphism. Here, focusing on a model pharmaceutical cocrystal of nicotinamide and adipic acid, we demonstrate that changes to the choice of milling media (i.e. number and material of milling balls) and/or the choice of milling assembly (i.e. jar material) can be used to direct polymorphism of mechanochemical cocrystallisation, enabling the selective synthesis, and even reversible and repeatable interconversion of cocrystal polymorphs. While real-time mechanistic studies of mechanochemical transformations of metal–organic materials have previously suggested that reactions follow a path described by Ostwald's rule of stages, i.e. from metastable to increasingly more stable product structures, the herein presented systematic study presents an exception to that rule, revealing that modification of energy input in the mechanochemical system, combined with a small energy difference between polymorphs, permits the selective synthesis of either the more stable room temperature form, or the new metastable high-temperature form, of the target cocrystal.

The choice of milling assembly (jar and ball material, number and size of balls) can be used to direct polymorphism in mechanochemical cocrystallisation, enabling the selective synthesis, and even reversible interconversion of cocrystal polymorphs.

Mechanochemical Metathesis between AgNO3 and NaX (X = Cl, Br, I) and Ag2XNO3 Double-Salt Formation

Here we describe real-time, in situ monitoring of mechanochemical solid-state metathesis between silver nitrate and the entire series of sodium halides, on the basis of tandem powder X-ray diffraction and Raman spectroscopy monitoring. The mechanistic monitoring reveals that reactions of AgNO₃ with NaX (X = Cl, Br, I) differ in reaction paths, with only the reaction with NaBr providing the NaNO₃ and AgX products directly. The reaction with NaI revealed the presence of a novel, short-lived intermediate phase, while the reaction with NaCl progressed the slowest through the well-defined Ag₂ClNO₃ intermediate double salt. While the corresponding iodide and bromide double salts were not observed as intermediates, all three are readily prepared as pure compounds by milling equimolar mixtures of AgX and AgNO₃. The in situ observation of reactive intermediates in these simple metathesis reactions reveals a surprising resemblance of reactions involving purely ionic components to those of molecular organic solids and cocrystals. This study demonstrates the potential of in situ reaction monitoring for mechanochemical reactions of ionic compounds as well as completes the application of these techniques to all major compound classes.

Mechanochemistry: an efficient and versatile toolbox for synthesis, transformation, and functionalization of porous metal–organic frameworks

Multiple ways in which the synergy of mechanochemistry and MOFs advances the field of materials chemistry are presented here.

Improving Predicted Nuclear Magnetic Resonance Chemical Shifts Using the Quasi-Harmonic Approximation

Ab initio nuclear magnetic resonance chemical shift prediction plays an important role in the determination or validation of crystal structures. The ability to predict chemical shifts more accurately can translate to increased confidence in the resulting chemical shift or structural assignments. Standard electronic structure predictions for molecular crystal structures neglect thermal expansion, which can lead to an appreciable underestimation of the molar volumes. This study examines this volume error and its impact on 68 ¹³C- and 28 ¹⁵N-predicted chemical shifts taken from 20 molecular crystals. It assesses the ability to recover more realistic room-temperature crystal structures using the quasi-harmonic approximation and how refining the structures impacts the chemical shifts. Several pharmaceutical molecular crystals are also examined in more detail. On the whole, accounting for quasi-harmonic expansion changes the ¹³C and ¹⁵N chemical shifts by 0.5 and 1.0 ppm on average. This, in turn, reduces the root-mean-square errors relative to experiment by 0.3 ppm for ¹³C and 0.7 ppm for ¹⁵N. Although the statistical impacts are modest, changes in individual chemical shifts can reach multiple ppm. Accounting for thermal expansion in molecular crystal chemical shift prediction may not be needed routinely, but the systematic trend toward improved accuracy with the experiment could be useful in cases where discrimination between structural candidates is challenging, as in the pharmaceutical theophylline.

Continuous-flow crystallisation in 3D-printed compact devices

A flexible and cost-effective methodology to develop compact flow devices with heat exchange ability is presented here.

Impact of Additives on Heterogeneous Crystallization of Acetaminophen

Introducing foreign substrates or additives is the common way to regulate polymorphism or kinetics of crystallization. Most present studies consider the substrate factor and additive factor separately. Here, the interplay between the additive, crystallizing molecules, and the substrate was investigated. Acetaminophen (APAP) was used as the model compound. 5 wt.% dioctyl sodium sulfosuccinate (AOT), poly(acrylic) acid (PAA), and hydroxypropyl methylcellulose (HPMC) were employed as additives. The interfacial crystal growth rate of APAP in the presence of additives was studied between slides coated with chitosan (CS) film. The crystallization kinetics of the additive/APAP mixture on CS substrate was also investigated. The additive/APAP was characterized by differential scanning calorimetry (DSC), and the interfacial molecular interaction was studied by Fourier transform infrared spectroscopy (FTIR). The results indicate that the additive-substrate interaction can change the interfacial growth behavior observed in the additive-APAP binary system. Nevertheless, crystallizing without confinement, the additive-APAP interaction is more effective at controlling the crystallization of APAP, and the substrates did not have much effect.

Bulk and Surface-Stabilized Structures of Paracetamol Revisited by Raman Confocal Microscopy

We revisit the polymorphism of paracetamol by means of a micro-Raman technique, which has proved to be a powerful tool for structure recognition. Distinct lattice phonon spectra clearly identified the pure phases. Confocality enabled us to detect phase mixing between form II and either I or III on a micrometric scale in the same crystallite. Following the most recent findings on surface-mediated structures, we also investigated spin-coated films grown on glass, gold, and polystyrene substrates, confirming the selectivity of these surfaces for the metastable form III, which shows an unprecedented stability over a time span of several months. A mechanism of its transformation to phase II, via a partially ordered intermediate state, is suggested by polarized Raman measurements.

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.

Tandem In Situ Monitoring for Quantitative Assessment of Mechanochemical Reactions Involving Structurally Unknown Phases

We report herein quantitative in situ monitoring by simultaneous PXRD and Raman spectroscopy of the mechanochemical reaction between benzoic acid and nicotinamide, affording a rich polymorphic system with four new cocrystal polymorphs, multiple phase transformations, and a variety of reaction pathways. After observing polymorphs by in situ monitoring, we were able to isolate and characterize three of the four polymorphs, most of which are not accessible from solution. Relative stabilities among the isolated polymorphs at ambient conditions were established by slurry experiments. Using two complementary methods for in situ monitoring enabled quantitative assessment and kinetic analysis of each studied mechanochemical reaction, even when involving unknown crystal structures, and short-lived intermediates. In situ Raman monitoring was introduced here also as a standalone laboratory technique for quantitative assessment of mechanochemical reactions and understanding of mechanochemical reactivity. Our results provide an important step toward a complete and high-throughput quantitative approach to mechanochemical reaction kinetics and mechanisms, necessary for the development of the mechanistic framework of milling reactions.

Enhanced NMR Discrimination of Pharmaceutically Relevant Molecular Crystal Forms through Fragment-Based Ab Initio Chemical Shift Predictions

Chemical shift prediction plays an important role in the determination or validation of crystal structures with solid-state nuclear magnetic resonance (NMR) spectroscopy. One of the fundamental theoretical challenges lies in discriminating variations in chemical shifts resulting from different crystallographic environments. Fragment-based electronic structure methods provide an alternative to the widely used plane wave gauge-including projector augmented wave (GIPAW) density functional technique for chemical shift prediction. Fragment methods allow hybrid density functionals to be employed routinely in chemical shift prediction, and we have recently demonstrated appreciable improvements in the accuracy of the predicted shifts when using the hybrid PBE0 functional instead of generalized gradient approximation (GGA) functionals like PBE. Here, we investigate the solid-state ¹³C and ¹⁵N NMR spectra for multiple crystal forms of acetaminophen, phenobarbital, and testosterone. We demonstrate that the use of the hybrid density functional instead of a GGA provides both higher accuracy in the chemical shifts and increased discrimination among the different crystallographic environments. Finally, these results also provide compelling evidence for the transferability of the linear regression parameters mapping predicted chemical shieldings to chemical shifts that were derived in an earlier study.

Monitoring model drug microencapsulation in PLGA scaffolds using X-ray powder diffraction

The microencapsulation of three model drugs; metronidazole, paracetamol and sulphapyridine into Poly (dl-Lactide-Co-Glycolide) (PLGA) scaffolds were probed using X-ray Powder Diffraction (XRPD). Changes in the diffraction patterns of the PLGA scaffolds after encapsulation was suggestive of a chemical interaction between the pure drugs and the scaffolds and not a physical intermixture.

Controlled production of the elusive metastable form II of acetaminophen (paracetamol): a fully scalable templating approach in a cooling environment

A designed multi-component templating approach to paracetamol form-II provides a reproducible route to scalable (>100 g in a single run) production of this previously elusive polymorph, with its enhanced compressibility and solubility.

Stabilisation of metastable polymorphs: the case of paracetamol form III

A unique approach is applied in preparing the first ‘air-stable’ formulation of form III of paracetamol, long regarded as a highly unstable form.

Real-Time and In Situ Monitoring of Mechanochemical Reactions: A New Playground for All Chemists

We provide a brief overview of the first techniques for direct, real-time observation of mechanochemical reactions by milling. Whereas mechanisms and kinetics of solid-state reactions induced by temperature or pressure have been extensively investigated, transformations of materials under continuous impact in a milling assembly remain largely unexplored and based on ex situ studies. The recent introduction and development of techniques for in situ monitoring of milling reactions by synchrotron X-ray powder diffraction and Raman spectroscopy has enabled the first direct insight into milling mechanochemistry, opening a new area for studies of chemical reactivity. So far, these techniques have revealed rapid, multistep reaction mechanisms and metastable intermediates that are impossible or difficult to observe or isolate in solution and have highlighted shortcomings of ex situ mechanistic studies. These pioneering advances also highlight the low level of mechanistic understanding and future challenges in developing a clear mechanistic picture of physicochemical transformations by milling.

Quantitative in situ and real-time monitoring of mechanochemical reactions

An experimental technique for in situ and real-time monitoring of mechanochemical reactions in a shaker ball mill was recently described, which utilises highly penetrating X-ray radiation available at the ID15B beamline of the European Synchrotron Radiation Facility. Herein, we describe the first attempts to perform such reaction monitoring in a quantitative fashion, by introducing an internal X-ray diffraction standard. The use of silicon as an internal standard resolved the issue with variations of the amount of the sample in the X-ray beam due to the non-uniform distribution of the sample in the reaction jar and allowed, via Rietveld analysis, the first quantitative estimate of the amorphous phase content in a mechanochemical reaction as it is being milled. We also highlight problems associated with the non-ideal mixing of the reaction mixture.

Understanding pharmaceutical polymorphic transformations I: influence of process variables and storage conditions

The active pharmaceutical ingredient (API) of a dosage form is affected by number of mechanical and environmental factors which have a tendency to alter its crystalline state. Polymorphic transitions have been observed to occur during various unit operations like granulation, milling and compression. Forces of pressure, shear and temperature have an ability to induce alterations in crystal habit. A conversion in polymorphic form during a unit operation is very likely to affect the handling of API in the subsequent unit operation. Transitions have also been observed during storage of formulations where the relative humidity and temperature play a major role. An increase in temperature during storage can dehydrate or desolvate the crystal and hence produce crystal defects, whilst, high humidity conditions produce higher molecular mobility leading to either crystallization of API or alteration of its crystalline form.

In situ Raman mapping for identifying transient solid forms

In situ Raman surface mapping on a trihydrate sample provides evidence for the existence of a transient phase X on the route to its stable anhydrous polymorph I. The phase transformation events are observed in the order trihydrate → transient phase → form I.

Surface Mediated Structures: Stabilization of Metastable Polymorphs on the Example of Paracetamol

The preparation of typically thermodynamically unstable polymorphic structures is a challenge. However, solid surfaces are well established aids for the formation and stabilization of polymorphic structures within, for instance, organic electronics. In this study, we report the stabilization of a pharmaceutically relevant substance via a solid surface at ambient conditions. Form III of paracetamol, which is typically unstable in the bulk at standard conditions, can be stabilized with a model silica surface by a standard spin coating procedure followed by rapid heat treatment. Such a preparation technique allows the use of atomic force microscopy and grazing incidence X-ray diffraction measurements revealing detailed information on the morphology and structure of the polymorph. Furthermore, the results exhibit that this polymorph is stable over a long period of time revealing surface mediated stabilization. These findings demonstrate a novel approach to provide thermodynamic stability when applied to similar molecules with specific applications.

Study of the single-crystal X-ray diffuse scattering in paracetamol polymorphs

Single-crystal diffuse X-ray scattering from paracetamol polymorphs is successfully calculated with Monte Carlo (MC) models that are used to simulate the crystals. In order to obtain the correct model appropriate force constants are required that describe the interatomic potentials used in the MC algorithm. Coefficients for an empirical `Buckingham'-type formula are used to determine these force constants. These coefficients are subsequently refined using the least-squares method and are found to converge on similar values for both polymorphic forms. An investigation of the correlation space generated from each model provides what would be expected given that strong displacive correlations exist between the molecules comprising the densely hydrogen-bonded layers. More disordered motions between these layers are present in the model for form (II) as opposed to form (I). An investigation into the peculiarities of librational disorder was also conducted, however, correlation values turn out to be so small that any structural information concerning librational correlation is inconclusive. The purpose of this experiment was to identify if the diffuse scattering features could provide further insight into understanding the physical reasoning behind the metastability of form (II). The form (II) → (I) phase transition is also not currently well understood and usually phase transitional information can be obtained from pronounced diffuse scattering features. Since the diffuse scattering is modelled adequately using harmonic potentials it is our conjecture that the `diffuse' is essentially thermal in origin and does not afford any extra information about the form (II) → (I) phase transition.