4-Aminoquinaldine monohydrate polymorphism: Prediction and impurity aided discovery of a difficult to access stable form

Frontiers of molecular crystal structure prediction for pharmaceuticals and functional organic materials

The reliability of organic molecular crystal structure prediction has improved tremendously in recent years. Crystal structure predictions for small, mostly rigid molecules are quickly becoming routine. Structure predictions for larger, highly flexible molecules are more challenging, but their crystal structures can also now be predicted with increasing rates of success. These advances are ushering in a new era where crystal structure prediction drives the experimental discovery of new solid forms. After briefly discussing the computational methods that enable successful crystal structure prediction, this perspective presents case studies from the literature that demonstrate how state-of-the-art crystal structure prediction can transform how scientists approach problems involving the organic solid state. Applications to pharmaceuticals, porous organic materials, photomechanical crystals, organic semi-conductors, and nuclear magnetic resonance crystallography are included. Finally, efforts to improve our understanding of which predicted crystal structures can actually be produced experimentally and other outstanding challenges are discussed.

Global analysis of the energy landscapes of molecular crystal structures by applying the threshold algorithm

Polymorphism in molecular crystals has important consequences for the control of materials properties and our understanding of crystallization. Computational methods, including crystal structure prediction, have provided important insight into polymorphism, but have usually been limited to assessing the relative energies of structures. We describe the implementation of the Monte Carlo threshold algorithm as a method to provide an estimate of the energy barriers separating crystal structures. By sampling the local energy minima accessible from multiple starting structures, the simulations yield a global picture of the crystal energy landscapes and provide valuable information on the depth of the energy minima associated with crystal structures. We present results from applying the threshold algorithm to four polymorphic organic molecular crystals, examine the influence of applying space group symmetry constraints during the simulations, and discuss the relationship between the structure of the energy landscape and the intermolecular interactions present in the crystals.

Progress in understanding crystallisation: a personal perspective

After this discussion meeting, most participants felt that we do not understand crystallisation. However, in the 1980s I believe that most scientists would have considered that crystallisation was adequately understood....

The trimorphism of 3-hydroxybenzoic acid: an experimental and computational study

A computationally driven experimental search for polymorphs of 3-hydroxybenzoic acid confirmed the third form and the small energy differences between the polymorphs.

Crystal forms in pharmaceutical applications: olanzapine, a gift to crystal chemistry that keeps on giving

This contribution reviews the efforts of many scientists around the world to discover and structurally characterize olanzapine crystal forms, clearing up inconsistencies in the scientific and patent literature and highlighting the challenges in identifying new forms amidst 60+ known polymorphs and solvates. Owing to its remarkable solid-state chemistry, olanzapine has emerged over the last three decades as a popular tool compound for developing new experimental and computational methods for enhanced molecular level understanding of solid-state structure, form diversity and crystallization outcomes. This article highlights the role of olanzapine in advancing the fundamental understanding of crystal forms, interactions within crystal structures, and growth units in molecular crystallization, as well as influencing the way in which drugs are developed today.

Pharmaceutical Hydrates Analysis-Overview of Methods and Recent Advances

This review discusses a set of instrumental and computational methods that are used to characterize hydrated forms of APIs (active pharmaceutical ingredients). The focus has been put on highlighting advantages as well as on presenting some limitations of the selected analytical approaches. This has been performed in order to facilitate the choice of an appropriate method depending on the type of the structural feature that is to be analyzed, that is, degree of hydration, crystal structure and dynamics, and (de)hydration kinetics. The presented techniques include X-ray diffraction (single crystal X-ray diffraction (SCXRD), powder X-ray diffraction (PXRD)), spectroscopic (solid state nuclear magnetic resonance spectroscopy (ssNMR), Fourier-transformed infrared spectroscopy (FT-IR), Raman spectroscopy), thermal (differential scanning calorimetry (DSC), thermogravimetric analysis (TGA)), gravimetric (dynamic vapour sorption (DVS)), and computational (molecular mechanics (MM), Quantum Mechanics (QM), molecular dynamics (MD)) methods. Further, the successful applications of the presented methods in the studies of hydrated APIs as well as studies on the excipients' influence on these processes have been described in many examples.

