Structure determination from powder data without prior indexing, using a similarity measure based on cross-correlation functions

Trace polymer coated clarithromycin spherulites: Formation mechanism, improvement in pharmaceutical properties and development of high-drug-loading direct compression tablets

Clarithromycin (CLA) is a high dose antibiotic drug exhibiting poor flowability and tabletability, making the tablet development challenging. This study aims to develop spherulitic CLA by introducing trace amount of polymer in crystallization solution. Its formation mechanism, physicochemical properties and potential for the direct compression (DC) tablets development were also investigated. Morphological analyses and the in situ observation on crystallization process revealed that the CLA spherulites are formed by fractal branching growth from both sides of the threadlike precursor fibers. 1H NMR analysis and nucleation time monitoring indicated that the existence of hydroxypropyl cellulose in solution slowed down the crystal nucleation and growth rate by forming hydrogen bonding interactions with CLA molecules, making the system maintain high supersaturation, providing high driving forces for CLA spherulitic growth. In comparison to commercial CLA, the CLA spherulites exhibit profoundly improved flowability, tabletability and dissolution behaviors. XPS, contact angle and Raman mapping analysis confirmed the presence of a thin HPC layer on the surfaces and interior of CLA spherulitic particles, resulting in increasing powder plasticity, interparticulate bonding strength and powder wettability, thus better tabletability and dissolution performances. The improved flowability and tabletability of CLA spherulites also enabled the successful development of DC tablet formulation with a high CLA loading (82.8 wt%) and similar dissolution profiles to reference listed drug. This study provides a novel solid form of CLA with superior manufacturability for further development.

Leucopterin, the white pigment in butterfly wings: structural analysis by PDF fit, FIDEL fit, Rietveld refinement, solid-state NMR and DFT-D

Leucopterin (C6H5N5O3) is the white pigment in the wings of Pieris brassicae butterflies, and other butterflies; it can also be found in wasps and other insects. Its crystal structure and its tautomeric form in the solid state were hitherto unknown. Leucopterin turned out to be a variable hydrate, with 0.5 to about 0.1 molecules of water per leucopterin molecule. Under ambient conditions, the preferred state is the hemihydrate. Initially, all attempts to grow single crystals suitable for X-ray diffraction were to no avail. Attempts to determine the crystal structure by powder diffraction using the direct-space method failed, because the trials did not include the correct, but rare, space group P2/c. Attempts were made to solve the crystal structure by a global fit to the pair distribution function (PDF-Global-Fit), as described by Prill and co-workers [Schlesinger et al. (2021). J. Appl. Cryst. 54, 776-786]. The approach worked well, but the correct structure was not found, because again the correct space group was not included. Finally, tiny single crystals of the hemihydrate could be obtained, which allowed at least the determination of the crystal symmetry and the positions of the C, N and O atoms. The tautomeric state of the hemihydrate was assessed by multinuclear solid-state NMR spectroscopy. 15N CPMAS spectra showed the presence of one NH2 and three NH groups, and one unprotonated N atom, which agreed with the 1H MAS and 13C CPMAS spectra. Independently, the tautomeric state was investigated by lattice-energy minimizations with dispersion-corrected density functional theory (DFT-D) on 17 different possible tautomers, which also included the prediction of the corresponding 1H, 13C and 15N chemical shifts in the solid. All methods showed the presence of the 2-amino-3,5,8-H tautomer. The DFT-D calculations also confirmed the crystal structure. Heating of the hemihydrate results in a slow release of water between 130 and 250 °C, as shown by differential thermal analysis and thermogravimetry (DTA-TG). Temperature-dependent powder X-ray diffraction (PXRD) showed an irreversible continuous shift of the reflections upon heating, which reveals that leucopterin is a variable hydrate. This observation was also confirmed by PXRD of samples obtained under various synthetic and drying conditions. The crystal structure of a sample with about 0.2 molecules of water per leucopterin was solved by a fit with deviating lattice parameters (FIDEL), as described by Habermehl et al. [Acta Cryst. (2022), B78, 195-213]. A local fit, starting from the structure of the hemihydrate, as well as a global fit, starting from random structures, were performed, followed by Rietveld refinements. Despite dehydration, the space group remains P2/c. In both structures (hemihydrate and variable hydrate), the leucopterin molecules are connected by 2-4 hydrogen bonds into chains, which are connected by further hydrogen bonds to neighbouring chains. The molecular packing is very efficient. The density of leucopterin hemihydrate is as high as 1.909 kg dm-3, which is one of the highest densities for organic compounds consisting of C, H, N and O only. The high density might explain the good light-scattering and opacity properties of the wings of Pieris brassicae and other butterflies.

