Combined crystal structure prediction and high-pressure crystallization in rational pharmaceutical polymorph screening

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

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

Exchange-Hole Dipole Dispersion Model for Accurate Energy Ranking in Molecular Crystal Structure Prediction II: Nonplanar Molecules

The crystal structure prediction (CSP) of a given compound from its molecular diagram is a fundamental challenge in computational chemistry with implications in relevant technological fields. A key component of CSP is the method to calculate the lattice energy of a crystal, which allows the ranking of candidate structures. This work is the second part of our investigation to assess the potential of the exchange-hole dipole moment (XDM) dispersion model for crystal structure prediction. In this article, we study the relatively large, nonplanar, mostly flexible molecules in the first five blind tests held by the Cambridge Crystallographic Data Centre. Four of the seven experimental structures are predicted as the energy minimum, and thermal effects are demonstrated to have a large impact on the ranking of at least another compound. As in the first part of this series, delocalization error affects the results for a single crystal (compound X), in this case by detrimentally overstabilizing the π-conjugated conformation of the monomer. Overall, B86bPBE-XDM correctly predicts 16 of the 21 compounds in the five blind tests, a result similar to the one obtained using the best CSP method available to date (dispersion-corrected PW91 by Neumann et al.). Perhaps more importantly, the systems for which B86bPBE-XDM fails to predict the experimental structure as the energy minimum are mostly the same as with Neumann's method, which suggests that similar difficulties (absence of vibrational free energy corrections, delocalization error,...) are not limited to B86bPBE-XDM but affect GGA-based DFT-methods in general. Our work confirms B86bPBE-XDM as an excellent option for crystal energy ranking in CSP and offers a guide to identify crystals (organic salts, conjugated flexible systems) where difficulties may appear.

Solvent Vapor Annealing of Amorphous Carbamazepine Films for Fast Polymorph Screening and Dissolution Alteration

Solubility enhancement and thus higher bioavailability are of great importance and a constant challenge in pharmaceutical research whereby polymorph screening and selection is one of the most important tasks. A very promising approach for polymorph screening is solvent vapor annealing where a sample is exposed to an atmosphere saturated with molecules of a specific chemical/solvent. In this work, amorphous carbamazepine thin films were prepared by spin coating, and the transformation into crystalline forms under exposure to solvent vapors was investigated. Employing grazing incidence X-ray diffraction, four distinct carbamazepine polymorphs, a solvate, and hydrates could be identified, while optical microscopy showed mainly spherulitic morphologies. In vitro dissolution experiments revealed different carbamazepine release from the various thin-film samples containing distinct polymorphic compositions: heat treatment of amorphous samples at 80 °C results in an immediate release; samples exposed to EtOH vapors show a drug release about 5 times slower than this immediate one; and all the others had intermediate release profiles. Noteworthy, even the sample of slowest release has a manifold faster release compared to a standard powder sample demonstrating the capabilities of thin-film preparation for faster drug release in general. Despite the small number of samples in this screening experiment, the results clearly show how solvent vapor annealing can assist in identifying potential polymorphs and allows for estimating their impact on properties like bioavailability.

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

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

Combined Chemoinformatics Approach to Solvent Library Design Using clusterSim and Multidimensional Scaling

Reported here is a rational approach for the selection of solvents intended for use in physical form screening based on a novel chemoinformatics analysis of solvent properties. A comprehensive assessment of eight clustering methods was carried out on a series of 94 solvents described by calculated molecular descriptors using the clusterSim package in R. The effectiveness of clustering methods was evaluated using a range of statistical measures as well as increasing efficiency of solid form discovery using a cluster-based solvent selection approach. Multidimensional scaling was used to illustrate cluster analysis on a two-dimensional solvent map. The map presented here is a valuable tool to aid efficient solvent selection in physical form screens. This tool is equally applicable to any scientific area which requires a solubility dependent decision on solvent choice.

Antitumoural Sulphur and Selenium Heteroaryl Compounds: Thermal Characterization and Stability Evaluation

The physicochemical properties of a compound play a crucial role in the cancer development process. In this context, polymorphism can become an important obstacle for the pharmaceutical industry because it frequently leads to the loss of therapeutic effectiveness of some drugs. Stability under manufacturing conditions is also critical to ensure no undesired degradations or transformations occur. In this study, the thermal behaviour of 40 derivatives of a series of sulphur and selenium heteroaryl compounds with potential antitumoural activity were studied. In addition, the most promising cytotoxic derivatives were analysed by a combination of differential scanning calorimetry, X-ray diffraction and thermogravimetric techniques in order to investigate their polymorphism and thermal stability. Moreover, stability under acid, alkaline and oxidative media was tested. Degradation under stress conditions as well as the presence of polymorphism was found for the compounds VA6E and VA7J, which might present a hurdle to carrying on with formulation. On the contrary, these obstacles were not found for derivative VA4J.

