Packing Preferences of Chalcones: A Model Conjugated Pharmaceutical Scaffold.

We sought the crystal packing preferences of the chalcone scaffold by analyzing 232 single-component crystal structures of chalcones with a small (six or fewer non-hydrogen atoms) substituent on either or both rings, including the unsubstituted molecule. This covers 216 molecules, as some are polymorphic, and 277 independent molecular conformations, as 16% of the crystal structures have more than one symmetry independent molecule. Quantum mechanical conformational profiles of the unsubstituted molecule and the almost 5000 crystal structures within 20 kJ mol-1 of the global minimum generated in a crystal structure prediction (CSP) study have been used to complement this analysis. Although π conjugation would be expected to favor a planar molecule, there are a significant number of crystal structures containing nonplanar molecules with an approximately 50° angle between the aromatic rings. The relative orientations of the molecules in the inversion-related dimers and translation-related dimers in the experimental crystal structures show the same trends as in the CSP-generated structures for the unsubstituted molecule, allowing for the substituent making the side-to-side distances larger. There is no type of dimer geometry associated with particularly favorable lattice energies for the chalcone core. Less than a third of the experimental structures show a face-to-face contact associated with π···π stacking. Analysis of the experimental crystal structures with XPac and Mercury finds various pairs of isostructural crystals, but the largest isostructural set had only 15 structures, with all substituents (mainly halogens) in the para position. The most common one-dimensional motif, found in half of the experimental crystal structures, is a translation-related side-to-side packing, which can be adopted by all the observed conformations. This close-packed motif can be adopted by chalcones with a particularly wide variety of substituents as the substituents are at the periphery. Thus, although the crystal structures of the substituted chalcones show thermodynamically plausible packings of the chalcone scaffold, there is little evidence for any crystal engineering principle of preferred chalcone scaffold packing beyond close packing of the specific molecule.

Visualising early-stage liquid phase organic crystal growth via liquid cell electron microscopy.

Here, we show that the development of nuclei and subsequent growth of a molecular organic crystal system can be induced by electron beam irradiation by exploiting the radiation chemistry of the carrier solvent. The technique of Liquid Cell Electron Microscopy was used to probe the crystal growth of flufenamic acid; a current commercialised active pharmaceutical ingredient. This work demonstrates liquid phase electron microscopy analysis as an essential tool for assessing pharmaceutical crystal growth in their native environment while giving insight into polymorph identification of nano-crystals at their very inception. Possible mechanisms of crystal nucleation due to the electron beam with a focus on radiolysis are discussed along with the innovations this technique offers to the study of pharmaceutical crystals and other low contrast materials.

Calculation of Diamagnetic Susceptibility Tensors of Organic Crystals: From Coronene to Pharmaceutical Polymorphs.

Understanding why crystallization in strong magnetic fields can lead to new polymorphs requires methods to calculate the diamagnetic response of organic molecular crystals. We develop the calculation of the macroscopic diamagnetic susceptibility tensor, χcryst, for organic molecular crystals using periodic density functional methods. The crystal magnetic susceptibility tensor, χcryst, for all experimentally known polymorphs, and its molecular counterpart, χmol, are calculated for flexible pharmaceuticals such as carbamazepine, flufenamic acid, and chalcones, and rigid molecules, such as benzene, pyridine, acridine, anthracene, and coronene, whose molecular magnetic properties have been traditionally studied. A tensor addition method is developed to approximate the crystal diamagnetic susceptibility tensor, χcryst, from the molecular one, χmol, giving good agreement with those calculated directly using the more costly periodic density functional method for χcryst. The response of pharmaceutical molecules and crystals to magnetic fields, as embodied by χcryst, is largely determined by the packing in the crystal, as well as the molecular conformation. The anisotropy of χcryst can vary considerably between polymorphs though the isotropic terms are fairly constant. The implications for developing a computational method for predicting whether crystallization in a magnetic field could produce a novel or different polymorph are discussed.

A Prolific Solvate Former, Galunisertib, under the Pressure of Crystal Structure Prediction, Produces Ten Diverse Polymorphs.