From Powders to Single Crystals: A Crystallographer's Toolbox for Small-Molecule Structure Determination

Although the crystal structures of small-molecule compounds are often determined from single-crystal X-ray diffraction (scXRD), recent advances in three-dimensional electron diffraction (3DED) and crystal structure prediction (CSP) methods promise to expand the structure elucidation toolbox available to the crystallographer. Herein, a comparative assessment of scXRD, 3DED, and CSP in combination with powder X-ray diffraction is carried out on two former drug candidate compounds and a multicomponent crystal of a key building block in the synthesis of gefapixant citrate.

Ambiguous structure determination from powder data: four different structural models of 4,11-di-fluoro-quinacridone with similar X-ray powder patterns, fit to the PDF, SSNMR and DFT-D

Four different structural models, which all fit the same X-ray powder pattern, were obtained in the structure determination of 4,11-di-fluoro-quinacridone (C20H10N2O2F2) from unindexed X-ray powder data by a global fit. The models differ in their lattice parameters, space groups, Z, Z', molecular packing and hydrogen bond patterns. The molecules form a criss-cross pattern in models A and B, a layer structure built from chains in model C and a criss-cross arrangement of dimers in model D. Nevertheless, all models give a good Rietveld fit to the experimental powder pattern with acceptable R-values. All molecular geometries are reliable, except for model D, which is slightly distorted. All structures are crystallochemically plausible, concerning density, hydrogen bonds, intermolecular distances etc. All models passed the checkCIF test without major problems; only in model A a missed symmetry was detected. All structures could have probably been published, although 3 of the 4 structures were wrong. The investigation, which of the four structures is actually the correct one, was challenging. Six methods were used: (1) Rietveld refinements, (2) fit of the crystal structures to the pair distribution function (PDF) including the refinement of lattice parameters and atomic coordinates, (3) evaluation of the colour, (4) lattice-energy minimizations with force fields, (5) lattice-energy minimizations by two dispersion-corrected density functional theory methods, and (6) multinuclear CPMAS solid-state NMR spectroscopy (1H, 13C, 19F) including the comparison of calculated and experimental chemical shifts. All in all, model B (perhaps with some disorder) can probably be considered to be the correct one. This work shows that a structure determination from limited-quality powder data may result in totally different structural models, which all may be correct or wrong, even if they are chemically sensible and give a good Rietveld refinement. Additionally, the work is an excellent example that the refinement of an organic crystal structure can be successfully performed by a fit to the PDF, and the combination of computed and experimental solid-state NMR chemical shifts can provide further information for the selection of the most reliable structure among several possibilities.

Structure determination from unindexed powder data from scratch by a global optimization approach using pattern comparison based on cross-correlation functions