Pharmaceutical cocrystals, salts and polymorphs: Advanced characterization techniques

The main goal of a novel drug development is to obtain it with optimal physiochemical, pharmaceutical and biological properties. Pharmaceutical companies and scientists modify active pharmaceutical ingredients (APIs), which often are cocrystals, salts or carefully selected polymorphs, to improve the properties of a parent drug. To find the best form of a drug, various advanced characterization methods should be used. In this review, we have described such analytical methods, dedicated to solid drug forms. Thus, diffraction, spectroscopic, thermal and also pharmaceutical characterization methods are discussed. They all are necessary to study a solid API in its intrinsic complexity from bulk down to the molecular level, gain information on its structure, properties, purity and possible transformations, and make the characterization efficient, comprehensive and complete. Furthermore, these methods can be used to monitor and investigate physical processes, involved in the drug development, in situ and in real time. The main aim of this paper is to gather information on the current advancements in the analytical methods and highlight their pharmaceutical relevance.

Dense Semi-Clathrates at High Pressure: A Study of the Water-tert-Butylamine System

In situ high-pressure crystallization and diffraction techniques have been applied to obtain two very structurally distinct semi-clathrates of the tert-butylamine-water system with hydration numbers 5.65 and 5.8, respectively, thereby considerably reducing a notable hydration gap between the monohydrate and the 71/4 -hydrate that results when crystallization space is explored by temperature alone. Both structures can be considered as an intriguing solid-state example of hydrophobic hydration, in which the water network creates wide tert-butylamine-filled channels stabilized by cross-linking hydrogen bonds. The existence of interconnected channels might also add low hydration structures to a list of potential targets for hydrogen storage. A detailed analysis of the topology of host water and host-guest interactions is reported and extended to those of other hydrates of the compound. This analysis offers new insight into properties of the tert-butylamine-water system and provides some clues as to the occurrence of the sizable number of hydrates of this compound.

The application of tailor-made force fields and molecular dynamics for NMR crystallography: a case study of free base cocaine

Motional averaging has been proven to be significant in predicting the chemical shifts in ab initio solid-state NMR calculations, and the applicability of motional averaging with molecular dynamics has been shown to depend on the accuracy of the molecular mechanical force field. The performance of a fully automatically generated tailor-made force field (TMFF) for the dynamic aspects of NMR crystallography is evaluated and compared with existing benchmarks, including static dispersion-corrected density functional theory calculations and the COMPASS force field. The crystal structure of free base cocaine is used as an example. The results reveal that, even though the TMFF outperforms the COMPASS force field for representing the energies and conformations of predicted structures, it does not give significant improvement in the accuracy of NMR calculations. Further studies should direct more attention to anisotropic chemical shifts and development of the method of solid-state NMR calculations.

Diabat Interpolation for Polymorph Free-Energy Differences

Existing methods to compute free-energy differences between polymorphs use harmonic approximations, advanced non-Boltzmann bias sampling techniques, and/or multistage free-energy perturbations. This work demonstrates how Bennett's diabat interpolation method ( J. Comput. Phys. 1976, 22, 245 ) can be combined with energy gaps from lattice-switch Monte Carlo techniques ( Phys. Rev. E 2000, 61, 906 ) to swiftly estimate polymorph free-energy differences. The new method requires only two unbiased molecular dynamics simulations, one for each polymorph. To illustrate the new method, we compute the free-energy difference between face-centered cubic and body-centered cubic polymorphs for a Gaussian core solid. We discuss the justification for parabolic models of the free-energy diabats and similarities to methods that have been used in studies of electron transfer.

Pressure-driven phase transition mechanisms revealed by quantum chemistry: l-serine polymorphs

The present study delivers a computational approach for the understanding of the mechanism of phase transitions between polymorphs of small organic molecules.

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.

Predicting Molecular Crystal Properties from First Principles: Finite-Temperature Thermochemistry to NMR Crystallography