Periodic DFT Calculations-Review of Applications in the Pharmaceutical Sciences

In the introduction to this review the complex chemistry of solid-state pharmaceutical compounds is summarized. It is also explained why the density functional theory (DFT) periodic calculations became recently so popular in studying the solid APIs (active pharmaceutical ingredients). Further, the most popular programs enabling DFT periodic calculations are presented and compared. Subsequently, on the large number of examples, the applications of such calculations in pharmaceutical sciences are discussed. The mentioned topics include, among others, validation of the experimentally obtained crystal structures and crystal structure prediction, insight into crystallization and solvation processes, development of new polymorph synthesis ways, and formulation techniques as well as application of the periodic DFT calculations in the drug analysis.

A novel lurasidone hydrochloride–shikimic acid co-amorphous system formed by hydrogen-bonding interaction with the retained pH-dependent solubility behavior

The current study was aimed at investigating the lurasidone hydrochloride–shikimic acid co-amorphous system using a new type of organic acid.

Surface Induced Phenytoin Polymorph. 1. Full Structure Solution by Combining Grazing Incidence X-ray Diffraction and Crystal Structure Prediction

Understanding the behavior and properties of molecules assembled in thin layers requires knowledge of their crystalline packing. The drug phenytoin (5,5-diphenylhydantoin) is one of the compounds that can be grown as a surface induced polymorph. By using grazing incidence X-ray diffraction, the monoclinic unit cell of the new form II can be determined, but, due to crystal size and the low amount of data, a full solution using conventional structure solving strategies fails. In this work, the full solution has been obtained by combining computational structure generation and experimental results. The comparison between the bulk and the new surface induced phase reveals significant packing differences of the hydrogen-bonding network, which might be the reason for the faster dissolution of form II with respect to form I. The results are very satisfactory, and the method might be adapted for other systems, where, due to the limited amount of experimental data, one must rely on additional approaches to gain access to more detailed information to understand the solid-state behavior.

Dapsone Form V: A Late Appearing Thermodynamic Polymorph of a Pharmaceutical

Five anhydrate polymorphs (forms I-V) and the isomorphic dehydrate (Hy_dehy) of dapsone (4,4'-diaminodiphenyl sulfone or DDS) were prepared and characterized in an interdisciplinary experimental and computational study, elucidating the kinetic and thermodynamic stabilities, solid form interrelationships, and structural features of the known forms I-IV, the novel polymorph form V, and Hy_dehy. Calorimetric measurements, solubility experiments, and lattice energy calculations revealed that form V is the thermodynamically stable polymorph from absolute zero to at least 90 °C. At higher temperatures, form II, and then form I, becomes the most stable DDS solid form. The computed 0 K stability order (lattice energy calculations) was confirmed with calorimetric measurements as follows, V (most stable) > III > Hy_dehy > II > I > IV (least stable). The discovery of form V was complicated by the fact that the metastable but kinetically stabilized form III shows a higher nucleation and growth rate. By combining laboratory powder X-ray diffraction data and ab initio calculations, the crystal structure of form V ( P2₁/ c, Z' = 4) was solved, with a high energy DDS conformation allowing a denser packing and more stable intermolecular interactions, rationalizing the formation of a high Z' structure. The structures of the forms I and IV, only observed from the melt and showing distinct packing features compared to the forms II, III, and V, were derived from the computed crystal energy landscapes. Dehydration modeling of the DDS hydrate led to the Hy_dehy structure. This study expands our understanding about the complex crystallization behavior of pharmaceuticals and highlights the big challenge in solid form screening, especially that there is no clear end point.

Insight into the role of piperazine in the thermodynamics and nucleation kinetics of the triethylenediamine–methyl tertiary butyl ether system

