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

Exploring the Solid-State Landscape of Carbamazepine during Dehydration: A Low Frequency Raman Spectroscopy Perspective

The solid-state landscape of carbamazepine during its dehydration was explored using Raman spectroscopy in the low- (-300 to -15, 15 to 300) and mid- (300 to 1800 cm^-1) frequency spectral regions. Carbamazepine dihydrate and forms I, III, and IV were also characterized using density functional theory with periodic boundary conditions and showed good agreement with experimental Raman spectra with mean average deviations less than 10 cm^-1. The dehydration of carbamazepine dihydrate was examined under different temperatures (40, 45, 50, 55, and 60 °C). Principal component analysis and multivariate curve resolution were used to explore the transformation pathways of different solid-state forms during the dehydration of carbamazepine dihydrate. The low-frequency Raman domain was able to detect the rapid growth and subsequent decline of carbamazepine form IV, which was not as effectively observed by mid-frequency Raman spectroscopy. These results showcased the potential benefits of low-frequency Raman spectroscopy for pharmaceutical process monitoring and control.

A data-driven and topological mapping approach for the a priori prediction of stable molecular crystalline hydrates

Predictions of the structures of stoichiometric, fractional, or nonstoichiometric hydrates of organic molecular crystals are immensely challenging due to the extensive search space of different water contents, host molecular placements throughout the crystal, and internal molecular conformations. However, the dry frameworks of these hydrates, especially for nonstoichiometric or isostructural dehydrates, can often be predicted from a standard anhydrous crystal structure prediction (CSP) protocol. Inspired by developments in the field of drug binding, we introduce an efficient data-driven and topologically aware approach for predicting organic molecular crystal hydrate structures through a mapping of water positions within the crystal structure. The method does not require a priori specification of water content and can, therefore, predict stoichiometric, fractional, and nonstoichiometric hydrate structures. This approach, which we term a mapping approach for crystal hydrates (MACH), establishes a set of rules for systematic determination of favorable positions for water insertion within predicted or experimental crystal structures based on considerations of the chemical features of local environments and void regions. The proposed approach is tested on hydrates of three pharmaceutically relevant compounds that exhibit diverse crystal packing motifs and void environments characteristic of hydrate structures. Overall, we show that our mapping approach introduces an advance in the efficient performance of hydrate CSP through generation of stable hydrate stoichiometries at low cost and should be considered an integral component for CSP workflows.

Solid-State Forms of Koumine Hydrochloride: Phase Transformations and the Crystal Structure and Properties of the Stable Form

To investigate the solid-state forms of koumine hydrochloride (KMY), solid form screening was performed, and one amorphous form and five crystalline forms (forms A, B, C, D, and E) were identified by powder X-ray diffraction. Form A was the dominant crystal product, and its crystal structure and packing pattern were determined by single-crystal X-ray diffraction. The crystals displayed an orthorhombic crystal system and symmetry of space group P2₁2₁2₁ with Z' = 1. The amorphous form transformed to form A at 105-120 °C or 75% RH, while forms B, C, D, and E could only be intermediate phases and readily transformed to form A at room temperature. Therefore, the phase transformations of KMY solid-state forms were established. The properties of the amorphous form and form A were further elucidated by applying vibrational spectroscopy, moisture sorption analysis, and thermal analysis. Accordingly, form A, the KMY anhydrate, was found to be the thermodynamically stable form with low hygroscopicity under ambient conditions. These characteristics are crucial in the manufacture and storage of active pharmaceutical ingredients.

Emerging Landscape of Computational Modeling in Pharmaceutical Development

Computational chemistry applications have become an integral part of the drug discovery workflow over the past 35 years. However, computational modeling in support of drug development has remained a relatively uncharted territory for a significant part of both academic and industrial communities. This review considers the computational modeling workflows for three key components of drug preclinical and clinical development, namely, process chemistry, analytical research and development, as well as drug product and formulation development. An overview of the computational support for each step of the respective workflows is presented. Additionally, in context of solid form design, special consideration is given to modern physics-based virtual screening methods. This covers rational approaches to polymorph, coformer, counterion, and solvent virtual screening in support of solid form selection and design.

Toward an Understanding of the Propensity for Crystalline Hydrate Formation by Molecular Compounds. Part 2

