Predicting crystal form stability under real-world conditions.

The physicochemical properties of molecular crystals, such as solubility, stability, compactability, melting behaviour and bioavailability, depend on their crystal form1. In silico crystal form selection has recently come much closer to realization because of the development of accurate and affordable free-energy calculations2-4. Here we redefine the state of the art, primarily by improving the accuracy of free-energy calculations, constructing a reliable experimental benchmark for solid-solid free-energy differences, quantifying statistical errors for the computed free energies and placing both hydrate crystal structures of different stoichiometries and anhydrate crystal structures on the same energy landscape, with defined error bars, as a function of temperature and relative humidity. The calculated free energies have standard errors of 1-2 kJ mol-1 for industrially relevant compounds, and the method to place crystal structures with different hydrate stoichiometries on the same energy landscape can be extended to other multi-component systems, including solvates. These contributions reduce the gap between the needs of the experimentalist and the capabilities of modern computational tools, transforming crystal structure prediction into a more reliable and actionable procedure that can be used in combination with experimental evidence to direct crystal form selection and establish control5.

Directing Crystallization Outcomes of Conformationally Flexible Molecules: Polymorphs, Solvates, and Desolvation Pathways of Fluconazole.

Control over polymorphism and solvatomorphism in API assisted by structural information, e.g., molecular conformation or associations via hydrogen bonds, is crucial for the industrial development of new drugs, as the crystallization products differ in solubility, dissolution profile, compressibility, or melting temperature. The stability of the final formulation and technological factors of the pharmaceutical powders further emphasize the importance of precise crystallization protocols. This is particularly important when working with highly flexible molecules with considerable conformational freedom and a large number of hydrogen bond donors or acceptors (e.g., fluconazole, FLU). Here, cooling and suspension crystallization were applied to access polymorphs and solvates of FLU, a widely used azole antifungal agent with high molecular flexibility and several reported polymorphs. Each of four polymorphic forms, FLU I, II, III, or IV, can be obtained from the same set of alcohols (MeOH, EtOH, isPrOH) and DMF via careful control of the crystallization conditions. For the first time, two types of isostructural channel solvates of FLU were obtained (nine new structures). Type I solvates were prepared by cooling crystallization in Tol, ACN, DMSO, BuOH, and BuON. Type II solvates formed in DCM, ACN, nPrOH, and BuOH during suspension experiments. We propose desolvation pathways for both types of solvates based on the structural analysis of the newly obtained solvates and their desolvation products. Type I solvates desolvate to FLU form I by hydrogen-bonded chain rearrangements. Type II solvates desolvation leads first to an isomorphic desolvate, followed by a phase transition to FLU form II through hydrogen-bonded dimer rearrangement. Combining solvent-mediated phase transformations with structural analysis and solid-state NMR, supported by periodic electronic structure calculations, allowed us to elucidate the interrelations and transformation pathways of FLU.

In Vitro Investigation of Thiolated Chitosan Derivatives as Mucoadhesive Coating Materials for Solid Lipid Nanoparticles.

In the present study, chitosan (CS) was thiolated by introducing l-cysteine via amide bond formation. Free thiol groups were protected with highly reactive 6-mercaptonicotinic acid (6-MNA) and less-reactive l-cysteine, respectively, via thiol/disulfide-exchange reactions. Unmodified CS, l-cysteine-modified thiolated CS (CS-Cys), 6-MNA-S-protected thiolated CS (CS-Cys-MNA), and l-cysteine-S-protected thiolated CS (CS-Cys-Cys) were applied as coating materials to solid lipid nanoparticles (SLN). The strength of mucus interaction followed the rank order plain < CS < CS-Cys-Cys < CS-Cys < CS-Cys-MNA, whereas mucus diffusion followed the rank order CS-Cys < CS-Cys-Cys < CS < CS-Cys-MNA < plain. In accordance with lower reactivity, CS-Cys-Cys-coated SLN were immobilized to a lower extent than CS-Cys-coated SLN, while CS-Cys-MNA-coated SLN dissociated from their coating material resulting in a similar diffusion behavior as plain SLN. Consequently, CS-Cys-Cys-coated SLN and CS-Cys-MNA-coated SLN showed the highest retention on porcine intestinal mucosa by enabling a synergism of efficient mucus diffusion and strong mucoadhesion.

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

