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

Molecular Insights into Water-Chloride and Water-Water Interactions in the Supramolecular Architecture of Aconine Hydrochloride Dihydrate

Despite the previous preparation of aconine hydrochloride monohydrate (AHM), accurate determination of the crystal's composition was hindered by severely disordered water molecules within the crystal. In this study, we successfully prepared a new dihydrate form of the aconine hydrochloride [C25H42NO9+Cl-·2(H2O), aconine hydrochloride dihydrate (AHD)] and accurately refined all water molecules within the AHD crystal. Our objective is to elucidate both water-chloride and water-water interactions in the AHD crystal. The crystal structure of AHD was determined at 136 K using X-ray diffraction and a multipolar atom model was constructed by transferring charge-density parameters to explore the topological features of key short contacts. By comparing the crystal structures of dihydrate and monohydrate forms, we have observed that both AHD and AHM exhibit identical aconine cations, except for variations in the number of water molecules present. In the AHD crystal, chloride anions and water molecules serve as pivotal connecting hubs to establish three-dimensional hydrogen bonding networks and one-dimensional hydrogen bonding chain; both water-chloride and water-water interactions assemble supramolecular architectures. The crystal packing of AHD exhibits a complete reversal in the stacking order compared to AHM, thereby emphasizing distinct disparities between them. Hirshfeld surface analysis reveals that H···Cl- and H···O contacts play a significant role in constructing the hydrogen bonding network and chain within these supramolecular architectures. Furthermore, topological analysis and electrostatic interaction energy confirm that both water-chloride and water-water interactions stabilize supramolecular architectures through electrostatic attraction facilitated by H···Cl- and H···O contacts. Importantly, these findings are strongly supported by the existing literature evidence. Consequently, navigating these water-chloride and water-water interactions is imperative for ensuring storage and safe processing of this pharmaceutical compound.

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.

Navigating the Complex Solid Form Landscape of the Quercetin Flavonoid Molecule

Quercetin, a naturally occurring bioflavonoid substance widely used in the nutraceutical and food industries, exists in various solid forms that can have different physicochemical properties, thus impacting this compound's performance in various applications. In this work, we will clarify the complex solid-form landscape of this molecule. Two elusive isostructural solvates of quercetin were obtained from ethanol and methanol. The obtained crystals were characterized experimentally, but the crystallographic structure could not be solved due to their high instability. Nevertheless, the desolvated structure resulting from a high-temperature treatment (or prolonged storage at ambient conditions) of both these two labile crystals was characterized and solved via powder X-ray diffraction and solid-state nuclear magnetic resonance (SSNMR). This anhydrous crystal structure was compared with another anhydrous quercetin form obtained in our previous work, indicating that, at least, two different anhydrous polymorphs of quercetin exist. Navigating the solid-form landscape of quercetin is essential to ensure accurate control of the functional properties of food, nutraceutical, or pharmaceutical products containing crystal forms of this substance.

A theoretical framework for the design of molecular crystal engines

Photomechanical molecular crystals have garnered attention for their ability to transform light into mechanical work, but difficulties in characterizing the structural changes and mechanical responses experimentally have hindered the development of practical organic crystal engines. This study proposes a new computational framework for predicting the solid-state crystal-to-crystal photochemical transformations entirely from first principles, and it establishes a photomechanical engine cycle that quantifies the anisotropic mechanical performance resulting from the transformation. The approach relies on crystal structure prediction, solid-state topochemical principles, and high-quality electronic structure methods. After validating the framework on the well-studied [4 + 4] cycloadditions in 9-methyl anthracene and 9-tert-butyl anthracene ester, the experimentally-unknown solid-state transformation of 9-carboxylic acid anthracene is predicted for the first time. The results illustrate how the mechanical work is done by relaxation of the crystal lattice to accommodate the photoproduct, rather than by the photochemistry itself. The large ∼10⁷ J m^-3 work densities computed for all three systems highlight the promise of photomechanical crystal engines. This study demonstrates the importance of crystal packing in determining molecular crystal engine performance and provides tools and insights to design improved materials in silico.

Efficient Crystal Structure Prediction for Structurally Related Molecules with Accurate and Transferable Tailor-Made Force Fields

Crystal structure prediction (CSP) his generally used to complement experimental solid form screening and applied to individual molecules in drug development. The fast development of algorithms and computing resources offers the opportunity to use CSP earlier and for a broader range of applications in the drug design cycle. This study presents a novel paradigm of CSP specifically designed for structurally related molecules, referred to as Quick-CSP. The approach prioritizes more accurate physics through robust and transferable tailor-made force fields (TMFFs), such that significant efficiency gains are achieved through the reduction of expensive ab initio calculations. The accuracy of the TMFF is increased by the introduction of electrostatic multipoles, and the fragment-based force field parameterization scheme is demonstrated to be transferable for a family of chemically related molecules. The protocol is benchmarked with structurally related compounds from the Bromodomain and Extraterminal (BET) domain inhibitors series. A new convergence criterion is introduced that aims at performing only as many ab initio optimizations of crystal structures as required to locate the bottom of the crystal energy landscape within a user-defined accuracy. The overall approach provides significant cost savings ranging from three- to eight-fold less than the full-CSP workflow. The reported advancements expand the scope and utility of the underlying CSP building blocks as well as their novel reassembly to other applications earlier in the drug design cycle to guide molecule design and selection.