Considering the role of piperazine in the TEDA–MTBE system, we investigate how to control the crystal size distribution.

Control and prediction of the organic solid state: a challenge to theory and experiment†

The ability of theoretical chemists to quantitatively model the weak forces between organic molecules is being exploited to predict their crystal structures and estimate their physical properties. Evolving crystal structure prediction methods are increasingly being used to aid the design of organic functional materials and provide information about thermodynamically plausible polymorphs of speciality organic materials to aid, for example, pharmaceutical development. However, the increasingly sophisticated experimental studies for detecting the range of organic solid-state behaviours provide many challenges for improving quantitative theories that form the basis for the computer modelling. It is challenging to calculate the relative thermodynamic stability of different organic crystal structures, let alone understand the kinetic effects that determine which polymorphs can be observed and are practically important. However, collaborations between experiment and theory are reaching the stage of devising experiments to target the first crystallization of new polymorphs or create novel organic molecular materials.

Supramolecular Organization of Nonstoichiometric Drug Hydrates: Dapsone

The observed moisture- and temperature dependent transformations of the dapsone (4,4'-diaminodiphenyl sulfone, DDS) 0. 33-hydrate were correlated to its structure and the number and strength of the water-DDS intermolecular interactions. A combination of characterization techniques was used, including thermal analysis (hot-stage microscopy, differential scanning calorimetry and thermogravimetric analysis), gravimetric moisture sorption/desorption studies and variable humidity powder X-ray diffraction, along with computational modeling (crystal structure prediction and pair-wise intermolecular energy calculations). Depending on the relative humidity the hydrate contains between 0 and 0.33 molecules of water per molecule DDS. The crystal structure is retained upon dehydration indicating that DDS hydrate shows a non-stoichiometric (de)hydration behavior. Unexpectedly, the water molecules are not located in structural channels but at isolated-sites of the host framework, which is counterintuitively for a hydrate with non-stoichiometric behavior. The water-DDS interactions were estimated to be weaker than water-host interactions that are commonly observed in stoichiometric hydrates and the lattice energies of the isomorphic dehydration product (hydrate structure without water molecules) and (form III) differ only by ~1 kJ mol-1. The computational generation of hypothetical monohydrates confirms that the hydrate with the unusual DDS:water ratio of 3:1 is more stable than a feasible monohydrate structure. Overall, this study highlights that a deeper understanding of the formation of hydrates with non-stoichiometric behavior requires a multidisciplinary approach including suitable experimental and computational methods providing a firm basis for the development and manufacturing of high quality drug products.

Is zeroth order crystal structure prediction (CSP_0) coming to maturity? What should we aim for in an ideal crystal structure prediction code?

Crystal structure prediction based on searching for the global minimum in the lattice energy (CSP_0) is growing in use for guiding the discovery of new materials, for example, new functional materials, new phases of interest to planetary scientists and new polymorphs relevant to pharmaceutical development. This Faraday Discussion can assess the progress of CSP_0 over the range of types of materials to which CSP is currently and could be applied, which depends on our ability to model the variety of interatomic forces in crystals. The basic hypothesis, that the outcome of crystallisation is determined by thermodynamics, needs examining by considering methods of modelling relative thermodynamic stability not only as a function of pressure and temperature, but also of size, solvent and the presence of heterogeneous templates or impurities (CSP_thd). Given that many important materials persist, and indeed may be formed, when they are not the most thermodynamically stable structure, we need to define what would be required of an ideal CSP code (CSP_aim).

New Insights into Solid Form Stability and Hydrate Formation: o-Phenanthroline HCl and Neocuproine HCl

The moisture- and temperature dependent stabilities and interrelation pathways of the practically relevant solid forms of o-phenanthroline HCl (1) and neocuproine HCl (2) were investigated using thermal analytical techniques (HSM, DSC and TGA) and gravimetric moisture sorption/desorption studies. The experimental stability data were correlated with the structural changes observed upon dehydration and the pairwise interaction and lattice energies calculated. For 1 the monohydrate was identified as the only stable form under conditions of RH typically found during production and storage, but at RH values >80% deliquescence occurs. The second compound, 2, forms an anhydrate and two different hydrates, mono- (2-Hy1) and trihydrate (2-Hy3). The 2-Hy1 structure was solved from SCXRD data and the anhydrate structure derived from a combination of PXRD and CSP. Depending on the environmental conditions (moisture) either 2-Hy1 or 2-Hy3 is the most sable solid form of 2 at RT. The monohydrates 1-Hy1 and 2-Hy1 show a high enthalpic stabilization (≥20 kJ mol⁻¹) relative to the anhydrates. The anhydrates are unstable at ambient conditions and readily transform to the monohydrates even in the presence of traces of moisture. This study demonstrates how the right combination of experiment and theory can unravel the properties and interconversion pathways of solid forms.

Understanding the Role of Water in 1,10-Phenanthroline Monohydrate

The solid forms emerging from an experimental screening programme of 1,10-phenanthroline (o-phen), a heavily used bidentate ligand, and interconversion pathways of its two neat forms, the monohdyrate (Hy1) and four solvates with acetone, chloroform, dichloromethane and 1,2-dichloroethane are described. The solvates, identified and characterised with thermoanalyical methods, are unstable when removed from the mother liquor and desolvate at room temperature depending on the relative humidity (RH) to anhydrate I° (AH I°) or transform to Hy1. At ambient conditions Hy1, a stoichiometric channel hydrate, is the thermodynaically most stable o-phen solid form. The enthalpically stabilised Hy1 melts at 102 °C or dehydrates to AH I° at RH < 10% at 25 °C. The potential energy difference between Hy1 and AH I° was calculated to be approx. 15 kJ mol-1. The second anhydrate polymorph (AH II) can be obatined from the quench cooled melt of o-phen, but is unstable at ambient conditions and transforms wihtin minutes to either AH I° or Hy1. The two neat polymorphs are enantiotropically related and water-free o-phen transforms to Hy1 at RH > 16%. The structural and stablity features of the solid forms, in paricular Hy1, are unravelled by a combination of experimental (thermal analysis, moisture sorption/desorption and storage experiments, infrared spectroscopy and powder X-ray diffraction) and computational modelling (crystal structure prediction and lattice energy calculations), providing a consistent picture why o-phen forms a very stable Z' = 3 channel hydrate.

Unraveling Complexity in the Solid Form Screening of a Pharmaceutical Salt: Why so Many Forms? Why so Few?

The solid form landscape of 5-HT2a antagonist 3-(4-(benzo[d]isoxazole-3-yl)piperazin-1-yl)-2,2-dimethylpropanoic acid hydrochloride (B5HCl) proved difficult to establish. Many crystalline materials were produced by solid form screening, but few forms readily grew high quality crystals to afford a clear picture or understanding of the solid form landscape. Careful control of crystallization conditions, a range of experimental methods, computational modeling of solvate structures, and crystal structure prediction were required to see potential arrangements of the salt in its crystal forms. Structural diversity in the solid form landscape of B5HCl was apparent in the layer structures for the anhydrate polymorphs (Forms I and II), dihydrate and a family of solvates with alcohols. The alcohol solvates, which provided a distinct packing from the neat forms and the dihydrate, form layers with conserved hydrogen bonding between B5HCl and the solvent, as well as stacking of the aromatic rings. The ability of the alcohol hydrocarbon moieties to efficiently pack between the layers accounted for the difficulty in growing some solvate crystals and the inability of other solvates to crystallize altogether. Through a combination of experiment and computation, the crystallization problems, form stability, and desolvation pathways of B5HCl have been rationalized at a molecular level.

Use of Crystal Structure Informatics for Defining the Conformational Space Needed for Predicting Crystal Structures of Pharmaceutical Molecules

Determining the range of conformations that a flexible pharmaceutical-like molecule could plausibly adopt in a crystal structure is a key to successful crystal structure prediction (CSP) studies. We aim to use conformational information from the crystal structures in the Cambridge Structural Database (CSD) to facilitate this task. The conformations produced by the CSD Conformer Generator are reduced in number by considering the underlying rotamer distributions, an analysis of changes in molecular shape, and a minimal number of molecular ab initio calculations. This method is tested for five pharmaceutical-like molecules where an extensive CSP study has already been performed. The CSD informatics-derived set of crystal structure searches generates almost all the low-energy crystal structures previously found, including all experimental structures. The workflow effectively combines information on individual torsion angles and then eliminates the combinations that are too high in energy to be found in the solid state, reducing the resources needed to cover the solid-state conformational space of a molecule. This provides insights into how the low-energy solid-state and isolated-molecule conformations are related to the properties of the individual flexible torsion angles.

Experimental and Computational Hydrate Screening: Cytosine, 5-Flucytosine and Their Solid Solution

The structural, temperature- and moisture dependent stability features of cytosine and 5-flucytosine monohydrates, two pharmaceutically important compounds, were rationalized using complementary experimental and computational approaches. Moisture sorption/desorption, water activity, thermal analysis and calorimetry were applied to determine the stability ranges of hydrate ↔ anhydrate systems, while X-ray diffraction, IR spectroscopy and crystal structure prediction provided the molecular level understanding. At 25 °C, the critical water activity for the cytosine hydrate ↔ anhydrate system is ~0.43 and for 5-flucytosine ~0.41. In 5-flucytosine the water molecules are arranged in open channels, therefore the kinetic desorption data, dehydration < 40% relative humidity (RH), conform with the thermodynamic data, whereas for the cytosine isolated site hydrate dehydration was observed at RH < 15%. Peritectic dissociation temperatures of the hydrates were measured to be 97 °C and 84.2 °C for cytosine and 5-flucytosine, respectively, and the monohydrate to anhydrate transition enthalpies to be around 10 kJ mol-1. Computed crystal energy landscapes not only revealed that the substitution of C5 (H or F) controls the packing and properties of cytosine/5-flucytosine solid forms, but also have enabled the finding of a monohydrate solid solution of the two substances which shows increased thermal- and moisture-dependent stability compared to 5-flucytosine monohydrate.

Prediction and experimental validation of solid solutions and isopolymorphs of cytosine/5-flucytosine

A computational search for polymorphs of cytosine, 5-flucytosine and a 1 : 1 mixture of the two substances not only rationalised the preferred packing arrangements but also enabled the finding and characterisation of cytosine/5-flucytosine solid solutions. The structures of the new solid forms were determined by combining laboratory powder X-ray diffraction data and computational modelling.

Why do Hydrates (Solvates) Form in Small Neutral Organic Molecules? Exploring the Crystal Form Landscapes of the Alkaloids Brucine and Strychnine

Computational methods were used to generate and explore the crystal structure landscapes of the two alkaloids strychnine and brucine. The computed structures were analyzed and rationalized by correlating the modelling results to a rich pool of available experimental data. Despite their structural similarity, the two compounds show marked differences in the formation of solid forms. For strychnine only one anhydrous form is reported in the literature and two new solvates from 1,4-dioxane were detected in the course of this work. In contrast, 22 solid forms are so far known to exist for brucine, comprising two anhydrates, four hydrates (HyA - HyC and a 5.25-hydrate), twelve solvates (alcohols and acetone) and four heterosolvates (mixed solvates with water and alcohols). For strychnine it is hard to produce any solid form other than the stable anhydrate while the formation of specific solid state forms of brucine is governed by a complex interplay between temperature and relative humidity/water activity and it is rather a challenging to avoid hydrate formation. Differences in crystal packing and the high tendency for brucine to form hydrates are not intuitive from the molecular structure alone, as both molecules have hydrogen bond acceptor groups but lack hydrogen bond donor groups. Only the evaluation of the crystal energy landscapes, in particular the close-packed crystal structures and high-energy open frameworks containing voids of molecular (water) dimensions, allowed us to unravel the diverse solid state behavior of the two alkaloids at a molecular level. In this study we demonstrate that expanding the analysis of anhydrate crystal energy landscapes to higher energy structures and calculating the solvent-accessible volume can be used to estimate non-stoichiometric or channel hydrate (solvate) formation, without explicitly computing the hydrate/solvate crystal energy landscapes.

Stoichiometric and Non-Stoichiometric Hydrates of Brucine

The complex interplay of temperature and water activity (a_w) / relative humidity (RH) on the solid form stability and transformation pathways of three hydrates (HyA, HyB and HyC), an isostructural dehydrate (HyA_dehy ), an anhydrate (AH) and amorphous brucine has been elucidated and the transformation enthalpies quantified. The dihydrate (HyA) shows a non-stoichimetric (de)hydration behavior at RH < 40% at 25 °C and the removal of the water molecules results in an isomorphic dehydrate structure. The metastable dehydration product converts to AH upon storage at driest conditions or to HyA if exposed to moisture. HyB is a stoichiometric tetrahydrate. The loss of the water molecules causes HyB to collapse to an amorphous phase. Amorphous brucine transforms to AH at RH < 40% RH and a mixture of hydrated phases at higher RH values. The third hyrdate (HyC) is only stable at RH ≥ 55% at 25 °C and contains 3.65 to 3.85 mole equivalent of water. Dehydration of HyC occurs in one step at RH < 55% at 25 °C or upon heating and AH is obtained. The AH is the thermodynamically most stable phase of brucine at RH < 40% at 25 °C. Depending on the conditions, temperature and a_w, each of the three hydrates becomes the thermodynamically most stable form. This study demonstrates the importance of applying complimentary analytical techniques and appropriate approaches for understanding the stability ranges and transition behavior between the solid forms of compounds with multiple hydrates.

Can computed crystal energy landscapes help understand pharmaceutical solids?

Computational crystal structure prediction (CSP) methods can now be applied to the smaller pharmaceutical molecules currently in drug development. We review the recent uses of computed crystal energy landscapes for pharmaceuticals, concentrating on examples where they have been used in collaboration with industrial-style experimental solid form screening. There is a strong complementarity in aiding experiment to find and characterise practically important solid forms and understanding the nature of the solid form landscape.

Computational and Experimental Characterization of Five Crystal Forms of Thymine: Packing Polymorphism, Polytypism/Disorder and Stoichiometric 0.8-Hydrate

New polymorphs of thymine emerged in an experimental search for solid forms, which was guided by the computationally generated crystal energy landscape. Three of the four anhydrates (AH) are homeoenergetic (A° - C) and their packing modes differ only in the location of oxygen and hydrogen atoms. AHs A° and B are ordered phases, whereas AH C shows disorder (X-ray diffuse scattering). Anhydrates AHs A° and B are ordered phases, whereas AH C shows disorder (X-ray diffuse scattering). Analysis of the crystal energy landscape for alternative AH C hydrogen bonded ribbon motifs identified a number of different packing modes, whose 3D structures were calculated to deviate by less than 0.24 kJ mol-1 in lattice energy. These structures provide models for stacking faults. The three anhydrates A° - C show strong similarity in their powder X-ray diffraction, thermoanalytical and spectroscopic (IR and Raman) characteristics. The already known anhydrate AH A° was identified as the thermodynamically most stable form at ambient conditions; AH B and AH C are metastable but show high kinetic stability. The hydrate of thymine is stable only at water activities (aw) > 0.95 at temperatures ≤ 25 °C. It was found to be a stoichiometric hydrate despite being a channel hydrate with an unusual water:thymine ratio of 0.8:1. Depending on the dehydration conditions, either AH C or AH D is obtained. The hydrate is the only known precursor to AH D. This study highlights the value and complementarity of simultaneous explorations of computationally and experimentally generated solid form landscapes of a small molecule anhydrate ↔ hydrate system.