Five anhydrate polymorphs (forms I-V) and the isomorphic dehydrate (Hydehy) 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 Hydehy. 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 > Hydehy > 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 ( P21/ 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 Hydehy 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.

Experimental and computational approaches to rationalise multicomponent supramolecular assemblies: dapsone monosolvates.

The monosolvate crystal energy landscapes of dapsone (DDS) including the solvents carbon tetrachloride, acetone, cyclohexanone, dimethyl formamide, tetrahydrofuran, methyl ethyl ketone, 1,2-dichloroethane, 1,4-dioxane, dichloromethane and chloroform were established using experimental and computational approaches. To rationalise and understand solvate formation, solvate stability and desolvation reactions a careful control of the experimental crystallisation and storage conditions, a range of thermoanalytical methods and crystal structure prediction were required. Six of the eight DDS monosolvates are reported and characterised for the first time. Structural similarity and diversity of the at ambient conditions unstable monosolvates were apparent from the computed crystal energy landscapes, which had the experimental packings as lowest energy structures. The computed structures were used as input for Rietveld refinements and isostructurality of four of the monosolvates was confirmed. Packing comparisons of the solvate structures and molecular properties of the solvent molecules indicated that both size/shape of the solvent molecule and the possible DDSsolvent interactions are the important factors for DDS solvate formation. Through the combination of experiment and theory solvate stability and structural features have been rationalised.

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.

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-1) 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.

Structural Properties, Order-Disorder Phenomena, and Phase Stability of Orotic Acid Crystal Forms.

Orotic acid (OTA) is reported to exist in the anhydrous (AH), monohydrate (Hy1), and dimethyl sulfoxide monosolvate (SDMSO) forms. In this study we investigate the (de)hydration/desolvation behavior, aiming at an understanding of the elusive structural features of anhydrous OTA by a combination of experimental and computational techniques, namely, thermal analytical methods, gravimetric moisture (de)sorption studies, water activity measurements, X-ray powder diffraction, spectroscopy (vibrational, solid-state NMR), crystal energy landscape, and chemical shift calculations. The Hy1 is a highly stable hydrate, which dissociates above 135 °C and loses only a small part of the water when stored over desiccants (25 °C) for more than one year. In Hy1, orotic acid and water molecules are linked by strong hydrogen bonds in nearly perfectly planar arranged stacked layers. The layers are spaced by 3.1 Å and not linked via hydrogen bonds. Upon dehydration the X-ray powder diffraction and solid-state NMR peaks become broader, indicating some disorder in the anhydrous form. The Hy1 stacking reflection (122) is maintained, suggesting that the OTA molecules are still arranged in stacked layers in the dehydration product. Desolvation of SDMSO, a nonlayer structure, results in the same AH phase as observed upon dehydrating Hy1. Depending on the desolvation conditions, different levels of order-disorder of layers present in anhydrous OTA are observed, which is also suggested by the computed low energy crystal structures. These structures provide models for stacking faults as intergrowth of different layers is possible. The variability in anhydrate crystals is of practical concern as it affects the moisture dependent stability of AH with respect to hydration.

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