The solid form screening of galunisertib produced many solvates, prompting an extensive investigation into possible risks to the development of the favored monohydrate form. Inspired by crystal structure prediction, the search for neat polymorphs was expanded to an unusual range of experiments, including melt crystallization under pressure, to work around solvate formation and the thermal instability of the molecule. Ten polymorphs of galunisertib were found; however, the structure predicted to be the most stable has yet to be obtained. We present the crystal structures of all ten unsolvated polymorphs of galunisertib, showing how state-of-the-art characterization methods can be combined with emerging computational modeling techniques to produce a complete structure landscape and assess the risk of late-appearing, more stable polymorphs. The exceptional conformational polymorphism of this prolific solvate former invites further development of methods, computational and experimental, that are applicable to larger, flexible molecules with complex solid form landscapes.

Towards crystal structure prediction of complex organic compounds--a report on the fifth blind test.

Following on from the success of the previous crystal structure prediction blind tests (CSP1999, CSP2001, CSP2004 and CSP2007), a fifth such collaborative project (CSP2010) was organized at the Cambridge Crystallographic Data Centre. A range of methodologies was used by the participating groups in order to evaluate the ability of the current computational methods to predict the crystal structures of the six organic molecules chosen as targets for this blind test. The first four targets, two rigid molecules, one semi-flexible molecule and a 1:1 salt, matched the criteria for the targets from CSP2007, while the last two targets belonged to two new challenging categories - a larger, much more flexible molecule and a hydrate with more than one polymorph. Each group submitted three predictions for each target it attempted. There was at least one successful prediction for each target, and two groups were able to successfully predict the structure of the large flexible molecule as their first place submission. The results show that while not as many groups successfully predicted the structures of the three smallest molecules as in CSP2007, there is now evidence that methodologies such as dispersion-corrected density functional theory (DFT-D) are able to reliably do so. The results also highlight the many challenges posed by more complex systems and show that there are still issues to be overcome.

A strategy for producing predicted polymorphs: catemeric carbamazepine form V.

A computationally assisted approach has enabled the first catemeric polymorph of carbamazepine (form V) to be selectively formed by templating the growth of carbamazepine from the vapour phase onto the surface of a crystal of dihydrocarbamazepine form II.

Modelling organic crystal structures using distributed multipole and polarizability-based model intermolecular potentials.

Crystal structure prediction for organic molecules requires both the fast assessment of thousands to millions of crystal structures and the greatest possible accuracy in their relative energies. We describe a crystal lattice simulation program, DMACRYS, emphasizing the features that make it suitable for use in crystal structure prediction for pharmaceutical molecules using accurate anisotropic atom-atom model intermolecular potentials based on the theory of intermolecular forces. DMACRYS can optimize the lattice energy of a crystal, calculate the second derivative properties, and reduce the symmetry of the spacegroup to move away from a transition state. The calculated terahertz frequency k = 0 rigid-body lattice modes and elastic tensor can be used to estimate free energies. The program uses a distributed multipole electrostatic model (Q, t = 00,...,44s) for the electrostatic fields, and can use anisotropic atom-atom repulsion models, damped isotropic dispersion up to R(-10), as well as a range of empirically fitted isotropic exp-6 atom-atom models with different definitions of atomic types. A new feature is that an accurate model for the induction energy contribution to the lattice energy has been implemented that uses atomic anisotropic dipole polarizability models (alpha, t = (10,10)...(11c,11s)) to evaluate the changes in the molecular charge density induced by the electrostatic field within the crystal. It is demonstrated, using the four polymorphs of the pharmaceutical carbamazepine C(15)H(12)N(2)O, that whilst reproducing crystal structures is relatively easy, calculating the polymorphic energy differences to the accuracy of a few kJ mol(-1) required for applications is very demanding of assumptions made in the modelling. Thus DMACRYS enables the comparison of both known and hypothetical crystal structures as an aid to the development of pharmaceuticals and other speciality organic materials, and provides a tool to develop the modelling of the intermolecular forces involved in molecular recognition processes.