Molecular crystals occur widely in pharmaceuticals, foods, explosives, organic semiconductors, and many other applications. Thanks to substantial progress in electronic structure modeling of molecular crystals, attention is now shifting from basic crystal structure prediction and lattice energy modeling toward the accurate prediction of experimentally observable properties at finite temperatures and pressures. This Account discusses how fragment-based electronic structure methods can be used to model a variety of experimentally relevant molecular crystal properties. First, it describes the coupling of fragment electronic structure models with quasi-harmonic techniques for modeling the thermal expansion of molecular crystals, and what effects this expansion has on thermochemical and mechanical properties. Excellent agreement with experiment is demonstrated for the molar volume, sublimation enthalpy, entropy, and free energy, and the bulk modulus of phase I carbon dioxide when large basis second-order Møller-Plesset perturbation theory (MP2) or coupled cluster theories (CCSD(T)) are used. In addition, physical insight is offered into how neglect of thermal expansion affects these properties. Zero-point vibrational motion leads to an appreciable expansion in the molar volume; in carbon dioxide, it accounts for around 30% of the overall volume expansion between the electronic structure energy minimum and the molar volume at the sublimation point. In addition, because thermal expansion typically weakens the intermolecular interactions, neglecting thermal expansion artificially stabilizes the solid and causes the sublimation enthalpy to be too large at higher temperatures. Thermal expansion also frequently weakens the lower-frequency lattice phonon modes; neglecting thermal expansion causes the entropy of sublimation to be overestimated. Interestingly, the sublimation free energy is less significantly affected by neglecting thermal expansion because the systematic errors in the enthalpy and entropy cancel somewhat. Second, because solid state nuclear magnetic resonance (NMR) plays an increasingly important role in molecular crystal studies, this Account discusses how fragment methods can be used to achieve higher-accuracy chemical shifts in molecular crystals. Whereas widely used plane wave density functional theory models are largely restricted to generalized gradient approximation (GGA) functionals like PBE in practice, fragment methods allow the routine use of hybrid density functionals with only modest increases in computational cost. In extensive molecular crystal benchmarks, hybrid functionals like PBE0 predict chemical shifts with 20-30% higher accuracy than GGAs, particularly for ¹H, ¹³C, and ¹⁵N nuclei. Due to their higher sensitivity to polarization effects, ¹⁷O chemical shifts prove slightly harder to predict with fragment methods. Nevertheless, the fragment model results are still competitive with those from GIPAW. The improved accuracy achievable with fragment approaches and hybrid density functionals increases discrimination between different potential assignments of individual shifts or crystal structures, which is critical in NMR crystallography applications. This higher accuracy and greater discrimination are highlighted in application to the solid state NMR of different acetaminophen and testosterone crystal forms.

Coarse-grained modeling of crystal growth and polymorphism of a model pharmaceutical molecule

We describe a systematic coarse-graining method to study crystallization and predict possible polymorphs of small organic molecules. In this method, a coarse-grained (CG) force field is obtained by inverse-Boltzmann iteration from the radial distribution function of atomistic simulations of the known crystal. With the force field obtained by this method, we show that CG simulations of the drug phenytoin predict growth of a crystalline slab from a melt of phenytoin, allowing determination of the fastest-growing surface, as well as giving the correct lattice parameters and crystal morphology. By applying meta-dynamics to the coarse-grained model, a new crystalline form of phenytoin (monoclinic, space group P2₁) was predicted which is different from the experimentally known crystal structure (orthorhombic, space group Pna2₁). Atomistic simulations and quantum calculations then showed the polymorph to be meta-stable at ambient temperature and pressure, and thermodynamically more stable than the conventional orthorhombic crystal at high pressure. The results suggest an efficient route to study crystal growth of small organic molecules that could also be useful for identification of possible polymorphs as well.

Ab initio structure determination from prion nanocrystals at atomic resolution by MicroED

Electrons, because of their strong interaction with matter, produce high-resolution diffraction patterns from tiny 3D crystals only a few hundred nanometers thick in a frozen-hydrated state. This discovery offers the prospect of facile structure determination of complex biological macromolecules, which cannot be coaxed to form crystals large enough for conventional crystallography or cannot easily be produced in sufficient quantities. Two potential obstacles stand in the way. The first is a phenomenon known as dynamical scattering, in which multiple scattering events scramble the recorded electron diffraction intensities so that they are no longer informative of the crystallized molecule. The second obstacle is the lack of a proven means of de novo phase determination, as is required if the molecule crystallized is insufficiently similar to one that has been previously determined. We show with four structures of the amyloid core of the Sup35 prion protein that, if the diffraction resolution is high enough, sufficiently accurate phases can be obtained by direct methods with the cryo-EM method microelectron diffraction (MicroED), just as in X-ray diffraction. The success of these four experiments dispels the concern that dynamical scattering is an obstacle to ab initio phasing by MicroED and suggests that structures of novel macromolecules can also be determined by direct methods.