The propensity of molecular organic compounds to form stoichiometric or nonstoichiometric crystalline hydrates remains a challenging aspect of crystal engineering and is of practical relevance to fields such as pharmaceutical science. In this work, we address the propensity for hydrate formation of a library of eight compounds comprised of 5- and 6-membered N-heterocyclic aromatics classified into three subgroups: linear dipyridyls, substituted Schiff bases, and tripodal molecules. Each molecular compound studied possesses strong hydrogen bond acceptors and is devoid of strong hydrogen bond donors. Four methods were used to screen for hydrate propensity using the anhydrate forms of the molecular compounds in our library: water slurry under ambient conditions, exposure to humidity, aqueous solvent drop grinding (SDG), and dynamic water vapor sorption (DVS). In addition, crystallization from mixed solvents was studied. Water slurry, aqueous SDG, and exposure to humidity were found to be the most effective methods for hydrate screening. Our study also involved a structural analysis using the Cambridge Structural Database, electrostatic potential (ESP) maps, full interaction maps (FIMs), and crystal packing motifs. The hydrate propensity of each compound studied was compared to a compound of the same type known to form a hydrate through a previous study of ours. Out of the eight newly studied compounds (herein numbered 4-11), three Schiff bases were observed to form hydrates. Three crystal structures (two hydrates and one anhydrate) were determined. Compound 6 crystallized as an isolated site hydrate in the monoclinic space group P2₁/a, while 7 and 10 crystallized in the monoclinic space group P2₁/c as a channel tetrahydrate and an anhydrate, respectively. Whereas we did not find any direct correlation between the number of H-bond acceptors and either hydrate propensity or the stoichiometry of the resulting hydrates, analysis of FIMs suggested that hydrates tend to form when the corresponding anhydrate structure does not facilitate intermolecular interactions.

Challenges and opportunities of pharmaceutical cocrystals: a focused review on non-steroidal anti-inflammatory drugs

The focus of the review is to discuss the relevant and essential aspects of pharmaceutical cocrystals in both academia and industry with an emphasis on non-steroidal anti-inflammatory drugs (NSAIDs). Although cocrystals have been prepared for a plethora of drugs, NSAID cocrystals are focused due to their humongous application in different fields of medication such as antipyretic, anti-inflammatory, analgesic, antiplatelet, antitumor, and anti-carcinogenic drugs. The highlights of the review are (a) background of cocrystals and other solid forms of an active pharmaceutical ingredient (API) based on the principles of crystal engineering, (b) why cocrystals are an excellent opportunity in the pharma industry, (c) common methods of preparation of cocrystals from the lab scale to bulk quantity, (d) some latest case studies of NSAIDs which have shown better physicochemical properties for example; mechanical properties (tabletability), hydration, solubility, bioavailability, and permeability, and (e) latest guidelines of the US FDA and EMA opening new opportunities and challenges.

Combined Use of Structure Analysis, Studies of Molecular Association in Solution, and Molecular Modelling to Understand the Different Propensities of Dihydroxybenzoic Acids to Form Solid Phases

The arrangement of hydroxyl groups in the benzene ring has a significant effect on the propensity of dihydroxybenzoic acids (diOHBAs) to form different solid phases when crystallized from solution. All six diOHBAs were categorized into distinctive groups according to the solid phases obtained when crystallized from selected solvents. A combined study using crystal structure and molecule electrostatic potential surface analysis, as well as an exploration of molecular association in solution using spectroscopic methods and molecular dynamics simulations were used to determine the possible mechanism of how the location of the phenolic hydroxyl groups affect the diversity of solid phases formed by the diOHBAs. The crystal structure analysis showed that classical carboxylic acid homodimers and ring-like hydrogen bond motifs consisting of six diOHBA molecules are prominently present in almost all analyzed crystal structures. Both experimental spectroscopic investigations and molecular dynamics simulations indicated that the extent of intramolecular bonding between carboxyl and hydroxyl groups in solution has the most significant impact on the solid phases formed by the diOHBAs. Additionally, the extent of hydrogen bonding with solvent molecules and the mean lifetime of solute-solvent associates formed by diOHBAs and 2-propanol were also investigated.

Insight into amoxicillin sodium heterosolvates and non-solvated form: crystal structures, phase transformation behaviors, and desolvation mechanism

Methanol and methyl acetate molecules played different roles in the formation of S_M–M amoxicillin sodium heterosolvate. A two-step desolvation mechanism of novel heterosolvates was proposed.

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.

Supramolecular organisation of sulphate salt hydrates exemplified with brucine sulphate

The frequency of hydrate formation among organic sulphate salts is unravelled. Interconversion of the hydrates of brucine sulphate occurs with small changes in the relative humidity.

Design, screening, and properties of novel solvates of azoxystrobin based on isomorphism

The molecular size of the solvent is important for the formation of isomorphic azoxystrobin solvates.

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.

Explaining crystallization preferences of two polyphenolic diastereoisomers by crystal structure prediction

Crystal structure prediction is used to understand the differences in crystallization of catechin and epicatechin, and to explore the predictability of solvate formation.

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.

From serendipity to supramolecular design: assessing the utility of computed crystal form landscapes in inferring the risks of crystal hydration in carboxylic acids

Calculated structural descriptors for predicted anhydrate polymorphs are used to assess the risks of crystal hydration in carboxylic acids.

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.