Absorbing a Little Water: The Structural, Thermodynamic, and Kinetic Relationship between Pyrogallol and Its Tetarto-Hydrate

Designer Gelators for the Crystallization of a Salt Active Pharmaceutical Ingredient-Mexiletine Hydrochloride

We report an approach to obtain drug-mimetic supramolecular gelators, which are capable of stabilizing metastable polymorphs of the pharmaceutical salt mexiletine hydrochloride, a highly polymorphic antiarrhythmic drug. Solution-phase screening led to the discovery of two new solvated solid forms of mexiletine, a type C 1,2,4-trichlorobenzene tetarto-solvate and a type D nitrobenzene solvate. Various metastable forms were crystallized within the gels under conditions which would not have been possible in solution. Despite typically crystallizing concomitantly with form 1, a pure sample of form 3 was crystallized within a gel of ethyl methyl ketone. Various type A channel solvates were crystallized from gels of toluene and ethyl acetate, in which the contents of the channels varied from those of solution-phase forms. Most strikingly, the high-temperature-stable form 2 was crystallized from a gel in 1,2-dibromoethane: the only known route to access this form at room temperature. These results exemplify the powerful stabilizing effect of drug-mimetic supramolecular gels, which can be exploited in pharmaceutical polymorph screens to access highly metastable or difficult-to-nucleate solid forms.

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 peculiar dehydration and solid–solid phase transition of the active pharmaceutical ingredient AZD9898 based on in situ single crystal-to-single crystal transformations

Isostructural dehydration from form A hydrate to form B, and solid–solid phase transition from form B to C of AZD9898 were revealed by in situ single crystal-to-single crystal transformations.

Infrared spectroscopy of flavones and flavonols. Reexamination of the hydroxyl and carbonyl vibrations in relation to the interactions of flavonoids with membrane lipids

Detailed vibrational assignments for twelve flavonoids (seven flavones (flavone, 3- and 5-hydroxyflavone, chrysin, apigenin, fisetin and luteolin) and five flavonols (galangin, kaempferol, quercetin, morin and myricetin)) have been made based on own and reported experimental data and calculations at the B3LYP/6-31+G(d,p) level of theory. All the molecules are treated in a uniform way by using the same set of redundancy-free set of internal coordinates. A generalized harmonic mode mixing is used to corroborate the vibrational characteristics of this important class of molecules. Each flavonoid molecule can be treated from the vibrational point of view as made of relatively weakly coupled chromone and phenyl part. It has been shown that the strongest band around 1600cm^-1 need not be attributable to the CO stretching. The way the vibrations of any of the hydroxyl groups are mixed with ring vibrations and vibrations of other neighboring hydroxyl groups is rather involved. This imposes severe limitations on any attempt to describe normal modes of a flavonol in terms of hydroxyl or carbonyl group vibrations. The role of water molecules in the appearance of flavonoid IR spectra is emphasized. Knowing for the great affinity of phosphate groups in lipids towards water, the immediate consequence is a reasonable assumption that flavonoid lipid interactions is mediated by water.

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.

Diversity of crystal structures and physicochemical properties of ciprofloxacin and norfloxacin salts with fumaric acid

The crystallization of norfloxacin and ciprofloxacin – antibacterial fluoroquinolone compounds – with fumaric acid resulted in the isolation of six distinct solid forms of the drugs with different stoichiometries and hydration levels.

Accurate and efficient representation of intramolecular energy in ab initio generation of crystal structures. I. Adaptive local approximate models

The global search stage of crystal structure prediction (CSP) methods requires a fine balance between accuracy and computational cost, particularly for the study of large flexible molecules. A major improvement in the accuracy and cost of the intramolecular energy function used in the CrystalPredictor II [Habgood et al. (2015). J. Chem. Theory Comput. 11, 1957-1969] program is presented, where the most efficient use of computational effort is ensured via the use of adaptive local approximate model (LAM) placement. The entire search space of the relevant molecule's conformations is initially evaluated using a coarse, low accuracy grid. Additional LAM points are then placed at appropriate points determined via an automated process, aiming to minimize the computational effort expended in high-energy regions whilst maximizing the accuracy in low-energy regions. As the size, complexity and flexibility of molecules increase, the reduction in computational cost becomes marked. This improvement is illustrated with energy calculations for benzoic acid and the ROY molecule, and a CSP study of molecule (XXVI) from the sixth blind test [Reilly et al. (2016). Acta Cryst. B72, 439-459], which is challenging due to its size and flexibility. Its known experimental form is successfully predicted as the global minimum. The computational cost of the study is tractable without the need to make unphysical simplifying assumptions.

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.

Solvent inclusion in the crystal structure of bis-[(adamantan-1-yl)methanaminium chloride] 1,4-dioxane hemisolvate monohydrate explained using the computed crystal energy landscape