Crystal structure prediction studies indicated the existence of an unknown high density monohydrate structure (Hy1B°) as global energy minimum for 4-aminoquinaldine (4-AQ). We thus performed an interdisciplinary experimental and computational study elucidating the crystal structures, solid form inter-relationships, kinetic and thermodynamic stabilities of the stable anhydrate (AH I°), the kinetic monohydrate (Hy1A ) and this novel monohydrate polymorph (Hy1B°) of 4-AQ. The crystal structure of Hy1B° was determined by combining laboratory powder X-ray diffraction data and ab initio calculations. Dehydration studies with differential scanning calorimetry and solubility measurements confirmed the result of the lattice energy calculations, which identified Hy1B° as the thermodynamically most stable hydrate form. At 25 °C the equilibrium of the 4-AQ hydrate/anhydrate system was observed at an aw (water activity) of 0.14. The finding of Hy1B° was complicated by the fact that the metastable but kinetically stable Hy1A shows a higher nucleation and growth rate. The presence of an impurity in an available 4-AQ sample facilitated the nucleation of Hy1B°, whose crystallisation is favored under hydrothermal conditions. The value of combining experimental with theoretical studies in hydrate screening and characterisation, as well as the reasons for hydrate formation in 4-AQ, are discussed.

Navigating the Waters of Unconventional Crystalline Hydrates.

Elucidating the crystal structures, transformations, and thermodynamics of the two zwitterionic hydrates (Hy2 and HyA) of 3-(4-dibenzo[b,f][1,4]oxepin-11-yl-piperazin-1-yl)-2,2-dimethylpropanoic acid (DB7) rationalizes the complex interplay of temperature, water activity, and pH on the solid form stability and transformation pathways to three neutral anhydrate polymorphs (Forms I, II°, and III). HyA contains 1.29 to 1.95 molecules of water per DB7 zwitterion (DB7z). Removal of the essential water stabilizing HyA causes it to collapse to an amorphous phase, frequently concomitantly nucleating the stable anhydrate Forms I and II°. Hy2 is a stoichiometric dihydrate and the only known precursor to Form III, a high energy disordered anhydrate, with the level of disorder depending on the drying conditions. X-ray crystallography, solid state NMR, and H/D exchange experiments on highly crystalline phase pure samples obtained by exquisite control over crystallization, filtration, and drying conditions, along with computational modeling, provided a molecular level understanding of this system. The slow rates of many transformations and sensitivity of equilibria to exact conditions, arising from its varying static and dynamic disorder and water mobility in different phases, meant that characterizing DB7 hydration in terms of simplified hydrate classifications was inappropriate for developing this pharmaceutical.

Solid state forms of 4-aminoquinaldine - From void structures with and without solvent inclusion to close packing.

Polymorphs of 4-aminoquinaldine (4-AQ) have been predicted in silico and experimentally identified and characterised. The two metastable forms, AH (anhydrate) II and AH III, crystallise in the trigonal space group [Formula: see text] and are less densely packed than the thermodynamically most stable phase AH I° (P21/c ). AH II can crystallise and exist both, as a solvent inclusion compound and as an unsolvated phase. The third polymorph, AH III, is exclusively obtained by desolvation of a carbon tetrachloride solvate. Theoretical calculations correctly estimated the experimental 0K stability order, confirmed that AH II can exist without solvents, gave access to the AH III structure, and identified that there exists a subtle balance between close packing and number of hydrogen bonding interactions in the solid state of anhydrous 4-AQ. Furthermore, the prevalence of void space and solvent inclusion in [Formula: see text] structures is discussed.

Insights into hydrate formation and stability of morphinanes from a combination of experimental and computational approaches.