Polymorphism in Ionic Cocrystals Comprising Lithium Salts and l-Proline

The occurrence of polymorphism in ionic cocrystals formed by two lithium salts, lithium salicylate (LIS) and lithium 4-methoxybenzoate (L4M), and l-proline (PRO) has been investigated. The previously reported monoclinic form of the 1:1 cocrystal of LIS and PRO, LISPRO(α), and a new thermodynamically stable orthorhombic polymorph, LISPRO(β), were prepared and characterized. The two polymorphs form square grid, sql, topology coordination networks and differ mainly in the conformation of the salicylate ions and positioning of the sql nets. LISPRO(α) was observed to transform to LISPRO(β) under slurry conditions. The 1:1 ionic cocrystal of L4M and PRO (L4MPRO) was found to form three polymorphs. Apart from the previously reported orthorhombic crystal form, L4MPRO(α), two new monoclinic crystal forms, L4MPRO(β) and L4MPRO(γ), were obtained by modifying crystallization conditions. The new polymorphs were found to be metastable, undergoing transformations to L4MPRO(α) upon exposure to humidity. Experimental conditions that induce transformations between the polymorphs of LISPRO and L4MPRO are detailed, and the structural differences between the polymorphs are discussed in the broader context of polymorphism.

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.

How many more polymorphs of ROY remain undiscovered

With 12 crystal forms, 5-methyl-2-[(2-nitrophenyl)amino]-3-thiophenecabonitrile (a.k.a. ROY) holds the current record for the largest number of fully characterized organic crystal polymorphs. Four of these polymorph structures have been reported since 2019, raising the question of how many more ROY polymorphs await future discovery. Employing crystal structure prediction and accurate energy rankings derived from conformational energy-corrected density functional theory, this study presents the first crystal energy landscape for ROY that agrees well with experiment. The lattice energies suggest that the seven most stable ROY polymorphs (and nine of the twelve lowest-energy forms) on the Z' = 1 landscape have already been discovered experimentally. Discovering any new polymorphs at ambient pressure will likely require specialized crystallization techniques capable of trapping metastable forms. At pressures above 10 GPa, however, a new crystal form is predicted to become enthalpically more stable than all known polymorphs, suggesting that further high-pressure experiments on ROY may be warranted. This work highlights the value of high-accuracy crystal structure prediction for solid-form screening and demonstrates how pragmatic conformational energy corrections can overcome the limitations of conventional density functionals for conformational polymorphs.

Modeling the α- and β-resorcinol phase boundary via combination of density functional theory and density functional tight-binding

The ability to predict not only what organic crystal structures might occur but also the thermodynamic conditions under which they are the most stable would be extremely useful for discovering and designing new organic materials. The present study takes a step in that direction by predicting the temperature- and pressure-dependent phase boundary between the α and β polymorphs of resorcinol using density functional theory (DFT) and the quasi-harmonic approximation. To circumvent the major computational bottleneck associated with computing a well-converged phonon density of states via the supercell approach, a recently developed approximation is employed, which combines a supercell phonon density of states from dispersion-corrected third-order density functional tight binding [DFTB3-D3(BJ)] with frequency corrections derived from a smaller B86bPBE-XDM functional DFT phonon calculation on the crystallographic unit cell. This mixed DFT/DFTB quasi-harmonic approach predicts the lattice constants and unit cell volumes to within 1%-2% at lower pressures. It predicts the thermodynamic phase boundary in almost perfect agreement with the experiment, although this excellent agreement does reflect fortuitous cancellation of errors between the enthalpy and entropy of transition.

The Eight Hydrates of Strychnine Sulfate

Commercial samples of strychnine sulfate were used as the starting material in crystallization experiments accompanied by stability studies. Eight hydrate forms (HyA-HyG), including five novel hydrates, were verified. The crystal structures of HyA ("pentahydrate") and HyF ("hexahydrate") were determined from single-crystal X-ray diffraction data. HyF was identified as the most stable hydrate at high water activities at room temperature (RT), and HyA and HyC were also found to be stable at ambient conditions. Long-time storage experiments over nearly two decades confirm that these three hydrates are stable at ambient conditions (20-60% relative humidity). The other five hydrates, HyB ("dihydrate"), HyD, HyE, HyG, and HyH, are only observable at the low(est) relative humidity (RH) levels at RT. Some of these latter forms can only exist within a very narrow RH range and are therefore intermediate phases. By applying a range of complementary experimental techniques such as gravimetric moisture sorption analysis, thermal analysis, moisture controlled PXRD measurements, and variable temperature IR spectroscopy in combination with principal component analysis, it was possible to identify the distinct hydrate phases and elucidate their stability and dehydration pathways. The observed (de)hydration routes, HyA ↔ HyB, HyC ↔ HyD ↔ HyE, HyF ↔ HyG ↔ HyH and HyF → HyA ↔ HyB, depended on the initial hydrate form, particle size, and atmospheric conditions. In addition, a transformation from HyC/HyA to HyF occurs at high RH values at RT. The specific moisture and temperature conditions of none of the applied drying regimes yielded a crystalline water-free form, which highlights the essential role of water molecules for the formation and stability of the crystalline strychnine sulfate phases.