Repeated attempts to crystallize 1-adamantane-methyl-amine hydro-chloride as an anhydrate failed but the salt was successfully crystallized as a solvate (2C11H20N+·2Cl-·0.5C4H8O2·H2O), with water and 1,4-dioxane playing a structural role in the crystal and engaging in hydrogen-bonding inter-actions with the cation and anion. Computational crystal-structure prediction was used to rationalize the solvent-inclusion behaviour of this salt by computing the solvent-accessible voids in the predicted low-energy structures for the anhydrate: the global lattice-energy minimum structure, which has the same packing of the ions as the solvate, has solvent-accessible voids that account for 3.71% of the total unit-cell volume and is 6 kJ mol-1 more stable than the next most stable predicted structure.

Using crystal structure prediction to rationalize the hydration propensities of substituted adamantane hydrochloride salts

The crystal energy landscapes of the salts of two rigid pharmaceutically active molecules reveal that the experimental structure of amantadine hydrochloride is the most stable structure with the majority of low-energy structures adopting a chain hydrogen-bond motif and packings that do not have solvent accessible voids. By contrast, memantine hydrochloride which differs in the substitution of two methyl groups on the adamantane ring has a crystal energy landscape where all structures within 10 kJ mol−1of the global minimum have solvent-accessible voids ranging from 3 to 14% of the unit-cell volume including the lattice energy minimum that was calculated after removing water from the hydrated memantine hydrochloride salt structure. The success in using crystal structure prediction (CSP) to rationalize the different hydration propensities of these substituted adamantane hydrochloride salts allowed us to extend the model to predict under blind test conditions the experimental crystal structures of the previously uncharacterized 1-(methylamino)adamantane base and its corresponding hydrochloride salt. Although the crystal structure of 1-(methylamino)adamantane was correctly predicted as the second ranked structure on the static lattice energy landscape, the crystallization of aZ′ = 3 structure of 1-(methylamino)adamantane hydrochloride reveals the limits of applying CSP when the contents of the crystallographic asymmetric unit are unknown.

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.

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.

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 (S_DMSO) 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 S_DMSO, 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.

Hydrous salts of 1-aminoethylidenediphosphonic acid and piperazidine: temperature induced reversible structural transformation in a humid environment

A series of hydrous salts constructed using AEDPH₄ and piperazidine (pip) were synthesized and structurally characterized. A reversible structural transformation circle is observed in single crystals of compounds 1–3.

Quantitative analysis of intermolecular interactions in orthorhombic rubrene

Rubrene is one of the most studied organic semiconductors to date due to its high charge carrier mobility which makes it a potentially applicable compound in modern electronic devices. Previous electronic device characterizations and first principles theoretical calculations assigned the semiconducting properties of rubrene to the presence of a large overlap of the extended π-conjugated core between molecules. We present here the electron density distribution in rubrene at 20 K and at 100 K obtained using a combination of high-resolution X-ray and neutron diffraction data. The topology of the electron density and energies of intermolecular interactions are studied quantitatively. Specifically, the presence of Cπ⋯Cπ interactions between neighbouring tetracene backbones of the rubrene molecules is experimentally confirmed from a topological analysis of the electron density, Non-Covalent Interaction (NCI) analysis and the calculated interaction energy of molecular dimers. A significant contribution to the lattice energy of the crystal is provided by H-H interactions. The electron density features of H-H bonding, and the interaction energy of molecular dimers connected by H-H interaction clearly demonstrate an importance of these weak interactions in the stabilization of the crystal structure. The quantitative nature of the intermolecular interactions is virtually unchanged between 20 K and 100 K suggesting that any changes in carrier transport at these low temperatures would have a different origin. The obtained experimental results are further supported by theoretical calculations.

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.

An NMR crystallography study of the hemihydrate of 2', 3'-O-isopropylidineguanosine

An NMR crystallography study of the hemihydrate of 2', 3'-O-isopropylidineguanosine (Gace) is presented, together with powder X-ray diffraction and thermogravimetric analysis. (1)H double-quantum and (14)N-(1)H HMQC spectra recorded at 850MHz and 75kHz MAS (using a JEOL 1mm probe) are presented together with a (1)H-(13)C refocused INEPT spectrum recorded at 500MHz and 12.5kHz MAS using eDUMBO-122(1)H homonuclear decoupling. NMR chemical shieldings are calculated using the GIPAW (gauge-including projector augmented wave) method; good two-dimensional agreement between calculation and experiment is observed for (13)C and (1)H chemical shifts for directly bonded CH and CH3 peaks. There are two Gace molecules in the asymmetric unit cell: differences in specific (1)H chemical shifts are rationalised in terms of the strength of CH-π and intermolecular hydrogen bonding interactions.

Creatine: Polymorphs Predicted and Found

Hydrate and anhydrate crystal structure prediction (CSP) of creatine (CTN), a heavily used, badly water soluble, zwitterionic compound, has enabled the finding and characterization of its anhydrate polymorphs, including the thermodynamic room temperature form. Crystal structures of the novel forms were determined by combining laboratory powder X-ray diffraction data and ab initio generated structures. The computational method not only revealed all experimental forms but predicted the correct stability order, which was experimentally confirmed by measurements of the heat of hydration.

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 P2₁2₁2₁. 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.