Morphine, codeine, and ethylmorphine are important drug compounds whose free bases and hydrochloride salts form stable hydrates. These compounds were used to systematically investigate the influence of the type of functional groups, the role of water molecules, and the Cl(-) counterion on molecular aggregation and solid state properties. Five new crystal structures have been determined. Additionally, structure models for anhydrous ethylmorphine and morphine hydrochloride dihydrate, two phases existing only in a very limited humidity range, are proposed on the basis of computational dehydration modeling. These match the experimental powder X-ray diffraction patterns and the structural information derived from infrared spectroscopy. All 12 structurally characterized morphinane forms (including structures from the Cambridge Structural Database) crystallize in the orthorhombic space group P212121. Hydrate formation results in higher dimensional hydrogen bond networks. The salt structures of the different compounds exhibit only little structural variation. Anhydrous polymorphs were detected for all compounds except ethylmorphine (one anhydrate) and its hydrochloride salt (no anhydrate). Morphine HCl forms a trihydrate and dihydrate. Differential scanning and isothermal calorimetry were employed to estimate the heat of the hydrate ↔ anhydrate phase transformations, indicating an enthalpic stabilization of the respective hydrate of 5.7 to 25.6 kJ mol(-1) relative to the most stable anhydrate. These results are in qualitative agreement with static 0 K lattice energy calculations for all systems except morphine hydrochloride, showing the need for further improvements in quantitative thermodynamic prediction of hydrates having water···water interactions. Thus, the combination of a variety of experimental techniques, covering temperature- and moisture-dependent stability, and computational modeling allowed us to generate sufficient kinetic, thermodynamic and structural information to understand the principles of hydrate formation of the model compounds. This approach also led to the detection of several new crystal forms of the investigated morphinanes.

Flavone Cocrystals: A Comprehensive Approach Integrating Experimental and Virtual Methods.

The dapsone/flavone cocrystal system served as a benchmark for both experimental and virtual screening methods. Expanding beyond this, two additional active pharmaceutical ingredients (APIs), sulfanilamide and sulfaguanidine, structurally related to dapsone were chosen to investigate the impact of substituents on cocrystal formation. The experimental screening involved mechanochemical methods, slurry experiments, hot-melt extrusion, and the contact preparation method. The virtual screening focused on crystal structure prediction (CSP), molecular complementarity, hydrogen-bond propensity, and molecular electrostatic potentials. The CSP studies not only indicated that each of the three APIs should form cocrystals with flavone but also reproduced the known single- and multicomponent phases. Experimentally, dapsone/flavone cocrystals ACC, BCC, CCC, and DCC were reproduced, CCC was identified as a nonstoichiometric hydrate, and a fifth cocrystal (ECC), a t-butanol solvate, was discovered. The cocrystal polymorphs ACC and BCC are enantiotripically related, and DCC, exhibiting a different stoichiometric ratio, is enthalpically stabilized over the other cocrystals. For the sulfaguanidine/flavone system, two novel, enantiotripically related cocrystals were identified. The crystal structures of two cocrystals and a flavone polymorph were solved from powder X-ray diffraction data, and the stability of all cocrystals was assessed through differential scanning calorimetry and lattice energy calculations. Despite computational indications, a diverse array of cocrystallization techniques did not result in a sulfanilamide/flavone cocrystal. The driving force behind dapsone's tendency to cocrystallize with flavone can be attributed to the overall strength of flavone interactions in the cocrystals. For sulfaguanidine, the potential to form strong API···API and API···coformer interactions in the cocrystal is a contributing factor. Furthermore, flavone was found to be trimorphic.

Comprehensive Insights into Sulfaguanidine in the Solid State: An Experimental and Computational Study.