Report on the sixth blind test of organic crystal structure prediction methods

The sixth blind test of organic crystal structure prediction (CSP) methods has been held, with five target systems: a small nearly rigid molecule, a polymorphic former drug candidate, a chloride salt hydrate, a co-crystal and a bulky flexible molecule. This blind test has seen substantial growth in the number of participants, with the broad range of prediction methods giving a unique insight into the state of the art in the field. Significant progress has been seen in treating flexible molecules, usage of hierarchical approaches to ranking structures, the application of density-functional approximations, and the establishment of new workflows and `best practices' for performing CSP calculations. All of the targets, apart from a single potentially disordered Z' = 2 polymorph of the drug candidate, were predicted by at least one submission. Despite many remaining challenges, it is clear that CSP methods are becoming more applicable to a wider range of real systems, including salts, hydrates and larger flexible molecules. The results also highlight the potential for CSP calculations to complement and augment experimental studies of organic solid forms.

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.

Can computed crystal energy landscapes help understand pharmaceutical solids?

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

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

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

Recent developments in the experimental investigations of relaxations in pharmaceuticals by dielectric techniques at ambient and elevated pressure

In recent years, there is a growing interest in improving the physicochemical stability of amorphous pharmaceutical solids due to their very promising applications to manufacture medicines characterized by a better water solubility, and consequently by a higher dissolution rate than those of their crystalline counterparts. In this review article, we show that the molecular mobility investigated both in the supercooled liquid and glassy states is the crucial factor required to understand molecular mechanisms that govern the physical stability of amorphous drugs. We demonstrate that pharmaceuticals can be thoroughly examined by means of the broadband dielectric spectroscopy, which is a very useful experimental technique to explore different relaxation processes and crystallization kinetics as well. Such studies conducted in the wide temperature and pressure ranges provide data needed in searching correlations between properties of molecular dynamics and crystallization process, which are aimed at developing effective and efficient methods for stabilizing amorphous drugs.

Modeling Polymorphic Molecular Crystals with Electronic Structure Theory

Interest in molecular crystals has grown thanks to their relevance to pharmaceuticals, organic semiconductor materials, foods, and many other applications. Electronic structure methods have become an increasingly important tool for modeling molecular crystals and polymorphism. This article reviews electronic structure techniques used to model molecular crystals, including periodic density functional theory, periodic second-order Møller-Plesset perturbation theory, fragment-based electronic structure methods, and diffusion Monte Carlo. It also discusses the use of these models for predicting a variety of crystal properties that are relevant to the study of polymorphism, including lattice energies, structures, crystal structure prediction, polymorphism, phase diagrams, vibrational spectroscopies, and nuclear magnetic resonance spectroscopy. Finally, tools for analyzing crystal structures and intermolecular interactions are briefly discussed.

Convergence Properties of Crystal Structure Prediction by Quasi-Random Sampling

Generating sets of trial structures that sample the configurational space of crystal packing possibilities is an essential step in the process of ab initio crystal structure prediction (CSP). One effective methodology for performing such a search relies on low-discrepancy, quasi-random sampling, and our implementation of such a search for molecular crystals is described in this paper. Herein we restrict ourselves to rigid organic molecules and, by considering their geometric properties, build trial crystal packings as starting points for local lattice energy minimization. We also describe a method to match instances of the same structure, which we use to measure the convergence of our packing search toward completeness. The use of these tools is demonstrated for a set of molecules with diverse molecular characteristics and as representative of areas of application where CSP has been applied. An important finding is that the lowest energy crystal structures are typically located early and frequently during a quasi-random search of phase space. It is usually the complete sampling of higher energy structures that requires extended sampling. We show how the procedure can first be refined, through targetting the volume of the generated crystal structures, and then extended across a range of space groups to make a full CSP search and locate experimentally observed and lists of hypothetical polymorphs. As the described method has also been created to lie at the base of more involved approaches to CSP, which are being developed within the Global Lattice Energy Explorer (Glee) software, a few of these extensions are briefly discussed.

The role of fluids in high-pressure polymorphism of drugs: different behaviour of β-chlorpropamide in different inert gas and liquid media

Compression of β-chlorpropamide gives different phases depending on the choice of non-dissolving pressure-transmitting fluid (paraffin, neon and helium).

Giant strain geared to transformable H-bonded network in compressed β-d-mannose

The clearly visible shape change of β-d-mannose crystals is associated with pressure-induced isostructural phase transition.

Pharmaceutical crystallization with nanocellulose organogels

A carboxylated nanocellulose organogel facilitates crystallization of pharmaceutical solid forms including a new sulfapyridine solvate.