A method of ab initio crystal structure determination from powder diffraction data for organic and metal-organic compounds, which does not require prior indexing of the powder pattern, has been developed. Only a reasonable molecular geometry is required, needing knowledge of neither unit-cell parameters nor space group. The structures are solved from scratch by a global fit to the powder data using the new program FIDEL-GO (`FIt with DEviating Lattice parameters - Global Optimization'). FIDEL-GO uses a similarity measure based on cross-correlation functions, which allows the comparison of simulated and experimental powder data even if the unit-cell parameters deviate strongly. The optimization starts from large sets of random structures in various space groups. The unit-cell parameters, molecular position and orientation, and selected internal degrees of freedom are fitted simultaneously to the powder pattern. The optimization proceeds in an elaborate multi-step procedure with built-in clustering of duplicate structures and iterative adaptation of parameter ranges. The best structures are selected for an automatic Rietveld refinement. Finally, a user-controlled Rietveld refinement is performed. The procedure aims for the analysis of a wide range of `problematic' powder patterns, in particular powders of low crystallinity. The method can also be used for the clustering and screening of a large number of possible structure candidates and other application scenarios. Examples are presented for structure determination from unindexed powder data of the previously unknown structures of the nanocrystalline phases of 4,11-difluoro-, 2,9-dichloro- and 2,9-dichloro-6,13-dihydro-quinacridone, which were solved from powder patterns with 14-20 peaks only, and of the coordination polymer dichloro-bis(pyridine-N)copper(II).

Structure determination of organic compounds by a fit to the pair distribution function from scratch without prior indexing

A method for the ab initio crystal structure determination of organic compounds by a fit to the pair distribution function (PDF), without prior knowledge of lattice parameters and space group, has been developed. The method is called 'PDF-Global-Fit' and is implemented by extension of the program FIDEL (fit with deviating lattice parameters). The structure solution is based on a global optimization approach starting from random structural models in selected space groups. No prior indexing of the powder data is needed. The new method requires only the molecular geometry and a carefully determined PDF. The generated random structures are compared with the experimental PDF and ranked by a similarity measure based on cross-correlation functions. The most promising structure candidates are fitted to the experimental PDF data using a restricted simulated annealing structure solution approach within the program TOPAS, followed by a structure refinement against the PDF to identify the correct crystal structure. With the PDF-Global-Fit it is possible to determine the local structure of crystalline and disordered organic materials, as well as to determine the local structure of unindexable powder patterns, such as nanocrystalline samples, by a fit to the PDF. The success of the method is demonstrated using barbituric acid as an example. The crystal structure of barbituric acid form IV solved and refined by the PDF-Global-Fit is in excellent agreement with the published crystal structure data.

Comparison and evaluation of pair distribution functions, using a similarity measure based on cross-correlation functions

An approach for the comparison of pair distribution functions (PDFs) has been developed using a similarity measure based on cross-correlation functions. The PDF is very sensitive to changes in the local structure, i.e. small deviations in the structure can cause large signal shifts and significant discrepancies between the PDFs. Therefore, a comparison based on pointwise differences (e.g. R values and difference curves) may lead to the assumption that the investigated PDFs as well as the corresponding structural models are not in agreement at all, whereas a careful visual inspection of the investigated structural models and corresponding PDFs may reveal a relatively good match. To quantify the agreement of different PDFs for those cases an alternative approach is introduced: the similarity measure based on cross-correlation functions. In this paper, the power of this application of the similarity measure to the analysis of PDFs is highlighted. The similarity measure is compared with the classical R wp values as representative of the comparison based on pointwise differences as well as with the Pearson product-moment correlation coefficient, using polymorph IV of barbituric acid as an example.

Crystal Structure Prediction Methods for Organic Molecules: State of the Art

The prediction of the crystal structures that a given organic molecule is likely to form is an important theoretical problem of significant interest for the pharmaceutical and agrochemical industries, among others. As evidenced by a series of six blind tests organized over the past 2 decades, methodologies for crystal structure prediction (CSP) have witnessed substantial progress and have now reached a stage of development where they can begin to be applied to systems of practical significance. This article reviews the state of the art in general-purpose methodologies for CSP, placing them within a common framework that highlights both their similarities and their differences. The review discusses specific areas that constitute the main focus of current research efforts toward improving the reliability and widening applicability of these methodologies, and offers some perspectives for the evolution of this technology over the next decade.

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.

The PO13 crystal structure of ROY

A novel pumpkin-orange polymorph of ROY was matched to a predicted structure and refined to PXRD data.

Accuracy and reproducibility in crystal structure prediction: the curious case of ROY

Because of excessive electron delocalization, the polymorphs of ROY constitute a surprisingly challenging system for crystal structure prediction.

Challenging structure determination from powder diffraction data: two pharmaceutical salts and one cocrystal with Z′ = 2

Structure solution of molecular crystals from powder diffraction data by real-space methods becomes challenging when the total number of degrees of freedom (DoF) for molecular position, orientation and intramolecular torsions exceeds a value of 20. Here we describe the structure determination from powder diffraction data of three pharmaceutical salts or cocrystals, each with four molecules per asymmetric unit on general position: Lamivudine camphorsulfonate (1, P 21, Z=4, Z′=2; 31 DoF), Theophylline benzamide (2, P 41, Z=8, Z′=2; 23 DoF) and Aminoglutethimide camphorsulfonate hemihydrate [3, P 21, Z=4, Z′=2; 31 DoF (if the H2O molecule is ignored)]. In the salts 1 and 3 the cations and anions have two intramolecular DoF each. The molecules in the cocrystal 2 are rigid. The structures of 1 and 2 could be solved without major problems by DASH using simulated annealing. For compound 3, indexing, space group determination and Pawley fit proceeded without problems, but the structure could not be solved by the real-space method, despite extensive trials. By chance, a single crystal of 3 was obtained and the structure was determined by single-crystal X-ray diffraction. A post-analysis revealed that the failure of the real-space method could neither be explained by common sources of error such as incorrect indexing, wrong space group, phase impurities, preferred orientation, spottiness or wrong assumptions on the molecular geometry or other user errors, nor by the real-space method itself. Finally, is turned out that the structure solution failed because of problems in the extraction of the integrated reflection intensities in the Pawley fit. With suitable extracted reflection intensities the structure of 3 could be determined in a routine way.

Control and prediction of the organic solid state: a challenge to theory and experiment†

The ability of theoretical chemists to quantitatively model the weak forces between organic molecules is being exploited to predict their crystal structures and estimate their physical properties. Evolving crystal structure prediction methods are increasingly being used to aid the design of organic functional materials and provide information about thermodynamically plausible polymorphs of speciality organic materials to aid, for example, pharmaceutical development. However, the increasingly sophisticated experimental studies for detecting the range of organic solid-state behaviours provide many challenges for improving quantitative theories that form the basis for the computer modelling. It is challenging to calculate the relative thermodynamic stability of different organic crystal structures, let alone understand the kinetic effects that determine which polymorphs can be observed and are practically important. However, collaborations between experiment and theory are reaching the stage of devising experiments to target the first crystallization of new polymorphs or create novel organic molecular materials.

Improved multivariate analysis for fast and selective monitoring of structural dynamics by in situ X-ray powder diffraction

The development of two solid-state reactions, Xe absorption into MFI and molecular complex formation, where samples are affected by changes of crystal lattice due to temperature or pressure variation was structurally monitored through in situ or in operando X-ray powder diffraction experiments. Consequent variations of the peak positions prevent collective analysis of measured patterns, aiming at investigating structural changes occurring within the crystal cell. Moreover, an intrinsic and variable error in peak position is unavoidable when using the Bragg-Brentano geometry and, in some cases (sticky, bulky, aggregate samples) the sample mounting can increase the error within a dataset. Here we present a general multivariate analysis method to process in a fast and automatic way in situ XRPD data collected on charge transfer complexes and porous materials, with the capacity of disentangling peak shifts from intensity and shape variations in diffraction signals, thus allowing an efficient separation of the contribution of crystal lattice changes from structural changes. The peak shift correction allowed an improved PCA analysis that turned out to be more sensible than the traditional single pattern Rietveld analysis. The developed algorithms allowed, with respect to the traditional approach, the location of two new Xe positions into MFI with a better interpretation of the experimental data, while a much faster and more efficient recovery of the reaction coordinate was achieved in the molecular complex formation reaction.

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

3-Cyanopyridine as a bridging and terminal ligand in coordination polymers

Crystal structures of [M^IIBr₂(3-CNpy)_x]_(n) with M = Mn, Fe, Co, Ni and x = 1, 2, 4 were determined by powder diffraction. For x = 1, the 3-cyanopyridine ligand is bridging two metal atoms.

Thermal Expansion of Carbamazepine: Systematic Crystallographic Measurements Challenge Quantum Chemical Calculations

We report systematic temperature-dependent X-ray measurements on the most stable carbamazepine polymorph. This active pharmaceutical ingredient is used to demonstrate how the thermal expansion can probe certain intermolecular interactions resulting in anisotropic expansion behavior. We show that most structural features can be captured by electronic structure calculations at the quasi-harmonic approximation (QHA) provided a dispersion-corrected density functional based method is employed. The impact of thermal expansion on the phonon modes and hence free energy contributions is large enough to impact the relative stability of different polymorphs.

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.