A thorough re-examination of sulfaguanidine's (SGD) solid-state behavior was conducted, 65 years after the initial report on SGD polymorphism. This investigation focuses on the polymorphic nature of the compound, the formation of hydrates and solvates, and the pivotal role of experimental and computational methods in screening, assessing stability, and understanding transformation processes. The findings confirm the presence of five anhydrates (AH-I-V), two monohydrate polymorphs (Hy1-I and Hy1-II), and nine solvates (with tetrahydrofuran, methanol, ethanol, t-butanol, acetone, cyclohexanone, dimethyl sulfoxide, dimethyl formamide, and dimethyl acetamide). Notably, nine novel structures-two anhydrates and seven solvates-are reported, solved from powder X-ray diffraction data. Calorimetric measurements have revealed that AH-II is the thermodynamically stable polymorph at room and low temperatures. In contrast, AH-I emerges as the stable polymorph at higher temperatures, yet it exhibits remarkable kinetic stability at RT and demonstrates greater stability in terms of hydration. The anhydrate forms exhibit distinctive packing arrangements, while the two hydrates share a close structural resemblance. Among the seven structurally characterized solvates, only the tetrahydrofuran and dimethyl sulfoxide solvates are isostructural. Controlled desolvation experiments enabled the formation of AH-I, AH-II, and, notably, AH-V for the first time. The anhydrate and monohydrate crystal structure prediction studies reveal that the computed lowest-energy structures correspond to experimentally observed forms and propose models for the elusive AH-IV structure. Overall, the exploration of SGD's solid-state landscape confirms a rich array of highly stable H-bonding motifs and packing arrangements, positioning this study as an ideal model for complex solid-state systems and shedding light on its intricate solid-state nature.

Suppression of Disorder in Benzamide and Thiobenzamide Crystals by Fluorine Substitution.

Disorder is a common feature of molecular crystals that complicates determination of structures and can potentially affect electric and mechanical properties. Suppression of disorder is observed in otherwise severely disordered benzamide and thiobenzamide crystals by substituting hydrogen with fluorine in the ortho-position of the phenyl ring. Fluorine occupancies of 20-30% are sufficient to suppress disorder without changing the packing motif. Crystal structure prediction calculations reveal a much denser lattice energy landscape for benzamide compared to 2-fluorobenzamide, suggesting that fluorine substitution makes disorder less likely.

The seventh blind test of crystal structure prediction: structure ranking methods.

A seventh blind test of crystal structure prediction has been organized by the Cambridge Crystallographic Data Centre. The results are presented in two parts, with this second part focusing on methods for ranking crystal structures in order of stability. The exercise involved standardized sets of structures seeded from a range of structure generation methods. Participants from 22 groups applied several periodic DFT-D methods, machine learned potentials, force fields derived from empirical data or quantum chemical calculations, and various combinations of the above. In addition, one non-energy-based scoring function was used. Results showed that periodic DFT-D methods overall agreed with experimental data within expected error margins, while one machine learned model, applying system-specific AIMnet potentials, agreed with experiment in many cases demonstrating promise as an efficient alternative to DFT-based methods. For target XXXII, a consensus was reached across periodic DFT methods, with consistently high predicted energies of experimental forms relative to the global minimum (above 4 kJ mol-1 at both low and ambient temperatures) suggesting a more stable polymorph is likely not yet observed. The calculation of free energies at ambient temperatures offered improvement of predictions only in some cases (for targets XXVII and XXXI). Several avenues for future research have been suggested, highlighting the need for greater efficiency considering the vast amounts of resources utilized in many cases.

The seventh blind test of crystal structure prediction: structure generation methods.

A seventh blind test of crystal structure prediction was organized by the Cambridge Crystallographic Data Centre featuring seven target systems of varying complexity: a silicon and iodine-containing molecule, a copper coordination complex, a near-rigid molecule, a cocrystal, a polymorphic small agrochemical, a highly flexible polymorphic drug candidate, and a polymorphic morpholine salt. In this first of two parts focusing on structure generation methods, many crystal structure prediction (CSP) methods performed well for the small but flexible agrochemical compound, successfully reproducing the experimentally observed crystal structures, while few groups were successful for the systems of higher complexity. A powder X-ray diffraction (PXRD) assisted exercise demonstrated the use of CSP in successfully determining a crystal structure from a low-quality PXRD pattern. The use of CSP in the prediction of likely cocrystal stoichiometry was also explored, demonstrating multiple possible approaches. Crystallographic disorder emerged as an important theme throughout the test as both a challenge for analysis and a major achievement where two groups blindly predicted the existence of disorder for the first time. Additionally, large-scale comparisons of the sets of predicted crystal structures also showed that some methods yield sets that largely contain the same crystal structures.