Overcoming the difficulties of predicting conformational polymorph energetics in molecular crystals via correlated wavefunction methods

Molecular crystal structure prediction is increasingly being applied to study the solid form landscapes of larger, more flexible pharmaceutical molecules. Despite many successes in crystal structure prediction, van der Waals-inclusive density functional theory (DFT) methods exhibit serious failures predicting the polymorph stabilities for a number of systems exhibiting conformational polymorphism, where changes in intramolecular conformation lead to different intermolecular crystal packings. Here, the stabilities of the conformational polymorphs of o-acetamidobenzamide, ROY, and oxalyl dihydrazide are examined in detail. DFT functionals that have previously been very successful in crystal structure prediction perform poorly in all three systems, due primarily to the poor intramolecular conformational energies, but also due to the intermolecular description in oxalyl dihydrazide. In all three cases, a fragment-based dispersion-corrected second-order Møller-Plesset perturbation theory (MP2D) treatment of the crystals overcomes these difficulties and predicts conformational polymorph stabilities in good agreement with experiment. These results highlight the need for methods which go beyond current-generation DFT functionals to make crystal polymorph stability predictions truly reliable.

Same or different – that is the question: identification of crystal forms from crystal structure data

An analysis of the CSD with structural comparison tools shows that differentiating between polymorphism and redeterminations is not always straight forward and requires of complementary tools at the hands of an expert practitioner.

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.

The solid state forms of the sex hormone 17-β-estradiol

The crystal structure of the female sex hormone has been established despite its high affinity for water.

Analysis of stacking disorder in ice I using pair distribution functions

Stacking-disordered materials display crystalline order in two dimensions but are disordered along the direction in which layered structural motifs are stacked. Countless examples of stacking disorder exist, ranging from close-packed metals, ice I and diamond to open-framework materials and small-molecule pharmaceuticals. In general, the presence of stacking disorder can have profound consequences for the physical and chemical properties of a material. Traditional analyses of powder diffraction data are often complicated by the presence of memory effects in the stacking sequences. Here it is shown that experimental pair distribution functions of stacking-disordered ice I can be used to determine local information on the fractions of cubic and hexagonal stacking. Ice is a particularly challenging material in this respect, since both the stacking disorder and the orientational disorder of the water molecules need to be described. Memory effects are found to contribute very little to the pair distribution functions, and consequently, the analysis of pair distribution functions is the method of choice for characterizing stacking-disordered samples with complicated and high-order memory effects. In the context of this work, the limitations of current structure-reconstruction approaches are also discussed.

Crystal structure prediction of flexible pharmaceutical-like molecules: density functional tight-binding as an intermediate optimisation method and for free energy estimation

Successful methodologies for theoretical crystal structure prediction (CSP) on flexible pharmaceutical-like organic molecules explore the lattice energy surface to find a set of plausible crystal structures. The initial search stages of CSP studies use relatively simple lattice energy approximations as hundreds of thousands of minima have to be considered. These generated crystal structures often have poor molecular geometries, as well as inaccurate lattice energy rankings, and performing reasonably accurate but computationally affordable optimisations of the crystal structures generated in a search would be highly desirable. Here, we seek to explore whether semi-empirical quantum-mechanical methods can perform this task. We employed the dispersion-corrected tight-binding Hamiltonian (DFTB3-D3) to relax all the inter- and intra-molecular degrees of freedom of several thousands of generated crystal structures of five pharmaceutical-like molecules, saving a large amount of computational effort compared to earlier studies. The computational cost scales better with molecular size and flexibility than other CSP methods, suggesting that it could be extended to even larger and more flexible molecules. On average, this optimisation improved the average reproduction of the eight experimental crystal structures (RMSD15) and experimental conformers (RMSD1) by 4% and 23%, respectively. The intermolecular interactions were then further optimised using distributed multipoles, derived from the molecular wave-functions, to accurately describe the electrostatic components of the intermolecular energies. In all cases, the experimental crystal structures are close to the top of the lattice energy ranking. Phonon calculations on some of the lowest energy structures were also performed with DFTB3-D3 methods to calculate the vibrational component of the Helmholtz free energy, providing further insights into the solid-state behaviour of the target molecules. We conclude that DFTB3-D3 is a cost-effective method for optimising flexible molecules, bridging the gap between the approximate methods used in CSP searches for generating crystal structures and more accurate methods required in the final energy ranking.

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